Search Results

4 - Nonlinear Analysis with Pore Water Pressure

1.0 Introduction

This tutorial is composed of a nonlinear analysis as a quasi-effective stress analysis, with generation and dissipation of pore water pressure.

In this analysis, we demonstrate the porewater generation option in RSSeismic. The user has the option to use one of several models. The porewater-pressure formulation used represents a quasi-coupling commonly used in 1-D nonlinear site response analysis platform, whereby the generated porewater pressure results in a degradation of soil stiffness. The user has the option to allow for porewater pressure dissipation during ground shaking as can control the model hydraulic boundary condition (permeable or impermeable)

Topics Covered in this Tutorial:

  • Nonlinear analysis time domain analysis
  • Pore water pressure generation and dissipation

Finished Product:

The finished product of this tutorial can be found in the Tutorial 4 folder. All tutorial files installed with RSSeismic can be accessed by selecting File > Recent Folders > Tutorials Folder from the RSSeismic main menu.

2.0 Project Setup

  1. Begin by opening the RSSeismic program and selecting New Project new project icon. A new project will open, and you will be taken to the Project Setup tab.
  2. In the Project Setup, set the Unit System to Imperial, stress as psf.
  3. For the Profile Generation, we will be leaving the default: Create Profiles Manually. Since we are only defining one soil profile, Create Multiple Profiles will remain turned off.
  4. For Analysis Method select Nonlinear.
    1. Notice that the Pore Pressure Options are no longer greyed out. Tick the checkbox for Generate excess pore water pressure.
    2. Select Enable Dissipation. Leave the permeability options top and bottom of the profile boundaries selected.
  5. For the Solution Type, Time Domain is selected by default for nonlinear analysis.
  6. For Default Soil Model, ensure General quadratic/hyperbolic model (GQ/H) is selected.
  7. For Default hysteretic re/unloading formulation ensure the default Non-masing re/unloading (recommended) is selected.
When conducting a Nonlinear analysis, users have the option of obtaining the site response results using the Equivalent Linear method automatically. This can be done by using the Complementary analyses option. We will not be using this option for the tutorial.
  1. For Complementary analyses, turn OFF Equivalent Linear – Frequency Domain.

project setup

3.0 Profile

  1. Go the Profile tab profile tab icon

For this tutorial, we will be creating a soil profile that consist of 5 soil layers.

  1. Layer 1 is created by default. Add 4 more layers to the soil profile. This can be done by clicking Add Layer Below 4 times, or by clicking Append Rows and entering 4. You should now how 5 layers total
  2. Change the Thickness of each layer as indicated in the table below:

Name

Thickness (ft)

Layer 1

10

Layer 2

15

Layer 3

15

Layer 4

20

Layer 5

20

  1. Tick the Water table at top of layer checkbox to add a water table above Layer 1.

profile tab

4.0 Properties

  1. Go to the Properties tab properties icon
  2. Select the Advanced Table View.
  3. Copy and paste the basic soil profile values from the table below.

Name

Thickness (ft)

Unit Weight (pcf)

Shear Wave Velocity (ft/s)

Shear Strength (psf)

Layer 1

10

125

1000

3288

Layer2

15

125

1500

7624

Layer 3

15

125

1500

8166

Layer 4

20

125

2000

14237

Layer 5

20

125

2000

14960

We will use the same curve fitting procedure we used in the previous tutorial.

  1. Under Reference Curve > Sand select the Darendeli, 2001 reference curve.
  2. Enter Ko = 0.5. Leave all others as the default values shown below:
    1. OCR = 1
    2. Ko = 0.5
    3. N = 10
    4. Frequency = 1
    5. PI = 0
  3. Under Curve Fitting, enter Fitting Procedure = MRDF with UIUC Reduction Factor.
  4. Click Fit.
  5. The Fitting Limits dialog will appear. Leave the default values (Max strain = 0.05% and Min Strength = 95%) selected and click OK.
  6. The Soil Model properties and Reduction factor formulation will be calculated. Click Use Fit to apply the fitted curve to Layer 1.

