Esope source files

   1 : C VISPLA    SOURCE    CHAT      05/01/13    04:07:29     5004           
   2 :       SUBROUTINE VISPLA(SIGR,SIGF,DSIG1,NSTRS,DPSTF1,DPSTF2,
   3 :      &                      SIGP,SIGRV,DSTRN,DP,BETJEF,VISCO)
   4 : C
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 :       DIMENSION SIGR(4),SIGF(4),DSIG(4),DSTRN(4),D(4,4),EPST(4)
   8 :       DIMENSION DSIG1(4),DSIG2(4),DSIG3(4),DSTRN1(4),DP(4,4)
   9 :       DIMENSION DSIGP(4),SIGP(4),SIGP0(4),SIGRV(4)
  10 : C
  11 : *
  12 :       SEGMENT BETJEF
  13 :        REAL*8 AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,CAR,ETA,TDEF,
  14 :      &          TCON,DPS1,DPS2,TETA,PDT,TP0
  15 :        INTEGER ICT,ICC,IMOD,IVIS,ITER,
  16 :      &                ISIM,IBB1,IGAU1,IZON
  17 :       ENDSEGMENT
  18 :       SEGMENT VISCO
  19 :        REAL*8 DPSTV1,DPSTV2,SIGV1,SIGV2,ENDV
  20 :       ENDSEGMENT
  21 : *
  22 : *     COMMON /DBETJEF/AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,CAR,ETA,TDEF,
  23 : *    &           TCON,DPS1,DPS2,TETA,PDT,ICT,ICC,IMOD,IVIS,ITER,
  24 : *    &                ISIM,IBB1,IGAU1,IZON
  25 : *     COMMON /VISCO/ DPSTV1,DPSTV2,SIGV1,SIGV2
  26 : C
  27 : C
  28 : C ********************************************************************
  29 : C **                 COMPORTEMENT VISCOPLASTIQUE                    **
  30 : C **                    DE TYPE DUVAUT LIONS                        **
  31 : C ********************************************************************
  32 : C
  33 :       IPAS=0
  33 :       CALL ZERO(DSIG2,4,1)
  33 :       CALL ZERO(DSIG3,4,1)
  33 :       CALL ZERO(DSTRN1,4,1)
  34 : C
  35 : C -------------- recalcul des increments de contrainte  -------------
  36 : C
  37 : C
  38 :       DO 2 I=1,NSTRS
  39 :         SIGP(I)=SIGF(I)
  40 :     2 CONTINUE
  41 : C
  42 :       DO 3 I=1,NSTRS
  43 :         DSIGP(I)=SIGF(I)-SIGR(I)
  44 :     3 CONTINUE
  45 : C
  46 : C ---------------- calcul des termes viscoplastiques -----------------
  47 : C
  48 :       GAMA = 0.5
  49 :       A = (ETA)/(ETA+GAMA*PDT)
  50 :       B = (GAMA*PDT)/ETA
  51 :       C = PDT
  52 :       G=EXP(-PDT/ETA)
  53 :       IF (C.EQ.0.D0) THEN
  54 :       WRITE(*,*)'l increment de temps est de longueur nulle'
  55 :       stop
  56 :       ENDIF
  57 : C
  58 : C -------------------- mise a jour des contraintes -------------------
  59 : C 
  60 :       IF (IMOD.EQ.1) THEN
  60 :         CALL CREMAT(D,YOUN,XNU,3,-2) 
  61 :       ENDIF 
  62 :       IF (IMOD.EQ.2) THEN
  62 :         CALL CREMAT(D,YOUN,XNU,3,-1) 
  63 :       ENDIF
  64 : C
  65 : C ------------ Calcul de la matrice tangente consistante --------------
  66 : C
  67 :       DO 10 I=1,NSTRS        
  68 :       DO 10 J=1,NSTRS
  69 :         DP(I,J)=0.D0
  70 :         DP(I,J)=D(I,J)*A
  71 :    10 CONTINUE 
  72 : C
  73 : C -------------------- algorithme type SLUYS -------------------
  74 : C 
  75 :       IF (ISIM.EQ.0) THEN
  76 : C
  77 :       DO 20 I=1,NSTRS
  78 :         DSIG2(I)=0.D0
  79 :       DO 20 J=1,NSTRS
  80 :         DSIG2(I)=DSIG2(I)+D(I,J)*EPST(J)
  81 :    20 CONTINUE
  82 : C
  83 :       DO 30 I=1,NSTRS
  84 :         DSIG3(I)=A*(DSIG1(I)-C*DSIG2(I)+B*DSIGP(I))
  85 :    30 CONTINUE 
  86 : C
  87 :       DO 40 I=1,NSTRS
  88 :         SIGF(I)=SIGR(I)+DSIG3(I)
  89 :    40 CONTINUE
  90 :       ENDIF
  91 : C
  92 : C -------------------- algorithme type SIMO -------------------
  93 : C 
  94 :       IF (ISIM.EQ.1) THEN
  94 :       CALL ZERO(SIGF,4,1)
  95 :       DO 35 I=1,NSTRS
  96 :         SIGF(I)=G*SIGRV(I)+(1.D0-G)*SIGP(I)+((1.D0-G)/(C/ETA))*DSIG1(I)
  97 :    35 CONTINUE 
  98 :       ENDIF
  99 : C
 100 : C ------- mise a jour des variables d'ecrouissages --------
 101 : C 
 102 :       DPSTV1= DPSTV1*G+(1.D0-G)*DPSTF1
 102 :       CALL UNICOU(DPSTV1,PAEC,1,SIGV1,BETJEF)
 103 :       DPSTV2= DPSTV2*G+(1.D0-G)*DPSTF2 
 103 :       CALL UNICOU(DPSTV2,PAEC,2,SIGV2,BETJEF)
 104 : C
 105 : C ------- mise a jour des taux de deformation --------
 106 : C 
 107 :       IF (IMOD.EQ.1) THEN
 107 :         CALL MATSOU(D,YOUN,XNU,3,-2) 
 108 :       ENDIF 
 109 :       IF (IMOD.EQ.2) THEN
 109 :         CALL MATSOU(D,YOUN,XNU,4,-1) 
 110 :       ENDIF
 111 : C
 112 : C -------  taux de deformation viscoplastique --------
 113 : C
 114 :       DO 50 I=1,NSTRS
 115 :         DSTRN1(I)=0.D0
 116 :       DO 50 J=1,NSTRS
 117 :         DSTRN1(I)=DSTRN1(I)+D(I,J)*(DSIG3(J)-DSIGP(J))/ETA
 118 :    50 CONTINUE
 119 : C
 120 :       DO 60 I=1,NSTRS
 121 :         EPST(I)=EPST(I)+DSTRN1(I)
 122 :    60 CONTINUE
 123 : C 
 124 : C
 125 : C -------  taux de deformation en norme --------
 126 : C
 127 :       TDEF=0.D0
 128 :       D1=DSTRN(1)
 129 :       D2=DSTRN(2)
 130 :       D3=DSTRN(3)
 131 :       D4=DSTRN(4)
 132 :       TDEF=SQRT(D1*D1+D2*D2+D3*D3+D4*D4)/PDT
 133 : C
 134 : C -------  taux de contrainte en norme --------
 135 : C
 136 :       TCON=0.D0
 137 :       DO 62 I=1,NSTRS
 138 :       DSIG(I) = SIGF(I) - SIGRV(I)
 139 :    62 CONTINUE
 140 :       S1=DSIG(1)
 141 :       S2=DSIG(2)
 142 :       S3=DSIG(3)
 143 :       S4=DSIG(4)
 144 :       TCON=SQRT(S1*S1+S2*S2+S3*S3+S4*S4)/PDT
 145 : C
 145 :       CALL INDICA(DPSTV1,DPSTV2,IFIS,IPLA,3,BETJEF)
 146 :       RETURN
 147 :       END
 148 : 
 149 :  
 150 :  
 151 :