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. opti surv 2125600 ;
  7.  
  8. p0 = 0. 0. 0. ;
  9. e1 = 1. 0. 0. ; e2 = 0. 1. 0. ; e3 = 0. 0. 1. ;
  10.  
  11. l1 = d 1 P0 (p0 plus e1) ;
  12. *l2 = l1 d 2 (p0 plus e2) d 3 p0 ;
  13. s1 = l1 tran 1 e2 ;
  14. *s1 = surf plan l2 ;
  15. v1 = s1 volu 1 tran e3 ;
  16. s11 = changer s1 tri3 ;
  17. v11 = s11 volu 1 tran e3 ;
  18. v12 = v1 et (v1 plus e1) ;
  19.  
  20. mth1 = mode v1 thermique isotrope ;
  21.  
  22. me1 = mode v1 mecanique elastique plastique parfait phase pha1;
  23. me2 = mode v1 mecanique elastique fluage norton phase pha2;
  24. me5 = mode v1 mecanique elastique plastique parfait phase nu ;
  25.  
  26. me11 = mode v11 mecanique elastique plastique parfait phase pha1;
  27. me12 = mode v12 mecanique elastique plastique parfait phase pha1;
  28.  
  29. *mm1 = mode v1 melange ztmax serie (me1 et me2) ;
  30. mm2 = mode v1 melange ztmax parallele (me1 et me2) ;
  31. mm3 = mode v1 melange parallele (me1 et me2) ;
  32. mm4 = mode v1 melange parallele me1 ;
  33. mm5 = mode v1 melange parallele me5 ;
  34.  
  35.  
  36.  
  37. v2 = v1 plus e2 ;
  38. ms1 = mode v2 mecanique elastique plastique parfait phase pha1
  39. statio me1 ;
  40. ms2 = mode v2 mecanique elastique fluage norton phase pha2
  41. statio me2 ;
  42. mms2 = mode v2 melange ztmax parallele (ms1 et ms2) statio mm2 ;
  43. *
  44. ca1= mate me1 youn 2.e11 nu 0.3 sigy 200.e6 ;
  45. ca5= mate me5 youn 2.e11 nu 0.3 sigy 200.e6 ;
  46. ca2 = mate me2 young 150.e9 nu 0.3 rho 7.8e9 alpha 1.5e-5
  47. af1 ((exp (-8.56 * (log 10.)))*(1.e-6 ** 5.9)*1.e-6/3600.)
  48. af2 5.90 af3 1. smax 0.e6 ;
  49.  
  50. *cam2 = mate mm2 pha1 1. pha2 1. ;
  51. cam2 = mate mm2 pha1 0.8 pha2 0.2 AC1 600. AC2 1000.
  52. 'AC3' -200. 'AC4' -500. ;
  53. cam20 = mate mm2 pha1 0. pha2 1. VIPH 0. AC1 600. AC2 1000. ;
  54. cam3 = mate mm3 pha1 0.5 pha2 -.6 ;
  55. cam4 = mate mm4 pha1 0.5 ;
  56.  
  57. ri3 = rigi mm3 (cam3 et ca1 et ca2) ;
  58. ri4 = rigi mm4 (cam4 et ca1 ) ;
  59. ri5 = rigi mm5 ca5 ;
  60.  
  61.  
  62. chp1 = manu chpo s1 3 ux 0. uy 0. uz -0.01 nature diffuse ;
  63. chg1 = epsi (me1 et me2) chp1 ;
  64.  
  65. cho3 = comp mm3 (chg1 * 0.) (chg1 et cam3 et ca1 et ca2) ;
  66. sig2 = exco cho3 (extr me1 contraintes) ;
  67. sig3 = redu stric sig2 mm3 ;
  68. sig5 = redu sig2 (me1 et me2) ;
  69. fi3 = bsig mm3 sig3 ;
  70. fi4 = bsig (extr mm3 form mecanique) sig2 ;
  71.  
  72. * CB215821 : fi5 ne peut pas fonctionner comme ca, komcha ne trouve pas les contraintes du MODELE de melange
  73. *fi5 = bsig mm3 sig5 cam3 ;
  74.  
  75. chg41 = epsi me1 chp1 ;
  76. cho43 = comp mm4 (chg41 * 0.) (chg41 et cam4 et ca1) ;
  77. sig42 = exco cho43 (extr me1 contraintes) ;
  78. sig43 = redu stric sig42 mm4 ;
  79. fi43 = bsig mm4 sig43 ;
  80. fi44 = bsig (extr mm4 form mecanique) sig42 ;
  81.  
  82. chg51 = epsi me5 chp1 ;
  83. cho53 = comp mm5 (chg51 * 0.) (chg51 et ca5) ;
  84. sig52 = exco cho53 (extr me5 contraintes) ;
  85. sig53 = redu stric sig52 mm5 ;
  86. fi53 = bsig mm5 sig53 ;
  87. fi54 = bsig (extr mm5 form mecanique) sig52 ;
  88.  
  89.  
  90. cht0 = manu chml v1 'T' 100. ;
  91. cht1 = manu chml v1 'T' -800. ;
  92.  
  93. cho4 = comp mm2 (cht0 et cam20) (cam2 et cht1 et ca1 et ca2);
  94.  
  95.  
  96. *err1 = abs ((mini fi3) - (mini fi5)) ;
  97. *
  98. *si (err1 < (1.e-6*(abs (mini fi3))) ) ;
  99. * erre 0 ;
  100. *sinon ;
  101. * erre 5 ;
  102. *finsi ;
  103.  
  104.  
  105. fin ;
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  

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