[CIG-SHORT] Question

[Li, Teng]

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



________________________________________
From: CIG-SHORT [cig-short-bounces at geodynamics.org] on behalf of Brad Aagaard [baagaard at usgs.gov]
Sent: Tuesday, April 18, 2017 9:58 PM
To: 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
http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short