/[escript]/trunk/modellib/test/python/drucker_prager.py
ViewVC logotype

Diff of /trunk/modellib/test/python/drucker_prager.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 831 by gross, Thu Aug 31 08:47:40 2006 UTC revision 832 by gross, Fri Sep 1 02:04:46 2006 UTC
# Line 31  dom.order=1 Line 31  dom.order=1
31    
32  sq=Sequencer(debug)  sq=Sequencer(debug)
33  sq.t=0  sq.t=0
34  sq.t_end=3.0  sq.t_end=4.0
35  sq.dt_max=0.01  sq.dt_max=0.049/4
36    
37  iob=InterpolateOverBox(debug)  iob=InterpolateOverBox(debug)
38  iob.domain=Link(dom,"domain")  iob.domain=Link(dom,"domain")
# Line 79  m.shear_modulus=1. Line 79  m.shear_modulus=1.
79  m.plastic_stress=0.  m.plastic_stress=0.
80  m.friction_parameter=0.  m.friction_parameter=0.
81  m.dilatancy_parameter=0.  m.dilatancy_parameter=0.
82  m.shear_length=m.shear_modulus*2000.  m.shear_length=m.shear_modulus/2.
83    
84    
85  ug=UpdateGeometry(debug)  ug=UpdateGeometry(debug)
# Line 88  ug.displacement=Link(m,"displacement") Line 88  ug.displacement=Link(m,"displacement")
88    
89  vis=WriteVTK()  vis=WriteVTK()
90  vis.t=Link(sq)  vis.t=Link(sq)
91    vis.scalar=Link(m,"plastic_stress")
92  vis.vector=Link(m,"velocity")  vis.vector=Link(m,"velocity")
 # vis.scalar=Link(m,"plastic_stress")  
93  vis.tensor=Link(m,"stress")  vis.tensor=Link(m,"stress")
94  vis.dt=0.2  vis.dt=0.2
95  vis.filename=WORKDIR+"/temp.xml"  vis.filename=WORKDIR+"/temp.xml"
Legend:
Removed from v.831  
changed lines
  Added in v.832
  ViewVC Help
Powered by ViewVC 1.1.26