/[escript]/trunk/modellib/py_src/crustal/newcastle.py
ViewVC logotype

Contents of /trunk/modellib/py_src/crustal/newcastle.py

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1384 - (show annotations)
Fri Jan 11 02:29:38 2008 UTC (11 years, 9 months ago) by phornby
Original Path: temp_trunk_copy/modellib/py_src/crustal/newcastle.py
File MIME type: text/x-python
File size: 2485 byte(s)
Make a temp copy of the trunk before checking in the windows changes


1 """
2 mining activities in modelframe
3
4 @var __author__: name of author
5 @var __licence__: licence agreement
6 @var __url__: url entry point on documentation
7 @var __version__: version
8 @var __date__: date of the version
9 """
10
11 __copyright__=""" Copyright (c) 2006 by ACcESS MNRF
12 http://www.access.edu.au
13 Primary Business: Queensland, Australia"""
14 __license__="""Licensed under the Open Software License version 3.0
15 http://www.opensource.org/licenses/osl-3.0.php"""
16 __author__="Lutz Gross, l.gross@uq.edu.au"
17 __url__="http://www.iservo.edu.au/esys/escript"
18 __version__="$Revision$"
19 __date__="$Date$"
20
21 d=True
22 from setups import MiningHistory, DensityChange, LinearElasticStressChange, CoulombFailureStress
23 from esys.modellib.geometry import FinleyReader,VectorConstrainerOverBox
24 from esys.modellib.input import Sequencer
25 from esys.escript.modelframe import Link,Simulation, DataSource
26 import numarray
27 from esys.modellib.visualization import WriteVTK
28
29 dom=FinleyReader(debug=d)
30 dom.source=DataSource("./newcastle_mines.msh","gmsh")
31 dom.tag_map_source=DataSource("./tags.xml", "ESysXML")
32
33 sq=Sequencer(debug=d)
34 sq.t=1840.
35 sq.t_end=2000.
36 sq.dt_max=100.
37
38 hist=MiningHistory(debug=d)
39 hist.history=DataSource("./newcastle_mining.xml")
40 hist.t=Link(sq,"t")
41
42 dens_dot=DensityChange(debug=d)
43 dens_dot.domain=Link(dom,"domain")
44 dens_dot.tag_map=Link(dom,"tag_map")
45 dens_dot.mass_rate=Link(hist,"mass_changes")
46
47 fix=VectorConstrainerOverBox(debug=d)
48 fix.domain=Link(dom,"domain")
49 fix.top=[False, False, False]
50 fix.bottom= [True, True, True]
51 fix.front= [False, False, False]
52 fix.back=[False, False, False]
53 fix.left=[False, False, False]
54 fix.right=[False, False, False]
55
56 el=LinearElasticStressChange(debug=d)
57 el.domain=Link(dom,"domain")
58 el.tag_map=Link(dom,"tag_map")
59 el.density=8e3*0
60 el.lame_lambda=1.7e11
61 el.lame_mu=1.7e11
62 el.location_of_fixed_displacement=Link(fix,"location_of_constraint")
63 el.density_rate=Link(dens_dot,"density_rate")
64
65 cfs=CoulombFailureStress(debug=d)
66 cfs.stress=Link(el,"stress")
67 cfs.friction_coefficient=0.
68 cfs.normal=numarray.array([-1,0,1])
69
70 vis=WriteVTK()
71 vis.t=Link(sq)
72 vis.data0=Link(el,"displacement")
73 vis.data1=Link(cfs,"cfs")
74 vis.dt=10.
75 vis.filename="out.xml"
76
77
78 s=Simulation([sq, hist, dens_dot, fix, el, vis], debug=d)
79 s.run()

  ViewVC Help
Powered by ViewVC 1.1.26