Télécharger fluage_maxwell_1.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : fluage_maxwell_1.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Test fluage_maxwell_1.dgibi: Jeux de données *
  5. * -------------------------------------------- *
  6. * *
  7. *
  8. **************************************************
  9. * *
  10. * TEST DE VERIFICATION D'UNE LOI DE COMPORTEMENT*
  11. * -------------------------------------------- *
  12. * DE MATERIAU *
  13. * ----------- *
  14. * *
  15. * LOI DE COMPORTEMENT DE FLUAGE DE : *
  16. * MAXWELL *
  17. * *
  18. * TEST POUR DES ELEMENTS MASSIFS *
  19. * ------------------------------ *
  20. * *
  21. * MAILLAGE: *
  22. * UNE PLAQUE CARREE DE COTE L=1 M *
  23. * *
  24. * CHARGEMENT: *
  25. * CISAILLEMENT PUR *
  26. * FORCES IMPOSEES *
  27. * DEFORMATIONS PLANES *
  28. * *
  29. **************************************************
  30.  
  31. opti dime 2 elem qua4 ;
  32.  
  33. *
  34. * Maillage
  35. *
  36. p1 = 0. 0. ;
  37. p2 = 1. 0. ;
  38. p3 = 1. 1. ;
  39. p4 = 0. 1. ;
  40. li1 = d p1 1 p2 ;
  41. li2 = d p2 1 p3 ;
  42. li3 = d p3 1 p4 ;
  43. li4 = d p4 1 p1 ;
  44. mail1 = 'SURF' ( li1 'ET' li2 'ET' li3 'ET' li4 ) 'PLAN' ;
  45. *
  46. * Modele
  47. *
  48. AF1 = 1.E-45 ;
  49. AF2 = 5. ;
  50. AF3 = 1.2 ;
  51. SMAX = 2.1E8;
  52. youn0 = 1.5E+4;
  53. nu0 = 0.3 ;
  54.  
  55. prog1=prog 0. 1.D12;
  56. prog20=prog 1.5E+4 1.5E+4;
  57. E0=evol manu prog1 prog20;
  58.  
  59. *prog1=prog 0. 1.D12;
  60. prog21=prog 1.11E+6 1.11E+6;
  61. E1=evol manu prog1 prog21;
  62.  
  63. *prog1=prog 0. 1.D12;
  64. prog22=prog 1.12E+6 1.12E+6;
  65. E2=evol manu prog1 prog22;
  66.  
  67. *prog1=prog 0. 1.D12;
  68. prog23=prog 3.43E+5 3.43E+5;
  69. E3=evol manu prog1 prog23;
  70.  
  71. *prog1=prog 0. 1.D12;
  72. prog24=prog 4.14E+5 4.14E+5;
  73. E4=evol manu prog1 prog24;
  74.  
  75. proyou= prog20 + prog21 + prog22 + prog23 + prog24 ;
  76. youn0 = evol manu 'TEMP' prog1 'YOUN' proyou;
  77.  
  78. TRR1 = 0.46;
  79. TRR2 = 0.45;
  80. TRR3 = 3.62;
  81. TRR4 = 45.1;
  82.  
  83. mod0 =MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'FLUAGE' 'MAXWELL';
  84. mat0 =MATE mod0 'YOUN' youn0 'NU' nu0
  85. 'EM0' E0 'EM1' E1 'EM2' E2 'EM3' E3 'EM4' E4
  86. 'TR1' TRR1 'TR2' TRR2 'TR3' TRR3 'TR4' TRR4;
  87.  
  88. *
  89. * Conditions aux limites
  90. *
  91. cl1 = 'BLOQ' mail1 'UY' ;
  92. cl2 = 'BLOQ' li1 'UX' ;
  93. cl0 = cl1 'ET' cl2 ;
  94. *
  95. * Chargement
  96. *
  97. depmax = 7.E-1 ;
  98. *chp1 = 'COOR' 2 mail1 ;
  99. *chp1 = chp1 * depmax / 1. ;
  100. *chp1 = 'NOMC' chp1 'UX' ;
  101. *titr ' Champ de forces impose' ;
  102. **trace chp1 mail1 ;
  103. *chp2 = 'MANU' 'CHPO' mail1 1 'UY' 0. ;
  104. *chp0 = chp1 + chp2 ;
  105. dep1 = 'FORCE' li1 'FX' ( -1. * depmax ) 'FY' 0. ;
  106. dep2 = 'FORCE' li3 'FX' depmax 'FY' 0. ;
  107. dep0 = dep1 'ET' dep2 ;
  108. ev0 = 'EVOL' 'MANU' temps ( 'PROG' 0. 1.E-4 1000. )
  109. y ( 'PROG' 0. 1. 1. ) ;
  110. cha0 = 'CHAR' 'MECA' dep0 ev0 ;
  111. *
  112. dt0 = 0.020 ;
  113. dt1 = 2. * dt0 ;
  114. dt2 = 5. * dt1 ;
  115. dt3 = 5. * dt2 ;
  116. tfin0 = 20. ;
  117. *
  118. * Appel à PASAPAS
  119. *
  120. ta1 = 'TABLE' ;
  121. ta1 .'MODELE' = mod0 ;
  122. ta1 .'CARACTERISTIQUES' = mat0 ;
  123. ta1 .'CHARGEMENT' = cha0 ;
  124. ta1 .'BLOCAGES_MECANIQUES' = cl0 ;
  125. ta1 .'TEMPS_CALCULES' = ( 'PROG' 0. 'PAS' 1.E-4 2.E-4 ) 'ET'
  126. ( 'PROG' 3.E-4 dt0 ) 'ET'
  127. ( 'PROG' dt1 'PAS' dt0 dt2 ) 'ET'
  128. ( 'PROG' ( dt2 + dt2 ) 'PAS' dt2 dt3 ) 'ET'
  129. ( 'PROG' ( dt3 + dt3 ) 'PAS' dt3 tfin0 ) ;
  130. *
  131. list ta1.'TEMPS_CALCULES';
  132. TMASAU=table;
  133. ta1. 'MES_SAUVEGARDES'=TMASAU;
  134. TMASAU.'DEFTO'=VRAI;
  135. TMASAU.'DEFIN'=VRAI;
  136. *
  137. PASAPAS ta1;
  138.  
  139. FIN;
  140.  
  141.  
  142.  

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