自重应力计算的问题

zhaolei1987

gen zone brick p0 0 -24 0 p1 add 25 0 0 p2 0 20 0 p3 add 0 0 40 size 20 20 30按这个简单的模型建立一个土体，然后做自重应力计算，得到的均匀的沉降云图。

问题：

gen zone radcyl p0 5 40 4.5 p1 add 5 0 0 p2 add 0 20 0 p3 add 0 0 6.5 dim 2.5 2.5 2.5 2.5 size 1 20 10 10 rat 1 1 1 1.1 group 围岩
gen zone cshell p0 5 40 4.5 p1 add 2.5 0 0 p2 add 0 20 0 p3 add 0 0 2.5 dim 2 2 2 2 size 1 20 10 1 group liner1 fill group 水管

gen zone reflect origin 5 0 4.5 normal 0 0 -1
gen zone reflect origin 5 0 0 normal -1 0 0

gen zone radtun p0 0 40 11 p1 add 0 0 -27 p2 add 0 20 0 p3 add 25 0 0 &
p4 add 0 20 -27 p5 add 25 20 0 p6 add 25 0 -27 p7 add 25 20 -27 p8 add 0 0 -13 &
p9 add 10 0 0 p10 add 0 20 -13 p11 add 10 20 0 p12 add 10 0 -13 p13 add 10 20 -13 &
dim 11 20 10 16 size 10 20 10 15 rat 1 1 1 1 group 围岩
这个模型的长宽高和上面的长方体一样，只不过是模拟的是在土体中要挖一个隧洞，图面在附件里面。


两种建模采取相同的边界和材料属性，第二种建模的自重应力计算，沉降云图就是错误的，不仅不是均匀的，而且还有些地方出现了隆起，请只点。谢谢！


我把两种计算的代码和计算后的自重应力云图都发上去大家看看。
谢谢指点！