不平衡

luntanren_lpc

大神可以看下我的参数有问题吗，建立个模型，50万网格，不能平衡啊，最大不平衡力起步50步直接变零，然后不断增大，这是哪里的问题啊
Prop density 2500 bulk 6.74e9   shear 4.52e9  coh 2.67e6 fric 31 ...
tension 1.2e6 range  group 'yanshi1'
Prop density 2500 bulk 6.74e9   shear 4.52e9  coh 2.67e6 fric 31 ...
tension 1.2e6 range  group 'yanshi2'
Prop density 1400 bulk 4.33e9  shear 1.34e9   coh 1.07e6 fric 18 ...
tension 6e5   range x -8 8 y 2 15 z  -4 4
Prop density 2550 bulk 18.37e9  shear 13.96e9  coh 2.8e6 fric 52 ...
tension 1.78e6 range  group  'jiati'
set grav 0 0 -10
;约束与边界,埋深500，侧压系数0.7
fix x y z range z -9.1 -8.9
apply szz -12.5e6 range z 8.9 9.1
apply sxx 8.75e6 range x -8.1 -7.9
apply sxx -8.75e6 range x 7.9 8.1
fix x y   range  y -0.1 0.1
fix x y   range  y 14.9 15.1