Télécharger la plaquette Cast3M


fr en

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

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