为什么这么简单的一个接触面的模型都平衡不了啊。
gen zone brick p0 0 0 0 p1 20 0 0 p2 0 20 0 p3 0 0 20 size 10 10 10 ratio 1 1 1
m m
property bulk 10e6 shear 3.5e6 friction 10 cohesion 10e3 ten 1e1
fix x range x -0.1 0.1
fix z range z -0.1 0.1
fix y range y -0.1 0.1
fix x range x 19.9 20.1
fix y range y 19.9 20.1
set gravity 0 0 -10
initial density 1800 range z 0 20
group inner range x 8 10 y 0 20 z 10 20
group outer range group inner not
gen separate inner
gen merge 1e-4 range x 8 10 y 0 20 z 9.9 10.1
int 1 face range group inner x 7.9 8.1 y 0 20 z 10 20
int 2 face range group inner x 9.9 10.1 y 0 20 z 10 20
interface 1 prop kn 1e8 ks 1e8 fric 20 cohesion 2e3
interface 2 prop kn 1e8 ks 1e8 fric 20 cohesion 2e3
solve