/[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 998 - (show annotations)
Mon Feb 26 07:09:37 2007 UTC (12 years, 6 months ago) by gross
File MIME type: text/x-python
File size: 2139 byte(s)
basic mine scenario
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
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
27 dom=FinleyReader(debug=d)
28 dom.source=DataSource("./newcastle_mines.msh","gmsh")
29 dom.tag_map_source=DataSource("./tags.xml", "ESysXML")
30
31 sq=Sequencer(debug=d)
32 sq.t=1840.
33 sq.t_end=2000.
34 sq.dt_max=100.
35
36 hist=MiningHistory(debug=d)
37 hist.history=DataSource("./newcastle_mining.xml")
38 hist.t=Link(sq,"t")
39
40 dens_dot=DensityChange(debug=d)
41 dens_dot.domain=Link(dom,"domain")
42 dens_dot.tag_map=Link(dom,"tag_map")
43 dens_dot.mass_rate=Link(hist,"mass_changes")
44
45 fix=VectorConstrainerOverBox(debug=d)
46 fix.domain=Link(dom,"domain")
47 fix.top=[False, False, False]
48 fix.bottom= [True, True, True]
49 fix.front= [False, False, False]
50 fix.back=[False, False, False]
51 fix.left=[False, False, False]
52 fix.right=[False, False, False]
53
54 el=LinearElasticStressChange(debug=d)
55 el.domain=Link(dom,"domain")
56 el.tag_map=Link(dom,"tag_map")
57 el.density=8e3*0
58 el.lame_lambda=1.7e11
59 el.lame_mu=1.7e11
60 el.location_of_fixed_displacement=Link(fix,"location_of_constraint")
61 el.density_rate=Link(dens_dot,"density_rate")
62
63 s=Simulation([sq, hist, dens_dot, fix, el], debug=d)
64 s.run()

  ViewVC Help
Powered by ViewVC 1.1.26