Télécharger vispl1.eso

Retour à la liste

Numérotation des lignes :

vispl1
  1. C VISPL1 SOURCE CHAT 05/01/13 04:07:26 5004
  2. C VISPL1 SOURCE AM 00/12/21 21:26:00 4062
  3. SUBROUTINE VISPL1(SIGP,DSIGT,NSTRS,SIGF,SIGRV,DSTRN,
  4. & DEPS1,DEPS2,BETJEF,VISCO,NECH0,NECH1)
  5. C
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. DIMENSION SIGP(4),DSIGT(4),DSTRN(4),D(4,4)
  9. DIMENSION DSIG(4),DSIG2(4),DSIG3(4),DP(4,4)
  10. DIMENSION SIGF(4),SIGRV(4),DF(4)
  11. C
  12. *
  13. SEGMENT BETJEF
  14. REAL*8 AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,CAR,ETA,TDEF,
  15. & TCON,DPS1,DPS2,TETA,PDT,TP0
  16. INTEGER ICT,ICC,IMOD,IVIS,ITER,
  17. & ISIM,IBB1,IGAU1,IZON
  18. ENDSEGMENT
  19. SEGMENT VISCO
  20. REAL*8 DPSITV1,DPSITV2,SIGV1,SIGV2,ENDV
  21. ENDSEGMENT
  22. SEGMENT NECH0
  23. REAL*8 DT,DC,ALFG,S0,ENDO
  24. ENDSEGMENT
  25. SEGMENT NECH1
  26. REAL*8 ENDL
  27. ENDSEGMENT
  28. C
  29. C
  30. C ********************************************************************
  31. C ** COMPORTEMENT VISCOPLASTIQUE_ENDOMMAGEABLE **
  32. C ** DE TYPE DUVAUT LIONS **
  33. C ********************************************************************
  34. C
  35. IPAS=0
  36. DAM = 0.D0
  37. CALL ZERO(DSIG2,4,1)
  38. CALL ZERO(DSIG3,4,1)
  39. DPSTV1 = DPSITV1
  40. DPSTV2 = DPSITV2
  41. ENDV0 = ENDV
  42. ENDO0 = ENDO
  43. C
  44. C ---------------- calcul des termes viscoplastiques -----------------
  45. C
  46. GAMA = 0.5
  47. A = (ETA)/(ETA+GAMA*PDT)
  48. B = (GAMA*PDT)/ETA
  49. C = PDT
  50. G=EXP(-PDT/ETA)
  51. IF (C.EQ.0.D0) THEN
  52. WRITE(*,*)'l increment de temps est de longueur nulle'
  53. stop
  54. ENDIF
  55. C
  56. C -------------------- creation matrice de young -------------------
  57. C
  58. IF (IMOD.EQ.1.OR.IMOD.EQ.3) THEN
  59. CALL CREMAT(D,YOUN,XNU,3,-2)
  60. ENDIF
  61. IF (IMOD.EQ.2.OR.IMOD.EQ.4) THEN
  62. CALL CREMAT(D,YOUN,XNU,4,-1)
  63. ENDIF
  64. C
  65. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  66. C MISE A JOUR DES VARIABLES D'ECROUISSAGES
  67. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  68. C
  69. DPSTV1=DPSTV1*G+(1.D0-G)*DEPS1
  70. CALL UNICO1(DPSTV1,PAEC,1,SIGV1,BETJEF,NECH0,NECH1)
  71. DPSTV2=DPSTV2*G+(1.D0-G)*DEPS2
  72. CALL UNICO1(DPSTV2,PAEC,2,SIGV2,BETJEF,NECH0,NECH1)
  73. C
  74. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  75. C MISE A JOUR DE LA VARIABLE D'ENDOMMAGEMENT VISQUEUSE
  76. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  77. C
  78. CALL DAMAG1(DPSTV1,DPSTV2,ENDC,
  79. & ENDT,IFIS,IPLA,3,BETJEF,NECH0,NECH1)
  80. ENDV = 1.D0-((1.D0-ENDC)*(1.D0-ENDT))
  81. C
  82. C ------------ mise a jour contrainte ----------------------
  83. C ------ methode explicite
  84. CC DD = (ENDV - ENDV0) / PDT
  85. CC FF = 1. - ENDV0
  86. CC FK = -((1./ETA)+(DD/FF))
  87. CC DO 15 I=1,NSTRS
  88. CC DF(I) = FK*SIGRV(I)+FF*(DSIGT(I)/PDT)+(SIGP(I)/ETA)
  89. CC 15 CONTINUE
  90. C ------------ resolution numerique de l equa diff ---------------
  91. CC CALL ZERO(SIGF,4,1)
  92. CC DO 36 I=1,NSTRS
  93. CC SIGF(I)= SIGRV(I) + PDT * DF(I)
  94. CC 36 CONTINUE
  95. C
  96. DD = (ENDV - ENDV0)/(1. - ENDV)
  97. FK = (PDT/ETA)
  98. CC = 1./(1 + FK + DD)
  99. C
  100. C ------------ resolution numerique de l equa diff ---------------
  101. CALL ZERO(SIGF,4,1)
  102. DO 37 I=1,NSTRS
  103. SIGF(I)= CC *(SIGRV(I) + (1.-ENDV)*DSIGT(I) + FK*SIGP(I))
  104. 37 CONTINUE
  105. C
  106. C
  107. C ------- taux de deformation en norme --------
  108. C
  109. TDEF=0.D0
  110. D1=DSTRN(1)
  111. D2=DSTRN(2)
  112. D3=DSTRN(3)
  113. D4=DSTRN(4)
  114. TDEF=SQRT(D1*D1+D2*D2+D3*D3+D4*D4)/PDT
  115. C
  116. C ------- taux de contrainte en norme --------
  117. C
  118. TCON=0.D0
  119. DO 62 I=1,NSTRS
  120. DSIG(I) = (SIGF(I) - SIGRV(I))
  121. 62 CONTINUE
  122. S1=DSIG(1)
  123. S2=DSIG(2)
  124. S3=DSIG(3)
  125. S4=DSIG(4)
  126. TCON=SQRT(S1*S1+S2*S2+S3*S3+S4*S4)/PDT
  127. C
  128. CALL INDIC1(DPSTV1,DPSTV2,IFIS,IPLA,3,BETJEF,NECH0,NECH1)
  129. C
  130. DPSITV1 = DPSTV1
  131. DPSITV2 = DPSTV2
  132. RETURN
  133. END
  134.  
  135.  
  136.  
  137.  
  138.  

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