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

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