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

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