[CIG-SHORT] Question
Li, Teng
tengli2 at illinois.edu
Thu Apr 20 17:19:03 PDT 2017
Hi Brad,
Thanks! I will try this.
Best,
Teng
Sent from my iPhone
> On Apr 20, 2017, at 5:12 PM, Brad Aagaard <baagaard at usgs.gov> wrote:
>
> 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
>>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
More information about the CIG-SHORT
mailing list