Télécharger equatn.eso

Retour à la liste

Numérotation des lignes :

  1. C EQUATN SOURCE CB215821 16/04/21 21:16:39 8920
  2. SUBROUTINE EQUATN(INPLAS,SIG,TIME,DT,XMAT,DEPSCP,DEPSCS)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL *8 (A-H,O-Z)
  5. DIMENSION XMAT(*)
  6. IF(INPLAS.EQ.20) THEN
  7. C
  8. C BLACKBURN EQUATION
  9. C
  10. F1 = XMAT(6)*EXP(XMAT(7)*SIG)
  11. F2 = XMAT(8)*SIG**XMAT(9)
  12. FS = F1 + F2
  13. R1 = XMAT(10)*EXP(XMAT(11)*SIG)
  14. R2 = XMAT(12)*SIG**XMAT(13)
  15. RS = R1 + R2
  16. T1 = XMAT(15)*SIG
  17. SINH = 0.5D0*(EXP(T1)-EXP(-T1))
  18. G1 = XMAT(14)*SINH**XMAT(16)
  19. G2 = XMAT(17)*EXP(XMAT(18)*SIG)
  20. GS = G1 + G2
  21. DEPSCP= FS * RS * EXP(-RS*TIME)
  22. DEPSCS= GS
  23. RETURN
  24. C
  25. ELSE IF(INPLAS.EQ.61) THEN
  26. C
  27. C BLACKBURN_2 EQUATION
  28. C
  29. F1 = XMAT(6)*EXP(XMAT(7)*SIG)
  30. F2 = XMAT(8)*SIG**XMAT(9)
  31. FS = F1 + F2
  32. R1 = XMAT(10)*EXP(XMAT(11)*SIG)
  33. R2 = XMAT(12)*SIG**XMAT(13)
  34. RS = R1 + R2
  35. T1 = XMAT(15)*SIG
  36. SINH = 0.5D0*(EXP(T1)-EXP(-T1))
  37. G1 = XMAT(14)*SINH**XMAT(16)
  38. G2 = XMAT(17)*SIG**XMAT(18)
  39. GS = G1 + G2
  40. DEPSCP= FS * RS * EXP(-RS*TIME)
  41. DEPSCS= GS
  42. RETURN
  43. C
  44. ELSE IF(INPLAS.EQ.22) THEN
  45. C
  46. C RCC-MR EQUATION FOR 1S MATERIAL (SS 316LN)
  47. C
  48. DEPSCP= 0.0D0
  49. DEPSCS= 0.0D0
  50. TFP = XMAT(9)*SIG**XMAT(10)
  51. IF(TIME.LT.DT) TIME=ABS(DT)
  52. IF(TIME.LE.TFP)
  53. & DEPSCP=(XMAT(6)*XMAT(7))*(TIME**(XMAT(7)-1.0D0))*(SIG**XMAT(8))
  54. IF(TIME.GT.TFP)
  55. & DEPSCS=XMAT(11)*SIG**XMAT(12)
  56. RETURN
  57. ELSE IF(INPLAS.EQ.23) THEN
  58. C
  59. C RCC-MR EQUATION FOR 2S MATERIAL (SS 304LN)
  60. C
  61. CON1 = XMAT(7)/XMAT(8)*SIG
  62. EPSM = XMAT(6)*(0.5D0*(EXP(CON1)-EXP(-CON1)))**XMAT(8)
  63. CON2 = XMAT(10)/XMAT(11)*SIG
  64. R = XMAT(9)*(0.5D0*(EXP(CON2)-EXP(-CON2)))**XMAT(11)
  65. EPSTR = XMAT(12)*EPSM
  66. S = XMAT(13)/XMAT(9)*R
  67. EPSX = 0.0D0
  68. IF(SIG.GT.XMAT(16)) EPSX=XMAT(14)+XMAT(15)*SIG
  69. T1 = R*TIME
  70. T2 = S*TIME
  71. F1 = 0.0D0
  72. F2 = 0.0D0
  73. IF(ABS(T1).LT.100.0) F1=EXP(-T1)
  74. IF(ABS(T2).LT.100.0) F2=EXP(-T2)
  75. DEPSCP= EPSTR*F1+EPSX*S*F2
  76. DEPSCS= EPSM
  77. C
  78. ENDIF
  79. RETURN
  80. END
  81.  
  82.  
  83.  
  84.  
  85.  

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