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. *
  9. d_bloq = vrai;
  10. d_graph = faux ;
  11.  
  12.  
  13. p0 = 0. 0. 0. ; e_z = 0. 0. 1. ;
  14. q2 = e_z ; q1 = p0 ;
  15. lig1 = d 20 q2 q1 ;
  16. qf1 = point lig1 proc (0. 0. 0.6) ;
  17.  
  18. MO1= MODE LIG1 MECANIQUE ELASTIQUE POUT ;
  19.  
  20. MATPL1=MATE MO1 YOUNG 2.E10 NU 0.4 RHO 7800;
  21. CARPL1= 'CARA' MO1 SECT 0.1 INRY 2.08E-8 INRZ 3.33E-8
  22. TORS 2.41E-7;
  23. MATE1=MATPL1 ET CARPL1;
  24. ***
  25.  
  26.  
  27. ***
  28.  
  29. RIGPL1 =RIGI MATE1 MO1 ;
  30. MASPLA1=MASS MATE1 MO1 ;
  31. BLOQ2=BLOQ DEPL ROTA Q2 ;
  32. bloq1r = (bloq rota q1) ; blq1uz = bloq q1 uz ;
  33. bloq1rz = bloq q1 rz ;
  34. blq1ux = bloq q1 ux ; blq1uy = bloq q1 uy ;
  35. bloq1rx = bloq q1 rx ; bloq1ry= bloq q1 ry ;
  36. *
  37. *----
  38. *
  39. PR=PROG 0.01 ; n_mv = lect 4 ;
  40.  
  41. tblsta = table 'LIAISONS_STATIQUES' ;
  42. tblsta . 1 = table ; tblsta . 1 . point_liaison = q1 ;
  43. si d_bloq ; tblsta . 1 . ddl_liaison = ux ;
  44. sinon ; tblsta . 1 . ddl_liaison = fx ; finsi ;
  45. tblsta . 2 = table ; tblsta . 2 . point_liaison = q1 ;
  46. si d_bloq ; tblsta . 2 . ddl_liaison = uy ;
  47. sinon ; tblsta . 2 . ddl_liaison = fy ; finsi ;
  48. tblsta . 3 = table ; tblsta . 3 . point_liaison = q1 ;
  49. si d_bloq ; tblsta . 3 . ddl_liaison = rx ;
  50. sinon ; tblsta . 3 . ddl_liaison = mx ; finsi ;
  51.  
  52. si d_bloq ;
  53. bliaiq1 = bloque tblsta ; depi tblsta ;
  54. sinon ;
  55. force tblsta ;
  56. finsi ;
  57. si d_bloq ;
  58. * modes vibrations
  59. bl1bloq = bloq1ry et bloq1rz et blq1uz et bliaiq1 ;
  60. MODPLA1=VIBR PROC PR n_mv (RIGPL1 et bl1bloq)
  61. MASPLA1 ;
  62. * solutions statiques
  63. tbsta = resout (RIGPL1 et bloq1ry et bloq1rz et
  64. blq1uz et bliaiq1) tblsta ;
  65. react bliaiq1 tbsta ;
  66. sinon ;
  67. * modes vibrations
  68. MODPLA1=VIBR PROC PR n_mv (RIGPL1 et bloq1rx et bloq1rz et
  69. blq1uz) MASPLA1 ;
  70. * solutions statiques
  71. tbsta = resout RIGPL1 tblsta ;
  72. finsi ;
  73.  
  74. tbsta maspla1 * ;
  75. opti mode freq ;
  76.  
  77. repeter bmod ((dime tbsta) - 1) ;
  78. ind1 = &bmod ;
  79. psta1 = manu poi1 tbsta . ind1 . point_repere ;
  80. mostai1 = mode psta1 mecanique elastique statique ;
  81. si d_bloq ;
  82. castai1 = mate mostai1 defo tbsta . ind1 . deformee
  83. made tbsta . ind1 . masse_deformee
  84. ride tbsta . ind1 . reaction ;
  85. sinon ;
  86. castai1 = mate mostai1 defo tbsta . ind1 . deformee
  87. made tbsta . ind1 . masse_deformee
  88. ride tbsta . ind1 . force ;
  89. finsi ;
  90. si (ega 1 ind1) ;
  91. mosta1 = mostai1 ; casta1 = castai1 ;
  92. sinon ;
  93. mosta1 = mosta1 et mostai1 ; casta1 = casta1 et castai1 ;
  94. finsi ;
  95. def1 = defo tbsta . ind1 . deformee lig1 1. vert ;
  96. def0 = defo tbsta . ind1 . deformee lig1 0. bleu ;
  97. titre ' sol statique ind = ' ind1 ;
  98. si d_graph ;
  99. trac (10. 10. 0.) (def0 et def1) ;
  100. finsi ;
  101. fin bmod ;
  102. masta1 = masse mosta1 casta1 ;
  103. rista1 = rigi mosta1 casta1 ;
  104.  
  105. modto1 = mosta1 ; casto1 = casta1 ; tmod = modpla1 . modes ;
  106. tmod (maspla1 et bl1bloq) * ;
  107.  
  108. maxfre5 = 0. ;
  109.  
  110. repeter bmod ((dime modpla1 . modes) - 2) ;
  111. ind1 = &bmod ;
  112. psta1 = manu poi1 tmod . ind1 . point_repere ;
  113. momoda1 = mode psta1 mecanique elastique modal ;
  114. camoda1 = mate momoda1 freq tmod . ind1 . frequence
  115. mass tmod . ind1 . masse_generalisee
  116. defo tmod . ind1 . deformee_modale
  117. made tmod . ind1 . masse_deformee ;
  118. maxfre5 = maxi (prog maxfre5 tmod . ind1 . frequence) ;
  119. si (ega 1 ind1) ;
  120. momod1 = momoda1 ; camod1 = camoda1 ;
  121. sinon ;
  122. momod1 = momod1 et momoda1 ; camod1 = camod1 et camoda1 ;
  123. finsi ;
  124. def1 = defo tmod . ind1 . deformee_modale lig1 1. vert ;
  125. def0 = defo tmod . ind1 . deformee_modale lig1 0. bleu ;
  126. titre ' mode propre ind = ' ind1 ;
  127. si d_graph ; trac (10. 10. 0.) (def0 et def1) ; finsi ;
  128. fin bmod ;
  129.  
  130. modto1 = modto1 et momod1 ; casto1 = casto1 et camod1 ;
  131.  
  132. * forces
  133.  
  134. fy1 = manu chpo qf1 fy 1. ; fy1 = fy1 * 1.e4 ;
  135. fy2 = manu chpo q1 fy 1. ; fy2 = fy2 * 1.e4 ;
  136.  
  137. fq2 = manu chpo q2 fy 1.e4 ;
  138. dq2 = resout (RIGPL1 et bloq1ry et bloq1rz et
  139. blq1uz et bliaiq1) fq2 ;
  140. maaq2 = manu poi1 (8. 8. 8.) ;
  141. mosq2 = modele maaq2 mecanique elastique statique ;
  142. casq2 = mate mosq2 defo dq2 ride fq2 made (maspla1 * dq2) ;
  143. def1 = defo dq2 lig1 1. vert ;
  144. def0 = defo dq2 lig1 0. bleu ;
  145. titre ' sol statique q2 fy 1.e4' ;
  146. si d_graph ;
  147. trac (10. 10. 0.) (def0 et def1) ;
  148. finsi ;
  149.  
  150. * frequentiel
  151. abs_f = prog 0. 1000. ; ord_f = prog 1. 1. ;
  152. *abs_f = prog 1. pas 0.3 4. pas 1. 10. ;
  153. *abs_f = l_zf ; ord_f = prog (dime abs_f)*1. ;
  154.  
  155. mimcoq2 = modele maaq2 mecanique elastique statique impedance complexe;
  156. evim2 = (evol manu 'FREQ' abs_f 'MOCO' (prog 1.e3 2.e3))
  157. et (evol manu 'FREQ' abs_f 'VISC' (prog 2.e3 1.e3)) ;
  158. cimcoq2 = mate mimcoq2 'MOCO' evim2 ;
  159. cimcoq3 = mate mimcoq2 'RAID' 1.e3 'VISC' 2.e4 ;
  160.  
  161. ev_xf = evol manu abs_f ord_f ;
  162. fy3 = pjba fy1 momod1 camod1 ;
  163. cha3 = char meca ev_xf fy3 ;
  164. fy4 = pjba fy1 mosta1 casta1 ;
  165. cha4 = char meca ev_xf fy4 ;
  166. fy5 = pjba fy1 mosq2 casq2 ;
  167. cha5 = char meca ev_xf fy5 ;
  168. fy6 = pjba fy2 (mosq2 et mosta1 et momod1) (casq2 et casta1 et camod1);
  169. cha6 = char meca ev_xf fy6 ;
  170. *list cha6 ;
  171.  
  172. t_freq = table 'PASAPAS' ;
  173. t_freq . modele = mosq2 et mosta1 et momod1 et mimcoq2 ;
  174. t_freq . caracteristiques = casq2 et casta1 et camod1 et cimcoq3 ;
  175. *t_freq . chargement = cha3 et cha5 et cha4 ;
  176. t_freq . chargement = cha6 ;
  177. t_freq .hypothese_deformations = 'LINEAIRE' ;
  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.  
  245.  
  246.  
  247.  
  248.  

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