求助,动力边坡分析,结果始终不对
本帖最后由 rogergoodbye 于 2009-9-1 16:09 编辑我的边坡动力分析结果是坡脚处的位移及加速度最大,但是查询各种已有资料均应该是边坡沿着坡面位移及加速度逐渐增大,到坡顶处最大。 实在查不出我错在了哪里,真心请教各位高手。
附命令流如下:
new
Gen zone brick &
P0 0 0 0 p1 15 0 0 p2 0 5 0 p3 10 0 8 &
p4 15 5 0 p5 10 5 8 p6 15 0 8 p7 15 5 8 &
size 15 5 8
Gen zone brick &
p0 -5 0 -1 p1 15 0 -1 p2 -5 5 -1 p3 -5 0 0 &
p4 15 5 -1 p5 -5 5 0 p6 15 0 0 p7 15 5 0 &
size 2 2 2
model m
prop bulk 6.8e7 shear 3.4e7 coh 45e3 fri 15 ten 8e6
ini dens 2100
;重力添加:
set grav 0 0 -9.8
fix x range x -0.01 0.01
fix x range x 14.99 15.01
fix y range y -0.01 0.01
fix y range y 4.99 5.01
fix z range z -0.01 0.01
solve
ini xd 0 yd 0 zd 0
ini xvel 0 yvel 0zvel 0
;动力:
Conf dyn
set large
Def setup
Freq=2.0
Omega=1.0*pi*freq
Old_time=clock
End
Setup
Def wave
Wave=sin(omega*dytime)
End
free x y z ran z -0.1 0.1
apply nquiet squiet dquiet ran z -0.1 0.1
Apply xacc=1 hist wave range z=-.1 .1
Hist dytime
hist unbal
Def tim
Tim=0.01*(clock-old_time)
End
Set dyn multi on
Solve age 2.0
PS:由于只想找地震动对边坡的破坏的反应规律,故未考虑设置自由场边界 怎么没人看啊,急盼高手回复 在静态边界上不能直接加加速度(你加的是加速度) 只能加应力
可以先将加速度时程转化成速度时程再转化成应力时程输入 apply nquiet squiet dquiet ran z -0.1 0.1
Apply xacc=1 hist wave range z=-.1 .1
静态边界上输入动荷载只能输入应力时程
可以将加速度、速度时程通过公式转换成应力时程施加在静力边界上。 楼上是正解,看看手册 apply nquiet squiet dquiet ran z -0.1 0.1
xacc=1 hist wave range z=-.1 .1
静态边界上输入动荷载只能输入应力时程
可以将加速度、速度时程通过公式转换成应力时程施加在静力边界上。 我采用静态边界+应力波输入,得出的加速度分布没什么特别好的规律,反而是坡脚处的加速度最大 模量比较大的刚性地基,可以直接在底部施加动力荷载,无需设置静态边界。 如果设置静态粘滞边界,应力要乘以2的,如果直接施加速度、加速度是不需要乘2的。 跟以上人说的没关吧,我也是设置的静态边界,施加应力时程,也是坡脚处加速度峰值最大 这数值模拟就是玄学
页:
[1]