/[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 1387 - (show annotations)
Fri Jan 11 07:45:26 2008 UTC (11 years, 9 months ago) by trankine
Original Path: temp/modellib/py_src/crustal/newcastle.py
File MIME type: text/x-python
File size: 2485 byte(s)
Restore the trunk that existed before the windows changes were committed to the (now moved to branches) old trunk.
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