Esope source files

   1 : C DEPO      SOURCE    CHAT      05/01/12    22:42:46     5004           
   2 : C DEPO     SOURCE    INSL       24/10/96    
   3 :       SUBROUTINE DEPO(S1,DEP,PAEC,SEQ,NSTRS,IFOU,D,DP,EX,PXY,AA,BB,
   4 :      1 DK1,DK2,RB,ILO)
   5 : C
   6 :       IMPLICIT INTEGER(I-N)
   7 :       IMPLICIT REAL*8(A-H,O-Z)
   8 :       DIMENSION S1(NSTRS),DEP(NSTRS,NSTRS),DP(NSTRS,NSTRS)
   9 :       DIMENSION D(NSTRS,NSTRS),AC(6),AB(6),DGSI(6),DFSI(6)
  10 : C
  11 : C             *************************************************
  12 : C             **  CALCUL DE LA MATRICE ELASTOPLASTIQUE [Dep] **
  13 : C             *************************************************
  14 : C-------------------------------------------------------------------
  14 :       CALL ZERO(DEP,NSTRS,NSTRS)
  14 :       CALL ZERO(DP,NSTRS,NSTRS)
  14 :       CALL ZERO(AB,6,1)
  14 :       CALL ZERO(AC,6,1)
  14 :       CALL ZERO(DFSI,6,1)
  14 :       CALL ZERO(DGSI,6,1)
  15 : C-------------------------------------------------------------------
  15 :       CALL MATHOO(D,EX,PXY,NSTRS,IFOU)
  15 :       CALL DFSIG(S1,DFSI,DGSI,SEQ,NSTRS,RB,AA,BB,DK1,DK2,ILO)
  16 : C
  17 :       DO 5 I=1,NSTRS
  18 :       DO 5 J=1,NSTRS
  19 :         AB(I)=AB(I)+D(I,J)*DGSI(J)
  20 :     5 CONTINUE
  21 : C
  22 :       DO 2 I=1,NSTRS
  23 :       DO 2 J=1,NSTRS
  24 :        AC(I)=AC(I)+D(J,I)*DFSI(J)
  25 :     2 CONTINUE
  26 : C
  27 :       DO 3 I=1,NSTRS
  28 :       DO 3 J=1,NSTRS
  29 :        DP(I,J)=DP(I,J)+AB(I)*AC(J)
  30 :     3 CONTINUE
  31 : C
  32 :       H1=0.D0
  33 :       DO 4 I=1,NSTRS
  34 :         H1=H1+DFSI(I)*AB(I)
  35 :     4 CONTINUE
  36 : C
  37 :       HA=H1+PAEC
  38 :       DO 1 I=1,NSTRS
  39 :       DO 1 J=1,NSTRS
  40 :        IF(DABS(HA).LT.1.D-10) THEN
  41 :          DP(I,J)=0.D0
  42 :        ELSE
  43 :          DP(I,J)=DP(I,J)/HA
  44 :        ENDIF
  45 :        DEP(I,J)=D(I,J)-DP(I,J)
  46 :     1 CONTINUE
  47 : C-------------------------------------------------------------------
  48 :       IECR=0
  49 :       IF(IECR.EQ.10) THEN
  50 :       WRITE(*,*) '***  S1 DANS DEPO  ***  PAEC=',PAEC,'SEQ=',SEQ
  51 :       WRITE(*,1991) (S1(IA),IA=1,NSTRS)
  52 :       WRITE(*,*) '***  DFSI DANS DEPO  ***  IFOU=',IFOU,'NSTRS=',NSTRS
  53 :       WRITE(*,1991) (DFSI(IA),IA=1,NSTRS)
  54 :       WRITE(*,*) '***  DGSI ***  H1=',H1,' HA=',HA
  55 :       WRITE(*,1991) (DGSI(IA),IA=1,NSTRS)
  56 :       WRITE(*,*) ' ***  D (cont) **   '
  57 :       DO 115 IA = 1,NSTRS
  58 :         WRITE(*,1991) (D(IA,IC),IC=1,NSTRS)
  59 :   115 CONTINUE
  60 :       WRITE(*,*) ' ***  DEP **  '
  61 :       DO 114 IA = 1,NSTRS
  62 :         WRITE(*,1991) (DEP(IA,IC),IC=1,NSTRS)
  63 :   114 CONTINUE
  64 :       ENDIF
  65 : 1991  FORMAT(18(1X,E12.5))
  66 : C-------------------------------------------------------------------
  67 :       RETURN
  68 :       END
  69 :  
  70 :