Télécharger spectral.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : spectral.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. *
  6. **********************************************************************
  7. * *
  8. * Mots-clés : Vibrations, calcul modal, modes statiques, poutre *
  9. * reponse frequentielle *
  10. * *
  11. * calcul spectral sans/avec amortissement, base modale seule. *
  12. * auteur : JK / a parfaire *
  13. * *
  14. **********************************************************************
  15. *
  16. opti dime 3 elem seg2 ;
  17. OPTI EPSI LINEAIRE;
  18. *
  19. d_bloq = vrai;
  20. d_graph = faux ;
  21.  
  22. p0 = 0. 0. 0. ; e_z = 0. 0. 1. ;
  23. q2 = e_z ; q1 = p0 ;
  24. lig1 = d 20 q2 q1 ;
  25. qf1 = point lig1 proc (0. 0. 0.6) ;
  26.  
  27. MO1= MODE LIG1 MECANIQUE ELASTIQUE POUT ;
  28.  
  29. MATPL1=MATE MO1 YOUNG 2.E10 NU 0.4 RHO 7800;
  30. CARPL1= 'CARA' MO1 SECT 0.1 INRY 2.08E-8 INRZ 3.33E-8
  31. TORS 2.41E-7;
  32. MATE1=MATPL1 ET CARPL1;
  33. ***
  34.  
  35.  
  36. ***
  37.  
  38. RIGPL1 =RIGI MATE1 MO1 ;
  39. MASPLA1=MASS MATE1 MO1 ;
  40. BLOQ2=BLOQ DEPL ROTA Q2 ;
  41. bloq1r = (bloq rota q1) ; blq1uz = bloq q1 uz ;
  42. bloq1rz = bloq q1 rz ;
  43. blq1ux = bloq q1 ux ; blq1uy = bloq q1 uy ;
  44. bloq1rx = bloq q1 rx ; bloq1ry= bloq q1 ry ;
  45. *
  46. *----
  47. *
  48. PR=PROG 0.01 ; n_mv = lect 4 ;
  49.  
  50. tblsta = table 'LIAISONS_STATIQUES' ;
  51. tblsta . 1 = table ; tblsta . 1 . point_liaison = q1 ;
  52. si d_bloq ; tblsta . 1 . ddl_liaison = ux ;
  53. sinon ; tblsta . 1 . ddl_liaison = fx ; finsi ;
  54. tblsta . 2 = table ; tblsta . 2 . point_liaison = q1 ;
  55. si d_bloq ; tblsta . 2 . ddl_liaison = uy ;
  56. sinon ; tblsta . 2 . ddl_liaison = fy ; finsi ;
  57. tblsta . 3 = table ; tblsta . 3 . point_liaison = q1 ;
  58. si d_bloq ; tblsta . 3 . ddl_liaison = rx ;
  59. sinon ; tblsta . 3 . ddl_liaison = mx ; finsi ;
  60.  
  61. si d_bloq ;
  62. bliaiq1 = bloque tblsta ; depi tblsta ;
  63. sinon ;
  64. force tblsta ;
  65. finsi ;
  66. si d_bloq ;
  67. * modes vibrations
  68. bl1bloq = bloq1ry et bloq1rz et blq1uz et bliaiq1 ;
  69. MODPLA1=VIBR PROC PR n_mv (RIGPL1 et bl1bloq)
  70. MASPLA1 tbas;
  71. * solutions statiques
  72. tbsta = resout (RIGPL1 et bloq1ry et bloq1rz et
  73. blq1uz et bliaiq1) tblsta ;
  74. react bliaiq1 tbsta ;
  75. sinon ;
  76. * modes vibrations
  77. MODPLA1=VIBR PROC PR n_mv (RIGPL1 et bloq1rx et bloq1rz et
  78. blq1uz) MASPLA1 tbas;
  79. * solutions statiques
  80. tbsta = resout RIGPL1 tblsta ;
  81. finsi ;
  82.  
  83. tmod = modpla1 . modes ;
  84. tmod (maspla1 et bl1bloq) * ;
  85.  
  86. repeter bmod ((dime modpla1 . modes) - 2) ;
  87. ind1 = &bmod ;
  88. psta1 = manu poi1 tmod . ind1 . point_repere ;
  89. momoda1 = mode psta1 mecanique elastique modal ;
  90. camoda1 = mate momoda1 freq tmod . ind1 . frequence
  91. mass tmod . ind1 . masse_generalisee
  92. defo tmod . ind1 . deformee_modale
  93. ;
  94. si (ega 1 ind1) ;
  95. momod1 = momoda1 ; camod1 = camoda1 ;
  96. sinon ;
  97. momod1 = momod1 et momoda1 ; camod1 = camod1 et camoda1 ;
  98. finsi ;
  99. def1 = defo tmod . ind1 . deformee_modale lig1 1. vert ;
  100. def0 = defo tmod . ind1 . deformee_modale lig1 0. bleu ;
  101. titre ' mode propre ind = ' ind1 ;
  102. si d_graph ; trac (10. 10. 0.) (def0 et def1) ; finsi ;
  103. fin bmod ;
  104.  
  105. riam1 = amor momod1 camod1 ;
  106.  
  107.  
  108. * forces
  109.  
  110. fy1 = manu chpo qf1 fy 1. ; fy1 = fy1 * 1.e4 ;
  111. fy2 = manu chpo q1 fy 1. ; fy2 = fy2 * 1.e4 ;
  112.  
  113. fq2 = manu chpo q2 fy 1.e4 ;
  114. dq2 = resout (RIGPL1 et bloq1ry et bloq1rz et
  115. blq1uz et bliaiq1) fq2 ;
  116. maaq2 = manu poi1 (8. 8. 8.) ;
  117.  
  118.  
  119. *
  120. abs_f = prog 0. 1000. ; ord_f = prog 1. 1. ;
  121. ev_xf = evol manu abs_f ord_f ;
  122. fy3 = pjba fy1 momod1 camod1 ;
  123. cha3 = char meca ev_xf fy3 ;
  124.  
  125. opti mode freq ;
  126. t_freq = table 'PASAPAS' ;
  127.  
  128. t_freq . modele = momod1 ;
  129. t_freq . caracteristiques = camod1 ;
  130. t_freq . chargement = cha3 ;
  131.  
  132. *t_freq . temps_calcules = prog 0. 0.1 ;
  133. pasapas t_freq ;
  134.  
  135. o_y = prog ; a_f = prog ; o_alp = prog ;
  136. repeter b_post ((dime t_freq.temps) - 1) ;
  137. ind1 = &b_post ;
  138. chdep1 = t_freq . deplacements . ind1 ;
  139. dep11 = reco t_freq . modele t_freq . caracteristiques chdep1 ;
  140.  
  141. uu_alp = xty chdep1 chdep1 (mots alfa beta ) (mots alfa beta) ;
  142. o_alp = o_alp et (prog uu_alp) ;
  143.  
  144. si faux ;
  145. deff1 = defo dep11 lig1 1. vert ;
  146. deff0 = defo dep11 lig1 0. bleu ;
  147. titre ' reponse frequentielle f = ' t_freq . temps . ind1 ;
  148. trac (10. 10. 0.) (deff0 et deff1) ;
  149. finsi ;
  150.  
  151. dep11p2 = redu dep11 qf1 ;
  152. uu_p2 = xty dep11p2 dep11p2
  153. (mots uy) (mots uy) ;
  154. o_y = o_y et (prog uu_p2);
  155. a_f = a_f et (prog t_freq . temps . ind1) ;
  156. fin b_post ;
  157. o_y = o_y ** 0.5 ;
  158. o_alp = o_alp ** 0.5 ;
  159.  
  160. ev_dsp0 = evol manu 'FREQ' a_f 'UY qf1' o_y coul vert;
  161. ev_alp0 = evol manu 'FREQ' a_f 'ALPH' o_alp coul vert;
  162.  
  163. si d_graph ;
  164. dess ev_dsp0 logx logy titre 'reponse frequentielle Q2' ;
  165. dess ev_alp0 logx logy titre 'reponse reel modal' ;
  166. finsi ;
  167. o_y0 = 2671.3 ; o_alp0 = 7973.7 ;
  168. err01 = maxi o_y - o_y0 ; err02 = maxi o_alp - o_alp0 ;
  169.  
  170. si (((abs err01) < 1.e-1) et ((abs err02) < 1.e-1)) ;
  171. erre 0 ; sinon ; erre 5 ; finsi ;
  172. *
  173. *
  174. *
  175. opti mode freq ;
  176. repeter bmod ((dime modpla1 . modes) - 2) ;
  177. ind1 = &bmod ;
  178. psta1 = manu poi1 tmod . ind1 . point_repere ;
  179. momoda1 = mode psta1 mecanique elastique modal ;
  180. camoda1 = mate momoda1 freq tmod . ind1 . frequence
  181. mass tmod . ind1 . masse_generalisee
  182. defo tmod . ind1 . deformee_modale
  183. amor 0.05 ;
  184. si (ega 1 ind1) ;
  185. momod1 = momoda1 ; camod1 = camoda1 ;
  186. sinon ;
  187. momod1 = momod1 et momoda1 ; camod1 = camod1 et camoda1 ;
  188. finsi ;
  189. def1 = defo tmod . ind1 . deformee_modale lig1 1. vert ;
  190. def0 = defo tmod . ind1 . deformee_modale lig1 0. bleu ;
  191. titre ' mode propre ind = ' ind1 ;
  192. si d_graph ; trac (10. 10. 0.) (def0 et def1) ; finsi ;
  193. fin bmod ;
  194.  
  195. riam1 = amor momod1 camod1 ;
  196.  
  197. fy6 = fy3 ;
  198. fiy6 = changer fy6 comp (mots falf fbet) (mots ifal ifbe) ;
  199. cha6 = char meca ev_xf (fy6 et fiy6) ;
  200. *list cha6 ;
  201.  
  202. t_freq = table 'PASAPAS' ;
  203.  
  204. t_freq . modele = momod1 ;
  205. t_freq . caracteristiques = camod1 ;
  206. t_freq . chargement = cha3 ;
  207.  
  208. *t_freq . temps_calcules = prog 0. 0.1 ;
  209. pasapas t_freq ;
  210.  
  211.  
  212. o_y = prog ; io_y = prog ; a_f = prog ; o_alp = prog ; o_ial = prog ;
  213. o_pha = prog ;
  214. repeter b_post ((dime t_freq.temps) - 1) ;
  215. ind1 = &b_post ;
  216. chdep1 = t_freq . deplacements . ind1 ;
  217. dep11 = reco t_freq . modele t_freq . caracteristiques chdep1 ;
  218.  
  219. uu_alp = xty chdep1 chdep1 (mots alfa beta ) (mots alfa beta) ;
  220. uu_ial = xty chdep1 chdep1 (mots ialf ibet) (mots ialf ibet) ;
  221. o_alp = o_alp et (prog uu_alp) ;o_ial = o_ial et (prog uu_ial) ;
  222.  
  223. si faux ;
  224. deff1 = defo dep11 lig1 1. vert ;
  225. deff0 = defo dep11 lig1 0. bleu ;
  226. titre ' reponse frequentielle f = ' t_freq . temps . ind1 ;
  227. trac (10. 10. 0.) (deff0 et deff1) ;
  228. finsi ;
  229.  
  230. dep11p2 = redu dep11 qf1 ;
  231. uu_p2 = xty dep11p2 dep11p2
  232. (mots uy iuy ) (mots uy iuy) ;
  233. im_p2 = xty dep11p2 dep11p2
  234. (mots iuy) (mots iuy) ;
  235. o_y = o_y et (prog uu_p2); io_y = io_y et (prog im_p2) ;
  236. o_pha = o_pha et
  237. (prog ((extr dep11p2 iuy qf1 )/ (extr dep11p2 uy qf1))) ;
  238. a_f = a_f et (prog t_freq . temps . ind1) ;
  239. fin b_post ;
  240.  
  241. o_y = o_y ** 0.5 ; io_y = io_y ** 0.5 ;
  242. o_alp = o_alp ** 0.5 ; o_ial = o_ial ** 0.5 ;
  243. o_pha = atg o_pha ;
  244.  
  245. ev_dsp = evol manu 'FREQ' a_f 'UY qf1' o_y ;
  246. ev_dsp1 = ev_dsp coul rouge ;
  247. ev_pha = evol manu 'FREQ' a_f 'PHASE' o_pha ;
  248. ev_pha1 = ev_pha coul rouge ;
  249. ev_idsp = evol manu 'FREQ' a_f 'iUY qf1' io_y ;
  250. ev_alp = evol manu 'FREQ' a_f 'ALPH' o_alp ;
  251. ev_alp1 = ev_alp coul rouge ;
  252. ev_ial = evol manu 'FREQ' a_f 'IALP' o_ial ;
  253. ev_ial1 = ev_ial coul rouge ;
  254.  
  255. si d_graph ;
  256. *dess ev_dsp logx logy titre 'reponse frequentielle qf1';
  257. dess (ev_dsp et ev_dsp0) logx logy titre 'reponse frequentielle qf1';
  258. dess ev_pha logx titre 'phase qf1'
  259. *dess ev_idsp logx logy titre 'reponse frequentielle ima qf1' ;
  260. dess (ev_alp et ev_alp0) logx logy titre 'reponse reel modal' ;
  261. dess ev_ial logx logy titre 'reponse ima modal' ;
  262. finsi ;
  263.  
  264. o_y3 = 20.666 ;
  265. o_alp3 = 30.549 ;
  266. o_pha3 = 89.434 ;
  267.  
  268. err1 = (maxi o_y) - o_y3 ; err2 = (maxi o_alp) - o_alp3 ;
  269. err3 = (maxi o_pha) - o_pha3 ;
  270.  
  271.  
  272. si (((abs err1) < 1.e-1) et ((abs err2) < 1.e-1) et
  273. ((abs err3) < 1.e-1)) ;
  274. erre 0 ; sinon ; erre 5 ; finsi ;
  275. *
  276. *
  277. t_freq = table 'PASAPAS' ;
  278. t_freq . chargement = cha6 ;
  279.  
  280. t_freq . modele = momod1 ;
  281. t_freq . caracteristiques = camod1 ;
  282.  
  283. *t_freq . temps_calcules = prog 0. 0.1 ;
  284. pasapas t_freq ;
  285.  
  286.  
  287. o_y = prog ; io_y = prog ; a_f = prog ; o_alp = prog ; o_ial = prog ;
  288. o_pha = prog ;
  289. repeter b_post ((dime t_freq.temps) - 1) ;
  290. ind1 = &b_post ;
  291. chdep1 = t_freq . deplacements . ind1 ;
  292. dep11 = reco t_freq . modele t_freq . caracteristiques chdep1 ;
  293.  
  294. uu_alp = xty chdep1 chdep1 (mots alfa beta ) (mots alfa beta) ;
  295. uu_ial = xty chdep1 chdep1 (mots ialf ibet) (mots ialf ibet) ;
  296. o_alp = o_alp et (prog uu_alp) ;o_ial = o_ial et (prog uu_ial) ;
  297.  
  298. si faux ;
  299. deff1 = defo dep11 lig1 1. vert ;
  300. deff0 = defo dep11 lig1 0. bleu ;
  301. titre ' reponse frequentielle f = ' t_freq . temps . ind1 ;
  302. trac (10. 10. 0.) (deff0 et deff1) ;
  303. finsi ;
  304.  
  305. dep11p2 = redu dep11 qf1 ;
  306. uu_p2 = xty dep11p2 dep11p2
  307. (mots uy iuy ) (mots uy iuy) ;
  308. im_p2 = xty dep11p2 dep11p2
  309. (mots iuy) (mots iuy) ;
  310. o_y = o_y et (prog uu_p2); io_y = io_y et (prog im_p2) ;
  311. o_pha = o_pha et
  312. (prog ((extr dep11p2 iuy qf1 )/ (extr dep11p2 uy qf1))) ;
  313. a_f = a_f et (prog t_freq . temps . ind1) ;
  314. fin b_post ;
  315.  
  316. o_y = o_y ** 0.5 ; io_y = io_y ** 0.5 ;
  317. o_alp = o_alp ** 0.5 ; o_ial = o_ial ** 0.5 ;
  318. * ceci n est pas le dephasage relatif --- ici ecart 45°
  319. o_pha = atg o_pha ;
  320.  
  321. ev_dsp = evol manu 'FREQ' a_f 'UY qf1' o_y ;
  322. ev_pha = evol manu 'FREQ' a_f 'PHASE' o_pha ;
  323. ev_idsp = evol manu 'FREQ' a_f 'iUY qf1' io_y ;
  324. ev_alp = evol manu 'FREQ' a_f 'ALPH' o_alp ;
  325. ev_ial = evol manu 'FREQ' a_f 'IALP' o_ial ;
  326.  
  327. si d_graph ;
  328. dess (ev_dsp et ev_dsp1 et ev_dsp0)
  329. logx logy titre 'reponse frequentielle qf1';
  330. dess (ev_pha et ev_pha1) logx titre 'phase qf1'
  331. *dess ev_idsp logx logy titre 'reponse frequentielle ima qf1' ;
  332. dess (ev_alp et ev_alp1 et ev_alp0)
  333. logx logy titre 'reponse reel modal' ;
  334. dess (ev_ial et ev_ial1) logx logy titre 'reponse ima modal' ;
  335. finsi ;
  336.  
  337. o_y5 = 29.226 ; o_alp5 = 69.570 ;
  338. *o_pha5 = 82.335 ; * kich 22/07/14
  339. o_pha5 = 47.503 ;
  340. err5 = (maxi o_y) - o_y5 ; err6 = (maxi o_alp) - o_alp5 ;
  341. err7 = (maxi o_pha) - o_pha5 ;
  342. si (((abs err5) < 1.e-1) et ((abs err6) < 1.e-1) et
  343. ((abs err7) < 1.e-1)) ;
  344. erre 0 ; sinon ; erre 5 ; finsi ;
  345.  
  346.  
  347.  
  348. fin ;
  349.  
  350.  
  351.  
  352.  
  353.  
  354.  
  355.  
  356.  
  357.  
  358.  
  359.  
  360.  
  361.  
  362.  

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