![]() ![]() |
![]() ![]() ![]() |
Manual for LS-DYNA Wood Material Model 143APPENDIX E. DERIVATION OF CONSISTENCY PARAMETER FOR PLASTICITY ALGORITHMPlasticity is modeled by enforcing separate consistency conditions for the parallel and perpendicular modes. The goal of each consistency condition is to partition the total strain increments into elastic and plastic components:
The superscripts e and p indicate the elastic and plastic components, respectively. The total strain increments (Deij) are calculated by the finite element code from the dynamic equations of motion and the time step. Once this partition is known, then the stress increments are updated from the elastic strain increments:
Here, n denotes the nth time step in the finite element analysis. E.1 CHECK FOR YIELDINGThe partition into elastic and plastic components requires two steps. The first step is to check for yielding. This is done by temporarily updating the stress components from the incremental strains by assuming that the entire strain increment is elastic:
These updated stresses are called the trial elastic stresses ( E.2 CALCULATE CONSISTENCY PARAMETERThe second step is to enforce the consistency condition if f *> 0. Enforcement of the consistency condition requires an assumption about the direction of plastic flow. It is assumed that the plastic strain increments are normal to the yield surface:
where Dl is a proportionality constant known as the consistency parameter. This assumption is known as an associated flow rule, or normality condition. Use of a potential function other than the yield function in equation 180 results in a nonassociated flow rule. Recent studies reported by Pucik(35) suggest that rate-independent models with nonassociated flow lead to spurious (non-unique) dynamic solutions, so only associated flow is proposed for the present model. The plasticity algorithm calculates Dl by enforcing the plastic consistency condition. This condition is expressed as:
where f is the yield surface function at time increments n to n + 1. The stress state at the beginning of the time step lies on the yield surface, thus, f n = 0. The stress state at the end of the time step is returned to the yield surface by the plasticity algorithm, thus, f n+1 = 0. Therefore, Df = 0. The solution of the consistency condition in equation 181 determines Dl, which, in turn, determines the partitioning of the total strain rate into elastic and plastic components. The stresses are updated from the elastic strain components. Separate Dl solutions are proposed for the parallel and perpendicular modes. Parallel Modes The consistency condition is derived in terms of the invariants rather than the stresses, because the parallel failure criterion in equation 13 is formulated in terms of two transversely isotropic stress invariants. For the purposes of this derivation, the parallel failure criterion from equation 13 is defined as f||(I1,I4) ≥ 0, with:
A first-order Taylor series expansion of the consistency condition
Expansion of the stress invariant increments gives:
where
The expression in the numerator is recognized as the first-order Taylor series expansion of
Perpendicular Modes The perpendicular failure criterion in equation 14 is derived in terms of the invariants rather than the stresses. The perpendicular failure criterion is defined as f^(I2,I3) = 0, with:
Expansion of the consistency condition Df^ = 0 gives:
The stress invariant updates are:
Substitution of these updates into the consistency condition gives the following express for Dl^.
E.3 UPDATE STRESSESThe third step is to update the stresses. There are two options: (1) a purely elastic update and (2) an elastoplastic update. If the trial elastic stress state lies inside the failure surface (f *< 0), then Dl = 0 and the stress state is purely elastic. In this case, the stress update from equations 178 and 179 is trivial:
If the trial elastic stress state lies outside the failure surface (f * > 0), then Dl ¹ 0 and the stress state is elastoplastic. In this case, equations 178, 179, and 180 combine to give:
FHWA-HRT-04-097 |
![]() |