Télécharger pecker_f.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : pecker_f.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * etude frequentielle
  5. * chargement structure avec impedance
  6. *
  7. opti dime 3 elem seg2 trac x ;
  8. OPTION epsilon lineaire;
  9. *
  10. d_bloq = vrai;
  11. d_graph = faux ;
  12.  
  13.  
  14. p0 = 0. 0. 0. ; e_z = 0. 0. 1. ;
  15. q2 = e_z ; q1 = p0 ;
  16. lig1 = d 20 q2 q1 ;
  17. qf1 = point lig1 proc (0. 0. 0.6) ;
  18.  
  19. MO1= MODE LIG1 MECANIQUE ELASTIQUE POUT ;
  20.  
  21. MATPL1=MATE MO1 YOUNG 2.E10 NU 0.4 RHO 7800;
  22. CARPL1= 'CARA' MO1 SECT 0.1 INRY 2.08E-8 INRZ 3.33E-8
  23. TORS 2.41E-7;
  24. MATE1=MATPL1 ET CARPL1;
  25. ***
  26.  
  27.  
  28. ***
  29.  
  30. RIGPL1 =RIGI MATE1 MO1 ;
  31. MASPLA1=MASS MATE1 MO1 ;
  32. BLOQ2=BLOQ DEPL ROTA Q2 ;
  33. bloq1r = (bloq rota q1) ; blq1uz = bloq q1 uz ;
  34. bloq1rz = bloq q1 rz ;
  35. blq1ux = bloq q1 ux ; blq1uy = bloq q1 uy ;
  36. bloq1rx = bloq q1 rx ; bloq1ry= bloq q1 ry ;
  37. *
  38. *----
  39. *
  40. PR=PROG 0.01 ; n_mv = lect 4 ;
  41.  
  42. tblsta = table 'LIAISONS_STATIQUES' ;
  43. tblsta . 1 = table ; tblsta . 1 . point_liaison = q1 ;
  44. si d_bloq ; tblsta . 1 . ddl_liaison = ux ;
  45. sinon ; tblsta . 1 . ddl_liaison = fx ; finsi ;
  46. tblsta . 2 = table ; tblsta . 2 . point_liaison = q1 ;
  47. si d_bloq ; tblsta . 2 . ddl_liaison = uy ;
  48. sinon ; tblsta . 2 . ddl_liaison = fy ; finsi ;
  49. tblsta . 3 = table ; tblsta . 3 . point_liaison = q1 ;
  50. si d_bloq ; tblsta . 3 . ddl_liaison = rx ;
  51. sinon ; tblsta . 3 . ddl_liaison = mx ; finsi ;
  52.  
  53. si d_bloq ;
  54. bliaiq1 = bloque tblsta ; depi tblsta ;
  55. sinon ;
  56. force tblsta ;
  57. finsi ;
  58. si d_bloq ;
  59. * modes vibrations
  60. bl1bloq = bloq1ry et bloq1rz et blq1uz et bliaiq1 ;
  61. MODPLA1=VIBR PROC PR n_mv (RIGPL1 et bl1bloq)
  62. MASPLA1 tbas;
  63. * solutions statiques
  64. tbsta = resout (RIGPL1 et bloq1ry et bloq1rz et
  65. blq1uz et bliaiq1) tblsta ;
  66. react bliaiq1 tbsta ;
  67. sinon ;
  68. * modes vibrations
  69. MODPLA1=VIBR PROC PR n_mv (RIGPL1 et bloq1rx et bloq1rz et
  70. blq1uz) MASPLA1 tbas;
  71. * solutions statiques
  72. tbsta = resout RIGPL1 tblsta ;
  73. finsi ;
  74.  
  75. tbsta maspla1 * ;
  76. opti mode freq ;
  77.  
  78. repeter bmod ((dime tbsta) - 1) ;
  79. ind1 = &bmod ;
  80. psta1 = manu poi1 tbsta . ind1 . point_repere ;
  81. mostai1 = mode psta1 mecanique elastique statique ;
  82. si d_bloq ;
  83. castai1 = mate mostai1 defo tbsta . ind1 . deformee
  84. made tbsta . ind1 . masse_deformee
  85. ride tbsta . ind1 . reaction ;
  86. sinon ;
  87. castai1 = mate mostai1 defo tbsta . ind1 . deformee
  88. made tbsta . ind1 . masse_deformee
  89. ride tbsta . ind1 . force ;
  90. finsi ;
  91. si (ega 1 ind1) ;
  92. mosta1 = mostai1 ; casta1 = castai1 ;
  93. sinon ;
  94. mosta1 = mosta1 et mostai1 ; casta1 = casta1 et castai1 ;
  95. finsi ;
  96. def1 = defo tbsta . ind1 . deformee lig1 1. vert ;
  97. def0 = defo tbsta . ind1 . deformee lig1 0. bleu ;
  98. titre ' sol statique ind = ' ind1 ;
  99. si d_graph ;
  100. trac (10. 10. 0.) (def0 et def1) ;
  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. tmod (maspla1 et bl1bloq) * ;
  108.  
  109. maxfre5 = 0. ;
  110.  
  111. repeter bmod ((dime modpla1 . modes) - 2) ;
  112. ind1 = &bmod ;
  113. psta1 = manu poi1 tmod . ind1 . point_repere ;
  114. momoda1 = mode psta1 mecanique elastique modal ;
  115. camoda1 = mate momoda1 freq tmod . ind1 . frequence
  116. mass tmod . ind1 . masse_generalisee
  117. defo tmod . ind1 . deformee_modale
  118. made tmod . ind1 . masse_deformee ;
  119. maxfre5 = maxi (prog maxfre5 tmod . ind1 . frequence) ;
  120. si (ega 1 ind1) ;
  121. momod1 = momoda1 ; camod1 = camoda1 ;
  122. sinon ;
  123. momod1 = momod1 et momoda1 ; camod1 = camod1 et camoda1 ;
  124. finsi ;
  125. def1 = defo tmod . ind1 . deformee_modale lig1 1. vert ;
  126. def0 = defo tmod . ind1 . deformee_modale lig1 0. bleu ;
  127. titre ' mode propre ind = ' ind1 ;
  128. si d_graph ; trac (10. 10. 0.) (def0 et def1) ; finsi ;
  129. fin bmod ;
  130.  
  131. modto1 = modto1 et momod1 ; casto1 = casto1 et camod1 ;
  132.  
  133. * forces
  134.  
  135. fy1 = manu chpo qf1 fy 1. ; fy1 = fy1 * 1.e4 ;
  136. fy2 = manu chpo q1 fy 1. ; fy2 = fy2 * 1.e4 ;
  137.  
  138. fq2 = manu chpo q2 fy 1.e4 ;
  139. dq2 = resout (RIGPL1 et bloq1ry et bloq1rz et
  140. blq1uz et bliaiq1) fq2 ;
  141. maaq2 = manu poi1 (8. 8. 8.) ;
  142. mosq2 = modele maaq2 mecanique elastique statique ;
  143. casq2 = mate mosq2 defo dq2 ride fq2 made (maspla1 * dq2) ;
  144. def1 = defo dq2 lig1 1. vert ;
  145. def0 = defo dq2 lig1 0. bleu ;
  146. titre ' sol statique q2 fy 1.e4' ;
  147. si d_graph ;
  148. trac (10. 10. 0.) (def0 et def1) ;
  149. finsi ;
  150.  
  151. * frequentiel
  152. abs_f = prog 0. 1000. ; ord_f = prog 1. 1. ;
  153. *abs_f = prog 1. pas 0.3 4. pas 1. 10. ;
  154. *abs_f = l_zf ; ord_f = prog (dime abs_f)*1. ;
  155.  
  156. mimcoq2 = modele maaq2 mecanique elastique statique impedance complexe;
  157. evim2 = (evol manu 'FREQ' abs_f 'MOCO' (prog 1.e3 2.e3))
  158. et (evol manu 'FREQ' abs_f 'VISC' (prog 2.e3 1.e3)) ;
  159. cimcoq2 = mate mimcoq2 'MOCO' evim2 ;
  160. cimcoq3 = mate mimcoq2 'RAID' 1.e3 'VISC' 2.e4 ;
  161.  
  162. ev_xf = evol manu abs_f ord_f ;
  163. fy3 = pjba fy1 momod1 camod1 ;
  164. cha3 = char meca ev_xf fy3 ;
  165. fy4 = pjba fy1 mosta1 casta1 ;
  166. cha4 = char meca ev_xf fy4 ;
  167. fy5 = pjba fy1 mosq2 casq2 ;
  168. cha5 = char meca ev_xf fy5 ;
  169. fy6 = pjba fy2 (mosq2 et mosta1 et momod1) (casq2 et casta1 et camod1);
  170. cha6 = char meca ev_xf fy6 ;
  171. *list cha6 ;
  172.  
  173. t_freq = table 'PASAPAS' ;
  174. t_freq . modele = mosq2 et mosta1 et momod1 et mimcoq2 ;
  175. t_freq . caracteristiques = casq2 et casta1 et camod1 et cimcoq3 ;
  176. *t_freq . chargement = cha3 et cha5 et cha4 ;
  177. t_freq . chargement = cha6 ;
  178. *t_freq . modele = mosq2 ;
  179. *t_freq . caracteristiques = casq2 ;
  180. *t_freq . chargement = cha5 ;
  181. *t_freq . modele = momod1 ;
  182. *t_freq . caracteristiques = camod1 ;
  183. *t_freq . chargement = cha3 ;
  184. *t_freq . temps_calcules = abs_f ;
  185. * par defaut pas_defa propose une liste ; ensuite on affine
  186.  
  187. pasapas t_freq ;
  188.  
  189. o_y = prog ; a_f = prog ; o_iy = prog ;
  190. repeter b_post ((dime t_freq.temps) - 1) ;
  191. ind1 = &b_post ;
  192. chdep1 = t_freq . deplacements . ind1 ;
  193. dep11 = reco t_freq . modele t_freq . caracteristiques chdep1 ;
  194.  
  195. si faux ;
  196. deff1 = defo dep11 lig1 1. vert ;
  197. deff0 = defo dep11 lig1 0. bleu ;
  198. titre ' reponse frequentielle f = ' t_freq . temps . ind1 ;
  199. trac (10. 10. 0.) (deff0 et deff1) ;
  200. finsi ;
  201.  
  202. dep11p2 = redu dep11 q2 ;
  203. uu_p2 = xty dep11p2 dep11p2
  204. (mots ux uy uz) (mots ux uy uz) ;
  205. ui_p2 = xty dep11p2 dep11p2
  206. (mots iux iuy iuz) (mots iux iuy iuz) ;
  207.  
  208. o_y = o_y et (prog uu_p2 ); o_iy = o_iy et (prog ui_p2 );
  209. a_f = a_f et (prog t_freq . temps . ind1) ;
  210. fin b_post ;
  211. o_y = o_y ** 0.5 ;o_iy = o_iy ** 0.5 ;
  212.  
  213. ev_dsp = evol manu 'FREQ' a_f 'DEPL' o_y ;
  214. ei_dsp = evol manu 'FREQ' a_f 'DEPL' o_iy ;
  215. si d_graph ;
  216. dess ev_dsp logx logy titre 'reponse frequentielle Q2' ;
  217. dess ei_dsp logx logy titre 'reponse frequentielle Q2 ima' ;
  218. finsi ;
  219.  
  220. * test t = 4. avec impedance
  221. o_y4 = .41974 ;
  222. si (abs((ipol 4. ev_dsp) - o_y4) > (o_y4*1.e-2)) ;
  223. erre 5 ;
  224. sinon ;
  225. erre 0 ;
  226. finsi ;
  227.  
  228. fin ;
  229.  
  230.  
  231.  
  232.  
  233.  
  234.  
  235.  
  236.  
  237.  
  238.  
  239.  
  240.  
  241.  
  242.  
  243.  
  244.  

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