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

1 ksteube 1809
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 gross 986 """
23     mining activities in modelframe
24 ksteube 1809
25 gross 986 @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 gross 1000 from setups import MiningHistory, DensityChange, LinearElasticStressChange, CoulombFailureStress
36 gross 986 from esys.modellib.geometry import FinleyReader,VectorConstrainerOverBox
37     from esys.modellib.input import Sequencer
38     from esys.escript.modelframe import Link,Simulation, DataSource
39 gross 1000 import numarray
40     from esys.modellib.visualization import WriteVTK
41 gross 986
42     dom=FinleyReader(debug=d)
43     dom.source=DataSource("./newcastle_mines.msh","gmsh")
44 gross 997 dom.tag_map_source=DataSource("./tags.xml", "ESysXML")
45 gross 986
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 gross 997 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 gross 986
60     fix=VectorConstrainerOverBox(debug=d)
61     fix.domain=Link(dom,"domain")
62 gross 998 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 gross 986
69 gross 997 el=LinearElasticStressChange(debug=d)
70     el.domain=Link(dom,"domain")
71     el.tag_map=Link(dom,"tag_map")
72 gross 998 el.density=8e3*0
73     el.lame_lambda=1.7e11
74     el.lame_mu=1.7e11
75 gross 997 el.location_of_fixed_displacement=Link(fix,"location_of_constraint")
76     el.density_rate=Link(dens_dot,"density_rate")
77 gross 986
78 gross 1000 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 gross 986 s.run()

  ViewVC Help
Powered by ViewVC 1.1.26