叶harry 发表于 2008-7-9 22:56:34

我做的一个简单边坡动力,位移总是算不对,帮忙看下

我的命令流
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





我感觉边界和阻尼都是没问题的,但是位移总是算不对,加速度也做了基线校正,请斑竹帮我看看到底那里出问题了

lpfvsfl 发表于 2009-6-1 20:55:01

可以指导下我做动力分析部?我也是做边坡!QQ:312169218
页: [1]
查看完整版本: 我做的一个简单边坡动力,位移总是算不对,帮忙看下