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  26/08/24    21:16:04     12622          
  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 1992 I=1,NSTRS
  26.       DO 5 J=1,NSTRS
  27.         AB(I)=AB(I)+D(I,J)*DGSI(J)
  28.     5 CONTINUE
  29.  1992 CONTINUE
  30. C
  31.       DO 1993 I=1,NSTRS
  32.       DO 2 J=1,NSTRS
  33.        AC(I)=AC(I)+D(J,I)*DFSI(J)
  34.     2 CONTINUE
  35.  1993 CONTINUE
  36. C
  37.       DO 1994 I=1,NSTRS
  38.       DO 3 J=1,NSTRS
  39.        DP(I,J)=DP(I,J)+AB(I)*AC(J)
  40.     3 CONTINUE
  41.  1994 CONTINUE
  42. C
  43.       H1=0.D0
  44.       DO 4 I=1,NSTRS
  45.         H1=H1+DFSI(I)*AB(I)
  46.     4 CONTINUE
  47. C
  48.       HA=H1+PAEC
  49.       DO 1995 I=1,NSTRS
  50.       DO 1 J=1,NSTRS
  51.        IF(ABS(HA).LT.1.D-10) THEN
  52.          DP(I,J)=0.D0
  53.        ELSE
  54.          DP(I,J)=DP(I,J)/HA
  55.        ENDIF
  56.        DEP(I,J)=D(I,J)-DP(I,J)
  57.     1 CONTINUE
  58.  1995 CONTINUE
  59. C-------------------------------------------------------------------
  60.       IECR=0
  61.       IF(IECR.EQ.10) THEN
  62.       WRITE(*,*) '***  S1 DANS DEPO  ***  PAEC=',PAEC,'SEQ=',SEQ
  63.       WRITE(*,1991) (S1(IA),IA=1,NSTRS)
  64.       WRITE(*,*) '***  DFSI DANS DEPO  ***  IFOU=',IFOU,'NSTRS=',NSTRS
  65.       WRITE(*,1991) (DFSI(IA),IA=1,NSTRS)
  66.       WRITE(*,*) '***  DGSI ***  H1=',H1,' HA=',HA
  67.       WRITE(*,1991) (DGSI(IA),IA=1,NSTRS)
  68.       WRITE(*,*) ' ***  D (cont) **   '
  69.       DO 115 IA = 1,NSTRS
  70.         WRITE(*,1991) (D(IA,IC),IC=1,NSTRS)
  71.   115 CONTINUE
  72.       WRITE(*,*) ' ***  DEP **  '
  73.       DO 114 IA = 1,NSTRS
  74.         WRITE(*,1991) (DEP(IA,IC),IC=1,NSTRS)
  75.   114 CONTINUE
  76.       ENDIF
  77. 1991  FORMAT(18(1X,E12.5))
  78. C-------------------------------------------------------------------
  79.       RETURN
  80.       END
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  

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