马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
*set,h,100 !板高
*set,b,560 !板宽
*set,l,3270 !板长
*set,lz,100 !支座长
*set,zg,20 !支座高
*set,a,25 !保护层厚度
*set,dis,120 !钢筋间距
*set,zj,8 !钢筋直径
*set,pi,acos(-1)
*set,sr,pi*(zj/2)**2 !钢筋面积
*set,cb,300 !碳纤维布宽度
*set,ch,0.103 !碳纤维布厚度
*set,cl,l-2*lz !碳纤维布长度
*set,f1,3.14e-3
*set,f2,5.58e-3
/prep7
et,1,solid65 !设置混凝土单元
et,2,link8 !设置钢筋单元
et,3,shell41 !设置碳纤维布单元
keyopt,3,1,2
et,4,solid45 !设置支座单元
mp,ex,1,3.1e4
mp,prxy,1,0.2
tb,kinh,1,1,7
tbpt,,0.0001,3.1
tbpt,,0.0004,11.477
tbpt,,0.0008,20.40
tbpt,,0.0012,26.79
tbpt,,0.0016,30.60
tbpt,,0.002,31.88
tbpt,,0.0033,31.88
tb,conc,1,1,9
tbdata,,0.3,0.5,1.75,-1
mp,ex,2,2e5
mp,prxy,2,0.3
tb,bkin,2,1,2,1
tbdata,,265,0
mp,ex,3,2.8e5
mp,prxy,3,0
mp,ex,4,2.1e5
mp,prxy,4,0.3
r,1,2,0,0,0,,,
rmore,,,,,,,
r,2,sr,,
r,3,sr/2,,
r,4,ch,,,,,,
rmore,,
block,0,-b/2,0,h,(l/2-lz),-(l/2-lz)
lgen,2,8,,,b/2-2*dis
lgen,2,13,,,dis
adrag,13,,,,,,9
adrag,14,,,,,,9
vsba,1,7
vsba,3,8
lsel,s,loc,z,(l/2-lz)
lsel,r,loc,y,0
lgen,2,all,,,,a
lsel,s,loc,z,(l/2-lz)
lsel,r,loc,y,a
adrag,all,,,,,,9
allsel,all
vsba,2,15
vsba,4,14
vsba,1,13
lsel,s,loc,y,0 !划分混凝土单元
lsel,a,loc,y,a
lsel,a,loc,y,h
lsel,u,loc,z,(l/2-lz)
lsel,u,loc,z,-(l/2-lz)
lesize,all,,,50
lsel,s,loc,y,0
lsel,a,loc,y,a
lsel,a,loc,y,h
lsel,u,loc,x,0
lsel,u,loc,x,-dis
lsel,u,loc,x,-2*dis
lsel,u,loc,x,-b/2
lsel,r,loc,x,-2*dis,0
lesize,all,30
allsel
lsel,s,loc,y,0
lsel,a,loc,y,a
lsel,a,loc,y,h
lsel,u,loc,x,0
lsel,u,loc,x,-dis
lsel,u,loc,x,-2*dis
lsel,u,loc,x,-b/2
lsel,r,loc,x,-b/2,-2*dis
lesize,all,40
lsel,s,loc,z,(l/2-lz)
lsel,a,loc,z,-(l/2-lz)
lsel,u,loc,y,0
lsel,u,loc,y,a
lsel,u,loc,y,h
lesize,all,25
allsel,all
type,1
mat,1
real,1
esys,0
mshape,0,3d
mshkey,1
vmesh,all
allsel,all
asel,s,loc,z,(l/2-lz) !生成左支座上混凝土单元
type,1
extopt,esize,3,0
extopt,aclear,1
extopt,attr,0,0,0
mat,1
real,1
esys,0
vext,all,,,0,0,lz
asel,s,loc,z,-(l/2-lz) !生成右支座上混凝土单元
type,1
extopt,esize,3,0
extopt,aclear,1
extopt,attr,0,0,0
mat,1
real,1
esys,0
vext,all,,,0,0,-lz
allsel,all
lsel,s,loc,x,-2*dis !划分板计算长度内钢筋单元
lsel,a,loc,x,-dis
lsel,r,loc,y,a
lsel,r,loc,z,-(l/2-lz),(l/2-lz)
type,2
mat,2
real,2
esys,0
lesize,all,,,,50
lmesh,all
allsel,all
lsel,s,loc,x,0
lsel,r,loc,y,a
lsel,r,loc,z,-(l/2-lz),(l/2-lz)
type,2
mat,2
real,3
esys,0
lesize,all,,,,50
lmesh,all
allsel,all
lsel,s,loc,x,-2*dis !划分板支座上钢筋单元
lsel,a,loc,x,-dis
lsel,r,loc,y,a
lsel,r,loc,z,(l/2-lz),l/2
type,2
mat,2
real,2
esys,0
lesize,all,,,,3
lmesh,all
allsel,all
lsel,s,loc,x,0
lsel,r,loc,y,a
lsel,r,loc,z,-l/2,-(l/2-lz)
type,2
mat,2
real,3
esys,0
lesize,all,,,,3
lmesh,all
allsel,all
asel,s,loc,y,0 !生成支座
asel,u,loc,z,-(l/2-lz),(l/2-lz)
type,4
extopt,esize,2,0
extopt,aclear,1
extopt,attr,0,0,0
mat,4
esys,0
vext,all,,,,-zg
allsel,all
allsel,all !碳纤维布
wpstyle,,,,,,,,1
wprota,,-90
rectng,0,-cb/2,(l/2-lz),-(l/2-lz)
wpstyle,,,,,,,,0
allsel,all
*get,mian,area,,num,max !划分碳纤维布
asel,s,area,,mian
lsla,s
lsel,u,loc,z,(l/2-lz)
lsel,u,loc,z,-(l/2-lz)
lesize,all,,,50
allsel,all
asel,s,area,,mian
lsla,s
lsel,u,loc,x,0
lsel,u,loc,x,-cb/2
lesize,all,30
allsel,all
asel,s,area,,mian
type,3
mat,3
real,4
esys,0
mshape,1,2d
amesh,all
nummrg,node
numcmp,node
/solu
nsel,s,loc,x,0
dsym,symm,x
allsel,all
nsel,s,loc,z,l/2-lz
nsel,r,loc,y,-zg
d,all,,,,,,uy,uz
allsel,all
nsel,s,loc,z,-(l/2-lz)
nsel,r,loc,y,-zg
d,all,,,,,,uy
allsel,all
asel,s,loc,y,h
asel,r,loc,z,-(l/2-lz),(l/2-lz)
sfa,all,1,pres,f1
antype,0
time,1
nlgeom,1
autots,1
nropt,full
eqslv,spar
outres,all,all
nsubst,20
cnvtol,f,,0.05,2
solve
把子布等改了好多次,就是不收敛,请各位大侠帮帮我 |