Télécharger ouglova_1D.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : ouglova_1D.dgibi
  2. opti epsi lineaire ;
  3. *******************************************************
  4. *******************************************************
  5. *******************************************************
  6. * Cas test de l'implantation numerique du modele
  7. * OUGLOVA 1D
  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 SEG2 TRAC PSC
  24. DENSITE 1;
  25.  
  26. ******************************************************
  27. *** Parametres geometriques
  28. ******************************************************
  29.  
  30. Ra = 1;
  31. sc1 = 3.141592/4;
  32. h = 1;
  33.  
  34. ******************************************************
  35. *** Definition de la geometrie
  36. ******************************************************
  37.  
  38. P1 = 0. 0.;
  39. P2 = h 0.;
  40.  
  41. LT = D P1 P2;
  42.  
  43. ******************************************************
  44. *** Définition du modele
  45. ******************************************************
  46.  
  47. MOD1 = MODE LT MECANIQUE ELASTIQUE ISOTROPE
  48. PLASTIQUE OUGLOVA BARRE;
  49. *MOD1 = MODE LT MECANIQUE ELASTIQUE ISOTROPE BARRE;
  50. ******************************************************
  51. *** Parametres materiaux
  52. ******************************************************
  53.  
  54. * Module d'Young
  55. ygn = 200.e9;
  56.  
  57. * Coefficient de Poisson
  58. nun = 0.33;
  59.  
  60. * Limite elastique
  61. sigym = 500.e6;
  62.  
  63. * Coefficient d'ecrouissage
  64. km = 500.e6;
  65.  
  66. * Exposant d'ecrouissage
  67. mm = 2.786;
  68.  
  69. * Taux de corrosion en percentage
  70. tcm = 15;
  71.  
  72. * Endommagement critique
  73. dcm = 0.2;
  74.  
  75.  
  76. MAT1 = MATE MOD1 'YOUN' ygn 'NU' nun 'SECT' sc1
  77. 'SIGY' sigym 'K ' km 'M ' mm
  78. 'TC ' tcm 'DC ' dcm;
  79. *MAT1 = MATE MOD1 YOUN ygn NU nun SECT sc1;
  80.  
  81. ******************************************************
  82. *** Definition des conditions aux limites
  83. *** et des déplacements imposés
  84. ******************************************************
  85.  
  86. CL1 = BLOQ P1 UX UY;
  87. CL2 = BLOQ LT UY;
  88. CL11 = BLOQ P2 UX;
  89. D1 = DEPI CL11 1;
  90.  
  91. CLTOT = CL1 ET CL2 ET CL11;
  92.  
  93. ******************************************************
  94. *** Definition de chargement
  95. ******************************************************
  96.  
  97. LI1 = PROG 0. 1.;
  98. LI2 = PROG 0. 2.5e-1;
  99. LIS1 = PROG 0. PAS 0.01 1;
  100.  
  101. EV = EVOL MANU LI1 LI2;
  102.  
  103. CHA1 = CHAR 'DIMP' D1 EV;
  104.  
  105. CHTOT = CHA1;
  106.  
  107. ******************************************************
  108. *** Resolution
  109. ******************************************************
  110.  
  111. TAB1 = TABLE ;
  112. TAB1.'BLOCAGES_MECANIQUES' = CLTOT;
  113. TAB1.'MODELE' = MOD1;
  114. TAB1.'CHARGEMENT' = CHA1;
  115. TAB1.'CARACTERISTIQUES' = MAT1;
  116. TAB1.'TEMPS_CALCULES' = LIS1;
  117. TAB1.'MAXITERATION' = 1000;
  118.  
  119. PASAPAS TAB1;
  120.  
  121. ******************************************************
  122. *** Verification
  123. ******************************************************
  124.  
  125. EV2 = @GLOBAL TAB1 cl1 EV FX;
  126. liste EV2;
  127. dess ev2;
  128.  
  129. vref = EXTR EV2 'ORDO';
  130. vref = EXTR 29 vref;
  131.  
  132. err = abs(vref + 5.68172e8) / 5.68172e8;
  133. message ' erreur relative ';
  134. list err;
  135.  
  136. si (err > 1.e-3);
  137. erreur (5);
  138. finsi;
  139.  
  140. FIN;
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  

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