Télécharger tristru.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : tristru.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. **********************************************************************
  6. * *
  7. * Mots-clés : Vibrations, calcul modal, modes statiques, poutre *
  8. * *
  9. * test sous-structuration 3 poutres *
  10. * auteur : JK *
  11. * *
  12. **********************************************************************
  13.  
  14. opti dime 3 elem seg2 ;
  15. OPTI EPSI LINEAIRE;
  16.  
  17.  
  18. d_graph = 'F' ;
  19. * premiere poutre bloquee
  20. d_bloq = faux;
  21.  
  22. p0 = 0. 0. 0. ; e_z = 0. 0. 1. ;
  23. q1 = 1. 0. 0. ; q2 = 2. 0. 0. ;
  24. lig1 = d 12 q1 q2 ;
  25.  
  26.  
  27. MO1= MODE LIG1 MECANIQUE ELASTIQUE POUT ;
  28.  
  29. MATPL1=MATE MO1 YOUNG 2.E11 NU 0.3 RHO 7800;
  30. CARPL1= 'CARA' MO1 SECT 0.1 INRY 2.08E-7 INRZ 3.33E-8
  31. TORS 2.41E-7;
  32. MATE1=MATPL1 ET CARPL1;
  33.  
  34. RIGPL1 =RIGI MATE1 MO1 ;
  35. MASPLA1=MASS MATE1 MO1 ;
  36. BLOQ2=BLOQ DEPL ROTA Q2 ;
  37. bloq1r = (bloq rota q1) ; blq1uz = bloq q1 uz ;
  38. blq1ux = bloq q1 ux ; blq1uy = bloq q1 uy ;
  39. *
  40. *---- Recherche des modes propres de la 1�re poutre ---
  41. *
  42. PR=PROG 0.4 ; n_mv = lect 2 ;
  43. MODPLA01=VIBR PROC PR (RIGPL1 ET BLOQ2 et bloq1r et
  44. blq1uz et blq1ux et blq1uy) MASPLA1 ;
  45. modpla0 = vibr proc pr (RIGPL1 ET BLOQ2 et bloq1r et blq1uz) maspla1 ;
  46.  
  47. pl1 = 0. 0. 0. ;
  48. vl1 = manu mult pl1 q1 ;
  49. bl_co1 = manu rigidite type rigidite vl1
  50. (mots lx) (mots ux xy uz) dual (mots flx) (mots fx fy fz)
  51. (prog 0. 1.5 0. -0.8666 0. 0. 0. 0. 0. 0. ) ;
  52.  
  53.  
  54. tblsta = table 'LIAISONS_STATIQUES' ;
  55. tblsta . 1 = table ; tblsta . 1 . point_liaison = q1 ;
  56. si d_bloq ; tblsta . 1 . ddl_liaison = ux ;
  57. sinon ; tblsta . 1 . ddl_liaison = fx ; finsi ;
  58. tblsta . 2 = table ; tblsta . 2 . point_liaison = q1 ;
  59. si d_bloq ; tblsta . 2 . ddl_liaison = uy ;
  60. sinon ; tblsta . 2 . ddl_liaison = fy ; finsi ;
  61.  
  62. si d_bloq ;
  63. bliaiq1 = bloque tblsta ; depi tblsta ;
  64. sinon ;
  65. force tblsta ;
  66. finsi ;
  67. si d_bloq ;
  68. * modes vibrations
  69. MODPLA1=VIBR PROC PR n_mv (RIGPL1 ET BLOQ2 et bloq1r et
  70. blq1uz et bliaiq1) MASPLA1 tbas;
  71. * solutions statiques
  72. tbsta = resout (RIGPL1 ET BLOQ2 et bloq1r et bliaiq1) tblsta ;
  73. react bliaiq1 tbsta ;
  74. sinon ;
  75. * modes vibrations
  76. MODPLA1=VIBR PROC PR n_mv (RIGPL1 ET BLOQ2 et bloq1r et
  77. blq1uz) MASPLA1 tbas;
  78. * solutions statiques
  79. tbsta = resout (RIGPL1 ET BLOQ2 et bloq1r ) tblsta ;
  80. finsi ;
  81.  
  82. tbsta maspla1 * ;
  83.  
  84. repeter bmod ((dime tbsta) - 1) ;
  85. ind1 = &bmod ;
  86. psta1 = manu poi1 tbsta . ind1 . point_repere ;
  87. mostai1 = mode psta1 mecanique elastique statique ;
  88. si d_bloq ;
  89. castai1 = mate mostai1 defo tbsta . ind1 . deformee
  90. made tbsta . ind1 . masse_deformee
  91. ride tbsta . ind1 . reaction ;
  92. sinon ;
  93. castai1 = mate mostai1 defo tbsta . ind1 . deformee
  94. made tbsta . ind1 . masse_deformee
  95. ride tbsta . ind1 . force ;
  96. finsi ;
  97. si (ega 1 ind1) ;
  98. mosta1 = mostai1 ; casta1 = castai1 ;
  99. sinon ;
  100. mosta1 = mosta1 et mostai1 ; casta1 = casta1 et castai1 ;
  101. finsi ;
  102. fin bmod ;
  103. masta1 = masse mosta1 casta1 ;
  104. rista1 = rigi mosta1 casta1 ;
  105.  
  106. modto1 = mosta1 ; casto1 = casta1 ; tmod = modpla1 . modes ;
  107.  
  108. repeter bmod ((dime modpla1 . modes) - 2) ;
  109. ind1 = &bmod ;
  110. psta1 = manu poi1 tmod . ind1 . point_repere ;
  111. momoda1 = mode psta1 mecanique elastique modal ;
  112. camoda1 = mate momoda1 freq tmod . ind1 . frequence
  113. mass tmod . ind1 . masse_generalisee
  114. defo tmod . ind1 . deformee_modale ;
  115. modto1 = modto1 et momoda1 ; casto1 = casto1 et camoda1 ;
  116. fin bmod ;
  117.  
  118. rigto1 = rigi modto1 casto1 ;
  119. masto1 = mass modto1 casto1 ;
  120.  
  121. * solution statique
  122. chq1ux = depi blq1ux 1. ; chq1uy = depi blq1uy 1. ;
  123. stq1ux = resout (RIGPL1 ET BLOQ2 et bloq1r et blq1uz et blq1ux) chq1ux ;
  124. stq1uy = resout (RIGPL1 ET BLOQ2 et bloq1r et blq1uz et blq1uy) chq1uy ;
  125.  
  126.  
  127. * combinaison
  128. evmoq0 = evol chpo (tire modpla0 depl) lig1 uy coul bleu ;
  129. evmoq1 = evol chpo (tire modpla01 depl) lig1 uy ;
  130. evstq1 = evol chpo stq1uy lig1 uy ;
  131. coe_1 = 1. ;
  132. repeter bdes 5 ;
  133. *dess (((coe_1 * evmoq1) + ((1.-coe_1)*evstq1)) et (evmoq0 * (-1.))) ;
  134. coe_1 = coe_1 - 0.25 ;
  135. fin bdes ;
  136. ***
  137.  
  138. modlig1 = vibr proc (prog 0.) (lect 2) rigto1 masto1 tbas ;
  139. dlig1 = reco (modlig1 . modes . 1 . deformee_modale) modpla1 tbsta ;
  140. dlig11 = reco (modlig1 . modes . 1 . deformee_modale) modto1 casto1 ;
  141. deflig1 = defo dlig1 lig1 1. vert ;
  142. deflig11 = defo dlig11 lig1 1. rouge ;
  143. d_lig1 = defo dlig1 lig1 0. ;
  144. si (neg d_graph 'F') ;
  145. trac (0. 0. 1000.) (deflig1 et deflig11 et d_lig1) ;
  146. finsi ;
  147. *
  148. * 2eme poutre
  149. *
  150. q21 q22 lig2 mo2 mate2 = q1 q2 lig1 mo1 mate1 tour 120. p0 e_z ;
  151.  
  152. RIGPL2 =RIGI MATE2 MO2 ;
  153. MASPLA2=MASS MATE2 MO2 ;
  154. BLOQ22 =BLOQ DEPL ROTA Q22 ;
  155. BLOQ21R= BLOQ ROTA Q21 ; blq21uz = bloq q21 uz ;
  156. blq21ux = bloq q21 ux ; blq21uy = bloq q21 uy ;
  157.  
  158. *
  159. PR=PROG 0.4 ;
  160.  
  161. tblsta = table 'LIAISONS_STATIQUES' ;
  162. tblsta . 1 = table ; tblsta . 1 . point_liaison = q21 ;
  163. tblsta . 1 . ddl_liaison = ux ;
  164. tblsta . 2 = table ; tblsta . 2 . point_liaison = q21 ;
  165. tblsta . 2 . ddl_liaison = uy ;
  166.  
  167. bliaiq21 = bloque tblsta ;
  168. depi tblsta ;
  169. * modes vibrations
  170. MODPLA2=VIBR PROC PR n_mv (RIGPL2 ET BLOQ22 et bloq21r et
  171. blq21uz et bliaiq21) MASPLA2 tbas;
  172. * solutions statiques
  173. tbsta2 = resout (RIGPL2 ET BLOQ22 et bloq21r et bliaiq21) tblsta ;
  174. react bliaiq21 tbsta2 ;
  175. tbsta2 maspla2 * ;
  176. repeter bmod ((dime tbsta2) - 1) ;
  177. ind1 = &bmod ;
  178. psta1 = manu poi1 tbsta2 . ind1 . point_repere ;
  179. mostai2 = mode psta1 mecanique elastique statique ;
  180. castai2 = mate mostai2 defo tbsta2 . ind1 . deformee
  181. made tbsta2 . ind1 . masse_deformee
  182. ride tbsta2 . ind1 . reaction ;
  183. si (ega 1 ind1) ;
  184. mosta2 = mostai2 ; casta2 = castai2 ;
  185. sinon ;
  186. mosta2 = mosta2 et mostai2 ; casta2 = casta2 et castai2 ;
  187. finsi ;
  188. fin bmod ;
  189. masta2 = masse mosta2 casta2 ;
  190. rista2 = rigi mosta2 casta2 ;
  191.  
  192. modto2 = mosta2 ; casto2 = casta2 ; tmod = modpla2 . modes ;
  193.  
  194. repeter bmod ((dime modpla2 . modes) - 2) ;
  195. ind1 = &bmod ;
  196. psta1 = manu poi1 tmod . ind1 . point_repere ;
  197. momoda2 = mode psta1 mecanique elastique modal ;
  198. camoda2 = mate momoda2 freq tmod . ind1 . frequence
  199. mass tmod . ind1 . masse_generalisee
  200. defo tmod . ind1 . deformee_modale ;
  201. modto2 = modto2 et momoda2 ; casto2 = casto2 et camoda2 ;
  202. fin bmod ;
  203.  
  204. rigto21 = rigi modto2 casto2 ;
  205. masto21 = mass modto2 casto2 ;
  206.  
  207. t_sta = tbsta et tbsta2 ;
  208. modlig2 = vibr proc (prog 0.) (lect 2) rigto21 masto21 tbas ;
  209. dlig2 = reco (modlig2 . modes . 1 . deformee_modale) modpla2 tbsta2 ;
  210. deflig2 = defo dlig2 lig2 1. vert ;
  211. d_lig2 = defo dlig2 lig2 0. ;
  212. si (neg d_graph 'F') ;
  213. trac (0. 0. 1000.) (deflig2 et d_lig2) ;
  214. finsi ;
  215. *
  216. * 3 eme poutre
  217. *
  218. q31 q32 lig3 mo3 mate3 = q21 q22 lig2 mo2 mate2 tour 120. p0 e_z ;
  219.  
  220. RIGPL3 =RIGI MATE3 MO3 ;
  221. MASPLA3=MASS MATE3 MO3 ;
  222. BLOQ32=BLOQ DEPL ROTA Q32 ;
  223. BLOQ31R = BLOQ ROTA Q31 ; blq31uz = bloq q31 uz ;
  224. blq31ux = bloq q31 ux ; blq31uy = bloq q31 uy ;
  225.  
  226. PR=PROG 0.4 ;
  227. tblsta = table 'LIAISONS_STATIQUES' ;
  228. tblsta . 1 = table ; tblsta . 1 . point_liaison = q31 ;
  229. tblsta . 1 . ddl_liaison = ux ;
  230. tblsta . 2 = table ; tblsta . 2 . point_liaison = q31 ;
  231. tblsta . 2 . ddl_liaison = uy ;
  232.  
  233. bliaiq31 = bloque tblsta ;
  234. depi tblsta ;
  235. * modes vibrations
  236. MODPLA3=VIBR PROC PR n_mv (RIGPL3 ET BLOQ32 et bloq31r et
  237. blq31uz et bliaiq31) MASPLA3 tbas;
  238. * solutions statiques
  239. tbsta3 = resout (RIGPL3 ET BLOQ32 et bloq31r et bliaiq31) tblsta ;
  240. ds3_1 = defo tbsta3 . 1 . deformee lig3 1. rouge ;
  241. ds3_2 = defo tbsta3 . 2 . deformee lig3 1. bleu ;
  242. ds3_0 = defo tbsta3 . 1 . deformee lig3 0. ;
  243. si (neg d_graph 'F') ;
  244. trac (0. 0. 1000.) (ds3_1 et ds3_2 et ds3_0) ;
  245. finsi ;
  246. react bliaiq31 tbsta3 ;
  247. tbsta3 maspla3 * ;
  248.  
  249. repeter bmod ((dime tbsta3) - 1) ;
  250. ind1 = &bmod ;
  251. psta1 = manu poi1 tbsta3 . ind1 . point_repere ;
  252. mostai3 = mode psta1 mecanique elastique statique ;
  253. castai3 = mate mostai3 defo tbsta3 . ind1 . deformee
  254. made tbsta3 . ind1 . masse_deformee
  255. ride tbsta3 . ind1 . reaction ;
  256. si (ega 1 ind1) ;
  257. mosta3 = mostai3 ; casta3 = castai3 ;
  258. sinon ;
  259. mosta3 = mosta3 et mostai3 ; casta3 = casta3 et castai3 ;
  260. finsi ;
  261. fin bmod ;
  262. masta3 = masse mosta3 casta3 ;
  263. rista3 = rigi mosta3 casta3 ;
  264.  
  265. modto3 = mosta3 ; casto3 = casta3 ; tmod = modpla3 . modes ;
  266.  
  267. repeter bmod ((dime modpla3 . modes) - 2) ;
  268. ind1 = &bmod ;
  269. psta1 = manu poi1 tmod . ind1 . point_repere ;
  270. momoda3 = mode psta1 mecanique elastique modal ;
  271. camoda3 = mate momoda3 freq tmod . ind1 . frequence
  272. mass tmod . ind1 . masse_generalisee
  273. defo tmod . ind1 . deformee_modale ;
  274. modto3 = modto3 et momoda3 ; casto3 = casto3 et camoda3 ;
  275. fin bmod ;
  276.  
  277. rigto31 = rigi modto3 casto3 ;
  278. masto31 = mass modto3 casto3 ;
  279.  
  280. modlig3 = vibr proc (prog 0.) (lect 2) rigto31 masto31 tbas ;
  281. dlig3 = reco (modlig3 . modes . 1 .deformee_modale) modpla3 tbsta3 ;
  282. deflig3 = defo dlig3 lig3 1. vert ;
  283. d_lig3 = defo dlig3 lig3 0. ;
  284.  
  285. si (neg d_graph 'F') ;
  286. trac (0. 0. 1000.) (deflig3 et d_lig3) ;
  287. finsi ;
  288. t_sta = t_sta et tbsta3 ;
  289. t_vib = modpla1 et modpla2 et modpla3 ;
  290. mo_to = modto1 et modto2 et modto3 ;
  291. ca_to = casto1 et casto2 et casto3 ;
  292.  
  293. n_te = 1 ;
  294. re_bar1 = ( rela q1 ux + q21 ux + q31 ux )
  295. et ( rela q1 uy + q21 uy + q31 uy ) ;
  296. re_bap1 = pjba re_bar1 mo_to ca_to ;
  297. re_bap0 = pjba re_bar1 t_sta ;
  298. tt1 = idli re_bar1 lig1 ;
  299. tt2 = idli re_bar1 lig2 ;
  300. tt3 = idli re_bar1 lig3 ;
  301. rig_t2 = rigto1 et rigto21 et rigto31 et re_bap1 ;
  302. mas_t2 = masto1 et masto21 et masto31 ;
  303. mod_t2 = vibr proc pr (lect n_te) rig_t2 mas_t2 tbas ;
  304.  
  305. rig_t1 = rigpl1 et bloq2 et rigpl2 et bloq22 et rigpl3 et bloq32
  306. et bloq1r et blq1uz et bloq21r et blq21uz et bloq31r et blq31uz
  307. et re_bar1 ;
  308. mas_t1 = maspla1 et maspla2 et maspla3 ;
  309. mod_t1 = vibr proc pr (lect n_te) rig_t1 mas_t1 tbas ;
  310.  
  311. i_te = 0 ;
  312. repeter b_te n_te ;
  313. i_te = i_te + 1 ;
  314. dp1_t2 = reco (mod_t2 . modes . i_te . deformee_modale)
  315. t_sta t_vib ;
  316. dp1_t2 = reco (mod_t2 . modes . i_te . deformee_modale) mo_to ca_to ;
  317. m_t2 = mod_t2 . modes . i_te . masse_generalisee ** 0.5 ;
  318. m_t2 = -1.*m_t2;
  319. df1_t2 = defo (dp1_t2/m_t2) (lig1 et lig2 et lig3) 10. vert ;
  320.  
  321. dp1_t1 = mod_t1 . modes . i_te . deformee_modale ;
  322. m_t1 = mod_t1 . modes . i_te . masse_generalisee ** 0.5 ;
  323. df1_t1 = defo (dp1_t1 / m_t1) (lig1 et lig2 et lig3) 10. bleu ;
  324. si (neg d_graph 'F') ;
  325. trac (0. 0. 1000.) (df1_t1 et df1_t2) ;
  326. finsi ;
  327. fin b_te ;
  328. err0 = (dp1_t2/m_t2) - (dp1_t1 / m_t1) ;
  329. nerr0= (maxi (abs err0) sans (mots lx)) ;
  330. err2 = (dp1_t2/m_t2) + (dp1_t1 / m_t1) ;
  331. nerr2= (maxi (abs err2) sans (mots lx)) ;
  332. *opti donn 5 ;
  333. si ((nerr0 < 1.e-3) ou (nerr2 < 1.e-3)) ;
  334. erre 0 ; sinon ; erre 5 ; finsi ;
  335. fin ;
  336.  
  337.  
  338.  
  339.  
  340.  
  341.  
  342.  
  343.  
  344.  
  345.  
  346.  
  347.  
  348.  
  349.  
  350.  
  351.  
  352.  
  353.  

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