Télécharger la plaquette Cast3M


fr en

Télécharger credep.eso

Retour à la liste

Numérotation des lignes : 

credep

  1. C CREDEP    SOURCE    CB215821  16/04/21    21:16:07     8920
  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-------------------------------------------------------------------
  25.       CALL ZERO(DEP,4,4)
  26.       CALL ZERO(DP,4,4)
  27.       CALL ZERO(AB,6,1)
  28.       CALL ZERO(AC,6,1)
  29. C-------------------------------------------------------------------
  30. C
  31.       DO 5 I=1,NSTRS
  32.       DO 5 J=1,NSTRS
  33.         AB(I)=AB(I)+AH(I,J)*DFSIG(J)
  34.     5 CONTINUE
  35. C
  36.       DO 6 I=1,NSTRS
  37.       DO 6 J=1,NSTRS
  38.        AC(I)=AC(I)+AH(J,I)*DFSIG(J)
  39.     6 CONTINUE
  40. C
  41.       DO 7 I=1,NSTRS
  42.       DO 7 J=1,NSTRS
  43.        DP(I,J)=DP(I,J)+AB(I)*AC(J)
  44.     7 CONTINUE
  45. C
  46.       H1=0.D0
  47.       DO 8 I=1,NSTRS
  48.         H1=H1+DFSIG(I)*AB(I)
  49.     8 CONTINUE
  50. C
  51.       HA=H1+PAEC
  52.       DO 9 I=1,NSTRS
  53.       DO 9 J=1,NSTRS
  54.        DP(I,J)=DP(I,J)/HA
  55.        DEP(I,J)=AH(I,J)-DP(I,J)
  56.     9 CONTINUE
  57. C
  58.       DO 10 I=1,NSTRS
  59.       DO 10 J=1,NSTRS
  60.         IF (ABS(DEP(I,J)).LT.1.D-5) THEN
  61.           DEP(I,J)=0.D0
  62.         ENDIF
  63.    10 CONTINUE
  64. C-------------------------------------------------------------------
  65.       RETURN
  66.       END
  67.  
  68.  
  69.  
  70.  
  71.  
  72.  

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