Rocscience International Conference 2025 is going to take place in Sydney, Australia Read more

Search Results

Dynamic Data Analysis

1.0 Introduction

This tutorial starts from introducing the newest features in RS2 Dynamic Data Analysis. In this tutorial, an embankment is subjected to a dynamic loading at the base. The model consists of two soil layers and four stages and groundwater analysis is involved. In the starting model file, project settings, geometry and field stress have been defined. The tutorial will proceed from defining and assigning material properties. The dynamic data filter feature is incorporated in the model simulation. The focus of this tutorial is to compare the pore pressure build-up under seismic waves using two different material models: Finn and Mohr-Coulomb.

  1. Select Files > Recent Folders > Tutorials Folder
  2. Open the Load Input Filter tutorial file in the Dynamic Data Analysis tutorials folder.

2.0 Dynamic Data Analysis

  1. Select: Dynamic > Dynamic Data Analysis Dynamic Data Analysis icon

Once the load function time histories are entered, amplitude spectrum, time-acceleration, time-velocity and time-displacement graphs are plotted based on the input data. It can be selected that either the amplitude spectrum is generated based on acceleration or velocity.

The load function time history is already inputted in this tutorial, it can be found from the first sheet of the Excel file titled “Mexico City 1985 - Mesa Bivradora C. U.” at: C:\Users\Public\Documents\Rocscience\RS2 Examples\tutorials\Dynamic Analysis\Dynamic Slope Analysis

Dynamic Data Analysis dialog box

Response spectrum and arias intensity are also automatically plotted.

Response Spectrum plot

Arias Intensity plot

2.1 Response Spectrum

  1. Select: Response Spectrum from the left graph.

The maximum period and damping ratio can be adjusted to compare with the original response spectrum.

Response Spectrum graph

2.2 Filter Spectrum

  1. Select: Filter Spectrum, second tab from the data.

The minimum and maximum frequencies can be adjusted for filtering. As they are adjusted, all graphs are also replotted based on the filtered data set. Both the original and filtered dataset are plotted on all graphs for easy comparison. To find the global minimum and maximum frequencies of the entire data set, select Reset button to reset the Min and Max.

It can be observed from the Amplitude Spectrum that the amplitude is very minimal for frequencies bigger than 10, therefore, we set maximum frequency to 10 for the filter spectrum.

Dynamic Data Analysis dialog

Response Spectrum graph

Arias Intensity

2.3 Baseline Correction

  1. Select: Baseline Correction, third tab from the data.

Two methods are available for baseline correction: Cut Off and Sinusoid. Both methods can be applied to velocity to correct the velocity drift. However, only Sinusoid method can be applied to displacement to correct the displacement drift. Once either velocity or displacement is corrected, new sets of acceleration and displacement or acceleration and velocity are calculated based on the corrected values and plotted on the graphs to compare with the original dataset. Please note that the baseline correction is applied to the original data, not the filtered data from Filter Spectrum.

For Cut Off method, an Applicable Cut Off Range is generated based on velocity data set, giving the users a reference on the range of Cut off value that outputs the optimal baseline correction. Click Apply to process the Cut Off. The Undo button is also available for users to take a step back.

Dynamic Function 1 dialog

Amplitude Spec dialog

Similarly, a default value for the Sinusoid method, which is taken as the displacement/velocity value at the final timestep, is generated based on either velocity or displacement data set as a reference for the users. After clicking Apply, the residual value changed to 0, which means the baseline has been corrected.

Dynamic Function 1 dialog

Amplitude Sepc dialog

Dynamic Function 1 dialog

Amplitude Spec dialog

3.0 Defining Material Properties

  1. Select: File > Recent Folders > Tutorial Folder
  2. In the Dynamic Analysis > Dynamic Data Analysis folder, open MC and Rayleigh (Initial).fez

This file has the project setting, model geometry, boundaries, field stress, mesh, and restraints created. This tutorial will proceed to define material properties, determine groundwater conditions, set up dynamic boundary conditions, filter and apply the dynamic loading. Natural frequencies will then be computed to obtain the Rayleigh Damping parameters and apply to the project settings.

Image of model

3.1 Material Properties

