Télécharger fsi3.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : fsi3.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * CAS TEST DU 91/10/04 PROVENANCE : PETI
  5.  
  6. *******************************************
  7. * Test fsi3.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 X ;
  21. SINO ;
  22. OPTI ECHO 0 ;
  23. FINSI ;
  24.  
  25. SAUT PAGE;
  26. *******************************************
  27. * *
  28. * TEST FSI3 *
  29. * *
  30. * CYLINDRICAL TANK CONTAINING WATER *
  31. * *
  32. * Calculation of the frequencies of *
  33. * cylindrical tank containing water, for *
  34. * m = 30. The radius is 1.43m and *
  35. * height is 1.039m. *
  36. * *
  37. * *
  38. * P2 *
  39. * + *
  40. * | *
  41. * | *
  42. * fluid | *
  43. * | *
  44. * | *
  45. * | *
  46. * --------------+ *
  47. * P4 P3 *
  48. * *
  49. * The boundary condition : p3 is fixed *
  50. * *
  51. * dp | *
  52. * ---- | = 0. *
  53. * dz | h = 0. , 1.43 *
  54. * *
  55. * dp | *
  56. * ---- | = 0. *
  57. * dr | r = R *
  58. *******************************************
  59.  
  60. TEMPS;
  61.  
  62. OPTI MODE FOUR 10 ;
  63. OPTI ELEM QUA4;
  64.  
  65. *=========================================*
  66. * geometry *
  67. *=========================================*
  68. * Dimensions en metres
  69.  
  70. * Points
  71.  
  72. N1 = 20 ; N2 = 30 ;
  73. P2 = 1.43 1.039 ;
  74. P3 = 1.43 0.0 ;
  75. P4= 0. 0. ;
  76. P5 = 1.43 1.039 ;
  77. P6 = 1.43 0.0 ;
  78. P7= 0. 0. ;
  79. P8 = 0. 1.039 ;
  80.  
  81. S2 = P3 D N1 P2 ; S3 =P3 D N2 P4 ;
  82. S4 = P6 D N1 P5 ; S5 =P7 D N2 P6 ;
  83. S6 = P7 D N1 P8 ; S7 =P8 D N2 P5 ;
  84.  
  85. STEEL = S2 ET S3 ;
  86. WATER = DALL S4 S5 S6 S7 QUELC ;
  87. RAC1 = RACC 0.002 WATER STEEL ;
  88.  
  89. * OPTIO FOR TRACE
  90.  
  91. SI (NEG GRAPH 'N');
  92. TITR ' FSI3 : MAILLAGE';
  93. TRAC QUAL (WATER ET (0 0));
  94. FINSI;
  95.  
  96. *=========================================*
  97. * MODE - materiau - rigidite - masse *
  98. *=========================================*
  99. *
  100. MODCOQ = MODE STEEL MECANIQUE COQ2 ;
  101. MODLIQ = MODE WATER LIQUIDE LQU4 ;
  102. MODRAC = MODE RAC1 MECANIQUE LIQUIDE RACO;
  103.  
  104. MATCOQ = MATE MODCOQ RHO 7.8E3 YOUN 2.E11
  105. NU 0.3 ;
  106. MATLIQ = MATE (MODLIQ ET MODRAC)
  107. RHO 1.E3 RORF 1.E3 CSON 1435.
  108. CREF 1435. LCAR 1. G 0.;
  109.  
  110. CARCOQ = CARA MODCOQ EPAI 0.002;
  111. CARRAC = CARA MODRAC LIQU WATER;
  112. MATCOQ = MATCOQ ET CARCOQ;
  113. MATLIQ = MATLIQ ET CARRAC;
  114.  
  115. RIG1 = RIGI (MODCOQ ET MODLIQ ET MODRAC)
  116. (MATCOQ ET MATLIQ) ;
  117. MAS1 = MASS (MODCOQ ET MODLIQ ET MODRAC)
  118. (MATCOQ ET MATLIQ) ;
  119.  
  120. *=========================================*
  121. * boundary conditions *
  122. *=========================================*
  123.  
  124. RES1 = BLOQ DEPL ROTA P3 ;
  125. *
  126. *=========================================*
  127. * theritical value of frequency
  128.  
  129. FRE1 = 106.64 ;
  130. LIST1 = PROG FRE1 ;
  131.  
  132. *=========================================*
  133. * calculation of the frequencies *
  134. * and *
  135. * extraction of some results *
  136. *=========================================*
  137.  
  138. * Use of the operator VIBR. (option PROC)
  139.  
  140. RESUL = VIBR PROC LIST1 (RES1 ET RIG1 ) MAS1 ;
  141.  
  142. *=========================================*
  143. * results *
  144. *=========================================*
  145.  
  146. MESS ' RESULTATS ';
  147. MESS ' --------- ';
  148.  
  149. FRE2 = RESUL . MODES . 1 . FREQUENCE;
  150. ERG = 100 *
  151. (ABS ((FRE1 - FRE2 ) / FRE1 ));
  152. MESS ' MODE ' 1 ;
  153. MESS ' ----------';
  154. MESS 'Frequence theorique :' FRE1 'Hz';
  155. MESS 'Frequence calculee :' FRE2 'Hz';
  156. MESS ' Soit un ecart de : ' ERG '%';
  157.  
  158. * mode shape extraction and optional trace
  159.  
  160. SI (NEG GRAPH 'N');
  161. MOD = RESUL . MODES . 1 . DEFORMEE_MODALE;
  162. DEF0 = DEFO STEEL MOD 0. ;
  163. DEF = DEFO STEEL MOD BLEU;
  164. TRAC (DEF ET DEF0);
  165. FINSI;
  166.  
  167. *=========================================*
  168. * code validation *
  169. *=========================================*
  170.  
  171. ERGMAX = ERG ;
  172.  
  173. SI (ERGMAX <EG 5.);
  174. ERRE 0;
  175. SINON;
  176. ERRE 5;
  177. FINSI;
  178.  
  179. TEMPS;
  180.  
  181. FIN;
  182.  
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  
  189.  
  190.  

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