Esope source files

   1 : C CREDEP    SOURCE    CHAT      05/01/12    22:29:10     5004           
   2 :       SUBROUTINE CREDEP(AH,DFSIG,PAEC,NSTRS,DEP,BETJEF)
   3 : C
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8 (A-H,O-Z)
   6 :       DIMENSION DEP(4,4),DP(4,4)
   7 :       DIMENSION AC(6),AB(6),DFSIG(4)
   8 :       DIMENSION AH(4,4)
   9 : *
  10 :       SEGMENT BETJEF
  11 :        REAL*8 AA,BETA,RB,ALPHA,YOUN,XNU,GFC,GFT,CAR,ETA,TDEF,
  12 :      &          TCON,DPSTF1,DPSTF2,TETA,PDT
  13 :        INTEGER ICT,ICC,IMOD,IVISS,ITER,
  14 :      &                ISIM,IBB1,IGAU1,IZON
  15 :       ENDSEGMENT
  16 : *
  17 : *     COMMON /DBETJEF/AA,BETA,RB,ALPHA,YOUN,XNU,GFC,GFT,CAR,ETA,TDEF,
  18 : *    &          TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVISS,ITER,
  19 : *    &                ISIM,IBB1,IGAU1,IZON
  20 : C
  21 : C             *************************************************
  22 : C             **  CALCUL DE LA MATRICE ELASTOPLASTIQUE [Dep] **
  23 : C             *************************************************
  24 : C-------------------------------------------------------------------
  24 :       CALL ZERO(DEP,4,4)
  24 :       CALL ZERO(DP,4,4)
  24 :       CALL ZERO(AB,6,1)
  24 :       CALL ZERO(AC,6,1)
  25 : C-------------------------------------------------------------------
  26 : C
  27 :       DO 5 I=1,NSTRS
  28 :       DO 5 J=1,NSTRS
  29 :         AB(I)=AB(I)+AH(I,J)*DFSIG(J)
  30 :     5 CONTINUE
  31 : C
  32 :       DO 6 I=1,NSTRS
  33 :       DO 6 J=1,NSTRS
  34 :        AC(I)=AC(I)+AH(J,I)*DFSIG(J)
  35 :     6 CONTINUE
  36 : C
  37 :       DO 7 I=1,NSTRS
  38 :       DO 7 J=1,NSTRS
  39 :        DP(I,J)=DP(I,J)+AB(I)*AC(J)
  40 :     7 CONTINUE
  41 : C
  42 :       H1=0.D0
  43 :       DO 8 I=1,NSTRS
  44 :         H1=H1+DFSIG(I)*AB(I)
  45 :     8 CONTINUE
  46 : C
  47 :       HA=H1+PAEC
  48 :       DO 9 I=1,NSTRS
  49 :       DO 9 J=1,NSTRS
  50 :        DP(I,J)=DP(I,J)/HA
  51 :        DEP(I,J)=AH(I,J)-DP(I,J)
  52 :     9 CONTINUE
  53 : C
  54 :       DO 10 I=1,NSTRS
  55 :       DO 10 J=1,NSTRS
  56 :         IF (DABS(DEP(I,J)).LT.1.D-5) THEN
  57 :           DEP(I,J)=0.D0
  58 :         ENDIF
  59 :    10 CONTINUE   
  60 : C-------------------------------------------------------------------
  61 :       RETURN
  62 :       END
  63 : 
  64 :  
  65 :