关于FLAC3D边界改变后计算结果的困惑

daihaifeng

请教大家一个问题，最近做了一个FLAC3D的例子，有一个现象自己无法理解。例子如下：
new
gen zone brick size 20 20 10
model mohr
prop bulk 3.50e7 shear 1.50e7 coh 1e10 fri 24.73 ten 1e10
ini dens 1700
fix z ran z -0.1 0.1
fix x ran x -0.1 0.1
fix x ran x 19.9 20.1
fix y ran y -0.1 0.1
fix y ran y 19.9 20.1
ini dens 1000
set grav 0 0 -10
solve
prop bulk 3.50e7 shear 1.50e7 coh 1.94e4 fri 24.73 ten 0
solve
save k_1.sav

rest k_1.sav
ini xd 0 yd 0 zd 0
ini xv 0 yv 0 zv 0
conf dy
set dy on
def setup
  omega=2.0*pi*freq/6
end
set freq=100
setup
def wave
   if dytime<1.0
   wave=(1-cos(pi*dytime/3))*sin(omega*dytime)/2
   else
   wave=0.0
   endif
end
apply dq sq nq range x -0.1 0.1
apply dq sq nq range x 19.9 20.1
apply dq sq nq range y -0.1 0.1
apply dq sq nq range x 19.9 20.1
apply dq sq nq range z -0.1 0.1
set dyn dt 2e-3
set large
;apply nstress -1e3 hist wave range x 9.5 10.5 y 9.5 10.5 z 9.9 10.1
hist gp zd 10 10 10
hist gp zv 5 6 5
hist gp zv 5 5 10
hist gp zv 5 5 0
hist gp zd 10 10 0
hist dytime
set dyn damp local 0.314
plot create hhh
plot add hist 1  5  v 6
plot set plane ori 10 10 5 norm 0 1 0
plot add con zd plane
;plot add hist 2 v 6
plot show
solve age 2.0

两端程序运行完后，发现顶部（10，10 ，20）和底部（10，10，0）两点的速度都不为零，最后整个模型将会直线向下匀速运动，应该是开始时受力不平衡，从而有个加速度对模型加速，一段时间后受力恢复平衡整个模型以加速后速度做匀速运动。但是我找不出不平衡的原因在哪里。
后面我又用一个弹性模型做了相似的例子，只不过没有固结沉降这一部分，也就是没有加set grav 0 0 -10，此时模型底部不会运动，或者说相对于其他点位移数值极小。
new
gen zone brick size 20 20 20
model elas
prop bulk 1067 shear 400
ini dens 1200
conf dy
def wave
  if dytime<4e-3
     wave=1.0
  else
     wave=0.0
  endif
end
def aa
    p_gp=gp_near(a1,a2,20.0)
    pz=gp_id(p_gp)
end
set a1=10.0
set a2=10.0
aa
apply dq sq nq range x -0.1 0.1
apply dq sq nq range x 19.9 20.1
apply dq sq nq range y -0.1 0.1
apply dq sq nq range x 19.9 20.1
apply nq sq dq range z -0.1 0.1
apply xf 0.0 yf 0.0 zf -1e3 hist wave range id pz
set dyn damp local 0.314
set dyn dt 2e-3
set large
hist gp zvel 10 10 20
hist gp zdis 10 10 20
hist gp zdis 10 10 0
hist dytime
hist wave
plot create hhh
plot add hist 2 3 v 4
plot set plane ori 10 10 20 norm 0 1 0
plot add con zd plane
;plot add hist 1 v 3
plot show
solve age 2
第一个程序表明即使不加任何荷载，模型也会出现整体下移的现象。现在自己十分困惑，希望群里各位高手不吝赐教啊，小弟感激不尽。