[CIG-SHORT] Make two independent rough surfaces contact each other

[Charles Williams]

Hi Ekaterina,

Brad means that you need matching nodes (vertices) on either side of the fault, so that the topology matches across the fault.  PyLith isn’t set up to handle contact problems.

Cheers,
Charles


> On 17/01/2018, at 10:44 AM, Ekaterina Bolotskaya <bolee at mit.edu> wrote:
> 
> Dear Brad,
> 
> Two contacting bodies would work for my purposes as well.
> Could you tell me how to implement that? I can't find anything in the manual.
> 
> Thanks!
> 
> 
> Best regards,
> Ekaterina Bolotskaya
> 
> PhD student in Geophysics,
> Earth, Atmospheric and Planetary Science Department,
> Massachusetts Institute of Technology
> E-mail. bolee at mit.edu <mailto:bolee at mit.edu>
> Mob. +1 (857) 284-2805 <tel:%2B1%20%28857%29%20284-2805>
>          +7 (963) 995-36-33 <tel:%2B7%20%28963%29%20995-36-33>
> From: CIG-SHORT [cig-short-bounces at geodynamics.org <mailto:cig-short-bounces at geodynamics.org>] on behalf of Brad Aagaard [baagaard at usgs.gov <mailto:baagaard at usgs.gov>]
> Sent: Tuesday, January 16, 2018 15:32
> To: cig-short at geodynamics.org <mailto:cig-short at geodynamics.org>
> Subject: Re: [CIG-SHORT] Make two independent rough surfaces contact each other
> 
> PyLith's fault implementation assumes that the geometry of the fault surface is the same on the two opposing sides of the fault. To have two opposing sides with different geometry would generally involve a contact condition between two bodies, not a single internal surface that PyLith uses for the fault.
> 
> Regards,
> Brad
> 
> 
> 
> On Jan 16, 2018 10:26 AM, Ekaterina Bolotskaya <bolee at mit.edu <mailto:bolee at mit.edu>> wrote:
> Dear Pylith developers,
> 
> I have a 2D simulation of a rough fracture. It consists of 2 separate surfaces (upper and lower), one edge of which has a rough profile (bottom edge of the upper surface and upper edge of the lower surface). The rough edges are not initially in contact. 
> I would like to load the upper and the lower sides of the model to bring the asperities in contact with each other and look at plastic deformation. 
> 
> When I'm trying to do it Pylith moves both parts towards each other as if they are allowed to overlap (no contact, no deformation at the asperities).
> Is there a way to let it know the two surfaces are supposed to contact each other.
> 
> Thanks!
> 
> 
> Best regards,
> Ekaterina Bolotskaya
> 
> PhD student in Geophysics,
> Earth, Atmospheric and Planetary Science Department,
> Massachusetts Institute of Technology
> E-mail. bolee at mit.edu <mailto:bolee at mit.edu>
> Mob. +1 (857) 284-2805 <tel:%2B1%20%28857%29%20284-2805>
>          +7 (963) 995-36-33 <tel:%2B7%20%28963%29%20995-36-33>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org <mailto:CIG-SHORT at geodynamics.org>
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short>
Charles Williams I Geodynamic Modeler
GNS Science I Te Pῡ Ao
1 Fairway Drive, Avalon 5010, PO Box 30368, Lower Hutt 5040, New Zealand
Ph 0064-4-570-4566 I Mob 0064-22-350-7326 I Fax 0064-4-570-4600
http://www.gns.cri.nz/ <http://www.gns.cri.nz/> I Email: C.Williams at gns.cri.nz <mailto:your.email at gns.cri.nz>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20180117/daed9da3/attachment-0001.html>