Axial and Shear Deformation (2D Box)#

The files are in the directory examples/box-2d. The files and directories for this set of examples includes:

README.md

README file containing a brief description of the various examples.

*.cfg

PyLith parameter files.

*.mesh

Finite-element mesh files generated manually using a text editor.

*.spatialdb

Spatial database filesFiles associated with the spatial databases.

viz

Directory containing ParaView Python scripts and other files for visualizing results.

output

Directory containing simulation output. It is created automatically when running the simulations.

Overview#

This suite of examples demonstrates some basic concepts of using PyLith to solve the static and quasistatic boundary elasticity equation in a 2D box (Fig. 24) with uniform material properties. This example incrementally adds complexity through a series of steps:

Step 1

Axial extension with Dirichlet (displacement) boundary conditions.

Step 2

Shear deformation with Dirichlet (displacement) boundary conditions.

Step 3

Shear deformation with Dirichlet (displacement) and Neumann (traction) boundary conditions.

Step 4

Same as Step 2 but with initial conditions equal to the analytical solution.

Step 5

Shear deformation with time-dependent Dirichlet (displacement) and Neumann (traction) boundary conditions.

Diagram of geometry for 2D box.

Fig. 24 Diagram of geometry for 2D box that extends from -6 km to +6 km in the x direction and from -16 km to 0 km in the y direction. We refer to the domain boundaries using the names shown in the diagram.#

Example Workflow#