Télécharger la plaquette Cast3M


fr en

Télécharger incre8.eso

Retour à la liste

Numérotation des lignes : 

incre8

  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