U.S. Department of Transportation
Federal Highway Administration
1200 New Jersey Avenue, SE
Washington, DC 20590
202-366-4000


Skip to content
Facebook iconYouTube iconTwitter iconFlickr iconLinkedInInstagram

Federal Highway Administration Research and Technology
Coordinating, Developing, and Delivering Highway Transportation Innovations

Report
This report is an archived publication and may contain dated technical, contact, and link information
Publication Number: FHWA-HRT-04-094
Date: November 2004

Evaluation of LS-DYNA Soil Material Model 147

PDF Version (2.87 MB)

PDF files can be viewed with the Acrobat® Reader®

CHAPTER 3. BASELINE MODEL: DIRECT SHEAR TEST SIMULATION

In this chapter, a baseline model of the direct shear test is established in order to investigate the effects of the various parameters associated with the soil model.

A finite element model using 10,600 elements was developed based on the large-scale direct shear testing device (see figure 4). This model included solid rigid elements to model the steel, since no deformation had been seen during physical testing, and solid elements to implement the Federal Highway Administration (FHWA) soil model, material type 147.

An overburden pressure of 18.5 kPa was applied and dynamic relaxation was implemented before the soil was sheared in the finite element model. A prescribed motion condition of 1 millimeter per millisecond (mm/ms) was applied to the lower containing cylinder in the device, similar to the quasi-static testing motion condition. Loads and displacements were measured so that simulation results could be directly compared to the physical test results.

View Alternative Text

Figure 4. Finite element model of large-scale direct shear test.

Baseline Parameters

Initial parameters for material type 147 were selected from three LS-DYNA models provided by the developer. These models included two triaxial compression tests (cylin.k and txc3-4pr0c.k) and a hydrostatic tension test (hydten1.k). The initial parameters implemented in the developer's models and the comparison to the baseline model used in the parameter study are shown in table 1. The material input parameters in LS-DYNA keyword format are shown in table 2.

E-mail correspondence from the developer later recommended values for the shear modulus, G, and the bulk modulus, K, as 1.3 megapascals (MPa) and 3.25 MPa, respectively. Because of the lateness of obtaining these values, these changes were not reflected in the material parameter study of chapter 4. However, chapter 5 of this report, "Developer -Recommended Parameters," investigates these recommended values.

Table 1 . Developer's and baseline material parameters.

 

Material Parameter

Input Deck

RO

Nplot

Spgrav

Rhowat

Vn

Gammar

cylin.k

2.350E-6

3

2.79

1.0E-6

1.1

0.0

txc3-4pr0c.k

2.350E-6

3

2.79

1.0E-6

1.1

0.0

hydten1.k

2.350E-6

3

2.79

1.0E-6

1.1

0.0

Baseline

2.350E-6

3

2.79

1.0E-6

1.1

0.0

 

Input Deck

Itermax

K

G

Phimax

Ahyp

Coh

cylin.k

10

0.465

0.186

1.1

1.0E-7

1.0E-6

txc3-4pr0c.k

10

0.465

0.186

1.1

1.0E-7

1.0E-6

hydten1.k

10

0.465

0.186

1.1

1.0E-7

6.2E-6

Baseline

10

0.465

0.186

1.1

1.0E-7

6.2E-6

 

Input Deck

Eccen

An

Et

Mcont

Pwd1

PwKsk

cylin.k

0.7

0.4

10

0.034

0.0

0.0

txc3-4pr0c.k

0.7

0.0

10

0.034

0.0

0.0

hydten1.k

0.7

0.0

0

0.034

0.0

0.0

Baseline

0.7

0.0

0

0.034

0.0

0.0

 

Input Deck

Pwd2

Phires

Dint

Vdfm

Damlev

Epsmax

cylin.k

0.0

1.0E-3

5.0E-5

6.0E-7

0.98

1.00

txc3-4pr0c.k

0.0

0.0E-0

5.0E-5

1.0E-9

0.80

0.03

hydten1.k

0.0

0.0E-0

2.5E-3

5.0E-0

1.00

1.00

Baseline

0.0

0.0E-0

2.5E-3

5.0E-0

1.00

1.00


Table 2 . LS-DYNA format: Baseline values for parameter study, model R3.

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$$$$ FHWA Soil Material

$

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$

$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8

$

$$$$ Material 147 Nebraska soil

$

*MAT_FHWA_SOIL

$ mid ro NPLOT SPGRAV RHOWAT VN GAMMAR ITERMAX

1 2.350E-6 3 2.79 1.0E-6 1.1 0.0 10

$ K G PHIMAX AHYP COH ECCEN AN ET

0.465000 0.186000 1.1 1.0E-7 6.2E-6 0.7 0.0 0.0

$ MCONT PWD1 PWKSK PWD2 PHIRES DINT VDFM DAMLEV

0.034 0.00 0.0 0.0 0.0 0.00250 5.00 1.0

$ Epsmax

1.0

$

Baseline Results

Deformation of the baseline model results (without the frame structure) is shown in figure 5, the force to perform the direct shear simulation is shown in figure 6, and the internal energy absorbed by the soil during the simulation is shown in figure 7. Unfortunately, it appears that the soil model is only useable in this configuration for about 25 mm of deformation (as shown in figure 6). The force builds up nicely and peaks at 718 kilonewtons (kN) at 18.2 mm of displacement. The soil begins to soften, as it should, until it reaches the valley at 24.6 mm of displacement, with a corresponding force of 385 kN. After that time, the force begins to rise, which never occurs in physical tests. The force should either maintain a relative constant value or begin to drop, as previously shown in figures 2 and 3.

One reason for the rise in forces after 25 mm of displacement is because of the high damage level specification. However, if lower values of damage are used, the model becomes unstable and the code terminates. This will be discussed further in chapter 5. For the parameter study in this report (chapter 4), the main results of the comparison will be the peak force and the valley force, and their corresponding displacements and internal energy absorbed by the soil.


View Alternative Text
(a) 0 ms (b) 49.992 ms
Figure 5 . Soil displacement: Baseline model.

 

View Alternative Text

Figure 6. Direct shear force: Baseline model.

 

View Alternative Text

Figure 7. Direct shear soil internal energy: Baseline model.


Previous | Table of Contents | Next

Federal Highway Administration | 1200 New Jersey Avenue, SE | Washington, DC 20590 | 202-366-4000
Turner-Fairbank Highway Research Center | 6300 Georgetown Pike | McLean, VA | 22101