求助大侠高手帮忙看下我写的这个命令 [复制链接]

网格都建好了，但是初始应力求解不了，很苦恼，帮忙分析下原因
new
;矿体
gen zone brick p0 0 0 0 p1 80 0 0 p2 0 30 0 p3 0 45 45 p4 80 30 0 p5 0 75 45 p6 80 45 45 p7 80 75 45 size 20 20 2  0 group KT
gen zone brick p0 80 0 0 p1 320 0 0 p2 80 30 0 p3 80 45 45 p4 320 30 0 p5 80 75 45 p6 320 45 45 p7 320 75 45 size 18 20 20 ratio 1.1 1 1 group KT
gen zone brick p0 -240 0 0 p1 0 0 0 p2 -240 30 0 p3 -240 45 45 p4 0 30 0 p5 -240 75 45 p6 0 45 45 p7 0 75 45 size 18 20 20 ratio 0.91 1 1 group KT
gen zone brick p0 0 45 45 p1 80 45 45 p2 0 75 45 p3 0 180 180 p4 80 75 45 p5 0 210 180 p6 80 180 180 p7 80 210 180 size 20 20 18 ratio 1 1 1.1 group KT
gen zone brick p0 0 -135 -135 p1 80 -135 -135 p2 0 -105 -135 p3 0 0 0 p4 80 -105 -135 p5 0 30 0 p6 80 0 0 p7 80 30 0 size 20 20 18 ratio 1 1 0.91 group KT
gen zone brick p0 -240 -135 -135 p1 0 -135 -135 p2 -240 -105 -135 p3 -240 0 0 size 18 20 18 ratio 0.91 1 0.91 group KT
gen zone brick p0 80 -135 -135 p1 320 -135 -135 p2 80 -105 -135 p3 80 0 0 size 18 20 18 ratio 1.1 1 0.91 group KT
gen zone brick p0 -240 45 45 p1 0 45 45 p2 -240 75 45 p3 -240 180 180 size 18 20 18 ratio 0.91 1 1.1 group KT
gen zone brick p0 80 45 45 p1 320 45 45 p2 80 75 45 p3 80 180 180 size 18 20 18 ratio 1.1 1 1.1 group KT
;下盘
gen zone brick p0 0 -105 -135 p1 80 -105 -135 p2 0 210 -135 p3 0 30 0 p4 80 210 -135 p5 0 210 0 p6 80 30 0 p7 80 210 0 size 20 18 18 ratio 1 1.1 0.91 &
group XP
gen zone brick p0 0 30 0 p1 80 30 0 p2 0 210 0 p3 0 75 45 p4 80 210 0 p5 0 210 45 p6 80 75 45 p7 80 210 45 size 20 18 20 ratio 1 1.1 1 group XP
gen zone w p0 0 210 45 p1 0 75 45 p2 80 210 45 p3 0 210 180 size 18 20 18 ratio 0.91 1 1.1 group XP
gen zone brick p0 -240 -105 -135 p1 0 -105 -135 p2 -240 210 -135 p3 -240 30 0 p4 0 210 -135 p5 -240 210 0 p6 0 30 0 p7 0 210 0 size 18 18 18 ratio 0.91 1.1 0.91 &
group XP
gen zone brick p0 80 -105 -135 p1 320 -105 -135 p2 80 210 -135 p3 80 30 0 p4 320 210 -135 p5 80 210 0 p6 320 30 0 p7 320 210 0 size 18 18 18 ratio 1.1 1.1 0.91 &
group XP
gen zone brick p0 -240 30 0 p1 0 30 0 p2 -240 210 0 p3 -240 75 45 p4 0 210 0 p5 -240 210 45 p6 0 75 45 p7 0 210 45 size 18 18 20 ratio 0.91 1.1 1 &
group XP
gen zone brick p0 80 30 0 p1 320 30 0 p2 80 210 0 p3 320 75 45 p4 320 210 0 p5 80 210 45 p6 320 75 45 p7 320 210 45 size 18 18 20 ratio 1.1 1.1 1 &
group XP
gen zone w p0 -240 210 45 p1 -240 75 45 p2 0 210 45 p3 -240 210 180 p4 0 75 45 p5 0 210 180 size 18 18 18 ratio 0.91 0.91 1.1 group XP
gen zone w p0 80 210 45 p1 80 75 45 p2 320 210 45 p3 80 210 180 p4 320 75 45 p5 320 210 180 size 18 18 18 ratio 0.91 1.1 1.1 group XP
;上盘
gen zone brick p0 0 -135 0 p1 80 -135 0 p2 0 0 0 p3 0 -135 45 p4 80 0 0 p5 0 45 45 p6 80 -135 45 p7 80 45 45 size 20 18 20 ratio 1 0.91 1 group SP
gen zone w p0 0 -135 0 p1 0 0 0 p2 80 -135 0 p3 0 -135 -135 size 18 20 18 ratio 0.91 1 1.1 group SP
gen zone brick p0 -240 -135 0 p1 0 -135 0 p2 -240 0 0 p3 -240 -135 45 p4 0 0 0 p5 -240 45 45 p6 0 -135 45 p7 0 45 45 size 18 18 20 ratio 0.91 0.91 1 group SP
gen zone brick p0 80 -135 0 p1 320 -135 0 p2 80 0 0 p3 80 -135 45 p4 320 0 0 p5 80 45 45 p6 320 -135 45 p7 320 45 45 size 18 18 20 ratio 1.1 0.91 1 group SP
gen zone w p0 -240 -135 0 p1 -240 0 0 p2 0 -135 0 p3 -240 -135 -135 size 18 18 18 ratio 0.91 0.91 1.1 group SP
gen zone w p0 80 -135 0 p1 80 0 0 p2 320 -135 0 p3 80 -135 -135 size 18 18 18 ratio 0.91 1.1 1.1 group SP
gen zone brick p0 0 -135 45 p1 80 -135 45 p2 0 45 45 p3 0 -135 180 p4 80 45 45 p5 0 180 180 p6 80 -135 180 p7 80 180 180 size 20 18 18 ratio 1 0.91 1.1 group SP
gen zone brick p0 80 -135 45 p1 320 -135 45 p2 80 45 45 p3 80 -135 180 p4 320 45 45 p5 80 180 180 p6 320 -135 180 p7 320 180 180 size 18 18 18 ratio 1.1 0.91 1.1 group SP
gen zone brick p0 -240 -135 45 p1 0 -135 45 p2 -240 45 45 p3 -240 -135 180 p4 0 45 45 p5 -240 180 180 p6 0 -135 180 p7 0 180 180 size 18 18 18 ratio 0.91 0.91 1.1 group SP
plot add surface yellow
plot add axes black
plot show
;
model mohr
prop bulk 2.24e11 shear 1.68e11 coh 5.72e6 fri 30.6 ten 3.18e6 range group SP
prop bulk 2.42e11 shear 1.89e11 coh 6.43e6 fri 32.6 ten 3.72e6 range group KT
prop bulk 3.29e11 shear 2.07e11 coh 10.7e6 fri 36.9 ten 4.31e6 range group XP
ini dens 2706 range group SP
ini dens 2710 range group KT
ini dens 2635 range group XP
;应力条件
apply sxx -16.04e6 grad 0 0 1.95e4 ran x -240.1 -239.9
apply sxx -16.04e6 grad 0 0 1.95e4 ran x 319.9 320.1
apply syy -30.16e6 grad 0 0 3.62e4 ran y -135.1 -134.9
apply syy -30.16e6 grad 0 0 3.62e4 ran y 209.9 210.1
apply szz -16.59e6 grad 0 0 1.86e4 ran z -135 180
set grav 0 0 -9.81
solve elas
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
我查看了好多次  都找不出原因

