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) MASPLA1 ;
  70. * solutions statiques
  71. tbsta = resout (RIGPL1 et bloq1ry et bloq1rz et
  72. blq1uz et bliaiq1) tblsta ;
  73. react bliaiq1 tbsta ;
  74. sinon ;
  75. * modes vibrations
  76. MODPLA1=VIBR PROC PR n_mv (RIGPL1 et bloq1rx et bloq1rz et
  77. blq1uz) MASPLA1 ;
  78. * solutions statiques
  79. tbsta = resout RIGPL1 tblsta ;
  80. finsi ;
  81.  
  82. tmod = modpla1 . modes ;
  83. tmod (maspla1 et bl1bloq) * ;
  84.  
  85. repeter bmod ((dime modpla1 . modes) - 2) ;
  86. ind1 = &bmod ;
  87. psta1 = manu poi1 tmod . ind1 . point_repere ;
  88. momoda1 = mode psta1 mecanique elastique modal ;
  89. camoda1 = mate momoda1 freq tmod . ind1 . frequence
  90. mass tmod . ind1 . masse_generalisee
  91. defo tmod . ind1 . deformee_modale
  92. ;
  93. si (ega 1 ind1) ;
  94. momod1 = momoda1 ; camod1 = camoda1 ;
  95. sinon ;
  96. momod1 = momod1 et momoda1 ; camod1 = camod1 et camoda1 ;
  97. finsi ;
  98. def1 = defo tmod . ind1 . deformee_modale lig1 1. vert ;
  99. def0 = defo tmod . ind1 . deformee_modale lig1 0. bleu ;
  100. titre ' mode propre ind = ' ind1 ;
  101. si d_graph ; trac (10. 10. 0.) (def0 et def1) ; finsi ;
  102. fin bmod ;
  103.  
  104. riam1 = amor momod1 camod1 ;
  105.  
  106.  
  107. * forces
  108.  
  109. fy1 = manu chpo qf1 fy 1. ; fy1 = fy1 * 1.e4 ;
  110. fy2 = manu chpo q1 fy 1. ; fy2 = fy2 * 1.e4 ;
  111.  
  112. fq2 = manu chpo q2 fy 1.e4 ;
  113. dq2 = resout (RIGPL1 et bloq1ry et bloq1rz et
  114. blq1uz et bliaiq1) fq2 ;
  115. maaq2 = manu poi1 (8. 8. 8.) ;
  116.  
  117.  
  118. *
  119. abs_f = prog 0. 1000. ; ord_f = prog 1. 1. ;
  120. ev_xf = evol manu abs_f ord_f ;
  121. fy3 = pjba fy1 momod1 camod1 ;
  122. cha3 = char meca ev_xf fy3 ;
  123.  
  124. opti mode freq ;
  125. t_freq = table 'PASAPAS' ;
  126.  
  127. t_freq . modele = momod1 ;
  128. t_freq . caracteristiques = camod1 ;
  129. t_freq . chargement = cha3 ;
  130.  
  131. *t_freq . temps_calcules = prog 0. 0.1 ;
  132. pasapas t_freq ;
  133.  
  134. o_y = prog ; a_f = prog ; o_alp = prog ;
  135. repeter b_post ((dime t_freq.temps) - 1) ;
  136. ind1 = &b_post ;
  137. chdep1 = t_freq . deplacements . ind1 ;
  138. dep11 = reco t_freq . modele t_freq . caracteristiques chdep1 ;
  139.  
  140. uu_alp = xty chdep1 chdep1 (mots alfa beta ) (mots alfa beta) ;
  141. o_alp = o_alp et (prog uu_alp) ;
  142.  
  143. si faux ;
  144. deff1 = defo dep11 lig1 1. vert ;
  145. deff0 = defo dep11 lig1 0. bleu ;
  146. titre ' reponse frequentielle f = ' t_freq . temps . ind1 ;
  147. trac (10. 10. 0.) (deff0 et deff1) ;
  148. finsi ;
  149.  
  150. dep11p2 = redu dep11 qf1 ;
  151. uu_p2 = xty dep11p2 dep11p2
  152. (mots uy) (mots uy) ;
  153. o_y = o_y et (prog uu_p2);
  154. a_f = a_f et (prog t_freq . temps . ind1) ;
  155. fin b_post ;
  156. o_y = o_y ** 0.5 ;
  157. o_alp = o_alp ** 0.5 ;
  158.  
  159. ev_dsp0 = evol manu 'FREQ' a_f 'UY qf1' o_y coul vert;
  160. ev_alp0 = evol manu 'FREQ' a_f 'ALPH' o_alp coul vert;
  161.  
  162. si d_graph ;
  163. dess ev_dsp0 logx logy titre 'reponse frequentielle Q2' ;
  164. dess ev_alp0 logx logy titre 'reponse reel modal' ;
  165. finsi ;
  166. o_y0 = 2671.3 ; o_alp0 = 7973.7 ;
  167. err01 = maxi o_y - o_y0 ; err02 = maxi o_alp - o_alp0 ;
  168.  
  169. si (((abs err01) < 1.e-1) et ((abs err02) < 1.e-1)) ;
  170. erre 0 ; sinon ; erre 5 ; finsi ;
  171. *
  172. *
  173. *
  174. opti mode freq ;
  175. repeter bmod ((dime modpla1 . modes) - 2) ;
  176. ind1 = &bmod ;
  177. psta1 = manu poi1 tmod . ind1 . point_repere ;
  178. momoda1 = mode psta1 mecanique elastique modal ;
  179. camoda1 = mate momoda1 freq tmod . ind1 . frequence
  180. mass tmod . ind1 . masse_generalisee
  181. defo tmod . ind1 . deformee_modale
  182. amor 0.05 ;
  183. si (ega 1 ind1) ;
  184. momod1 = momoda1 ; camod1 = camoda1 ;
  185. sinon ;
  186. momod1 = momod1 et momoda1 ; camod1 = camod1 et camoda1 ;
  187. finsi ;
  188. def1 = defo tmod . ind1 . deformee_modale lig1 1. vert ;
  189. def0 = defo tmod . ind1 . deformee_modale lig1 0. bleu ;
  190. titre ' mode propre ind = ' ind1 ;
  191. si d_graph ; trac (10. 10. 0.) (def0 et def1) ; finsi ;
  192. fin bmod ;
  193.  
  194. riam1 = amor momod1 camod1 ;
  195.  
  196. fy6 = fy3 ;
  197. fiy6 = changer fy6 comp (mots falf fbet) (mots ifal ifbe) ;
  198. cha6 = char meca ev_xf (fy6 et fiy6) ;
  199. *list cha6 ;
  200.  
  201. t_freq = table 'PASAPAS' ;
  202.  
  203. t_freq . modele = momod1 ;
  204. t_freq . caracteristiques = camod1 ;
  205. t_freq . chargement = cha3 ;
  206.  
  207. *t_freq . temps_calcules = prog 0. 0.1 ;
  208. pasapas t_freq ;
  209.  
  210.  
  211. o_y = prog ; io_y = prog ; a_f = prog ; o_alp = prog ; o_ial = prog ;
  212. o_pha = prog ;
  213. repeter b_post ((dime t_freq.temps) - 1) ;
  214. ind1 = &b_post ;
  215. chdep1 = t_freq . deplacements . ind1 ;
  216. dep11 = reco t_freq . modele t_freq . caracteristiques chdep1 ;
  217.  
  218. uu_alp = xty chdep1 chdep1 (mots alfa beta ) (mots alfa beta) ;
  219. uu_ial = xty chdep1 chdep1 (mots ialf ibet) (mots ialf ibet) ;
  220. o_alp = o_alp et (prog uu_alp) ;o_ial = o_ial et (prog uu_ial) ;
  221.  
  222. si faux ;
  223. deff1 = defo dep11 lig1 1. vert ;
  224. deff0 = defo dep11 lig1 0. bleu ;
  225. titre ' reponse frequentielle f = ' t_freq . temps . ind1 ;
  226. trac (10. 10. 0.) (deff0 et deff1) ;
  227. finsi ;
  228.  
  229. dep11p2 = redu dep11 qf1 ;
  230. uu_p2 = xty dep11p2 dep11p2
  231. (mots uy iuy ) (mots uy iuy) ;
  232. im_p2 = xty dep11p2 dep11p2
  233. (mots iuy) (mots iuy) ;
  234. o_y = o_y et (prog uu_p2); io_y = io_y et (prog im_p2) ;
  235. o_pha = o_pha et
  236. (prog ((extr dep11p2 iuy qf1 )/ (extr dep11p2 uy qf1))) ;
  237. a_f = a_f et (prog t_freq . temps . ind1) ;
  238. fin b_post ;
  239.  
  240. o_y = o_y ** 0.5 ; io_y = io_y ** 0.5 ;
  241. o_alp = o_alp ** 0.5 ; o_ial = o_ial ** 0.5 ;
  242. o_pha = atg o_pha ;
  243.  
  244. ev_dsp = evol manu 'FREQ' a_f 'UY qf1' o_y ;
  245. ev_dsp1 = ev_dsp coul rouge ;
  246. ev_pha = evol manu 'FREQ' a_f 'PHASE' o_pha ;
  247. ev_pha1 = ev_pha coul rouge ;
  248. ev_idsp = evol manu 'FREQ' a_f 'iUY qf1' io_y ;
  249. ev_alp = evol manu 'FREQ' a_f 'ALPH' o_alp ;
  250. ev_alp1 = ev_alp coul rouge ;
  251. ev_ial = evol manu 'FREQ' a_f 'IALP' o_ial ;
  252. ev_ial1 = ev_ial coul rouge ;
  253.  
  254. si d_graph ;
  255. *dess ev_dsp logx logy titre 'reponse frequentielle qf1';
  256. dess (ev_dsp et ev_dsp0) logx logy titre 'reponse frequentielle qf1';
  257. dess ev_pha logx titre 'phase qf1'
  258. *dess ev_idsp logx logy titre 'reponse frequentielle ima qf1' ;
  259. dess (ev_alp et ev_alp0) logx logy titre 'reponse reel modal' ;
  260. dess ev_ial logx logy titre 'reponse ima modal' ;
  261. finsi ;
  262.  
  263. o_y3 = 20.666 ;
  264. o_alp3 = 30.549 ;
  265. o_pha3 = 89.434 ;
  266.  
  267. err1 = (maxi o_y) - o_y3 ; err2 = (maxi o_alp) - o_alp3 ;
  268. err3 = (maxi o_pha) - o_pha3 ;
  269.  
  270.  
  271. si (((abs err1) < 1.e-1) et ((abs err2) < 1.e-1) et
  272. ((abs err3) < 1.e-1)) ;
  273. erre 0 ; sinon ; erre 5 ; finsi ;
  274. *
  275. *
  276. t_freq = table 'PASAPAS' ;
  277. t_freq . chargement = cha6 ;
  278.  
  279. t_freq . modele = momod1 ;
  280. t_freq . caracteristiques = camod1 ;
  281.  
  282. *t_freq . temps_calcules = prog 0. 0.1 ;
  283. pasapas t_freq ;
  284.  
  285.  
  286. o_y = prog ; io_y = prog ; a_f = prog ; o_alp = prog ; o_ial = prog ;
  287. o_pha = prog ;
  288. repeter b_post ((dime t_freq.temps) - 1) ;
  289. ind1 = &b_post ;
  290. chdep1 = t_freq . deplacements . ind1 ;
  291. dep11 = reco t_freq . modele t_freq . caracteristiques chdep1 ;
  292.  
  293. uu_alp = xty chdep1 chdep1 (mots alfa beta ) (mots alfa beta) ;
  294. uu_ial = xty chdep1 chdep1 (mots ialf ibet) (mots ialf ibet) ;
  295. o_alp = o_alp et (prog uu_alp) ;o_ial = o_ial et (prog uu_ial) ;
  296.  
  297. si faux ;
  298. deff1 = defo dep11 lig1 1. vert ;
  299. deff0 = defo dep11 lig1 0. bleu ;
  300. titre ' reponse frequentielle f = ' t_freq . temps . ind1 ;
  301. trac (10. 10. 0.) (deff0 et deff1) ;
  302. finsi ;
  303.  
  304. dep11p2 = redu dep11 qf1 ;
  305. uu_p2 = xty dep11p2 dep11p2
  306. (mots uy iuy ) (mots uy iuy) ;
  307. im_p2 = xty dep11p2 dep11p2
  308. (mots iuy) (mots iuy) ;
  309. o_y = o_y et (prog uu_p2); io_y = io_y et (prog im_p2) ;
  310. o_pha = o_pha et
  311. (prog ((extr dep11p2 iuy qf1 )/ (extr dep11p2 uy qf1))) ;
  312. a_f = a_f et (prog t_freq . temps . ind1) ;
  313. fin b_post ;
  314.  
  315. o_y = o_y ** 0.5 ; io_y = io_y ** 0.5 ;
  316. o_alp = o_alp ** 0.5 ; o_ial = o_ial ** 0.5 ;
  317. * ceci n est pas le dephasage relatif --- ici ecart 45°
  318. o_pha = atg o_pha ;
  319.  
  320. ev_dsp = evol manu 'FREQ' a_f 'UY qf1' o_y ;
  321. ev_pha = evol manu 'FREQ' a_f 'PHASE' o_pha ;
  322. ev_idsp = evol manu 'FREQ' a_f 'iUY qf1' io_y ;
  323. ev_alp = evol manu 'FREQ' a_f 'ALPH' o_alp ;
  324. ev_ial = evol manu 'FREQ' a_f 'IALP' o_ial ;
  325.  
  326. si d_graph ;
  327. dess (ev_dsp et ev_dsp1 et ev_dsp0)
  328. logx logy titre 'reponse frequentielle qf1';
  329. dess (ev_pha et ev_pha1) logx titre 'phase qf1'
  330. *dess ev_idsp logx logy titre 'reponse frequentielle ima qf1' ;
  331. dess (ev_alp et ev_alp1 et ev_alp0)
  332. logx logy titre 'reponse reel modal' ;
  333. dess (ev_ial et ev_ial1) logx logy titre 'reponse ima modal' ;
  334. finsi ;
  335.  
  336. o_y5 = 29.226 ; o_alp5 = 69.570 ;
  337. *o_pha5 = 82.335 ; * kich 22/07/14
  338. o_pha5 = 47.503 ;
  339. err5 = (maxi o_y) - o_y5 ; err6 = (maxi o_alp) - o_alp5 ;
  340. err7 = (maxi o_pha) - o_pha5 ;
  341. si (((abs err5) < 1.e-1) et ((abs err6) < 1.e-1) et
  342. ((abs err7) < 1.e-1)) ;
  343. erre 0 ; sinon ; erre 5 ; finsi ;
  344.  
  345.  
  346.  
  347. fin ;
  348.  
  349.  
  350.  
  351.  
  352.  
  353.  
  354.  
  355.  
  356.  
  357.  
  358.  
  359.  
  360.  
  361.  
  362.  
  363.  

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