/[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 997 - (show annotations)
Mon Feb 26 06:31:45 2007 UTC (13 years, 6 months ago) by gross
File MIME type: text/x-python
File size: 2130 byte(s)
small small fixes in modellib
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.value=0.
48 fix.top=False
49 fix.bottom=True
50 fix.front=False
51 fix.back=False
52 fix.left=False
53 fix.right=False
54 # dens_dot.density_rate=Link(hist,"mass_changes")
55
56 el=LinearElasticStressChange(debug=d)
57 el.domain=Link(dom,"domain")
58 el.tag_map=Link(dom,"tag_map")
59 el.density=1.
60 el.lame_lambda=2.
61 el.lame_mu=1.
62 el.location_of_fixed_displacement=Link(fix,"location_of_constraint")
63 el.density_rate=Link(dens_dot,"density_rate")
64
65 # hist.domain=Link(dom,"domain")
66
67
68 s=Simulation([sq, hist, dens_dot, fix, el], debug=d)
69 s.run()

  ViewVC Help
Powered by ViewVC 1.1.26