Télécharger dyna_nl1.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dyna_nl1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test Dyna_nl1.dgibi: Jeux de données *
  5. * ------------------------------------ *
  6. * *
  7. complet = vrai;
  8. *
  9. *************************************************
  10. * pour calcul complet mettre complet à : vrai;
  11. *************************************************
  12. *******************************************************
  13. *
  14. * dynamique non lineaire geometrique
  15. *
  16. * oscillations libres d'un oscillateur de type Duffing
  17. *
  18. * reference :J. ARGYRIS H.P. MLEJNEK
  19. * DYNAMICS OF STRUCTURES (p. 494)
  20. *
  21. *******************************************************
  22. opti echo 1 ;
  23. graph = faux ;
  24.  
  25.  
  26. OPTI DIME 2 ELEM SEG2 MODE PLAN CONT ;
  27.  
  28. P1 = 0. 0. ; P2 = 1. 0. ;
  29. ST = P1 D 1 P2 ;
  30.  
  31. MO = MODE ST MECANIQUE BARR ;
  32. MA1 = MATE MO YOUN 1E7 NU 0.3 RHO 0.;
  33. CAR1 = CARA MO SECT 1. ;
  34. MACA= MA1 ET CAR1;
  35. MASS1 = MASS UY 250. P2 ;
  36.  
  37. CL1 = BLOQ UX uy P1 ;
  38. CL2 = BLOQ UX P2 ;
  39. CL3 = BLOQ UY P2 ;
  40. CL = CL1 ET CL2 ;
  41.  
  42.  
  43. * precontrainte de 500 N
  44.  
  45. RI1 = RIGI MACA MO ;
  46. f1 = force fx 500 p2 ;
  47. dep1 = reso (ri1 et cl1 et cl3) f1 ;
  48. sig1 = sigm dep1 mo maca ;
  49. *sig1 = manu 'CHML' mo 'EFFX' 500. 'TYPE' 'CONTRAIN'
  50. * 'STRESSES' ;
  51.  
  52. * deplacement impose en statique
  53.  
  54. t1 = 1e-10 ;
  55. FP11 = DEPI CL3 .2;
  56. LIX1 = PROG 0. t1 100 ;
  57. LIY1 = PROG 0. 1 1 ;
  58. EV1 = EVOL MANU T LIX1 F(T) LIY1 ;
  59. CHA1 = CHAR MECA FP11 EV1 ;
  60. tab2 = table ;
  61. tab2.'GRANDS_DEPLACEMENTS'=VRAI;
  62. TAB2.'BLOCAGES_MECANIQUES' = CL et cl3 ;
  63. TAB2.'MODELE' = MO;
  64. Tab2.'K_SIGMA' = VRAI ;
  65. TAB2.'CARACTERISTIQUES' = MACA;
  66. TAB2.'CHARGEMENT' = CHA1;
  67. TAB2.'CONTRAINTES' = TABLE;
  68. TAB2.'CONTRAINTES' . 0 = SIG1 ;
  69. LIS11 = PROG 0. (.33*t1) (.66*t1) t1 ;
  70. TAB2.'TEMPS_CALCULES' = LIS11;
  71. PASAPAS TAB2 ;
  72.  
  73. * lacher et oscillations libres en dynamique
  74. si complet;
  75. t2 = .27 ;
  76. sinon;
  77. t2 = .02;
  78. finsi;
  79. dt = .002 ;
  80. f1 = force fy 0. p2 ;
  81. cha2 = char meca f1 ev1 ;
  82. tab2.dynamique = vrai ;
  83. tab2.masse_constante = mass1 ;
  84. tab2.ftol = 1e-4 ;
  85. vv = manu chpo ST 1 'UX' 0.;
  86. tab2.'CONTINUATION'.'VITESSES'= vv;
  87. nso= dime (tab2.'TEMPS');tab2.'VITESSES'= table;
  88. tab2.'VITESSES'.(nso + 1) = vv;
  89. tab2.'CONTINUATION'.'ACCELERATIONS'= vv;
  90. tab2.'ACCELERATIONS'= table;
  91. tab2.'ACCELERATIONS'.(nso + 1) = vv;
  92. TAB2.'BLOCAGES_MECANIQUES' = CL ;
  93. LIS22 = PROG (t1 + dt) pas dt (t1 + t2) ;
  94. TAB2.'TEMPS_CALCULES' = LIS22;
  95. TAB2.'CHARGEMENT' = CHA2;
  96. PASAPAS TAB2 ;
  97.  
  98.  
  99. NDIM = (DIME ( TAB2 . DEPLACEMENTS )) ;
  100. PGY = PROG ;
  101. REPETER TBOU2 NDIM;
  102. LEDEP = TAB2 . DEPLACEMENTS. (&TBOU2 - 1);
  103. V = EXTR LEDEP UY P2 ;
  104. PGY = PGY ET ( PROG V ) ;
  105. FIN TBOU2 ;
  106.  
  107. ev1 = evol manu 'temps' (lis11 et lis22) 'depl' pgy ;
  108. si graph ;
  109. dess ev1 ;
  110. finsi ;
  111.  
  112. *recherche de l'amplitude et de la periode
  113. si complet;
  114. dim1 = enti ((t2 - .25)/dt) ;
  115. maxuy = 0 ;
  116. tt = .25 ;
  117. tp=0.;
  118. repeter bou1 dim1 ;
  119. tt = tt + dt ;
  120. uy1 = ipol tt (lis11 et lis22) pgy ;
  121. si (uy1 > maxuy) ;
  122. maxuy = uy1 ;
  123. tp = tt ;
  124. finsi ;
  125. fin bou1 ;
  126. valres = .2;
  127. sinon;
  128. tp =.2648;
  129. valres = .171 ;maxuy = pgy extr ( dime pgy);
  130. finsi;
  131.  
  132.  
  133. ** periode theorique .264791269s
  134. ** amplitude theorique .2 ;
  135. err1 = abs ((maxuy - valres)/valres * 100) ;
  136. err2 = abs ((tp - .2648)/.2648 *100) ;
  137. mess 'erreur en amplitude'err1'%';
  138. mess 'erreur en periode'err2'%' ;
  139. si ((err1 > .1) ou (err2 > .5));
  140. erre 5 ;
  141. sinon ;
  142. erre 0 ;
  143. finsi ;
  144. fin ;
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  
  177.  
  178.  
  179.  

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