Skip to article

 

Seismic Wave Propagation (SEG-Y Results)

CAE Fidesys allows you to upload solution results in SEG-Y format. This example considers the propagation of seismic waves in the ground based on the Lamb problem for a 2D case. The procedures for setting receivers, saving and subsequent analysis of data in the SEG-Y format are demonstrated.

The model is a part of the plane (xy), a point force is applied to vertex. Non-reflective boundary conditions are applied.

Geometry creation

1. Create a square plate.

On the command bar, select the module for constructing surface geometry (Mode - Geometry, Entity - Surface, Action - Create).

From the list of geometric primitives, select Rectangle.

Set block sizes:

Click Apply.

2. Due to symmetry, we consider half of the model.

On the command bar, select the module for constructing surface geometry (Mode - Geometry, Entity - Surface, Action - Webcut).

From the list of possible kind of webcuts, select Coordinate Plane.

Set the following parameters:

Click Apply.

Do the same, but in the ZX plane.

Click Apply.

As a result, the original Body 1 in the Model Tree will be divided into three bodies (Body 1, Body 2 and Body 3).

3. Delete Surface 3.

On the command bar, select the module for constructing surface geometry (Mode - Geometry, Entity - Surface, Action - Delete).

Set the parameters:

Click Apply.

4. Rotate the model.

On the command bar, select the module for constructing surface geometry (Mode - Geometry, Entity - Surface, Action - Transform).

Select Rotate in the list of operations.  

Set the following parameters:

5. Cut the model.

On the command bar, select the module for constructing surface geometry (Mode - Geometry, Entity - Surface, Action - Webcut).

From the list of possible kind of webcuts, select Coordinate Plane.

Set the following parameters:

Click Apply.

6. Imprint and merge the following commands on the command line:

Meshing

1. Make sure that in the main menu settings (Tools - Options — Mesh Defaults) Default Element Type Hex/Quad is selected.

2. On the command bar, select the surface mesh module (Mode - Mesh, Entity - Surface, Action - Intervals).

Specify the degree of mesh refinement:

Click Apply Size.

Click Mesh.

Specifying the material and block properties

1. In the command bar, select the module for specifying material properties (Mode - Material, Entity - Materials Management).

In the opened widget Material Management in the middle column specify the material name. In the properties column, open the Elasticity list and drag Isotropic Material to the Material Properties column, then in the left column go to the section and drag General → Density and Plasticity → Second Drucker-Prager Criterion. Set the following parameters:

Click Apply.

Close the Materials management window.

2. Create a block of  one material type.

On the command bar, select the block management module (Mode - Blocks, Entity - Block, Action - Add).

Set the following parameters:

Click Apply.

3. Set the block parameters.

On the command bar, select the block management module (Mode - Blocks, Entity - Block, Action - Block properties/parameters).

Set the following parameters:

Click Apply.

Setting boundary conditions

1. Set non-reflective boundary conditions.

On the command panel, select the boundary conditions module (Mode - Boundary Conditions, Entity – Absorbing Condition, Action - Create).

Set the following parameters:

Click Apply.

2. Set the force.

On the command panel, select the boundary conditions module  (Mode - Boundary Conditions, Entity - Force, Action - Create).

Set the following parameters:

Click Apply.

Set the BC dependency on time and / or coordinates

1. Create a formula 1 for strength 1.

On the command panel, select the module for specifying Dependency BC (Mode - BC Dependency).

Click Force 1 and choose Formula panel in the right. Then choose Berlage and set the following parameters:

Click Apply. Close the BC Dependency window.

Receivers

1. Create receivers on curve 16 along all directions.

On the command panel, select module for setting receivers (Mode - Receivers, Operation - Create).

From the drop-down list, select the fields whose data you want to save in SEG-Y format.

Set the following parameters:

Click Apply.

Repeat all the steps with the same parameters for each field in the drop-down list (velocity, principal stress).

The receiver lines are highlighted on the model in yellow when clicked in the corresponding section of the Model Tree.

Starting calculation

1. Set the type of problem you want to solve.

On the command bar, select the calculation settings module (Mode - Calculation Settings, Calculation Settings - Transient Analysis, Transient - General).

Set the following calculation parameters:

Click Apply.

2. Go to the output fields settings section.

On the command bar, select the calculation settings module (Mode - Calculation Settings, Calculation Settings - Transient Analysis, Transient - Output Fields).

Specify:

Click Apply.

Click Start Calculation.

3. In the window that appears, select the directory in which the result will be saved, and enter the file name.

4. In the case of a successful calculation, the console displays the message: Calculation finished successfully at "date" "time".

Results analysis

1. Open the file with the results. There are three ways to do that.

2. To analyze the results, go to the FidesysViewer.

3. On the top bar, select the required result data to display. From the first drop-down list, select Stress, from the second - Mises.

4. Set the step 1 in the step viewer panel. You should see the plate image in the initial state. Next, click on Play  . You should see the propagation of stress over time.

5. Open the saved data in SEG-Y format.

To do this, go to Menu - File - Open. In the drop-down list of file types, select SEG-Y Files (* .sgy, * .segy). Specify the file to view test_Vy.sgs

 

Set the viewing direction along the Y axis.

The calculation results for displacement Uy in the SEG-Y format are visualized in the field of visualization.

Using Console Interface

Geometry creating, meshing, setting boundary conditions and materials can be performed using the console interface. Below is a link to the program code that allows you to perform the steps described above manual, you only need to specify the full path and name.

boussinesq_problem_segy