Télécharger incre8.eso

Retour à la liste

Numérotation des lignes :

  1. C INCRE8 SOURCE CHAT 05/01/13 00:35:34 5004
  2.  
  3. SUBROUTINE INCRE8(SIG,VAR,T0,TF,EPSVPT,VARPT,XMAT,NSTRS0,
  4. & MFR,NVARI,NCOMAT)
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. DIMENSION SIG(*),VAR(*),EPSVPT(*)
  8. DIMENSION VARPT(*),XMAT(*)
  9. DIMENSION XX(6)
  10. DIMENSION SIG0(6),EPS0(6)
  11. PARAMETER (AMAX = 1.0D20 , AMIN = 1.D-10 ,ZER=0.D0 )
  12. DETIER = 2.0D0/3.0D0
  13. ROOT = SQRT(3.0D0/2.0D0)
  14. CALL ZDANUL(VARPT,NVARI)
  15. IF (MFR.EQ.5) THEN
  16. NSTRS=6
  17. SIG0(1)=SIG(1)
  18. SIG0(2)=SIG(2)
  19. SIG0(3)=0.D0
  20. SIG0(4)=SIG(3)
  21. SIG0(5)=SIG(4)
  22. SIG0(6)=SIG(5)
  23. ELSE
  24. NSTRS=NSTRS0
  25. DO 10 I=1,NSTRS
  26. SIG0(I)=SIG(I)
  27. 10 CONTINUE
  28. ENDIF
  29. DO 20 I=1,NSTRS
  30. EPS0(I) = ZER
  31. 20 CONTINUE
  32. C-------------------------------------------------------------------|
  33. C COEFFICIENTS MATERIAU |
  34. C-------------------------------------------------------------------|
  35. YOU =XMAT( 1)
  36. XNU =XMAT( 2)
  37. AA=XMAT( 5)
  38. B =XMAT( 6)
  39. RM =XMAT( 7)
  40. Q =XMAT( 8)
  41. R =XMAT( 9)
  42. H0 =XMAT(10)
  43. AP =XMAT(11)
  44. SB =XMAT(12)
  45. RN =XMAT(13)
  46. S0 =XMAT(14)
  47. C
  48. C ------------------------------------------------------------------|
  49. C PARTIE VISCOPLASTIQUE |
  50. C-------------------------------------------------------------------I
  51. C
  52. C
  53. TRACE =(SIG0(1)+SIG0(2)+SIG0(3))/3.0D0
  54. DO 71 I=1,NSTRS
  55. A = 0.0D0
  56. IF (I.LE.3) A=1.0D0
  57. XX(I) = SIG0(I)-A*TRACE
  58. 71 CONTINUE
  59. C
  60. C *** CRITERE VISCOPLASTIQUE
  61. C
  62. AJ2 = PROCON(XX,XX,NSTRS)
  63. AJ2 = SQRT(1.5D0*AJ2)
  64. IF (VAR(2).LT.AMIN) VAR(2)= S0
  65. RR = SINH (B * AJ2/VAR(2))
  66. C
  67. C temperature au debut du pas
  68. C
  69. QSRT = Q/(R*T0)
  70. Z = AA*XPUISS(RR,1.D0/RM,AMAX)
  71. VARPT(1) = Z*EXP(-QSRT)
  72. SSAT = SB * XPUISS(Z/AA,RN,AMAX)
  73. DSAT = SSAT - S0
  74. SSIG = ABS(DSAT )/DSAT
  75. RK = ABS((SSAT - VAR(2))/ DSAT)
  76. VARPT(2)= H0*XPUISS(RK,AP,AMAX)*SSIG*VARPT(1)
  77. IF (AJ2.LT.AMIN) GOTO 100
  78. DO 70 I=1,NSTRS
  79. EPS0(I) =XX(I)/AJ2*VARPT(1)/DETIER
  80. 70 CONTINUE
  81. C
  82. 100 CONTINUE
  83. C
  84. IF (MFR.EQ.5) THEN
  85. EPSVPT(1)=EPS0(1)
  86. EPSVPT(2)=EPS0(2)
  87. EPSVPT(3)=EPS0(4)
  88. EPSVPT(4)=EPS0(5)
  89. EPSVPT(5)=EPS0(6)
  90. ELSE
  91. DO 11 I=1,NSTRS
  92. EPSVPT(I)=EPS0(I)
  93. 11 CONTINUE
  94. ENDIF
  95. RETURN
  96. END
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  

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