[CIG-SHORT] Question

Wed Apr 19 10:13:12 PDT 2017

Hi Brad,
I double checked the script again, I have set interfaces=[fault,fault1] in the previous part.
And also, I have

[pylithapp.timedependent.interfaces]
# Change fault to dynamic fault interface.
fault = pylith.faults.FaultCohesiveDyn
fault1 = pylith.faults.FaultCohesiveDyn

As we can see in the attachment, I think the problem is just in this part. I wonder whether there should be some modification if the fault and fault1 have a shared nodes?
Best,
Teng

Teng Li

Master Candidate in Structures

Department of Civil and Environmental Engineering

University of Illinois at Urbana-Champaign

205 North Mathews Ave, Urbana, IL. 61801

Phone:(217)8196210, Email: tengli2 at illinois.edu

________________________________
From: CIG-SHORT [cig-short-bounces at geodynamics.org] on behalf of Charles Williams [willic3 at gmail.com]
Sent: Wednesday, April 19, 2017 3:01 AM
To: cig-short at geodynamics.org
Subject: Re: [CIG-SHORT] Question

Hi Teng Li,

I can’t be sure, but are you defining your interfaces correctly?  Looking at one of your previous .cfg files I didn’t see where you were doing that.  You should have something like:

[pylithapp.timedependent]
interfaces = [fault,fault1]

Cheers,
Charles

On 19/04/2017, at 6:20 PM, Li, Teng <tengli2 at illinois.edu<mailto:tengli2 at illinois.edu>> wrote:

Hi Brad,
After testing, I find the error only occurs in fault1 part as follows:

# ----------------------------------------------------------------------
# fault1
# ----------------------------------------------------------------------
[pylithapp.timedependent.interfaces.fault1]

# Material id
id = 101

# Name of nodeset from CUBIT
label = fault1
edge= fault_edge

# Specify the basis functions and quadrature:
# 1-D Lagrange cell in 2-D space with 2nd order quadrature
quadrature.cell = pylith.feassemble.FIATLagrange
quadrature.cell.dimension = 1
quadrature.cell.quad_order = 2

# slip-weakening friction
friction = pylith.friction.SlipWeakening
friction.label = Slip weakening

# Set the slip-weakening friction model parameters.
#  static coefficient of friction: 0.6
#  dynamic coefficient of friction: 0.3
#  slip-weakening parameter: 0.2 m
#  cohesion: 0 Pa
friction.db_properties = spatialdata.spatialdb.UniformDB
friction.db_properties.label = Slip weakening
friction.db_properties.values = [static-coefficient,dynamic-coefficient,slip-weakening-parameter,cohesion]
friction.db_properties.data = [0.6,0.3,0.2*m,0.0*Pa]
# Specify the initial tractions and a temporal perturbation on the
# fault using a uniform DB.
#
# Initial
#   shear: 75.0 MPa (right-lateral)
#   normal 120 MPa (compressive)
# Perturbation
#   shear: 25.0 MPa (right-lateral)
traction_perturbation = pylith.faults.TractPerturbation

[pylithapp.timedependent.interfaces.fault1.traction_perturbation]
#db_initial = spatialdata.spatialdb.UniformDB
#db_initial.label = Initial fault1 tractions
#db_initial.values = [traction-shear,traction-normal]
#db_initial.data = [-75.0*MPa, -120.0*MPa]
db_initial = spatialdata.spatialdb.SimpleDB
db_initial.label = Initial fault1 tractions

# Spatial database with physical properties for elastic material
db_initial.iohandler.filename = traction1.spatialdb

And I find that the ultimate reason is in the quadrature.cell part.

{command line}::
pyre.inventory(error)
-- 'quadrature.cell'
pylithapp: configuration error(s)

It seems in fault1, I can not use:

quadrature.cell = pylith.feassemble.FIATLagrange

Do you have any suggestions on this?
Best,
Teng

Teng Li

Master Candidate in Structures

Department of Civil and Environmental Engineering

University of Illinois at Urbana-Champaign

205 North Mathews Ave, Urbana, IL. 61801

Phone:(217)8196210, Email: tengli2 at illinois.edu<mailto:tengli2 at illinois.edu>

________________________________________
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, April 18, 2017 9:58 PM
To: cig-short at geodynamics.org<mailto:cig-short at geodynamics.org>
Subject: Re: [CIG-SHORT] Question

On 4/18/17 8:50 PM, Li, Teng wrote:
Hi,
I have a question about quadrature.cell.quad_order. In my case, I have a
main fault and a secondary fault. And these 2 faults meet at one node. I
create a nodeset named fault_edge for this edge. And in the .cfg file. I
write         edge=fault_edge only in the fault1 part, which is the
secondary fault. Should I also write edge=fault_edge in the fault part?

You should use edge = NODESET_NAME with each fault that has buried edges
(where NODESET_NAME is the name of the nodeset with the nodes on all
buried edges for that fault), that is, edges that end inside the domain.
Faults in which the all edges are at boundaries of the domain are the
only faults that do not need to use edge = NODESET.

And also, after running the code, I get the *_error:_*

{command line}::

pyre.inventory(error)

-- 'quadrature.cell.quad_order'

pylithapp: configuration error(s)

I am not sure why you are seeing this. The error message seems to
indicate you have specified the quad_oder on the command line.

Regards,
Brad

_______________________________________________
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
_______________________________________________
CIG-SHORT mailing list
CIG-SHORT at geodynamics.org
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/ 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/20170419/535d94a1/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Problem.png
Type: image/png
Size: 1218599 bytes
Desc: Problem.png
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20170419/535d94a1/attachment-0001.png>