1 : C EQUATN SOURCE CHAT 05/01/12 23:44:36 5004 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)*DEXP(XMAT(7)*SIG) 11 : F2 = XMAT(8)*SIG**XMAT(9) 12 : FS = F1 + F2 13 : R1 = XMAT(10)*DEXP(XMAT(11)*SIG) 14 : R2 = XMAT(12)*SIG**XMAT(13) 15 : RS = R1 + R2 16 : T1 = XMAT(15)*SIG 17 : SINH = 0.5D0*(DEXP(T1)-DEXP(-T1)) 18 : G1 = XMAT(14)*SINH**XMAT(16) 19 : G2 = XMAT(17)*DEXP(XMAT(18)*SIG) 20 : GS = G1 + G2 21 : DEPSCP= FS * RS * DEXP(-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)*DEXP(XMAT(7)*SIG) 30 : F2 = XMAT(8)*SIG**XMAT(9) 31 : FS = F1 + F2 32 : R1 = XMAT(10)*DEXP(XMAT(11)*SIG) 33 : R2 = XMAT(12)*SIG**XMAT(13) 34 : RS = R1 + R2 35 : T1 = XMAT(15)*SIG 36 : SINH = 0.5D0*(DEXP(T1)-DEXP(-T1)) 37 : G1 = XMAT(14)*SINH**XMAT(16) 38 : G2 = XMAT(17)*SIG**XMAT(18) 39 : GS = G1 + G2 40 : DEPSCP= FS * RS * DEXP(-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=DABS(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*(DEXP(CON1)-DEXP(-CON1)))**XMAT(8) 63 : CON2 = XMAT(10)/XMAT(11)*SIG 64 : R = XMAT(9)*(0.5D0*(DEXP(CON2)-DEXP(-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(DABS(T1).LT.100.0) F1=DEXP(-T1) 74 : IF(DABS(T2).LT.100.0) F2=DEXP(-T2) 75 : DEPSCP= EPSTR*F1+EPSX*S*F2 76 : DEPSCS= EPSM 77 : C 78 : ENDIF 79 : RETURN 80 : END 81 : 82 :
© Cast3M 2003 - All rights reserved.
Disclaimer