Télécharger intfic.eso

Retour à la liste

Numérotation des lignes :

  1. C INTFIC SOURCE CB215821 16/04/21 21:17:12 8920
  2. C sub intfibc
  3. C====&===1=========2=========3=========4=========5=========6=========7==
  4. C Commentaires : Subroutine permettant de mettre en oeuvre le
  5. C modele RICINT pour une formulation à fibre
  6. C Auteurs : B. Richard (doctorant)
  7. C Date : Février 2008
  8. C====&===1=========2=========3=========4=========5=========6=========7==
  9. SUBROUTINE INTFIC(XMAT,DEPS,SIG0,VAR0,SIGF,VARF)
  10. C
  11. C-----DECLARATION GENERALE
  12. C
  13. IMPLICIT INTEGER(I-N)
  14. IMPLICIT REAL*8(A-H,O-Z)
  15. C
  16. C-----DECLARATION DES VARIABLES
  17. C
  18. REAL*8 XMAT(18),VAR0(25),VARF(25),DEPST
  19. REAL*8 SIG0(3),SIGF(3),EINT,DEPS(3)
  20. REAL*8 GCOEFF,AD,ACOEFF,GAMCOE,LCCOEF,EPS0CO,TCI
  21. REAL*8 ENDI,ALPH,EPSF,XVAR,ZVAR,XLAM,EPTO,SIGMA
  22. INTEGER INDITER,I1,ITERMAX
  23. REAL*8 XENER,SIGMP,CRIT1,XENRR,CRIT2,DLAMPI,CRIT20
  24. REAL*8 DFSIGF,DFXVAR,DGSIGF,DGXVAR,XENE0
  25. C
  26. C-----LES PARAMETRES
  27. C
  28. ITERMAX = 100000
  29. C
  30. C-----NOMENCLATURE DES PARAMETRES MATERIAUX (XMAT) ET VARIABLES INTERNES
  31. C (VAR0 ET VARF)
  32. C
  33.  
  34. C ---------------------------
  35. C | INTERFACE A/B CORRODEE |
  36. C ---------------------------
  37. C XMAT(11) GCOEFF
  38. C XMAT(12) AD
  39. C XMAT(13) ACOEFF
  40. C XMAT(14) GAMCOE
  41. C XMAT(15) LCCOEF
  42. C XMAT(16) EPS0CO
  43. C XMAT(17) TCI
  44. C XMAT(18) CALA
  45. C
  46. C ---------------------------
  47. C | INTERFACE A/B CORRODEE |
  48. C ---------------------------
  49. C VAR0(11) VARF(11) ENDI
  50. C VAR0(12) VARF(12) ALPH
  51. C VAR0(13) VARF(13) EPSF
  52. C VAR0(14) VARF(14) XVAR
  53. C VAR0(15) VARF(15) ZVAR
  54. C VAR0(16) VARF(16) EINT
  55. C VAR0(17) VARF(17) SIGF
  56. C VAR0(18) VARF(18)
  57. C VAR0(19) VARF(19)
  58. C VAR0(20) VARF(20)
  59. C
  60. C ---------------------------
  61. C | AUTRES |
  62. C ---------------------------
  63. C VAR0(21) VARF(21) XLAM
  64. C VAR0(22) VARF(22) EPTO
  65. C VAR0(23) VARF(23) INDITER
  66.  
  67. C
  68. C-----LES PARAMETRES MATERIAUX
  69. C
  70. GCOEFF = XMAT(11)
  71. AD = XMAT(12)
  72. ACOEFF = XMAT(13)
  73. GAMCOE = XMAT(14)
  74. LCCOEF = XMAT(15)
  75. EPS0CO = XMAT(16)
  76. TCI = XMAT(17)
  77. C
  78. C-----LES VARIABLES INTERNES
  79. C
  80. ENDI = VAR0(11)
  81. ALPH = VAR0(12)
  82. EPSF = VAR0(13)
  83. XVAR = VAR0(14)
  84. ZVAR = VAR0(15)
  85. EINT = VAR0(16)
  86. SIGF(1) = VAR0(17)
  87.  
  88. XLAM = 0.0D0
  89. EPTO = VAR0(22)
  90. INDITER = nint(VAR0(23))
  91. C
  92. C-----AVANCE EN DEFORMATION
  93. C
  94. DEPST= DEPS(1)
  95. EPTO = (EPTO + DEPST)*LCCOEF*(1.0D0-XLAM)
  96. C
  97. C-----PREDICTION ELASTIQUE
  98. C
  99. SIGMA = 2.0D0*GCOEFF*EPTO*(1.0D0-ENDI) + 2.0D0*GCOEFF*(EPTO-EPSF)*
  100. & ENDI
  101. XENER = GCOEFF*EPTO*EPTO
  102. XENE0 = GCOEFF*EPS0CO*EPS0CO
  103. SIGMP = 2.0D0*GCOEFF*ENDI*(EPTO-EPSF)
  104. C
  105. C-----CALCUL DE YR POUR L ACCROISSEMENT DE LA ROUILLE
  106. C
  107. CALL APPFIB(TCI,XENRR,15)
  108. C
  109. C-----CALCUL DU CRITERE D ENDOMMAGEMENT
  110. C
  111. CRIT1 = XENER - (XENE0 + ZVAR + XENRR)
  112.  
  113. IF (CRIT1.GT.0.0D0) THEN
  114. ENDI = 1.0D0-1.0D0/(1.0D0+AD*(XENER-(XENE0+XENRR)))
  115. ZVAR = XENER-XENE0-XENRR
  116. END IF
  117. C
  118. C-----CALCUL DU CRITERE POUR LE FROTTEMENT
  119. C
  120. CRIT2 = ABS(SIGMP-XVAR)
  121.  
  122. IF (CRIT2.GT.0.0D0) THEN
  123.  
  124. CRIT20 = ABS(CRIT2)
  125.  
  126. DO I1=1,ITERMAX
  127.  
  128. CRIT2 = ABS(SIGMP-XVAR)
  129.  
  130. IF (ABS(CRIT2/CRIT20).LE.1.0D-5) THEN
  131. GOTO 10
  132. END IF
  133.  
  134. DFSIGF = ABS(SIGMP-XVAR)/(SIGMP-XVAR)
  135. DFXVAR = -ABS(SIGMP-XVAR)/(SIGMP-XVAR)
  136. DGSIGF = ABS(SIGMP-XVAR)/(SIGMP-XVAR)
  137. DGXVAR = -ABS(SIGMP-XVAR)/(SIGMP-XVAR) + ACOEFF*XVAR
  138.  
  139. DLAMPI = CRIT2/(ENDI*2.0D0*GCOEFF*DFSIGF*DGSIGF+
  140. & GAMCOE*DFXVAR*DGXVAR)
  141.  
  142. XVAR = XVAR - GAMCOE*DLAMPI*DGXVAR
  143. SIGMP = SIGMP - 2.0D0*GCOEFF*ENDI*DGSIGF*DLAMPI
  144.  
  145. END DO
  146.  
  147. 10 CONTINUE
  148.  
  149. ALPH = XVAR/GAMCOE
  150. EPSF = EPTO - SIGMP/(2.0D0*ENDI*GCOEFF)
  151.  
  152. END IF
  153. C
  154. C-----ON REAJUSTE LES CONTRAINTES
  155. C
  156. SIGMA = 2.0D0*GCOEFF*EPTO*(1.0D0-ENDI) + 2.0D0*GCOEFF*(EPTO-EPSF)*
  157. & ENDI
  158. C
  159. C-----LES SORTIES
  160. C
  161. VARF(11) = ENDI
  162. VARF(12) = ALPH
  163. VARF(13) = EPSF
  164. VARF(14) = XVAR
  165. VARF(15) = ZVAR
  166. VARF(16) = EINT
  167. VARF(17) = SIGMA
  168.  
  169. VARF( 7) = EPTO
  170. VARF( 8) = SIGMA
  171. VARF(22) = EPTO
  172. VARF(23) = INDITER
  173.  
  174. C
  175. C-----LES CONTRAINTES EN SORTIE
  176. C
  177. SIGF(1) = SIGMA
  178.  
  179. RETURN
  180. END
  181.  
  182.  
  183.  
  184.  
  185.  
  186.  
  187.  

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