/[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 1000 - (show annotations)
Wed Feb 28 00:34:42 2007 UTC (12 years, 6 months ago) by gross
File MIME type: text/x-python
File size: 2485 byte(s)
coulomb stress failure added
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