请教大家一个问题,最近做了一个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
def aa
p_gp=gp_near(a1,a2,10.0)
pz=gp_id(p_gp)
end
set a1=10.0
set a2=10.0
aa
;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
现在自己十分困惑,希望群里各位高手不吝赐教啊,小弟感激不尽。