shell43 做的箱涵分析，请拍砖

xugangseu

主要目的是得到截面剪力和弯距
1，我想考虑纵向的影响，并且想真实的模拟出车轮的作用，用的是shell43和combin14建的模型，但是对于弹簧单元的约速不知如何约速（是只约速底边还是四边都约速，如果四边都约束，弹簧的弹性系数如何取），
2，两边的土压力不是对称的，两边都是梯形施加，但是一边大一边小（考虑车辆的影响），我施加的是面荷载，单位是Kpa，弹簧我取的是地基反应模量，单位是N/M3,不知有没有什么问题，我看过用beam3建的隧道平面模型，好像是这么取得。
现附上我得命令流，请多多指教，多谢！
QQ：461068444  
fini
/cle
/nopr
keyw,pr_struc,1
/prep7
et,1,shell43
et,2,combin14
mp,ex,1,2.8e10
mp,nuxy,1,0.2
mp,dens,1,2500
r,1,0.5
r,2,30e6   !是否需要换算，用的是地基反应模量，单位N/m3
!!
k,1
k,2,5
k,3,5,5
k,4,,5
lstr,1,2
lstr,2,3
lstr,2,3
lstr,3,4
lstr,1,4
a,1,2,3,4
VEXT,1,,,,,10
vdele,1
adele,1,2,1
!!
lsel,s,line,,9,12,1
lesize,all,0.5
lsel,s,line,,1,8,1
lesize,all,0.5
allsel,all
type,1
mat,1
real,1
amesh,all
/pnum,elem,1
!!
nsel,s,loc,y,0
ngen,2,1000,all,,,,-1
nplot
allsel,all
/rep
!!
type,2
mat,1
real,2
*do,i,1,231,1
e,i,i+1000
*enddo
eplot
!!
/sol
nsel,s,loc,y,-1 !约束弹簧底端的所有自由度
d,all,all
allsel,all
nsel,s,loc,y,0  !我约束了箱涵底面的ux，uz两个方向，感觉多余了
d,all,ux
d,all,uz
allsel,all
!!
acel,0,10
sfa,5,,pres,,-63060 !箱涵顶面荷载，暂时按均布考虑，不过这个值
                    !是否包含顶板的混凝土自重？
allsel,all
sfgrad,pres,0,y,0,3720
sfa,4,,pres,-40000    !对这边施加的方向有点疑惑，怎么看到施加的效果呢？
sfgrad,pres,0,y,0,3720
sfa,6,,pres,-23400
sfa,3,,pres,-85560   !底板的反力，加了弹簧的话，是否还要加反力呢？
/psf,pres,norm,2,0,1
!!
nropt,full
allsel
outres,all,all
solve
fini
/post1
pldisp,1
etable,,smisc,4
pletab,smis4,-1,0
etable,,smisc,7
pletab,smis7,-1,0 !弯距的显示结果好像有很大问题
etable,,smisc,1
pletab,smis1,1,0
etable,,smisc,2
etable,,smisc,7
etable,,smisc,8
plls,smis1,smis2,1,0  ！轴力
plls,smis27,smis8,1   ！剪力
finish   !下面考虑的是杀死所有弹簧单元分析的结果
/solu
allsel,all
esel,s,type,,2
ekill,all
allsel,all
/post1
pldisp,1
etable,,smisc,7
pletab,smis7,-1,0 ! 感觉没多大差别
etable,,smisc,1
pletab,smis1,1,0
etable,,smisc,2
etable,,smisc,7
etable,,smisc,8
plls,smis1,smis2,1,0  ！轴力
plls,smis27,smis8,1   ！剪力

xugangseu

没人关注箱涵吗？