Télécharger dyne03.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : dyne03.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ****************************************************
  5. *
  6. * EXEMPLE D'UTILISATION DE L'OPERATEUR DYNE
  7. *
  8. * Rupture de tuyauterie avec impact
  9. * Calcul sur base modale avec raideur de choc
  10. * Opérateur Dyne
  11. * D. Combescure
  12. * Aout 2006
  13. *
  14. *****************************************************
  15. * GRAPH = VRAI;
  16. GRAPH = FAUX;
  17. *
  18. opti dime 3 elem seg2;
  19. OPTI EPSI LINEAIRE;
  20. *
  21. L1 = 3.;
  22. n1 = 20;
  23. L = 6.;
  24. n2 = 20;
  25. H = 2.;
  26. nh = 20;
  27. *
  28. p1 = 0. 0. 0.;
  29. p2 = L1 0. 0.;
  30. p2b = L1 0. 0.;
  31. p3 = L 0. 0.;
  32. p3b = L 0. 0.;
  33. lig1 = d n1 p1 p2;
  34. lig2 = d n2 p2 p3;
  35. mesh = lig1 et lig2;
  36. *
  37. P4 = L 0. ((-1.)*H);
  38. *
  39. meshV = d nh P4 P3b;
  40. *
  41. R = 0.30;
  42. P0 = 3.6D6;
  43. sigy = 219.D6;
  44. e = (P0/sigy)*R;
  45. I1 = 0.5*2.*pi*(R**3)*e;
  46. S1 = 2.0*pi*R*e;
  47. S0 = pi* (R**2);
  48. *
  49. R2 = 0.10;
  50. ep2 = 0.01;
  51. I2 = 0.5*2.*pi*(R2**3)*ep2;
  52. S2 = 2.0*pi*R2*ep2;
  53. *
  54. F0 = 1.26*P0*S0;
  55. *
  56. E1 = 200000.D6;
  57. E2= E1;
  58. *
  59. nu1 =0.2;
  60. ro1 = 7.8D3;
  61. *
  62. ro2 = ro1;
  63. nu2 = nu1;
  64. *
  65. *
  66. ustat = ((L**3)/(3.*E1*I1))*F0;
  67. mul = 4.*2.*sigy*e*(R**2);
  68. *
  69. mod1 = 'MODELI' mesh mecanique elastique POUT;
  70. mat1 = mate mod1 YOUN E1 NU nu1 RHO ro1
  71. INRY I1 INRZ I1 TORS (2.*I1)
  72. SECT S1;
  73. K1 = (3.*E1*I1/((L**3)));
  74. *
  75. mod2 = 'MODELI' meshV mecanique elastique POUT;
  76. mat2 = mate mod2 YOUN E2 NU nu2 RHO ro2
  77. INRY I2 INRZ I2 TORS (2.*I2)
  78. SECT S2;
  79. K2 = (3.*E2*I2/((H**3)));
  80. *
  81. modT = mod1 et mod2;
  82. Matt = mat1 et mat2;
  83. MeshT = mesh et MeshV;
  84. *
  85. rig1 = RIGI modT matT;
  86. mas1 = mass modT matT;
  87. bl1 = (BLOQ p1 RY)
  88. et (BLOQ p1 UY) et (BLOQ p1 RZ);
  89. bl2b = BLOQ p2b DEPL;
  90. bl3b = (BLOQ P4 DEPL ROTA) et (BLOQ p3b UX);
  91. blM = (BLOQ mesh UZ) et (BLOQ mesh RX)
  92. et (BLOQ mesh UX);
  93. *
  94. For1 = FORCE (0. (F0) 0.) p3;
  95. *
  96. *
  97. * Caractéristiques liaisons
  98. *
  99. umax = 1.17;
  100. Kliai = 15.6D6;
  101. Jeuliai = 0.500;
  102. ******************************
  103. ******************************
  104. *
  105. * Cas 1: 2 modes avec choc mou
  106. *
  107. ******************************
  108. ******************************
  109. nmode = 2;
  110. FR = VIBR PROCHE (PROG 11. ) (LECT nmode)
  111. (RIG1 et bl1 et blm et bl3b) MAS1;
  112. *
  113. Mequi3 = TABLE;
  114. *
  115. *
  116. lmot2 = MOTS UX UY UZ RX RY RZ;
  117. lmot1 = MOTS FX FY FZ MX MY MZ;
  118. lmot3 = lmot2;
  119. IdFYT = (manu chpo meshT 6 UX 0. UY (-1.0) UZ 0. RX 0. RY 0. RZ 0.
  120. NATURE DISCRET);
  121. *
  122. oeil = 10. 10. 2.;
  123. Repeter bou1 nmode;
  124. FIFI = FR . MODES . &Bou1 . FREQUENCE ;
  125. TT = FR . MODES . &Bou1 . DEFORMEE_MODALE ;
  126. DEf1 = DEFO meshT TT;
  127. SI GRAPH;
  128. titre (CHAIN 'Frequence=' Fifi);
  129. trac oeil def1;
  130. FINSI;
  131. Mm1 = resu ((mas1*(TT)) '*' (TT) lmot1 lmot2 lmot3);
  132. Mm1V = ( ((maxi (exco Mm1 UX))) +
  133. ((maxi (exco Mm1 UY))) +
  134. ((maxi (exco Mm1 UZ))) +
  135. ((maxi (exco Mm1 RX))) +
  136. ((maxi (exco Mm1 RY))) +
  137. ((maxi (exco Mm1 RZ))) );
  138. Mequi3.&bou1 = MM1V/((extr TT P3 UY)**2);
  139. fin bou1;
  140. *
  141. Cliai = 2.*((Kliai*(Mequi3.1))**0.5);
  142. *
  143. * Appel à DYNE
  144. *
  145. ************************************************
  146. * TEMPS
  147. *************************
  148. *
  149. DT = MINI (PROG (1.E-5) (0.1/(2.*pi*fifi))) ;
  150. NTT = ENTIER (0.1/DT) ;
  151. *
  152. *************************
  153. * TABLE AMORTISSEMENT
  154. *************************
  155. *
  156. TAMOR = TABLE 'AMORTISSEMENT' ;
  157. ValAmor = 0.;
  158. AMO = PROG NMODE*ValAmor ;
  159. TAMOR . AMORTISSEMENT = AMOR FR AMO ;
  160. *
  161. *************************
  162. * TABLE CHARGEMENT (CHARG1 CHARG2 CHARG3)
  163. *************************
  164. TCHAR = TABLE 'CHARGEMENT' ;
  165. *
  166. Ampl = 1.0;
  167. Evotime = evol manu (prog 0. 0.01 100.)
  168. (Ampl*(prog 0. 1. 1.));
  169. CHA1 = CHAR For1 Evotime;
  170. CHAP1 = PJBA FR CHA1;
  171. TCHAR.'BASE_A' = CHAP1;
  172. ***********************************************************
  173. *
  174. * LIAISONS POINT-POINT
  175. *
  176. ***********************************************************
  177. *
  178. TILIA = TABLE 'LIAISON' ;
  179. TTLB = TABLE 'LIAISON_B' ;
  180. TILIA.'LIAISON_B' = TTLB ;
  181. *
  182. TLA12 = TABLE 'LIAISON_ELEMENTAIRE' ;
  183. TTLB . 1 = TLA12 ;
  184. *
  185. TLA12 . 'TYPE_LIAISON' = MOT 'POINT_POINT';
  186. TLA12 . 'POINT_A' = P3 ;
  187. TLA12 . 'POINT_B' = P3b ;
  188. TLA12. 'NORMALE' = ( 0. 1. 0.);
  189. TLA12 . 'RAIDEUR' = Kliai ;
  190. TLA12 . 'JEU' = JeuLiai ;
  191. TLA12 . 'AMORTISSEMENT' = Cliai;
  192. *************************
  193. * TABLE SORTIE
  194. *************************
  195. *
  196. TABSOR = TABLE 'SORTIE';
  197. *
  198. TABSOR2 = TABLE 'LIAISON_B' ;
  199. TABSOR.'LIAISON_B' = TABSOR2 ;
  200. *
  201. TVARPP = TABLE 'VARIABLE' ;
  202. TVARPP.'FORCE_DE_CHOC_POINT_A' = VRAI ;
  203. TVARPP.'FORCE_DE_CHOC_POINT_B' = VRAI ;
  204. TVARPP.'UY_POINT_A' = VRAI ;
  205. TVARPP.'UY_POINT_B' = VRAI ;
  206. *
  207. TVARGD = TABLE 'VARIABLE' ;
  208. TVARGD.'FORCE_DE_CHOC' = VRAI ;
  209. TVARGD.'UY' = VRAI ;
  210. *
  211. TABSOR2. TLA12 = TVARPP;
  212. *
  213. ************************************************
  214. * CALCUL TEMPOREL PROPREMENT DIT
  215. *************************
  216. *
  217. TRESU = DYNE DE_VOGELAERE FR TAMOR TCHAR TILIA
  218. NTT DT 1 TABSOR ;
  219. *
  220. ************************************************
  221. * TRACE DES DEPLACEMNTS EN FONCTION DU TEMPS
  222. *************************
  223. *
  224. UXP2 = EVOL VERT RECO TRESU FR 'DEPL' P2 UY ;
  225. UXP3 = EVOL BLEU RECO TRESU FR 'DEPL' P3 UY ;
  226. UXPB3 = EVOL BLEU RECO TRESU FR 'DEPL' P3B UY ;
  227. UXPB3 = evol manu (extr UXPB3 absc)
  228. ((extr UXPB3 ordo) + (prog (dime (extr UXPB3 ordo))*jeuliai));
  229. *
  230. TDS = TRESU . 'TEMPS_DE_SORTIE' ;
  231. RESCHOC = TRESU . TLA12;
  232. FchocB = RESCHOC. FORCE_DE_CHOC_POINT_B;
  233. FchocA = RESCHOC. FORCE_DE_CHOC_POINT_A;
  234. *
  235. EVFchoA = 'EVOL' 'MANU' TDS FChocA;
  236. EVFchoB = 'EVOL' 'MANU' TDS FChocB;
  237. *
  238. *
  239. Ncont = 10.;
  240. REPETER bo3 (ENTIER ((NTT + 1)/Ncont));
  241. dep1 = RECO TRESU FR ((&bo3 - 1 )*DT*Ncont) 'DEPL';
  242. con1 = SIGMA MOD1 MAT1 dep1;
  243. con2 = SIGMA MOD2 MAT2 dep1;
  244. SI (&bo3 EGA 1);
  245. time = PROG ((&bo3 - 1 )*DT*Ncont);
  246. pfyp1 = prog (extr con1 1 1 2 EFFY);
  247. pmzp1 = prog (extr con1 1 1 2 MOMZ);
  248. deft = defo meshT dep1 1.;
  249. SINON;
  250. time = time et (PROG ((&bo3 - 1 )*DT*Ncont));
  251. pfyp1 = pfyp1 et (prog (extr con1 1 1 2 EFFY));
  252. pmzp1 = pmzp1 et (prog (extr con1 1 1 2 MOMZ));
  253. deft = deft et (defo meshT dep1 1.);
  254. FINSI;
  255. FIN bo3;
  256. tifyp1 = evol manu Time pfyp1;
  257. timzp1 = evol manu Time pmzp1;
  258. *
  259. SI GRAPH;
  260. TRAC oeil DEFT ANIME;
  261. TITRE 'Deplacement avec 20 modes';
  262. dess (UXP2 et UXP3 et UXPB3) xbord 0. 0.10 ybord 0. 1.40;
  263. Titre 'Force de choc';
  264. DESS EVFchoA xbord 0. 0.10;
  265. FINSI;
  266. *
  267. *
  268. FIN;
  269. *
  270.  
  271.  
  272.  
  273.  
  274.  
  275.  
  276.  
  277.  
  278.  
  279.  

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