Télécharger dyne02.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dyne02.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ******************************
  5. *
  6. * Impact d'une barre sur un plan rigide
  7. * Calcul sur base modale
  8. * Opérateur Dyne
  9. * D. Combescure
  10. * Aout 2006
  11. *
  12. ********************************
  13. *
  14. GRAPH = VRAI; opti trac PSC EPTR 5 POTR HELVETICA_16;
  15. *GRAPH = FAUX;
  16. *
  17. opti dime 3 elem seg2;
  18. OPTI EPSI LINEAIRE;
  19. *
  20. L1 = 1.;
  21. n1 = 10;
  22. L = 2.;
  23. n2 = 10;
  24. *
  25. p1 = 0. 0. 0.;
  26. p2 = L1 0. 0.;
  27. p2b = L1 0. 0.;
  28. p3 = L 0. 0.;
  29. lig1 = d n1 p1 p2;
  30. lig2 = d n2 p2b p3;
  31. mesh = lig1;
  32. *
  33. S1 = 0.1;
  34. *
  35. E1 = 2000.;
  36. nu1 =0.2;
  37. ro1 = 7.8D-3;
  38. *
  39. mod1 = 'MODE' mesh mecanique elastique BARR;
  40. mat1 = mate mod1 YOUN E1 NU nu1 RHO ro1
  41. SECT S1;
  42. *
  43. rig1 = RIGI mod1 mat1;
  44. mas1 = mass mod1 mat1;
  45. bl1 = (BLOQ lig1 UY) et (bloq lig1 UZ);
  46. bl2 = (BLOQ lig2 UY) et (bloq lig2 UZ);
  47. bl2b = (BLOQ p2b DEPL);
  48. *
  49. For1 = MAS1 * (MANU CHPO mesh 1 UX (-9.81));
  50. *
  51. nmode = 2;
  52. FR = VIBR PROCHE (PROG nmode) (LECT nmode)
  53. (RIG1 et bl1 et bl2b) MAS1 ;
  54. *
  55. SI GRAPH;
  56. oeil = 0. 0. 100.;
  57. Repeter bou1 nmode;
  58. FIFI = FR . MODES . &Bou1 . FREQUENCE ;
  59. TT = FR . MODES . &Bou1 . DEFORMEE_MODALE ;
  60. DEf1 = DEFO mesh TT;
  61. titre (CHAIN 'Frequence=' Fifi);
  62. trac oeil def1;
  63. fin bou1;
  64. FINSI;
  65. *
  66. *
  67. * Appel à DYNE
  68. *
  69. ************************************************
  70. * TEMPS
  71. *************************
  72. *
  73. DT = 1.E-4 ;
  74. NTT = 5000 ;
  75. *
  76. *************************
  77. * TABLE AMORTISSEMENT
  78. *************************
  79. *
  80. TAMOR = TABLE 'AMORTISSEMENT' ;
  81. AMO = PROG NMODE*0.0 ;
  82. TAMOR . AMORTISSEMENT = AMOR FR AMO ;
  83. *
  84. *************************
  85. * TABLE CHARGEMENT (CHARG1 CHARG2 CHARG3)
  86. *************************
  87. TCHAR = TABLE 'CHARGEMENT' ;
  88. *
  89. Ampl = 1.0;
  90. Evotime = evol manu (prog 0. 0.01 100.)
  91. (Ampl*(prog 0. 1. 1.));
  92. CHA1 = CHAR For1 Evotime;
  93. CHAP1 = PJBA FR CHA1;
  94. TCHAR.'BASE_A' = CHAP1;
  95. ***********************************************************
  96. *
  97. * LIAISONS
  98. *
  99. ***********************************************************
  100. *
  101. TILIA = TABLE 'LIAISON' ;
  102. TTLB = TABLE 'LIAISON_B' ;
  103. TILIA.'LIAISON_B' = TTLB ;
  104. ***********************************************************
  105. * POINT POINT
  106. ***********************************************************
  107. *
  108. Kliai = 30.0*E1*S1/L1;
  109. Jeuliai = .01;
  110. TLA12 = TABLE 'LIAISON_ELEMENTAIRE' ;
  111. TTLB . 1 = TLA12 ;
  112. TLA12 . 'TYPE_LIAISON' = MOT 'POINT_POINT';
  113. TLA12 . 'POINT_A' = P2 ;
  114. TLA12 . 'POINT_B' = P2b ;
  115. TLA12. 'NORMALE' = ( -1. 0. 0.);
  116. TLA12 . 'RAIDEUR' = Kliai ;
  117. TLA12 . 'JEU' = JeuLiai ;
  118. TLA12 . 'AMORTISSEMENT' = 0. ;
  119. *************************
  120. * TABLE SORTIE
  121. *************************
  122. *
  123. TABSOR = TABLE 'SORTIE';
  124. *
  125. TABSOR2 = TABLE 'LIAISON_B' ;
  126. TABSOR.'LIAISON_B' = TABSOR2 ;
  127. *
  128. TVARPP = TABLE 'VARIABLE' ;
  129. TVARPP.'FORCE_DE_CHOC_POINT_A' = VRAI ;
  130. TVARPP.'FORCE_DE_CHOC_POINT_B' = VRAI ;
  131. TVARPP.'UY_POINT_A' = VRAI ;
  132. TVARPP.'UY_POINT_B' = VRAI ;
  133. TVARPP.'VITESSE' = VRAI ;
  134. *
  135. TVARGD = TABLE 'VARIABLE' ;
  136. TVARGD.'FORCE_DE_CHOC' = VRAI ;
  137. TVARGD.'VITESSE' = VRAI ;
  138. TVARGD.'UY' = VRAI ;
  139. *
  140. TABSOR2. TLA12 = TVARPP;
  141. *
  142. ************************************************
  143. * CALCUL TEMPOREL PROPREMENT DIT
  144. ***********************************************
  145. *
  146. *TRESU = DYNE DE_VOGELAERE FR TAMOR TCHAR TILIA
  147. * NTT DT 1 TABSOR ;
  148. TRESU = DYNE DIFFERENCES_CENTREES FR TAMOR TCHAR TILIA
  149. NTT DT 1 TABSOR ;
  150. *
  151. ************************************************
  152. * TRACE DES DEPLACEMENTS EN FONCTION DU TEMPS
  153. **********************************************
  154. *
  155. UXP2 = EVOL VERT RECO TRESU FR 'DEPL' P2 UX ;
  156. UXP3 = EVOL BLEU RECO TRESU FR 'DEPL' P3 UX ;
  157. UXP3 = evol bleu manu (extr UXP3 absc)
  158. ((extr UXP3 ordo) + (PROG (dime (extr UXP3 ordo) )*(-0.01)));
  159. VXP2 = EVOL VERT RECO TRESU FR 'VITE' P2 UX ;
  160. VXP3 = EVOL BLEU RECO TRESU FR 'VITE' P3 UX ;
  161. *
  162. *
  163. *
  164. TDS = TRESU . 'TEMPS_DE_SORTIE' ;
  165. RESCHOC = TRESU . TLA12;
  166. FchocB = RESCHOC. FORCE_DE_CHOC_POINT_B;
  167. FchocA = RESCHOC. FORCE_DE_CHOC_POINT_A;
  168. *
  169. EVFchoA = 'EVOL' 'MANU' TDS FChocA;
  170. EVFchoB = 'EVOL' 'MANU' TDS FChocB;
  171. *
  172. SI GRAPH;
  173. ta = table;
  174. ta.2= 'TIRR';
  175. DESS (UXP2 ET UXP3) ta XBORD 0. 0.20 TITRE 'DEPL BASE MODALE' ;
  176. DESS (VXP2 ET VXP3) ta XBORD 0. 0.20 TITRE 'VITE BASE MODALE' ;
  177. DESS EVFchoA XBORD 0. 0.20 ;
  178. DESS EVFchoB XBORD 0. 0.20 ;
  179. FINSI;
  180. *
  181.  
  182. * test du fonctionnement de RECO
  183. tt = 0.1;
  184. * list 'RESUME' TRESU;
  185. uphys = RECO TRESU FR 'DEPL' tt;
  186.  
  187. ***********************************************
  188. * TEST DE NON REGRESSION
  189. * par rapport a une execution du cas test
  190. ***********************************************
  191. uref = -6.71056E-03;
  192. ureco = EXTR uphys P2 'UX';
  193. uevol = IPOL UXP2 tt;
  194. MESS uref ureco uevol;
  195. ERR1 = ABS ((ureco - uref) / uref);
  196. ERR2 = ABS ((uevol - uref) / uref);
  197. ERR3 = ABS ((ureco - uevol) / uref);
  198. MESS ERR1 ERR2 ERR3;
  199.  
  200. SI ((ERR1 > 0.005) OU (ERR2 > 0.005) OU (ERR3 > 0.0001));
  201. ERRE 5;
  202. FINSI;
  203.  
  204. FIN ;
  205.  
  206.  
  207.  
  208.  
  209.  
  210.  
  211.  
  212.  
  213.  
  214.  
  215.  
  216.  
  217.  

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