[CIG-SHORT] Initial stress cell_data_fields settings and interpolation
Chunyu Liu (cliu5)
cliu5 at memphis.edu
Mon Nov 9 07:29:04 PST 2015
Hi Brad,
Thank you for replying.
Actually my initial stress data dimension is 3.
When I change data-dim=3 and use linear interpolation, it works on SimpleGridDB but not in SimpleDB.
However, SimpleDB works by using nearest interpolation.
Any idea about it ?
Thanks again.
This is the initial stress file:
// -*- C++ -*- (tell Emacs to use C++ mode for syntax highlighting)
//
// This spatial database specifies the initial stresses for the model
// The stresses increase linearly with depth.
//
#SPATIAL_GRID.ascii 1
SimpleGridDB {
num-values = 6 // number of stress components
// names of the stress components
value-names = stress-xx stress-yy stress-zz stress-xy stress-yz stress-xz
value-units = Pa Pa Pa Pa Pa Pa // units
// num-locs = 125
// data-dim = 3
num-x = 5
num-y = 5
num-z = 5
space-dim = 3
cs-data = cartesian {
to-meters = 1.0
space-dim = 3
}
}
//columns are
// (1) x coordinates
-222000.000000 -111000.000000 0.000000 111000.000000 222000.000000
// (2) y coordinates
-222000.000000 -111000.000000 0.000000 111000.000000 222000.000000
// (3) z coordinates
0.000000 -50000.000000 -100000.000000 -120000.000000 -160000.000000
// Columns are
// (1) x coordinates(m)
// (2) y coordinates(m)
// (3) z coordinates(m)
// (4) stress-xx (Pa)
// (5) stress-yy (Pa)
// (6) stress-zz (Pa)
// (7) stress-xy (Pa)
// (8) stress-yz (Pa)
// (9) stress-xz (Pa)
-222000.000000 -222000.000000 0.000000 -45681491.000000 -2434233.000000 48115724.000000 1073136.000000 0.000000 0.000000
-222000.000000 -111000.000000 0.000000 -48460100.000000 -1221094.000000 49681194.000000 1061051.000000 0.000000 0.000000
-222000.000000 0.000000 0.000000 -44910160.000000 2325826.400000 42584333.599999 2885474.000000 0.000000 0.000000
-222000.000000 111000.000000 0.000000 -45189450.000000 4476817.000000 40712633.000000 4592233.400000 0.000000 0.000000
-222000.000000 222000.000000 0.000000 -50044152.000000 2600954.400000 47443197.599999 5230918.000000 0.000000 0.000000
-111000.000000 -222000.000000 0.000000 -43048840.000000 -5970312.000000 49019152.000000 13113931.000000 0.000000 0.000000
-111000.000000 -111000.000000 0.000000 -47813031.000000 -2786867.000000 50599898.000000 15048702.000000 0.000000 0.000000
-111000.000000 0.000000 0.000000 -46725583.000000 3007694.400000 43717888.599999 11741930.000000 0.000000 0.000000
-111000.000000 111000.000000 0.000000 -45840312.000000 4800151.000000 41040161.000000 8963863.000000 0.000000 0.000000
-111000.000000 222000.000000 0.000000 -47930592.000000 7399439.000000 40531153.000000 8105111.000000 0.000000 0.000000
0.000000 -222000.000000 0.000000 -64667330.000000 4102524.000000 60564806.000000 29007560.000000 0.000000 0.000000
0.000000 -111000.000000 0.000000 -55029184.000000 8561188.000000 46467996.000000 18754142.000000 0.000000 0.000000
0.000000 0.000000 0.000000 -53050243.000000 5777051.400000 47273191.599999 19112570.000000 0.000000 0.000000
0.000000 111000.000000 0.000000 -49783672.000000 10516800.000000 39266872.000000 14613130.000000 0.000000 0.000000
0.000000 222000.000000 0.000000 -48126370.000000 10993830.000000 37132540.000000 11684103.000000 0.000000 0.000000
111000.000000 -222000.000000 0.000000 -73447340.000000 11476001.000000 61971339.000000 -7343758.000000 0.000000 0.000000
111000.000000 -111000.000000 0.000000 -67824624.000000 6208069.000000 61616555.000000 13030550.000000 0.000000 0.000000
111000.000000 0.000000 0.000000 -53048111.000000 8156690.000000 44891421.000000 16181100.000000 0.000000 0.000000
111000.000000 111000.000000 0.000000 -53548432.000000 6474293.400000 47074138.599999 19395603.000000 0.000000 0.000000
111000.000000 222000.000000 0.000000 -49868980.000000 11999960.000000 37869020.000000 15272010.000000 0.000000 0.000000
222000.000000 -222000.000000 0.000000 -82692303.000000 4455240.000000 78237063.000000 -21022010.000000 0.000000 0.000000
222000.000000 -111000.000000 0.000000 -67696911.000000 7180925.000000 60515986.000000 -14294154.000000 0.000000 0.000000
222000.000000 0.000000 0.000000 -58599052.000000 1094036.000000 57505016.000000 4826587.000000 0.000000 0.000000
222000.000000 111000.000000 0.000000 -48208090.000000 3318901.400000 44889188.599999 12947510.000000 0.000000 0.000000
222000.000000 222000.000000 0.000000 -46904143.000000 2044234.000000 44859909.000000 19000964.000000 0.000000 0.000000
-222000.000000 -222000.000000 -50000.000000 -1415266784.560000 -1372019526.560000 -1321469569.560000 1073136.000000 0.000000 0.000000
-222000.000000 -111000.000000 -50000.000000 -1418734200.360000 -1371495194.360000 -1320592906.360000 1061051.000000 0.000000 0.000000
-222000.000000 0.000000 -50000.000000 -1412061641.784000 -1364825655.384000 -1324567148.184000 2885474.000000 0.000000 0.000000
-222000.000000 111000.000000 -50000.000000 -1411517383.520000 -1361851116.520000 -1325615300.520000 4592233.400000 0.000000 0.000000
-222000.000000 222000.000000 -50000.000000 -1419333533.944000 -1366688427.544000 -1321846184.344000 5230918.000000 0.000000 0.000000
-111000.000000 -222000.000000 -50000.000000 -1413031641.880000 -1375953113.880000 -1320963649.880000 13113931.000000 0.000000 0.000000
-111000.000000 -111000.000000 -50000.000000 -1418491361.120000 -1373465197.120000 -1320078432.120000 15048702.000000 0.000000 0.000000
-111000.000000 0.000000 -50000.000000 -1414375828.984000 -1364642551.584000 -1323932357.384000 11741930.000000 0.000000 0.000000
-111000.000000 111000.000000 -50000.000000 -1412312357.840000 -1361671894.840000 -1325431884.840000 8963863.000000 0.000000 0.000000
-111000.000000 222000.000000 -50000.000000 -1414178674.320000 -1358848643.320000 -1325716929.320000 8105111.000000 0.000000 0.000000
0.000000 -222000.000000 -50000.000000 -1439730219.639999 -1370960365.640000 -1314498083.640000 29007560.000000 0.000000 0.000000
0.000000 -111000.000000 -50000.000000 -1423889477.240000 -1360299105.240000 -1322392297.240000 18754142.000000 0.000000 0.000000
0.000000 0.000000 -50000.000000 -1422264822.303999 -1363437527.904000 -1321941387.704000 19112570.000000 0.000000 0.000000
0.000000 111000.000000 -50000.000000 -1415475470.680000 -1355174998.680000 -1326424926.680000 14613130.000000 0.000000 0.000000
0.000000 222000.000000 -50000.000000 -1412879062.600000 -1353758862.600000 -1327620152.600000 11684103.000000 0.000000 0.000000
111000.000000 -222000.000000 -50000.000000 -1449129104.159999 -1364205763.160000 -1313710425.160000 -7343758.000000 0.000000 0.000000
111000.000000 -111000.000000 -50000.000000 -1443350283.199999 -1369317590.200000 -1313909104.200000 13030550.000000 0.000000 0.000000
111000.000000 0.000000 -50000.000000 -1421214711.240000 -1360009910.240000 -1323275179.240000 16181100.000000 0.000000 0.000000
111000.000000 111000.000000 -50000.000000 -1422675427.983999 -1362652702.584000 -1322052857.384000 19395603.000000 0.000000 0.000000
111000.000000 222000.000000 -50000.000000 -1414945723.800000 -1353076783.800000 -1327207723.800000 15272010.000000 0.000000 0.000000
222000.000000 -222000.000000 -50000.000000 -1465530985.720000 -1378383442.720000 -1304601619.720000 -21022010.000000 0.000000 0.000000
222000.000000 -111000.000000 -50000.000000 -1442738319.839999 -1367860483.840000 -1314525422.840000 -14294154.000000 0.000000 0.000000
222000.000000 0.000000 -50000.000000 -1432315634.039999 -1372622546.040000 -1316211566.040000 4826587.000000 0.000000 0.000000
222000.000000 111000.000000 -50000.000000 -1416373707.984000 -1364846716.584000 -1323276429.384000 12947510.000000 0.000000 0.000000
222000.000000 222000.000000 -50000.000000 -1415056877.960000 -1366108500.960000 -1323292825.960000 19000964.000000 0.000000 0.000000
-222000.000000 -222000.000000 -100000.000000 -3158375340.000000 -3115128082.000000 -3064578125.000000 1073136.000000 0.000000 0.000000
-222000.000000 -111000.000000 -100000.000000 -3162719419.000000 -3115480413.000000 -3064578125.000000 1061051.000000 0.000000 0.000000
-222000.000000 0.000000 -100000.000000 -3152072618.600000 -3104836632.200000 -3064578125.000000 2885474.000000 0.000000 0.000000
-222000.000000 111000.000000 -100000.000000 -3150480208.000000 -3100813941.000000 -3064578125.000000 4592233.400000 0.000000 0.000000
-222000.000000 222000.000000 -100000.000000 -3162065474.600000 -3109420368.200000 -3064578125.000000 5230918.000000 0.000000 0.000000
-111000.000000 -222000.000000 -100000.000000 -3156646117.000000 -3119567589.000000 -3064578125.000000 13113931.000000 0.000000 0.000000
-111000.000000 -111000.000000 -100000.000000 -3162991054.000000 -3117964890.000000 -3064578125.000000 15048702.000000 0.000000 0.000000
-111000.000000 0.000000 -100000.000000 -3155021596.600000 -3105288319.200000 -3064578125.000000 11741930.000000 0.000000 0.000000
-111000.000000 111000.000000 -100000.000000 -3151458598.000000 -3100818135.000000 -3064578125.000000 8963863.000000 0.000000 0.000000
-111000.000000 222000.000000 -100000.000000 -3153039870.000000 -3097709839.000000 -3064578125.000000 8105111.000000 0.000000 0.000000
0.000000 -222000.000000 -100000.000000 -3189810261.000000 -3121040407.000000 -3064578125.000000 29007560.000000 0.000000 0.000000
0.000000 -111000.000000 -100000.000000 -3166075305.000000 -3102484933.000000 -3064578125.000000 18754142.000000 0.000000 0.000000
0.000000 0.000000 -100000.000000 -3164901559.599999 -3106074265.200000 -3064578125.000000 19112570.000000 0.000000 0.000000
0.000000 111000.000000 -100000.000000 -3153628669.000000 -3093328197.000000 -3064578125.000000 14613130.000000 0.000000 0.000000
0.000000 222000.000000 -100000.000000 -3149837035.000000 -3090716835.000000 -3064578125.000000 11684103.000000 0.000000 0.000000
111000.000000 -222000.000000 -100000.000000 -3199996804.000000 -3115073463.000000 -3064578125.000000 -7343758.000000 0.000000 0.000000
111000.000000 -111000.000000 -100000.000000 -3194019304.000000 -3119986611.000000 -3064578125.000000 13030550.000000 0.000000 0.000000
111000.000000 0.000000 -100000.000000 -3162517657.000000 -3101312856.000000 -3064578125.000000 16181100.000000 0.000000 0.000000
111000.000000 111000.000000 -100000.000000 -3165200695.599999 -3105177970.200000 -3064578125.000000 19395603.000000 0.000000 0.000000
111000.000000 222000.000000 -100000.000000 -3152316125.000000 -3090447185.000000 -3064578125.000000 15272010.000000 0.000000 0.000000
222000.000000 -222000.000000 -100000.000000 -3225507491.000000 -3138359948.000000 -3064578125.000000 -21022010.000000 0.000000 0.000000
222000.000000 -111000.000000 -100000.000000 -3192791022.000000 -3117913186.000000 -3064578125.000000 -14294154.000000 0.000000 0.000000
222000.000000 0.000000 -100000.000000 -3180682193.000000 -3120989105.000000 -3064578125.000000 4826587.000000 0.000000 0.000000
222000.000000 111000.000000 -100000.000000 -3157675403.600000 -3106148412.200000 -3064578125.000000 12947510.000000 0.000000 0.000000
222000.000000 222000.000000 -100000.000000 -3156342177.000000 -3107393800.000000 -3064578125.000000 19000964.000000 0.000000 0.000000
-222000.000000 -222000.000000 -120000.000000 -3855618762.176000 -3812371504.176000 -3761821547.176000 1073136.000000 0.000000 0.000000
-222000.000000 -111000.000000 -120000.000000 -3860313506.456000 -3813074500.456000 -3762172212.456000 1061051.000000 0.000000 0.000000
-222000.000000 0.000000 -120000.000000 -3848077009.326400 -3800841022.926400 -3760582515.726400 2885474.000000 0.000000 0.000000
-222000.000000 111000.000000 -120000.000000 -3846065337.792000 -3796399070.792000 -3760163254.792000 4592233.400000 0.000000 0.000000
-222000.000000 222000.000000 -120000.000000 -3859158250.862400 -3806513144.462399 -3761670901.262400 5230918.000000 0.000000 0.000000
-111000.000000 -222000.000000 -120000.000000 -3854091907.048000 -3817013379.048000 -3762023915.048000 13113931.000000 0.000000 0.000000
-111000.000000 -111000.000000 -120000.000000 -3860790931.152000 -3815764767.152000 -3762378002.152000 15048702.000000 0.000000 0.000000
-111000.000000 0.000000 -120000.000000 -3851279903.646399 -3801546626.246400 -3760836432.046400 11741930.000000 0.000000 0.000000
-111000.000000 111000.000000 -120000.000000 -3847117094.064000 -3796476631.064000 -3760236621.064000 8963863.000000 0.000000 0.000000
-111000.000000 222000.000000 -120000.000000 -3848584348.272000 -3793254317.272000 -3760122603.272000 8105111.000000 0.000000 0.000000
0.000000 -222000.000000 -120000.000000 -3889842277.544000 -3821072423.544000 -3764610141.544000 29007560.000000 0.000000 0.000000
0.000000 -111000.000000 -120000.000000 -3862949636.104000 -3799359264.104000 -3761452456.104000 18754142.000000 0.000000 0.000000
0.000000 0.000000 -120000.000000 -3861956254.518399 -3803128960.118400 -3761632819.918400 19112570.000000 0.000000 0.000000
0.000000 111000.000000 -120000.000000 -3848889948.328000 -3788589476.328000 -3759839404.328000 14613130.000000 0.000000 0.000000
0.000000 222000.000000 -120000.000000 -3844620223.960000 -3785500023.960000 -3759361313.960000 11684103.000000 0.000000 0.000000
111000.000000 -222000.000000 -120000.000000 -3900343883.935999 -3815420542.936000 -3764925204.936000 -7343758.000000 0.000000 0.000000
111000.000000 -111000.000000 -120000.000000 -3894286912.320000 -3820254219.320000 -3764845733.320000 13030550.000000 0.000000 0.000000
111000.000000 0.000000 -120000.000000 -3859038835.303999 -3797834034.304000 -3761099303.304000 16181100.000000 0.000000 0.000000
111000.000000 111000.000000 -120000.000000 -3862210802.646399 -3802188077.246400 -3761588232.046400 19395603.000000 0.000000 0.000000
111000.000000 222000.000000 -120000.000000 -3847264285.480000 -3785395345.480000 -3759526285.480000 15272010.000000 0.000000 0.000000
222000.000000 -222000.000000 -120000.000000 -3929498093.112000 -3842350550.112000 -3768568727.112000 -21022010.000000 0.000000 0.000000
222000.000000 -111000.000000 -120000.000000 -3892812102.863999 -3817934266.864000 -3764599205.864000 -14294154.000000 0.000000 0.000000
222000.000000 0.000000 -120000.000000 -3880028816.584000 -3820335728.584000 -3763924748.584000 4826587.000000 0.000000 0.000000
222000.000000 111000.000000 -120000.000000 -3854196081.846400 -3802669090.446400 -3761098803.246400 12947510.000000 0.000000 0.000000
222000.000000 222000.000000 -120000.000000 -3852856296.616000 -3803907919.616000 -3761092244.616000 19000964.000000 0.000000 0.000000
-222000.000000 -222000.000000 -160000.000000 -5250105606.528000 -5206858348.528000 -5156308391.528000 1073136.000000 0.000000 0.000000
-222000.000000 -111000.000000 -160000.000000 -5255501681.368000 -5208262675.368000 -5157360387.368000 1061051.000000 0.000000 0.000000
-222000.000000 0.000000 -160000.000000 -5240085790.779200 -5192849804.379200 -5152591297.179200 2885474.000000 0.000000 0.000000
-222000.000000 111000.000000 -160000.000000 -5237235597.375999 -5187569330.375999 -5151333514.375999 4592233.400000 0.000000 0.000000
-222000.000000 222000.000000 -160000.000000 -5253343803.387199 -5200698696.987200 -5155856453.787200 5230918.000000 0.000000 0.000000
-111000.000000 -222000.000000 -160000.000000 -5248983487.144000 -5211904959.144000 -5156915495.144000 13113931.000000 0.000000 0.000000
-111000.000000 -111000.000000 -160000.000000 -5256390685.456000 -5211364521.456000 -5157977756.456000 15048702.000000 0.000000 0.000000
-111000.000000 0.000000 -160000.000000 -5243796517.739200 -5194063240.339200 -5153353046.139200 11741930.000000 0.000000 0.000000
-111000.000000 111000.000000 -160000.000000 -5238434086.192000 -5187793623.192000 -5151553613.192000 8963863.000000 0.000000 0.000000
-111000.000000 222000.000000 -160000.000000 -5239673304.816000 -5184343273.816000 -5151211559.816000 8105111.000000 0.000000 0.000000
0.000000 -222000.000000 -160000.000000 -5289906310.632000 -5221136456.632000 -5164674174.632000 29007560.000000 0.000000 0.000000
0.000000 -111000.000000 -160000.000000 -5256698298.311999 -5193107926.312000 -5155201118.312000 18754142.000000 0.000000 0.000000
0.000000 0.000000 -160000.000000 -5256065644.355200 -5197238349.955200 -5155742209.755199 19112570.000000 0.000000 0.000000
0.000000 111000.000000 -160000.000000 -5239412506.984000 -5179112034.984000 -5150361962.984000 14613130.000000 0.000000 0.000000
0.000000 222000.000000 -160000.000000 -5234186601.880000 -5175066401.880000 -5148927691.880000 11684103.000000 0.000000 0.000000
111000.000000 -222000.000000 -160000.000000 -5301038043.808000 -5216114702.808001 -5165619364.808001 -7343758.000000 0.000000 0.000000
111000.000000 -111000.000000 -160000.000000 -5294822128.959999 -5220789435.960000 -5165380949.960000 13030550.000000 0.000000 0.000000
111000.000000 0.000000 -160000.000000 -5252081191.912000 -5190876390.912000 -5154141659.912000 16181100.000000 0.000000 0.000000
111000.000000 111000.000000 -160000.000000 -5256231016.739200 -5196208291.339200 -5155608446.139200 19395603.000000 0.000000 0.000000
111000.000000 222000.000000 -160000.000000 -5237160606.440001 -5175291666.440001 -5149422606.440001 15272010.000000 0.000000 0.000000
222000.000000 -222000.000000 -160000.000000 -5337479297.335999 -5250331754.336000 -5176549931.336000 -21022010.000000 0.000000 0.000000
222000.000000 -111000.000000 -160000.000000 -5292854264.591999 -5217976428.592000 -5164641367.592000 -14294154.000000 0.000000 0.000000
222000.000000 0.000000 -160000.000000 -5278722063.752000 -5219028975.752000 -5162617995.752000 4826587.000000 0.000000 0.000000
222000.000000 111000.000000 -160000.000000 -5247237438.339200 -5195710446.939199 -5154140159.739200 12947510.000000 0.000000 0.000000
222000.000000 222000.000000 -160000.000000 -5245884535.848000 -5196936158.848000 -5154120483.848000 19000964.000000 0.000000 0.000000
On Nov 6, 2015, at 18:13, Brad Aagaard <baagaard at usgs.gov<mailto:baagaard at usgs.gov>> wrote:
On 11/06/2015 11:06 AM, Chunyu Liu (cliu5) wrote:
Hi all,
Is there any way I can get initial stress data( not that resulting stress) in .VTK file?
I do not find any settings for the initial stress in “ cell_info_fields” and “cell-data_fields”.
For example:
cell_data_fields =[ stress, “stress_initial”]
This is an oversight on our part. It is not currently available. A workaround is to run a static simulation (single time step) without any deformation to get the "initial" stress.
Q2: Interpolation problem
1. I use linear interpolation.
2. several initial stress lines. For example:
// Columns are
25 // (1) x coordinates(m)
26 // (2) y coordinates(m)
27 // (3) z coordinates(m)
28 // (4) stress-xx (Pa)
29 // (5) stress-yy (Pa)
30 // (6) stress-zz (Pa)
31 // (7) stress-xy (Pa)
32 // (8) stress-yz (Pa)
33 // (9) stress-xz (Pa)
34 -222000.000000 -222000.000000 0.000000 -37264580.000000 -131340346.000000 0.000000 7981201.000000 0.000000 0.000000
35 -111000.000000 -222000.000000 0.000000 -36546640.000000 -129611870.000000 0.000000 8891624.000000 0.000000 0.000000
36 0.000000 -222000.000000 0.000000 -36008091.000000 -127917852.000000 0.000000 9781642.000000 0.000000 0.000000
37 111000.000000 -222000.000000 0.000000 -35531603.000000 -126206836.000000 0.000000 10653143.000000 0.000000 0.000000
38 222000.000000 -222000.000000 0.000000 -34957300.000000 -124402730.000000 0.000000 11503441.000000 0.000000 0.000000
39 -222000.000000 -111000.000000 0.000000 -22838144.000000 -116776288.000000 0.000000 11259071.000000 0.000000 0.000000
40 -111000.000000 -111000.000000 0.000000 -22551422.000000 -115173226.000000 0.000000 11932060.000000 0.000000 0.000000
41 0.000000 -111000.000000 0.000000 -22180910.000000 -113548630.000000 0.000000 12538850.000000 0.000000 0.000000
42 111000.000000 -111000.000000 0.000000 -21762040.000000 -111904430.000000 0.000000 13105260.000000 0.000000 0.000000
43 222000.000000 -111000.000000 0.000000 -21621862.000000 -110393807.000000 0.000000 13642334.000000 0.000000 0.000000
44 -222000.000000 0.000000 0.000000 -43329544.000000 -117636548.000000 0.000000 7273745.000000 0.000000 0.000000
45 -111000.000000 0.000000 0.000000 -42661628.000000 -115014364.000000 0.000000 8733876.000000 0.000000 0.000000
46 0.000000 0.000000 0.000000 -40696706.000000 -111989218.000000 0.000000 10180400.000000 0.000000 0.000000
47 111000.000000 0.000000 0.000000 -39202640.000000 -109388170.000000 0.000000 11592650.000000 0.000000 0.000000
48 222000.000000 0.000000 0.000000 -38123380.000000 -107149580.000000 0.000000 12943260.000000 0.000000 0.000000
Your points are not on a line. It looks like they are on a plane, so you should use area interpolation (dim = 2).
Regards,
Brad
_______________________________________________
CIG-SHORT mailing list
CIG-SHORT at geodynamics.org<mailto:CIG-SHORT at geodynamics.org>
http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/cig-short/attachments/20151109/4acd5502/attachment-0001.html>
More information about the CIG-SHORT
mailing list