/[escript]/trunk/finley/test/python/slip_stress.py
ViewVC logotype

Diff of /trunk/finley/test/python/slip_stress.py

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

revision 896 by gross, Thu Nov 9 05:46:31 2006 UTC revision 897 by gross, Thu Nov 9 07:37:40 2006 UTC
# Line 49  print "total ne = ",ne*ne_w*ne_w Line 49  print "total ne = ",ne*ne_w*ne_w
49    
50  fstart =  [lockToGrid(50000.0,width,ne_w), lockToGrid(40000.0,width,ne_w), lockToGrid(8000.,height,ne)]  fstart =  [lockToGrid(50000.0,width,ne_w), lockToGrid(40000.0,width,ne_w), lockToGrid(8000.,height,ne)]
51  fend =  [lockToGrid(50000.0,width,ne_w), lockToGrid(60000.0,width,ne_w), lockToGrid(20000.,height,ne)]  fend =  [lockToGrid(50000.0,width,ne_w), lockToGrid(60000.0,width,ne_w), lockToGrid(20000.,height,ne)]
52  s=[0.,0.,1.]  s=[0.,1.,0.]
53    
54  fstart =  [lockToGrid(30000.0,width,ne_w), lockToGrid(30000.0,width,ne_w), lockToGrid(20000.,height,ne)]  # fstart =  [lockToGrid(30000.0,width,ne_w), lockToGrid(30000.0,width,ne_w), lockToGrid(20000.,height,ne)]
55  fend =  [lockToGrid(70000.0,width,ne_w), lockToGrid(70000.0,width,ne_w), lockToGrid(20000.,height,ne)]  # fend =  [lockToGrid(70000.0,width,ne_w), lockToGrid(70000.0,width,ne_w), lockToGrid(20000.,height,ne)]
56  s=[1.,0.1,0.]  # s=[1.,0.,0.]
57    
58  dom=Brick(l0=width, l1=width, l2=height, n0=ne_w, n1=ne_w, n2=ne)  dom=Brick(l0=width, l1=width, l2=height, n0=ne_w, n1=ne_w, n2=ne)
59    
# Line 89  for i in range(d): Line 89  for i in range(d):
89        A[i,j,j,i] += mu        A[i,j,j,i] += mu
90        A[i,j,i,j] += mu        A[i,j,i,j] += mu
91        A[i,i,j,j] += lmbd        A[i,i,j,j] += lmbd
92  pde.setValue(A=A,q=mask,Y=-kronecker(Function(dom))[d-1]*g*rho,X=-sigma0)  pde.setValue(A=A,q=mask,Y=-kronecker(Function(dom))[d-1]*g*rho,X=sigma0)
93  u=pde.getSolution(verbose=True)  u=pde.getSolution(verbose=True)
94  g_s=grad(u)  g_s=grad(u)
95  saveVTK("dis.xml",disp=u,sigma=mu*symmetric(g_s)+lmbd*trace(g_s)*kronecker(d))  sigma=mu*symmetric(g_s)+lmbd*trace(g_s)*kronecker(d)
96    saveVTK("dis.xml",disp=u,sigma=sigma,cfs=sigma[0,1]-0.4*sigma[0,0])

Legend:
Removed from v.896  
changed lines
  Added in v.897

  ViewVC Help
Powered by ViewVC 1.1.26