[CIG-SHORT] fault opening near tip problem

[Brad Aagaard]

Eren,

In order to answer your question, it would help to have more information.

What is the loading?

What do you mean by "fracture"? Are you prescribing slip or using a 
frictional interface? In either case, what are the interface conditions?

By "tip nodes" do you mean you have marked the buried edges so that they 
do not slip?

Do the linear and nonlinear solves converge?

Regards,
Brad


On 5/10/16 8:04 AM, Eren Ucar wrote:
> Hi,
>
> I am currently working on fracture deformation modeling.  I am trying to
> learn how to use Pylith. I am willing to implement relative motion
> across the fault surface(cohesive cells). I have started with 2D
> implementation, 10mx10m square. I added the picture of the domain
> (domain.pdf).  There is a fracture in the middle of domain (green line).
> I wanted to implement fault opening (0.01m). The fracture doesn’t cut
> through the domain so I introduced the tip nodes as fault edge. The
> boundary conditions everywhere are zero with constrained in both x and y
> direction. I attached the figure of the deformation in the domain
> (deformation_pylith.pdf). It is a vector drawing (precisely; matlab
> quiver). I drew the fracture again with green line. My question is; as
> you can see from the picture, there are different behaviors near the
> tips. I circled the regions with yellow. Tips are not displacing which
> is something I expected. But there is a different behavior near the tip
> that I couldn’t understand.
>
> Is this something normal? Maybe I made an implementation mistake.
>
> Ps1: When I fine the mesh I still get the same behavior.
>
> Ps2: When the fracture cuts through the domain, solution seems as I
> expected. No problem.
>
> Ps3: I write my own ascii file for mesh.
>
> I would appreciate any help.
>
> Regards,
>
> Eren
>
>
>
>
>
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>