Material & Staging workflow tab

The only difference between the models using Finn model and Mohr-Coulomb model is the material properties, everything else is identical. Two material layers are used in this tutorial. This tutorial will start from the Mohr-Coulomb model.

  1. Select: Properties > Define Material
  2. Enter material properties for the Mohr-Coulomb model as follows.

MATERIAL 1

Initial Conditions

Define Material Properties- Initial Conditions dialog box

Strength Parameters

Define Material Properties- Strength Parameters dialog

Stiffness

Define Material Properties- Stiffness dialog

Stage Factors

Define Material Properties- Stage Factors dialog

Hydraulic Properties

Define Material Properties- Hydraulic Properties dialog

EMBANKMENT

For the embankment material, only initial conditions, stiffness and strength are different as shown in the following screenshots, the hydraulic properties, datum dependency and stage factors remain the same as material 1.

Initial Conditions

For Initial Conditions, the Initial Element Loading = Body Force Only:

Define Material Properties- Initial Conditions dialog

Strength

For Strength, Residual Cohesion (MPa) = 0.1:

Define Material Properties- Strength dialog

Stiffness

For Stiffness, the Young’s Modulus (MPa) = 500:

Define Material Properties- Stiffness dialog

Stage Factors

Stage Factors for Embankment is the same as Material 1:

Define Material Properties- Stage Factors dialog

Hydraulic Properties

Hydraulic Properties for Embankment is the same as Material 1:

Define Properties Properties

Select OK to save and close the Define Material Properties dialog.

3.2 Assign Material Properties

Click through all stages, ensure in Stage 1, Material 1 is assigned to the bottom layer while the top layer is excavated. From Stage 2 to 6, Material 1 is assigned to the bottom layer, and Embankment to the top layer.

4.0 Assigning Groundwater Conditions

Make sure the Groundwater workflow tab is selected throughout this section.

4.1 Set Boundary Conditions

  1. Select: Groundwater > Set Boundary Conditions
  2. Set BC Type to Total Head and input a Total head value of 50m.
    Set Boundary Conditions dialog box
  3. Make sure Stage 1 is selected.
  4. Assign this boundary condition to the left and right sides at stage 1 as the follows:
    Image of model
  5. Click Close to exit the Set Boundary Condition dialog.
  6. Select Stage 3 (Pressure), right click at the left side. Select Remove Boundary Condition as the follows:
    Image of model
  7. Repeat the same to the right side. Now the 50 m total head boundary condition should be applied only at stage 1 and 2 and removed at stage 3 to 6.
  8. Select: Groundwater > Set Boundary Conditions Set Boundary Conditions icon
  9. Set BC Type to Total Head and input a Total head value of 100m.
    Set Boundary Conditions dialog box
  10. Make sure Stage 3 (Pressure) is selected. Assign this boundary condition to the left boundary of embankment, and top left boundary of Material 1 as the follows:
    Image of model
    Click through stages, now the 100 m total head should be applied to stages 3 to 6.
  11. In the Set Boundary Conditions dialog, set BC Type to Unknown (P=0 or Q=0) and select Apply at this stage only (do not propagate).
    Set Boundary Conditions dialog box
  12. Make sure Stage 3 (Pressure) is selected. Assign this boundary condition to the top and right boundary of the embankment, and the right top boundary of Material 1 as the follows:
    Image of model
    Click through all stages, now the unknown boundary condition should apply only to stage 3.
  13. In the Set Boundary Conditions dialog, set BC Type to Zero Pressure (P=0) and Selection mode to Boundary Nodes as the follows:
    Set Boundary Conditions dialog
  14. Make sure stage 4 (10 s) is selected. Click and hold the left mouse button, and drag a selection window which encloses a portion of top right boundary at the ground surface. Hit Enter or click Apply to assign the Zero Pressure boundary condition to the selected boundary nodes as the follows:
    Image of model
    Click through stages, now the Zero Pressure boundary condition should be applied to stage 4 – 6.
  15. Click Close to exit the Set Boundary Conditions dialog.

4.2 Add Ponded Water Loads

