This report is an archived publication and may contain dated technical, contact, and link information 

Publication Number: FHWAHRT04097
Date: August 2007 

Measured Variability Of Southern Yellow Pine  Manual for LSDYNA Wood Material Model 143PDF Version (2.92 MB)
PDF files can be viewed with the Acrobat® Reader® 1.5 PLASTIC FLOWThe plasticity algorithms limit the stress components once the failure criterion in equations 13 or 14 is satisfied. This is done by returning the stress state back to the yield surface.^{1} Our traditional approach for modeling plasticity is to partition the stress and strain tensors into elastic and plastic parts. Partitioning is done with return mapping algorithms that enforce the plastic consistency conditions. Such algorithms allow one to control plastic strain generation. In addition, return mapping algorithms with associated flow satisfy the second law of thermodynamics. Associated flow is discussed in appendix D. ^{1}Two simple types of plasticity algorithms are those that reduce the moduli directly and those that scale back the stresses directly. Although simple to implement, such methods do not allow one to control plastic strain generation and do not necessarily satisfy the second law of thermodynamics. 1.5.1 Consistency Parameter UpdatesSeparate plasticity algorithms are modeled for the parallel and perpendicular modes by enforcing separate consistency conditions. Each consistency condition is derived in appendix D. The solution of each consistency condition determines the consistency parameters Dl_{úú} and Dl_{^}. The Dl solutions, in turn, determine the stress updates. Parallel Modes The parallel failure criterion from equation 13 is redefined as f_{} (I_{1}, I_{4}) ≥ 0, with: where I_{1} and I_{4} are two of the five invariants of a transversely isotropic material. In this case, the expression Dl_{} is: Dl_{} is calculated from specification of the total strain increments and the yield function . The trial elastic stress invariants, and , are calculated from the trial elastic stresses (see section 1.5.2). For the failure surface function in equation 15, the partial derivatives in equation 16 are: Perpendicular Modes The perpendicular failure criterion from equation 14 is redefined as f_{^} (I_{3}, I_{4}) ≥ 0, with: where I_{2} and I_{3} are two of the five invariants of a transversely isotropic material. In this case, the expression for Dl_{^} is: Dl_{^}is calculated from specification of the total strain increments and the yield function , where and are the trial elastic stress invariants. For the failure surface function in equation 21, the partial derivatives in equation 22 are: 1.5.2 Elastoplastic Stress UpdatesThe stresses are readily updated from the total strain increments and the consistency parameters: where: Total strain increments are calculated by the finite element code from the dynamic equations of motion and the time step. Each normal stress update depends on the consistency parameters and failure surface functions for both the parallel (Dl = Dl_{} and f = f_{}) and perpendicular (Dl = Dl_{^} and f = f_{^}) modes. Each shear stress update depends on just one consistency parameter and failure surface function. If neither parallel nor perpendicular yielding occurs and , then Dl = 0 and the stress update is trivial: . For the stress updates in equation 27, the partial derivatives are: 
Topics: research, safety, infrastructure, structures, materials Keywords: research, safety, infrastructure, structures, materials, wood, LSDYNA, orthotropic, material model, damage, rate effects, guardrail TRT Terms: guardrails, wood Updated: 04/12/2012
