Télécharger estit1.eso

Retour à la liste

Numérotation des lignes :

  1. C ESTIT1 SOURCE AM 08/12/19 21:16:36 6237
  2. SUBROUTINE ESTIT1(SIG,NSTRS,VAR,NVARI,XMAT,NCOMAT,T,TAUX)
  3. C
  4. C----------------------------------------------------------------------
  5. C Objet: Estimation du pas de temps par integration de la loi
  6. C d'evolution du dommage a contr. cte
  7. C----------------------------------------------------------------------
  8. C Entree: SIG(NSTRS) tenseur des contraintes
  9. C VAR(NVARI) tableau contenant les variables internes
  10. C pilotant les equations
  11. C VAR(1)=p ; VAR(2)=D11 ; VAR(3)=D22; VAR(4)=D33;
  12. C VAR(5)=D12; VAR(6)=D13; VAR(7)=D23; VAR(8)=DMAX
  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)=q ;XMAT(11)=ALPHAT
  19. C XMAT(12)=RHO; XMAT(13)=SIGY
  20. C T temperature
  21. C------------------------------------------------------------------------
  22. C
  23. C------------------------------------------------------------------------
  24. C Sortie: TAUX pas de temps estime
  25. C----------------------------------------------------------------------
  26. C
  27. IMPLICIT INTEGER(I-N)
  28. IMPLICIT REAL*8(A-H,O-Z)
  29. DIMENSION SIG(*),VAR(*),XMAT(*)
  30. DIMENSION XX(6),A(3,3),D(3),S(3,3)
  31. C
  32. C----------------------------------------------------------------------
  33. C CALCUL DE J1,J2
  34. C J1: CONTRAINTE MOYENNE
  35. C J2: CONTRAINTE DE VON MISES
  36. C----------------------------------------------------------------------
  37. C
  38. AJ1 = SIG(1)+SIG(2)+SIG(3)
  39. AJ1 = AJ1/3.0D0
  40. DO I=1,3
  41. XX(I) = SIG(I)-AJ1
  42. ENDDO
  43. DO I=4,NSTRS
  44. XX(I) = SIG(I)
  45. ENDDO
  46. AJ2 = PROCON(XX,XX,NSTRS)
  47. AJ2 = SQRT(1.5D0*AJ2)
  48. C
  49. C---------------------------------------------------------------------
  50. C CALCUL DE LA VALEUR PROPRE D'ENDOMMAGEMENT MAX. : D(*)
  51. C---------------------------------------------------------------------
  52. C
  53. DMAX = VAR(8)
  54. C
  55. C---------------------------------------------------------------------
  56. C CALCUL DE LA VARIABLE X
  57. C---------------------------------------------------------------------
  58. C
  59. ALF2 = XMAT(6)
  60. X = (ALF2*AJ1)+(AJ2)
  61. C
  62. C--------------------------------------------------------------------
  63. C Estimation du pas de temps par integration de la loi
  64. C d'evolution du dommage a contr. cte
  65. C---------------------------------------------------------------------
  66. C
  67. EA=XMAT(9)
  68. ER=XMAT(8)
  69. EQ=XMAT(10)
  70. TAUX = 1.0D10
  71. IF(DMAX.GE.0.3.AND.X.GT.0.0D0) THEN
  72. EQ1 = EQ + 1
  73. TAUX = (1.0D0-DMAX)**EQ1/EQ1/(X/EA)**ER
  74. ENDIF
  75. C
  76. RETURN
  77. END
  78.  
  79.  
  80.  
  81.  

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