Télécharger la plaquette Cast3M


fr en

Télécharger gurson.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  gurson.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ***************************************************
  5. ****        Dilatation uniforme d'un cub8       ***
  6. ****    Test sur l'implementation de GURSON     ***
  7. ***************************************************
  8. *
  9. * Date 16/02/94
  10. * Auteur A de Gayffier
  11. * Reference rapport DMT/94-108
  12. *
  13. * Le test consiste à imposer une dillatation uniforme
  14. * d'un cube à 8 noeuds et verifier la cohérence du champ 
  15. * de contrainte obtenu avec un modele Plastique de Gurson
  16. *
  17. * En effet l'etat de contrainte est hydrostatique et verifie la
  18. * relation ln(poro) +1.5P/sbar =< 0
  19. *
  20. option dime 3 elem cub8;
  21. opti echo 0 ;
  22. *
  23. *maillage
  24. *
  25. densite  1. ;
  26. p1 = 0. 0. 0.  ;
  27. p2 = 1. 0. 0.  ;
  28. p3 = 1. 1. 0.  ;
  29. p4 = 0. 1. 0.  ;
  30. p5 = 0. 0. 0.5 ;
  31. p6 = 1. 0. 0.5 ;
  32. p7 = 1. 1. 0.5 ;
  33. p8 = 0. 1. 0.5 ;
  34. vol = manu cub8 p1 p2 p3 p4 p5 p6 p7 p8 ;
  35. *
  36. *modele de calcul
  37. *
  38. K = 139.7e9 ;
  39. G = 46.6e9 ;
  40. poisson = (1.-(G/K*2./3.) ) / (2./3.*G/K + 2.) ;
  41. E = 3.*K*(1.-(2.*poisson)) ;
  42. list poisson ;
  43. list E ;
  44. phi0 = 3.e-4 ;
  45. sigbar = 3.e8 ;
  46. *
  47. mo = MODE vol mecanique elastique plastique gurson cub8 ;
  48. ma = MATE mo 'YOUNG' E 'NU' poisson 'SIGY' 2.6E8 'H' 2.5e8 'SBAR'
  49. sigbar 'PORO' phi0 ;
  50. *
  51. *rigidite
  52. *
  53. ri = rigi mo ma;
  54. *
  55. *conditions aux limites
  56. lambda = 10. ;
  57. cl1 = bloqu ux ( p2 et p3 et p6 et p7 ) ;
  58. cl2 = bloqu uy ( p3 et p4 et p7 et p8 )  ;
  59. cl3 = bloqu uz ( p5 et p6 et p7 et p8 )  ;
  60. cl4 = bloqu ux ( p1 et p4 et p5 et p8 ) ;
  61. cl5 = bloqu uy ( p1 et p2 et p5 et p6 ) ;
  62. cl6 = bloqu uz ( p1 et p2 et p3 et p4 ) ;
  63. stf = ri et cl1 et cl2 et cl3 et cl4 et cl5 et cl6 ;
  64. cl =  cl1 et cl2 et cl3 et cl4 et cl5 et cl6 ;
  65. *
  66. fo1 = depimp cl1 lambda ;
  67. fo2 = depimp cl2 lambda ;
  68. fo3 = depimp cl3 (lambda*0.5) ;
  69. for = fo1 et fo2 et fo3 ;
  70. dep = reso stf for ;
  71. *
  72. * table pour non lineaire
  73. *
  74. xtab = table ;
  75. *
  76. *   chronologie du chargement
  77. ltime = prog 0. pas 1.  10. 11. ;
  78. letir = prog  0. pas  0.1 1. 1.;
  79. evau = evol manu time ltime etir letir ;
  80. char1 = char dimp for evau;
  81. *
  82. * temps pour la sortie des resultats 
  83. tres = prog  1. pas 1. 10. ;
  84. *
  85. * Il est nécessaire de definir les variables internes au debut car elle
  86. * ne sont pas nulles (l'operateur ZERO initialise par défaut à 0)
  87. xtab.'VARIABLES_INTERNES' = TABLE;
  88. xtab.'VARIABLES_INTERNES' . 0 = manu chml mo 'EPSE' 0. 'VPOR' 3.e-4
  89.      'VHNS' 1.  'TYPE' 'VARIABLES INTERNES' 'STRESSES';
  90. * 
  91. * iterations plastiques
  92. xtab.blocages_mecaniques = cl;
  93. xtab.caracteristiques = ma;
  94. xtab.modele = mo;
  95. xtab.chargement = char1;
  96. xtab.temps_calcules = tres;
  97. pasapas xtab ;
  98. * 
  99. *  controle des resultats
  100. *
  101. i = 1 ;
  102. repeter boucle 10 ;
  103. var1 =  xtab.variables_internes.i ;
  104. phi = extr var1 vpor 1 1 1 ;
  105. rho = extr var1 vhns  1 1 1 ;
  106. sig1 = xtab.contraintes.i ;
  107. sxx = extr sig1 SMXX 1 1 1 ;
  108. syy = extr sig1 SMYY 1 1 1 ;
  109. szz = extr sig1 SMZZ 1 1 1 ;
  110. p = ( sxx + szz + syy ) /(-3.) ;
  111. pp = K*(1.-phi)*(1.-((1-phi)/(1-phi0)/rho)) ;
  112. conv1 =( 1.5 * p / sigbar) - (log phi) ;
  113. conv2 =( 1.5 * pp / sigbar) - (log phi) ;
  114. *
  115. *  Tests 
  116. err1 = abs(conv1) > 1.e-6 ;
  117. err2 = abs(conv2) > 1.e-6 ;
  118. si err1 ;
  119.   si (ou (phi > phi0) (conv1 < 0.));
  120.     ERREUR 5 ;
  121.   finsi ;
  122. finsi;
  123. si err2 ;
  124.   si (ou (phi > phi0) (conv2 < 0.));
  125.     ERREUR 5 ;
  126.   finsi ;
  127. finsi ;
  128. i = i + 1 ;
  129. fin boucle ;
  130. fin ;
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  

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