Télécharger poudrb.eso

Retour à la liste

Numérotation des lignes :

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

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