Télécharger ouglova_fibre.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : ouglova_fibre.dgibi
  2. *******************************************************
  3. *******************************************************
  4. *******************************************************
  5. * Cas test de l'implantation numerique du modele
  6. * OUGLOVA 1D MULTIFIBRES
  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. OPTION DIME 2 ELEM QUA4 TRAC PSC;
  22.  
  23. * Nombre d'elements
  24. nelem = 1;
  25.  
  26. ******************************************************
  27. *** Definition de la geometrie
  28. ******************************************************
  29. SACIER1 = 0. 0.;
  30. SACIER2 = 1. 0.;
  31. SACIER3 = 1. 1.;
  32. SACIER4 = 0. 1.;
  33.  
  34. LL1 = D 1 SACIER1 SACIER2;
  35. LL2 = D 1 SACIER2 SACIER3;
  36. LL3 = D 1 SACIER3 SACIER4;
  37. LL4 = D 1 SACIER4 SACIER1;
  38.  
  39. SACIER = DALL LL1 LL2 LL3 LL4;
  40.  
  41. ******************************************************
  42. *** Définition du modele
  43. ******************************************************
  44.  
  45. OPTI DIME 3 ELEM SEG2 TRAC PSC;
  46.  
  47. MOD1 = MODE SACIER MECANIQUE ELASTIQUE PLASTIQUE OUGLOVA
  48. QUAS;
  49.  
  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
  77. 'SIGY' sigym 'K ' km 'M ' mm
  78. 'TC ' tcm 'DC ' dcm;
  79.  
  80. CAR1 = CARA MOD1 ALPY .66 ALPZ .66;
  81.  
  82. MODTOT = MOD1;
  83. MATTOT = MAT1 et CAR1;
  84.  
  85. ******************************************************
  86. *** Element poutre
  87. ******************************************************
  88.  
  89. OPTI ELEM SEG2;
  90.  
  91. P1 = 0. 0. 0.;
  92. P2 = 0. 0. 1.;
  93.  
  94. Lpoutre = D nelem P1 P2;
  95.  
  96. MODPOU = MODE Lpoutre MECANIQUE ELASTIQUE SECTION
  97. PLASTIQUE SECTION TIMO;
  98. MATPOU = MATE MODPOU MODS MODTOT MATS MATTOT
  99. VECT (0 1 0);
  100.  
  101. ******************************************************
  102. *** Définition des conditions aux limites
  103. *** et des déplacements imposés
  104. ******************************************************
  105.  
  106. cl1 = BLOQ DEPL ROTA P1;
  107. cl2 = BLOQ UX UY RX RY RZ P2;
  108. cl3 = BLOQ UZ P2;
  109.  
  110. d1 = DEPI cl3 1.0;
  111.  
  112. ******************************************************
  113. *** Définition de chargement
  114. ******************************************************
  115.  
  116. LI1 = PROG 0. 1.;
  117. LI2 = PROG 0. 2.5e-1;
  118. LIS1 = PROG 0. PAS 0.01 1.;
  119.  
  120. EV1 = EVOL MANU LI1 LI2;
  121. CHA1 = CHAR DIMP d1 EV1;
  122.  
  123. ******************************************************
  124. *** Resolution
  125. ******************************************************
  126.  
  127. TAB1 = TABLE;
  128. TAB1 . BLOCAGES_MECANIQUES = cl1 et cl2 et cl3;
  129. TAB1 . MODELE = MODPOU;
  130. TAB1. 'HYPOTHESE_DEFORMATIONS' = MOT 'LINEAIRE';
  131. TAB1 . CARACTERISTIQUES = MATPOU;
  132. TAB1 . CHARGEMENT = CHA1;
  133. TAB1 . TEMPS_CALCULES = LIS1;
  134. TAB1 . MAXDEFOR = 1.E-5;
  135. TAB1 . MOVA = RIEN;
  136.  
  137. PASAPAS TAB1;
  138.  
  139. ******************************************************
  140. *** Verification
  141. ******************************************************
  142.  
  143. EV2 = @GLOBAL TAB1 cl1 EV1 FZ;
  144. liste EV2;
  145.  
  146. vref = EXTR EV2 'ORDO';
  147. vref = EXTR 29 vref;
  148.  
  149. err = abs(vref + 5.68172e8) / 5.68172e8;
  150. message ' erreur relative ';
  151. list err;
  152.  
  153. si (err > 1.e-3);
  154. erreur (5);
  155. finsi;
  156.  
  157. fin;
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  
  177.  
  178.  
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  
  185.  
  186.  

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