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 ;
  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.  
  50. t_sta = table 'STATIONNAIRE' ;
  51. t_sta . maillage = table ;
  52. t_sta . maillage . 1  = v1 ;
  53. t_sta . maillage . 2 = v2 ;
  54.  
  55. mstae1 = mode t_sta mecanique elastique plastique parfait phase pha1;
  56. mstae2 = mode t_sta mecanique elastique fluage norton phase pha2; 
  57.  
  58. mstam1 = mode t_sta melange ztmax ;
  59. mstam2 = mode t_sta melange parallele (mstae1 et mstae2) ;
  60.  
  61. mstam11 = redu mstam1 v1 ;
  62. mstam21 = redu mstam2 v1 ;
  63.  
  64.  
  65. ca1= mate me1 youn 2.e11 nu 0.3 sigy 200.e6 ; 
  66. ca2 = mate me2 young 150.e9 nu 0.3 rho 7.8e9 alpha 1.5e-5
  67.   af1 ((exp (-8.56 * (log 10.)))*(1.e-6 ** 5.9)*1.e-6/3600.)
  68.   af2 5.90 af3 1.  smax 0.e6 ;
  69.  
  70. cam1 = mate mm1   AC1 300. AC2 750.  'AC3' 800. 'AC4' 250. 
  71.        VIPH 1000. VDEH 1000. VPAR 1000. ;
  72.  
  73. cafa1= mate mefa1 youn 2.e11 nu 0.3 sigy 200.e6 ; 
  74. cafa2= mate mefa2 youn 3.e11 nu 0.4 sigy 300.e6 ; 
  75.  
  76. cas1= mate ms1 youn 2.e11 nu 0.3 sigy 200.e6 ; 
  77. cas2 = mate ms2 young 150.e9 nu 0.3 rho 7.8e9 alpha 1.5e-5
  78.   af1 ((exp (-8.56 * (log 10.)))*(1.e-6 ** 5.9)*1.e-6/3600.)
  79.   af2 5.90 af3 1.  smax 0.e6 ;
  80.  
  81. casm1 = mate mms1   AC1 300. AC2 750.  'AC3' 800. 'AC4' 250. 
  82.        VIPH 1000. VDEH 1000. VPAR 1000. ;
  83.  
  84. castam1 = mate mstam1   AC1 300. AC2 750.  'AC3' 800. 'AC4' 250. 
  85.        VIPH 1000. VDEH 1000. VPAR 1000. ;
  86.  
  87. castae1 = mate mstae1 youn 2.e11 nu 0.3 sigy 200.e6 ;
  88. castae2 = mate mstae2 young 150.e9 nu 0.3 rho 7.8e9 alpha 1.5e-5
  89.   af1 ((exp (-8.56 * (log 10.)))*(1.e-6 ** 5.9)*1.e-6/3600.)
  90.   af2 5.90 af3 1.  smax 0.e6 ;
  91.  
  92. ************************************************************
  93.  
  94.  chpha21 = manu chml mm2 pha1 0.8 pha2 0.2 'RIGIDITE' ;
  95.  chpha20 = manu chml mm2 pha1 0. pha2 1. 'RIGIDITE' ;
  96.  chpha10 = manu chml mm1 pha1 0.1 pha2 .9 'RIGIDITE' ;
  97.  
  98.  
  99. chfa12 = manu chml mmfa  fa1 0.5 fa2 -.6 'RIGIDITE' ;
  100. chpha4 = manu chml mm4  pha1 0.5  ;
  101.  
  102. ri20 = rigi mm2 (chpha20 et ca1 et ca2) ;
  103. ri21 = rigi mm2 (chpha21 et ca1 et ca2) ;
  104. ri4 = rigi mm4 (chpha4 et ca1 ) ;
  105.  
  106. rifa1 = rigi mmfa (chfa12 et cafa1 et cafa2) ; 
  107.  
  108. chp1 = manu chpo s1 3 ux 0. uy 0.  uz -0.01 nature diffuse ;
  109.  chg1 = epsi (me1 et me2)  chp1 ;
  110.  
  111. cho3 = comp mm2 (chg1 * 0.) (chg1 et chpha21 et ca1 et ca2) ;
  112.  sig2 = exco cho3 (extr me1 contraintes) ;
  113.  sig3 = redu stric sig2 mm2 ;
  114.  sig5 = redu sig2 (me1 et me2) ;
  115.  
  116. fi3 = bsig mm2 sig3 ; 
  117. * fi4 = fi44
  118. fi4 = bsig (extr mm2 form mecanique) sig2 ;
  119. fi44 = bsig (extr mm2 form mecanique) sig5 ;
  120.  
  121. * CB215821 : fi5 ne peut pas fonctionner comme ca, komcha ne trouve pas les contraintes du MODELE de melange / oui  fi5 = fi3 
  122. fi5 = bsig mm2  sig5 chpha21 ;
  123. err1 = maxi (abs ( fi3 - fi5)) ; mess err1 ;
  124. si (err1 < 1.e-6) ;
  125.   erre 0 ; 
  126. sinon ;
  127.   erre 5 ;
  128. finsi ;
  129.  
  130.  chg41 = epsi me1  chp1 ;
  131. cho43 = comp mm4 (chg41 * 0.) (chg41 et chpha4 et ca1) ;
  132.  sig42 = exco cho43 (extr me1 contraintes) ;
  133.  sig43 = redu stric sig42 mm4 ;
  134. fi43 = bsig mm4 sig43 ; 
  135. fi44 = bsig (extr mm4 form mecanique) sig42 ;
  136.  
  137.  
  138.  
  139. cht0 = manu chml v1 'T' 1000. ;
  140. cht1 = manu chml v1 'T' 200. ;
  141.  
  142.  chg2 = epsi mm2  chp1 ;
  143.  
  144. cho4 = comp (mm1 et mm4) (cht0 et cam1 et chpha10 et (chg41*0.)) 
  145.     (cht1 et ca1 et ca2 et chg41);
  146.  
  147.  
  148. cht0 = manu chml v1 'T' 1000. ;
  149. cht1 = manu chml v1 'T' 350. ;
  150. cho5 = comp (mm1 et mm2) (cht0 et cam1 et chpha10 et (chg2*0.)) 
  151.     (cht1 et ca1 et ca2 et chg2);
  152.  
  153.  
  154. chpha5 = redu cho5 mm2 ;
  155. sigs5 = exco cho5 (extr me1 contraintes) ;
  156. sigs12 = redu sigs5 (me1 et me2) ;
  157. sig5m2 = redu stric sigs5 mm2 ;
  158. fi5m2 = bsig mm2 sig5m2 ;
  159. fi5s12 = bsig mm2 sigs12 chpha5 ;
  160.  
  161. err2 = maxi (abs ( fi5s12 - fi5m2)) ; mess err2 ;
  162. si (err1 < 1.e-6) ;
  163.   erre 0 ; 
  164. sinon ;
  165.   erre 5 ;
  166. finsi ;
  167. ************************************************************
  168.  
  169. chts0 = manu chml v2 'T' 1000. ;
  170. chts1 = manu chml v2 'T' 500. ;
  171.  
  172. s2 = face 2 v2 ;
  173. chp2 = manu chpo s2 3 ux 0. uy 0.  uz -0.02 nature diffuse ;
  174.  chgs2 = epsi (ms1 et ms2)  chp2 ;
  175.  
  176. cho6 = comp (mm1 et mm2 et mms1 et mms2)
  177.       (cht0 et cam1 et chpha10 et (chg2*0.))
  178.       (cht1 et chts1 et ca1 et cas1 et ca2 et cas2 et chg2 et chgs2); 
  179.  
  180. sig6m2s = redu stric (exco cho6 (extr me1 contraintes)) (mm2 et mms2) ;
  181. fi6m2s = bsig (mm2 et mms2) sig6m2s ;
  182.  
  183.  
  184.  chpha11 = manu chml mstam11 pha1 0.1 pha2 .9 'RIGIDITE' ;
  185.  chg3 = epsi mstam21  chp1 ;
  186.  chgs3 = epsi (redu mstam2 v2) chp2 ;
  187.  
  188. cho8 = comp (mstam1 et mstam2)
  189.       (cht0 et castam1 et chpha11 et (chg3*0.))
  190.       (cht1 et chts1 et castae1 et castae2 et chg3 et chgs3); 
  191. sig8m2 = redu stric (exco cho8 (extr me1 contraintes)) mstam2 ;
  192. fi8m2 = bsig mstam2 sig8m2 ;
  193.  
  194. err3 = maxi (abs ( fi6m2s - fi8m2)) ; mess err3 ;
  195. mx82 = maxi fi8m2 ; mess 'maxi fi8m2 ' mx82 ;
  196.  
  197. si ((err3 < 1.e-6) et (mx82 > 1.e8)) ;
  198.   erre 0 ; 
  199. sinon ;
  200.   erre 5 ;
  201. finsi ;
  202. *opti donn 5 ;
  203.  
  204. fin ;
  205.  
  206.  
  207.  
  208.  
  209.  
  210.  
  211.  
  212.  
  213.  
  214.  
  215.  
  216.  
  217.  
  218.  

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