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

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

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 999 by gross, Mon Feb 26 07:09:37 2007 UTC revision 1000 by gross, Wed Feb 28 00:34:42 2007 UTC
# Line 19  __version__="$Revision$" Line 19  __version__="$Revision$"
19  __date__="$Date$"  __date__="$Date$"
20    
21  d=True  d=True
22  from setups import MiningHistory, DensityChange, LinearElasticStressChange  from setups import MiningHistory, DensityChange, LinearElasticStressChange, CoulombFailureStress
23  from esys.modellib.geometry import FinleyReader,VectorConstrainerOverBox  from esys.modellib.geometry import FinleyReader,VectorConstrainerOverBox
24  from esys.modellib.input import Sequencer  from esys.modellib.input import Sequencer
25  from esys.escript.modelframe import Link,Simulation, DataSource  from esys.escript.modelframe import Link,Simulation, DataSource
26    import numarray
27    from esys.modellib.visualization import WriteVTK
28    
29  dom=FinleyReader(debug=d)  dom=FinleyReader(debug=d)
30  dom.source=DataSource("./newcastle_mines.msh","gmsh")  dom.source=DataSource("./newcastle_mines.msh","gmsh")
# Line 60  el.lame_mu=1.7e11 Line 62  el.lame_mu=1.7e11
62  el.location_of_fixed_displacement=Link(fix,"location_of_constraint")  el.location_of_fixed_displacement=Link(fix,"location_of_constraint")
63  el.density_rate=Link(dens_dot,"density_rate")  el.density_rate=Link(dens_dot,"density_rate")
64    
65  s=Simulation([sq, hist, dens_dot, fix, el], debug=d)  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()  s.run()

Legend:
Removed from v.999  
changed lines
  Added in v.1000

  ViewVC Help
Powered by ViewVC 1.1.26