Télécharger la plaquette Cast3M


fr en

Télécharger ouglova_1D.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : ouglova_1D.dgibi
  2. *******************************************************
  3. *******************************************************
  4. *******************************************************
  5. *  Cas test de l'implantation numerique du modele
  6. *  OUGLOVA 1D
  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 SEG2 TRAC PSC
  23. DENSITE 1;
  24.  
  25. ******************************************************
  26. *** Parametres geometriques
  27. ******************************************************
  28.  
  29. Ra  = 1;
  30. sc1 = 3.141592/4;
  31. h   = 1; 
  32.  
  33. ******************************************************
  34. *** Definition de la geometrie 
  35. ******************************************************
  36.  
  37. P1 = 0.  0.;
  38. P2 = h   0.;
  39.  
  40. LT = D P1 P2;
  41.  
  42. ******************************************************
  43. *** Définition du modele
  44. ******************************************************
  45.  
  46. MOD1 = MODE LT MECANIQUE ELASTIQUE ISOTROPE 
  47.                          PLASTIQUE OUGLOVA BARRE;
  48. *MOD1 = MODE LT MECANIQUE ELASTIQUE ISOTROPE BARRE;
  49. ******************************************************
  50. *** Parametres materiaux
  51. ****************************************************** 
  52.  
  53. * Module d'Young
  54. ygn = 200.e9;
  55.  
  56. * Coefficient de Poisson
  57. nun = 0.33;
  58.  
  59. * Limite elastique
  60. sigym = 500.e6;
  61.  
  62. * Coefficient d'ecrouissage
  63. km = 500.e6;
  64.  
  65. * Exposant d'ecrouissage
  66. mm = 2.786;
  67.  
  68. * Taux de corrosion en percentage
  69. tcm = 15;
  70.  
  71. * Endommagement critique
  72. dcm = 0.2;
  73.  
  74.  
  75. MAT1 = MATE MOD1 'YOUN' ygn 'NU' nun 'SECT' sc1
  76.                  'SIGY' sigym 'K   ' km 'M   ' mm
  77.                  'TC  ' tcm 'DC  ' dcm;
  78. *MAT1 = MATE MOD1 YOUN ygn NU nun SECT sc1;  
  79.  
  80. ******************************************************
  81. *** Definition des conditions aux limites
  82. *** et des déplacements imposés
  83. ******************************************************
  84.  
  85. CL1 = BLOQ P1 UX UY;
  86. CL2 = BLOQ LT UY;
  87. CL11 = BLOQ P2 UX;
  88. D1   = DEPI CL11 1;
  89.  
  90. CLTOT = CL1 ET CL2 ET CL11;
  91.  
  92. ******************************************************
  93. *** Definition de chargement
  94. ******************************************************
  95.  
  96. LI1  = PROG 0. 1.;
  97. LI2  = PROG 0. 2.5e-1;
  98. LIS1 = PROG 0. PAS 0.01 1;
  99.  
  100. EV   = EVOL MANU LI1 LI2;
  101.  
  102. CHA1 = CHAR 'DIMP' D1 EV;
  103.  
  104. CHTOT = CHA1;
  105.  
  106. ******************************************************
  107. *** Resolution
  108. ******************************************************
  109.  
  110. TAB1 = TABLE ;
  111. TAB1.'BLOCAGES_MECANIQUES' = CLTOT;
  112. TAB1.'MODELE' = MOD1;
  113. TAB1.'CHARGEMENT' = CHA1;
  114. TAB1.'CARACTERISTIQUES' = MAT1;
  115. TAB1.'TEMPS_CALCULES' = LIS1;
  116. TAB1.'MAXITERATION' = 1000;
  117. TAB1.'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE' ;
  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.  
  165.  
  166.  

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