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.  
  62. *=========================================*
  63. * geometry *
  64. *=========================================*
  65.  
  66. * Dimensions en metres
  67.  
  68. * Points
  69.  
  70. P5 = 1.43 1.039 ;
  71. P6 = 1.43 0.0 ;
  72. P7= 0. 0. ;
  73. P8 = 0. 1.039 ;
  74.  
  75. N1 = 10 ;
  76. S4 = P6 D N1 P5 ; S5 =P7 D N1 P6 ;
  77. S6 = P7 D N1 P8 ; S7 =P8 D N1 P5 ;
  78.  
  79. WATER = DALL S4 S5 S6 S7 QUELC ;
  80.  
  81. * OPTIO FOR TRACE
  82.  
  83. SI (NEG GRAPH 'N');
  84. TITR ' FSI2 : MAILLAGE';
  85. TRAC QUAL (WATER ET (0 0));
  86. FINSI;
  87.  
  88. *=========================================*
  89. * MODE - materiau - rigidite - masse *
  90. *=========================================*
  91.  
  92. MODLIQ = MODE WATER LIQUIDE LQU4 ;
  93.  
  94. MATLIQ = MATE MODLIQ
  95. RHO 1.E3 RORF 1.E3 CSON 1435.
  96. CREF 1435. LCAR 1 G 0.;
  97.  
  98. RIG1 = RIGI MODLIQ MATLIQ ;
  99. MAS1 = MASS MODLIQ MATLIQ ;
  100.  
  101. *=========================================*
  102. * boundary conditions *
  103. *=========================================*
  104.  
  105. * No explicit boundary condition
  106. * the boundary conditions are natural .
  107. *=========================================*
  108. * calculation of the frequencies *
  109. * and *
  110. * extraction of some results *
  111. *=========================================*
  112.  
  113. * Use the operator VIBR. (option PROC)
  114.  
  115. FRE1 = TABLE;
  116. FRE1.1 = 294.06;
  117. FRE1.2 = 750.56;
  118. LIST1 = PROG FRE1.1 FRE1.2 ;
  119. *
  120. RESUL = VIBR PROC LIST1 RIG1 MAS1 ;
  121.  
  122. *=========================================*
  123. * results *
  124. *=========================================*
  125.  
  126. MESS ' RESULTATS ';
  127. MESS ' --------- ';
  128.  
  129. FRE2 = TABL;
  130. MOD = TABL;
  131. DEF = TABL;
  132. ERG = TABL;
  133.  
  134. I = 0;
  135. REPETER BLOC1 2;
  136. I = I + 1;
  137. FRE2.I = RESUL . MODES . I . FREQUENCE;
  138. ERG.I = 100 *
  139. (ABS ((FRE1.I - FRE2.I) / FRE1.I));
  140. MESS ' MODE ' I ;
  141. MESS ' ----------';
  142. MESS 'Frequence theorique :' FRE1.I 'Hz';
  143. MESS 'Frequence calculee :' FRE2.I 'Hz';
  144. MESS ' Soit un ecart de : ' ERG.I '%';
  145. FIN BLOC1;
  146.  
  147. *=========================================*
  148. * code validation *
  149. *=========================================*
  150.  
  151. ERGMAX = MAXI (PROG ERG.1 ERG.2 );
  152.  
  153. SI (ERGMAX <EG 5.);
  154. ERRE 0;
  155. SINON;
  156. ERRE 5;
  157. FINSI;
  158.  
  159. TEMPS;
  160.  
  161. FIN;
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  

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