马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
做了一个试验,120mm*250mm截面,跨度为1800mm,纵筋为不锈刚,直径为13mm,参数为:E=1.9e11pa,u=0.3,
抗拉强度为3.25e8pa.架立钢筋为普通钢筋,直径8mm,箍筋为6mm◎150mm,混凝土为C30。在离梁端面150mm处简支,荷载离支座500mm处对称荷载。
各位如果有兴趣的话帮我看看为什么不收敛。
命令流如下:
/filname,001,on
/title,001
/prep7
!单元选择
et,1,link8
et,2,solid65
et,3,plane42
!定义实常数
r,1,132.73e-6
r,2,50.27e-6
r,3,28.27e-6
r,4
!材料定义
mp,ex,1,1.9e11
mp,prxy,1,0.3
tb,bkin,1,1,2,1
tbdata,,3.25e8,0.0
mp,ex,2,2.1e11
mp,prxy,2,0.25
tb,bkin,2,1,2,1
tbdata,,2.35e8,0.0
mp,ex,3,3.0e10
mp,prxy,3,0.2
tb,kinh,3,1,8,
tbtemp,0
tbpt,,0.0001,3.0e6
tbpt,,0.0003,3.9683e6
tbpt,,0.0006,7.2930e6
tbpt,,0.0009,9.9743e6
tbpt,,0.0012,12.012e6
tbpt,,0.0016,13.728e6
tbpt,,0.002,14.3e6
tbpt,,0.0033,14.3e6
tb,conc,3,1,9,
tbdata,,0.5,0.55,2.0e6,-1
!创建节点
n,1,0,0,0
n,13,0.12,0,0
fill,1,13,11
ngen,26,20,1,13,1,0,0.01,0
ngen,37,1000,all,0,0,0,0,-0.05
!创建箍筋
type,1
mat,2
real,3
*do,i,0,36,3
*do,j,64,69,1
e,1000*i+j,1000*i+j+1
*enddo
*do,j,70,430,20
e,1000*i+j,1000*i+j+20
*enddo
*do,j,450,445,-1
e,1000*i+j,1000*i+j-1
*enddo
*do,j,444,84,-20
e,1000*i+j,1000*i+j-20
*enddo
*enddo
!创建架立钢筋
type,1
mat,2
real,2
*do,i,444,35444,1000
e,i,i+1000
*enddo
*do,i,450,35450,1000
e,i,i+1000
*enddo
!创建纵筋
type,1
mat,1
real,1
*do,i,64,35064,1000
e,i,i+1000
*enddo
*do,i,70,35070,1000
e,i,i+1000
*enddo
!创建拉伸平面
k,1,0,0,0
k,2,0.12,0,0
k,3,0.12,0.25,0
k,4,0,0.25,0
a,1,2,3,4
lplot
/pnum,line,1
/replot
type,3
lsel,s,line,,1,4,1
lesize,all,0.01
lsel,all
amesh,1
extopt,attr,0,0,0
mat,3
real,4
esys,0
extopt,esize,36,1.0
extopt,aclear,1
vext,all, , ,0,0,-1.8,,,,
/view,1,1,1,1
/replot
!节点合并与压缩
nummrg,all,,,,low
numcmp,all
!施加约束
nsel,s,loc,z,-0.15
nsel,r,loc,y,0
d,all,all
allsel,all
nsel,s,loc,z,-1.65
nsel,r,loc,y,0
d,all,ux,uy
allsel,all
!施加荷载
nsel,s,loc,z,-0.65
nsel,r,loc,y,0.25
f,all,fy,-1000
allsel,all
nsel,s,loc,z,-1.15
nsel,r,loc,y,0.25
f,all,fy,-1000
allsel,all
fini
/solu
!分析选项设置
antype,0
nlgeom,1
nropt,full,,
eqslv,spar,,0,
time,2.0
autots,1
nsubst,100,200,50,1
kbc,0
cnvtol,u,,0.03,0,,
neqit,50,
pred,on,,on
outress,all,all
欢迎各位参与讨论。 |