我的solid65与shell41的接触有问题，请大神帮助

!角钢与混凝土接触 直接合并节点 约束200处的一条线 厚度为10
/title,TJ-4
finish
/clear
/prep7
/title,TJ-4
!定义----单元及材料
et,1,65
!混凝土单元
et,2,8
!箍筋单元
et,3,8
!纵筋单元
et,4,41
!纤
et,5,45
!垫块
et,6,45
!角钢
et,7,45
!螺栓
et,8,8
!U形箍筋
!材料属性定义
!材料1-混凝土
MP,EX,1,30000
MP,PRXY,1,0.2
*SET,fc,23.4
*SET,ft,2.24
tb,concr,1
!保证良好收敛，关闭压碎选项
tbdata,,0.5,0.95,ft,fc
TB,KINH,1,1,12,0
!mat=1，temperature=1,12个数据点，tbopt=0升高温度不会应力松弛
TBPT,,0.00005,1.5
tbpt,,0.0002,fc*0.19
tbpt,,0.0004,fc*0.36
tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64
tbpt,,0.0010,fc*0.75
tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0.91
tbpt,,0.0016,fc*0.96
tbpt,,0.0018,fc*0.99
tbpt,,0.0020,fc
tbpt,,0.0033,fc*0.85
/axlab,x,strain
/axlab,y,stress
tbplot,kinh,1,,,yes
!查看本构
!材料2-箍筋
MP,EX,2,2.1e5
MP,PRXY,2,0.3
TB,BKIN,2,1,2,1
TBDATA,,315
/axlab,x,strain
/axlab,y,stress
tbplot,bkin,2,,,yes
!查看本构
!材料3-纵筋
MP,EX,3,2.0e5
MP,PRXY,3,0.3
TB,BKIN,3,1,2,1
TBDATA,,378.6
/axlab,x,strain
/axlab,y,stress
tbplot,bkin,3,,,yes
!查看本构
!材料4-纤维布
!纤
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,4,,231012
MPDATA,EY,4,,0.76
MPDATA,EZ,4,,0.76
MPDATA,PRXY,4,,0.34
MPDATA,PRYZ,4,,0.34
MPDATA,PRXZ,4,,0.34
MPDATA,GXY,4,,2070
MPDATA,GYZ,4,,2070
MPDATA,GXZ,4,,2070
TBDE,KINH,4,,,
TB,KINH,4,1,2,0
TBTEMP,0
TBPT,,0.0158,3650 !极限应变，抗拉强度
TBPT,,0.01580001,0
KEYOPT,4,1,2
/axlab,x,strain
/axlab,y,stress
tbplot,kinh,4,,,yes
!查看本构
!材料5-混凝土垫块
MP,EX,5,2.0e5
MP,PRXY,5,0.3
!材料6-角钢
mp,ex,6,206e3
mp,prxy,6,0.28
MP,DENS,6,7.85E-9
tb,biso,6,1,2,1
tbdata,2,235
tbdata,1,235
/axlab,x,strain
/axlab,y,stress
tbplot,biso,6,,,yes
!查看本构
!材料7-螺栓
mp,ex,7,206e3
mp,prxy,7,0.3
MP,DENS,7,7.85E-9
tb,biso,7,1,2,1
tbdata,1,640
tbdata,2,640
/axlab,x,strain
/axlab,y,stress
tbplot,biso,7,,,yes
!查看本构
!材料8-单肢箍筋
MP,EX,8,2.1e5
MP,PRXY,8,0.3
TB,BKIN,8,1,2,1
TBDATA,,315
/axlab,x,strain
/axlab,y,stress
tbplot,bkin,8,,,yes
!纤维布厚度
r0=16
h1=250
h2=150
l1=25
l2=12
*SET,RC,0.1
R,13,RC*2,RC*2,RC*2,RC*2
R,14,RC*2.001,RC*2.001,RC*2.001,RC*2.001
R,15,RC*3,RC*3,RC*3,RC*3
R,16,RC*3.001,RC*3.001,RC*3.001,RC*3.001
R,17,RC*2.0001,RC*2.0001,RC*2.0001,RC*2.0001
LOCAL,11,0,1000,1000,1000,90,
LOCAL,12,0,1200,1200,1200,,,90
LOCAL,13,0,1400,1400,1400
!创建实体单元（混凝土）
wpcsys,-1
csys,0
wpcsys,-1
block,0,1000,0,250,0,125
wpoff,0,250,0
blc4,0,0,150,1200,100
vadd,all
wpcsys,-1
wpoff,1000,0,0
blc4,0,0,50,250,125
cm,lz,volu
!创建
wpcsys,-1
wpoff,150,250,0
block,0,200,0,l2,0,100
wpcsys,-1
wpoff,150,250+l2,0
block,0,l2,0,200-l2,0,100
cmsel,u,lz
vadd,all
cm,jg,volu
!定义配筋率
*SET,pi,acos(-1)
RZ=16
RG=6
*SET,v11,0.25*pi*RZ*RZ*1.5/(125*50)!柱拉压区纵筋配筋率
*SET,v122,0.25*pi*RG*RG*4/(200*125)!柱箍筋加密区竖向
*SET,v133,0.25*pi*RG*RG*4/(200*50) !柱箍筋加密区水平向
*SET,v12,0.25*pi*RG*RG*6/(650*125)!柱箍筋非加密区竖向
*SET,v13,0.25*pi*RG*RG*6/(650*50) !柱箍筋非加密区水平向
*SET,v21,0.25*pi*RZ*RZ*1.5/(100*50)!梁拉压区纵筋配筋率
*SET,v222,0.25*pi*RG*RG*5/(250*125)!梁箍筋加密区竖向
*SET,v233,0.25*pi*RG*RG*5/(250*50) !梁箍筋加密区水平向
*SET,v22,0.25*pi*RG*RG*9/(950*125)!梁箍筋非加密区竖向
*SET,v23,0.25*pi*RG*RG*9/(950*50) !梁箍筋非加密区水平向
*SET,v0,0
r,1,3,v11,0,0,2,v133
rmore,0,90,2,v122,90,0
r,2,3,v11,0,0,2,v13
rmore,0,90,2,v12,90,0
r,3,3,v11,0,0,2,v21
rmore,90,0
r,4,2,v122,90,0
r,5,2,v12,90,0
r,6,3,V0,90,0
r,7,2,v222,0,0,2,v233
rmore,0,90,3,v21,90,0
r,8,2,v22,0,0,2,v23
rmore,0,90,3,v21,90,0
r,9,2,v222,0,0
r,10,2,v22,0,0
R,11,28.26
R,12, , , , , , ,
r,101,3,v11,0,0,2,v13 !与R2定义一样啊
rmore,0,90,2,v12,90,0
r,606
r,707
!划分梁柱网格
allsel,all
wpcsys,-1
wpoff,0,0,50
VSBW,ALL
wpcsys,-1
wpoff,0,0,100
VSBW,ALL
wpcsys,-1
wpro,,-90
wpoff,0,0,50
VSBW,ALL
wpoff,0,0,150
VSBW,ALL
wpoff,0,0,50
VSBW,ALL
wpoff,0,0,l2
vsbw,all
wpoff,0,0,50-l2
vsbw,all
wpoff,0,0,50
VSBW,ALL
wpoff,0,0,50
VSBW,ALL
wpoff,0,0,50
VSBW,ALL
*do,i,1,7
wpoff,0,0,100
VSBW,ALL
*enddo
wpcsys,-1
wpoffs,150,,
wprota,0,0,90
vsbw,all
wpof,,,l2
vsbw,all
wpof,,,50-l2
vsbw,all
wpof,,,50
vsbw,all
wpof,,,50
vsbw,all
wpof,,,50
vsbw,all
*do,i,1,5
wpof,,,100
vsbw,all
*enddo
wpcsys,-1
wpro,,,90
wpof,,,50
vsbw,all
wpof,,,50
vsbw,all
!创建
wpcsys,-1
wpof,250,0,50
wprot,,-90
numstr,volu,1000
cyl4,,,r0,,,,h1+l2
vsel,s,volu,,1000
cm,mdd1,volu
wpcsys,-1
wpof,,350,50
wpro,,,90
numstr,volu,2000
cyl4,,,r0,,,,h2+l2
vsel,s,volu,,2000
cm,mdd2,volu
cmsel,a,mdd1
cm,mdd,volu
allsel,all
allsel,all
cmsel,u,mdd
cm,lz,volu
alls
vsbv,lz,mdd
numstr,volu,1
!赋值
allsel,all
vsel,s,loc,y,0,250
vsel,r,loc,x,850,1000
vsel,u,loc,y,50,200
vatt,1,1,1
vsel,s,loc,y,0,250
vsel,r,loc,x,150,850
vsel,u,loc,y,50,200
vatt,1,2,1
vsel,s,loc,y,0,250
vsel,r,loc,x,0,150
vsel,u,loc,y,50,200
vatt,1,3,1
vsel,s,loc,y,50,200
vsel,r,loc,x,850,1000
vatt,1,4,1
vsel,s,loc,y,50,200
vsel,r,loc,x,150,850
vatt,1,5,1
vsel,s,loc,y,50,200
vsel,r,loc,x,0,150
vatt,1,6,1
vsel,s,loc,y,1150,1450
vsel,r,loc,x,100,150
vatt,1,7,1
vsel,s,loc,y,250,1150
vsel,r,loc,x,100,150
vatt,1,8,1
vsel,s,loc,y,1150,1450
vsel,r,loc,x,0,100
vatt,1,9,1
vsel,s,loc,y,250,1150
vsel,r,loc,x,0,100
vatt,1,10,1
ALLSEL,ALL
VSEL,S,LOC,X,1000,1050
VATT,5,12,5
lSEL,S,LOC,X,0
lSEL,R,LOC,Y,50,200
lSEL,R,LOC,Z,0,100
lSEL,U,LOC,Z,0
lSEL,U,LOC,Z,50
LATT,8,11,8
allsel,all
vsel,s,loc,y,250,250+l2
vsel,r,loc,x,150,350
vatt,6,606,6
vsel,s,loc,y,250+l2,450
vsel,r,loc,x,150,150+l2
vatt,6,606,6
!网格
allsel,all
wpcsys,-1
vsel,s,loc,x,200,300
vsel,r,loc,y,200,250
vsel,r,loc,z,0,100
cm,y301,volu
aslv
asel,u,loc,y,200
asel,u,loc,y,250
accat,271,417
accat,376,438
accat,480,453
accat,450,476
cmsel,s,y301
aslv
lsla
lsel,r,radius,,r0
lesize,all,,,4
cmsel,s,y301,volu
aslv
lsla
lsel,u,radius,,r0
lsel,r,loc,y,200
cm,line1,line
lesize,all,,,2
cmsel,s,y301,volu
aslv
lsla
lsel,u,radius,,r0
lsel,r,loc,y,250
cm,line2,line
lesize,all,,,2
cmsel,s,y301,volu
aslv
lsla
lsel,r,tan1,x
lsel,r,tan1,z
lesize,all,50
cmsel,s,y301,volu
aslv
lsla
vplot
mshape,0,3d
mshkey,1
vmesh,all
asel,s,acca
adele,all
lsel,s,lcca
ldele,all
allsel,all
vsel,s,loc,y,200,250
vsel,r,loc,x,200,300
vsel,r,loc,z,100,125
cm,z101,volu
aslv
lsla
lsel,r,tan1,x
lsel,r,tan1,y
lesize,all,l1
vsweep,z101
allsel,all
vsel,s,loc,y,200,250
vsel,r,loc,x,0,1000
vsel,u,loc,x,200,300
cm,z303,volu
aslv
lsla
lsel,r,tan1,y
lsel,r,tan1,z
lesize,all,50
vsweep,z303
allsel,all
vsel,s,loc,x,0,1000
vsel,r,loc,y,0,200
cm,z404,volu
aslv
lsla
lsel,r,tan1,x
lsel,r,tan1,z
lesize,all,50
vsweep,z404
vsel,s,type,,5
cm,z505,volu
aslv
lsla
lsel,r,tan1,x
lsel,r,tan1,z
lesize,all,50
cm,z505,volu
aslv
lsla
lsel,r,tan1,x
lsel,r,tan1,y
lesize,all,25
cm,z505,volu
aslv
lsla
lsel,r,tan1,z
lsel,r,tan1,y
lesize,all,50
vsweep,z505
!梁划分网格
alls
vsel,s,loc,y,300,400
vsel,r,loc,x,100,150
cm,y301,volu
aslv
asel,u,loc,x,100
asel,u,loc,x,150
accat,202,203
accat,200,192
accat,96,164
accat,152,144
cmsel,s,y301,volu
aslv
lsla
lsel,u,radius,,r0
lsel,r,loc,x,100
cm,line1,line
lesize,all,,,2
cmsel,s,y301,volu
aslv
lsla
lsel,u,radius,,r0
lsel,r,loc,x,150
cm,line2,line
lesize,all,,,2
cmsel,s,y301
aslv
lsla
lsel,r,tan1,y
lsel,r,tan1,z
lesize,all,,,1
cmsel,s,y301
mshape,0,3d
mshkey,1
vmesh,all
asel,s,acca
adele,all
lsel,s,lcca
ldele,all
allsel,all
wpcsys,-1
vsel,s,loc,y,250,1450
vsel,r,loc,x,100,150
vsel,u,loc,y,300,400
cm,z101,volu
aslv
lsla
lsel,r,tan1,x
lsel,r,tan1,z
lesize,all,50
vsweep,z101
wpcsys,-1
allsel,all
vsel,s,loc,x,0,100
vsel,r,loc,y,250,1450
cm,z404,volu
aslv
lsla
lsel,r,tan1,z
lsel,r,tan1,y
lesize,all,50
vsweep,z404
!角钢网格划分
allsel,all
wpcsys,-1
vsel,s,loc,x,200,300
vsel,r,loc,y,250,250+l2
vsel,r,loc,z,0,100
cm,y302,volu
aslv
asel,u,loc,y,250+l2
asel,u,loc,y,250
accat,122,54
accat,44,51
accat,365,48
accat,454,273
cmsel,s,y302
aslv
lsla
lsel,r,radius,,r0
lesize,all,,,4
cmsel,s,y302,volu
aslv
lsla
lsel,u,radius,,r0
lsel,r,loc,y,250
cm,line1,line
lesize,all,,,2
cmsel,s,y302,volu
aslv
lsla
lsel,u,radius,,r0
lsel,r,loc,y,250+l2
cm,line2,line
lesize,all,,,2
cmsel,s,y302,volu
aslv
lsla
lsel,r,tan1,x
lsel,r,tan1,z
lesize,all,l2
cmsel,s,y302,volu
aslv
lsla
vplot
mshape,0,3d
mshkey,1
vmesh,all
asel,s,acca
adele,all
lsel,s,lcca
ldele,all
allsel,all
wpcsys,-1
vsel,s,loc,x,150,350
vsel,r,loc,y,250,250+l2
vsel,r,loc,z,0,100
vsel,u,loc,x,200,300
cm,y303,volu
aslv
lsla
lsel,r,tan1,y
lsel,r,tan1,z
lesize,all,50
vsweep,y303
!角钢网格
allsel,all
wpcsys,-1
vsel,s,loc,x,150,150+l2
vsel,r,loc,y,300,400
vsel,r,loc,z,0,100
cm,y304,volu
aslv
asel,u,loc,x,150+l2
asel,u,loc,x,150
accat,180,189
accat,138,185
accat,167,136
accat,171,140
cmsel,s,y304
aslv
lsla
lsel,r,radius,,r0
lesize,all,,,4
cmsel,s,y304,volu
aslv
lsla
lsel,u,radius,,r0
lsel,r,loc,x,150
cm,line1,line
lesize,all,,,2
cmsel,s,y304,volu
aslv
lsla
lsel,u,radius,,r0
lsel,r,loc,x,150+l2
cm,line2,line
lesize,all,,,2
cmsel,s,y304,volu
aslv
lsla
lsel,r,tan1,y
lsel,r,tan1,z
lesize,all,l2
cmsel,s,y304,volu
aslv
lsla
vplot
mshape,0,3d
mshkey,1
vmesh,all
asel,s,acca
adele,all
lsel,s,lcca
ldele,all
wpcsys,-1
allsel,all
vsel,s,loc,x,150,150+l2
vsel,r,loc,y,250+l2,1450
vsel,u,loc,y,300,400
vsel,r,loc,z,0,100
cm,z405,volu
aslv
lsla
lsel,r,tan1,z
lsel,r,tan1,x
lesize,all,50
vsweep,z405
!创建铆钉
allsel,all
wpcsys,-1
wpof,250,0,50
wpro,,-90
numstr,volu,3000
cyl4,,,r0,,,,h1
vsel,s,volu,,3000
cm,md1,volu
aslv
asel,r,loc,y,h1
cm,md_vext,area
esize,,1
vext,md_vext,,,,l2
cm,md2,volu
wpcsys,-1
wpof,,,50
cmsel,s,md2
vsbw,all
wpcsys,-1
wpof,250,0,50
wpro,,,90
vsbw,all
wpcsys,-1
wpro,,-90
wpof,,,50
vsbw,all
cm,md3,volu
cmsel,s,md3
aslv
lsla
lsel,u,length,,200
lsel,u,length,,50
lsel,u,length,,l2
cm,md_line,line
lesize,all,,,4
cmsel,s,md3
aslv
lsla
cmsel,u,md_line
lsel,r,length,,200
lesize,all,,,4
cmsel,s,md3
aslv
lsla
cmsel,u,md_line
lsel,r,length,,50
lesize,all,,,1
cmsel,s,md3
aslv
lsla
cmsel,u,md_line
lsel,r,length,,l2
lesize,all,,,1
cmsel,s,md3
wpcsys,-1
vsel,r,loc,y,0,50
vatt,1,101,1
cmsel,s,md3
vsel,u,loc,y,0,50
vatt,7,707,7
cmsel,s,md3
mshape,0,3d
mshkey,1
vmesh,all
!水平铆钉
wpcsys,-1
allsel,all
wpof,,350,50
wpro,,,90
numstr,volu,4000
cyl4,,,r0,,,,h2
vsel,s,volu,,4000
cm,md1,volu
wpcsys,-1
wpof,100
wpro,,,90
vsbw,all
wpof,,,-50
vsbw,all
wpcsys,-1
wpof,,350,50
vsbw,all
wpro,,-90
vsbw,all
cm,md1,volu
aslv
asel,r,loc,x,h2
cm,md_vext,area
esize,,1
vext,md_vext,,,l2
cm,md2,volu
cmsel,s,md2
aslv
lsla
lsel,u,length,,50
lsel,u,length,,l2
cm,md_line,line
lesize,all,,,4
cmsel,s,md2
aslv
lsla
cmsel,u,md_line
lsel,r,length,,50
lesize,all,,,1
cmsel,s,md2
aslv
lsla
cmsel,u,md_line
lsel,r,length,,l2
lesize,all,,,1
cmsel,s,md2
wpcsys,-1
vsel,r,loc,x,0,50
vatt,1,10,1
cmsel,s,md2
vsel,u,loc,x,0,50
vatt,7,707,7
cmsel,s,md2
mshape,0,3d
mshkey,1
vmesh,all
!纤维布网格划分
ASEL,S,LOC,X,0,150
ASEL,R,LOC,Z,125
aatt,4,13,4,11
ASEL,S,LOC,Y,250,1050
ASEL,R,LOC,Z,100
asel,r,loc,x,0,150
aatt,4,13,4,11
ASEL,S,LOC,Y,250
ASEL,R,LOC,X,0,150
ASEL,R,LOC,Z,100,125
aatt,4,13,4,12
ASEL,S,LOC,X,0,150
ASEL,R,LOC,Y,0
aatt,4,13,4,12
asel,none
k,6500,150,250,0
k,6501,150,250,100
k,6502,150,1050,100
k,6503,150,1050,0
a,6500,6501,6502,6503
wpcsys,-1
wpoff,150,350,50
wpro,,,90
cyl4,0,0,r0
asba,1020,1021
wpcsys,-1
wpro,,-90
wpoff,0,0,l2+250
asbw,all
wpoff,0,0,50-l2
asbw,all
wpoff,0,0,50
aSBW,ALL
wpoff,0,0,50
aSBW,ALL
wpoff,0,0,50
aSBW,ALL
wpcsys,-1
wpoff,0,0,50
aSBW,ALL
aatt,4,14,4,11
asel,none
k,6504,150,250,0
k,6505,150,250,100
k,6506,750,250,100
k,6507,750,250,0
a,6504,6505,6506,6507
wpcsys,-1
wpoff,250,250,50
wpro,,-90
cyl4,0,0,r0
asba,1020,1021
wpcsys,-1
wpof,,,50
asbw,all
wpcsys,-1
wpro,,,90
wpof,,,l2+150
asbw,all
wpof,,,50-l2
asbw,all
wpof,,,50
asbw,all
wpof,,,50
asbw,all
wpof,,,50
asbw,all
aatt,4,14,4,13
vsel,s,type,,1
aslv
asel,r,loc,x,150,350
asel,r,loc,y,250
AGEN,2,ALL,,,,100
asel,r,loc,y,350
aatt,4,15,4,12
AGEN,2,ALL,,,,-100
asel,r,loc,y,350
adele,all
vsel,s,type,,1
aslv
asel,r,loc,x,150,350
asel,r,loc,z,125
AGEN,2,ALL,,,,,100
asel,r,loc,z,225
aatt,4,15,4,11
AGEN,2,ALL,,,,,-100
asel,r,loc,z,225
adele,all
vsel,s,type,,1
aslv
asel,r,loc,x,150,350
asel,r,loc,y,0
AGEN,2,ALL,,,,-100
asel,r,loc,y,-100
aatt,4,15,4,12
AGEN,2,ALL,,,,100
asel,r,loc,y,-100
adele,all
ASEL,S,LOC,X,750,850
ASEL,R,LOC,Y,250
AGEN,2,ALL,,,-100
ASEL,R,LOC,X,650,750
aatt,4,15,4,12
ASEL,S,LOC,X,750,850
ASEL,R,LOC,Y,0,250
ASEL,R,LOC,Z,125
AGEN,2,ALL,,,-100
ASEL,R,LOC,X,650,750
ASEL,U,TYPE,,4
aatt,4,15,4,11
ASEL,S,LOC,X,750,850
ASEL,R,LOC,Y,0
AGEN,2,ALL,,,-100
ASEL,r,LOC,X,650,750
aatt,4,15,4,12
ASEL,S,REAL,,15
ASEL,R,LOC,X,650,750
AGEN,2,ALL,,,-200
vsel,s,type,,1
aslv
asel,r,loc,y,250,450
asel,r,loc,x,0,150
asel,r,loc,z,100
AGEN,2,ALL,,,,,100
asel,r,loc,z,200
aatt,4,16,4,13
AGEN,2,ALL,,,,,-100
asel,r,loc,z,200
adele,all
vsel,s,type,,1
aslv
asel,r,loc,y,250,450
asel,r,loc,x,150
AGEN,2,ALL,,,100
asel,r,loc,x,250
aatt,4,16,4,12
AGEN,2,ALL,,,-100
asel,r,loc,x,250
adele,all
ASEL,S,LOC,Y,1050,1150
asel,r,loc,x,0,150
ASEL,R,LOC,Z,100
AGEN,2,ALL,,,,-100
ASEL,R,LOC,Y,950,1050
aatt,4,16,4,13
ASEL,S,LOC,Y,1050,1150
ASEL,R,LOC,X,150
AGEN,2,ALL,,,,-100
ASEL,R,LOC,Y,950,1050
aatt,4,16,4,12
ASEL,S,REAL,,16
asel,r,loc,y,950,1050
AGEN,3,ALL,,,,-200
asel,none
KSEL,S,LOC,X,0,150
KSEL,R,LOC,Y,0,250
KSEL,R,LOC,Z,125
A,8,87,405,407
A,405,407,488,487
A,488,487,196,197
A,87,69,389,405
A,389,405,487,481
A,487,481,190,196
A,69,5,387,389
A,387,389,481,479
A,481,479,188,190
ASEL,S,LOC,X,0,750
ASEL,R,LOC,Z,125
ASEL,U,TYPE,,4
aatt,4,17,4,13
!纤维布划分网格
wpcsys,-1
asel,s,real,,14
asel,r,loc,x,200,300
cm,a1,area
lsla
lccat,1374,1391
lccat,1386,1384
lccat,1397,1394
lccat,1393,1398
cmsel,s,a1
lsla
lsel,r,radius,,r0
lesize,all,,,4
alls
cmsel,s,a1,area
lsla
lsel,u,radius,,r0
lesize,all,,,2
cmsel,s,a1,area
aplot
mshape,0,2d
mshkey,1
amesh,all
asel,s,real,,15
asel,r,loc,x,200,300
asel,r,loc,y,250
asel,r,loc,z,0,100
cm,a2,area
lsla
lccat,1449,1460
lccat,1459,1478
lccat,1483,1473
lccat,1470,1463
cmsel,s,a2
lsla
lsel,r,radius,,r0
lesize,all,,,4
alls
cmsel,s,a2,area
lsla
lsel,u,radius,,r0
lesize,all,,,2
cmsel,s,a2,area
aplot
mshape,0,2d
mshkey,1
amesh,all
asel,s,real,,15
asel,r,loc,x,200,300
asel,r,loc,y,0
asel,r,loc,z,0,100
cm,a4,area
lsla
lccat,1633,1636
lccat,1635,1654
lccat,1659,1650
lccat,1647,1639
cmsel,s,a4
lsla
lsel,r,radius,,r0
lesize,all,,,4
alls
cmsel,s,a4,area
lsla
lsel,u,radius,,r0
lesize,all,,,2
cmsel,s,a4,area
aplot
mshape,0,2d
mshkey,1
amesh,all
wpcsys,-1
asel,s,real,,14
asel,r,loc,y,300,400
cm,a3,area
lsla
lccat,1351,1366
lccat,1368,1350
lccat,1347,1344
lccat,1348,1357
cmsel,s,a3
lsla
lsel,r,radius,,r0
lesize,all,,,4
alls
cmsel,s,a3,area
lsla
lsel,u,radius,,r0
lesize,all,,,2
cmsel,s,a3,area
aplot
mshape,0,2d
mshkey,1
amesh,all
asel,s,real,,16
asel,r,loc,y,300,400
asel,r,loc,x,150
cm,a5,area
lsla
lccat,1845,1857
lccat,1831,1847
lccat,1851,1840
lccat,1854,1843
cmsel,s,a5
lsla
lsel,r,radius,,r0
lesize,all,,,4
alls
cmsel,s,a5,area
lsla
lsel,u,radius,,r0
lesize,all,,,2
cmsel,s,a5,area
aplot
mshape,0,2d
mshkey,1
amesh,all
asel,s,type,,4
asel,u,area,,a1
asel,u,area,,a2
asel,u,area,,a3
asel,u,area,,a4
asel,u,area,,a5
cm,a6,area
cmsel,s,a6
aplot
lsla,s
LSEL,r,TAN1,y
lsel,r,tan1,x
lesize,all,25
cmsel,s,a6
aplot
lsla,s
LSEL,r,TAN1,y
lsel,r,tan1,z
lesize,all,50
cmsel,s,a6
aplot
lsla,s
LSEL,r,TAN1,x
lsel,r,tan1,z
lesize,all,50
cmsel,s,a6
MSHKEY,1 !映射网格划分
AMESH,ALL !生成纤维布单元
alls
asel,s,mat,,4
arsym,x,all
ALLS
vsym,x,all
alls
allsel,all
esel,s,type,,1
esel,a,type,,4
esel,a,type,,5
esel,a,type,,8
nummrg,elem
nsle
nummrg,node
ksln
nummrg,kp
NUMCMP,ALL
!混凝土—柱上角钢
wpcsys,-1
allsel,all
mp,mu,1,0.2 !定义接触摩擦系数
mp,mu,6,0.2 !定义接触摩擦系数
r,80
real,80
et,11,targe170 !定义接触单元类型
et,61,conta173 !接触面为3D4节点面面元
r,80,,,0.1,0.1,,1.0 !第一个0.15=fkn罚刚度，第二个0.1=ftoln
rmore,,,10.86,,,, !taumax=10.86
rmore,1.84,,,,,, !cohe=1.84
keyopt,61,4,0 !推荐用高斯点即=0
keyopt,61,5,3 !闭合间隙或者减小初始侵入
keyopt,61,7,2 !自动进行时间载荷步检查
keyopt,61,8,1 !避免伪接触
keyopt,61,9,3 !忽略模型引起的初始侵入
keyopt,61,10,2 !程序自行进行刚度修正
keyopt,61,12,2 !不分离接触，允许滑移
VSEL,S,MAT,,6
vsel,r,loc,x,150,350
vsel,r,loc,y,250,250+l2
vsel,r,loc,z,0,100
ASLV,S
asel,r,loc,y,250
cm,target,area
type,11
real,80
nsla,s,1 !选择与面相关的所有节点
esln,s !选择面内的单元
esurf,,top !生成目标单元
!*Generate the contact surface 下面创建接触面
alls
vsel,s,type,,1
aslv
asel,r,loc,x,150,350
asel,r,loc,z,0,100
asel,r,loc,y,250
cm,contact,area
type,61
real,80
nsla,s,1
esln,s,0
esurf,,bottom
eplot
allsel,all
VSEL,S,MAT,,6
vsel,r,loc,x,-350,-150
vsel,r,loc,y,250,250+l2
vsel,r,loc,z,0,100
ASLV,S
asel,r,loc,y,250
cm,target,area
type,11
real,80
nsla,s,1 !选择与面相关的所有节点
esln,s !选择面内的单元
esurf,,top !生成目标单元
!*Generate the contact surface 下面创建接触面
alls
vsel,s,mat,,1
aslv
asel,r,loc,x,-350,-150
asel,r,loc,z,0,100
asel,r,loc,y,250
cm,contact,area
type,61
real,80
nsla,s,1
esln,s,0
esurf,,bottom
eplot
allsel,all
!纤维布—柱上角钢
VSEL,S,MAT,,6
vsel,r,loc,x,150,150+l2
vsel,r,loc,y,250,450
vsel,r,loc,z,0,100
ASLV,S
asel,r,loc,x,150
cm,target,area
type,11
real,80
nsla,s,1 !选择与面相关的所有节点
esln,s !选择面内的单元
esurf,,top !生成目标单元
!*Generate the contact surface 下面创建接触面
alls
vsel,s,mat,,1
aslv
asel,r,loc,y,250,450
asel,r,loc,z,0,100
asel,r,loc,x,150
cm,contact,area
type,61
real,80
nsla,s,1
esln,s,0
esurf,,bottom
eplot
allsel,all
VSEL,S,MAT,,6
vsel,r,loc,x,-150-l2,-150
vsel,r,loc,y,250,450
vsel,r,loc,z,0,100
ASLV,S
asel,r,loc,x,-150
cm,target,area
type,11
real,80
nsla,s,1 !选择与面相关的所有节点
esln,s !选择面内的单元
esurf,,top !生成目标单元
!*Generate the contact surface 下面创建接触面
alls
vsel,s,mat,,1
aslv
asel,r,loc,y,250,450
asel,r,loc,z,0,100
asel,r,loc,x,-150
cm,contact,area
type,61
real,80
nsla,s,1
esln,s,0
esurf,,bottom
eplot
allsel,all
/COM,CONTACT PAIR CREATION - START
!进入电脑设置
!铆钉—混凝土
mp,mu,1,0.15 !定义接触摩擦系数
mp,mu,7,0.15 !定义接触摩擦系数
r,81
real,81
et,12,targe170 !定义接触单元类型
et,91,conta173 !接触面为3D4节点面面元
r,81,,,0.15,0.1,, !第一个0.15=fkn罚刚度，第二个0.1=ftoln
rmore,,,10.86,0.001,,, !taumax=10.86
rmore,1.84,,,,,, !cohe=1.84
keyopt,91,4,0 !推荐用高斯点即=0
keyopt,91,5,3 !闭合间隙或者减小初始侵入
keyopt,91,7,2 !自动进行时间载荷步检查
keyopt,91,8,1 !避免伪接触
keyopt,91,9,3 !
keyopt,91,10,2 !程序自行进行刚度修正
keyopt,91,12,2 !不分离接触，允许滑移
!*Generate the target surface 创建目标面1
vsel,s,mat,,1
aslv,s
aplot
lsla,s
lsel,r,radius,,r0
lsel,r,loc,y,0,250
asll,s,0
asel,u,loc,x,-250,250
asel,u,loc,y,0,50
asel,u,loc,y,200
asel,u,loc,y,250
asel,r,loc,y,50,250
cm,target,area
type,12
real,81
nsla,s,1
esln,s
esurf
!*Generate the contact surface 创建接触面
alls
VSEL,S,MAT,,7
vsel,r,loc,y,0,250
vsel,u,loc,x,-250,250
ASLV,S
asel,u,loc,y,50
asel,u,loc,y,250
asel,u,loc,z,50
asel,u,loc,x,250
asel,u,loc,x,-250
cm,contact,area
nsla,s,1
esln,s
type,91
real,81
esurf
eplot
allsel,all
!*Generate the target surface 创建目标面2
vsel,s,mat,,1
aslv,s
aplot
lsla,s
lsel,r,radius,,r0
lsel,r,loc,y,0,250
asll,s,0
asel,r,loc,x,-250,250
asel,u,loc,y,0,50
asel,u,loc,y,200
asel,u,loc,y,250
asel,r,loc,y,50,250
cm,target,area
type,12
real,81
nsla,s,1
esln,s
esurf
!*Generate the contact surface 创建接触面
alls
VSEL,S,MAT,,7
vsel,r,loc,y,0,250
vsel,r,loc,x,-250,250
ASLV,S
asel,u,loc,y,50
asel,u,loc,y,250
asel,u,loc,z,50
asel,u,loc,x,250
asel,u,loc,x,-250
cm,contact,area
nsla,s,1
esln,s
type,91
real,81
esurf
eplot
!*Generate the target surface 创建目标面3
vsel,s,mat,,1
aslv,s
aplot
lsla,s
lsel,r,radius,,r0
lsel,r,loc,y,270,350
asll,s,0
asel,u,loc,x,-50,50
asel,u,loc,x,150
asel,u,loc,x,-150
asel,u,loc,x,100
asel,u,loc,x,-100
cm,target,area
type,12
real,81
nsla,s,1
esln,s
esurf
!*Generate the contact surface 创建接触面
alls
VSEL,S,MAT,,7
vsel,r,loc,y,270,350
vsel,r,loc,x,-150,150
ASLV,S
asel,u,loc,x,150
asel,u,loc,x,-150
asel,u,loc,x,100
asel,u,loc,x,-100
asel,u,loc,x,50
asel,u,loc,x,-50
asel,u,loc,z,50
asel,u,loc,y,350
cm,contact,area
nsla,s,1
esln,s
type,91
real,81
esurf
eplot
!*Generate the target surface 创建目标面4
vsel,s,mat,,1
aslv,s
aplot
lsla,s
lsel,r,radius,,r0
lsel,r,loc,y,350,430
asll,s,0
asel,u,loc,x,-50,50
asel,u,loc,x,150
asel,u,loc,x,-150
asel,u,loc,x,100
asel,u,loc,x,-100
cm,target,area
type,12
real,81
nsla,s,1
esln,s
esurf
!*Generate the contact surface 创建接触面
alls
VSEL,S,MAT,,7
vsel,r,loc,y,350,430
vsel,r,loc,x,-150,150
ASLV,S
asel,u,loc,x,150
asel,u,loc,x,-150
asel,u,loc,x,100
asel,u,loc,x,-100
asel,u,loc,x,50
asel,u,loc,x,-50
asel,u,loc,z,50
asel,u,loc,y,350
cm,contact,area
nsla,s,1
esln,s
type,91
real,81
esurf
eplot
!合并螺栓与角钢的节点
allsel,all
vsel,s,mat,,7
nslv,s,1
nsel,r,loc,y,250
cm,ls_1,node
nplot
vsel,s,mat,,6
nslv,s,1
nsel,r,loc,y,250
cm,jg_node,node
nplot
cmsel,s,ls_1,node
cmsel,a,jg_node,node
nplot
NUMMRG,NODE,,,,LOW !合并节点
allsel,all
allsel,all
vsel,s,mat,,7
nslv,s,1
nsel,r,loc,y,250+l2
cm,ls_2,node
nplot
vsel,s,mat,,6
nslv,s,1
nsel,r,loc,y,250+l2
cm,jg_node2,node
nplot
cmsel,s,ls_2,node
cmsel,a,jg_node2,node
nplot
NUMMRG,NODE,,,,LOW !合并节点
allsel,all
allsel,all
vsel,s,mat,,7
nslv,s,1
nsel,r,loc,x,150
cm,ls_3,node
nplot
vsel,s,mat,,6
nslv,s,1
nsel,r,loc,x,150
cm,jg_node3,node
nplot
cmsel,s,ls_3,node
cmsel,a,jg_node3,node
nplot
NUMMRG,NODE,,,,LOW !合并节点
allsel,all
allsel,all
vsel,s,mat,,7
nslv,s,1
nsel,r,loc,x,150+l2
cm,ls_4,node
nplot
vsel,s,mat,,6
nslv,s,1
nsel,r,loc,x,150+l2
cm,jg_node4,node
nplot
cmsel,s,ls_4,node
cmsel,a,jg_node4,node
nplot
NUMMRG,NODE,,,,LOW !合并节点
allsel,all
allsel,all
vsel,s,mat,,7
nslv,s,1
nsel,r,loc,x,-150
cm,ls_5,node
nplot
vsel,s,mat,,6
nslv,s,1
nsel,r,loc,x,-150
cm,jg_node5,node
nplot
cmsel,s,ls_5,node
cmsel,a,jg_node5,node
nplot
NUMMRG,NODE,,,,LOW !合并节点
allsel,all
allsel,all
vsel,s,mat,,7
nslv,s,1
nsel,r,loc,x,-150-l2
cm,ls_6,node
nplot
vsel,s,mat,,6
nslv,s,1
nsel,r,loc,x,-150-l2
cm,jg_node6,node
nplot
cmsel,s,ls_6,node
cmsel,a,jg_node6,node
nplot
NUMMRG,NODE,,,,LOW !合并节点
allsel,all
nsel,s,loc,X,-150
nsel,r,loc,y,1300,1450 !从左侧加载
CP,1,UX,ALL !!!2850
FINISH
/SOL
nsel,s,loc,X,1050
D,ALL,UY,,,,,UZ,,,
lsel,s,loc,x,1050
lsel,r,loc,y,200
DL,all,,UX
nsel,s,loc,X,-1050
D,ALL,UY,,,,,UZ,,,
ALLS
ASEL,S,LOC,Z,0$DA,ALL,SYMM
nsel,s,loc,X,-1050
F,ALL,FX,3704
alls
EPLOT
/SOL !262页
anty,static !定义分析类型为静态分析
eqslv,spar,,0, !求解器选择 稀疏矩阵直接法
time,1 !荷载步结束时的时间值
outres,all,all !输出控制
LNSRCH,on !线性搜索
nsubst,20 !子步数和时间步长
NEQIT,50
cnvtol,u,,0.1,0,, !收敛准则
lswrite,1 !生成荷载步文件 第1荷载步
TIME,2
d,3129,ux,21 !对节点2456的UX施加约束，且约束位移值为23
nropt,full,,
nsubst,50
NEQIT,50 !子步数和时间步长
OUTRES,ALL,ALL !输出控制
LNSRCH,on !线性搜索
lswrite,2 !生成荷载步文件 第2荷载步
lssolve,1,2 !多荷载步文件求解
/POST26
nsol,2,3129,u,x !定义变量2为节点2456的X向位移
rforce,3,3129,f,x !定义变量3为节点2456的X向反力
/axlab,x,ux of 3129(mm)
/axlab,y,fx(KN)
prod,2,2,,,,,,1.37
prod,3,3,,,,,,0.002
xvar,2 !定义图形显示的X轴
plvar,3 !图形显示变量
prvar,2,3