Télécharger ouglova_DP.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : ouglova_DP.dgibi
  2. opti epsi lineaire ;
  3. *******************************************************
  4. *******************************************************
  5. *******************************************************
  6. * Cas test de l'implantation numerique du modele
  7. * OUGLOVA 2D en deformations planes
  8. *******************************************************
  9. * Developpe par :
  10. * Romili PAREDES
  11. * Benjamin RICHARD
  12. * Contact :
  13. * Romili.Paredes@cea.fr
  14. * Benjamin.Richard@cea.fr
  15. * Institution :
  16. * CEA\DEN\DANS\DM2S\SEMT\EMSI
  17. ******************************************************
  18.  
  19. *******************************************************
  20. * Options
  21. *******************************************************
  22.  
  23. OPTI DIME 2 ELEM TRI3 TRAC PSC MODE PLAN DEFO;
  24. DENSITE 1;
  25.  
  26. ******************************************************
  27. *** Definition de la geometrie
  28. ******************************************************
  29.  
  30. P1 = 0. 0.;
  31. P2 = 1. 0.;
  32. P3 = 1. 1.;
  33. P4 = 0. 1.;
  34.  
  35. L1 = D P1 P2;
  36. L2 = D P2 P3;
  37. L3 = D P3 P4;
  38. L4 = D P4 P1;
  39.  
  40. ST = SURF 'PLAN' (L1 ET L2 ET L3 ET L4);
  41.  
  42. ******************************************************
  43. *** Définition du modele
  44. ******************************************************
  45.  
  46. MOD1 = MODE ST MECANIQUE ELASTIQUE ISOTROPE PLASTIQUE OUGLOVA;
  47.  
  48. ******************************************************
  49. *** Parametres materiaux
  50. ******************************************************
  51.  
  52. * Module d'Young
  53. ygn = 200.e9;
  54.  
  55. * Coefficient de Poisson
  56. nun = 0.33;
  57.  
  58. * Limite elastique
  59. sigym = 500.e6;
  60.  
  61. * Coefficient d'ecrouissage
  62. km = 500.e6;
  63.  
  64. * Exposant d'ecrouissage
  65. mm = 2.786;
  66.  
  67. * Taux de corrosion en percentage
  68. tcm = 15;
  69.  
  70. * Endommagement critique
  71. dcm = 0.2;
  72.  
  73.  
  74. MAT1 = MATE MOD1 'YOUN' ygn 'NU' nun
  75. 'SIGY' sigym 'K ' km 'M ' mm
  76. 'TC ' tcm 'DC ' dcm;
  77.  
  78. ******************************************************
  79. *** Definition des conditions aux limites
  80. *** et des déplacements imposés
  81. ******************************************************
  82.  
  83. CL1 = BLOQ L1 UY;
  84. CL2 = BLOQ P1 UX;
  85.  
  86. CL11 = BLOQ L3 UY;
  87. D1 = DEPI CL11 1;
  88.  
  89. CLTOT = CL1 ET CL2 ET CL11;
  90.  
  91. ******************************************************
  92. *** Definition de chargement
  93. ******************************************************
  94.  
  95. LI1 = PROG 0. 1.;
  96. LI2 = PROG 0. 3.5e-1;
  97. LIS1 = PROG 0. PAS 0.01 1;
  98.  
  99. EV = EVOL MANU LI1 LI2;
  100.  
  101. CHA1 = CHAR 'DIMP' D1 EV;
  102.  
  103. CHTOT = CHA1;
  104.  
  105. ******************************************************
  106. *** Resolution
  107. ******************************************************
  108.  
  109. TAB1 = TABLE ;
  110. TAB1.'BLOCAGES_MECANIQUES' = CLTOT;
  111. TAB1.'MODELE' = MOD1;
  112. TAB1.'CHARGEMENT' = CHA1;
  113. TAB1.'CARACTERISTIQUES' = MAT1;
  114. TAB1.'TEMPS_CALCULES' = LIS1;
  115. TAB1.'MAXITERATION' = 1000;
  116.  
  117. PASAPAS TAB1;
  118.  
  119. ******************************************************
  120. *** Verification
  121. ******************************************************
  122.  
  123. EV2 = @GLOBAL TAB1 cl1 EV FY;
  124. liste EV2;
  125.  
  126. vref = EXTR EV2 'ORDO';
  127. vref = EXTR 19 vref;
  128.  
  129. err = abs(vref + 6.59512e8) / 6.59512e8;
  130. message ' erreur relative ';
  131. list err;
  132.  
  133. si (err > 1.e-3);
  134. erreur (5);
  135. finsi;
  136.  
  137. FIN;
  138.  
  139.  
  140.  
  141.  
  142.  

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