Télécharger depo.eso

Retour à la liste

Numérotation des lignes :

  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