Télécharger estito.eso

Retour à la liste

Numérotation des lignes :

  1. C ESTITO SOURCE CHAT 05/01/12 23:45:42 5004
  2. SUBROUTINE ESTITO(SIG,NSTRS,VAR,NVARI,FKX,NKX,NPTEUR,NC,
  3. & XMAT,NCOMAT,T,TAUX,TRUC,NCOURB)
  4. C
  5. C----------------------------------------------------------------------
  6. C Objet: Estimation du pas de temps par integration de la loi
  7. C d'evolution du dommage a contr. cte
  8. C----------------------------------------------------------------------
  9. C Entree: SIG(NSTRS) tenseur des contraintes
  10. C VAR(NVARI) tableau contenant les variables internes
  11. C pilotant les equations
  12. C VAR(1)=p ; VAR(2)=r ;VAR(3)=D
  13. C XMAT(NCOMAT) tableau des parametres scalaires du materiau
  14. C pouvant dependre de la temperature
  15. C XMAT(1)=YOUNG ;XMAT(2)=XNU ;XMAT(3)=n
  16. C XMAT(4)=m ;XMAT(5)=KK; XMAT(6)=ALPHA
  17. C XMAT(7)=BETA ;XMAT(8)=r; XMAT(9)=A ;
  18. C XMAT(10)=EVOL ;XMAT(11)=ALPHAT
  19. C XMAT(12)=RHO; XMAT(13)=SIGY
  20. C FKX(NKX) tableau de la fonction k(X,T) contenant
  21. C les courbes de k a differentes temperatures T1,T2,T3,...
  22. C NPTEUR(NC) tableau de pointeurs sur les differentes courbes
  23. C de la fonction k
  24. C T temperature
  25. C------------------------------------------------------------------------
  26. C
  27. C------------------------------------------------------------------------
  28. C Sortie: TAUX pas de temps estime
  29. C----------------------------------------------------------------------
  30. C
  31. IMPLICIT INTEGER(I-N)
  32. IMPLICIT REAL*8(A-H,O-Z)
  33. DIMENSION SIG(*),VAR(*),XMAT(*)
  34. DIMENSION FKX(*),NPTEUR(*),TRUC(*)
  35. DIMENSION XX(6),A(3,3),D(3),S(3,3)
  36. C----------------------------------------------------------------------
  37. C CALCULATE J0,J1,J2
  38. C J0: THE BIGEST PRINCIPAL STRESS
  39. C J1: MEAN STRESS
  40. C J2: CONTR. DE VON MISES
  41. C----------------------------------------------------------------------
  42. AJ1 = SIG(1)+SIG(2)+SIG(3)
  43. AJ1 = AJ1/3.0D0
  44. DO 30 I=1,3
  45. 30 XX(I) = SIG(I)-AJ1
  46. DO 40 I=4,NSTRS
  47. 40 XX(I) = SIG(I)
  48. AJ2 = PROCON(XX,XX,NSTRS)
  49. AJ2 = SQRT(1.5D0*AJ2)
  50.  
  51. CALL ZERO(A,3,3)
  52. DO 50 I=1,3
  53. 50 A(I,I) = SIG(I)
  54. A(1,2) = SIG(4)
  55. IF(NSTRS.NE.6) GOTO 54
  56. A(1,3) = SIG(5)
  57. A(2,3) = SIG(6)
  58. 54 CONTINUE
  59. A(2,1) = A(1,2)
  60. A(3,1) = A(1,3)
  61. A(3,2) = A(2,3)
  62. IF (NSTRS.EQ.4) CALL JACOB3(A,2,D,S)
  63. IF (NSTRS.EQ.6) CALL JACOB3(A,3,D,S)
  64. C
  65. C---------------------------------------------------------------------
  66. C CALCULATE PRINCIPLE STRESS : D(*)
  67. C---------------------------------------------------------------------
  68. AJ0 = DMAX1(D(1),D(2),D(3))
  69. C
  70. C---------------------------------------------------------------------
  71. C CALCULATE VARIABLE X
  72. C---------------------------------------------------------------------
  73. C
  74. ALFA = XMAT(6)
  75. BETA = XMAT(7)
  76. X = (ALFA*AJ0)+(BETA*AJ1)+((1.0D0-ALFA-BETA)*AJ2)
  77. C
  78. C---------------------------------------------------------------------
  79. C CALCULATE PARAMETER k(X,T)=AK
  80. C---------------------------------------------------------------------
  81. C
  82. CALL VALPAR(FKX,NPTEUR,NKX,NC,X,AK,AKV,T,TIN,TSU,TRUC,NCOURB)
  83. C--------------------------------------------------------------------
  84. C Estimation du pas de temps par integration de la loi
  85. C d'evolution du dommage a contr. cte
  86. C---------------------------------------------------------------------
  87. EA=XMAT(9)
  88. ER=XMAT(8)
  89. TAUX = 1.0D10
  90. IF(VAR(3).GE.0.3.AND.X.GT.0.0D0) THEN
  91. AK1 = AK + 1
  92. TAUX = (1.0D0-VAR(3))**AK1/AK1/(X/EA)**ER
  93. ENDIF
  94. C
  95. RETURN
  96. END
  97.  
  98.  
  99.  

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