ansys 箱涵结构处理完全版

http://www.simwe.com/forum/viewthread.php?tid=806364&pid=1184459&page=1&extra=page%3D1#pid1184459

ansys 箱涵结构处理完全版


说明：
1、箱涵顶部填土厚度为0
2、单位制KN，m
3，箱涵顶部车辆荷载为局部荷载，作用面积0.4*0.8m，荷载700KN
4，shell43模型，取涵长度为10m，两侧土压力通长施加荷载，车辆产生的侧向荷载也通长施加

finish
/cle
/CWD,’F:\ansys shujubaoceng’
/nopr
keyw,pr_struc,1
/prep7
*set,l1,3.4 !箱涵跨径
*set,h1,2.8 !箱涵高度
*set,d1,0.3 ! 箱涵厚度
*set,k1,50000 !地基反应模量 取密实土壤
et,1,shell43
et,2,combin14
mp,ex,1,2.8e7 !混凝土弹性模量
mp,nuxy,1,0.2 !混凝土泊松比
mp,dens,1,2.5 !混凝土密度，单位取KN,m制
r,1,d1
!!
k,1
k,2,l1
k,3,l1,h1
k,4,,h1
lstr,1,2
lstr,2,3
lstr,3,4
lstr,4,1
a,1,2,3,4
VEXT,1,,,,,10
vdele,1
adele,1,2,1
!!
lsel,s,line,,9,12,1
lesize,all,0.1
lsel,s,line,,1,8,1
lesize,all,0.1
allsel,all
type,1
mat,1
real,1
local,11,0
csys,11
esys,11
amesh,all
/pnum,elem,1
!!
nsel,s,loc,y,0
ngen,2,20000,all,,,,-1 !复制产生弹簧底部节点
allsel,all
nplot
nummrg,all
numcmp,all

nsel,s,loc,y,0
*get,n1,node,,count
nsel,s,loc,y,-1
*get,n2,node,,num,min
allsel
*dim,rr,array,n1 !设定地基的弹性模量
*do,i,1,n1,1
rr(i)=arnode(i)
r,i+1,rr(i)*k1
type,2
mat,1
real,i+1
e,i,i+n2-1
*enddo
allsel
/pnum,node,0
/pnum,elem,0
/view,1,1,2,3
eplot !!模型
k,9,0,0,5 ! 定义关键点 便于后边路径选取
k,10,3.4,0,5
k,11,3.4,2.8,5
k,12,0,2.8,5
nkpt,,9
nkpt,,10
nkpt,,11
nkpt,,12
save
!!


/solve !!荷载和约束
nsel,s,loc,y,-1
d,all,all
nsel,s,loc,x,l1
nsel,r,loc,y,0
d,all,ux,0
allsel,all
sfcum,pres,add !!荷载累加
/psf,pres,norm,2,1
! acel,0,10
!!顶部汽车荷载
*Set,pv,700/0.8/0.4 !汽车作用面积0.8*0.4，荷载800KN
wpoffs,l1/2,h1,5
csys,4
nsel,s,loc,y,0
nsel,r,loc,x,-0.3,0.3
nsel,r,loc,z,-0.1,0.1
esln,s
sfe,all,2,pres,,pv
allsel
csys,11
wpave
!! 右侧土三角形荷载
*afun,deg
*Set,fai,30 ! 内摩擦角30度
*Set,r1,20 !土的饱和重度
sfgrad,pres,0,y,0,-r1*h1*(1-sin(fai))/h1
asel,s,area,,2,,,1
esel,r,elem,,all
sfe,all,2,pres,,r1*h1*(1-sin(fai))
allsel
!! 左侧土三?