savcr3
C SAVCR3 SOURCE FD218221 26/02/13 21:15:30 12462 SUBROUTINE SAVCR3(VARF,NVARI,CMAT,NMAT,SIGF, $ SIGEF,NSTRS,IFOUR,EPSFLU,AN,BN,CN,SIGFI) C C_______________________________________________________________________ C C SOUS PROGRAMME APPELE PAR SAVCR QUI CALCUL LES MATRICES DE C DEFORMATION DE FLUAGE DE DESSICCATION A LA FIN DU PAS DE TEMPS C D APRES LE MODELE DE F. BENBOUDJEMA (similaire au modele de Bazant) C_______________________________________________________________________ C IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) DIMENSION VARF(NVARI),CMAT(NMAT),SIGF(NSTRS),SIGEF(NSTRS) DIMENSION EPSFLU(6),SIGFI(NSTRS) DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS) C C RECUPERATION DES VARIABLES INTERNES INITIALES C C FLUAGE PROPRE C EPSFLU(1) = VARF(21) EPSFLU(2) = VARF(22) EPSFLU(3) = VARF(23) EPSFLU(4) = VARF(24) EPSFLU(5) = VARF(25) EPSFLU(6) = VARF(26) C C TEST SI LE TEMPS EST NUL C TDT = CMAT(28) IF (TDT.EQ.0.D0) THEN GOTO 30 ENDIF C C INITIALISATION DES MATRICES C C DO 10 I=1,NSTRS SIGFI(I) = SIGF(I) 10 CONTINUE C C_______________________________________________________________________ C C FLUAGE PROPRE AMORTISSEUR C NDI=NSTRS CALL DEFCR3(VARF,NVARI,CMAT,NMAT,IFOUR,NSTRS, $ AN,BN,CN,NDI) C DO 20 I=1,NSTRS EPSFLU(I) = EPSFLU(I) + AN(I) DO 21 J=1,NSTRS EPSFLU(I) = EPSFLU(I) + BN(I,J)*SIGEF(J) + CN(I,J)*SIGFI(J) 21 CONTINUE 20 CONTINUE C 30 CONTINUE C C SAUVEGARDE DES DEFORMATIONS C C FLUAGE DESSICCATION C VARF(21) = EPSFLU(1) VARF(22) = EPSFLU(2) VARF(23) = EPSFLU(3) VARF(24) = EPSFLU(4) VARF(25) = EPSFLU(5) VARF(26) = EPSFLU(6) C RETURN END
© Cast3M 2003 - Tous droits réservés.
Mentions légales