Télécharger intfib.eso

Retour à la liste

Numérotation des lignes :

  1. C INTFIB SOURCE CB215821 16/04/21 21:17:12 8920
  2. C sub
  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 INTFIB(XMAT,XSECT,DEPST,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),XSECT,DEPST
  19. REAL*8 SIG0(3),SIGF,EINT
  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 = VAR0(17)
  87.  
  88. XLAM = VAR0(21)
  89. EPTO = VAR0(22)
  90. INDITER = nint(VAR0(23))
  91. C
  92. C-----AVANCE EN DEFORMATION
  93. C
  94. EPTO = (EPTO + DEPST)*LCCOEF*(1.0D0-XLAM)
  95. C
  96. C-----PREDICTION ELASTIQUE
  97. C
  98. SIGMA = 2.0D0*GCOEFF*EPTO*(1.0D0-ENDI) + 2.0D0*GCOEFF*(EPTO-EPSF)*
  99. & ENDI
  100. XENER = GCOEFF*EPTO*EPTO
  101. XENE0 = GCOEFF*EPS0CO*EPS0CO
  102. SIGMP = 2.0D0*GCOEFF*ENDI*(EPTO-EPSF)
  103. C
  104. C-----CALCUL DE YR POUR L ACCROISSEMENT DE LA ROUILLE
  105. C
  106. CALL APPFIB(TCI,XENRR,15)
  107. C
  108. C-----CALCUL DU CRITERE D ENDOMMAGEMENT
  109. C
  110. CRIT1 = XENER - (XENE0 + ZVAR + XENRR)
  111.  
  112. IF (CRIT1.GT.0.0D0) THEN
  113. ENDI = 1.0D0-1.0D0/(1.0D0+AD*(XENER-(XENE0+XENRR)))
  114. ZVAR = XENER-XENE0-XENRR
  115. END IF
  116. C
  117. C-----CALCUL DU CRITERE POUR LE FROTTEMENT
  118. C
  119. CRIT2 = ABS(SIGMP-XVAR)
  120.  
  121. IF (CRIT2.GT.0.0D0) THEN
  122.  
  123. CRIT20 = ABS(CRIT2)
  124.  
  125. DO I1=1,ITERMAX
  126.  
  127. CRIT2 = ABS(SIGMP-XVAR)
  128.  
  129. IF (ABS(CRIT2/CRIT20).LE.1.0D-5) THEN
  130. GOTO 10
  131. END IF
  132.  
  133. DFSIGF = ABS(SIGMP-XVAR)/(SIGMP-XVAR)
  134. DFXVAR = -ABS(SIGMP-XVAR)/(SIGMP-XVAR)
  135. DGSIGF = ABS(SIGMP-XVAR)/(SIGMP-XVAR)
  136. DGXVAR = -ABS(SIGMP-XVAR)/(SIGMP-XVAR) + ACOEFF*XVAR
  137.  
  138. DLAMPI = CRIT2/(ENDI*2.0D0*GCOEFF*DFSIGF*DGSIGF+
  139. & GAMCOE*DFXVAR*DGXVAR)
  140.  
  141. XVAR = XVAR - GAMCOE*DLAMPI*DGXVAR
  142. SIGMP = SIGMP - 2.0D0*GCOEFF*ENDI*DGSIGF*DLAMPI
  143.  
  144. END DO
  145.  
  146. 10 CONTINUE
  147.  
  148. ALPH = XVAR/GAMCOE
  149. EPSF = EPTO - SIGMP/(2.0D0*ENDI*GCOEFF)
  150.  
  151. END IF
  152. C
  153. C-----ON REAJUSTE LES CONTRAINTES
  154. C
  155. SIGMA = 2.0D0*GCOEFF*EPTO*(1.0D0-ENDI) + 2.0D0*GCOEFF*(EPTO-EPSF)*
  156. & ENDI
  157. C
  158. C-----LES SORTIES
  159. C
  160. VARF(11) = ENDI
  161. VARF(12) = ALPH
  162. VARF(13) = EPSF
  163. VARF(14) = XVAR
  164. VARF(15) = ZVAR
  165. VARF(16) = EINT
  166. VARF(17) = SIGMA
  167.  
  168. VARF( 7) = EPTO
  169. VARF( 8) = SIGMA
  170.  
  171. C
  172. C-----LES CONTRAINTES EN SORTIE
  173. C
  174. SIGF = SIGMA
  175.  
  176. RETURN
  177. END
  178.  
  179.  
  180.  
  181.  
  182.  
  183.  
  184.  

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