/[escript]/trunk/doc/examples/usersguide/heatedblock.py
ViewVC logotype

Diff of /trunk/doc/examples/usersguide/heatedblock.py

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

revision 2502 by gross, Tue Jun 30 05:49:22 2009 UTC revision 2513 by gross, Thu Jul 2 08:42:11 2009 UTC
# Line 57  u=mypde.getSolution() Line 57  u=mypde.getSolution()
57  g=grad(u)  g=grad(u)
58  sigma=mu*(g+transpose(g))+lam*trace(g)*kronecker(mydomain)-sigma0  sigma=mu*(g+transpose(g))+lam*trace(g)*kronecker(mydomain)-sigma0
59  sigma_mises=sqrt(((sigma[0,0]-sigma[1,1])**2+(sigma[1,1]-sigma[2,2])**2+ \  sigma_mises=sqrt(((sigma[0,0]-sigma[1,1])**2+(sigma[1,1]-sigma[2,2])**2+ \
60                    (sigma[2,2]-sigma[0,0])**2)/6. \                    (sigma[2,2]-sigma[0,0])**2)/2. \
61                     +sigma[0,1]**2 + sigma[1,2]**2 + sigma[2,0]**2)                     +3*(sigma[0,1]**2 + sigma[1,2]**2 + sigma[2,0]**2))
62  #... output ...  #... output ...
63  saveVTK("deform.xml",disp=u,stress=sigma_mises)  saveVTK("deform.xml",disp=u,stress=sigma_mises)
64    

Legend:
Removed from v.2502  
changed lines
  Added in v.2513

  ViewVC Help
Powered by ViewVC 1.1.26