Make sure Stage 1 is selected.

  1. Select: Loading > Ponded Water Loads > Add Ponded Water Loads
  2. For Ponded Water Load 1, select the Use Calculated Water Pressure option as shown below:
    Add Ponded Water Load dialog box
  3. Click OK.
  4. Enter coordinate (-50, 50) as the first point, and (75, 50) as the second point. Hit Enter to complete selection. Ponded Water Load 1 will be added to the boundary as the follows:
    Image of model
  5. Select Stage 2.
  6. Select: Loading > Ponded Water Loads > Add Ponded Water Loads
  7. For Ponded Water Load 2, select the Use Calculated Water Pressure option. Check the Stage Load checkbox as shown below:
    Add Ponded Water Load dialog box
  8. Click the Stage Application button. Check the Apply option for stages 2 to 6 as the follows:
    Stage Factors dialog box
  9. Click OK to close both the Stage Factors dialog and the Add Ponded Water Load dialog.
  10. Enter coordinate (75, 50) as the first point, and (140, 100) as the second point. Hit Enter to complete selection. Ponded Water Load 2 will be added to the sloped boundary as follows:

Image of model

Click through all stages. Ponded Water Load 1 should be applied to all stages, and Ponded Water Load 2 to Stages 2 - 6.

5.0 Assigning Dynamic Conditions

Dynamic workflow tab

5.1 Dynamic Data Analysis

In this tutorial, we use the seismic acceleration data from the 1985 Mexico City Earthquake. The data can be found in the Excel Sheet in the installation folder.

C:\Users\Public\Documents\Rocscience\RS2 Examples\tutorials\Dynamic Analysis\Dynamic Slope Analysis

  1. Open the Excel file titled “Mexico City 1985 - Mesa Bivradora C. U.
  2. Ensure the first sheet, titled “Mexico City 1985” is selected.
  3. Select: Dynamic > Dynamic Data Analysis
  4. For Dynamic Function 1, under Data Input tab, select Input = Acceleration.
  5. Copy and paste the time and acceleration data from Excel Sheet, sheet 1 “Mexico City 1985” (Column A and B) to the data grid as shown below:
    Dynamic Data Analysis dialog
  6. In the middle section of the dialog, make sure Input Type = Acceleration for the Amplitude Spec tab.
  7. Select the Filter Spectrum tab.
  8. As observed from the Amplitude Spectrum, little amplitude occurs beyond 10 Hz frequency. Thus, we can set Max = 10 and click Apply. Your dialog should look as shown below:
    Dynamic Data Analysis dialog
  9. Select the Baseline Correction tab.
  10. Choose Process data = Displacement, and the Method will be set to Sinusoid. By default, Residual Value = 0.2914 m, and Manual Correction is set to the same.
    Dynamic Data Analysis dialog
  11. Click Apply to apply a Manual Correction of 0.2914 m. The dialog will display as shown:
    Dynamic Data Analysis dialog
    Now the Residual Value has changed to 0 m. The change is not significant on Acceleration (m/s^2) vs. Time (s) or Velocity (m/s) vs. Time plots, but it is noticeable on Displacement (m) vs. Time plot.
  12. Select the Filtered Data tab.

Output data type can be selected from the dropdown bar. We will use Velocity output as the dynamic load input for the next step.

Dynamic Data Analysis dialog

5.2 Define Dynamic Loads

  1. Select: Dynamic > Define Dynamic Loads
  2. Rename “Dynamic Load 1” to “Mexico City Velocity Data”. Select Type = Velocity. Check the Compliant Base checkbox. Your dialog should look as follows:
    Define Dynamic Loads dialog
  3. Click the Define button.
  4. In the Velocity vs. Time dialog, check the “Use filtered data” option at the bottom, and select Dynamic Function 1. The filtered velocity data earlier in Dynamic Data Analysis should be inputted as shown below:
    Velocity vs Time inputs
  5. Click OK in both the Velocity vs. Time and Define Dynamic Loads dialogs.

5.3 Set Dynamic Boundary Conditions

  1. Go to tab Stage 4 (the dynamic stage), and select Dynamic > Set Dynamic Boundary Conditions
  2. Select BC Type > Absorb for the bottom boundary and select Transmit for the boundaries on the sides:
    Image of model
  3. Click Close to exit the Set Dynamic BCs dialog. Now the dynamic boundary conditions should be applied to stage 4 - 6.

