Télécharger dy_devo5.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dy_devo5.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Dy_devo5.dgibi: Jeux de donn�es *
  5. * ------------------------------------ *
  6. * *
  7. GRAPH='N';
  8. SAUT PAGE ;
  9. d_pasa = 'O' ;
  10. *******************************************************
  11. * REPONSE TRANSITOIRE D'UNE POUTRE -METHODE MODALE
  12. *******************************************************
  13. *
  14. * DESCRIPTION DU PROBLEME
  15. * -----------------------
  16. *
  17. * UNE POUTRE EST MONTEE SUR DES APPUIS ELASTIQUES .
  18. * SUPPOSANT UNE CHARGE REPARTIE TRIANGULAIRE AVCE UNE
  19. * INTENSITE VARIABLE DANS LE TEMPS, ON CHERCHE A
  20. * CALCULER LA FLECHE MAXIMALE EN DIFFERENTS INSTANTS.
  21. *
  22. * REFERENCE :
  23. * ----------
  24. *
  25. * WALTER C. HURY AND MOSHE F. RUBINSTEIN,DYNAMICS OF
  26. * STRUCTURES, (ENGLEWOOD CLIFFS: PRENTICE-HALL,
  27. * INC.,1964),P.294
  28. *******************************************************
  29. OPTION ECHO 1 ;
  30. OPTION DIME 3 ELEM SEG2 ;
  31. OPTI EPSI LINEAIRE;
  32. *
  33. * MAILLAGE
  34. *
  35. P1 = 0. 0. 0. ; P2 = 50. 0. 0. ;
  36. L1 = D 4 P1 P2 ;
  37. *
  38. * MODELE
  39. *
  40. MOD1 = MODELE L1 MECANIQUE ELASTIQUE POUT ;
  41. *
  42. * CARACTERISTIQUES MATERIELLES ET GEOMETRIQUES
  43. *
  44. MAT1 = MATER MOD1 YOUN 1.E6 NU 0.3 RHO 0.1 ;
  45. CAR1 = CARAC MOD1 SECT 1 INRY 10. INRZ 10. TORS 1.
  46. VECT (0. 1. 0.) ;
  47. *
  48. * CONDITIONS AUX LIMITES : APPUI ELASTIQUE
  49. * ET CONDITIONS DE SYMETRIE
  50. *
  51. K = 974.1 / 2 ;
  52. CL1 = APPUI UY K P1 ;
  53. CL10 = bloq uy p1 ;
  54. CL2 = BLOQ L1 RX RY UX UZ ;
  55. CL3 = BLOQ P2 RZ ;
  56. CLTOT = CL1 ET CL2 ET CL3 ;
  57. CLTOT0 = CL10 ET CL2 ET CL3 ;
  58.  
  59. * MATRICE DE RIGIDITE
  60. *
  61. RIG1 = RIGIDITE MOD1 (MAT1 ET CAR1 ) ;
  62. *
  63. * MATRICE DE MASSE
  64. *
  65. MAS1 = MASSE MOD1 (MAT1 ET CAR1 ) ;
  66. *
  67. * CALCUL DES MODES
  68. *
  69. TAB1 = VIBR 'INTERVALLE' 0. 5. (RIG1 ET CLTOT) MAS1 ;
  70. TAB2 = TAB1.MODES ;
  71. NNOR TAB2 AVEC (MOTS UY ) ;
  72. NMOD = DIME TAB2 - 2 ;
  73. IB = 0 ;
  74. REPETE BOUC NMOD ;
  75. IB = IB + 1 ;
  76. F1 = (TAB2.IB).'FREQUENCE' ;
  77. OMEG1 = (F1 * 2.) * PI ;
  78. SI (IB 'EGA' 1 ) ; OMEG1TH = 7.155 ;
  79. SINON ; OMEG1TH = 28.23 ; FINSI ;
  80. ERR1 = ((OMEG1 - OMEG1TH) / (ABS OMEG1TH)) * 100. ;
  81. 'MODE 'IB 'THEORIE 'OMEG1TH'CASTEM 'OMEG1 '%ERREUR '
  82. ERR1 ;
  83. DEP1 = (TAB2.IB).DEFORMEE_MODALE ;
  84. DEF1 = DEFORM DEP1 L1 1. ROUG ;
  85. DEF0 = DEFORM DEP1 L1 0. ;
  86. SI (NEG GRAPH 'N');
  87. TRAC (0 0 10000)(DEF0 ET DEF1) ;
  88. FINSI;
  89.  
  90. ma1 = manu poi1 tab2 . ib . point_repere ;
  91. mo1 = mode ma1 mecanique elastique modal ;
  92. ca1 = mate mo1 freq tab2 . ib . frequence
  93. mass tab2 . ib . masse_generalisee defo tab2 . ib . deformee_modale;
  94. * ride ((rig1 et cltot) * tab2 . ib . deformee_modale) ;
  95. si (&bouc ega 1) ;
  96. ma_moda = tab2. ib . 'POINT_REPERE' ;
  97. mo_moda = mo1 ; ca_moda = ca1 ;
  98. sinon ;
  99. ma_moda = ma_moda et tab2. ib . 'POINT_REPERE' ;
  100. mo_moda = mo_moda et mo1 ; ca_moda = ca_moda et ca1 ;
  101. finsi ;
  102. FIN BOUC ;
  103. *
  104.  
  105. *
  106. TABM = TAB1 ;
  107. *
  108. * -- PREPARATION DES DONNEES POUR L'OPERATEUR DYNE
  109. *
  110. *
  111. * CHARGEMENT
  112. *
  113. *
  114. * FORCE REPARTIE
  115. *
  116. FC1 = COOR 1 L1 ;
  117. FC2 = (FC1 / -50.) * 2. ;
  118. FC3 = NOMC FC2 UY ;
  119. MODB = MODEL L1 MECANIQUE ELASTIQUE BARR ;
  120. MATB = MATER MODB SECT 1.`RHO 1. ;
  121. MASB = MASSE MODB MATB ;
  122. FC4 = MASB * FC3 ;
  123. FC5 = EXCO FC4 FY FY ;
  124. LI1 = PROG 0. 1. 1.025 3.0 ;
  125. LI2 = PROG 0. 1. 0. 0. ;
  126. EV = EVOL MANU T LI1 F(T) LI2 ;
  127. CHA1 = CHAR 'FORC' FC5 EV ;
  128. TABCHAR = TABLE 'CHARGEMENT' ;
  129. *
  130. * PROJECTION SUR LA BASE MODALE
  131. *
  132. CHA2 = PJBA CHA1 TAB1 ;
  133. TABCHAR.'BASE_A' = CHA2 ;
  134. cha3 = pjba cha1 mo_moda ca_moda ;
  135. *
  136. * TABLE DE SORTIE
  137. *
  138. TABSOR = TABLE 'SORTIE' ;
  139. TABVAR = TABLE 'VARIABLE' ;
  140. TABSOR.'VARIABLE' = TABVAR ;
  141. TABVAR.'DEPLACEMENT' = VRAI ;
  142. TABVAR.'VITESSE' = FAUX ;
  143. TABVAR.'ACCELERATION' = FAUX ;
  144. *
  145. * APPEL A DYNE DE_VOGELAERE
  146. *
  147. n2n = 20 ;
  148. tabdyn = table 'PASAPAS' ;
  149. tabdyn . modele = mo_moda ;
  150. tabdyn . caracteristiques = ca_moda ;
  151. tabdyn . chargement = cha3 ;
  152. * tabdyn . modele = mod1 ;
  153. * tabdyn . caracteristiques = mat1 et car1 ;
  154. * tabdyn . chargement = cha1 ;
  155. * tabdyn . blocages_mecaniques = cl2 et cl3 ;
  156. * tabdyn . rigidite_constante = cl1 ;
  157.  
  158. si ('EGA' d_pasa 'O') ;
  159. tabdyn . temps_calcules = prog 0 pas 0.025 (4*n2n*0.025);
  160. tabdyn.dynamique = vrai ;
  161. tabdyn . convergence_forcee = faux ;
  162. pasapas tabdyn ;
  163. sinon ;
  164. tabdyn . nombre_pas = 4*n2n ;
  165. tabdyn . pas_de_temps = 0.025 ;
  166. tabdyn . pas_de_sortie = 1 ;
  167. dyne de_vogelaere tabdyn ;
  168. finsi ;
  169.  
  170. *
  171. TABRES =
  172. DYNE DE_VOGELAERE TABM TABCHAR TABSOR 80 0.025 1 ;
  173. *
  174. ou1 = prog 0. ; ou11 = prog 0. ;
  175. *d_pasa = 'O';
  176. repeter b_post n2n ;
  177. DEPU1 = RECO TABRES TABM
  178. (0. + (&b_post* 0.1)) 'DEPL' ;
  179. chdepu1 = peche tabdyn deplacements
  180. (0. + (&b_post* 0.1)) ;
  181. si ('EGA' d_pasa 'O') ; depu11 = reco mo_moda ca_moda chdepu1 ;
  182. sinon ; depu11 = chdepu1 ; finsi ;
  183. FLECHu1 = EXTR DEPu1 P2 UY ;
  184. FLECHu11 = EXTR DEPU11 P2 UY ;
  185. ou1 = ou1 et (prog flechu1) ;
  186. ou11 = ou11 et (prog flechu11) ;
  187. fin b_post ;
  188. ab1 = prog 0. pas 0.1 2. ;
  189. e1 = evol vert manu ab1 ou1 ;
  190. e2 = evol rouge manu ab1 ou11 ;
  191. SI (NEG GRAPH 'N');
  192. dess (e1 et e2 et ev) ;
  193. FINSI ;
  194. *
  195. * RECOMPOSITION POUR LES TEMPS 1.45 ET 1.9
  196. *
  197. DEP1 = RECO TABRES TABM 1.45 'DEPL' ;
  198. DEP2 = RECO TABRES TABM 1.9 'DEPL' ;
  199. *
  200. chdep1 = peche tabdyn deplacements 1.45 ;
  201. si ('EGA' d_pasa 'O') ; dep11 = reco mo_moda ca_moda chdep1 ;
  202. sinon ; dep11 = chdep1 ; finsi ;
  203. chdep2 = peche tabdyn deplacements 1.9 ;
  204. si ('EGA' d_pasa 'O') ; dep22 = reco mo_moda ca_moda chdep2 ;
  205. sinon ; dep22 = chdep2 ; finsi ;
  206.  
  207. *
  208. * COMPARAISON AVEC DES RESULTATS ANALYTIQUES
  209. *
  210. FLECH1 = EXTR DEP1 P2 UY ; FLECH11 = EXTR DEP11 P2 UY ;
  211.  
  212. FLECH1TH = 0.23382 ;
  213. ERR1 = (( FLECH1 - FLECH1TH) / (ABS FLECH1TH)) * 100.;
  214. ERR11 = (( FLECH11 - FLECH1TH) / (ABS FLECH1TH)) * 100.;
  215. MESS ' FLECHE MAXIMALE T = 1.45 ' ;
  216. MESS'THEORIE 'FLECH1TH'CASTEM 'FLECH1'%ERREUR 'ERR1;
  217. MESS'THEORIE 'FLECH1TH'CASTEM 'FLECH11'%ERREUR 'ERR11;
  218. *
  219. FLECH2 = EXTR DEP2 P2 UY ; FLECH22 = EXTR DEP22 P2 UY ;
  220. FLECH2TH = -0.24283 ;
  221. ERR2 = (( FLECH2 - FLECH2TH) / (ABS FLECH2TH)) * 100.;
  222. ERR22 = (( FLECH22 - FLECH2TH) / (ABS FLECH2TH)) * 100.;
  223. MESS ' FELECHE MAXIMALE T = 1.9 ' ;
  224. MESS'THEORIE 'FLECH2TH'CASTEM 'FLECH2'%ERREUR 'ERR2;
  225. MESS'THEORIE 'FLECH2TH'CASTEM 'FLECH22'%ERREUR 'ERR22;
  226. *
  227. SI ((MAXI ABS (PROG ERR1 ERR2 ERR11 ERR22)) < 5);
  228. ERRE 0;
  229. SINON;
  230. ERRE 5;
  231. FINSI;
  232. FIN ;
  233.  
  234.  
  235.  
  236.  
  237.  
  238.  
  239.  
  240.  
  241.  
  242.  
  243.  
  244.  
  245.  
  246.  
  247.  
  248.  

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