我对巷道不是太了解，由于你的模型是三维的，我运行了一下，很慢。如果程序出了问题，我建议你把三维的先改成二维的在去查看。
程序中我有两点看法，不知道对不对了，你可以参考下。一：我看你后面没有定义FIX，这个边界条件在隧道里头都是要定义的，不然应该是没法计算的，你看你这个要不要加上；二：你定义的是摩尔库伦的模型，但是在最后是运算的是弹性，不知道这个对不对；三：你的程序应该还没有完全，为什么运算了过后再清零没有在做初支或者什么加固。第一个清零是初始应力产生的位移的清零，这个在隧道里头常用，第二个清零不太清楚，但感觉程序还没有运算完全。就这些了

恩  谢谢你。
我本想建好模型，有矿体 下盘 上盘，按照德是分阶段求初始应力，所以才给他设了摩尔库伦模型，最后运用弹性求解，就是赋予模型摩尔库伦模型的时候，出现了单元体有不合法的规格的警告语句，是否建模的时候出现了问题  才导致求解不了

刚看过了，就是你的网格有问题。
因为你是分块建立的网格，各个分块之间的网格节点没有对齐，导致模型内部存在不耦合面（自由面）。
你可以在建模时每执行一次gen 句命令，就显示网格看下，看看到底那部分的坐标没有对齐。
以下为错误提示，以及显示模型表面网格图，从图2 中明显可见你的模型内部的不耦合面

再者meizu29兄说的不错，没有加边界（fix）约束，等等吧

恩   我现在一步一步查看下   非常谢谢你！

网格建错了  有个坐标出错了   现在可以求解了  谢谢你们

是我的有个坐标弄错了  所以才不行   谢谢你们帮忙

我建的矿体是45度角的，我想定义个空模型，采完矿体的一部分，比如矿体采高  x  和 z 方向的好定  可是y方向我不知道怎么去定义

把矿体定义为一个group，这样在采矿时，只用z和group来控制就可以了，
命令格式如下:
model null range group 矿体 z ? ?
;这个range后面控制的采矿部分为“group 矿体”和“z ? ?”两个范围相交的部分，即采掉的部分为“group 矿体”中的“z ? ?”部分

恩  这下明白了  头一次做  遇到了不少问题   谢谢哈

还有个疑问  从昨天晚上我开始求解运行了30000多步  还没平衡  结果SZZ  应力云图竟然都是0，什么原因？  是否把网格建成宽一点  把size  20  20  20 改成size 10 10  10 是否求解快点？