Télécharger ouglova_CP.dgibi

Retour à la liste

Numérotation des lignes :

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

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