/[escript]/trunk/modellib/py_src/crustal/newcastle.py
ViewVC logotype

Annotation of /trunk/modellib/py_src/crustal/newcastle.py

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1388 - (hide annotations)
Fri Jan 11 07:45:58 2008 UTC (11 years, 9 months ago) by trankine
File MIME type: text/x-python
File size: 2485 byte(s)
And get the *(&(*&(* name right
1 gross 986 """
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 gross 1000 from setups import MiningHistory, DensityChange, LinearElasticStressChange, CoulombFailureStress
23 gross 986 from esys.modellib.geometry import FinleyReader,VectorConstrainerOverBox
24     from esys.modellib.input import Sequencer
25     from esys.escript.modelframe import Link,Simulation, DataSource
26 gross 1000 import numarray
27     from esys.modellib.visualization import WriteVTK
28 gross 986
29     dom=FinleyReader(debug=d)
30     dom.source=DataSource("./newcastle_mines.msh","gmsh")
31 gross 997 dom.tag_map_source=DataSource("./tags.xml", "ESysXML")
32 gross 986
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 gross 997 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 gross 986
47     fix=VectorConstrainerOverBox(debug=d)
48     fix.domain=Link(dom,"domain")
49 gross 998 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 gross 986
56 gross 997 el=LinearElasticStressChange(debug=d)
57     el.domain=Link(dom,"domain")
58     el.tag_map=Link(dom,"tag_map")
59 gross 998 el.density=8e3*0
60     el.lame_lambda=1.7e11
61     el.lame_mu=1.7e11
62 gross 997 el.location_of_fixed_displacement=Link(fix,"location_of_constraint")
63     el.density_rate=Link(dens_dot,"density_rate")
64 gross 986
65 gross 1000 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 gross 986 s.run()

  ViewVC Help
Powered by ViewVC 1.1.26