tufadv
C TUFADV SOURCE CHAT 05/01/13 03:54:58 5004 C TUFADV SOUCE BROC 94/04/09 21:19:18 C---------------------------------------------------------------- & SIGPT,EPSVPT,VARPT,XMAT,XCAR, & NSTRS,NVARI,IFOUR,INPLAS,NCOMAT ) C ----------------------------------------------------------------- IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) C Include contenant quelques constantes dont XPI : -INC CCREEL DIMENSION SIG(*),EPSV(*),VAR(*),SIG1(*),EPSV1(*),VAR1(*) DIMENSION SIGPT(*),EPSVPT(*),VARPT(*),XMAT(*),XCAR(*) DIMENSION DTUFI(2,2),CTUFI(2,2) C---------------------------------------------------------------- C CETTE ROUTINE CALCULE LES VALEURS FINALES DE SIG EPSV VAR A C PARTIR DES VALEURS INITIALES DE L'INCREMENT DE TEMPS ET DES C DERIVEES SIGPT,EPSVPT,VARPT C---------------------------------------------------------------- EPAIS = XCAR(2) RAYOM = XCAR(1)-(EPAIS/2) YOUN = XMAT(1) C EVALUATION DE L INCREMENT D ELONGATION DE FLUAGE DF = DTUFI(1,1) DEF= DTUFI(1,2) DE = DTUFI(2,2) CM = CTUFI(1,1) CMP= CTUFI(1,2) CP = CTUFI(2,2) PHIF = CM*SIG(6) + CMP*SIG(1) ELF = CMP*SIG(6) + CP*SIG(1) RAPF=ELF/PHIF DO 10 I=1,6 SIG1(I)=SIG(I) 10 CONTINUE C________________________________________________________________ C CALCUL DES DEFORMATIONS VISQUEUSES ELONGATION & ROTATION C DUES AU FLUAGE C________________________________________________________________ DO 20 I =2,6 EPSV1(I) =EPSV(I)+TAU*EPSVPT(I) 20 CONTINUE EPSV1(1)=RAPF*EPSV1(6) EPSVPT(1)=RAPF*EPSVPT(6) C________________________________________________________________ C MISE A JOUR DES VARIABLES INTERNES C________________________________________________________________ DO 30 I=1,5 VAR1(I)=VARPT(I) 30 CONTINUE VAR1(6)=VAR(6)+TAU*VARPT(6) VAR1(7)=TAU*VARPT(7) SIG1(1) = SIG(1) + TAU*(SIGPT(1) - DEF*EPSVPT(6)-DE*EPSVPT(1)) SIG1(6) = SIG(6) + TAU*(SIGPT(6) - DF*EPSVPT(6)-DEF*EPSVPT(1)) 99 CONTINUE RETURN END
© Cast3M 2003 - Tous droits réservés.
Mentions légales