自由场边界条件的实例

zxf613

new
;第一步：静力计算阶段
config dyn
set dyn off
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
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
hist unbal
solve
save 11-3_1.sav
;第二步：动力计算阶段
set dyn on
def iniwave
  per = 0.01   
end
iniwave
def wave
  wave = 0.5 * (1.0 - cos (2*pi*dytime/per))
end
free x y z  ran z -0.1 0.1                ;去掉模型底部原有的静力条件
apply nquiet squiet dquiet ran z -0.1 0.1  ;静态边界条件
apply dstress 1.0 hist wave ran z -0.1 0.1 ;加动力荷载
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
set dyn time = 0                         ;设置动力计算从0s开始
hist reset                                ;清空已有的历史信息
hist  unbal
hist  dytime
; 主体网格
hist  gp xvel 2 1 0
hist  gp xvel 2 1 5.0
; 柱体网格
hist  gp xvel -1 -1 0
hist  gp xvel -1 -1 5.0
; 平行于y方向的二维自由场网格
hist  gp xvel -1 0 0
hist  gp xvel -1 0 5.0
; 平行于x方向的二维自由场网格
hist  gp xvel 2 -1 0
hist  gp xvel 2 -1 5.0
solve age 0.015
save 11-3_2.sav
这是陈育民老师书上的一个例题，我有几个问题想问问各位。
1.静力状态下在z=0处只固定z方向，为什么后面动力计算free x y z  ran z -0.1 0.1呢？
2.最后面主体网格是怎么形成的？
hist  gp xvel 2 1 5.0是什么意思 啊？

syx2006

释放约束条件，让此边界吸收能量。

zxf613

请问楼上：动力计算怎么没有设置阻尼啊？例如上面这个程序

baby20023402

就是，还是为什么他的密度设置这么小，原因是什么呢？？