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