finish
/clear
antype,static
et,1,solid65
r,1,
TBTEMP,0
TBTEMP,0
TBDATA,,0,1,0.4e6,-1,,
lesize,i,50
*do,i,9,12,1
lesize,i,100
vmesh,1
et,2,link8 !是否应用PIPE 单元
r,2,0.000576
mp,ex,2,2.1e11
k,11,150,150
lmesh,13
keyopt,3,3,3 !定义Z方向弹簧
r,3,0.012,340.98,0.023144013,676.822,0.041007667,1009.222
rmore,0.060508355,1329.6,0.078506,1669.288667,0.13106,2006.533333
rmore,0.200204,2332.511333,0.297857,2678.866667,0.381036,3004.466667
rmore,0.491887,3340.466667,0.634276,3667.711333,2.961887,38.733333
rmore,2.9726,3982.733333,3.05323,4143.8,3.1414,4204.355333
rmore,3.16088,4294.266667,3.232973,4385.4667,3.244492,4473.044667
rmore,3.333723,4047.6,
e,37,736
e,86,737
save
et,4,shell63
r,4,0.015,0.015,0.015,0.015
mp,ex,4,2.5e4
mp,prxy,4,0.3
type,4
asel,,,,3,6,1
aatt,4,4,4,
amesh,3,6,1
allsel
nsel,s,loc,z,1400
nsel,r,loc,x,0
d,all,all
allsel
nsel,s,loc,z,1400
nsel,r,loc,x,300
d,all,all
allsel
time,0
*do,i,1,11,1
time,i
f,737,fz,5000*i
cnvtol,f,,0.1,2
nropt,full
fdele,all
*do,i,12,25,1
f,737,fz,55000+2500*(i-11)
time,i
cnvtol,f,,0.05,2
nropt,full
fdele,all
*enddo
ANSYS中使用COMBIN39单元的实例2
/PREP7
R,1,0.0,0.0,.1,.204,.2,.432 ! SPRING DATA
RMORE,.3,.708,.4,1.056,.5,1.5
RMORE,.6,2.0,.7,2.772,.8,3.8
RMORE,.9,4.716,1.0,6.0
N,1
N,2
E,1,2
TYPE,2
REAL,2
E,2
D,1,ALL
D,2,UX
SAVE
FINISH
SOLCONTROL,0
OUTRES,NSOL,1
NSUBST,5
OUTPR,BASIC,NONE
NSUBST,40
OUTPR,BASIC,LAST
FINISH
/POST26
TIMERANGE,.003,.18
*GET,PER,VARI,2,EXTREM,TMIN
*status,parm
FINISH
frp钢筋混凝土加固梁的粘结滑移问题
我最近在做frp钢筋混凝土加固梁的剥离问题，采用ansys建立模型。solid65模拟混凝土，combin39单元来模拟frp与混凝土之间的粘结滑移，有一些错误，希望高手们有时间帮我看看。多谢了
很多人说，采用ansys算粘结滑移很难算到收敛，但是不知道有没有算出来的，只要有可能，我就不放弃了。
命令流如下：
!第一组
/prep7
!定义参数和单元类型
!混凝土参数
*SET,hntra,20.5
*SET,ec,28800
*SET,possible,0.3
*SET,hntrl,2.35
!FRP参数
*SET,ff,1800 !抗拉强度
*SET,ef,220000 !弹性模量
!FRP纤维布的极限拉应变
*SET,ckf,ff/ ef
!FRP厚度
*SET,tf,0.121*2
et,1,solid65 !混凝土单元
KEYOPT,1,1,1 !不考虑大变形
KEYOPT,1,7,1 !考虑应力松弛，有助于计算收敛
et,2,link8 !钢筋单元
et,3,shell43 !模拟frp布
et,4,solid45 !刚性垫块
!定义实常数
r,1,50.3 !纵向受拉钢筋实常数
r,2,28.3 !纵向受压钢筋，箍筋实常数
r,3 !混凝土
R,4,tf,tf,tf,tf,,, !frp的实常数
RMORE,,
r,5 !垫块
ET,5,COMBIN39 !法向弹簧
keyopt,5,1,0
keyopt,5,2,1
KEYOPT,5,3,3
r,6,0,0,0.0005,29340,0.001,41577 !粘结滑移关系
rmore,0.005,92970,0.01,131478,0.015,161028
rmore,0.02,185939,0.025,207886,0.029214,224725
!定义材料属性
mp,ex,1,ec !混凝土材料属性
mp,prxy,1,possible
tb,conc,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.
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,1.94e5 !纵向受拉钢筋
mp,prxy,2,0.3
tb,bkin,2,1,2,1
tbdata,,368.3,2000
mp,ex,3,1.97e5 !横向箍筋，受压钢筋材料属性
mp,prxy,3,0.3
tb,bkin,3,1,2,1
tbdata,,410.6,2000
!定义FRP
mp,ex,4,ef
mp,prxy,4,0.3
mp,ex,5,ec !垫块材料属性
mp,prxy,5,possible
!建模
blc4,,,100,200,-1100
k,9,20,40,0 !钢筋
k,10,20,40,-1100
l,9,10
lgen,2,13,,,60
lsel,s,line,,13,14
LPLOT
CM,zlj,LINE
allsel,all
lgen,2,13,,,,120
lgen,2,14,,,,120
lsel,s,line,,15,16
lplot
cm,zyj,line
allsel,all
k,17,20,40,-500 !箍筋
k,18,80,40,-500
k,19,20,160,-500
l,17,18
l,17,19
lgen,2,17,,,,120
lgen,2,18,,,60
lsel,s,line,,17,20,1
LPLOT
*do,ii,1,5
lgen,2,17,20,1,,,-ii*100
ii=ii+1
*enddo
LPLOT
cm,gj,line
WPAVE,0,0,0 !CFRP
wprota,,90 !旋转工作平面yz，90°
wpoffs,,,0.835 !偏移工作平面z cfrp面与混凝土面之间距离为0.835mm
blc4,,,100,-1000
wpoffs,,,-0.835 !支座面与cfp下表面距离为0.835mm
wpoffs,,-1100
blc4,,,100,100,20 !支座处混凝土垫块100*100*20
wpoffs,,600
wpoffs,,,-200
wprota,,-90
blc4,,,100,20,-100 !荷载处混凝土垫块是100*100*20
allsel,all
!划分单元
VSEL,s,volu,,1 !混凝土 网格划分较细为25*50*25
VATT,1,3,1
LSEL,s,line,,12
LSEL,a,line,,7
LSEL,a,line,,8
LESIZE,12,50,,,,,,,0
LESIZE,7,20,,,,,,,0
LESIZE,8,20,,,,,,,0
MSHAPE,0,3D
MSHKEY,1
VMESH,1
cmsel,s,zlj,line !主la筋划分为50mm
lATT,2,1,2
LESIZE,all,50,,,,,,,0
LMESH,all
cmsel,s,zyj,line !主ya筋划分为50mm
lATT,3,2,2
LESIZE,all,50,,,,,,,0
LMESH,all
cmsel,s,gj,line !箍筋为20mm
lplot
lATT,3,2,2
LESIZE,all,20,,,,,,,0
LMESH,all
asel,s,area,,7 !CFRP划分单元为四边形20*50
aplot
AATT,3,4,3
lsel,s,,,42
lsel,a,,,43
LESIZE,42,50,,,,,,,0
LESIZE,43,20,,,,,,,0
mshape,0,2D
MSHKEY,1
AMESH,7
allsel,all !划分支座垫块为25*50*20
VSEL,s,,,2 !选择体2
VATT,5,3,4
LESIZE,50,50,,,,,,,0
LESIZE,51,20,,,,,,,0
LESIZE,55,20,,,,,,,0
LESIZE,47,50,,,,,,,0
LESIZE,49,20,,,,,,,0
LESIZE,53,20,,,,,,,0
MSHAPE,0,3D
MSHKEY,1
VMESH,all
/VIEW,1,1,1,1
allsel,all !划分支座垫块为25*50*20
VSEL,s,,,3 !选择体3
VATT,5,3,4
LESIZE,67,50,,,,,,,0
LESIZE,62,20,,,,,,,0
LESIZE,61,20,,,,,,,0
LESIZE,66,50,,,,,,,0
LESIZE,57,20,,,,,,,0
LESIZE,58,20,,,,,,,0
MSHAPE,0,3D
MSHKEY,1
VMESH,all
allsel,all
nsel,s,loc,y,-0.835,0 !弹簧单元0
nplot
TYPE,5
REAL,6
ESYS,0
eintf,0.84,,low !结点号码从最低点到最高点编号
allsel,all
nsel,s,loc,y,-1.67,0 !混凝土和CFRP界面重合节点 约束UX UY方向
CPINTF,UX,1.67
CPINTF,Uy,1.67
nsel,s,loc,y,-1.67,-0.835 ！cfrp内部结点耦合
CPINTF,Uz,0.835
ALLSEL,ALL !其他节点耦合
CPINTF,UX
CPINTF,Uy
CPINTF,Uz
allsel,all
!求解设置
/solu
antype,0
allsel,all
asel,s,,,9 !右边界约束
da,all,ux
da,all,uy
allsel,all
asel,s,loc,z,0 !对称面约束
da,all,symm
allsel,all !施加外部载荷
asel,s,,,17
sfa,all,,pres,2.456
allsel,all
autots,on
time,1
nsubst,100
lnsrch,on
pred,on
outres,basic,all
ncnv,2
allsel,all
SOLVE
FINISH
有一些错误，请高手们帮我看看
*** ERROR *** SUPPRESSED MESSAGE CP = 7.375 TIME= 21:16:15
Node 154 direction UY included in coupled sets 696 and 275. 下载本文
