Télécharger poudra.eso

Retour à la liste

Numérotation des lignes :

  1. C POUDRA SOURCE PV 09/12/22 21:15:02 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. C Modification ChD nouvelles equations de c et f selon
  10. C les equations de DA pour sa these.
  11. C
  12. C Modification LS nouvelle loi viscoplastique a trois parametres
  13. C Lois viscoplastiques d Abouaf et 3S sont disponibles ...
  14. C
  15. SUBROUTINE POUDRA(SIG,EPSV,VAR,EPSVPT,VARPT,XMAT,NSTRS,
  16. & NVARI,NCOMAT,KERRE)
  17. C
  18. IMPLICIT INTEGER(I-N)
  19. IMPLICIT REAL*8(A-H,O-Z)
  20. -INC CCOPTIO
  21. DIMENSION SIG(*),EPSV(*),VAR(*),EPSVPT(*),VARPT(*)
  22. DIMENSION XMAT(*)
  23. DIMENSION XX(6)
  24. C
  25. C---------------------------------------------------------------------
  26. C ====> Connaissant l'etat, c'est a dire SIG, EPSV et VAR on
  27. C calcule les derivees EPSVPT et VARPT(1) en utilisant les
  28. C constantes XMAT
  29. C----- calcul des invariants S1=Tr(sigma)
  30. C S2C=S2*S2=Tr(Dev(sigma):Dev(sigma))
  31. C avec Dev(sigma)=sigma-1/3Tr(sigma)
  32. C
  33. S1 = SIG(1)+SIG(2)+SIG(3)
  34. C
  35. C Modif Ch.D du 20/10/94 : Cf DEA Canel page A18
  36. C assure tr(XX) = 0
  37. DO I=1,2
  38. XX(I) = SIG(I)-S1/3.0D0
  39. END DO
  40. XX(3) = -XX(1)-XX(2)
  41. C Fin de modification Ch.D du 20/10/94
  42. C
  43. DO I=4,NSTRS,1
  44. XX(I) = SIG(I)
  45. END DO
  46. S2C = PROCON(XX,XX,NSTRS)
  47. C
  48. C----- calcul de F et C : fonctions de rho
  49. RHOR = VAR(NVARI)
  50. A = XMAT(6)
  51. RN = XMAT(5)
  52. QSRT = XMAT(7)
  53. F0 = XMAT(8)
  54. F1 = XMAT(9)
  55. F2 = XMAT(10)
  56. F3 = XMAT(11)
  57. F4 = XMAT(12)
  58. F5 = XMAT(13)
  59. C0 = XMAT(14)
  60. C1 = XMAT(15)
  61. C2 = XMAT(16)
  62. C3 = XMAT(17)
  63. C4 = XMAT(18)
  64. C5 = XMAT(19)
  65. RMIN = F5+1.0D-03
  66. IF (RHOR.LT.RMIN) THEN
  67. RHOR1 = RMIN
  68. ELSE
  69. RHOR1 = RHOR
  70. ENDIF
  71. C
  72. C----- Modification L3S du PROGRAMME :
  73. C Integration de la loi a 3 parametres
  74. C Le 15/10/1997 par Laurent SANCHEZ ...
  75. C COEF : Coefficient F(Rhor)
  76. C COEC : Coefficient C(Rhor)
  77. C COEE : Coefficient E(Rhor)
  78. C COSS : Coefficient S(Rhor)
  79. C
  80. C COF1 : Coefficient f1(Rhor)
  81. C COF2 : Coefficient f2(Rhor)
  82. C COF3 : Coefficient f3(Rhor)
  83. C
  84. IF ((1.0D0-RHOR1).LT.1.0D-06) THEN
  85. COEF = 0.0D0
  86. COEC = 1.0D0
  87. ELSE
  88. COEF = F0*(((1.D0-RHOR1)/(RHOR1-F5))**(F1*RHOR1+F2))
  89. COEC = 1.0D0+C0*(((1.D0-RHOR1)/(RHOR1-F5))**(C1*RHOR1+C2))
  90. ENDIF
  91. IF (COEF.GT.(C5)) COEF = C5
  92. IF (COEC.GT.(C5)) COEC = C5
  93. C
  94. C----- calcul de sigma equivalent :
  95. SIGEQ = SQRT(COEF*S1*S1+1.5D0*COEC*S2C)
  96. C
  97. C----- calcul de EPSVPT en fonction du parametre A :
  98. C A : INT(F4) = 0
  99. C Log10(A) : INT(F4) = 1
  100. C loi a 3 parametres : INT(F4) = 2
  101. C
  102. *pv pour être sur d'avoir la bonne valeur, on augmente un peu f4
  103. NF4=INT(F4+0.1)
  104. C
  105. C----- Loi d'Abouaf :
  106. IF (NF4.EQ.1) THEN
  107. DO I=1,NSTRS,1
  108. B = 0.D0
  109. IF (I.LE.3) B=1.D0
  110. EPSVPT(I) = 1.D0*(10.D0**A)*(SIGEQ**(RN-1.D0))*RHOR
  111. & *(COEF*S1*B+1.5D0*COEC*XX(I))
  112. END DO
  113. C
  114. C----- Loi viscoplastique 3S a trois parametres :
  115. ELSE IF (NF4.EQ.2) THEN
  116. IF ((1.D0-RHOR1).LT.1.D-6) THEN
  117. COEE = 0.0D0
  118. COES = 1.0D0
  119. COF1 = 0.D0
  120. COF2 = 1.D0
  121. COF3 = 0.D0
  122. ELSE
  123. COEE = (3.D0*COEF)/(COEF+COEC)
  124. COSS = (COEF+COEC)**(-(RN+1.D0)/(2.D0*RN))
  125. COF1 = (2.D0*COEF/3.D0)**((RN+1.D0)/2.D0)
  126. COF2 = (1.D0-COEE)*(COSS**(-RN))
  127. & + 2D0*( (1.5D0)**((RN+1.D0)/2.D0) )*COF1
  128. COF3 = COEE*(COSS**(-RN))*((2.D0/3.D0)**((RN+5.D0)/4.D0))
  129. & - 2.D0*( (1.5D0)**((RN+1.D0)/4.D0) )*COF1
  130. ENDIF
  131. DO I=1,NSTRS,1
  132. B = 0.D0
  133. IF (I.LE.3) B=1.D0
  134. EPSVPT(I) = (3.D0/2.D0)**((RN+1.D0)/2.D0)*EXP(-QSRT)
  135. & *(10.D0**A)*RHOR *((COF1*S1*(S1**2.D0)**((RN-1.D0)/2.D0)
  136. & +0.5D0*COF3*S1*S2C*(S1**2.D0*S2C)**((RN-3.D0)/4.D0))
  137. & *B+(COF2*S2C**(RN/2.D0)+0.5D0*COF3*S1**2.D0*S2C**0.5D0
  138. & *(S1**2.D0*S2C)**((RN-3.D0)/4.D0))*XX(I))
  139. END DO
  140. C
  141. C----- Loi viscoplastique par defaut
  142. ELSE
  143. DO I=1,NSTRS,1
  144. B = 0.D0
  145. IF (I.LE.3) B=1.D0
  146. EPSVPT(I) = 1.D0*A*(SIGEQ**(RN-1.D0))*EXP(-QSRT)*RHOR
  147. & *(COEF*S1*B+1.5D0*COEC*XX(I))
  148. END DO
  149. ENDIF
  150. C
  151. C----- calcul de VARPT(1) derivee de EPSE
  152. PPT = SQRT(2.D0/3.D0*PROCON(EPSVPT,EPSVPT,NSTRS))
  153. VARPT(1) = PPT
  154. C---------------------------------------------------------------------
  155. C
  156. RETURN
  157. END
  158.  
  159.  
  160.  
  161.  
  162.  

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