/[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 997 - (hide annotations)
Mon Feb 26 06:31:45 2007 UTC (14 years, 9 months ago) by gross
File MIME type: text/x-python
File size: 2130 byte(s)
small small fixes in modellib
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 997 from setups import MiningHistory, DensityChange, LinearElasticStressChange
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    
27     dom=FinleyReader(debug=d)
28     dom.source=DataSource("./newcastle_mines.msh","gmsh")
29 gross 997 dom.tag_map_source=DataSource("./tags.xml", "ESysXML")
30 gross 986
31     sq=Sequencer(debug=d)
32     sq.t=1840.
33     sq.t_end=2000.
34     sq.dt_max=100.
35    
36     hist=MiningHistory(debug=d)
37     hist.history=DataSource("./newcastle_mining.xml")
38     hist.t=Link(sq,"t")
39    
40 gross 997 dens_dot=DensityChange(debug=d)
41     dens_dot.domain=Link(dom,"domain")
42     dens_dot.tag_map=Link(dom,"tag_map")
43     dens_dot.mass_rate=Link(hist,"mass_changes")
44 gross 986
45     fix=VectorConstrainerOverBox(debug=d)
46     fix.domain=Link(dom,"domain")
47     fix.value=0.
48     fix.top=False
49     fix.bottom=True
50     fix.front=False
51     fix.back=False
52     fix.left=False
53     fix.right=False
54 gross 997 # dens_dot.density_rate=Link(hist,"mass_changes")
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     el.density=1.
60     el.lame_lambda=2.
61     el.lame_mu=1.
62     el.location_of_fixed_displacement=Link(fix,"location_of_constraint")
63     el.density_rate=Link(dens_dot,"density_rate")
64 gross 986
65 gross 997 # hist.domain=Link(dom,"domain")
66    
67    
68     s=Simulation([sq, hist, dens_dot, fix, el], debug=d)
69 gross 986 s.run()

  ViewVC Help
Powered by ViewVC 1.1.26