Télécharger ouglova_fibre.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : ouglova_fibre.dgibi
  2. opti epsi lineaire ;
  3. *******************************************************
  4. *******************************************************
  5. *******************************************************
  6. * Cas test de l'implantation numerique du modele
  7. * OUGLOVA 1D MULTIFIBRES
  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. OPTION DIME 2 ELEM QUA4 TRAC PSC;
  23.  
  24. * Nombre d'elements
  25. nelem = 1;
  26.  
  27. ******************************************************
  28. *** Definition de la geometrie
  29. ******************************************************
  30. SACIER1 = 0. 0.;
  31. SACIER2 = 1. 0.;
  32. SACIER3 = 1. 1.;
  33. SACIER4 = 0. 1.;
  34.  
  35. LL1 = D 1 SACIER1 SACIER2;
  36. LL2 = D 1 SACIER2 SACIER3;
  37. LL3 = D 1 SACIER3 SACIER4;
  38. LL4 = D 1 SACIER4 SACIER1;
  39.  
  40. SACIER = DALL LL1 LL2 LL3 LL4;
  41.  
  42. ******************************************************
  43. *** Définition du modele
  44. ******************************************************
  45.  
  46. OPTI DIME 3 ELEM SEG2 TRAC PSC;
  47.  
  48. MOD1 = MODE SACIER MECANIQUE ELASTIQUE PLASTIQUE OUGLOVA
  49. QUAS;
  50.  
  51. ******************************************************
  52. *** Parametres materiaux
  53. ******************************************************
  54.  
  55. * Module d'Young
  56. ygn = 200.e9;
  57.  
  58. * Coefficient de Poisson
  59. nun = 0.33;
  60.  
  61. * Limite elastique
  62. sigym = 500.e6;
  63.  
  64. * Coefficient d'ecrouissage
  65. km = 500.e6;
  66.  
  67. * Exposant d'ecrouissage
  68. mm = 2.786;
  69.  
  70. * Taux de corrosion en percentage
  71. tcm = 15;
  72.  
  73. * Endommagement critique
  74. dcm = 0.2;
  75.  
  76.  
  77. MAT1 = MATE MOD1 'YOUN' ygn 'NU' nun
  78. 'SIGY' sigym 'K ' km 'M ' mm
  79. 'TC ' tcm 'DC ' dcm;
  80.  
  81. CAR1 = CARA MOD1 ALPY .66 ALPZ .66;
  82.  
  83. MODTOT = MOD1;
  84. MATTOT = MAT1 et CAR1;
  85.  
  86. ******************************************************
  87. *** Element poutre
  88. ******************************************************
  89.  
  90. OPTI ELEM SEG2;
  91.  
  92. P1 = 0. 0. 0.;
  93. P2 = 0. 0. 1.;
  94.  
  95. Lpoutre = D nelem P1 P2;
  96.  
  97. MODPOU = MODE Lpoutre MECANIQUE ELASTIQUE SECTION
  98. PLASTIQUE SECTION TIMO;
  99. MATPOU = MATE MODPOU MODS MODTOT MATS MATTOT
  100. VECT (0 1 0);
  101.  
  102. ******************************************************
  103. *** Définition des conditions aux limites
  104. *** et des déplacements imposés
  105. ******************************************************
  106.  
  107. cl1 = BLOQ DEPL ROTA P1;
  108. cl2 = BLOQ UX UY RX RY RZ P2;
  109. cl3 = BLOQ UZ P2;
  110.  
  111. d1 = DEPI cl3 1.0;
  112.  
  113. ******************************************************
  114. *** Définition de chargement
  115. ******************************************************
  116.  
  117. LI1 = PROG 0. 1.;
  118. LI2 = PROG 0. 2.5e-1;
  119. LIS1 = PROG 0. PAS 0.01 1.;
  120.  
  121. EV1 = EVOL MANU LI1 LI2;
  122. CHA1 = CHAR DIMP d1 EV1;
  123.  
  124. ******************************************************
  125. *** Resolution
  126. ******************************************************
  127.  
  128. TAB1 = TABLE;
  129. TAB1 . BLOCAGES_MECANIQUES = cl1 et cl2 et cl3;
  130. TAB1 . MODELE = MODPOU;
  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.  

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