Télécharger la plaquette Cast3M


fr en

Télécharger fluage_maxwell_thve.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  fluage_maxwell_1.dgibi
  2. ************************************************************************
  3.  
  4. opti dime 2 elem qua8 ; 
  5.  
  6. * 
  7. * Maillage 
  8. * 
  9. p1 = 0. 0. ; 
  10. p2 = 1. 0. ; 
  11. p3 = 1. 1. ; 
  12. p4 = 0. 1. ; 
  13. li1 = d p1 1 p2 ; 
  14. li2 = d p2 1 p3 ; 
  15. li3 = d p3 1 p4 ; 
  16. li4 = d p4 1 p1 ; 
  17. mail1 = 'SURF' ( li1 'ET' li2 'ET' li3 'ET' li4 ) 'PLAN' ; 
  18. * 
  19. * Modele 
  20. * 
  21. AF1 = 1.E-45 ; 
  22. AF2 = 5. ; 
  23. AF3 = 1.2 ; 
  24. SMAX = 2.1E8; 
  25. youn0 = 1.5E+4; 
  26. nu0 = 0.3 ; 
  27.  
  28. prog1=prog 0.    1.D12; 
  29. prog20=prog 1.5E+4 1.5E+4; 
  30. E0=evol manu prog1 prog20; 
  31.  
  32. *prog1=prog 0. 1.D12; 
  33. prog21=prog 1.11E+6 1.11E+6; 
  34. E1=evol manu prog1 prog21; 
  35.  
  36. *prog1=prog 0.    1.D12; 
  37. prog22=prog 1.12E+6 1.12E+6; 
  38. E2=evol manu prog1 prog22; 
  39.  
  40. *prog1=prog 0. 1.D12; 
  41. prog23=prog 3.43E+5 3.43E+5; 
  42. E3=evol manu prog1 prog23; 
  43.  
  44. *prog1=prog 0. 1.D12; 
  45. prog24=prog 4.14E+5 4.14E+5; 
  46. E4=evol manu prog1 prog24; 
  47.  
  48. proyou= prog20 + prog21 + prog22 + prog23 + prog24 ;
  49. youn0 = evol manu 'TEMP' prog1 'YOUN' proyou; 
  50.  
  51. TRR1 = 0.46; 
  52. TRR2 = 0.45; 
  53. TRR3 = 3.62; 
  54. TRR4 = 45.1; 
  55.  
  56. mod0 =MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'FLUAGE' 'MAXWELL'; 
  57. mat0 =MATE mod0 'YOUN' youn0 'NU' nu0 'RHO' 2350. 'TREF' 15. 'TALP' 15.
  58. 'ALPHA' 29.e-6
  59. 'EM0' E0 'EM1' E1 'EM2' E2 'EM3' E3 'EM4' E4
  60. 'TR1' TRR1 'TR2' TRR2 'TR3' TRR3 'TR4' TRR4;
  61.  
  62. K1 = 2. ; C1 = 869. ; RHO1 = 2350. ;
  63. K1 = evol manu T (prog 0. 1.) K (prog 2. 2.) ;
  64. mod1 = MODE mail1 'THERMIQUE' 'ISOTROPE' 'CONDUCTION' ;
  65. mat1 = MATE mod1 'K' K1 'C' C1 'RHO' RHO1 'TINI' 15. ;
  66.  
  67. mod2 = mod0 ET mod1;
  68. mat2 = mat0 ET mat1;
  69.  
  70.  
  71. * Conditions aux limites 
  72. * 
  73. cl1 = 'BLOQ' mail1 'UY' ; 
  74. cl2 = 'BLOQ' li1 'UX' ; 
  75. cl0 = cl1 'ET' cl2 ; 
  76. * 
  77. * Chargement 
  78. * 
  79. depmax = 7.E-1 ; 
  80. *chp1 = 'COOR' 2 mail1 ; 
  81. *chp1 = chp1 * depmax / 1. ; 
  82. *chp1 = 'NOMC' chp1 'UX' ; 
  83. *titr ' Champ de forces impose' ; 
  84. **trace chp1 mail1 ; 
  85. *chp2 = 'MANU' 'CHPO' mail1 1 'UY' 0. ; 
  86. *chp0 = chp1 + chp2 ; 
  87. dep1 = 'FORCE' li1 'FX' ( -1. * depmax ) 'FY' 0. ; 
  88. dep2 = 'FORCE' li3 'FX' depmax 'FY' 0. ; 
  89. dep0 = dep1 'ET' dep2 ; 
  90. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1.E-4 1000. ) 
  91. y ( 'PROG' 0. 1. 1. ) ; 
  92. cha0 = 'CHAR' 'MECA' dep0 ev0 ; 
  93. * 
  94. dt0 = 0.020 ; 
  95. dt1 = 2. * dt0 ; 
  96. dt2 = 5. * dt1 ; 
  97. dt3 = 5. * dt2 ; 
  98. tfin0 = 20. ; 
  99. * 
  100. * Appel à PASAPAS 
  101. * 
  102. ta1 = 'TABLE' ; 
  103. ta1 .'MODELE' = mod2 ; 
  104. ta1 .'CARACTERISTIQUES' = mat2 ; 
  105. ta1 .'CHARGEMENT' = cha0 ; 
  106. ta1 .'BLOCAGES_MECANIQUES' = cl0 ; 
  107. ta1 .'TEMPS_CALCULES' = ( 'PROG' 0. 'PAS' 1.E-4 2.E-4 ) 'ET' 
  108. ( 'PROG' 3.E-4 dt0 ) 'ET' 
  109. ( 'PROG' dt1 'PAS' dt0 dt2 ) 'ET' 
  110. ( 'PROG' ( dt2 + dt2 ) 'PAS' dt2 dt3 ) 'ET' 
  111. ( 'PROG' ( dt3 + dt3 ) 'PAS' dt3 tfin0 ) ; 
  112. * 
  113. PASAPAS ta1; 
  114.  
  115. FIN;
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  

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