Télécharger credep.eso

Retour à la liste

Numérotation des lignes :

  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