Télécharger la plaquette Cast3M


fr en

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. *
  32. *  MAILLAGE 
  33. *
  34.  P1 = 0. 0. 0. ; P2 = 50. 0. 0. ;
  35.  L1 = D 4 P1 P2  ;
  36. * 
  37. * MODELE
  38. *
  39.  MOD1 = MODELE L1 MECANIQUE ELASTIQUE POUT ;
  40. *
  41. *  CARACTERISTIQUES MATERIELLES ET GEOMETRIQUES 
  42. *
  43.  MAT1 = MATER MOD1 YOUN 1.E6 NU 0.3  RHO 0.1 ;
  44.  CAR1 = CARAC MOD1 SECT 1 INRY 10. INRZ 10. TORS 1. 
  45.   VECT (0. 1. 0.)  ;
  46. *
  47. *   CONDITIONS AUX LIMITES : APPUI ELASTIQUE 
  48. *          ET CONDITIONS DE SYMETRIE
  49. *
  50.   K = 974.1 / 2  ;
  51.  CL1 = APPUI UY K  P1  ;
  52.   CL10 = bloq uy p1 ;
  53.  CL2 = BLOQ L1 RX RY UX UZ ;
  54.  CL3 = BLOQ P2 RZ    ;
  55.  CLTOT = CL1  ET CL2 ET CL3 ;
  56.   CLTOT0 = CL10  ET CL2 ET CL3 ;
  57.  
  58. *  MATRICE DE  RIGIDITE 
  59. * 
  60.  RIG1 = RIGIDITE MOD1 (MAT1 ET CAR1 ) ;
  61. *
  62. *  MATRICE DE MASSE 
  63. *
  64.  MAS1 = MASSE MOD1 (MAT1 ET CAR1 ) ;
  65. *
  66. *  CALCUL DES MODES 
  67. * 
  68.  TAB1 = VIBR  'INTERVALLE' 0. 5. (RIG1 ET CLTOT) MAS1     ;
  69. TAB1.MODES = NNOR TAB1.MODES AVEC (MOTS UY ) ;
  70. TAB2 = TAB1.MODES ;
  71. NMOD = DIME TAB2  - 2 ;
  72. IB = 0 ;
  73.  REPETE BOUC NMOD ;
  74.   IB = IB + 1   ;
  75.   F1 = (TAB2.IB).'FREQUENCE'  ;
  76.   OMEG1 = (F1 * 2.) * PI ;
  77.   SI (IB 'EGA' 1 ) ; OMEG1TH = 7.155 ;
  78.   SINON ; OMEG1TH = 28.23   ; FINSI ;
  79.  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.  tabdyn . hypothese_deformations = 'LINEAIRE' ;
  163.  pasapas tabdyn ;
  164. sinon ;
  165.  tabdyn . nombre_pas = 4*n2n ;
  166.  tabdyn . pas_de_temps = 0.025 ;
  167.  tabdyn . pas_de_sortie = 1 ;
  168.  dyne de_vogelaere tabdyn ;
  169. finsi ;
  170.  
  171. *
  172. TABRES = 
  173. DYNE DE_VOGELAERE TABM TABCHAR TABSOR 80 0.025 1 ;
  174. *
  175.  ou1 = prog 0. ; ou11 = prog 0. ;
  176. *d_pasa = 'O';
  177. repeter b_post n2n ;
  178.   DEPU1 = RECO TABRES TABM 
  179.          (0. + (&b_post* 0.1)) 'DEPL'  ;
  180.   chdepu1 = peche tabdyn deplacements
  181.           (0. + (&b_post* 0.1)) ;
  182.   si ('EGA' d_pasa 'O')  ; depu11 = reco mo_moda ca_moda chdepu1 ;
  183.     sinon ; depu11 = chdepu1 ; finsi ;
  184.   FLECHu1 = EXTR DEPu1 P2 UY  ;  
  185.   FLECHu11 = EXTR DEPU11 P2 UY  ; 
  186.    ou1 = ou1 et (prog flechu1) ; 
  187.    ou11 = ou11 et (prog flechu11) ;
  188. fin b_post ;
  189.  ab1 = prog 0. pas 0.1  2. ;
  190.  e1 = evol vert manu ab1 ou1  ;
  191.  e2 = evol rouge manu ab1 ou11  ;
  192.   SI (NEG GRAPH 'N');
  193.  dess (e1 et e2 et ev) ;
  194.   FINSI ;
  195. *
  196. *  RECOMPOSITION POUR LES TEMPS 1.45 ET 1.9 
  197. *
  198. DEP1 = RECO TABRES TABM 1.45 'DEPL'  ;
  199. DEP2 = RECO TABRES TABM 1.9 'DEPL'  ;
  200. *
  201. chdep1 = peche tabdyn deplacements 1.45 ;
  202. si ('EGA' d_pasa 'O')  ; dep11 = reco mo_moda ca_moda chdep1 ;
  203.            sinon ; dep11 = chdep1 ; finsi ;
  204. chdep2 = peche tabdyn deplacements 1.9 ;
  205. si ('EGA' d_pasa 'O')  ; dep22 = reco mo_moda ca_moda chdep2 ;
  206.           sinon ; dep22 = chdep2 ; finsi ;
  207.  
  208. *
  209. *  COMPARAISON AVEC DES RESULTATS ANALYTIQUES 
  210. *
  211.  FLECH1 = EXTR DEP1 P2 UY  ;  FLECH11 = EXTR DEP11 P2 UY  ; 
  212.  
  213.  FLECH1TH = 0.23382   ;
  214.  ERR1 = (( FLECH1 - FLECH1TH) / (ABS FLECH1TH)) * 100.;
  215.  ERR11 = (( FLECH11 - FLECH1TH) / (ABS FLECH1TH)) * 100.;
  216. MESS '    FLECHE MAXIMALE  T = 1.45       '  ;
  217. MESS'THEORIE   'FLECH1TH'CASTEM  'FLECH1'%ERREUR 'ERR1; 
  218. MESS'THEORIE   'FLECH1TH'CASTEM  'FLECH11'%ERREUR 'ERR11; 
  219. *
  220.  FLECH2 = EXTR DEP2 P2 UY  ; FLECH22 = EXTR DEP22 P2 UY  ;
  221. FLECH2TH = -0.24283   ;
  222.  ERR2 = (( FLECH2 - FLECH2TH) / (ABS FLECH2TH)) * 100.;
  223.  ERR22 = (( FLECH22 - FLECH2TH) / (ABS FLECH2TH)) * 100.;
  224. MESS '    FELECHE MAXIMALE  T = 1.9      '  ;
  225. MESS'THEORIE   'FLECH2TH'CASTEM  'FLECH2'%ERREUR 'ERR2; 
  226. MESS'THEORIE   'FLECH2TH'CASTEM  'FLECH22'%ERREUR 'ERR22; 
  227. *
  228. SI ((MAXI ABS (PROG ERR1 ERR2 ERR11 ERR22)) < 5);
  229. ERRE 0;
  230. SINON;
  231. ERRE 5;
  232. FINSI;
  233. FIN ;
  234.  
  235.  
  236.  
  237.  
  238.  
  239.  
  240.  
  241.  
  242.  
  243.  
  244.  
  245.  
  246.  
  247.  
  248.  
  249.  
  250.  
  251.  
  252.  
  253.  

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