search

geodynamics / pylith

PyLith is a finite element code for the solution of dynamic and quasi-static tectonic deformation problems.
Other
150 stars 96 forks source link

Problem in fault vertices (mesh) in a curved fault model #756

Open DianKusumawati opened 3 weeks ago

DianKusumawati commented 3 weeks ago

Dear Sir,

I tried to simulate the deformation of a non-planar (curved) fault (shown in the orange dots in the figures) due to a shearing boundary condition, using FaultCohesiveDyn in Pylith 2.2.1. I made the mesh using Cubit, and set the mesh to hexahedron, through the pave scheme, to match the FIATLagrange settings. But when I ran the simulation the code showed an error message: Fault vertices cannot be constrained.

I tried to edit the meshing by setting the scheme into the map scheme, but I found the same error: Fault vertices cannot be constrained.

Kindly suggest or point out what should I do for the solution. I attached the mesh and error message, as well as the corresponding files. Thank you.

The files pylithfiles.zip

The mesh with pave scheme

image image

The mesh with map scheme

image image
baagaard-usgs commented 2 weeks ago

The error is associated with the overlap in your fault and Dirichlet boundary conditions. As the error message states, you cannot constrain the displacement for a vertex on the fault. Either you need to adjust your geometry so that the fault does not extend all of the way to the boundary of your domain, or you need to exclude all vertices on the fault from your boundary condition.

DianKusumawati commented 2 weeks ago

Thank you for the solution. We'll try it. Thanks!

The error is associated with the overlap in your fault and Dirichlet boundary conditions. As the error message states, you cannot constrain the displacement for a vertex on the fault. Either you need to adjust your geometry so that the fault does not extend all of the way to the boundary of your domain, or you need to exclude all vertices on the fault from your boundary condition.