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#