/[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 1809 - (show annotations)
Thu Sep 25 06:43:44 2008 UTC (11 years ago) by ksteube
File MIME type: text/x-python
File size: 2674 byte(s)
Copyright updated in all python files

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__="http://www.uq.edu.au/esscc/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