Télécharger la plaquette Cast3M


fr en

Télécharger equatn.eso

Retour à la liste

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