Télécharger la plaquette Cast3M


fr en

Télécharger four9.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * calcul coq2 elasto-plastique mode FOURIER
  2. * exploitation : flambage 
  3. * sol statique elastique initiale et modes de flambage
  4. * sol elastoplastique et modes de flambage
  5.  
  6. graph = 'EGA' 'SANS' 'VRAI' ;
  7. opti dime 3 elem seg2 ;
  8. p_ori = 0. 0. 0. ;
  9. e_x = 1. 0. 0. ; e_y = 0. 1. 0. ;  e_z = 0. 0. 1. ;
  10. e_r = e_x plus p_ori ;
  11.  
  12. R_hb = 2.5d0 ; L_m = 750.d-3 ; H_c = 200.d-3 ; L_s = 100.d-3 ;
  13. e_b = 25.d-3 ; e_a = 10.d-3 ; e_s = 25.d-3 ;
  14. d1 = L_m / 3 ;
  15. opti dens d1 ;
  16.  
  17. px_0 = p_ori plus p_ori ;
  18. pr_0 = px_0 plus (R_hb * e_r) ;
  19. N_m = 34 ;
  20. mess ' demi-longueur principale ' (N_m * L_m) ;
  21. px_1 = px_0 plus (N_m * L_m * e_z) ;
  22. * 
  23. pr_v1 = px_1 plus (R_hb * e_r) ;
  24. pd_1 = px_1 plus (R_hb * e_z) ;
  25. s_d = c pr_v1 px_1 pd_1 ;
  26. s_b = d N_m pr_0 pr_v1 ;
  27. s_0 = s_b et s_d ;
  28.  
  29. elim 1.e-4 s_0 ;
  30. 'SI' graph ;
  31. titre ' maillage 2D initial ';
  32. trac (s_0 et px_0 et px_1) 'QUAL' nclk ;
  33. 'FINSI' ;
  34. ********************************
  35. s_d2 s_b2 = s_d s_b syme plan px_0 pr_0 (px_0 plus e_y) ;
  36. pd_2 = pd_1 syme plan px_0 pr_0 (px_0 plus e_y) ;
  37. elim 1.e-4 (pd_2 et s_d2) ;
  38. s_d3 = s_d et (inve s_d2) ; s_b3 = s_b et (inve s_b2) ; elim 1.e-4 s_b3 ;
  39.  s_2 = s_d2 et s_b2 ; s_2 = s_2 coul vert ;
  40.  s_1 = s_0 ;
  41. s_0 = s_d3 et s_b3 ;
  42. 'SI' graph ;
  43. titre ' maillage 2D total ';
  44. trac (s_0 et px_0 et px_1) 'QUAL' nclk ;
  45. 'FINSI' ;
  46. ********************************
  47. E_1 = 200.d9 ; nu_1 = 0.3d0 ; rho_1 = 7.9e3 ; cf_1 = 2./3 ;
  48. * 10 m de profondeur
  49. pres0 = 1.e5 ;
  50. ******************************************************
  51. *---------- CONSTRUCTION BASE MODES de FLAMBAGE / FOURIER / 3D ------------------
  52. t23F = table 'BASE_MODALE' ; t23F . modes = table 'BASE_DE_MODES' ;
  53.  
  54. opti mode four noharm ;
  55.  
  56. mod_b = mode s_b3 mecanique elastique isotrope plastique isotrope coq2 ;
  57. mod_d = mode s_d3 mecanique elastique isotrope plastique isotrope coq2 ;
  58.  
  59. mod_t = mod_b et mod_d ;
  60.  
  61. ev_ecro = evol vert manu ep ((prog 0. 10.)*1.e-2) sig ((prog 400. 800.)*1.e6) ;
  62. ev_ecro = evol vert manu ep ((prog 0. 10.)*1.e-2) sig ((prog 100. 800.)*1.e6) ;
  63. car_b = mate mod_b youn E_1 nu nu_1 epai e_b rho rho_1 ecrou ev_ecro calf 0.7;
  64. car_d = mate mod_d youn E_1 nu nu_1 epai (e_b*1.5) rho rho_1 ecrou ev_ecro calf 0.7 ;
  65. car_t = car_b et car_d ;
  66.  
  67. *-------- CONDITIONS AUX LIMITES  ----------
  68. bl_0 = bloque pr_0  'UZ' 'UT' ; 
  69. blr_0z =  bloque pr_0  'UZ' ;
  70. bls_0 = symt s_0 'DEPL' 'ROTA'  px_0 pr_0 (px_0 plus e_y);
  71. bl_d3 = bloque (pd_1 et pd_2) 'UR' 'RT' ; 
  72. ********************************
  73.  
  74. npmo = 2 ;
  75. ****************************************************
  76. n_four = 0 ;
  77. opti mode four n_four ;
  78. mess '***** FOURIER ***** n :' n_four ' **********' ;
  79. rig_n = rigi mod_t car_t ;
  80. ** solution statique
  81. ** pression 
  82. mpr_bd = mode s_0 chargement pression coq2 ;
  83. **
  84. cpr_bd = pression mpr_bd 'PRES' pres0;
  85. fpr_bd = bsig mpr_bd cpr_bd ;
  86.  
  87. ri_0 = rig_n ;
  88. dep_m0 = reso (ri_0 et bl_d3 et bl_0) fpr_bd ;
  89. 'SI' graph ;
  90. titre ' deplacement radial (m)  ' ;
  91. trac (exco dep_m0 ur) (s_0 et px_0 et px_1) ;
  92. titre ' deplacement axial (m)  ' ;
  93. trac (exco dep_m0 uz) (s_0 et px_0 et px_1) ;
  94. 'FINSI' ;
  95. sig_m0 = sigma dep_m0 mod_t car_t ;
  96. 'SI' graph ;
  97. titre 'FOURIER 0 - contraintes' ;
  98. trac sig_m0 mod_t car_t ;
  99. 'FINSI' ;
  100.  
  101. e_m0 = epsi dep_m0 mod_t car_t ;
  102. s_m0 = elas e_m0 mod_t car_t ;
  103.  
  104. n_four = 2 ;
  105. opti mode four n_four ;
  106. mess '***** FOURIER ***** n :' n_four ' **********' ;
  107.    rig_n = rigi mod_t car_t ;
  108.    MKSI= KSIGMA mod_t car_t sig_m0 FLAM;
  109.    MKP = KP (manu chpo s_0 1 P (pres0 * -1.)) mod_t FLAM ;
  110.  
  111.   t_flamb0 = table ;
  112.   t_flamb0 . 'CLIM' = bl_d3 et bl_0 ;
  113.   t_flamb0 . 'OBJM' = mod_t         ;
  114.   t_flamb0 . 'LAM1' = 0.001         ;
  115.   t_flamb0 . 'LAM2' = 10000.        ;
  116.   t_flamb0 . 'NMOD' = npmo          ;
  117.   t_flamb0 . 'RIGI' = rig_n         ;
  118.   t_flamb0 . 'KSIG' = MKSI ET MKP   ;
  119.   t_fo = flambage t_flamb0 ;
  120.  
  121.   defla0 = (defo t_fo . 1 . depl s_0 0.) coul bleu ;
  122.   defla1 = (defo t_fo . 1 . depl s_0 1.5) coul oran ;
  123. titre 'MODE FOURIER #' 2 ' ORDRE #' t_fo . 1 . modm ' MULTIPLICATEUR ' t_fo . 1 . lamb ;
  124.   si graph ; trac (defla0 et defla1) ; finsi ;
  125.  
  126.  err1 = (7.23260E+00 - t_fo . 1 . lamb)  / 7.23260E+00 ; err1 = abs err1 ;
  127. ****************************************************
  128. * ou bien
  129. *
  130. *   MKSI= KSIGMA mod_t car_t  (sig_m0 * -1.d0) FLAM;
  131. *   MKP = KP (manu chpo s_0 1 P pres0) mod_t FLAM ;
  132. * 
  133. *
  134. ** t2   = VIBR 'PROCH' (PROG 0.) (rig_n ET bl_d3 et bl_0) MKSI;
  135. * t2   = VIBR 'INTERVALLE' 0.001 1.e2 'BASSE' npmo (rig_n ET bl_d3 et bl_0) (MKSI et MKP);
  136. * npm2 = dime t2 . modes - 2 ; 
  137. * it2 = 1 ;
  138. * def0 = (defo t2 . modes . it2 . deformee_modale s_0 0.) coul jaune ;
  139. * def2 = defo t2 . modes . it2 . deformee_modale s_0 1.5 ;
  140. * w1 = t2 . modes . it2 . frequence ; 
  141. *
  142. * LAMBDA1=(W1 * PI * 2) ** 2 * (SIGN W1);
  143. * 'SI' graph ;
  144. * titre 'ordre Fourier = ' 2 ' - mode = ' it2 ' - lambda1 = '  lambda1;
  145. * trac (def0 et (def2 coul vert)) ;
  146. * 'FINSI' ;
  147. ****************************************************
  148. si  (non (existe t23F . modes 'MAILLAGE')) ; 
  149.   mai_3d = madef3d s_0 t_fo . 1 . depl n_four 'MAILLAGE' ;
  150. *  mai_3d = mai_3d oriente point px_0 ;
  151.   elim 1.e-4 (mai_3d et s_0) ; 
  152.  
  153.   t23F . modes . maillage = mai_3d ;
  154.   opti mode trid ;
  155.  'SI' graph ;
  156.   titre 'maillage 3D noeuds = ' (nbno mai_3d) ' elements = ' (nbel mai_3d) ;
  157.   trac mai_3d nclk ; 
  158.  'FINSI' ;
  159. finsi ; 
  160.  
  161. it3d = 0 ;
  162.  repeter b_dm npmo ;
  163.  
  164.  it3d = it3d + 1 ;
  165.  t23F . modes . it3d = table ;
  166.  t23F . modes . it3d . deformee_modale = madef3d s_0 t_fo . it3d . depl n_four 'DEFORMEE' ;
  167.     mai_3d = t23F . modes . maillage ;
  168.  elim 1.e-4 (mai_3d et (extr t23F . modes . it3d . deformee_modale mail)) ;
  169.  t23F . modes . it3d . lambda = lambda1 ;
  170.   opti mode trid ;
  171.   oeil_1 = 1000. 1000. 1000. ;     
  172.  def3d0 = (defo t23F . modes . it3d . deformee_modale mai_3d 0.) coul bleu ;
  173.  lambda1 = t_fo . it3d . lamb ;
  174.  def3d2 = defo t23F . modes . it3d . deformee_modale mai_3d 2. ;
  175.   'SI' graph ;
  176.   titre '3D ordre Fourier = ' n_four ' - numero mode = ' it3d ' - lambda1 = ' lambda1 ;
  177.      trac ((def3d2 coul CORA) et def3d0) oeil_1 ;
  178.      trac (def3d2 coul CORA) oeil_1 ;
  179. *     trac t23d . modes . it3d . deformee_modale mai_3d ;
  180.   'FINSI' ;
  181.  
  182.   fin b_dm ;
  183.  
  184. **********************************************************
  185.  
  186. l_10 = prog 0. pas 1. 20. ;
  187. ev_pres = evol manu 'TEMP' l_10 'PRES' l_10 ;
  188. * chargement pression
  189. ch_pres1 = char 'MECA'  ev_pres fpr_bd ;
  190.  
  191. opti mode four 0 ;
  192.  
  193.  
  194. * calcul COQ2 elasto-plastique
  195. ta_2f = table 'PASAPAS' ;
  196. ta_2f . modele =           mod_t ;
  197. ta_2f . caracteristiques = car_t ;
  198. ta_2f . blocages_mecaniques = bl_d3 et bl_0 ;
  199. ta_2f. chargement = ch_pres1 ;
  200. ta_2f . temps_calcules = l_10 ;
  201.  
  202. pasapas ta_2f ;
  203.  
  204. **
  205. n_four = 2 ;
  206. opti mode four n_four ;
  207. mess '***** FOURIER ***** n :' n_four ' **********' ;
  208.    rig_n = rigi mod_t car_t ;
  209.    MKSI= KSIGMA mod_t car_t ta_2f . contraintes . 10 FLAM ;
  210.    MKP = KP (manu chpo s_0 1 P (pres0 * -10.)) mod_t FLAM ;
  211.  
  212.   t_flamb1 = table ;
  213.   t_flamb1 . 'CLIM' = bl_d3 et bl_0 ;
  214.   t_flamb1 . 'OBJM' = mod_t         ;
  215.   t_flamb1 . 'LAM1' = 0.001         ;
  216.   t_flamb1 . 'LAM2' = 10000.        ;
  217.   t_flamb1 . 'NMOD' = npmo          ;
  218.   t_flamb1 . 'RIGI' = rig_n         ;
  219.   t_flamb1 . 'KSIG' = MKSI ET MKP   ;
  220.   t_fo1 = flambage t_flamb1 ;
  221.  
  222.   defla10 = (defo t_fo1 . 1 . depl s_0 0.) coul bleu ;
  223.   defla11 = (defo t_fo1 . 1 . depl s_0 1.5) coul oran ;
  224. titre 'MODE FOURIER #' 2 ' ORDRE #' t_fo1 . 1 . modm ' MULTIPLICATEUR ' t_fo1 . 1 . lamb ;
  225.   si graph ; trac (defla10 et defla11) ; finsi ;
  226.  
  227.  err2 = (7.23260E-01 - t_fo1 . 1 . lamb) / 7.23260E-01 ; err2 = abs err2 ;
  228. **
  229. n_four = 2 ;
  230. opti mode four n_four ;
  231. mess '***** FOURIER ***** n :' n_four ' **********' ;
  232.    car_13 = mota mod_t ta_2f . contraintes . 13 ta_2F . variables_internes . 13 car_t ;
  233.    car_13 = (nomc car_13 'YOUN') et (exco car_t nu) et (exco car_t epai) ;
  234.    rig_n = rigi mod_t car_13 ;
  235.    MKSI= KSIGMA mod_t car_t ta_2f . contraintes . 13 FLAM ;
  236.    MKP = KP (manu chpo s_0 1 P (pres0 * -13.)) mod_t FLAM ;
  237.  
  238.   t_flam13 = table ;
  239.   t_flam13 . 'CLIM' = bl_d3 et bl_0 ;
  240.   t_flam13 . 'OBJM' = mod_t         ;
  241.   t_flam13 . 'LAM1' = 0.001         ;
  242.   t_flam13 . 'LAM2' = 10000.        ;
  243.   t_flam13 . 'NMOD' = npmo          ;
  244.   t_flam13 . 'RIGI' = rig_n         ;
  245.   t_flam13 . 'KSIG' = MKSI ET MKP   ;
  246.   t_fo13 = flambage t_flam13 ;
  247.  
  248.   defl130 = (defo t_fo13 . 1 . depl s_0 0.) coul bleu ;
  249.   defl131 = (defo t_fo13 . 1 . depl s_0 1.5) coul oran ;
  250. titre 'MODE FOURIER #' 2 ' ORDRE #' t_fo13 . 1 . modm ' MULTIPLICATEUR ' t_fo13 . 1 . lamb ;
  251.   si graph ; trac (defl130 et defl131) ; finsi ;
  252.  
  253.  
  254.  err3 = (1.94492E-02 - t_fo13 . 1 . lamb) / 1.94492E-02 ; err3 = abs err3 ;
  255.  
  256.  si ((err1 > 0.01) ou (err2 > 0.01) ou (err3 > 0.01)) ; erre 5 ; sinon ; erre 0 ; finsi ;
  257.  
  258. fin;
  259.  
  260.  
  261.  
  262.  

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