layer 1 properties with fit applied

  1. Repeat steps 4-9 for Layers 2-5.
  2. Return to Layer 1 and select the PWP Properties tab.
  3. For the PWP Model select Sand-Vucetic Dobry.
  4. Enter the parameters below:

Max ru

Cv (ft2/sec)

Cv–exponent

f

p

F

s

y_tvp

ν

0.95

0.1

0.0

1.0

1.0

0.73

1.0

0.02

3.8

layer1 pwp properties

  1. Enter the same Pore Water Pressure Properties for Layers 2-5. (This can be done more quickly by copy and pasting the values into the Advanced Table View tab. A CSV of the Advanced Table View values has been included in your tutorial folder to make this easier.)

advanced table view of pwp properties

  1. Return to the Layer Properties tab and select the Bedrock layer. Enter:
    1. Halfspace option = Elastic Halfspace
    2. Shear wave velocity (ft/s) = 5000
    3. Unit Weight (pcf) = 160
    4. Damping ratio (%) = 2
    5. For Halfspace porewater pressure, tick the checkbox for Use Cv of last layer.

bedrock layer properties

  1. Select the Profile Plots tab to review your data.

profile plots

5.0 Motions

  1. Select the Motions tab motions tab icon
  2. Under Resources\Input Motions select the Kobe motion.

The time-histories (Acceleration, Velocity, Displacement, Arias and Housner Intensity), FAS, and 5% damped spectral acceleration for Kobe motion can be viewed below.

motions tab kobe

6.0 Damping Formulation

  1. Select the Damping Formulation tab damping formulation icon
  2. For Damping Matrix Type we will use the default selection of Frequency Independent.
  3. For Damping matrix update, use the default Do not update matrix.

damping formulation tab

7.0 Compute Options

  1. Go to the Compute Options tab compute icon
  2. For the nonlinear Time Domain analysis use the default values:
    1. Step Control = Flexible
    2. Maximum Strain Increment (%) = 0.005 %
    3. Integration scheme = Implicit: Newmark Beta Method
    4. Time History Interpolation Method = Linear in time domain
  3. For Output settings, select Layers = All Layers, since we will be comparing the results of the different layers.
  4. Save your project. Select File > Save save icon and name the file Tutorial 4 Nonlinear with PWP.
  5. Click Compute.

compute options

8.0 Results

  1. Go to the Results tab results tab icon
  2. In the Time History Plots tab, tick the checkboxes for Layer 1 and Layer 3 to compare the results.
  3. To better view the layers, use the Color Mode drop down in the bottom right of the screen and select blend motion and layer colour.

The Acceleration, Velocity, Displacement, and Arias Intensity time histories are displayed.

time history plots

  1. Select the Profile Plots tab. The results are shown below.

profile plots results

  1. Select the Stress-strain plots tab. Maximum strain of approximately 0.4% in site profile occurs in Layer 3 due to development of significant level of pore water pressure.

results stress strain plot

  1. Select the Response Spectra Summary tab. 5% damped spectral acceleration, Fourier Amplitude Spectrum and Ratio for Layer 1 and Layer 3 are presented.

response spectra summary

Rocscience logo, click here to return to the homepage Portal Account Portal Account Log In Log Out Home Shopping Cart icon Click here to search our site Click here to close Learning Tech Support Documentation Info Chevron Delete Back to Top View More" Previous Next Edit PDF File Calendar Location Language External Link Apply to ACC External Link Fees Video Click here to visit Rocscience's LinkedIn page Click here to visit Rocscience's YouTube page Click here to visit Rocscience's X page Click here to visit Rocscience's Facebook page Click here to visit Rocscience's Instagram page Click here to visit Rocscience's Reddit page Bookmark Network Scroll down for more Checkmark Download Print Back to top Single User Multiple Users RSLog RocFall3 CPillar Dips EX3 RocFall RocPlane RocSlope3 RocSupport RocTopple RS2 RS3 RSData RSPile RSWall RSSeismic Settle3 Slide2 Slide3 SWedge UnWedge RocTunnel3 RocSlope2 BlastMetrix ShapeMetriX FragMetriX TestLicense Commercial License Education License Trial License Shop safe & secure Money-back guarantee