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 :
© Cast3M 2003 - All rights reserved.
Disclaimer