[cig-commits] r19604 - short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d
brad at geodynamics.org
brad at geodynamics.org
Wed Feb 8 18:07:38 PST 2012
Author: brad
Date: 2012-02-08 18:07:38 -0800 (Wed, 08 Feb 2012)
New Revision: 19604
Modified:
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/createbc.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/gradient.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_050m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.jou
short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.jou
Log:
Added interface to geometry associated with discontinuity in stress field in tpv13.
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/createbc.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/createbc.jou 2012-02-09 01:15:50 UTC (rev 19603)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/createbc.jou 2012-02-09 02:07:38 UTC (rev 19604)
@@ -1,44 +1,43 @@
# ----------------------------------------------------------------------
# Create blocks for materials
# ----------------------------------------------------------------------
-block 1 surface 8 9
+block 1 surface 10 11 12 13
block 1 name "elastic"
# ----------------------------------------------------------------------
# Create nodeset for fault
# ----------------------------------------------------------------------
-group "fault" add node in curve fault_surface
+group "fault" add node in curve fault_surface fault_surface at A
nodeset 20 group fault
nodeset 20 name "fault"
# ----------------------------------------------------------------------
# Create nodeset for +x face
# ----------------------------------------------------------------------
-group "face_xpos" add node in curve 14
+group "face_xpos" add node in curve 30 31
nodeset 10 group face_xpos
nodeset 10 name "face_xpos"
# ----------------------------------------------------------------------
# Create nodeset for -x face
# ----------------------------------------------------------------------
-group "face_xneg" add node in curve 16
+group "face_xneg" add node in curve 26 28
nodeset 11 group face_xneg
nodeset 11 name "face_xneg"
# ----------------------------------------------------------------------
# Create nodeset for +y face
# ----------------------------------------------------------------------
-group "face_ypos" add node in curve 19
-group "face_ypos" add node in curve 20
+group "face_ypos" add node in curve 19 20
nodeset 12 group face_ypos
nodeset 12 name "face_ypos"
# ----------------------------------------------------------------------
# Create nodeset for -y face
# ----------------------------------------------------------------------
-group "face_yneg" add node in curve 18
-group "face_yneg" add node in curve 21
+group "face_yneg" add node in curve 18 21
nodeset 13 group face_yneg
nodeset 13 name "face_yneg"
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry.jou 2012-02-09 01:15:50 UTC (rev 19603)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/geometry.jou 2012-02-09 02:07:38 UTC (rev 19604)
@@ -34,10 +34,12 @@
#
#{faultDipAngle=60.0}
#{faultWidth=15.0*km}
+#
+#{xoffset=0.5*faultWidth*cosd(faultDipAngle)}
brick x {blockWidth} y {blockHeight} z {blockWidth}
# Translate block so the top is at z=0
-volume 1 move x {-0.5*faultWidth*cosd(faultDipAngle)} y {-0.5*blockHeight} z 0
+volume 1 move x {-xoffset} y {-0.5*blockHeight} z 0
create midsurface volume 1 surface 1 2
delete volume 1
@@ -51,6 +53,11 @@
split curve 17 at position {-faultWidth*cosd(faultDipAngle)} {-faultWidth*sind(faultDipAngle)} 0
curve 23 name "fault_surface"
+# Create interface at 13.8 km down-dip
+# to honor discontinuity in stress in TPV13
+split surface 9 across location position {-blockWidth/2-xoffset} {sind(faultDipAngle)*-13.8*km} 0 location position {cosd(faultDipAngle)*-13.8*km} {sind(faultDipAngle)*-13.8*km} 0
+split surface 8 across location position {+blockWidth/2-xoffset} {sind(faultDipAngle)*-13.8*km} 0 location position {cosd(faultDipAngle)*-13.8*km} {sind(faultDipAngle)*-13.8*km} 0
+
# ----------------------------------------------------------------------
# Imprint all volumes, merging surfaces
# ----------------------------------------------------------------------
@@ -58,3 +65,4 @@
merge all
# End of file
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/gradient.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/gradient.jou 2012-02-09 01:15:50 UTC (rev 19603)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/gradient.jou 2012-02-09 02:07:38 UTC (rev 19604)
@@ -19,12 +19,18 @@
# dxC - size at vertex 13
#{dxStart=dx}{curveL=Length(22)}{execute("dxC="//fun_dx)}
-# dxD - size at vertex 9
-#{dxStart=dxA}{curveL=Length(14)}{execute("dxD="//fun_dx)}
+# dxD - size at vertex 18
+#{dxStart=dx}{curveL=Length(29)}{execute("dxD="//fun_dx)}
-# dxE - size at vertex 10
-#{dxStart=dxC}{curveL=Length(21)}{execute("dxE="//fun_dx)}
+# dxE - size at vertex 17
+#{dxStart=dx}{curveL=Length(24)}{execute("dxE="//fun_dx)}
+# dxF - size at vertex 9
+#{dxStart=dxD}{curveL=Length(30)}{execute("dxF="//fun_dx)}
+
+# dxG - size at vertex 10
+#{dxStart=dxC}{curveL=Length(21)}{execute("dxG="//fun_dx)}
+
# ----------------------------------------------------------------------
# Reset sizes
# ----------------------------------------------------------------------
@@ -32,27 +38,37 @@
surface all sizing function none
# Set size on faults
-curve fault_surface size {dx}
+curve fault_surface fault_surface at A size {dx}
# Set bias on curves extending from faults
curve 19 scheme bias fine size {dx} factor {bias_factor} start vertex 14
curve 20 scheme bias fine size {dx} factor {bias_factor} start vertex 14
curve 22 scheme bias fine size {dx} factor {bias_factor} start vertex 15
+curve 29 scheme bias fine size {dx} factor {bias_factor} start vertex 16
+curve 24 scheme bias fine size {dx} factor {bias_factor} start vertex 16
# A to D
-curve 14 scheme bias fine size {dxA} factor {bias_factor} start vertex 11
+curve 31 scheme bias fine size {dxA} factor {bias_factor} start vertex 11
-# C to E
+# E to B
+curve 26 scheme bias fine size {dxE} coarse size {dxB} start vertex 17
+
+# C to G
curve 21 scheme bias fine size {dxC} factor {bias_factor} start vertex 13
-# C to D
-curve 18 scheme bias fine size {dxC} coarse size {dxD} start vertex 13
+# D to F
+curve 30 scheme bias fine size {dxD} coarse size {dxF} start vertex 18
-# B to E
-curve 16 scheme bias fine size {dxB} coarse size {dxE} start vertex 12
+# C to F
+curve 18 scheme bias fine size {dxC} coarse size {dxF} start vertex 13
-surface 8 9 sizing function type bias start curve fault_surface factor {bias_factor}
+# E to G
+curve 28 scheme bias fine size {dxE} coarse size {dxG} start vertex 17
+surface 10 13 sizing function type bias start curve fault_surface factor {bias_factor}
+surface 11 12 sizing function type bias start vertex 16 factor {bias_factor}
+
# End of file
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_050m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_050m.jou 2012-02-09 01:15:50 UTC (rev 19603)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_050m.jou 2012-02-09 02:07:38 UTC (rev 19604)
@@ -15,8 +15,11 @@
playback 'gradient.jou'
# A to D
-curve 14 scheme bias fine size {dxA*1.01} factor {bias_factor} start vertex 11
+curve 31 scheme bias fine size {dxA*0.98} factor {bias_factor} start vertex 11
+# C to F
+curve 18 scheme bias fine size {dxC*0.95} coarse size {dxF} start vertex 13
+
# ----------------------------------------------------------------------
# Generate the mesh
# ----------------------------------------------------------------------
@@ -42,3 +45,4 @@
# End of file
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.jou 2012-02-09 01:15:50 UTC (rev 19603)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_100m.jou 2012-02-09 02:07:38 UTC (rev 19604)
@@ -14,8 +14,8 @@
surface all scheme pave
playback 'gradient.jou'
-# A to D
-curve 14 scheme bias fine size {dxA*1.01} factor {bias_factor} start vertex 11
+# D to F
+curve 30 scheme bias fine size {dxD*1.02} coarse size {dxF} start vertex 18
# ----------------------------------------------------------------------
# Generate the mesh
@@ -42,3 +42,4 @@
# End of file
+
Modified: short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.jou
===================================================================
--- short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.jou 2012-02-09 01:15:50 UTC (rev 19603)
+++ short/3D/PyLith/benchmarks/trunk/dynamic/scecdynrup/tpv210-2d/quad4_200m.jou 2012-02-09 02:07:38 UTC (rev 19604)
@@ -14,8 +14,8 @@
surface all scheme pave
playback 'gradient.jou'
-# A to D
-curve 14 scheme bias fine size {dxA*1.01} factor {bias_factor} start vertex 11
+# D to F
+curve 30 scheme bias fine size {dxD*1.02} coarse size {dxF} start vertex 18
# ----------------------------------------------------------------------
# Generate the mesh
@@ -41,3 +41,4 @@
# End of file
+
More information about the CIG-COMMITS
mailing list