/[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 871 by gross, Sun Sep 17 23:27:00 2006 UTC revision 872 by gross, Sun Oct 15 23:33:14 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=[1.,1.,1.]  dom.dim=3
28  dom.n=[50,50,2]  dom.l=[0.5,1.,1.]
29  dom.order=2  dom.n=[30,6,6]
30  dom.integrationOrder=2  dom.order=1
31    dom.integrationOrder=-1
32    
33    
34  sq=Sequencer(debug)  sq=Sequencer(debug)
35  sq.t=0  sq.t=0
36  sq.t_end=3  sq.t_end=0.8
37  sq.dt_max=0.3  sq.dt_max=100.
38    
39  iob=InterpolateOverBox(debug)  iob=InterpolateOverBox(debug)
40  iob.domain=Link(dom,"domain")  iob.domain=Link(dom,"domain")
41  iob.value_left_bottom_front=[1.,0.,0.]  iob.value_left_bottom_front=[-1.,0.,0.]
42  iob.value_right_bottom_front=[0.,0.,0.]  iob.value_right_bottom_front=[0.,0.,0.]
43  iob.value_left_bottom_back=[1.,0.,0.]  iob.value_left_bottom_back=[-1.,0.,0.]
44  iob.value_right_bottom_back=[0.,0.,0.]  iob.value_right_bottom_back=[0.,0.,0.]
45  iob.value_left_top_front=[1.,0.,0.]  iob.value_left_top_front=[-1.,0.,0.]
46  iob.value_right_top_front=[0.,0.,0.]  iob.value_right_top_front=[0.,0.,0.]
47  iob.value_left_top_back=[1.,0.,0.]  iob.value_left_top_back=[-1.,0.,0.]
48  iob.value_right_top_back=[0.,0.,0.]  iob.value_right_top_back=[0.,0.,0.]
49    
50  iob.value_left_bottom_front=[-1.,0.]  # iob.value_left_bottom_front=[-1.,0.]
51  iob.value_right_bottom_front=[0.,0.]  # iob.value_right_bottom_front=[0.,0.]
52  iob.value_left_bottom_back=[-1.,0.]  # iob.value_left_bottom_back=[-1.,0.]
53  iob.value_right_bottom_back=[0.,0.]  # iob.value_right_bottom_back=[0.,0.]
54  iob.value_left_top_front=[-1.,0.]  # iob.value_left_top_front=[-1.,0.]
55  iob.value_right_top_front=[0.,0.]  # iob.value_right_top_front=[0.,0.]
56  iob.value_left_top_back=[-1.,0.]  # iob.value_left_top_back=[-1.,0.]
57  iob.value_right_top_back=[0.,0.]  # iob.value_right_top_back=[0.,0.]
58    
59  m=DruckerPrager(debug)  m=DruckerPrager(debug)
60  m.domain=Link(dom,"domain")  m.domain=Link(dom,"domain")
# Line 63  cv.domain=Link(dom,"domain") Line 64  cv.domain=Link(dom,"domain")
64  cv.value=Link(iob,"out")  cv.value=Link(iob,"out")
65    
66  cv.left=[True, False, False]  cv.left=[True, False, False]
67  cv.right= [True, False, False]  cv.right= [True, True, True]
68  cv.bottom= [False, True, False]  cv.bottom= [False, False, False]
69  cv.top= [False, False, False]  cv.top= [False, False, False]
70  cv.front= [False, False, True]  cv.front= [False, False, False]
71  cv.back= [False, False, False]  cv.back= [False, False, False]
72    
73  m.velocity=Link(iob,"out")  m.velocity=Link(iob,"out")
# Line 80  m.shear_modulus=1. Line 81  m.shear_modulus=1.
81  m.plastic_stress=0.  m.plastic_stress=0.
82  m.friction_parameter=0.  m.friction_parameter=0.
83  m.dilatancy_parameter=0.  m.dilatancy_parameter=0.
84  m.shear_length=m.shear_modulus*0.75*1.  m.shear_length=m.shear_modulus*0.9
85    
86    
87  ug=UpdateGeometry(debug)  ug=UpdateGeometry(debug)
# Line 92  vis.t=Link(sq) Line 93  vis.t=Link(sq)
93  vis.scalar=Link(m,"plastic_stress")  vis.scalar=Link(m,"plastic_stress")
94  vis.vector=Link(m,"velocity")  vis.vector=Link(m,"velocity")
95  vis.tensor=Link(m,"stress")  vis.tensor=Link(m,"stress")
96  vis.dt=0.5  vis.dt=0.01
97  vis.filename=WORKDIR+"/temp.xml"  vis.filename=WORKDIR+"/temp.xml"
98    
99  s=Simulation([sq,m,vis],debug=True)  s=Simulation([sq,m,vis],debug=True)
Legend:
Removed from v.871  
changed lines
  Added in v.872
  ViewVC Help
Powered by ViewVC 1.1.26