UVSCPL
UVSCPL
UVSCPL
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.
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.
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.