Télécharger la plaquette Cast3M


fr en

Télécharger poudrb.eso

Retour à la liste

Numérotation des lignes : 

poudrb

  1. C POUDRB    SOURCE    PV        09/12/22    21:15:17     6589
  2. C
  3. C  PLGDebut : implantation des equations relatives au modele
  4. C  viscoplastique poudre.
  5. C
  6. C  Modification ChD division par des rÅels et non pas des entiers
  7. C  COEC doit etre plus grand que 1.0 et non pas 0.0
  8. C
  9.        SUBROUTINE POUDRB(SIG,EPSV,VAR,EPSVPT,VARPT,XMAT,NSTRS,
  10.      +                  NVARI,NCOMAT,KERRE)
  11. C
  12.       IMPLICIT INTEGER(I-N)
  13.        IMPLICIT REAL*8(A-H,O-Z)
  14. -INC PPARAM
  15. -INC CCOPTIO
  16.        DIMENSION SIG(*),EPSV(*),VAR(*),EPSVPT(*),VARPT(*)
  17.        DIMENSION XMAT(*)
  18.        DIMENSION XX(6)
  19. C----------------------------------------------------------------
  20. C  ====> Connaissant l'etat, c'est a dire SIG, EPSV et VAR on
  21. C    calcule les derivees EPSVPT et VARPT(1) en utilisant les
  22. C    constantes XMAT
  23. C----------------------------------------------------------------
  24. C
  25. C-------calcul des invariants S1=Tr(sigma)
  26. C                             S2C=S2*S2=Tr(Dev(sigma):Dev(sigma))
  27. C                        avec Dev(sigma)=sigma-1/3Tr(sigma)
  28.        S1 = SIG(1)+SIG(2)+SIG(3)
  29. C Modif Ch.D du 20/10/94 : Cf DEA Canel page A18 assure tr(XX) = 0
  30.        DO 30 I=1,2
  31.        XX(I) = SIG(I)-S1/3.0D0
  32.  30    CONTINUE
  33.        XX(3) = -XX(1)-XX(2)
  34. C Fin de modification Ch.D du 20/10/94
  35.        DO 40 I=4,NSTRS,1
  36.        XX(I) = SIG(I)
  37.  40    CONTINUE
  38.        S2C = PROCON(XX,XX,NSTRS)
  39. C
  40.        RHOR = VAR(NVARI)
  41.        A    = XMAT(6)
  42.        RN   = XMAT(5)
  43.        QSRT = XMAT(7)
  44.        F0 = XMAT(8)
  45.        F1 = XMAT(9)
  46.        F2 = XMAT(10)
  47.        F3 = XMAT(11)
  48.        F4 = XMAT(12)
  49.        F5 = XMAT(13)
  50.        C0 = XMAT(14)
  51.        C1 = XMAT(15)
  52.        C2 = XMAT(16)
  53.        C3 = XMAT(17)
  54.        C4 = XMAT(18)
  55.        C5 = XMAT(19)
  56.        G0 = XMAT(20)
  57.        G1 = XMAT(21)
  58.        G2 = XMAT(22)
  59.        G3 = XMAT(23)
  60.        G4 = XMAT(24)
  61.        G5 = XMAT(25)
  62.        G6 = XMAT(26)
  63. C
  64.        RMIN = F5+1.0D-03
  65.        IF (RHOR.LT.RMIN) THEN
  66.         RHOR1 = RMIN
  67.        ELSE
  68.         RHOR1 = RHOR
  69.        ENDIF
  70. C
  71. C
  72. C-------calcul de EPSVPT en fonction du parametre F4 :
  73. C------- F4 = 0      : loi ABOUAF avec taille de grain
  74. C------- F4 = 1      : loi de Stutz
  75. C
  76. *pv  pour être sur d'avoir la bonne valeur
  77.        NF4=INT(F4+0.1)+1
  78.        GOTO(10,20),NF4
  79. C
  80. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  81. C
  82. C        Loi Abouaf avec calcul de la taille de grain
  83. C
  84. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  85. C
  86. 10     CONTINUE
  87. C
  88. C-------calcul de F et C : fonctions de rho
  89. C
  90.        IF ((1.0D0-RHOR1).LT.(1.D-06)) THEN
  91.           COEF = 0.0
  92.           COEC = 1.0
  93.        ELSE
  94.           COEF = F0*(((1-RHOR1)/(RHOR1-F5))**(F1*RHOR1+F2))
  95.           COEC = 1.0D0+C0*(((1-RHOR1)/(RHOR1-F5))**(C1*RHOR1+C2))
  96.        ENDIF
  97.        IF (COEF.GT.(C5)) COEF = C5
  98.        IF (COEC.GT.(C5)) COEC = C5
  99. C
  100. C-------calcul de sigma equivalent :
  101. C
  102.        SIGEQ = SQRT(COEF*S1*S1+1.5D0*COEC*S2C)
  103.        DO 500 I=1,NSTRS,1
  104.         B = 0.D0
  105.         IF (I.LE.3) B=1.D0
  106.         EPSVPT(I) = 1.D0*(10**A)*(SIGEQ**(RN-1.D0))*RHOR
  107.      &              *(COEF*S1*B+1.5D0*COEC*XX(I))
  108. 500    CONTINUE
  109.        GOTO 1000
  110. C
  111. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  112. C
  113. C        Loi de Stutz
  114. C
  115. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  116. C
  117.  20    CONTINUE
  118.        COEH = F0*(((1.-RHOR)/(RHOR-F5))**F1)
  119.        COEE = C0*(RHOR-RHOR**2)
  120.        COES = G0-(G1*(1.-RHOR))+(G2*((1.-RHOR)**2))
  121.        COF1 = (2./3.)**((RN+1)/2)*COEH/3.
  122.        COF2 = (1.-COEE)*(COES**(-RN))+2.*((3./2.)**((RN+1)/2))*COF1
  123.        COF3 = COEE*(COES**(-RN))*((2./3.)**((RN+5)/4))
  124.      &        -2.*COF1*((3./2.)**((RN+1)/4))
  125. C
  126. C-------calcul de EPSVPT :
  127. C
  128. C
  129.        DO 600 I=1,NSTRS,1
  130.        B = 0.0
  131.          IF (I.LE.3) B=1.0
  132.          SIGEQ = (S1**2*S2C)**((RN-3)/4)
  133.          EPSVPT(I) = (3./2.)**((RN+1)/2)*A*EXP(-QSRT)*RHOR
  134.      &   *(3.*S1*(COF1*S1**(RN-1)+COF3/2.*S2C*SIGEQ)*B
  135.      &   +(COF2*S2C**((RN-1)/2)+COF3/2.*S1**2*SIGEQ)*XX(I))
  136.  600    CONTINUE
  137. C
  138. C-------calcul de VARPT(1) derivee de EPSE                                         C
  139. 1000   CONTINUE
  140.        PPT = SQRT(2.0/3.0*PROCON(EPSVPT,EPSVPT,NSTRS))
  141.        VARPT(1) = PPT
  142. C
  143. C--------on rentre
  144. C
  145.        RETURN
  146.        END
  147. C
  148.  
  149.  
  150.  
  151.  
  152.  

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