Télécharger la plaquette Cast3M


fr en

Télécharger depo.eso

Retour à la liste

Numérotation des lignes : 

depo

  1. C DEPO      SOURCE    CB215821  16/04/21    21:16:23     8920
  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-------------------------------------------------------------------
  15.       CALL ZERO(DEP,NSTRS,NSTRS)
  16.       CALL ZERO(DP,NSTRS,NSTRS)
  17.       CALL ZERO(AB,6,1)
  18.       CALL ZERO(AC,6,1)
  19.       CALL ZERO(DFSI,6,1)
  20.       CALL ZERO(DGSI,6,1)
  21. C-------------------------------------------------------------------
  22.       CALL MATHOO(D,EX,PXY,NSTRS,IFOU)
  23.       CALL DFSIG(S1,DFSI,DGSI,SEQ,NSTRS,RB,AA,BB,DK1,DK2,ILO)
  24. C
  25.       DO 5 I=1,NSTRS
  26.       DO 5 J=1,NSTRS
  27.         AB(I)=AB(I)+D(I,J)*DGSI(J)
  28.     5 CONTINUE
  29. C
  30.       DO 2 I=1,NSTRS
  31.       DO 2 J=1,NSTRS
  32.        AC(I)=AC(I)+D(J,I)*DFSI(J)
  33.     2 CONTINUE
  34. C
  35.       DO 3 I=1,NSTRS
  36.       DO 3 J=1,NSTRS
  37.        DP(I,J)=DP(I,J)+AB(I)*AC(J)
  38.     3 CONTINUE
  39. C
  40.       H1=0.D0
  41.       DO 4 I=1,NSTRS
  42.         H1=H1+DFSI(I)*AB(I)
  43.     4 CONTINUE
  44. C
  45.       HA=H1+PAEC
  46.       DO 1 I=1,NSTRS
  47.       DO 1 J=1,NSTRS
  48.        IF(ABS(HA).LT.1.D-10) THEN
  49.          DP(I,J)=0.D0
  50.        ELSE
  51.          DP(I,J)=DP(I,J)/HA
  52.        ENDIF
  53.        DEP(I,J)=D(I,J)-DP(I,J)
  54.     1 CONTINUE
  55. C-------------------------------------------------------------------
  56.       IECR=0
  57.       IF(IECR.EQ.10) THEN
  58.       WRITE(*,*) '***  S1 DANS DEPO  ***  PAEC=',PAEC,'SEQ=',SEQ
  59.       WRITE(*,1991) (S1(IA),IA=1,NSTRS)
  60.       WRITE(*,*) '***  DFSI DANS DEPO  ***  IFOU=',IFOU,'NSTRS=',NSTRS
  61.       WRITE(*,1991) (DFSI(IA),IA=1,NSTRS)
  62.       WRITE(*,*) '***  DGSI ***  H1=',H1,' HA=',HA
  63.       WRITE(*,1991) (DGSI(IA),IA=1,NSTRS)
  64.       WRITE(*,*) ' ***  D (cont) **   '
  65.       DO 115 IA = 1,NSTRS
  66.         WRITE(*,1991) (D(IA,IC),IC=1,NSTRS)
  67.   115 CONTINUE
  68.       WRITE(*,*) ' ***  DEP **  '
  69.       DO 114 IA = 1,NSTRS
  70.         WRITE(*,1991) (DEP(IA,IC),IC=1,NSTRS)
  71.   114 CONTINUE
  72.       ENDIF
  73. 1991  FORMAT(18(1X,E12.5))
  74. C-------------------------------------------------------------------
  75.       RETURN
  76.       END
  77.  
  78.  
  79.  
  80.  
  81.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales