Télécharger la plaquette Cast3M


fr en

Télécharger estit1.eso

Retour à la liste

Numérotation des lignes : 

estit1

  1. C ESTIT1    SOURCE    CB215821  17/11/30    21:16:04     9639           
  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.  
  82.  

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