| 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() |
Parent Directory
| 
Revision Log