Télécharger dyna10.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dyna10.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Dyna10.dgibi: Jeux de données *
  5. * --------------------------------- *
  6. * *
  7.  
  8. *******************************************************
  9. * Test dyna10.dgibi: jeux de données *
  10. * --------------------------------- *
  11. *******************************************************
  12.  
  13. * SI GRAPH = N PAS DE GRAPHIQUE AFFICHE
  14. * SINON SI GRAPH DIFFERENT DE N TOUS
  15. * LES GRAPHIQUES SONT AFFICHES
  16.  
  17. GRAPH = 'N' ;
  18.  
  19. SAUT PAGE;
  20. SI (NEG GRAPH 'N') ;
  21. OPTI ECHO 1 ;
  22. OPTI TRAC PSC ;
  23. SINO ;
  24. OPTI ECHO 0 ;
  25. FINSI ;
  26.  
  27. SAUT PAGE;
  28. *******************************************************************
  29. * DYNA10
  30. *******************************************************************
  31. * MODES PROPES D'UNE PLAQUE CIRCULAIRE ,SYMETRIE CYCLIQUE
  32.  
  33. * DESCRIPTION DU PROBLEME
  34. * ------------------------
  35. * IL S'AGIT DE CALCULER LES FREQUENCES NATURELLES D'UNE PLAQUE
  36. * CIRCULAIRE LIBRE
  37. *
  38. * REFERENCE :
  39. * ----------
  40. * S. TIMOSHENKO,D. H. YOUNG AND WEAVER,JR.,VIBRATION PROBLEMES IN
  41. * ENGINEERING ,4TH ED. (NEW YORK :WILEY & SONS,1974) P. 501
  42. *
  43. ***********************************************************************
  44. OPTI DIME 2 ELEM QUA4 ;
  45. OPTI EPSI LINEAIRE;
  46. *
  47. * MAILLAGE
  48. P1 = 0. 0. ; P2 = 11. 0. ;
  49. L1 = D 11 P1 P2 ;
  50.  
  51. SI (NEG GRAPH 'N');
  52. TRAC 'QUAL' L1 ;
  53. FINSI;
  54.  
  55. OPTI MODE FOUR 0 ;
  56. *
  57. * MODELE
  58. MOD1 = MODELE L1 MECANIQUE ELASTIQUE COQ2 ;
  59. *
  60. * PROPRIETES DE MATERIAUX
  61. MAT1 = MATE MOD1 YOUN 1.E7 NU 0.3 RHO 0.0002588 ;
  62. CAR1 = CARA MOD1 EPAI 0.1 ;
  63. *
  64. * MATRICE DE RIGIDITE
  65. RIG1 = RIGIDITE MOD1 (MAT1 ET CAR1) ;
  66. *
  67. * MATRICE DE MASSE
  68. MAS1 = MASSE MOD1 (MAT1 ET CAR1) ;
  69. *
  70. * CALCUL DES MODES
  71. TAB1 = VIBR 'INTERVALE' 1. 1000. (RIG1 ) MAS1 ;
  72. *
  73. * COMPARAISON AVEC LES RESULTATS ANALYTIQUES
  74. TAB2 = TAB1.MODES ;
  75. OMEG1 = (((TAB2.1).FREQUENCE) * 2. ) * PI ;
  76. OMEG1TH = 446.18 ;ER01 =((OMEG1 - OMEG1TH) / (ABS OMEG1TH))*100.;
  77. OMEG2 = (((TAB2.2).FREQUENCE) * 2. ) * PI ;
  78. OMEG2TH = 1893.68 ;ER02 =((OMEG2 - OMEG2TH) / (ABS OMEG2TH))*100.;
  79. MESS ' PREMIER MODE ' ;
  80. MESS 'THEORIQUE ' OMEG1TH 'CASTEM 'OMEG1 '%ERREUR ' ER01 ;
  81. MESS ' DEUXIEME MODE ' ;
  82. MESS 'THEORIQUE ' OMEG2TH 'CASTEM 'OMEG2 '%ERREUR ' ER02 ;
  83. *
  84. * DEFORMEES
  85. DEF1 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 1. ROUG ;
  86. DEF2 = DEFORM (TAB2.2).DEFORMEE_MODALE L1 1. ROUG ;
  87. DEF0 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 0. VERT ;
  88. SI (NEG GRAPH 'N');
  89. TRAC (DEF0 ET DEF1) ;
  90. TRAC (DEF0 ET DEF2) ;
  91. FINSI;
  92. *
  93. * CALCUL DES MODES (MODE FOURIER 1)
  94. OPTIO MODE FOUR 1 ;
  95. RIG1 = RIGIDITE MOD1 (MAT1 ET CAR1) ;
  96. MAS1 = MASSE MOD1 (MAT1 ET CAR1) ;
  97. TAB1 = VIBR 'INTERVALE' 1. 1000. (RIG1 ) MAS1 ;
  98. *
  99. * COMPARAISON AVEC LES RESULTATS ANALYTIQUES
  100. TAB2 = TAB1.MODES ;
  101. OMEG1 = (((TAB2.1).FREQUENCE) * 2. ) * PI ;
  102. OMEG1TH = 1008.78 ;ER11 =((OMEG1 - OMEG1TH) / (ABS OMEG1TH))*100 ;
  103. OMEG2 = (((TAB2.2).FREQUENCE) * 2. ) * PI ;
  104. OMEG2TH = 2937.86 ;ER12 =((OMEG2 - OMEG2TH) / (ABS OMEG2TH))*100 ;
  105. MESS ' PREMIER MODE ' ;
  106. MESS 'THEORIQUE ' OMEG1TH 'CASTEM 'OMEG1 '%ERREUR ' ER11 ;
  107. MESS ' DEUXIEME MODE ' ;
  108. MESS 'THEORIQUE ' OMEG2TH 'CASTEM 'OMEG2 '%ERREUR ' ER12 ;
  109. *
  110. * DEFORMEES
  111. DEF1 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 1. ROUG ;
  112. DEF2 = DEFORM (TAB2.2).DEFORMEE_MODALE L1 1. ROUG ;
  113. DEF0 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 0. VERT ;
  114. SI (NEG GRAPH 'N');
  115. TRAC (DEF0 ET DEF1) ;
  116. TRAC (DEF0 ET DEF2) ;
  117. FINSI;
  118. *
  119. * CALCUL DES MODES (MODE FOURIER 2)
  120. OPTIO MODE FOUR 2 ;
  121. RIG1 = RIGIDITE MOD1 (MAT1 ET CAR1) ;
  122. MAS1 = MASSE MOD1 (MAT1 ET CAR1) ;
  123. TAB1 = VIBR 'INTERVALE' 1. 1000. (RIG1 ) MAS1 ;
  124. *
  125. * COMPARAISON AVEC LES RESULTATS ANALYTIQUES
  126. TAB2 = TAB1.MODES ;
  127. OMEG1 = (((TAB2.1).FREQUENCE) * 2. ) * PI ;
  128. OMEG1TH = 258.14 ;ER21 =((OMEG1 - OMEG1TH) / (ABS OMEG1TH))*100 ;
  129. OMEG2 = (((TAB2.2).FREQUENCE) * 2. ) * PI ;
  130. OMEG2TH = 1732.43 ;ER22 =((OMEG2 - OMEG2TH) / (ABS OMEG2TH))*100 ;
  131. MESS ' PREMIER MODE ' ;
  132. MESS 'THEORIQUE ' OMEG1TH 'CASTEM 'OMEG1 '%ERREUR ' ER21 ;
  133. MESS ' DEUXIEME MODE ' ;
  134. MESS 'THEORIQUE ' OMEG2TH 'CASTEM 'OMEG2 '%ERREUR ' ER22 ;
  135. *
  136. * DEFORMEES
  137. DEF1 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 1. ROUG ;
  138. DEF2 = DEFORM (TAB2.2).DEFORMEE_MODALE L1 1. ROUG ;
  139. DEF0 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 0. VERT ;
  140. SI (NEG GRAPH 'N');
  141. TRAC (DEF0 ET DEF1) ;
  142. TRAC (DEF0 ET DEF2) ;
  143. FINSI;
  144. *
  145. * CALCUL DES MODES (MODE FOURIER 3)
  146. OPTIO MODE FOUR 3 ;
  147. RIG1 = RIGIDITE MOD1 (MAT1 ET CAR1) ;
  148. MAS1 = MASSE MOD1 (MAT1 ET CAR1) ;
  149. TAB1 = VIBR 'INTERVALE' 1. 1000. (RIG1 ) MAS1 ;
  150. *
  151. * COMPARAISON AVEC LES RESULTATS ANALYTIQUES
  152. TAB2 = TAB1.MODES ;
  153. OMEG1 = (((TAB2.1).FREQUENCE) * 2. ) * PI ;
  154. OMEG1TH = 601.23 ;ER31 =((OMEG1 - OMEG1TH) / (ABS OMEG1TH))*100 ;
  155. OMEG2 = (((TAB2.2).FREQUENCE) * 2. ) * PI ;
  156. OMEG2TH = 2601.11 ;ER32 =((OMEG2 - OMEG2TH) / (ABS OMEG2TH))*100;
  157. MESS ' PREMIER MODE ' ;
  158. MESS 'THEORIQUE ' OMEG1TH 'CASTEM 'OMEG1 '%ERREUR ' ER31 ;
  159. MESS ' DEUXIEME MODE ' ;
  160. MESS 'THEORIQUE ' OMEG2TH 'CASTEM 'OMEG2 '%ERREUR ' ER32 ;
  161. *
  162. * DEFORMEES
  163. DEF1 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 1. ROUG ;
  164. DEF2 = DEFORM (TAB2.2).DEFORMEE_MODALE L1 1. ROUG ;
  165. DEF0 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 0. VERT ;
  166. SI (NEG GRAPH 'N');
  167. TRAC (DEF0 ET DEF1) ;
  168. TRAC (DEF0 ET DEF2) ;
  169. FINSI;
  170.  
  171.  
  172. SI ((MAXI ABS (PROG ER01 ER02 ER11 ER12 ER21 ER22 ER31 ER32)) < 2.5);
  173. ERRE 0;
  174. SINON;
  175. ERRE 5;
  176. FINSI;
  177.  
  178. FIN ;
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  
  185.  

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