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  26/08/24    21:15:57     12622          
  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 11 I=1,NSTRS
  32.       DO 5 J=1,NSTRS
  33.         AB(I)=AB(I)+AH(I,J)*DFSIG(J)
  34.     5 CONTINUE
  35.    11 CONTINUE
  36. C
  37.       DO 12 I=1,NSTRS
  38.       DO 6 J=1,NSTRS
  39.        AC(I)=AC(I)+AH(J,I)*DFSIG(J)
  40.     6 CONTINUE
  41.    12 CONTINUE
  42. C
  43.       DO 13 I=1,NSTRS
  44.       DO 7 J=1,NSTRS
  45.        DP(I,J)=DP(I,J)+AB(I)*AC(J)
  46.     7 CONTINUE
  47.    13 CONTINUE
  48. C
  49.       H1=0.D0
  50.       DO 8 I=1,NSTRS
  51.         H1=H1+DFSIG(I)*AB(I)
  52.     8 CONTINUE
  53. C
  54.       HA=H1+PAEC
  55.       DO 14 I=1,NSTRS
  56.       DO 9 J=1,NSTRS
  57.        DP(I,J)=DP(I,J)/HA
  58.        DEP(I,J)=AH(I,J)-DP(I,J)
  59.     9 CONTINUE
  60.    14 CONTINUE
  61. C
  62.       DO 15 I=1,NSTRS
  63.       DO 10 J=1,NSTRS
  64.         IF (ABS(DEP(I,J)).LT.1.D-5) THEN
  65.           DEP(I,J)=0.D0
  66.         ENDIF
  67.    10 CONTINUE
  68.    15 CONTINUE
  69. C-------------------------------------------------------------------
  70.       RETURN
  71.       END
  72.  
  73.  
  74.  
  75.  
  76.  
  77.  
  78.  

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