combin39粘结滑移问题 ，急
我在做FRP加固混凝土梁剥离破坏模拟时，采用combin39单元模拟粘结滑移，单元keyopt(3)=UZ,模拟Z方向滑移，弹簧单元是在重合的节点上建立的，把x，y向进行耦合，为什么计算的时候 老是警告约束不够呢，而且计算很慢，但是当我也耦合节点UZ方向的时候，就没有那个警告了，速度也快多了，但是我这样做好像没有起到弹簧的作用，到底是怎么回事啊啊？我的QQ：50343975
/config,nres,50000
/config,fsplit,750
...........................设置参数
/prep7
ET,1,SOLID65 混凝土
ET,2,LINK8
ET,3,SHELL41 CFRP
ET,4,SOLID45 垫块
et,5,combin39 弹簧单元
keyopt,5,1,0
keyopt,5,2,1
KEYOPT,5,3,3
keyopt,5,4,0
keyopt,1,1,1
keyopt,1,7,1
keyopt,3,1,2
*SET,rd0,6
*SET,rd1,12
*SET,hntra,20.1 钢筋参数
*SET,hntrl,2.01
r,1
r,2,acos(-1)*0.25*rd0*rd0
r,3,acos(-1)*0.25*rd1*rd1
r,5,0.0005,0.2934,0.001,0.41577,0.005,0.9297, 粘结滑移关系
rmore,0.01,1.31478,0.015,1.61028,0.02,1.85939,
rmore,0.025,2.07886,0.029214,2.24725
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,5,,3240
MPDATA,PRXY,5,,0.3
mp,ex,1,3e4
mp,prxy,1,0.2 混凝土
tb,concr,1
tbdata,,0.5,0.9,hntrl,-1
tb,miso,1,,23
tbpt,,0.000005,hntra*0.0074627
tbpt,,0.00001,hntra*0.00992537
tbpt,,0.00005,hntra*0.0491293
tbpt,,0.0001,hntra* 0.0970149
tbpt,,0.0002,hntra*0.19
tbpt,,0.0004,hntra*0.36
tbpt,,0.0006,hntra*0.51
tbpt,,0.0008,hntra*0.64
tbpt,,0.0010,hntra*0.75
tbpt,,0.0012,hntra*0.84
tbpt,,0.0014,hntra*0.91
tbpt,,0.0016,hntra*0.96
tbpt,,0.0018,hntra*0.99
tbpt,,0.0020,hntra*1.00
tbpt,,0.0022,hntra*0.98
tbpt,,0.0024,hntra*0.97
tbpt,,0.0026,hntra*0.95
tbpt,,0.0028,hntra*0.93
tbpt,,0.003,hntra*0.92
tbpt,,0.0032,hntra*0.9
tbpt,,0.0034,hntra*0.88
tbpt,,0.0036,hntra*0.87
tbpt,,0.0038,hntra*0.85
mp,ex,2,2e5 钢筋
mp,prxy,2,0.3
TB,BKIN,2,1,2,1
TBDATA,,335,2000,,,,
mp,ex,3,2.3e5 CFRP
mp,prxy,3,0.3
R,4,0.167, , , , ,
MPTEMP,,,,,,,,
MPTEMP,1,0 弹簧
MPDATA,EX,5,,3240
MPDATA,PRXY,5,,0.3
mp,ex,4,2e5
mp,prxy,4,0.3
tb,bkin,4,1,2 垫块
tbdata,,335,2000
建模
blc4,,,100,200,1500 混凝土
/VIEW,1,1,1,1
/ANG,1
k,9,25,25,50 钢筋
k,10,25,25,1450
l,9,10
lgen,2,13,,,50,,
lgen,2,13,,,,150,
lgen,2,14,,,,150,
lsel,s,line,,13,16,1
LPLOT
CM,zj,LINE
allsel,all
k,17,25,25,120 箍筋
k,18,75,25,120
k,19,25,175,120
l,17,18
l,17,19
lgen,2,17,,,,150
lgen,2,18,,,50
lsel,s,line,,17,20,1
LPLOT
*do,ii,1,7
lgen,2,17,20,1,,,ii*180
ii=ii+1
*enddo
LPLOT
cm,gj,line
wpoffs,,,75 垫块
blc4,,,100,-20,50
vgen,2,2,,,,,1300
vgen,2,2,,,,220,450
vgen,2,4,,,,,400
WPAVE,0,0,0 CFRP
wpoffs,,,150
wprota,,90
blc4,,,100,1200
allsel,all
aplot
划分单元
/PREP7
VSEL,s,volu,,1 混凝土
VATT,1,1,1,0
LSEL,s,line,,12
LESIZE,12,50,,,,,,,0
LSEL,s,line,,7
LESIZE,7,25,,,,,,,0
LSEL,s,line,,8
LESIZE,8,25,,,,,,,0
MSHAPE,0,3D
MSHKEY,1
VMESH,1
VSEL,s,volu,,2,5,1 垫块
VATT,4,,4,0
LSEL,s,line,,55,91,12
LSEL,a,line,,59,95,12
LESIZE,all,25,,,,,,,0
LSEL,s,line,,54,90,12
LESIZE,all,10,,,,,,,0
MSHAPE,0,3D
MSHKEY,1
VMESH,all
LSEL,s,line,,13,14,1 纵筋
lATT,2,3,2,0
LSEL,s,line,,15,16,1
lATT,2,2,2,0
cmsel,s,zj,line
LESIZE,all,50,,,,,,,0
LMESH,all
cmsel,s,gj,line 箍筋
lplot
lATT,2,2,2,0
LESIZE,all,25,,,,,,,0
LMESH,all
asel,s,area,,31 CFRP
aplot
AATT,3,4,3,0,
LSEL,s,line,,98
LESIZE,98,50,,,,,,,0
LSEL,s,line,,99
LESIZE,99,25,,,,,,,0
mshape,1,2D
MSHKEY,1
aMESH,31
allsel,all
eplot
nsel,s,loc,y,0 弹簧单元
nsel,r,loc,z,150,1350
nplot
TYPE,5
MAT,5
REAL,5
ESYS,0
SECNUM,none,
TSHAP,LINE
eintf,0.0001,,low,,,,,
allsel,all
eplot
nsel,s,loc,y,0 混凝土和CFRP界面重合节点 约束UX UY方向
nsel,r,loc,z,150,1350
nplot
CPINTF,UX,0.0001,
CPINTF,Uy,0.0001,
ALLSEL,ALL 其他节点耦合
NSEL,U,LOC,Y,0
NPLOT
CPINTF,ALL,0.0001,
nsel,s,loc,y,0
nsel,u,loc,z,150,1350
NPLOT
CPINTF,ALL,0.0001,
VSEL,s,volu,,1
ESLV,S
VSEL,s,volu,,2,5,1
NSLV,S,1
nsel,r,loc,y,0,200
nsel,u,loc,z,100
nsel,u,loc,z,1400
nsel,u,loc,z,550
nsel,u,loc,z,950
CEINTF,0.25,ux,uy,uz, , , ,0,
eplot
VSEL,s,volu,,1
ESLV,S
cmsel,s,zj,line
ESLL,a
cmsel,s,gj,line
NSLL,S,1
nsel,u,loc,z,300
nsel,u,loc,z,1200
CEINTF,0.25,ux,uy,uz, , , ,0,
allsel,all
eplot
求解设置
/solu
antype,0
nsel,s,loc,z,100
nsel,a,loc,z,1400
nsel,r,loc,y,-20
d,all,uy
d,all,ux
nsel,r,loc,z,1400
d,all,uz
*do,ii,0,60
allsel,all
nsel,s,loc,z,550
nsel,a,loc,z,950
nsel,r,loc,y,220
f,all,fy,-(500+ii*200)
allsel,all
autots,on
time,ii+1
nsubst,200
lnsrch,on
pred,on
neqit,3000
outres,basic,all
ncnv,2
lswrite,ii+1
ii=ii+1
*enddo
allsel,all
EPLOT
/solu
allsel,all
LSSOLVE,1,61,1
页:
[1]