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

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

Parent Directory Parent Directory | Revision Log Revision Log


Revision 998 - (hide annotations)
Mon Feb 26 07:09:37 2007 UTC (12 years, 7 months ago) by gross
File MIME type: text/x-python
File size: 2139 byte(s)
basic mine scenario
1 gross 986 """
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 gross 997 from setups import MiningHistory, DensityChange, LinearElasticStressChange
23 gross 986 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 gross 997 dom.tag_map_source=DataSource("./tags.xml", "ESysXML")
30 gross 986
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 gross 997 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 gross 986
45     fix=VectorConstrainerOverBox(debug=d)
46     fix.domain=Link(dom,"domain")
47 gross 998 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 gross 986
54 gross 997 el=LinearElasticStressChange(debug=d)
55     el.domain=Link(dom,"domain")
56     el.tag_map=Link(dom,"tag_map")
57 gross 998 el.density=8e3*0
58     el.lame_lambda=1.7e11
59     el.lame_mu=1.7e11
60 gross 997 el.location_of_fixed_displacement=Link(fix,"location_of_constraint")
61     el.density_rate=Link(dens_dot,"density_rate")
62 gross 986
63 gross 997 s=Simulation([sq, hist, dens_dot, fix, el], debug=d)
64 gross 986 s.run()

  ViewVC Help
Powered by ViewVC 1.1.26