Télécharger elas_hook_endom.dgibi

Retour à la liste

Numérotation des lignes :

  1. ************************************************************************
  2. * Petit test des operateurs ELAS et HOOK dans les cas des modeles *
  3. * endommageables avec fourniture d'un MCHAML de variables internes. *
  4. * On verifie que la matrice de Hooke utilisee tiens bien compte de *
  5. * l'endommagement. *
  6. ************************************************************************
  7.  
  8. ** Options generales
  9. OPTI 'DIME' 3 'ELEM' 'CUB8' 'ECHO' 1 ;
  10.  
  11.  
  12. ** Valeur d'endommagement
  13. d1 = 0.4815162342 ;
  14.  
  15.  
  16. ** Un petit maillage
  17. mail1 = ((DROI 1 (0. 0. 0.) (1. 0. 0.)) TRAN 1 (0. 1. 0.))
  18. VOLU 'TRAN' 1 (0. 0. 1.) ;
  19.  
  20.  
  21. ** Modeles et parametres materiaux (on boucle sur les modeles avec
  22. * endommagement geres par ELAS et HOOK, sauf le FLUAGE CERAMIQUE)
  23. ct = EVOL 'MANU' (PROG 0. 1. 2.) (PROG 0. 1. 2.) ;
  24. * conversion corube traction > ecrouissage
  25. lep1 = ct extr absc ;
  26. lsm1 = ct extr ordo ;
  27. k1 = (extr lsm1 2) / (extr lep1 2) ;
  28. lsm1 = lsm1 enle 1 ;
  29. lep1 = (lep1 enle 1) - (lsm1 / k1) ;
  30. ec = evol manu eps lep1 sig lsm1 ;
  31. REPE b1 6 ;
  32. MESS ;
  33. MESS '---------------------------------------------------------' ;
  34. icas = &b1 ;
  35. SI (EGA icas 1) ;
  36. MESS 'Modele PLASTIQUE ENDOMMAGEABLE' ;
  37. mo = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE' 'ENDOMMAGEABLE';
  38. ma = MATE mo 'YOUN' 1. 'NU' 0.3 'TRAC' ct 'EPSD' 1. 'DC' 1.
  39. 'EPSR' 1. ;
  40. mot1 = MOT 'VHWD' ;
  41. FINSI ;
  42. SI (EGA icas 2) ;
  43. MESS 'Modele VISCOPLASTIQUE VISCODOMMAGE' ;
  44. mo = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'VISCOPLASTIQUE'
  45. 'VISCODOMMAGE' ;
  46. ma = MATE mo 'YOUN' 1. 'NU' 0.3 'N' 1. 'M' 1. 'KK' 1. 'ALP1' 1.
  47. 'BLP1' 1. 'R' 1. 'A' 1. 'EVOL' 1. ;
  48. mot1 = MOT 'VHWD' ;
  49. FINSI ;
  50. SI (EGA icas 3) ;
  51. MESS 'Modele ENDOMMAGEMENT MAZARS' ;
  52. mo = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ENDOMMAGEMENT' 'MAZARS' ;
  53. ma = MATE mo 'YOUN' 1. 'NU' 0.3 'BTRA' 1. 'KTR0' 1. 'ATRA' 1.
  54. 'ACOM' 1. 'BCOM' 1. 'BETA' 1. ;
  55. mot1 = MOT 'D' ;
  56. FINSI ;
  57. SI (EGA icas 4) ;
  58. MESS 'Modele ENDOMMAGEMENT MVM' ;
  59. mo = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'ENDOMMAGEMENT' 'MVM' ;
  60. ma = MATE mo 'YOUN' 1. 'NU' 0.3 'B1' 1. 'B2' 1. 'RATI' 1. 'LOI' 1. ;
  61. mot1 = MOT 'D' ;
  62. FINSI ;
  63. SI (EGA icas 5) ;
  64. MESS 'Modele PLASTIQUE_ENDOM ROUSSELIER' ;
  65. mo = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE_ENDOM'
  66. 'ROUSSELIER' ;
  67. ma = MATE mo 'YOUN' 1. 'NU' 0.3 'ECRO' ec 'F' 0.1 'D' 1. 'SIG1' 1.
  68. 'FC' 0.9 ;
  69. mot1 = MOT 'VHOT' ;
  70. FINSI ;
  71. SI (EGA icas 6) ;
  72. MESS 'Modele PLASTIQUE_ENDOM GURSON2' ;
  73. mo = MODE mail1 'MECANIQUE' 'ELASTIQUE' 'PLASTIQUE_ENDOM'
  74. 'GURSON2' ;
  75. ma = MATE mo 'YOUN' 1. 'NU' 0.3 'ECRO' ec 'Q' 1. 'FU' 1. 'FF' 1.
  76. 'FC' 1. 'FNS0' 1. 'FNE0' 1. 'SNS' 1. 'SNE' 1. 'SIGN' 1.
  77. 'EPSN' 1. 'F0' 0. ;
  78. mot1 = MOT 'VHWF' ;
  79. FINSI ;
  80.  
  81.  
  82. ** Champ de variables internes avec la composante D non nulle
  83. vi0 = ZERO mo 'VARINTER' ;
  84. vi1 = vi0 + (MANU 'CHML' mo mot1 d1 'STRESSES') ;
  85. SI (EGA icas 5) ;
  86. vi1 = vi0 + (MANU 'CHML' mo mot1 (1. - d1) 'STRESSES') ;
  87. FINSI ;
  88.  
  89.  
  90. ** Test du passage deformation --> contrainte avec ELAS et HOOK
  91. ep0 = ZERO mo 'DEFORMAT' ;
  92. ep1 = ep0 + (MANU 'CHML' mo 'EPXX' 1. 'STRESSES') ;
  93. val0 = 1. * (1.-d1) * 1. * (1.-0.3)/((1.+0.3)*(1.-(2.*0.3))) ;
  94. sg1 = ELAS mo ep1 ma vi1 ;
  95. val1 = EXTR sg1 'SMXX' 1 1 1 ;
  96. ec1 = ABS ((val1 - val0) / val0) ;
  97. ho1 = HOOK mo ma vi1 ;
  98. sg2 = ho1 * ep1 mo ;
  99. val2 = EXTR sg2 'SMXX' 1 1 1 ;
  100. ec2 = ABS ((val2 - val0) / val0) ;
  101.  
  102.  
  103. ** Test du passage contrainte --> deformation avec ELAS seulement
  104. sg0 = ZERO mo 'CONTRAIN' ;
  105. sg1 = sg0 + (MANU 'CHML' mo 'SMXX' 1. 'STRESSES') ;
  106. val00 = 1. / (1. * (1.-d1)) ;
  107. ep1 = ELAS mo sg1 ma vi1 ;
  108. val3 = EXTR ep1 'EPXX' 1 1 1 ;
  109. ec3 = ABS ((val3 - val00) / val00) ;
  110.  
  111.  
  112. ** Test des valeurs calculees
  113. MESS 'Contrainte [theorique] : ' val0 ;
  114. MESS 'Contrainte [ELAS] et ecart (%) : ' val1 ec1 ;
  115. MESS 'Contrainte [HOOK] et ecart (%) : ' val2 ec2 ;
  116. MESS 'Deformation [theorique] : ' val00 ;
  117. MESS 'Deformation [ELAS] et ecart (%) : ' val3 ec3 ;
  118. MESS ;
  119.  
  120. ** Erreur si ecart trop grand
  121. SI (ec1 > 1.E-15) ;
  122. MESS '***** Erreur dans l operateur ELAS lors du calcul de '
  123. 'la matrice de Hooke';
  124. ERRE 5 ;
  125. FINSI ;
  126. SI (ec2 > 1.E-15) ;
  127. MESS '***** Erreur dans l operateur HOOK lors du calcul de '
  128. 'la matrice de Hooke';
  129. ERRE 5 ;
  130. FINSI ;
  131. SI (ec3 > 1.E-15) ;
  132. MESS '***** Erreur dans l operateur ELAS lors du calcul de '
  133. 'la matrice de Hooke';
  134. ERRE 5 ;
  135. FINSI ;
  136. FIN b1 ;
  137.  
  138. FIN ;
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  

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