对动力时步的考虑

wangguobo4220

我的算例：
config dynamic
set log on
gen zone brick size 30 25 10 p0 0 0 0 p1 3 0 0 p2 0 2.5 0 p3 0 0 1 0 model e
prop bulk 1e8 shear 4.5e8
ini dens 1000
;固定左、右及下边界
fix x range x -0.1 0.1
fix x range x 29.9 30.1
fix y range y -0.1 0.1
fix y range y 24.9 25.1
fix z range z -0.1 0.1
set dyna off
set grav 0 0 -10
hist unbal
hist gp zdis 15 12 5
cy 500
save static.sav
;------------------------以下开始动力计算-----------------------
;动力计算前，应该先清零
ini xdisp=0 ydisp=0 zdisp=0
ini xvel=0 yvel=0 zvel=0
;------------定义输入波――――――――
def wave
wave=ampl*sin(2.0*pi*freq*dytime)
end
set freq=25.0 ampl=0.5
wave
set dyn on
；set dynamic dt 1.25e-4   ；设置动力时步
;边界设置自由场边界
free x y z range z -0.1 0.1    ；释放底面，地震动输入
app xacc=10.0 his wave range z -0.1 0.1 ；地震波输入
apply ff                  ；加自由场边界
；set dyn damp rayleigh 1 25 ；先试算，未考虑阻尼
;选取的几个比较典型的节点，来进行记录。
hist reset
hist unbal
hist dytime
;-------x direction acceleration-------
his gp xacc 0 12 0      ；模型底部         
his gp xacc 0 12 5      ；模型中部
his gp xacc 0 12 10
solve age 5
save dpdynamic2.sav

手册里在动力分析之后有个例子来说明动力分析的过程，其中对动力时步的考虑有两个条件见手册（这里贴不上来，抱歉），一个是临界动力时步，一个是要保证波能够在介质中传播。我模型很规则，对第一个条件，我的计算模型有：l＝0.1m，A＝0.01m2，时步为1.25e-4，而系统默认值为5.556e－5，还是比我计算值小，就用了系统默认值计算，对于第二个条件，我有不同的理解，这里贴不上，烦请各位看附件！就不多写了

lakewater

静力边界条件错误，修改后请继续分析。

freddie_1999

lakewater兄能不能说说这里的静止边界哪里错了?希望能得到你的指导,谢谢

xjiang01

1、
gen zone brick size 30 25 10 p0 0 0 0 p1 3 0 0 p2 0 2.5 0 p3 0 0 1 0 model e
其中p3坐标输入错误，应为p3 0 0 1；
2、
fix x range x -0.1 0.1
fix x range x 29.9 30.1
fix y range y -0.1 0.1
fix y range y 24.9 25.1
fix z range z -0.1 0.1
静力边界错误，模型尺寸为３＊２.5＊１，而不是３０＊２５＊１０。

freddie_1999

哦,搞错了,还以为是静止边界

liucheng83

好像计算compression velocity和dt有错吧
还有请教一个问题：
  如果在底边施加了加速度波，就不用加quiet boudary，那么好像真的会回弹波，所以在上面会出现很混乱的加速度波形

zxf613

好例子……

spwxxn

是的，不要把网格数当作边长了呀，4楼分析的太对了，顶·

qiufenghust

如果要把静力边界换成粘滞边界，命令是怎样的？