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)
  112. MAS1 MULT TBAS;
  113. *
  114. Mequi3 = TABLE;
  115. *
  116. *
  117. lmot2 = MOTS UX UY UZ RX RY RZ;
  118. lmot1 = MOTS FX FY FZ MX MY MZ;
  119. lmot3 = lmot2;
  120. IdFYT = (manu chpo meshT 6 UX 0. UY (-1.0) UZ 0. RX 0. RY 0. RZ 0.
  121. NATURE DISCRET);
  122. *
  123. oeil = 10. 10. 2.;
  124. Repeter bou1 nmode;
  125. FIFI = FR . MODES . &Bou1 . FREQUENCE ;
  126. TT = FR . MODES . &Bou1 . DEFORMEE_MODALE ;
  127. DEf1 = DEFO meshT TT;
  128. SI GRAPH;
  129. titre (CHAIN 'Frequence=' Fifi);
  130. trac oeil def1;
  131. FINSI;
  132. Mm1 = resu ((mas1*(TT)) '*' (TT) lmot1 lmot2 lmot3);
  133. Mm1V = ( ((maxi (exco Mm1 UX))) +
  134. ((maxi (exco Mm1 UY))) +
  135. ((maxi (exco Mm1 UZ))) +
  136. ((maxi (exco Mm1 RX))) +
  137. ((maxi (exco Mm1 RY))) +
  138. ((maxi (exco Mm1 RZ))) );
  139. Mequi3.&bou1 = MM1V/((extr TT P3 UY)**2);
  140. fin bou1;
  141. *
  142. Cliai = 2.*((Kliai*(Mequi3.1))**0.5);
  143. *
  144. * Appel à DYNE
  145. *
  146. ************************************************
  147. * TEMPS
  148. *************************
  149. *
  150. DT = MINI (PROG (1.E-5) (0.1/(2.*pi*fifi))) ;
  151. NTT = ENTIER (0.1/DT) ;
  152. *
  153. *************************
  154. * TABLE AMORTISSEMENT
  155. *************************
  156. *
  157. TAMOR = TABLE 'AMORTISSEMENT' ;
  158. ValAmor = 0.;
  159. AMO = PROG NMODE*ValAmor ;
  160. TAMOR . AMORTISSEMENT = AMOR FR AMO ;
  161. *
  162. *************************
  163. * TABLE CHARGEMENT (CHARG1 CHARG2 CHARG3)
  164. *************************
  165. TCHAR = TABLE 'CHARGEMENT' ;
  166. *
  167. Ampl = 1.0;
  168. Evotime = evol manu (prog 0. 0.01 100.)
  169. (Ampl*(prog 0. 1. 1.));
  170. CHA1 = CHAR For1 Evotime;
  171. CHAP1 = PJBA FR CHA1;
  172. TCHAR.'BASE_A' = CHAP1;
  173. ***********************************************************
  174. *
  175. * LIAISONS POINT-POINT
  176. *
  177. ***********************************************************
  178. *
  179. TILIA = TABLE 'LIAISON' ;
  180. TTLB = TABLE 'LIAISON_B' ;
  181. TILIA.'LIAISON_B' = TTLB ;
  182. *
  183. TLA12 = TABLE 'LIAISON_ELEMENTAIRE' ;
  184. TTLB . 1 = TLA12 ;
  185. *
  186. TLA12 . 'TYPE_LIAISON' = MOT 'POINT_POINT';
  187. TLA12 . 'POINT_A' = P3 ;
  188. TLA12 . 'POINT_B' = P3b ;
  189. TLA12. 'NORMALE' = ( 0. 1. 0.);
  190. TLA12 . 'RAIDEUR' = Kliai ;
  191. TLA12 . 'JEU' = JeuLiai ;
  192. TLA12 . 'AMORTISSEMENT' = Cliai;
  193. *************************
  194. * TABLE SORTIE
  195. *************************
  196. *
  197. TABSOR = TABLE 'SORTIE';
  198. *
  199. TABSOR2 = TABLE 'LIAISON_B' ;
  200. TABSOR.'LIAISON_B' = TABSOR2 ;
  201. *
  202. TVARPP = TABLE 'VARIABLE' ;
  203. TVARPP.'FORCE_DE_CHOC_POINT_A' = VRAI ;
  204. TVARPP.'FORCE_DE_CHOC_POINT_B' = VRAI ;
  205. TVARPP.'UY_POINT_A' = VRAI ;
  206. TVARPP.'UY_POINT_B' = VRAI ;
  207. *
  208. TVARGD = TABLE 'VARIABLE' ;
  209. TVARGD.'FORCE_DE_CHOC' = VRAI ;
  210. TVARGD.'UY' = VRAI ;
  211. *
  212. TABSOR2. TLA12 = TVARPP;
  213. *
  214. ************************************************
  215. * CALCUL TEMPOREL PROPREMENT DIT
  216. *************************
  217. *
  218. TRESU = DYNE DE_VOGELAERE FR TAMOR TCHAR TILIA
  219. NTT DT 1 TABSOR ;
  220. *
  221. ************************************************
  222. * TRACE DES DEPLACEMNTS EN FONCTION DU TEMPS
  223. *************************
  224. *
  225. UXP2 = EVOL VERT RECO TRESU FR 'DEPL' P2 UY ;
  226. UXP3 = EVOL BLEU RECO TRESU FR 'DEPL' P3 UY ;
  227. UXPB3 = EVOL BLEU RECO TRESU FR 'DEPL' P3B UY ;
  228. UXPB3 = evol manu (extr UXPB3 absc)
  229. ((extr UXPB3 ordo) + (prog (dime (extr UXPB3 ordo))*jeuliai));
  230. *
  231. TDS = TRESU . 'TEMPS_DE_SORTIE' ;
  232. RESCHOC = TRESU . TLA12;
  233. FchocB = RESCHOC. FORCE_DE_CHOC_POINT_B;
  234. FchocA = RESCHOC. FORCE_DE_CHOC_POINT_A;
  235. *
  236. EVFchoA = 'EVOL' 'MANU' TDS FChocA;
  237. EVFchoB = 'EVOL' 'MANU' TDS FChocB;
  238. *
  239. *
  240. Ncont = 10.;
  241. REPETER bo3 (ENTIER ((NTT + 1)/Ncont));
  242. dep1 = RECO TRESU FR ((&bo3 - 1 )*DT*Ncont) 'DEPL';
  243. con1 = SIGMA MOD1 MAT1 dep1;
  244. con2 = SIGMA MOD2 MAT2 dep1;
  245. SI (&bo3 EGA 1);
  246. time = PROG ((&bo3 - 1 )*DT*Ncont);
  247. pfyp1 = prog (extr con1 1 1 2 EFFY);
  248. pmzp1 = prog (extr con1 1 1 2 MOMZ);
  249. deft = defo meshT dep1 1.;
  250. SINON;
  251. time = time et (PROG ((&bo3 - 1 )*DT*Ncont));
  252. pfyp1 = pfyp1 et (prog (extr con1 1 1 2 EFFY));
  253. pmzp1 = pmzp1 et (prog (extr con1 1 1 2 MOMZ));
  254. deft = deft et (defo meshT dep1 1.);
  255. FINSI;
  256. FIN bo3;
  257. tifyp1 = evol manu Time pfyp1;
  258. timzp1 = evol manu Time pmzp1;
  259. *
  260. SI GRAPH;
  261. TRAC oeil DEFT ANIME;
  262. TITRE 'Deplacement avec 20 modes';
  263. dess (UXP2 et UXP3 et UXPB3) xbord 0. 0.10 ybord 0. 1.40;
  264. Titre 'Force de choc';
  265. DESS EVFchoA xbord 0. 0.10;
  266. FINSI;
  267. *
  268. *
  269. FIN;
  270. *
  271.  
  272.  
  273.  
  274.  
  275.  
  276.  
  277.  
  278.  

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