Télécharger avanp.eso

Retour à la liste

Numérotation des lignes :

  1. C AVANP SOURCE PV 15/04/10 21:15:01 8474
  2. C
  3. C PLG93Debut
  4. C
  5. SUBROUTINE AVANP(TAU,SIG,EPSV,VAR,SIG1,EPSV1,VAR1,DSPT,
  6. & EPSVPT,VARPT,XMAT,NSTRS,NVARI,IFOUR,NCOMAT)
  7. IMPLICIT INTEGER(I-N)
  8. IMPLICIT REAL*8(A-H,O-Z)
  9. DIMENSION SIG(*),EPSV(*),VAR(*),SIG1(*),EPSV1(*),VAR1(*)
  10. DIMENSION DSPT(*),EPSVPT(*),VARPT(*),XMAT(*)
  11. C ____________________________________________________________________
  12. C
  13. C THIS ROUTINE UPDATES THE STRAINS ,STRESSES & INTERNAL VARIABLES
  14. C Version modifiee de la procedure advan1 pour le modele
  15. C viscoplastique poudre.
  16. C ____________________________________________________________________
  17. C
  18. C connaissant l'etat c'est a dire SIG, EPSV et VAR ainsi que les
  19. C derivees EPSVPT et DSPT on calcule l'etat SIG1, EPSV1 et VAR1
  20. C au temps T + TAU
  21. C ____________________________________________________________________
  22. C
  23. -INC CCREEL
  24.  
  25. C STEP-1 : UPDATING STRAINS
  26. C ----------------------------
  27. YOU = XMAT(1)
  28. XNU = XMAT(2)
  29. DO 10 I=1,NSTRS
  30. EPSV1(I)=EPSV(I)+TAU*EPSVPT(I)
  31. 10 CONTINUE
  32. C
  33. C STEP-2 : UPDATING STESSES
  34. C ----------------------------
  35. IF(IFOUR.EQ.-2) THEN
  36. SIG(3)=0.0D0
  37. X2MU= YOU/(1.0D0+XNU)
  38. XCO = X2MU/(1.0D0-XNU)
  39. SIG1(1)=SIG(1)+TAU*(DSPT(1)-XCO*(EPSVPT(1)+XNU*EPSVPT(2)))
  40. SIG1(2)=SIG(2)+TAU*(DSPT(2)-XCO*(EPSVPT(2)+XNU*EPSVPT(1)))
  41. SIG1(4)=SIG(4)+TAU*(DSPT(4)-X2MU*EPSVPT(4))
  42. SIG1(3)=0.0D0
  43. ELSEIF (IFOUR.EQ.2.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-1.OR.
  44. & IFOUR.EQ.1.OR.IFOUR.EQ.-3) THEN
  45. TRACE=EPSVPT(1)+EPSVPT(2)+EPSVPT(3)
  46. X2MU = YOU / (1.0D0+XNU)
  47. ALMB = X2MU*XNU/(1.0D0-2.0D0*XNU)
  48. DO 20 I=1,NSTRS,1
  49. A = 0.0D0
  50. IF (I.LE.3) A=1.0D0
  51. SIG1(I) = SIG(I) + TAU*(DSPT(I)-ALMB*A*TRACE-X2MU*EPSVPT(I))
  52. 20 CONTINUE
  53. C= Modes de calcul 1D contraintes planes suivant z (DYCZ et GYCZ)
  54. ELSE IF (IFOUR.EQ.4.OR.IFOUR.EQ.8) THEN
  55. SIG(3) = 0.D0
  56. XCO = YOU/(1.D0-XNU*XNU)
  57. SIG1(1) = SIG(1) + TAU*(DSPT(1)-XCO*(EPSVPT(1)+XNU*EPSVPT(2)))
  58. SIG1(2) = SIG(2) + TAU*(DSPT(2)-XCO*(EPSVPT(2)+XNU*EPSVPT(1)))
  59. SIG1(3) = 0.D0
  60. C= Modes de calcul 1D contraintes planes suivant y (CYDZ et CYGZ) et
  61. C= 1D axisymetrique contraintes axiales nulles (AXCZ)
  62. ELSE IF (IFOUR.EQ.5.OR.IFOUR.EQ.10.OR.IFOUR.EQ.13) THEN
  63. SIG(2) = 0.D0
  64. XCO = YOU/(1.D0-XNU*XNU)
  65. SIG1(1) = SIG(1) + TAU*(DSPT(1)-XCO*(EPSVPT(1)+XNU*EPSVPT(3)))
  66. SIG1(2) = 0.D0
  67. SIG1(3) = SIG(3) + TAU*(DSPT(3)-XCO*(EPSVPT(3)+XNU*EPSVPT(1)))
  68. C= Modes de calcul 1D contraintes planes suivant y et z (CYCZ)
  69. ELSE IF (IFOUR.EQ.6) THEN
  70. SIG(2) = 0.D0
  71. SIG(3) = 0.D0
  72. SIG1(1) = SIG(1) + TAU*(DSPT(1)-YOU*EPSVPT(1))
  73. SIG1(2) = 0.D0
  74. SIG1(3) = 0.D0
  75. C= Autres modes de calcul 1D deformations planes (DYDZ GYDZ DYGZ GYGZ)
  76. C= 1D axisymetrique (AXDZ AXGZ) et 1D spherique
  77. ELSE IF (IFOUR.EQ. 3.OR.IFOUR.EQ. 7.OR.IFOUR.EQ. 9.OR.IFOUR.EQ.11
  78. . .OR.IFOUR.EQ.12.OR.IFOUR.EQ.14.OR.IFOUR.EQ.15) THEN
  79. X2MU = YOU/(1.D0+XNU)
  80. TRACE = EPSVPT(1)+EPSVPT(2)+EPSVPT(3)
  81. XCO = X2MU*XNU*TRACE/(1.D0-XNU)
  82. SIG1(1) = SIG(1) + TAU*(DSPT(1)-X2MU*EPSVPT(1)-XCO)
  83. SIG1(2) = SIG(2) + TAU*(DSPT(2)-X2MU*EPSVPT(2)-XCO)
  84. SIG1(3) = SIG(3) + TAU*(DSPT(3)-X2MU*EPSVPT(3)-XCO)
  85. ENDIF
  86. C
  87. C STEP-3 : UPDATING THE INTERNAL VARIABLES
  88. C------------------------------------------
  89. DO 30 I=1,NSTRS
  90. VAR1(1+I)=EPSV1(I)
  91. 30 CONTINUE
  92. C
  93. C majoration de la densite relative : RHOR=RHOR-TraceEPSV*RHOR
  94. C
  95. TRE = EPSVPT(1)+EPSVPT(2)+EPSVPT(3)
  96. taux=-tre*tau
  97. taux=max(min(log(xgrand),taux),log(xpetit))
  98. VAR1(NVARI) = VAR(NVARI) * EXP(taux)
  99. IF (VAR1(NVARI).GT.1D0) VAR1(NVARI) = 1.D0
  100. IF (VAR1(NVARI).ne.var1(nvari)) VAR1(NVARI) = 1.D0
  101. C
  102. C majoration de EPSE :
  103. C
  104. VAR1(1) = VAR(1)+TAU*VARPT(1)
  105. C
  106. C _____________________________________________________________________
  107. RETURN
  108. END
  109. C
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  

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