Télécharger ouglova_3D.dgibi

Retour à la liste

Numérotation des lignes :

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

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