[CIG-SHORT] fault opening near tip problem

[Charles Williams]

It sounds like he is doing prescribed slip, so there will just be a linear solve.  I think it would also be useful to see his mesh.

Cheers,
Charles


> On 11/05/2016, at 3:46 AM, Brad Aagaard <baagaard at usgs.gov> wrote:
> 
> 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
>> 
> 
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short

Charles A. Williams
Scientist
GNS Science
1 Fairway Drive, Avalon
PO Box 30368
Lower Hutt  5040
New Zealand
ph (office): 0064-4570-4566
fax (office): 0064-4570-4600
C.Williams at gns.cri.nz