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.
  72. (RIG1 ) MAS1 'TBAS';
  73. *
  74. * COMPARAISON AVEC LES RESULTATS ANALYTIQUES
  75. TAB2 = TAB1.MODES ;
  76. OMEG1 = (((TAB2.1).FREQUENCE) * 2. ) * PI ;
  77. OMEG1TH = 446.18 ;ER01 =((OMEG1 - OMEG1TH) / (ABS OMEG1TH))*100.;
  78. OMEG2 = (((TAB2.2).FREQUENCE) * 2. ) * PI ;
  79. OMEG2TH = 1893.68 ;ER02 =((OMEG2 - OMEG2TH) / (ABS OMEG2TH))*100.;
  80. MESS ' PREMIER MODE ' ;
  81. MESS 'THEORIQUE ' OMEG1TH 'CASTEM 'OMEG1 '%ERREUR ' ER01 ;
  82. MESS ' DEUXIEME MODE ' ;
  83. MESS 'THEORIQUE ' OMEG2TH 'CASTEM 'OMEG2 '%ERREUR ' ER02 ;
  84. *
  85. * DEFORMEES
  86. DEF1 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 1. ROUG ;
  87. DEF2 = DEFORM (TAB2.2).DEFORMEE_MODALE L1 1. ROUG ;
  88. DEF0 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 0. VERT ;
  89. SI (NEG GRAPH 'N');
  90. TRAC (DEF0 ET DEF1) ;
  91. TRAC (DEF0 ET DEF2) ;
  92. FINSI;
  93. *
  94. * CALCUL DES MODES (MODE FOURIER 1)
  95. OPTIO MODE FOUR 1 ;
  96. RIG1 = RIGIDITE MOD1 (MAT1 ET CAR1) ;
  97. MAS1 = MASSE MOD1 (MAT1 ET CAR1) ;
  98. TAB1 = VIBR 'INTERVALE' 1. 1000.
  99. (RIG1 ) MAS1 'TBAS';
  100. *
  101. * COMPARAISON AVEC LES RESULTATS ANALYTIQUES
  102. TAB2 = TAB1.MODES ;
  103. OMEG1 = (((TAB2.1).FREQUENCE) * 2. ) * PI ;
  104. OMEG1TH = 1008.78 ;ER11 =((OMEG1 - OMEG1TH) / (ABS OMEG1TH))*100 ;
  105. OMEG2 = (((TAB2.2).FREQUENCE) * 2. ) * PI ;
  106. OMEG2TH = 2937.86 ;ER12 =((OMEG2 - OMEG2TH) / (ABS OMEG2TH))*100 ;
  107. MESS ' PREMIER MODE ' ;
  108. MESS 'THEORIQUE ' OMEG1TH 'CASTEM 'OMEG1 '%ERREUR ' ER11 ;
  109. MESS ' DEUXIEME MODE ' ;
  110. MESS 'THEORIQUE ' OMEG2TH 'CASTEM 'OMEG2 '%ERREUR ' ER12 ;
  111. *
  112. * DEFORMEES
  113. DEF1 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 1. ROUG ;
  114. DEF2 = DEFORM (TAB2.2).DEFORMEE_MODALE L1 1. ROUG ;
  115. DEF0 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 0. VERT ;
  116. SI (NEG GRAPH 'N');
  117. TRAC (DEF0 ET DEF1) ;
  118. TRAC (DEF0 ET DEF2) ;
  119. FINSI;
  120. *
  121. * CALCUL DES MODES (MODE FOURIER 2)
  122. OPTIO MODE FOUR 2 ;
  123. RIG1 = RIGIDITE MOD1 (MAT1 ET CAR1) ;
  124. MAS1 = MASSE MOD1 (MAT1 ET CAR1) ;
  125. TAB1 = VIBR 'INTERVALE' 1. 1000.
  126. (RIG1 ) MAS1 'TBAS';
  127. *
  128. * COMPARAISON AVEC LES RESULTATS ANALYTIQUES
  129. TAB2 = TAB1.MODES ;
  130. OMEG1 = (((TAB2.1).FREQUENCE) * 2. ) * PI ;
  131. OMEG1TH = 258.14 ;ER21 =((OMEG1 - OMEG1TH) / (ABS OMEG1TH))*100 ;
  132. OMEG2 = (((TAB2.2).FREQUENCE) * 2. ) * PI ;
  133. OMEG2TH = 1732.43 ;ER22 =((OMEG2 - OMEG2TH) / (ABS OMEG2TH))*100 ;
  134. MESS ' PREMIER MODE ' ;
  135. MESS 'THEORIQUE ' OMEG1TH 'CASTEM 'OMEG1 '%ERREUR ' ER21 ;
  136. MESS ' DEUXIEME MODE ' ;
  137. MESS 'THEORIQUE ' OMEG2TH 'CASTEM 'OMEG2 '%ERREUR ' ER22 ;
  138. *
  139. * DEFORMEES
  140. DEF1 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 1. ROUG ;
  141. DEF2 = DEFORM (TAB2.2).DEFORMEE_MODALE L1 1. ROUG ;
  142. DEF0 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 0. VERT ;
  143. SI (NEG GRAPH 'N');
  144. TRAC (DEF0 ET DEF1) ;
  145. TRAC (DEF0 ET DEF2) ;
  146. FINSI;
  147. *
  148. * CALCUL DES MODES (MODE FOURIER 3)
  149. OPTIO MODE FOUR 3 ;
  150. RIG1 = RIGIDITE MOD1 (MAT1 ET CAR1) ;
  151. MAS1 = MASSE MOD1 (MAT1 ET CAR1) ;
  152. TAB1 = VIBR 'INTERVALE' 1. 1000.
  153. (RIG1 ) MAS1 'TBAS';
  154. *
  155. * COMPARAISON AVEC LES RESULTATS ANALYTIQUES
  156. TAB2 = TAB1.MODES ;
  157. OMEG1 = (((TAB2.1).FREQUENCE) * 2. ) * PI ;
  158. OMEG1TH = 601.23 ;ER31 =((OMEG1 - OMEG1TH) / (ABS OMEG1TH))*100 ;
  159. OMEG2 = (((TAB2.2).FREQUENCE) * 2. ) * PI ;
  160. OMEG2TH = 2601.11 ;ER32 =((OMEG2 - OMEG2TH) / (ABS OMEG2TH))*100;
  161. MESS ' PREMIER MODE ' ;
  162. MESS 'THEORIQUE ' OMEG1TH 'CASTEM 'OMEG1 '%ERREUR ' ER31 ;
  163. MESS ' DEUXIEME MODE ' ;
  164. MESS 'THEORIQUE ' OMEG2TH 'CASTEM 'OMEG2 '%ERREUR ' ER32 ;
  165. *
  166. * DEFORMEES
  167. DEF1 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 1. ROUG ;
  168. DEF2 = DEFORM (TAB2.2).DEFORMEE_MODALE L1 1. ROUG ;
  169. DEF0 = DEFORM (TAB2.1).DEFORMEE_MODALE L1 0. VERT ;
  170. SI (NEG GRAPH 'N');
  171. TRAC (DEF0 ET DEF1) ;
  172. TRAC (DEF0 ET DEF2) ;
  173. FINSI;
  174.  
  175.  
  176. SI ((MAXI ABS (PROG ER01 ER02 ER11 ER12 ER21 ER22 ER31 ER32)) < 2.5);
  177. ERRE 0;
  178. SINON;
  179. ERRE 5;
  180. FINSI;
  181.  
  182. FIN ;
  183.  
  184.  
  185.  
  186.  
  187.  

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