5.4 Assign Dynamic Load

  1. Select: Dynamic > Add Dynamic Load…
  2. Select Load Function = Mexico City Velocity Data (Velocity) and set Installed at stage to 4 as follows:
    Add Dynamic Load dialog
  3. Click OK to close the dialog.
  4. Use the selection tool to select the bottom boundary. Press Enter to complete selection. The Dynamic Load will apply as follows:

Image of model

5.5 Define Rayleigh Damping

To decide the Rayleigh Damping parameters, we need to compute natural frequency based on the filtered velocity.

  1. Select: Dynamic > Compute Natural Frequencies…
  2. Enter Minimum Frequency = 1e-05 Hz and Maximum Frequency = 20 Hz. Then click Compute…
    Natural Frequencies dialog box
  3. Once the computation is done, select save and view results. The Natural Frequencies Results – Rayleigh Damping dialog should look as follows:

Natural Frequencies Results- Rayleigh Damping dialog

The fourth and fifth frequency numbers on the table (0.535 Hz and 1.054 Hz) need to be inputted in project settings in the next step.

  1. Select: Analysis > Project Settings and go to the Dynamic tab.
  2. Under Rayleigh Damping, select the second option.
  3. Enter:
    • Frequency 1 = 1.054, Frequency 2 = 0.535
    • Damping Ratio 1 = 0.05, and Damping Ratio 2 = 0.05.

    Now the defined dynamic analysis parameters are inputted. Your dialog should look as follows:

    Project Settings dialog box

  4. Click OK.

5.6 Add Time Queries

  1. Select: Dynamic > Time Query > Add Solid Time History Query
  2. Enter the coordinates (113, 50). Hit Enter to add the point and hit Enter again to complete selection. Now a time query point on the material boundary between Embankment and Material 1 has been added.

This step concludes the completion of the Mohr-Coulomb model construction.

  1. Save the model as “MC and Rayleigh.fez”.

6.0 Change to Finn Model

To compare the results between Mohr-Coulomb and Finn model, save the model we just created as a separate file, “Finn and Rayleigh.fez”. In this new saved Finn model, change the material properties as following:

  1. Select: Properties > Define Materials
  2. For Material 1, under the Strength tab, change Failure Criterion to Dynamic > Finn. Under the Dynamic Parameters section, select Finn Formulation = Byrne, Byrne Definition = N1, and (N1)60 = 10.
  3. For Embankment, under the Strength tab, change Failure Criterion to Dynamic > Finn. Under the Dynamic Parameters section, select Finn Formulation = Byrne, Byrne Definition = N1, and (N1)60 = 6.

    The initial conditions, hydraulic properties, datum dependency and stage factors remain the same as the Mohr-Coulomb model. Your dialog should look as follows:
    Define Material Properties dialog
    Define Material Properties dialog
  4. Click OK to close the dialog.

    Other settings remain unchanged. Now, we are ready to run analysis for both models.
  5. Select: Analysis > Compute

7.0 Results

When the models finish running, open Interpret to see the results.

  1. Select: Analysis > Interpret
  2. Click the tab for Stage 6 to observe the final pore pressure
  3. Change the contours to plot pore pressure using the drop-down menu in the toolbar.

7.1 Pore Pressure Contour

MOHR-COULOMB MODEL

Mohr-Coulomb Model

FINN MODEL

Finn Model

7.2 Pore Pressure Graphs at Time Query

Right click the time query and select pore pressure in the menu to plot the pore pressure changes through all stages over time for both models.

MOHR-COULOMB MODEL

TIme query pore pressure graph

FINN MODEL

The query pore pressure graph

Notice that Finn model results in higher pore pressure at the end of the time period than the Mohr-Coulomb model.

This concludes the tutorial.

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 PDF File Calendar Location Language 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 Settle3 Slide2 Slide3 SWedge UnWedge RocTunnel3 RocSlope2 BlastMetrix ShapeMetriX Fragmenter Commercial License Education License Trial License Shop safe & secure Money-back guarantee