/[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 840 by gross, Tue Sep 5 22:45:55 2006 UTC revision 841 by gross, Wed Sep 6 09:48:37 2006 UTC
# Line 25  debug=True Line 25  debug=True
25    
26  dom=RectangularDomain(debug)  dom=RectangularDomain(debug)
27  dom.l=[1.,1.,1.]  dom.l=[1.,1.,1.]
28  dom.n=[50,30,2]  dom.n=[50,50,2]
29  dom.order=2  dom.order=2
30  dom.integrationOrder=2  dom.integrationOrder=2
31    
# Line 33  dom.integrationOrder=2 Line 33  dom.integrationOrder=2
33  sq=Sequencer(debug)  sq=Sequencer(debug)
34  sq.t=0  sq.t=0
35  sq.t_end=0.8  sq.t_end=0.8
36  sq.dt_max=0.02  sq.dt_max=0.03
37    
38  iob=InterpolateOverBox(debug)  iob=InterpolateOverBox(debug)
39  iob.domain=Link(dom,"domain")  iob.domain=Link(dom,"domain")
# Line 72  cv.back= [False, False, False] Line 72  cv.back= [False, False, False]
72  m.velocity=Link(iob,"out")  m.velocity=Link(iob,"out")
73  m.prescribed_velocity=Link(cv,"value_of_constraint")  m.prescribed_velocity=Link(cv,"value_of_constraint")
74  m.location_prescribed_velocity=Link(cv,"location_of_constraint")  m.location_prescribed_velocity=Link(cv,"location_of_constraint")
75  m.rel_tol=0.0001  m.rel_tol=0.01
76    
77  m.expansion_coefficient= 0.  m.expansion_coefficient= 0.
78  m.bulk_modulus=1000.  m.bulk_modulus=1000.
# Line 80  m.shear_modulus=1. Line 80  m.shear_modulus=1.
80  m.plastic_stress=0.  m.plastic_stress=0.
81  m.friction_parameter=0.  m.friction_parameter=0.
82  m.dilatancy_parameter=0.  m.dilatancy_parameter=0.
83  m.shear_length=m.shear_modulus*0.75*10000.  m.shear_length=m.shear_modulus*0.75*1.
84    
85    
86  ug=UpdateGeometry(debug)  ug=UpdateGeometry(debug)
87  ug.domain=Link(dom,"domain")  ug.domain=Link(dom,"domain")
88  ug.displacement=Link(m,"displacement")  ug.displacement=Link(m,"displacement")
89    
90  vis=WriteVTK()  vis=WriteVTK(debug)
91  vis.t=Link(sq)  vis.t=Link(sq)
92  vis.scalar=Link(m,"plastic_stress")  vis.scalar=Link(m,"plastic_stress")
93  vis.vector=Link(m,"velocity")  vis.vector=Link(m,"velocity")

Legend:
Removed from v.840  
changed lines
  Added in v.841

  ViewVC Help
Powered by ViewVC 1.1.26