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. OPTI EPSI LINEAIRE;
  65.  
  66. *=========================================*
  67. * geometry *
  68. *=========================================*
  69. * Dimensions en metres
  70.  
  71. * Points
  72.  
  73. N1 = 20 ; N2 = 30 ;
  74. P2 = 1.43 1.039 ;
  75. P3 = 1.43 0.0 ;
  76. P4= 0. 0. ;
  77. P5 = 1.43 1.039 ;
  78. P6 = 1.43 0.0 ;
  79. P7= 0. 0. ;
  80. P8 = 0. 1.039 ;
  81.  
  82. S2 = P3 D N1 P2 ; S3 =P3 D N2 P4 ;
  83. S4 = P6 D N1 P5 ; S5 =P7 D N2 P6 ;
  84. S6 = P7 D N1 P8 ; S7 =P8 D N2 P5 ;
  85.  
  86. STEEL = S2 ET S3 ;
  87. WATER = DALL S4 S5 S6 S7 QUELC ;
  88. RAC1 = RACC 0.002 WATER STEEL ;
  89.  
  90. * OPTIO FOR TRACE
  91.  
  92. SI (NEG GRAPH 'N');
  93. TITR ' FSI3 : MAILLAGE';
  94. TRAC QUAL (WATER ET (0 0));
  95. FINSI;
  96.  
  97. *=========================================*
  98. * MODE - materiau - rigidite - masse *
  99. *=========================================*
  100. *
  101. MODCOQ = MODE STEEL MECANIQUE COQ2 ;
  102. MODLIQ = MODE WATER LIQUIDE LQU4 ;
  103. MODRAC = MODE RAC1 MECANIQUE LIQUIDE RACO;
  104.  
  105. MATCOQ = MATE MODCOQ RHO 7.8E3 YOUN 2.E11
  106. NU 0.3 ;
  107. MATLIQ = MATE (MODLIQ ET MODRAC)
  108. RHO 1.E3 RORF 1.E3 CSON 1435.
  109. CREF 1435. LCAR 1. G 0.;
  110.  
  111. CARCOQ = CARA MODCOQ EPAI 0.002;
  112. CARRAC = CARA MODRAC LIQU WATER;
  113. MATCOQ = MATCOQ ET CARCOQ;
  114. MATLIQ = MATLIQ ET CARRAC;
  115.  
  116. RIG1 = RIGI (MODCOQ ET MODLIQ ET MODRAC)
  117. (MATCOQ ET MATLIQ) ;
  118. MAS1 = MASS (MODCOQ ET MODLIQ ET MODRAC)
  119. (MATCOQ ET MATLIQ) ;
  120.  
  121. *=========================================*
  122. * boundary conditions *
  123. *=========================================*
  124.  
  125. RES1 = BLOQ DEPL ROTA P3 ;
  126. *
  127. *=========================================*
  128. * theritical value of frequency
  129.  
  130. FRE1 = 106.64 ;
  131. LIST1 = PROG FRE1 ;
  132.  
  133. *=========================================*
  134. * calculation of the frequencies *
  135. * and *
  136. * extraction of some results *
  137. *=========================================*
  138.  
  139. * Use of the operator VIBR. (option PROC)
  140.  
  141. RESUL = VIBR PROC LIST1 (RES1 ET RIG1 ) MAS1 ;
  142.  
  143. *=========================================*
  144. * results *
  145. *=========================================*
  146.  
  147. MESS ' RESULTATS ';
  148. MESS ' --------- ';
  149.  
  150. FRE2 = RESUL . MODES . 1 . FREQUENCE;
  151. ERG = 100 *
  152. (ABS ((FRE1 - FRE2 ) / FRE1 ));
  153. MESS ' MODE ' 1 ;
  154. MESS ' ----------';
  155. MESS 'Frequence theorique :' FRE1 'Hz';
  156. MESS 'Frequence calculee :' FRE2 'Hz';
  157. MESS ' Soit un ecart de : ' ERG '%';
  158.  
  159. * mode shape extraction and optional trace
  160.  
  161. SI (NEG GRAPH 'N');
  162. MOD = RESUL . MODES . 1 . DEFORMEE_MODALE;
  163. DEF0 = DEFO STEEL MOD 0. ;
  164. DEF = DEFO STEEL MOD BLEU;
  165. TRAC (DEF ET DEF0);
  166. FINSI;
  167.  
  168. *=========================================*
  169. * code validation *
  170. *=========================================*
  171.  
  172. ERGMAX = ERG ;
  173.  
  174. SI (ERGMAX <EG 5.);
  175. ERRE 0;
  176. SINON;
  177. ERRE 5;
  178. FINSI;
  179.  
  180. TEMPS;
  181.  
  182. FIN;
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  
  189.  

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