UVSCPL

Download as pdf or txt
Download as pdf or txt
You are on page 1of 3

CHAPTER 4 337 Viscoplasticity and Generalized Plasticity User Subroutines

UVSCPL
Definition of the Inelastic Strain Rate
Description This user subroutine is used for computing the inelastic strain increment for an elastic-viscoplastic material. This routine allows very general material laws to be entered. The user must define the inelastic strain and the stress increment. This user subroutine is activated when the implicit creep procedure is used, and VISCO PLAS material is selected on the ISOTROPIC or ORTHOTROPIC option. Format User subroutine UVSCPL is written with the following headers:
SUBROUTINE UVSCPL(YOUNG,POISS,SHEAR,B,USTRRT,ETOT,E, 1 THMSTI,EELAS,S,SINC,GF,EPL,AVGINE,EQCRP,EQCPNC,YD,YD1, 2 VSCPAR,DT,DTDL,CPTIM,TIMINC,XINTP,NGENS,M,NN,KCUS,MATUS, 3 NDI,NSHEAR,NCRD,IANISO,NSTATS,INC,NCYCLE,LOVL,NVSPLM) IMPLICIT REAL *8 (A-H, O-Z) DIMENSION POISS(3,2),YOUNG(3,2),B(NGENS,NGENS), 1 USTRRT(NGENS),ETOT(NGENS),E(NGENS),THMSTI(NGENS), 2 EELAS(NGENS),S(NGENS),SINC(NGENS),GF(NGENS),EPL(NGENS), 3 AVGINE(NGENS),DT(NSTATS),DTDL(NSTATS),XINTP(NCRD), SHEAR(3,2),VSCPAR(NVSPLM),MATUS(2),KCUS(2) user coding RETURN END

where:
Input: YOUNG POISS SHEAR B ETOT E is the Youngs modulus. is the Poissons modulus. is the shear modulus. is the tangent elastic matrix. is the accumulated total strain at beginning of increment. is the current strain increment.

338 Marc Volume D: User Subroutines and Special Routines

THMSTI EELAS S EPL EQRCP EQCPNC YD YDL VSCPAR DT DTDL CPTIM TIMINC XINTP NGENS M NN KCUS(1) KCUS(2) MATUS(1) MATUS(2) NDI NSHEAR NCRD IANISO NSTATS INC NCYCLE LOVL NVSPLM

is the thermal strain increment. is the accumulated elastic strain at beginning of increment. is the accumulated stress at beginning of increment. is the accumulated inelastic strain at beginning of increment. is the equivalent inelastic strain at beginning of increment. is the increment equivalent inelastic strain. is the flow stress at temperature t. is the flow stress at temperature t + dt. is the viscoplastic data read off isotropic or orthotropic option. is the state variables at beginning of increment. is the incremental state variables. is the elapsed time at beginning of increment. is the time increment. is the integration point coordinates. is the number of strain components. is the element number. is the integration point number. is your layer number (always 1 for continuum elements). is the internal layer number (always 1 for continuum element). is the user material identifier. is the internal material identifier. is the number of direct components. is the number of shear components. is the number of coordinate directions. is the flag to indicate nonisotropic elasticity. is the number of state variables. is the increment number. is the cycle number. = 4 during stiffness formation. = 6 during residual calculation. is the number of viscoplastic data read from input.

Required Output: USTRRT SINC is the inelastic strain rate. is the stress increment.

CHAPTER 4 339 Viscoplasticity and Generalized Plasticity User Subroutines

GF AVGINE

is the change in stress due to change in elastic material properties associated with DT. is the inelastic strain increment.

Note: To ensure convergence, it should be noted that the returned values of these quantities must be mutually compatible; that is, they simultaneously must satisfy within tolerance: 1. SINC = B*(E - AVGINE - THMSTI) + GF 2. the creep law employed.

The tolerance should be at least one order of magnitude smaller than the global Newton-Raphson tolerance. The values of USTRRT, AVGINE, and SINC are expected to be returned from the routine for both LOVL=4 and LOVL=6.

You might also like