Télécharger la plaquette Cast3M


fr en

Télécharger fsi2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  fsi2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * CAS TEST DU 91/10/04     PROVENANCE : PETI
  5.  
  6. *******************************************
  7. *    Test fsi2.dgibi: jeux de données     *
  8. *    ---------------------------------    *
  9. *******************************************
  10.  
  11. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  12. * SINON SI GRAPH DIFFERENT DE N TOUS
  13. * LES GRAPHIQUES SONT AFFICHES
  14.  
  15. GRAPH = 'N' ;
  16.  
  17. SAUT PAGE;
  18. SI (NEG GRAPH 'N') ;
  19.   OPTI ECHO 1   ;
  20.   OPTI TRAC PSC   ;
  21. SINO ;
  22.   OPTI ECHO 0   ;
  23. FINSI ;
  24.  
  25. SAUT PAGE;
  26. *******************************************
  27. *                                         *
  28. *              TEST FSI2                  *
  29. *                                         *
  30. *     CYLINDRICAL FLUID CAVITY WITHOUT    *
  31. *                FREE SURFACE             *
  32. *                                         *
  33. * Calculation of the acoustic frequencies,*
  34. * for m = 1, of cylindrical water volume  *
  35. * of radius 1.43m and height 1.039m       *
  36. *                                         *
  37. *                                         *
  38. *      P8+----------+P5                   *
  39. *        |          |                     *
  40. *        |          |                     *
  41. *        |          |                     *
  42. *        |          |                     *
  43. *        |          |                     *
  44. *      P7+----------+P6                   *
  45. *                                          *
  46. *                                         *
  47. *    The boundary conditions are          *
  48. *                                         *
  49. *      dp  |                              *
  50. *     ---- |         =  0                 *
  51. *      dz  | h = 0.                       *
  52. *                                         *
  53. *      dp  |                              *
  54. *     ---- |         =  0                 *
  55. *      dr  | r = R                        *
  56. *******************************************
  57. TEMPS;
  58.  
  60. OPTI MODE FOUR 1;
  61. OPTI ELEM QUA4;
  62.  
  63. *=========================================*
  64. *               geometry                  *
  65. *=========================================*
  66.  
  67. * Dimensions en metres
  68.  
  69. * Points
  70.  
  71. P5 = 1.43 1.039  ;
  72. P6 = 1.43 0.0    ;
  73. P7= 0. 0.        ;
  74. P8 = 0. 1.039    ;
  75.  
  76. N1 = 10       ;
  77. S4 = P6 D N1 P5 ; S5 =P7 D N1 P6 ;
  78. S6 = P7 D N1 P8 ; S7 =P8 D N1 P5 ;
  79.  
  80. WATER = DALL S4 S5 S6 S7  QUELC   ;
  81.  
  82. * OPTIO FOR TRACE
  83.  
  84. SI (NEG GRAPH 'N');
  85.   TITR ' FSI2 : MAILLAGE';
  86.   TRAC QUAL (WATER ET (0 0));
  87. FINSI;
  88.  
  89. *=========================================*
  90. *      MODE - materiau - rigidite - masse *
  91. *=========================================*
  92.  
  93. MODLIQ = MODE WATER  LIQUIDE LQU4 ;
  94.  
  95. MATLIQ = MATE  MODLIQ
  96.          RHO 1.E3 RORF 1.E3 CSON 1435. 
  97.          CREF 1435. LCAR 1 G 0.;
  98.  
  99. RIG1 = RIGI  MODLIQ MATLIQ  ;
  100. MAS1 = MASS  MODLIQ MATLIQ  ;
  101.  
  102. *=========================================*
  103. *                boundary conditions      *
  104. *=========================================*
  105.  
  106. * No explicit boundary condition
  107. * the boundary conditions are natural .
  108. *=========================================*
  109. *     calculation of the frequencies      *
  110. *                   and                   *
  111. *       extraction of some results        *
  112. *=========================================*
  113.  
  114. * Use the operator VIBR. (option PROC)
  115.  
  116. FRE1 = TABLE;
  117. FRE1.1 = 294.06;
  118. FRE1.2 = 750.56;
  119. LIST1 = PROG FRE1.1 FRE1.2 ;
  120. *
  121. RESUL = VIBR PROC LIST1 RIG1  MAS1 ;
  122.  
  123. *=========================================*
  124. *                   results               *
  125. *=========================================*
  126.  
  127. MESS ' RESULTATS ';
  128. MESS ' --------- ';
  130.  
  131. FRE2 = TABL;
  132. MOD  = TABL;
  133. DEF  = TABL;
  134. ERG  = TABL;
  135.  
  136. I = 0;
  137. REPETER BLOC1 2;
  138.   I = I + 1;
  139.   FRE2.I = RESUL . MODES . I . FREQUENCE;
  140.   ERG.I = 100 * 
  141.         (ABS ((FRE1.I - FRE2.I) / FRE1.I));
  142.   MESS ' MODE '  I     ;
  143.   MESS ' ----------';
  144.   MESS 'Frequence theorique :' FRE1.I 'Hz';
  145.   MESS 'Frequence calculee  :' FRE2.I 'Hz';
  146.   MESS '    Soit un ecart de : ' ERG.I '%';
  148. FIN BLOC1;
  149.  
  150. *=========================================*
  151. *            code validation              *
  152. *=========================================*
  153.  
  154. ERGMAX = MAXI (PROG ERG.1 ERG.2 );
  155.  
  156. SI (ERGMAX <EG 5.);
  157.    ERRE 0;
  158. SINON;
  159.    ERRE 5;
  160. FINSI;
  161.  
  163. TEMPS;
  165.  
  166. FIN;
  167.  
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales