/[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 2344 - (show annotations)
Mon Mar 30 02:13:58 2009 UTC (10 years, 6 months ago) by jfenwick
File MIME type: text/x-python
File size: 2669 byte(s)
Change __url__ to launchpad site

1
2 ########################################################
3 #
4 # Copyright (c) 2003-2008 by University of Queensland
5 # Earth Systems Science Computational Center (ESSCC)
6 # http://www.uq.edu.au/esscc
7 #
8 # Primary Business: Queensland, Australia
9 # Licensed under the Open Software License version 3.0
10 # http://www.opensource.org/licenses/osl-3.0.php
11 #
12 ########################################################
13
14 __copyright__="""Copyright (c) 2003-2008 by University of Queensland
15 Earth Systems Science Computational Center (ESSCC)
16 http://www.uq.edu.au/esscc
17 Primary Business: Queensland, Australia"""
18 __license__="""Licensed under the Open Software License version 3.0
19 http://www.opensource.org/licenses/osl-3.0.php"""
20 __url__="https://launchpad.net/escript-finley"
21
22 """
23 mining activities in modelframe
24
25 @var __author__: name of author
26 @var __licence__: licence agreement
27 @var __url__: url entry point on documentation
28 @var __version__: version
29 @var __date__: date of the version
30 """
31
32 __author__="Lutz Gross, l.gross@uq.edu.au"
33
34 d=True
35 from setups import MiningHistory, DensityChange, LinearElasticStressChange, CoulombFailureStress
36 from esys.modellib.geometry import FinleyReader,VectorConstrainerOverBox
37 from esys.modellib.input import Sequencer
38 from esys.escript.modelframe import Link,Simulation, DataSource
39 import numarray
40 from esys.modellib.visualization import WriteVTK
41
42 dom=FinleyReader(debug=d)
43 dom.source=DataSource("./newcastle_mines.msh","gmsh")
44 dom.tag_map_source=DataSource("./tags.xml", "ESysXML")
45
46 sq=Sequencer(debug=d)
47 sq.t=1840.
48 sq.t_end=2000.
49 sq.dt_max=100.
50
51 hist=MiningHistory(debug=d)
52 hist.history=DataSource("./newcastle_mining.xml")
53 hist.t=Link(sq,"t")
54
55 dens_dot=DensityChange(debug=d)
56 dens_dot.domain=Link(dom,"domain")
57 dens_dot.tag_map=Link(dom,"tag_map")
58 dens_dot.mass_rate=Link(hist,"mass_changes")
59
60 fix=VectorConstrainerOverBox(debug=d)
61 fix.domain=Link(dom,"domain")
62 fix.top=[False, False, False]
63 fix.bottom= [True, True, True]
64 fix.front= [False, False, False]
65 fix.back=[False, False, False]
66 fix.left=[False, False, False]
67 fix.right=[False, False, False]
68
69 el=LinearElasticStressChange(debug=d)
70 el.domain=Link(dom,"domain")
71 el.tag_map=Link(dom,"tag_map")
72 el.density=8e3*0
73 el.lame_lambda=1.7e11
74 el.lame_mu=1.7e11
75 el.location_of_fixed_displacement=Link(fix,"location_of_constraint")
76 el.density_rate=Link(dens_dot,"density_rate")
77
78 cfs=CoulombFailureStress(debug=d)
79 cfs.stress=Link(el,"stress")
80 cfs.friction_coefficient=0.
81 cfs.normal=numarray.array([-1,0,1])
82
83 vis=WriteVTK()
84 vis.t=Link(sq)
85 vis.data0=Link(el,"displacement")
86 vis.data1=Link(cfs,"cfs")
87 vis.dt=10.
88 vis.filename="out.xml"
89
90
91 s=Simulation([sq, hist, dens_dot, fix, el, vis], debug=d)
92 s.run()

  ViewVC Help
Powered by ViewVC 1.1.26