Télécharger fsi6.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : fsi6.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * CAS TEST DU 91/10/04 PROVENANCE : PETI
  5.  
  6. *******************************************
  7. * Test fsi6.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 FSI6 *
  29. * *
  30. * PARALLELEPIPEDIQUE FLUID VOLUME WITH *
  31. * FREE SURFACE *
  32. * *
  33. * Calculation of the sloshing frequency *
  34. * of a parallelpipedique fluid volume *
  35. * Height (h) 0.1m length (a) 2m *
  36. * Breadth (b) 1.58m *
  37. * *
  38. * P2 *
  39. * _________________ *
  40. * / /| *
  41. * / / | *
  42. * / b / | *
  43. * / / | *
  44. * / a / | *
  45. * P1 /_________________/ / *
  46. * | | / *
  47. * | | / *
  48. * | h | / *
  49. * |_________________|/ *
  50. * *
  51. * *
  52. * *
  53. * *
  54. * The boundary conditions are *
  55. * *
  56. * dp | *
  57. * ---- | = 0. *
  58. * dz | h = 0. *
  59. * *
  60. * dp | *
  61. * ---- | = 0. *
  62. * dx | x = 0 , a *
  63. * *
  64. * dp | *
  65. * ---- | = 0. *
  66. * dy | y = 0 , b *
  67. *******************************************
  68.  
  69. TEMPS;
  70.  
  71. OPTI ELEM CUB8;
  72.  
  73. *=========================================*
  74. * geometry *
  75. *=========================================*
  76.  
  77. * Dimensions en metres
  78.  
  79. * Points
  80.  
  81. P1 = 0. 0. 0. ;
  82. P2 = 0. 1.58114 0. ;
  83. EYE = 5. -3. 3. ;
  84.  
  85. S1 =P1 D 15 P2 ;
  86.  
  87. FRSUR = S1 TRAN 15 (2. 0. 0.) ;
  88. WATER = FRSUR VOLU TRAN 2 (0. 0. -0.1) ;
  89.  
  90. * OPTIO FOR TRACE
  91.  
  92. SI (NEG GRAPH 'N');
  93. TITR ' FSI6 : MAILLAGE';
  94. TRAC QUAL WATER CACH EYE ;
  95. FINSI;
  96.  
  97. *=========================================*
  98. * MODE - materiau - rigidite - masse *
  99. *=========================================*
  100.  
  101. MODLIQ1=MODE WATER LIQUIDE LCU8;
  102. MODLIQ2=MODE FRSUR LIQUIDE LSU4;
  103.  
  104. MATLIQ = MATE (MODLIQ1 ET MODLIQ2)
  105. RHO 1.E3 RORF 1.E3
  106. CSON 1435. CREF 1435.
  107. LCAR 1.0 G 9.81 ;
  108.  
  109. RIG1 = RIGI (MODLIQ1 ET MODLIQ2) MATLIQ ;
  110. MAS1 = MASS (MODLIQ1 ET MODLIQ2) MATLIQ ;
  111.  
  112. *=========================================*
  113. * boundary conditions *
  114. *=========================================*
  115.  
  116. * No explicit boundary condition
  117. * the boundary conditions are natural
  118.  
  119. * Theoritical value of frequency
  120.  
  121. FRE1 = TABL ;
  122. FRE1.1 = 0.2476;
  123. FRE1.2 = 0.3132;
  124. FRE1.3 = 0.3993;
  125. FRE1.4 = 0.4952;
  126. LIST1 = PROG FRE1.1 FRE1.2 FRE1.3 FRE1.4 ;
  127.  
  128. *=========================================*
  129. * calculation of the frequencies *
  130. * and *
  131. * extraction of some results *
  132. *=========================================*
  133.  
  134. * Use of the operator VIBR. (option PROC)
  135. RESUL = VIBR PROC LIST1 RIG1 MAS1 ;
  136.  
  137. *=========================================*
  138. * results *
  139. *=========================================*
  140.  
  141. MESS ' RESULTATS ';
  142. MESS ' --------- ';
  143.  
  144. FRE2 = TABL;
  145. MOD = TABL;
  146. DEF = TABL;
  147. ERG = TABL;
  148.  
  149. I = 0;
  150. REPETER BLOC1 (dime LIST1);
  151. I = I + 1;
  152. FRE2.I = RESUL . MODES . I . FREQUENCE;
  153. ERG.I = 100 *
  154. (ABS ((FRE1.I - FRE2.I) / FRE1.I));
  155. MESS ' MODE ' I ;
  156. MESS ' ----------';
  157. MESS 'Frequence theorique :' FRE1.I 'Hz';
  158. MESS 'Frequence calculee :' FRE2.I 'Hz';
  159. MESS ' Soit un ecart de : ' ERG.I '%';
  160.  
  161. * mode shape extraction and optional trace
  162. SI (NEG GRAPH 'N');
  163. MOD.I= RESUL . MODES . I . DEFORMEE_MODALE;
  164. DEF0 = DEFO FRSUR MOD.I 0.;
  165. DEF.I= DEFO FRSUR MOD.I VERT;
  166. DEPZ = EXCO MOD.I 'UZ';
  167. TRAC CACH (DEF.I ET DEF0) EYE CACH ;
  168. TRAC DEPZ FRSUR
  169. (DEFO FRSUR MOD.I NOIR) ;
  170. FINSI;
  171. FIN BLOC1;
  172.  
  173. *=========================================*
  174. * code validation *
  175. *=========================================*
  176.  
  177. ERGMAX = MAXI
  178. (PROG ERG.1 ERG.2 ERG.3 ERG.4 ) ;
  179.  
  180. SI (ERGMAX <EG 5.0);
  181. ERRE 0;
  182. SINON;
  183. ERRE 5;
  184. FINSI;
  185.  
  186. TEMPS;
  187.  
  188. FIN;
  189.  
  190.  
  191.  
  192.  
  193.  
  194.  
  195.  
  196.  
  197.  

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