我做的一个简单边坡动力,位移总是算不对,帮忙看下
我的命令流conf dyn
gen zone brick p0 -50 0 -50 p1 0 0 -50 p2 -50 1 -50 p3 -50 0 0 size 10 1 10
gen zone brick p0 0 0 -50 p1 20 0 -50 p2 0 1 -50 p3 0 0 0 size 4 1 10
gen zone brick p0 20 0 -50 p1 100 0 -50 p2 20 1 -50 p3 20 0 0 size 16 1 10
gen zone brick p0 20 0 0 p1 100 0 0 p2 20 1 0 p3 20 0 20 size 16 1 4
gen zone wedge p0 20 0 0 p1 20 0 20 p2 20 1 0 p3 0 0 0 size 4 1 4
range name bottom z=-50.1 -49.9
model mohr
prop shea 1e10 bulk 1.2e10 fric=30 cohes=1e6 tens=7e5
ini dens 2000
fix x range x=-50.1 -49.1 z=-50.1 0.1
fix x range x=99.9,100.1 z=-50.1 20.1
fix y
fix x y z range bottom
set dyn=off grav=0,0,-10 small
his unbal
solve
ini xv=0 yv=0 zv=0 xdis=0 ydis=0 zdis=0
save slope1.sav
rest slope1.sav
set dyn=on large
set dynamic dt 1e-4
free x y z ran bottom
table 1 read 2eicentrol.dat
apply xacc 1 his table 1 yacc=0 zacc=0 range bottom
apply ff
hist gp xvel 0,0,0
hist gp xvel 5,0,5
hist gp xvel 10,0,10
hist gp xvel 15,0,15
hist gp xvel 20,0,20
hist gp xacc 0,0,0
hist gp xacc 5,0,5
hist gp xacc 10,0,10
hist gp xacc 15,0,15
hist gp xacc 20,0,20
hist gp xdisp 0,0,0
hist gp xdisp 5,0,5
hist gp xdisp 10,0,10
hist gp xdisp 15,0,15
hist gp xdisp 20,0,20
hist dytime
hist gp xacc 0,0,-50
hist gp xdisp 100,0,-50
hist gp xdisp -50,0,-50
set dyn damp local 0.15 ;(approx 5% damping)
solve age 55
save slope2.sav
我感觉边界和阻尼都是没问题的,但是位移总是算不对,加速度也做了基线校正,请斑竹帮我看看到底那里出问题了 可以指导下我做动力分析部?我也是做边坡!QQ:312169218
页:
[1]