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.  
  59. OPTI MODE FOUR 1;
  60. OPTI ELEM QUA4;
  61. OPTI EPSI LINEAIRE;
  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 ' --------- ';
  129.  
  130. FRE2 = TABL;
  131. MOD = TABL;
  132. DEF = TABL;
  133. ERG = TABL;
  134.  
  135. I = 0;
  136. REPETER BLOC1 2;
  137. I = I + 1;
  138. FRE2.I = TIRE RESUL FREQ RANG I;
  139. ERG.I = 100 *
  140. (ABS ((FRE1.I - FRE2.I) / FRE1.I));
  141. MESS ' MODE ' I ;
  142. MESS ' ----------';
  143. MESS 'Frequence theorique :' FRE1.I 'Hz';
  144. MESS 'Frequence calculee :' FRE2.I 'Hz';
  145. MESS ' Soit un ecart de : ' ERG.I '%';
  146. FIN BLOC1;
  147.  
  148. *=========================================*
  149. * code validation *
  150. *=========================================*
  151.  
  152. ERGMAX = MAXI (PROG ERG.1 ERG.2 );
  153.  
  154. SI (ERGMAX <EG 5.);
  155. ERRE 0;
  156. SINON;
  157. ERRE 5;
  158. FINSI;
  159.  
  160. TEMPS;
  161.  
  162. FIN;
  163.  
  164.  
  165.  
  166.  
  167.  

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