Télécharger fsi5.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : fsi5.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * CAS TEST DU 91/10/04 PROVENANCE : PETI
  5.  
  6. *******************************************
  7. * Test fsi5.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 TRAC PSC ;
  20. FINSI ;
  21.  
  22. SAUT PAGE;
  23. *******************************************
  24. * TEST FSI5 *
  25. * *
  26. * ANNULAR FLUID LAYER WITH FREE SURFACE *
  27. * *
  28. * Calculation of the frequencies of an *
  29. * annular water layer with free surface, *
  30. * for m= 1 *
  31. * The mean radius is 0.5m and the *
  32. * thickness of annulus is 5mm *
  33. * *
  34. * P5_________ *
  35. * | | *
  36. * | | *
  37. * | | *
  38. * | fluid | *
  39. * | | *
  40. * | | *
  41. * | | *
  42. * P6 ----------- *
  43. * *
  44. * The boundary conditions are *
  45. * *
  46. * dp | *
  47. * ---- | = 0. *
  48. * dz | h = 0. *
  49. * *
  50. * dp | *
  51. * ---- | = 0. *
  52. * dr | r = Ri , Ro *
  53. *******************************************
  54.  
  55. OPTI MODE FOUR 1 ;
  56. OPTI ELEM QUA4;
  57. OPTI EPSI LINEAIRE;
  58. *=========================================*
  59. * geometry *
  60. *=========================================*
  61.  
  62. * Dimensions en metres
  63.  
  64. * Points
  65.  
  66. THICK = 5.0 ;
  67. N1 = 20 ;
  68. R = 500. ;
  69. P5 = (R - (THICK * 0.5 )) 1000.0 ;
  70. P6 = (R - (THICK * 0.5 )) 0.0 ;
  71.  
  72. S5 =P5 D N1 P6 ;
  73.  
  74. WATER = S5 TRAN 1 (THICK 0. ) ;
  75. FRSUR = COTE WATER 4 ;
  76.  
  77. * OPTIO FOR TRACE
  78.  
  79. SI (NEG GRAPH 'N');
  80. TITR ' FSI5 : MAILLAGE';
  81. TRAC QUAL (WATER ET (0 0));
  82. FINSI;
  83.  
  84. *=========================================*
  85. * MODE - materiau - rigidite - masse *
  86. *=========================================*
  87.  
  88. MODLIQ =MODE WATER LIQUIDE LQU4 ;
  89. MODSUR =MODE FRSUR LIQUIDE LSU2;
  90.  
  91.  
  92. MATLIQ = MATE (MODLIQ ET MODSUR)
  93. RHO 1.E-9 RORF 1.E-9
  94. CSON 1435.E3 CREF 1435.E3
  95. LCAR 1000. G 9810.;
  96.  
  97. RIG1 = RIGI (MODLIQ ET MODSUR) MATLIQ ;
  98. MAS1 = MASS (MODLIQ ET MODSUR) MATLIQ ;
  99.  
  100. *=========================================*
  101. * boundary conditions *
  102. *=========================================*
  103. * No explicit boundary condition
  104. * the boundary conditions are natural
  105.  
  106. * theoretical value of frequency
  107. FRE1 = 0.6922;
  108. LIST1 = PROG FRE1 ;
  109.  
  110. *=========================================*
  111. * calculation of the frequencies *
  112. * and *
  113. * extraction of some results *
  114. *=========================================*
  115.  
  116. * Use of the operator VIBR (option PROC)
  117. RESUL2 = VIBR 'PROC' LIST1 RIG1 MAS1 ;
  118.  
  119. * Use of the operator VIBR (option IRAM)
  120. * opti impi 5;
  121. RESUL3 = VIBR 'IRAM' 0.001 1 RIG1 MAS1 ;
  122. * opti impi 0;
  123.  
  124. *=========================================*
  125. * results *
  126. *=========================================*
  127.  
  128. OPTI ECHO 0;
  129. MESS ' RESULTATS ';
  130. MESS ' ----------------------------------------------------------';
  131. MESS ' | Frequence propre en Hz (erreur relative en %) ';
  132. MESS ' Mode | Theorie | VIBR PROCH | VIBR IRAM ';
  133. MESS ' ----------------------------------------------------------';
  134.  
  135. FRE2 = RESUL2 . MODES . 1 . FREQUENCE;
  136. ERG2 = 100 * (ABS ((FRE1 - FRE2 ) / FRE1 ));
  137. FRE3 = RESUL3 . MODES . 1 . FREQUENCE;
  138. ERG3 = 100 * (ABS ((FRE1 - FRE3 ) / FRE1 ));
  139. chacha = chai 1 *5 ' |' FORMAT '(F9.6)' FRE1
  140. ' |' FORMAT '(F9.6)' FRE2 ' (' FORMAT '(F5.3)' ERG2 '%)'
  141. ' |' FORMAT '(F9.6)' FRE3 ' (' FORMAT '(F5.3)' ERG3 '%)';
  142. MESS chacha;
  143.  
  144. * mode shape extraction and optional trace
  145.  
  146. SI (NEG GRAPH 'N');
  147. MOD2 = RESUL2 . MODES . I . DEFORMEE_MODALE;
  148. DEF0 = DEFO FRSUR MOD2 0. ;
  149. DEF2 = DEFO FRSUR MOD2 'BLEU';
  150. TRAC (DEF2 ET DEF0);
  151. DEF3 = DEFO FRSUR MOD3 'ROSE';
  152. TRAC (DEF3 ET DEF0);
  153. FINSI;
  154.  
  155. MESS ' ----------------------------------------------------------';
  156. OPTI ECHO 1;
  157.  
  158. *=========================================*
  159. * code validation *
  160. *=========================================*
  161.  
  162. ERG2MAX = ERG2 ;
  163. ERG3MAX = ERG3 ;
  164.  
  165. SI (ERG2MAX > 4.5);
  166. MESS 'ERREUR VIBR OPTION PROCHE';
  167. ERRE 5;
  168. FINSI;
  169. SI (ERG3MAX > 4.5);
  170. MESS 'ERREUR VIBR OPTION IRAM';
  171. ERRE 5;
  172. FINSI;
  173. ERRE 0;
  174.  
  175.  
  176. FIN;
  177.  
  178.  
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  
  185.  

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