Télécharger la plaquette Cast3M


fr en

Télécharger avanp.eso

Retour à la liste

Numérotation des lignes : 

avanp

  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