Numérotation des lignes :

equatn
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