有关钢筋混凝土粘结滑移的问题
附件是我做的一根混凝土梁的粘结滑移，但是计算一直不收敛，跪求各位大神帮忙啊！！！完整命令流如下：
Finish
/clear
/prep7
!定义单元类型
et,1,solid65 !混凝土
et,2,link8 !钢筋
et,3,solid45 !垫块
et,4,combin39 ! 弹簧单元 （沿切向Z向）
keyopt,1,1,1 ! 不考虑大变形
keyopt,1,5,1 !混凝土的线性解输出控制：给出每个积分点的解
keyopt,1,6,3 !混凝土的非线性解输出控制：给出每个积分点的解
keyopt,1,7,1 ! 开裂后考虑应力松弛
keyopt,4,3,3 !沿Z 向平动
keyopt,4,1,0 ! 卸载路径与加载路径相同
keyopt,4,2,1 !不能承受压力
keyopt,4,4,0 ! 控制三唯轴向单元
!定义实常数
r,1 ! 混凝土实常数
r,2,5.026548E-5 !箍筋d8的截面积
r,3,3.1415926e-4 !纵筋d20的截面积
r,4,0,0,0.5e-3,3.575e6,0.8e-3,4.29e6,11e-3,1.43e6
!定义材料属性
!混凝土
mp,ex,1,2.99e10 !混凝土弹模
mp,prxy,1,0.2 !混凝土泊松比
tb,concr,1 !定义非线性材料特性
tbdata,,0.5,0.9,1.43e6,-1!张开系数，闭合，单轴抗拉强度，单轴抗压
tb,miso,1,1,8 !多线性等向强化
tbpt,,0.0001,2.99e6 !应变应力关系
tbpt,,0.0003,3.9683e6
tbpt,,0.0006,7.2930e6
tbpt,,0.0009,9.9743e6
tbpt,,0.0012,12.0120e6
tbpt,,0.0016,13.7280e6
tbpt,,0.0020,14.3e6
tbpt,,0.0033,14.3e6
!箍筋
mp,ex,2,2.1e11 !箍筋弹模
mp,prxy,2,0.25 !箍筋泊松比
tb,bkin,2 !双线性随动强化
tbdata,,300e6 !d8的屈服应力
!纵向钢筋
mp,ex,3,2e11 !纵筋弹模
mp,prxy,3,0.3 !纵筋泊松比
tb,bkin,3 !双线性随动强化
tbdata,,400e6 !d20的屈服应力
! 垫块
mp,ex,4,2e20
mp,prxy,4,0.3
!建立所有节点(200*400*3000)
!!!创建受压纵筋和箍筋
k,1,0.02,0.38,0
k,2,0.02,0.38,3
k,3,0.18,0.38,0
k,4,0.18,0.38,3
k,5,0.02,0.02
k,6,0.18,0.02
k,7,0.02,0.02,3
k,8,0.18,0.02,3
!!!箍筋，不考虑粘结滑移
l,1,3
l,1,5
l,5,6
l,6,3
lsel,all
lgen,16,all,,,,,0.2
lsel,all
latt,2,2,2
lsel,s,tan1,y
lesize,all,,,4
lmesh,all
lsel,s,tan1,x
lesize,all,,,6
lmesh,all
!!!受压纵筋，不考虑粘结滑移
l,1,2
l,3,4
latt,3,3,2
lsel,s,tan1,x
lsel,r,tan1,y
latt,3,3,2
lesize,all,,,30
lmesh,all
allsel,all
!!!创建混凝土
blc4,0,0,0.2,0.4,3
vatt,1,1,1
vsel,s,,,1
wprota,,,90
wpoffs,,,0.02
vsbw,all
wpoffs,,,0.16
vsbw,all
wprota,,90
wpoffs,,,-0.02
vsbw,all
wpoffs,,,-0.36
vsbw,all
lsel,s,loc,z,1.5
lesize,all,,,30
lsel,s,loc,z,0
lsel,r,loc,x,0.1
lesize,all,,,4
lsel,s,loc,z,0
lsel,r,loc,y,0.2
lesize,all,,,6
lsel,s,length,,0.02
lesize,all,,,1
vmesh,all
nummrg,node
nummrg,kp
numcmp,node
numcmp,kp
!!!创建受拉纵筋,考虑粘结滑移
k,100,0.02,0.02
k,101,0.18,0.02
k,102,0.02,0.02,3
k,103,0.18,0.02,3
l,100,102
l,101,103
lsel,s,loc,y,0.02
lsel,r,tan1,x
lsel,r,tan1,y
latt,3,3,2
lesize,all,,,30
lmesh,all
!!!定义弹簧单元
allsel,all
cpintf,ux,0.001
cpintf,uy,0.001
type,4
real,4
eintf,0.0001,,low,,,,,
allsel,all
!!!创建垫块
wpcsys,-1
blc4,0,0,0.2,-0.1,0.2
vsel,s,loc,y,-0.05
vgen,2,all,,,,,2.8
vsel,s,loc,y,-0.05
vatt,4,,3
esize,0.1
vmesh,all
!!!垫块与混凝土的约束
nsel,s,loc,y,0
cpintf,ux,0.0001
cpintf,uy,0.0001
cpintf,uz,0.0001
nsel,s,loc,y,-0.11
d,all,all
!!!加载点设置
blc4,0,0.4,0.2,0.1,0.2
vatt,4,,3
vsel,s,loc,y,0.45
vgen,2,all,,,,,1.40,,,1
esize,0.1
vmesh,all
!!加载点约束
nsel,s,loc,y,0.4
cpintf,ux,0.0001
cpintf,uy,0.0001
cpintf,uz,0.0001
!!!施加荷载
asel,s,loc,y,0.5
sfa,all,,pres,2000000
dtran
sftran
allsel,all
/solu
antype,0
autos,on
neqit,50 !设置每个荷载子步迭代次数为50
cnvtol,u,,0.015!采用位移收敛
nsubst,800
outres,all,all
solve
/post26
xx=node(0.1,0,1.5)
nsol,2,xx,u,y
prod,3,2,,,,,,-1
prod,4,1,,,,,,2000000*0.04
xvar,3
plvar,4 大神们，跪求啊，快来大神。 本人已算出，需要命令流的可以问我要。
页:
[1]