Télécharger la plaquette Cast3M


fr en

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. *  pas de controle de positivite
  137. ta1.'STABILITE' = faux;
  138. * 
  139. PASAPAS ta1; 
  140.  
  141. FIN;
  142.  
  143.  
  144.  
  145.  
  146.  

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