/[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 837 by gross, Mon Sep 4 22:37:25 2006 UTC revision 838 by gross, Tue Sep 5 22:45:55 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=[20,20,2]  dom.n=[50,30,2]
29  dom.order=2  dom.order=2
30    dom.integrationOrder=2
31    
32    
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.1  sq.dt_max=0.02
37    
38  iob=InterpolateOverBox(debug)  iob=InterpolateOverBox(debug)
39  iob.domain=Link(dom,"domain")  iob.domain=Link(dom,"domain")
# Line 71  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.00001  m.rel_tol=0.0001
76    
77  m.expansion_coefficient= 0.  m.expansion_coefficient= 0.
78  m.bulk_modulus=1000.  m.bulk_modulus=1000.
# Line 79  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/5.  m.shear_length=m.shear_modulus*0.75*10000.
84    
85    
86  ug=UpdateGeometry(debug)  ug=UpdateGeometry(debug)

Legend:
Removed from v.837  
changed lines
  Added in v.838

  ViewVC Help
Powered by ViewVC 1.1.26