最新整理ansys实例混凝土加弹簧Word文档下载推荐.docx
《最新整理ansys实例混凝土加弹簧Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《最新整理ansys实例混凝土加弹簧Word文档下载推荐.docx(14页珍藏版)》请在冰豆网上搜索。
8
*SET,QX,100!
水平荷载
*SET,QZ,100
*SET,b0,1
*SET,h0,0.5
!
进入前处理器
/prep7
et,1,solid45!
桩身单元
et,2,combin14!
弹簧单元
et,3,surf154!
表面效应单元
keyopt,1,1,0
keyopt,2,3,0!
作为纵向三维弹簧
keyopt,3,4,1!
不包含中间节点
r,1
mp,ex,1,3。
25e10!
桩的材料属性
mp,nuxy,1,0。
3
mp,dens,1,2700
建立整个桩模型的节点
n,1
n,6,1
fill,1,6
ngen,6,6,1,6,1,,,0。
2
ngen,21,36,1,36,1,,0.5
/view,1,1,1,1!
以正等侧显示
/ang,1
/REP,FAST
GPLOT
建立桩的几何实体模型
k,1
k,2,1
k,3,1,,1
k,4,,,1
k,5,,10,
k,6,1,10,
k,7,1,10,1
k,8,,10,1
/PNUM,kp,1
kplot
v,1,2,3,4,5,6,7,8
lsel,s,loc,y,0
lsel,a,loc,y,10
lesize,all,,,5
allsel
lsel,s,loc,y,5
lesize,all,,,20
vmesh,1
nummrg,all!
压缩重复元素
numcmp,all!
合并重复元素
建立第一组弹簧模型,共7*16=102个弹簧
nsel,s,loc,x,1
nsel,r,loc,y,—0。
1,8.1
nplot
ngen,2,1000,all,,,2
type,2
*do,i,1,17,1
*if,i,le,4,then!
判断属于哪一层土
kk=1
*elseif,i,le,8,then
kk=2
*elseif,i,le,12,then
kk=3
*elseif,i,le,17,then
kk=4
*endif
*do,j,6,36,6
*get,vy,node,j+(i-1)*36,loc,y!
获取节点的Y坐标值
r,kk+1,mk(1,kk,1)*abs(vy—8)*b0*h0/6!
根据Y值计算水平抗力系数,即弹簧的刚度
real,kk+1
e,j+(i—1)*36,1000+j+(i—1)*36
*enddo
最大单元号602
第二组弹簧单元,共7*16=102个弹簧
nsel,s,loc,x,0
nsel,r,loc,y,-0。
nsel,r,loc,z,-0.1,1。
1
ngen,2,2000,all,,,-2
*do,i,1,17,1
*if,i,le,4,then
*elseif,i,le,8,then
*elseif,i,le,12,then
*elseif,i,le,17,then
*do,j,1,31,6
*get,vy,node,j+(i-1)*36,loc,y
r,kk+10,mk(1,kk,1)*abs(vy—8)*b0*h0/6
real,kk+10
e,j+(i—1)*36,2000+j+(i-1)*36
*enddo
最大单元号704
第三组弹簧单元,共7*16=102个弹簧
local,11,0,,,,,90
csys,0
nsel,s,loc,z,1
nsel,r,loc,y,—0.1,8.1
nsel,r,loc,x,-0.1,1.1
ngen,2,3000,all,,,,,2
*do,i,1,17,1
*if,i,le,4,then
*elseif,i,le,8,then
*elseif,i,le,12,then
*elseif,i,le,17,then
*endif
*do,j,31,36,1
*get,vy,node,j+(i—1)*36,loc,y
r,kk+20,mk(1,kk,1)*abs(vy—8)*b0*h0/6
real,kk+20
esys,0
e,j+(i-1)*36,3000+j+(i—1)*36
最大单元号码806
第四组弹簧单元,共7*16=102个弹簧
nsel,s,loc,z,0
nsel,r,loc,y,—0.1,8.1
nsel,r,loc,x,—0。
1,1。
ngen,2,4000,all,,,,,-2
*do,i,1,17,1
*if,i,le,4,then
*elseif,i,le,8,then
*elseif,i,le,12,then
*elseif,i,le,17,then
*do,j,1,6,1
*get,vy,node,j+(i—1)*36,loc,y
r,kk+20,mk(1,kk,1)*abs(vy-8)*b0*h0/6
real,kk+20
esys,0
e,j+(i-1)*36,4000+j+(i—1)*36
最大单元号码为908
底层弹簧
csys,0
ngen,2,5000,1,36,1,,—1
r,1000,mk(1,1,1)*8*b0*h0/6*20
real,1000
*do,i,1,36,1
e,i,5000+i
gplot
第一组表面效应单元
r,1001
real,1001
local,13,0,,,,90
type,3
*do,i,1,16,1
*do,j,6,30,6
nsel,s,node,,j+(i—1)*36
nsel,a,node,,j+6+(i-1)*36
nsel,a,node,,j+42+(i—1)*36
nsel,a,node,,j+36+(i—1)*36
esys,13
esurf
第二组表面效应单元
local,13,0,,,,90
*do,i,1,16,1
*do,j,1,29,6
nsel,a,node,,j+6+(i-1)*36
nsel,a,node,,j+42+(i-1)*36
nsel,a,node,,j+36+(i—1)*36
esys,13
第三组表面效应单元
local,14,0,,,,90
*do,i,1,16,1
*do,j,1,5,1
nsel,s,node,,j+(i-1)*36
nsel,a,node,,j+1+(i—1)*36
nsel,a,node,,j+37+(i-1)*36
nsel,a,node,,j+36+(i—1)*36
esys,14
第四组表面效应单元
local,15,0,,,,90
*do,i,1,16,1
*do,j,31,35,1
nsel,s,node,,j+(i—1)*36
nsel,a,node,,j+1+(i—1)*36
nsel,a,node,,j+37+(i—1)*36
nsel,a,node,,j+36+(i—1)*36
esys,15
save
finish
进入求解器
/solu
nsel,s,loc,x,3
nsel,a,loc,x,-2
nsel,a,loc,z,3
nsel,a,loc,z,—2
d,all,all!
约束边界条件
nsel,a,loc,y,-1
d,all,all
nsel,a,loc,y,0
d,all,ux
d,all,uz
ACEL,0,9.81,0,!
施加重力加速度
桩顶表面平均分布载荷
nsel,s,loc,y,10
nsel,r,loc,x,—0。
nsel,r,loc,z,-0.1,1。
esln,s
sfe,all,6,pres,,5000!
桩顶表面分布荷载
nsel,r,loc,x,1
f,all,fx,QX!
水平载荷
nsel,s,loc,y,10
nsel,r,loc,z,1
f,all,fz,QZ!
第四层土的表面正摩擦力
nsel,s,loc,y,0,1。
5
esln,s
esel,r,type,,3
sfe,all,2,pres,,FQ(1,1,1)
/PSF,PRES,TANx,2,0,1
第三层土的表面负摩擦力
nsel,s,loc,y,2。
5,3。
esel,r,type,,3
sfe,all,2,pres,,FQ(1,2,1)
/PSF,PRES,TANx,2,0,1
第二层土的表面负摩擦力
nsel,s,loc,y,4。
5,5。
esel,r,type,,3
sfe,all,2,pres,,FQ(1,3,1)
/PSF,PRES,TANx,2,0,1
第一层土的表面负摩擦力
nsel,s,loc,y,6.5,7。
sfe,all,2,pres,,FQ(1,4,1)
/PSF,PRES,TANx,2,0,1
time,1
autot,on
nlgeom,on
nropt,full
solve
/POST1
PLNSOL,U,Y,1,1
PLNSOL,S,Z,1,1
PATH,path1,2,30,20,
PPATH,P51X,1
PATH,STAT
/PBC,PATH,1
/REPLOT
/PBC,PATH,0
AVPRIN,0,0,
PDEF,SY,S,Y,AVG
/PBC,PATH,,0
PRANGE,1,1,10,YG,
PLPATH,SY