sste
C SSTE SOURCE JB251061 23/05/10 21:15:17 11667 SUBROUTINE SSTE *--------------------------------------------------------------------- * ECOULEMENT INELASTIQUE *--------------------------------------------------------------------- * * SYNTAXE POUR LES LOIS PLASTIQUES: * * SIGF VARF DEPIN RI1 = 'SUBSTEP' MODL SIG0 VAR0 DEPST CARAC * (PRECIS) (NMAXSSTEPS) (NITMAX) * (NNUMER) (DELTAX); * IN= * MMODEL | MODL OBJET MODELE * MCHAML | SIG0 CONTRAINTES AU DEBUT DU PAS * | VAR0 VARIABLES INTERNES AU DEBUT DU PAS * | DEPST INCREMENT DE DEFORMATIONS TOTALES * | CARAC CONSTANTES DU MATERIAU * FLOTTANT | PRECIS PRECISION POUR ITERATIONS INTERNES * INTEGER | NMAXSSTEPS Maximum number of substeps * INTEGER | NITMAX Maximum number of iterations * INTEGER | NNUMER Numerical differentiation type (MRSLade) * FLOTTANT | DELTAX Numerical differentiation parameter (MRSLade) * * OUT= * MCHAML | SIGF CONTRAINTES A LA FIN DU PAS * | VARF VARIABLES INTERNES A LA FIN DU PAS * | DEPIN INCREMENT DE DEFORMATIONS INELASTIQUES * | RI1 CONSISTENT TANGENT MATRIX * *--------------------------------------------------------------------- IMPLICIT INTEGER(I-N) IMPLICIT REAL*8(A-H,O-Z) -INC PPARAM -INC CCOPTIO -INC SMCHAML * * Lectura de datos * * modelo IF (IERR.NE.0) RETURN * tensiones iniciales IF (IERR.NE.0) RETURN IF(IERR .NE. 0) RETURN MCHELM=IPCHE1 SEGACT MCHELM c! Usa L1, N1, N3 IF(MCHELM.TITCHE.NE.'CONTRAINTES')THEN MOTERR(1:32)='CONTRAINTES' RETURN ENDIF * variables internas iniciales IF (IERR.NE.0) RETURN IF(IERR .NE. 0) RETURN MCHELM=IPCHE2 SEGACT MCHELM IF(TITCHE.NE.'VARIABLES INTERNES')THEN MOTERR(1:32)='VARIABLES INTERNES' RETURN ENDIF * incremento de deformacion IF (IERR.NE.0) RETURN IF(IERR .NE. 0) RETURN MCHELM=IPCHE4 SEGACT MCHELM IF(TITCHE.NE.'DEFORMATIONS')THEN MOTERR(1:32)='DEFORMATIONS' RETURN ENDIF * caracteristicas materiales y geometricas IF (IERR.NE.0) RETURN IF(IERR .NE. 0) RETURN MCHELM=IPCAR SEGACT MCHELM IF(TITCHE.NE.'CARACTERISTIQUES')THEN MOTERR(1:32)='CARACTERISTIQUES' RETURN ENDIF * precision para iteraciones internas IF(IRT.EQ.0) PRECIS=1.D-10 * num max de substeps IF (IRT.EQ.0) NMAXSSTEPS=1 * num max de iteraciones internas IF (IRT.EQ.0) NITMAX=25 * tipo de derivacion numerica IF (IRT.EQ.0) NNUMER=0 * delta de x relativo para las derivadas numericas IF (IRT.EQ.0) DELTAX=0.D0 * comprovar que es def plana o axisimetrico IF ((IFOUR.gt.0).or.(IFOUR.lt.-1)) then write(*,*) ' Formulacion no disponible' return endif . PRECIS,NITMAX,NMAXSSTEPS,NNUMER,DELTAX, . IPCHE7,IPCHE8,IPCHE9,IPRIGI) IF(IERR.EQ.0)THEN ENDIF END
© Cast3M 2003 - Tous droits réservés.
Mentions légales