Télécharger la plaquette Cast3M


fr en

Télécharger melange.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  melange.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. opti dime 3 elem cub8 echo 1 ;
  5.  
  6. p0 = 0. 0. 0. ; 
  7. e1 = 1. 0. 0. ; e2 = 0. 1. 0. ; e3 = 0. 0. 1. ; 
  8.  
  9. l1 = d 1 P0  (p0 plus e1) ;
  10. *l2 = l1 d 2 (p0 plus e2) d 3 p0 ;
  11. s1 = l1 tran 1 e2 ;
  12. *s1 = surf plan l2 ;
  13. v1 = s1 volu 1 tran e3 ;
  14. s11 = changer s1 tri3 ;
  15. v11 = s11 volu 1 tran e3 ;
  16. v12 = v1 et (v1 plus e1) ;
  17.  
  18. mth1 = mode v1 thermique isotrope ;
  19.  
  20.  me1 = mode v1 mecanique elastique plastique parfait phase pha1;
  21.  me2 = mode v1 mecanique elastique fluage norton phase pha2; 
  22.  
  23. mefa1 =  mode v1 mecanique elastique plastique parfait phase fa1;
  24. mefa2 =  mode v1 mecanique elastique plastique parfait phase fa2;
  25.  
  26. me12 =  mode v12 mecanique elastique plastique parfait phase fa1;
  27.  
  28. * cohérence noms des phases mecaniques et noms phases modele ztmax
  29. mm1 = mode v1 melange ztmax ;
  30.  
  31. mm2 = mode v1 melange parallele (me1 et me2) ;
  32.  
  33. mmfa = mode v1 melange parallele (mefa1 et mefa2) ;
  34.  
  35. mm4 = mode v1 melange parallele me1  ;
  36.  
  37. *
  38. *
  39. v2 = v1 plus e2 ; v3 = v2 plus e2 ;
  40. ms1 =  mode v2 mecanique elastique plastique parfait phase pha1
  41.                   statio me1 ;
  42. ms2 =  mode v2 mecanique  elastique fluage norton phase pha2 
  43.                  statio me2 ;
  44.  
  45. mms1 = mode v2 melange ztmax statio mm1 ;
  46. mms2 = mode v2 melange parallele (ms1 et ms2) statio mm2 ;
  47. *
  48. * 
  49. t_mai = table 'MAILLAGE' ;
  50. t_mai . 1 = v1 ;
  51. t_mai . 2 = v2 ;
  52. t_mai . 3 = v3 ;
  53.  
  54. mstazt = mode t_mai melange statio mm1 ;
  55. mstapa = mode t_mai melange statio mm2 ;
  56.  
  57.  
  58. t_sta = table 'STATIONNAIRE' ;
  59. t_sta . maillage = table ;
  60. t_sta . maillage . 1  = v1 ;
  61. t_sta . maillage . 2 = v2 ;
  62.  
  63. mstae1 = mode t_sta mecanique elastique plastique parfait phase pha1 ;
  64. mstae2 = mode t_sta mecanique elastique fluage norton phase pha2; 
  65.  
  66. mstam1 = mode t_sta melange ztmax ;
  67. mstam2 = mode t_sta melange parallele (mstae1 et mstae2) ;
  68.  
  69. mstam11 = redu mstam1 v1 ; 
  70. mstam21 = redu mstam2 v1 ; 
  71.  
  72.  
  73. ca1= mate me1 youn 2.e11 nu 0.3 sigy 200.e6 ; 
  74. ca2 = mate me2 young 150.e9 nu 0.3 rho 7.8e9 alpha 1.5e-5
  75.   af1 ((exp (-8.56 * (log 10.)))*(1.e-6 ** 5.9)*1.e-6/3600.)
  76.   af2 5.90 af3 1.  smax 0.e6 ;
  77.  
  78. cam1 = mate mm1   AC1 300. AC2 750.  'AC3' 800. 'AC4' 250. 
  79.        VIPH 1000. VDEH 1000. VPAR 1000. ;
  80.  
  81. cafa1= mate mefa1 youn 2.e11 nu 0.3 sigy 200.e6 ; 
  82. cafa2= mate mefa2 youn 3.e11 nu 0.4 sigy 300.e6 ; 
  83.  
  84. cas1= mate ms1 youn 2.e11 nu 0.3 sigy 200.e6 ; 
  85. cas2 = mate ms2 young 150.e9 nu 0.3 rho 7.8e9 alpha 1.5e-5
  86.   af1 ((exp (-8.56 * (log 10.)))*(1.e-6 ** 5.9)*1.e-6/3600.)
  87.   af2 5.90 af3 1.  smax 0.e6 ;
  88.  
  89. casm1 = mate mms1   AC1 300. AC2 750.  'AC3' 800. 'AC4' 250. 
  90.        VIPH 1000. VDEH 1000. VPAR 1000. ;
  91.  
  92. castam11 = mate mstam11   AC1 300. AC2 750.  'AC3' 800. 'AC4' 250. 
  93.        VIPH 1000. VDEH 1000. VPAR 1000. ;
  94.  
  95. castae1 = mate mstae1 youn 2.e11 nu 0.3 sigy 200.e6 ;
  96. castae2 = mate mstae2 young 150.e9 nu 0.3 rho 7.8e9 alpha 1.5e-5
  97.   af1 ((exp (-8.56 * (log 10.)))*(1.e-6 ** 5.9)*1.e-6/3600.)
  98.   af2 5.90 af3 1.  smax 0.e6 ;
  99.  
  100. ************************************************************
  101.  
  102.  chpha21 = manu chml mm2 pha1 0.8 pha2 0.2 'RIGIDITE' ;
  103.  chpha20 = manu chml mm2 pha1 0. pha2 1. 'RIGIDITE' ;
  104.  chpha10 = manu chml mm1 pha1 0.1 pha2 .9 'RIGIDITE' ;
  105.  
  106.  
  107. chfa12 = manu chml mmfa  fa1 0.5 fa2 -.6 'RIGIDITE' ;
  108. chpha4 = manu chml mm4  pha1 0.5  ;
  109.  
  110. ri20 = rigi mm2 (chpha20 et ca1 et ca2) ;
  111. ri21 = rigi mm2 (chpha21 et ca1 et ca2) ;
  112. ri4 = rigi mm4 (chpha4 et ca1 ) ;
  113.  
  114. rifa1 = rigi mmfa (chfa12 et cafa1 et cafa2) ; 
  115.  
  116. chp1 = manu chpo s1 3 ux 0. uy 0.  uz -0.01 nature diffuse ;
  117.  chg1 = epsi (me1 et me2)  chp1 ;
  118.  
  119. cho3 = comp mm2 (chg1 * 0.) (chg1 et chpha21 et ca1 et ca2) ;
  120.  sig2 = exco cho3 (extr me1 contraintes) ;
  121.  sig3 = redu stric sig2 mm2 ;
  122.  sig5 = redu sig2 (me1 et me2) ;
  123.  
  124. fi3 = bsig mm2 sig3 ; 
  125. * fi4 = fi44
  126. fi4 = bsig (extr mm2 form mecanique) sig2 ;
  127. fi44 = bsig (extr mm2 form mecanique) sig5 ;
  128.  
  129. * CB215821 : fi5 ne peut pas fonctionner comme ca, komcha ne trouve pas les contraintes du MODELE de melange / oui  fi5 = fi3 
  130. fi5 = bsig mm2  sig5 chpha21 ;
  131. err1 = maxi (abs ( fi3 - fi5)) ; mess err1 ;
  132. si (err1 < 1.e-6) ;
  133.   erre 0 ; 
  134. sinon ;
  135.   erre 5 ;
  136. finsi ;
  137.  
  138.  chg41 = epsi me1  chp1 ;
  139. cho43 = comp mm4 (chg41 * 0.) (chg41 et chpha4 et ca1) ;
  140.  sig42 = exco cho43 (extr me1 contraintes) ;
  141.  sig43 = redu stric sig42 mm4 ;
  142. fi43 = bsig mm4 sig43 ; 
  143. fi44 = bsig (extr mm4 form mecanique) sig42 ;
  144.  
  145.  
  146.  
  147. cht0 = manu chml v1 'T' 1000. ;
  148. cht1 = manu chml v1 'T' 200. ;
  149.  
  150.  chg2 = epsi mm2  chp1 ;
  151.  
  152. cho4 = comp (mm1 et mm4) (cht0 et cam1 et chpha10 et (chg41*0.)) 
  153.     (cht1 et ca1 et ca2 et chg41);
  154.  
  155.  
  156. cht0 = manu chml v1 'T' 1000. ;
  157. cht1 = manu chml v1 'T' 350. ;
  158. cho5 = comp (mm1 et mm2) (cht0 et cam1 et chpha10 et (chg2*0.)) 
  159.     (cht1 et ca1 et ca2 et chg2);
  160.  
  161.  
  162. chpha5 = redu cho5 mm2 ;
  163. sigs5 = exco cho5 (extr me1 contraintes) ;
  164. sigs12 = redu sigs5 (me1 et me2) ;
  165. sig5m2 = redu stric sigs5 mm2 ;
  166. fi5m2 = bsig mm2 sig5m2 ;
  167. fi5s12 = bsig mm2 sigs12 chpha5 ;
  168.  
  169. err2 = maxi (abs ( fi5s12 - fi5m2)) ; mess err2 ;
  170. si (err1 < 1.e-6) ;
  171.   erre 0 ; 
  172. sinon ;
  173.   erre 5 ;
  174. finsi ;
  175. ************************************************************
  176.  
  177. chts0 = manu chml v2 'T' 1000. ;
  178. chts1 = manu chml v2 'T' 500. ;
  179.  
  180. s2 = face 2 v2 ;
  181. chp2 = manu chpo s2 3 ux 0. uy 0.  uz -0.02 nature diffuse ;
  182.  chgs2 = epsi (ms1 et ms2)  chp2 ;
  183.  
  184. cho6 = comp (mm1 et mm2 et mms1 et mms2)
  185.       (cht0 et cam1 et chpha10 et (chg2*0.))
  186.       (cht1 et chts1 et ca1 et cas1 et ca2 et cas2 et chg2 et chgs2); 
  187.  
  188. sig6m2s = redu stric (exco cho6 (extr me1 contraintes)) (mm2 et mms2) ;
  189. fi6m2s = bsig (mm2 et mms2) sig6m2s ;
  190.  
  191.  
  192.  chpha11 = manu chml mstam11 pha1 0.1 pha2 .9 'RIGIDITE' ;
  193.  chg3 = epsi mstam21  chp1 ;
  194.  chgs3 = epsi (redu mstam2 v2) chp2 ;
  195.  
  196. cho8 = comp (mstam1 et mstam2)
  197.       (cht0 et castam11 et chpha11 et (chg3*0.))
  198.       (cht1 et chts1 et castae1 et castae2 et chg3 et chgs3); 
  199. sig8m2 = redu stric (exco cho8 (extr me1 contraintes)) mstam2 ;
  200. fi8m2 = bsig mstam2 sig8m2 ;
  201.  
  202. err3 = maxi (abs ( fi6m2s - fi8m2)) ; mess err3 ;
  203. mx82 = maxi fi8m2 ; mess 'maxi fi8m2 ' mx82 ;
  204.  
  205.  
  206. si ((err3 < 1.e-6) et (mx82 > 1.e8)) ;
  207.   erre 0 ; 
  208. sinon ;
  209.   erre 5 ;
  210. finsi ;
  211.  
  212. fin ;
  213.  
  214.  
  215.  
  216.  
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  
  228.  

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