边界条件的疑惑，请求解答。

zy7788250

手册上的动力分析的例子，但是对于边界条件的设置，很疑惑，请达人指教。
;
; --- Free-field test ---
;
config dyn
;
def iniwave
  per = 0.01
end
def wave
  wave = 0.5 * (1.0 - cos (2*pi*dytime/per))
end
iniwave
;
gen zone brick size 6 3 2
gen zone brick size 2 3 2 p0 0 0 2
gen zone brick size 2 3 2 p0 4 0 2
gen zone wedge size 1 3 2 p0 2 0 2
gen zone wedge size 1 3 2 p0 4 3 2 p1 3 3 2 p2 4 0 2 p3 4 3 4  &
                                   p4 3 0 2 p5 4 0 4
;
model elastic
prop  bulk 66667 shear 40000
ini  dens 0.0025
;
; --- static state ---
;
set grav 0 0 -10
;
fix x  range x -0.01 0.01
fix x  range x  5.99 6.01
fix y  range y -0.01 0.01
fix y  range y  2.99 3.01
fix z  range z -0.1 0.1
;
set dyn off
;
hist  unbal
hist  dytime
hist  gp zdis 2 1 0
hist  gp zdis 2 1 2.0
hist  gp zdis 2 1 5.0
;
solve
;
save ff0.sav
;
; --- dynamic run A: applied velocity at base ---
;
set dyn on
;
; --- set dyn b.c. at base ---
free x y z  ran z -0.1 0.1  ？？？静力状态下 range z -0.1 0.1 只是fix z，为什么进行动力计算的时候free x y z？free究竟是什么意思？
apply nquiet squiet dquiet ran z -0.1 0.1
apply dstress 1.0 hist wave ran z -0.1 0.1
;
; --- FF ---
apply ff
;
group ff_corner
group ff_side   ran x 0 6
group ff_side   ran y 0 3
group main_grid ran x 0 6 y 0 3
;
;pause
set dyn time = 0
hist reset
;
hist  unbal
hist  dytime
;
; main grid
hist  gp xvel 2 1 0
hist  gp xvel 2 1 5.0
;
; corner FF
hist  gp xvel -1 -1 0
hist  gp xvel -1 -1 5.0
;
; side FF parallel to y
hist  gp xvel -1 0 0
hist  gp xvel -1 0 5.0
;
; side FF parallel to x
hist  gp xvel 2 -1 0
hist  gp xvel 2 -1 5.0
;
solve age 0.015
save ff1.sav
ret

zy7788250

自己来解答一下
分别试算了一下几种情况：
1. free z  ran z -0.1 0.1；
2. free x y  ran z -0.1 0.1；
3. 去掉free。
这三种情况的计算结果是一样的。
解释如下：
手册上关于quiet boundary有如下描述：Dynamic analysis starts from some in-situ condition. If a fixed boundary is used while generating the static stress state, this boundary condition can be replaced by quiet boundaries; the boundary gridpoints will be freed, and the boundary reaction forces will be automatically calculated and maintained throughout the dynamic loading phase.
我是这样理解，动力计算中的quiet boundary会将静力计算中的固定边界取代，固定边界将会释放。因此，free和去掉free得出的结果是一样的。但是，如果施加，速度和加速度荷载的话，此时不能用queit boundary，就必须free固定边界。
不知道理解的对不对，请达人指点。

[ 本帖最后由 zy7788250 于 2007-8-6 10:52 编辑 ]

zy7788250

我理解的对吗？
恳求达人指点。

liucheng83

原帖由 zy7788250 于 2007-8-6 10:51 发表
自己来解答一下
分别试算了一下几种情况：
1. free z  ran z -0.1 0.1；
2. free x y  ran z -0.1 0.1；
3. 去掉free。
这三种情况的计算结果是一样的。
解释如下：
手册上关于quiet boundary有如下描述 ...

跟我的理解一样，谢谢！

菜鸟鸟

有人能解释一下在FLAC3D中怎样加粘弹性边界吗？

婴儿乐

郁闷，仿真币又不够了。