U.S. Department of Transportation
Federal Highway Administration
1200 New Jersey Avenue, SE
Washington, DC 20590
2023664000
Federal Highway Administration Research and Technology
Coordinating, Developing, and Delivering Highway Transportation Innovations
This report is an archived publication and may contain dated technical, contact, and link information 

Publication Number: FHWAHRT04095
Date: November 2004 

Manual for LSDYNA Soil Material Model 147PDF Version (1.96 MB)
PDF files can be viewed with the Acrobat® Reader® CHAPTER 2. USER'S MANUALThe user's manual was written as the model was being implemented, verified, and validated. The user's manual consists of a user input guide (much like material model sections in the LSDYNA user's manual); a brief theory manual (LSTC theory manual), which is a condensed version of the first section of this report; and a discussion of the use of the model. Both manuals will be added to an updated LSDYNA manual. The user's manual addresses the basics of the model, input parameters, and basic equations. Table 2 contains a brief description of the user input variables for the soil model, along with the corresponding symbols used in the LSTC theory manual. The bold text is the LSTC theory manual symbol, which is typically followed by a brief description and then the user input value symbol. The parameters that need to be specified are dependent on the soil and the specific application. Table 2. Input parameters for soil model.
USER INPUT GUIDEFHWA Soil Material Model Input *MAT_FHWA_SOIL_OPTION
This is material type 147. This is an isotropic material with damage and is available for solid elements in LSDYNA. The model has a modified MohrCoulomb surface to determine the pressuredependent peak shear strength. It was developed for applications involving roadbase soils. *MAT_FHWA_SOIL_NEBRASKA It is an option to use the default properties determined for soils used at the University of Nebraska at Lincoln. The default units used for this material are millimeter (mm), millisecond (ms), and kilogram (kg). If different units are desired, the conversion factors must be input. Card Format
*MAT_NCHRP_SOIL_blank Define the following cards: Card Format
THEORY MANUALMAT_FHWA_SOIL A brief discussion of the FHWA soil model is given. The elastic properties of the soil are isotropic. The implementation of the modified MohrCoulomb plasticity surface is based on the work of Abbo and Sloan. ^{(9)} The model is extended to include excess porewater effects, strain softening, kinematic hardening, strainrate effects, and element deletion. The modified yield surface is a hyperbola fitted to the MohrCoulomb surface. At the crossing of the pressure axis (zero shear strength), the modified surface is a smooth surface and it is perpendicular to the pressure axis. The yield surface is given as
where: P=pressure = internal friction angle K ()= function of the angle in deviatoric plane = square root of the second invariant of the stress deviator c = amount of cohesion, J_{3} = third invariant of the stress deviator ahyp = parameter for determining how close to the standard MohrCoulomb yield surface the modified surface is fitted If ahyp is input as zero, the standard MohrCoulomb surface is recovered. The input parameter ahyp should be set close to zero, based on numerical considerations, but always less than . It is best not to set the cohesion, c, cot to very small values since this causes excessive iterations in the plasticity routines. To generalize the shape in the deviatoric plane, the standard MohrCoulomb K () function was changed to a function used by Klisinski:^{(10)}
where: J_{3} = third invariant of the stress deviator e = material parameter describing the ratio of triaxial extension strength to triaxial compression strength If e is set to 1, then a circular cone surface is formed. If e is set to 0.55, then a triangular surface is formed. K () is defined for 0.5. < e 1.0. To simulate nonlinear strain hardening behavior, the friction angle is increased as a function of the effective plastic strain:
where: _{eff plas} = effective plastic strain A_{n} = fraction of the peak strength internal friction angle where nonlinear behavior begins, 0 < A_{n}  1 The input parameter E_{1} determines the rate of the nonlinear hardening. If there is no strain hardening, then _{max} = _{init= } . To simulate the effects of moisture and air voids, including excess porewater pressure, both the elastic and plastic behaviors can be modified. The bulk modulus is:
where: K_{i} = nonporous bulk modulus n _{cur} = current porosity = w = volumetric strain corresponding to the volumne of the air voids = n (1S) _{v} = total volumetric strain D _{1} = material constant controlling the stiffness before the air voids are collapsed n = porosity of the soil = e = void ratio = S = degree of saturation = , _{sp}, m _{c}, P _{w} = soil density, specific gravity, moisture content, and water density.Figure 24 shows the effect of the D _{1} parameter on the pressurevolumetric strain relationship (bulk modulus). The bulk modulus will always be a monotonically increasing value, that is:
Note that the model is following the standard practice of assuming that compressive stresses and strains are positive. If the input parameter is zero, then the standard linear elastic bulk modulus behavior is used. Figure 24. Pressure versus volumetric strain showing the effects of the D subscript 1 parameter. If D_{1} is not set to zero, the bulk modulus input should be the fully collapsed bulk modulus. To simulate the loss of shear strength caused by excess porewater effects, the model uses a standard soil mechanics technique^{(11)} of reducing the total pressure, P, by the excess porewater pressure, u, to get an "effective pressure," P' :
Figure 25. Effects on pressure caused by porewater pressure. Figure 25 shows how porewater pressure affects the algorithm for the plasticity surface. The excess porewater pressure reduces the total pressure, which lowers the shear strength, . Significant excess porewater pressure can cause the effective pressure to become zero. To calculate the porewater pressure, u, the model uses an equation similar to the equation used for the moisture effects on the bulk modulus:
where: K^{sk} = bulk modulus for soil without air voids (skeletal bulk modulus) n^{cur} = current porosity = Max[0, (we_{v})] w = volumetric strain corresponding to the volume of air voids =n(1 S) ^{v} = total volumetric strain D_{2}= material constant controlling the porewater pressure before the air voids are collapsed D_{2} 0 n = porosity of the soil = e = void ratio = S = degree of saturation = = soil density, specific gravity, and moisture content, respectively The increment porewater pressure is zero if the incremental mean strain is negative (tensile). Figure 26 is a plot of the pore pressure versus volumetric strain for different parameter values. With the D_{2} parameter set relatively high compared to K_{sk}, there is no pore pressure until the volumetric strain is greater than the strains associated with the air voids. However, as D_{2} is lowered, the pore pressure starts to increase before the air voids are totally collapsed. The K_{sk} parameter affects the slope of the postvoid collapse pressurevolumetric strain behavior The parameter D_{2} is found from Skempton porewater pressure parameter B, where B is defined as:^{(7)}
Figure 26. Effects of D subscript 2 and K subscript SK parameters on porewater pressure. To simulate strain softening behavior, the FHWA soil model uses a continuum damage algorithm. The strainbased damage algorithm is based on the work of J.W. Ju and J.C. Simo. They proposed a strainbased damage criterionthat is uncoupled from the plasticity algorithm.^{(12,13)} For the damage criterion , , where = pressure and _{pv} = plastic volumetric strain, the damaged stress is found from the undamaged stresses:
where: d = isotropic damage parameter (diso) The damage parameter is found at step j + 1 as:
where = damage threshold surface The meshsensitivity parameter, , is described below. Typically, the damage, d, varies from 0 to a maximum of 1. However, some soils can have a residual strength that is pressuredependent. The residual strength is represented by _{res}, the minimum internal friction angle. The maximum damage allowed is related to the internal friction angle of residual strength by:
If _{res} > 0 , then d _{max} , the maximum damage, will not reach 1 and the soil will have residual strength. When material models include strain softening, special techniques must be used to prevent mesh sensitivity. Mesh sensitivity is the tendency of the finite element model/analysis to produce significantly different results as the element size is reduced. Mesh sensitivity occurs because the softening in the model is concentrated in one element. As the element size is reduced, the failure becomes localized in smaller volumes, which causes less energy to be dissipated by the softening. This can lead to instabilities or, at least, meshsensitive behavior. To eliminate or reduce the effects of strain softening mesh sensitivity, the softening parameter, (the strain at full damage), must be modified as the element size changes. The FHWA soil model uses an input parameter, "void formation," G _{f} , that is like the fracture energy material property for metals. The void formation parameter is the area under the softening region of the pressurevolumetric strain curve times the cube root of the element volume , V ^{1/3} :
with _{0} as the volumetric strain at peak pressure (strain at initial damage (Dint)). Then, can be found as a function of the volume of the element V:
If G_{f} is made very small relative to , then the softening behavior will be brittle. Strainrateenhanced strength is simulated by a twoparameter DevautLions viscoplastic update algorithm developed by Y. Murray.^{(15)} This algorithm interpolates between the elastic trial stress (beyond the plasticity surface) and the inviscid stress. The inviscid stresses () are on the plasticity surface , with and As becomes 1, then the viscoplastic stress becomes the elastic trial stress. Setting the input value _{r} = 0 (gamma) eliminates any strainrateenhanced strength effects. The model allows element deletion if needed. As the strain softening (damage) increases, the effective stiffness of the element can become very small, causing severe element distortion and "hourglassing." The element can be "deleted" to remedy this behavior. There are two input parameters that affect the point of element deletion. Damlev is the damage threshold where element deletion will be considered. Epsmax is the maximum principal strain where the element will be deleted. Both d Damlev and _{pr max} Epsmax are required for element deletion to occur. If Damlev is set to zero, there is no element deletion. Care must be taken when employing element deletion to ensure that the internal forces are very small (element stiffness is zero) or significant errors may be introduced into the analysis. MAT_FHWA_SOIL_NEBRASKA This option gives the soil parameters that were used to validate the material model with experiments performed at the University of Nebraska at Lincoln. The units of these default inputs are milliseconds, kilograms, and millimeters. There are no required input parameters except for material ID (MID). If different units are desired, the appropriate unit conversion factors can be input. DISCUSSION OF SOIL MODEL USEMaterial models for geomaterials (soils, concrete, rock, etc.) tend to be complex. The determination of the input parameters for the models is complicated. In addition, modeling different loading conditions and accurate simulation of boundary conditions add to the complexity involved in using these material models. There are two methods that are typically used to determine the material input variables for soils. The most accurate method is to perform laboratory tests that include both triaxial compression and uniaxial strain tests. These tests can be used to determine the elastic moduli, yield surface parameters, and softening parameters. Typically, these tests use drained soil conditions. Laboratory tests with undrained soil conditions can be used to determine the porewater effects. A second method is to use fullscale testing of the specific application (e.g., a bogie impacting a steel post) to fit the parameters in a trialanderror method. This method requires more time by the analyst. Since the soil model is nonlinear, there may not be a set of unique input parameters that can be determined. Compaction of the soil is typically used to remove some of the air voids that exist in disturbed soils. However, the density, porewater effects, stiffness, and strength are also changed upon compacting the soil. To simulate compaction in highway safety applications where the soil is exposed, we recommend that the values for the soil density, porewater effects, stiffness, and strength be modified. Applying pressure to the ground surface to account for the effects of compaction is a less accurate method that will incorrectly simulate how the soil is deformed at the surface. In fullscale testing or applications, the soil typically extends to infinity. Analyses typically do not extend to infinity, so some type of boundary condition must be applied to the exterior surfaces of a soil analysis model (except for soil surfaces exposed to atmospheric pressure). Standard boundaries reflect dynamic disturbances (stress waves), which does not happen in the real applications. Such reflections can cause serious contamination of the analysis results. Exterior boundaries for analyses involving soil need a nonreflecting boundary. A partial nonreflecting boundary exists in LSDYNA. This boundary is an impedancematching boundary, which is only good for highfrequency (highly transient) behavior. At this time, there is no nonreflecting boundary that matches both low (quasistatic) and high (highly transient) frequency behaviors. Also, only linear behavior is assumed. Thus, to use the current nonreflecting boundary, the material near the boundaries must only behave linearly. Also, the nonreflecting boundary should only experience highfrequency behavior. 