test/python/drucker_prager.py

# $Id$

__copyright__="""  Copyright (c) 2006 by ACcESS MNRF
                    http://www.access.edu.au
                Primary Business: Queensland, Australia"""
__license__="""Licensed under the Open Software License version 3.0
             http://www.opensource.org/licenses/osl-3.0.php"""

# run by scons build/posix/modellib/test/python/drucker_prager.passed from check out cd.

from esys.modellib.mechanics import DruckerPrager
from esys.escript.modelframe import Link,Simulation
from esys.modellib.geometry import RectangularDomain, VectorConstrainer
from esys.modellib.input import Sequencer, InterpolateOverBox

debug=True

dom=RectangularDomain(debug)
dom.order=1


sq=Sequencer(debug)
sq.t=0
sq.t_end=1.
sq.dt_max=0.1

iob=InterpolateOverBox(debug)
iob.domain=Link(dom,"domain")
iob.value_left_bottom_front=[1.,0.,0.]
iob.value_right_bottom_front=[0.,0.,0.]
iob.value_left_bottom_back=[1.,0.,0.]
iob.value_right_bottom_back=[0.,0.,0.]
iob.value_left_top_front=[1.,0.,0.]
iob.value_right_top_front=[0.,0.,0.]
iob.value_left_top_back=[1.,0.,0.]
iob.value_right_top_back=[0.,0.,0.]

m=DruckerPrager(debug)
m.domain=Link(dom,"domain")

cv=VectorConstrainer(debug)
cv.domain=Link(dom,"domain")
cv.value=Link(iob,"out")
cv.left=[True, False, False] 
cv.top= [False, False, False]
cv.bottom= [False, False, False]
cv.front= [False, False, False]
cv.back= [False, False, False]
cv.right= [True, True, True]

m.prescribed_velocity=Link(cv,"value_of_constraint")
m.location_prescribed_velocity=Link(cv,"location_of_constraint")

dom.displacement=Link(m,"displacement")

s=Simulation([sq,dom,cv,m],debug=True)
s.run()
1	# $Id$
2
3	__copyright__=""" Copyright (c) 2006 by ACcESS MNRF
4	http://www.access.edu.au
5	Primary Business: Queensland, Australia"""
6	__license__="""Licensed under the Open Software License version 3.0
7	http://www.opensource.org/licenses/osl-3.0.php"""
8
9	# run by scons build/posix/modellib/test/python/drucker_prager.passed from check out cd.
10
11	from esys.modellib.mechanics import DruckerPrager
12	from esys.escript.modelframe import Link,Simulation
13	from esys.modellib.geometry import RectangularDomain, VectorConstrainer
14	from esys.modellib.input import Sequencer, InterpolateOverBox
15
16	debug=True
17
18	dom=RectangularDomain(debug)
19	dom.order=1
20
21
22	sq=Sequencer(debug)
23	sq.t=0
24	sq.t_end=1.
25	sq.dt_max=0.1
26
27	iob=InterpolateOverBox(debug)
28	iob.domain=Link(dom,"domain")
29	iob.value_left_bottom_front=[1.,0.,0.]
30	iob.value_right_bottom_front=[0.,0.,0.]
31	iob.value_left_bottom_back=[1.,0.,0.]
32	iob.value_right_bottom_back=[0.,0.,0.]
33	iob.value_left_top_front=[1.,0.,0.]
34	iob.value_right_top_front=[0.,0.,0.]
35	iob.value_left_top_back=[1.,0.,0.]
36	iob.value_right_top_back=[0.,0.,0.]
37
38	m=DruckerPrager(debug)
39	m.domain=Link(dom,"domain")
40
41	cv=VectorConstrainer(debug)
42	cv.domain=Link(dom,"domain")
43	cv.value=Link(iob,"out")
44	cv.left=[True, False, False]
45	cv.top= [False, False, False]
46	cv.bottom= [False, False, False]
47	cv.front= [False, False, False]
48	cv.back= [False, False, False]
49	cv.right= [True, True, True]
50
51	m.prescribed_velocity=Link(cv,"value_of_constraint")
52	m.location_prescribed_velocity=Link(cv,"location_of_constraint")
53
54	dom.displacement=Link(m,"displacement")
55
56	s=Simulation([sq,dom,cv,m],debug=True)
57	s.run()