Télécharger gurson.eso

Retour à la liste

Numérotation des lignes :

gurson
  1. C GURSON SOURCE CHAT 05/01/13 00:22:25 5004
  2. c*****************************************************************************
  3. c sous rouine utilisee dans l'ecoulement plastique de la loi
  4. c de gurson
  5. c**********************************************************************
  6. subroutine gurson (s,p,depsp,gamma,poro,densit,epsi,dtime,
  7. & wrkgur)
  8. c
  9. c this subroutine is the interface between the finite element program
  10. c and the rest of the material update program
  11. c given the history variables and the starin rate it updates the
  12. c state
  13. c by arnaud de gayffier 06/13/93 in san diego california
  14. c
  15. c--- arguments
  16. c s........deviatoric stress in gpa
  17. c s(1)=sxx|
  18. c s(2)=syy| s(1)+s(2)+s(3)=0
  19. c s(3)=szz|
  20. c s(4)=sxy
  21. c s(5)=sxz
  22. c s(6)=syz
  23. c p........pressure (if >0 compression) in gpa
  24. c gamma....cumulative plastic strain
  25. c poro.....porosity or void fraction
  26. c densit...material density in gr/cm3
  27. c epsi.....total strain rate same conventions as "s"
  28. c depsp....total plasti strain rate same convention
  29. c dtime....time increment
  30. c
  31. c warning!! the argument epsi might be modified
  32. c
  33. c--- variables
  34. c
  35. IMPLICIT INTEGER(I-N)
  36. IMPLICIT REAL*8 (A-H,O-Z)
  37. segment wrkgur
  38. real*8 sigbar, sy0,phi0,rho0,g,b,h
  39. real*8 epn,phin,sqrtj2,rho,sig(6)
  40. real*8 e(7),dt
  41. real*8 conv,tol1,tol2
  42. endsegment
  43.  
  44. * COMMON/PROP/SIGBAR,SY0,PHI0,RHO0,G,B,H
  45. * COMMON/STATE/EPN,PHIN,SQRTJ2,RHO,SIG(6)
  46. c e is the strain rate:
  47. c e(1-6) is the deviatoric part
  48. c e(7) is the trace/3
  49. c dt ist the time increment
  50. * COMMON /DELTA/ E(7),DT
  51. dimension s(6),epsi(6),depsp(6)
  52. c
  53. c--- identification
  54. DPHI=0.D0
  55. SQRTJ2=0.D0
  56. epn=gamma
  57. phin=poro
  58. rho=densit
  59. do 100 i=1,6
  60. e(i)=epsi(i)
  61. SIG(I)=S(I)
  62. 100 CONTINUE
  63. dt=dtime
  64. c
  65. c--- partition in deviatoric and hydrostatic
  66. E(7)=(E(1)+E(2)+E(3))/3.D0
  67. E(1)=E(1)-E(7)
  68. E(2)=E(2)-E(7)
  69. E(3)=E(3)-E(7)
  70. E(4)=E(4)/2.D0
  71. E(5)=E(5)/2.D0
  72. E(6)=E(6)/2.D0
  73. c
  74. 101 FORMAT (A,7(1PE12.5,1X))
  75. 102 FORMAT (5(A,1PE12.5))
  76. c
  77. c--- material update
  78. call newsta(wrkgur)
  79. c
  80. c--- arguments update
  81. DPHI = PHIN-PORO
  82. DEPN = EPN-GAMMA
  83. c
  84. P=(1.D0-PHIN)*B*(1.D0-(1.D0-PHIN)*RHO0/((1.D0-PHI0)*RHO))
  85. gamma=epn
  86. poro=phin
  87. densit=rho
  88. do 200 i=1,6
  89. 200 S(I)=SIG(I)
  90. c plastic strain
  91. DO 300 I=1,6
  92. IF (SQRTJ2 .GE. 1.D-6 ) THEN
  93. DEPSP(I)=1.5D0*SIG(I)/SQRTJ2*DEPN
  94. ELSE
  95. DEPSP(I) =0.D0
  96. ENDIF
  97. 300 CONTINUE
  98. c
  99. DEPSP(4)=2.D0*DEPSP(4)
  100. DEPSP(5)=2.D0*DEPSP(5)
  101. DEPSP(6)=2.D0*DEPSP(6)
  102. c
  103. DO 400 I=1,3
  104. DEPSP(I)=DEPSP(I)+DPHI/(1.D0-PHIN)/3.D0
  105. 400 CONTINUE
  106. c
  107. RETURN
  108. end
  109. c------------------------------------------------------------------------------
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  

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