frp钢筋混凝土加固梁的粘结滑移问题

sarina

我最近在做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.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,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        !选择体2
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.            

*** ERROR ***     SUPPRESSED MESSAGE    CP =       7.375   TIME= 21:16:15
Node 1872 direction UY included in coupled sets 696 and 275.            

*** ERROR ***     SUPPRESSED MESSAGE    CP =       7.375   TIME= 21:16:15
Node 155 direction UY included in coupled sets 697 and 276.            

*** ERROR ***     SUPPRESSED MESSAGE    CP =       7.375   TIME= 21:16:15
Node 1866 direction UY included in coupled sets 697 and 276.            

*** ERROR ***     SUPPRESSED MESSAGE    CP =       7.391   TIME= 21:16:15
The step data was checked and there were errors found.                  
  Please check output or errors file ( G:\ansys calculating              
result\instances of zhaotong\BM1-4\my\9-12.err ) for messages.

sarina

先自己顶一下了。
高手们，给指点一下了！
谢谢

Szncu

sarina 发表于 2008-9-13 17:37
先自己顶一下了。
高手们，给指点一下了！
谢谢

求教一下，FRP与混凝土是如何接触的，或者说如何通过GUI使两者粘结上，粘结强度如何设置？