粘结滑移本构
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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!箍筋
!对(NC-2)个中间节点施加粘结单元
*DO,i,2,NC-1
NSEL,S,,,nCrR(i)
NSEL,A,,,nNStR3T(i)
TYPE,5 !COMBIN39单元，只有x向自由度
REAL,i+nRDne!每个节点有一个实常数
MAT,5
E,nCrR(i),nNStR3T(i) !单元
NSEL,S,,,nCrR(i)
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
!定义Combin39实常数及粘结滑移本构模型
NC=10 !需要设置Combin39的节点数
*dim,nCrR,,NC !需耦合的砼节点号
*dim,nNStR3T,,NC !需耦合的钢筋节点号
!耦合首尾节点
cp,next,all,nCrR(1),nNStR3T(1)
cp,next,all,nCrR(NC),nNStR3T(NC)
*IF,xNStR3T(i),lt,xDfTu,THEN
TuC39N(i)= 9.52
*ELSE
TuC39N(i)= 15.2
*ENDIF
*ENDDO
*get,nRDne,ACTIVE,0,REAL, !得到已定义的实常数个数
!定义每个节点的新常数
*DO,i,2,cNSt-1
Ts= 3.16
Tcr=百度文库8.231
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
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
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
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
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
*DIM,RFR3T,ARRAY,11 !力
RDR3T(1)= -12, -9.72, -0.6624, -0.432, -0.0144, 0, 0.0144, 0.432, 0.6624, 9.72, 12
xDfTu=-141 !决定锚固Tu的值
*DIM,LC39R3T,,cNSt !每个节点的作用的有效钢筋长度
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
RFR3T(11)= Tr * CR3T * LC39R3T(i)
!定义实常数,一个节点对应一个常数
R,nRDne+i,RDR3T(1),RFR3T(1),RDR3T(2),RFR3T(2),RDR3T(3),RFR3T(3),RDR3T(4),RFR3T(4),RDR3T(5),RFR3T(5),RDR3T(6),RFR3T(6)
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
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,
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,,,,,
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
RFR3T(5)= -1* Ts * CR3T * LC39R3T(i)
RFR3T(6)= 0
RFR3T(7)= Ts * CR3T * LC39R3T(i)
RFR3T(8)= Tcr * CR3T * LC39R3T(i)
RFR3T(9)= Tcu * CR3T * LC39R3T(i)
RFR3T(10)= Tr * CR3T * LC39R3T(i)
NSEL,A,,,nNStR3T(i)
CP,NEXT,UY,nCrR(i),nNStR3T(i) !只在x向加弹簧,yz方向耦合
CP,NEXT,Uz,nCrR(i),nNStR3T(i)
*ENDDO
!****定义钢筋和混凝土之间的粘结滑移本构关系(徐有邻公式)
*DIM,RDR3T,ARRAY,11 !滑移
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
我的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!弹簧单元
*ENDDO
我在做FRP加固混凝土梁剥离破坏模拟时，采用combin39单元模拟粘结滑移，单元keyopt(3)=UZ,模拟Z方向滑移，弹簧单元是在重合的节点上建立的，把x，y向进行耦合，为什么计算的时候老是警告约束不够呢，而且计算很慢，但是当我也耦合节点UZ方向的时候，就没有那个警告了，速度也快多了，但是我这样做好像没有起到弹簧的作用，到底是怎么回事啊啊？
!划分单元
/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
RMORE,RDR3T(7),RFR3T(7),RDR3T(8),RFR3T(8),RDR3T(9),RFR3T(9),RDR3T(10),RFR3T(10),RDR3T(11),RFR3T(11)
!记录实常数，节点力
RFR3TBK(1,i)= RFR3T(1),RFR3T(2),RFR3T(3),RFR3T(4),RFR3T(5),RFR3T(6),RFR3T(7),RFR3T(8),RFR3T(9),RFR3T(10),RFR3T(11)
*DIM,TuC39N,,cNSt !每个节点的Tu
*DIM,RFR3TBk,,11,cNSt !记录全部节点实常数
!求每个节点的Tu
*DO,i,2,cNSt-1
LC39R3T(i)=(nx(nNStR3T(i+1)) - nx(nNStR3T(i-1))) / 2 !每个节点的作用长度
CR3T=18*3.1415/2 !钢筋周长
Tcu= TuC39N(i)
Tr= 3.136
RFR3T(1)= -1* Tr * CR3T * LC39R3T(i) !求曲线上的节点力
RFR3T(2)= -1* Tr * CR3T * LC39R3T(i)
RFR3T(3)= -1* Tcu * CR3T * LC39R3T(i)
RFR3T(4)= -1* Tcr * CR3T * LC39R3T(i)
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
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
/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!箍筋
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!箍筋
!对(NC-2)个中间节点施加粘结单元
*DO,i,2,NC-1
NSEL,S,,,nCrR(i)
NSEL,A,,,nNStR3T(i)
TYPE,5 !COMBIN39单元，只有x向自由度
REAL,i+nRDne!每个节点有一个实常数
MAT,5
E,nCrR(i),nNStR3T(i) !单元
NSEL,S,,,nCrR(i)
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
!定义Combin39实常数及粘结滑移本构模型
NC=10 !需要设置Combin39的节点数
*dim,nCrR,,NC !需耦合的砼节点号
*dim,nNStR3T,,NC !需耦合的钢筋节点号
!耦合首尾节点
cp,next,all,nCrR(1),nNStR3T(1)
cp,next,all,nCrR(NC),nNStR3T(NC)
*IF,xNStR3T(i),lt,xDfTu,THEN
TuC39N(i)= 9.52
*ELSE
TuC39N(i)= 15.2
*ENDIF
*ENDDO
*get,nRDne,ACTIVE,0,REAL, !得到已定义的实常数个数
!定义每个节点的新常数
*DO,i,2,cNSt-1
Ts= 3.16
Tcr=百度文库8.231
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
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
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
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
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
*DIM,RFR3T,ARRAY,11 !力
RDR3T(1)= -12, -9.72, -0.6624, -0.432, -0.0144, 0, 0.0144, 0.432, 0.6624, 9.72, 12
xDfTu=-141 !决定锚固Tu的值
*DIM,LC39R3T,,cNSt !每个节点的作用的有效钢筋长度
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
RFR3T(11)= Tr * CR3T * LC39R3T(i)
!定义实常数,一个节点对应一个常数
R,nRDne+i,RDR3T(1),RFR3T(1),RDR3T(2),RFR3T(2),RDR3T(3),RFR3T(3),RDR3T(4),RFR3T(4),RDR3T(5),RFR3T(5),RDR3T(6),RFR3T(6)
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
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,
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,,,,,
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
RFR3T(5)= -1* Ts * CR3T * LC39R3T(i)
RFR3T(6)= 0
RFR3T(7)= Ts * CR3T * LC39R3T(i)
RFR3T(8)= Tcr * CR3T * LC39R3T(i)
RFR3T(9)= Tcu * CR3T * LC39R3T(i)
RFR3T(10)= Tr * CR3T * LC39R3T(i)
NSEL,A,,,nNStR3T(i)
CP,NEXT,UY,nCrR(i),nNStR3T(i) !只在x向加弹簧,yz方向耦合
CP,NEXT,Uz,nCrR(i),nNStR3T(i)
*ENDDO
!****定义钢筋和混凝土之间的粘结滑移本构关系(徐有邻公式)
*DIM,RDR3T,ARRAY,11 !滑移
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
我的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!弹簧单元
*ENDDO
我在做FRP加固混凝土梁剥离破坏模拟时，采用combin39单元模拟粘结滑移，单元keyopt(3)=UZ,模拟Z方向滑移，弹簧单元是在重合的节点上建立的，把x，y向进行耦合，为什么计算的时候老是警告约束不够呢，而且计算很慢，但是当我也耦合节点UZ方向的时候，就没有那个警告了，速度也快多了，但是我这样做好像没有起到弹簧的作用，到底是怎么回事啊啊？
!划分单元
/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
RMORE,RDR3T(7),RFR3T(7),RDR3T(8),RFR3T(8),RDR3T(9),RFR3T(9),RDR3T(10),RFR3T(10),RDR3T(11),RFR3T(11)
!记录实常数，节点力
RFR3TBK(1,i)= RFR3T(1),RFR3T(2),RFR3T(3),RFR3T(4),RFR3T(5),RFR3T(6),RFR3T(7),RFR3T(8),RFR3T(9),RFR3T(10),RFR3T(11)
*DIM,TuC39N,,cNSt !每个节点的Tu
*DIM,RFR3TBk,,11,cNSt !记录全部节点实常数
!求每个节点的Tu
*DO,i,2,cNSt-1
LC39R3T(i)=(nx(nNStR3T(i+1)) - nx(nNStR3T(i-1))) / 2 !每个节点的作用长度
CR3T=18*3.1415/2 !钢筋周长
Tcu= TuC39N(i)
Tr= 3.136
RFR3T(1)= -1* Tr * CR3T * LC39R3T(i) !求曲线上的节点力
RFR3T(2)= -1* Tr * CR3T * LC39R3T(i)
RFR3T(3)= -1* Tcu * CR3T * LC39R3T(i)
RFR3T(4)= -1* Tcr * CR3T * LC39R3T(i)
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
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
/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!箍筋