[CIG-SHORT] Question

Thu Apr 20 15:11:58 PDT 2017

Teng,

If you are not using PyLith v2.1.4 or v2.2.0, please use the latest 
version as you will likely get a better error message.

Some suggestions:

1. Look to see if you have extra stuff in your pylithapp.cfg file that 
is irrelevant once it reads in the secondary_fault.cfg file.

2. If you are using one of these versions, user --help-properties with 
the interface components to try to zero in on the problem.

For example,
pylith secondary_fault.cfg 
--timedependent.interfaces.fault1.quadrature.cell.help-properties

3. If you are using v2.2.0, then use the pylith parameter viewer to look 
at the hierarchy of components and search for quad_order in the fully 
expanded tree to look for issues.

Regards,
Brad

On 4/19/17 11:13 AM, Li, Teng wrote:
> 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>
>
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>