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 : 50 : CALL ZERO(A,3,3) 51 : DO 50 I=1,3 52 : 50 A(I,I) = SIG(I) 53 : A(1,2) = SIG(4) 54 : IF(NSTRS.NE.6) GOTO 54 55 : A(1,3) = SIG(5) 56 : A(2,3) = SIG(6) 57 : 54 CONTINUE 58 : A(2,1) = A(1,2) 59 : A(3,1) = A(1,3) 60 : A(3,2) = A(2,3) 60 : IF (NSTRS.EQ.4) CALL JACOB3(A,2,D,S) 60 : IF (NSTRS.EQ.6) CALL JACOB3(A,3,D,S) 61 : C 62 : C--------------------------------------------------------------------- 63 : C CALCULATE PRINCIPLE STRESS : D(*) 64 : C--------------------------------------------------------------------- 65 : AJ0 = DMAX1(D(1),D(2),D(3)) 66 : C 67 : C--------------------------------------------------------------------- 68 : C CALCULATE VARIABLE X 69 : C--------------------------------------------------------------------- 70 : C 71 : ALFA = XMAT(6) 72 : BETA = XMAT(7) 73 : X = (ALFA*AJ0)+(BETA*AJ1)+((1.0D0-ALFA-BETA)*AJ2) 74 : C 75 : C--------------------------------------------------------------------- 76 : C CALCULATE PARAMETER k(X,T)=AK 77 : C--------------------------------------------------------------------- 78 : C 78 : CALL VALPAR(FKX,NPTEUR,NKX,NC,X,AK,AKV,T,TIN,TSU,TRUC,NCOURB) 79 : C-------------------------------------------------------------------- 80 : C Estimation du pas de temps par integration de la loi 81 : C d'evolution du dommage a contr. cte 82 : C--------------------------------------------------------------------- 83 : EA=XMAT(9) 84 : ER=XMAT(8) 85 : TAUX = 1.0D10 86 : IF(VAR(3).GE.0.3.AND.X.GT.0.0D0) THEN 87 : AK1 = AK + 1 88 : TAUX = (1.0D0-VAR(3))**AK1/AK1/(X/EA)**ER 89 : ENDIF 90 : C 91 : RETURN 92 : END 93 : 94 :
© Cast3M 2003 - All rights reserved.
Disclaimer