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