/[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 833 by gross, Fri Sep 1 02:04:46 2006 UTC revision 834 by gross, Fri Sep 1 09:05:16 2006 UTC
# Line 24  except KeyError: Line 24  except KeyError:
24  debug=True  debug=True
25    
26  dom=RectangularDomain(debug)  dom=RectangularDomain(debug)
27  dom.l=[5.,1.,1.]  dom.l=[1.,1.,1.]
28  dom.n=[25,5,5]  dom.n=[20,20,2]
29  dom.order=1  dom.order=2
30    
31    
32  sq=Sequencer(debug)  sq=Sequencer(debug)
33  sq.t=0  sq.t=0
34  sq.t_end=4.0  sq.t_end=0.8
35  sq.dt_max=0.049/4  sq.dt_max=0.049/10.
36    
37  iob=InterpolateOverBox(debug)  iob=InterpolateOverBox(debug)
38  iob.domain=Link(dom,"domain")  iob.domain=Link(dom,"domain")
# Line 71  cv.back= [False, False, False] Line 71  cv.back= [False, False, False]
71  m.velocity=Link(iob,"out")  m.velocity=Link(iob,"out")
72  m.prescribed_velocity=Link(cv,"value_of_constraint")  m.prescribed_velocity=Link(cv,"value_of_constraint")
73  m.location_prescribed_velocity=Link(cv,"location_of_constraint")  m.location_prescribed_velocity=Link(cv,"location_of_constraint")
74  m.rel_tol=0.005  m.rel_tol=0.00001
75    
76  m.expansion_coefficient= 0.  m.expansion_coefficient= 0.
77  m.bulk_modulus=10000.  m.bulk_modulus=1000.
78  m.shear_modulus=1.  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/2.  m.shear_length=m.shear_modulus*1000.
83    
84    
85  ug=UpdateGeometry(debug)  ug=UpdateGeometry(debug)

Legend:
Removed from v.833  
changed lines
  Added in v.834

  ViewVC Help
Powered by ViewVC 1.1.26