[cig-commits] r19413 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv18

[brad at geodynamics.org]

Author: brad
Date: 2012-01-21 09:43:49 -0800 (Sat, 21 Jan 2012)
New Revision: 19413

Modified:
   short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv18/create_stresses.py
Log:
Correct stress setup.

Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv18/create_stresses.py
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv18/create_stresses.py	2012-01-21 02:11:00 UTC (rev 19412)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv18/create_stresses.py	2012-01-21 17:43:49 UTC (rev 19413)
@@ -50,8 +50,8 @@
 Pf = -densityW*gacc*z
 Szz = density*gacc*z
 
-Syy = mask1*(b22*(Szz+Pf)-Pf) + ~mask1*Szz
-Sxx = mask1*(b33*(Szz+Pf)-Pf) + ~mask1*Szz
+Syy = mask1*(b22*(Szz+Pf)-0*Pf) + ~mask1*Szz
+Sxx = mask1*(b33*(Szz+Pf)-0*Pf) + ~mask1*Szz
 Sxy = mask1*(-b23*(Szz+Pf)) + ~mask1*0.0
 
 Syz = 0.0*z
@@ -99,15 +99,20 @@
 ny = 0.0
 nz = 0.0
 
-# normal traction is Tx
-tractionNormal = Sxx*nx + Sxy*ny + Sxz*nz
+Tx = Sxx*nx + Sxy*ny + Sxz*nz
+Ty = Sxy*nx + Syy*ny + Syz*nz
+Tz = Sxz*nx + Syz*ny + Szz*nz
 
-# LL shear traction is Ty
-tractionShearLL = Sxy*nx + Syy*ny + Syz*nz
+tractionNormal = Tx*nx + Ty*ny + Tz*nz
+sx = 0.0
+sy = 1.0
+sz = 0.0
+tractionShearLL = Tx*sx + Ty*sy + Tz*sz
+sx = 0.0
+sy = 0.0
+sz = 1.0
+tractionShearUD = Tx*sx + Ty*sy + Tz*sz
 
-# Up-dip shear traction is Tz
-tractionShearUD = Sxz*nx + Syz*ny + Szz*nz
-
 writer.filename("%s_traction_main.spatialdb" % sim)
 dataOut = {'points': points,
            'coordsys': cs,
@@ -136,15 +141,20 @@
 ny = -sin(branchAngle)
 nz = 0.0
 
-# normal traction is Tx
-tractionNormal = Sxx*nx + Sxy*ny + Sxz*nz
+Tx = Sxx*nx + Sxy*ny + Sxz*nz
+Ty = Sxy*nx + Syy*ny + Syz*nz
+Tz = Sxz*nx + Syz*ny + Szz*nz
 
-# LL shear traction is Ty
-tractionShearLL = Sxy*nx + Syy*ny + Syz*nz
+tractionNormal = Tx*nx + Ty*ny + Tz*nz
+sx = sin(branchAngle)
+sy = cos(branchAngle)
+sz = 0.0
+tractionShearLL = Tx*sx + Ty*sy + Tz*sz
+sx = 0.0
+sy = 0.0
+sz = 1.0
+tractionShearUD = Tx*sx + Ty*sy + Tz*sz
 
-# Up-dip shear traction is Tz
-tractionShearUD = Sxz*nx + Syz*ny + Szz*nz
-
 writer.filename("%s_traction_branch.spatialdb" % sim)
 dataOut = {'points': points,
            'coordsys': cs,