刚学FLAC3D，就碰到了问题，请教一下论坛的高手们，希望能解决一下！
自己建立了一个模拟地下岩体应力平衡的模型，是一个200*100*50的长方条模型，模型由五层岩性组成，由中间向两头呈对称性分布，
由上至下分别为砂岩层、软弱带过渡层（砂岩夹糜棱岩）、软弱带层（断层带物质）、软弱带过渡层、砂岩层，然后施加应力使其平衡，计算平衡后得到的SZZ向应力结果图如下所示。我的问题是：1、为什么计算平衡后的模型，SSZ应力不呈对称性分布呢？2、模型右下角（蓝色块）为什么会出现应力集中区？3、右下部分（蓝色块左侧）的黄色应力区域为什么会穿越到模型底部去？4、这个计算平衡后得到的SZZ应力分布结果图合理吗？
模型的具体代码如下所示：
gen zone brick size 20 10 1 p0 0,0,25 p1 200,0,25 p2 0,100,25 p3 0,0,25.5
gen zone brick size 20 10 3 p0 0,0,25.5 p1 200,0,25.5 p2 0,100,25.5 p3 0,0,28.5
gen zone brick size 20 10 8 p0 0,0,28.5 p1 200,0,28.5 p2 0,100,28.5 p3
0,0,50 ratio 1 1 1.2
gen zone reflect dip 0 dd 90 ori 0 0 25
;先用宏对象分别定义三种不同的材料，然后通过range调用宏来给模型的不同区域赋上不同的材料参数--------------
model mohr
macro sandstone 'bulk 99.6e8 shear 93.5e8 fric 44 coh 5e6 ten 4e6'
macro weakstone 'bulk 3.64e8 shear 2.07e8 fric 23 coh 0.3e6 ten 0.1e6'
macro flow 'bulk 0.89e8 shear 0.8e8 fric 18 coh 0.013e6 ten 1e6'
prop sandstone range z 0 21.5
prop weakstone range z 21.5 24.5
prop flow range z 24.5 25.5
prop weakstone range z 25.5 28.5
prop sandstone range z 28.5 50
initial density 2800 range z 0 21.5
initial density 2600 range z 21.5 24.5
initial density 2000 range z 24.5 25.5
initial density 2600 range z 25.5 28.5
initial density 2600 range z 28.5 50
set grav 0 0 -10
fix x range x -0.1 0.1
fix x range x 199.9 200.1
fix y range y -0.1 0.1
fix y range y 99.9 100.1
fix z range z -0.1 0.1
apply nstress -26.1e6 range z 49.9 50.1
;用26.1MPa法向应力模拟岩层上覆的自重应力
ini szz -27.5e6 grad 0,0,5.22e5
ini sxx -22.5e6 grad 0,0,2.22e5
ini syy -22.5e6 grad 0,0,2.22e5
set grav 0 0 -10
hist unbal
hist gp zdisp 100 25 25
solve
plot create szz_contour
plot add cont szz out on shade on
plot add axes
plot show
