Télécharger la plaquette Cast3M


fr en

Télécharger ouglova_3D.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : ouglova_3D.dgibi
  2. *******************************************************
  3. *******************************************************
  4. *******************************************************
  5. *  Cas test de l'implantation numerique du modele
  6. *  OUGLOVA 3D
  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 3 ELEM CUB8 TRAC PSC;
  23.  
  24. ******************************************************
  25. *** Parametres geometriques
  26. ******************************************************
  27.  
  28. b = 1.0;
  29. h = 1.0;
  30. e = 1.0;
  31.  
  32. nbe = 3;
  33. nbh = 3;
  34.  
  35. ******************************************************
  36. *** Definition de la geometrie 
  37. ******************************************************
  38.  
  39. P1 = 0.  0.  0.;
  40. P2 = h   0.  0.;
  41. P3 = b   h   0.;
  42. P4 = 0.  b   0.;
  43. P5 = 0.  0.  e;
  44.  
  45. L1  = D  nbe  P1  P2;
  46. L2  = D  nbe  P2  P3;
  47. L3  = D  nbe  P3  P4;
  48. L4  = D  nbe  P4  P1;
  49.  
  50. CT1 = L1 ET L2  ET L3 ET L4 ;
  51.  
  52. S1 = SURF CT1 'PLAN';
  53.  
  54. VT = VOLU S1 nbh TRAN P5;
  55.  
  56. S2 = FACE VT 2;
  57.  
  58. ******************************************************
  59. *** Définition du modele
  60. ******************************************************
  61.  
  62. MOD1 = MODE VT MECANIQUE ELASTIQUE ISOTROPE PLASTIQUE OUGLOVA;
  63.  
  64. ******************************************************
  65. *** Parametres materiaux
  66. ****************************************************** 
  67.  
  68. * Module d'Young
  69. ygn = 200.e9;
  70.  
  71. * Coefficient de Poisson
  72. nun = 0.33;
  73.  
  74. * Limite elastique
  75. sigym = 500.e6;
  76.  
  77. * Coefficient d'ecrouissage
  78. km = 500.e6;
  79.  
  80. * Exposant d'ecrouissage
  81. mm = 2.786;
  82.  
  83. * Taux de corrosion en percentage
  84. tcm = 15;
  85.  
  86. * Endommagement critique
  87. dcm = 0.2;
  88.  
  89.  
  90. MAT1 = MATE MOD1 'YOUN' ygn 'NU' nun 
  91.                  'SIGY' sigym 'K   ' km 'M   ' mm
  92.                  'TC  ' tcm 'DC  ' dcm;
  93.  
  94. ******************************************************
  95. *** Definition des conditions aux limites
  96. *** et des déplacements imposés
  97. ******************************************************
  98.  
  99. CL1 = BLOQ S1 UZ;
  100. CL2 = BLOQ P1 UX UY;
  101. CL2 = CL2 ET (BLOQ P2 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. TAB1.'HYPOTHESE_DEFORMATIONS' = 'LINEAIRE' ;
  135.  
  136. PASAPAS TAB1;
  137.  
  138. ******************************************************
  139. *** Verification
  140. ******************************************************
  141.  
  142. EV2 = @GLOBAL TAB1 cl1 EV FZ;
  143. liste EV2;
  144.  
  145. vref = EXTR EV2 'ORDO';
  146. vref = EXTR 16 vref;
  147.  
  148. err =  abs(vref  + 5.87513e8) / 5.87513e8;
  149. message ' erreur relative ';
  150. list err;
  151.  
  152. si (err > 1.e-3);
  153.    erreur (5);
  154. finsi; 
  155.  
  156. FIN;
  157.  
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  

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