ANSYS 分析钢筋混凝土请指教
该模型为钢筋混凝土玻璃纤维加固模型,为线弹性,但分析下来钢筋和纤维布不受力,请各位指教。


/prep7
!单位长度:毫米,密度:千克/mm3  例:牛顿
!沿桥长1200毫米处每隔150mm建立截面上的关键点
*do,i,0,8
k,1+i*16,0,0,i*150
k,2+i*16,30,0,i*150
k,3+i*16,30,40,i*150
k,4+i*16,0,40,i*150
k,5+i*16,0,220,i*150
k,6+i*16,30,220,i*150
k,7+i*16,30,250,i*150
k,8+i*16,0,250,i*150
k,9+i*16,120,250,i*150
k,10+i*16,120,220,i*150
k,11+i*16,150,220,i*150
k,12+i*16,150,250,i*150
k,13+i*16,150,40,i*150
k,14+i*16,120,40,i*150
k,15+i*16,120,0,i*150
k,16+i*16,150,0,i*150
*enddo

!沿梁长1250毫米处的截面上的关键点
*do,i,9,9
k,1+i*16,0,0,1250
k,2+i*16,30,0,1250
k,3+i*16,30,40,1250
k,4+i*16,0,40,1250
k,5+i*16,0,220,1250
k,6+i*16,30,220,1250
k,7+i*16,30,250,1250
k,8+i*16,0,250,1250
k,9+i*16,120,250,1250
k,10+i*16,120,220,1250
k,11+i*16,150,220,1250
k,12+i*16,150,250,1250
k,13+i*16,150,40,1250
k,14+i*16,120,40,1250
k,15+i*16,120,0,1250
k,16+i*16,150,0,1250
*enddo

!将个关键点连成体
*do,i,1,1
v,1+16*(i-1),2+16*(i-1),2+16*i,1+16*i,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i
v,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i
v,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i,8+16*(i-1),7+16*(i-1),7+16*i,8+16*i

v,2+16*(i-1),15+16*(i-1),15+16*i,2+16*i,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i
v,15+16*(i-1),16+16*(i-1),16+16*i,15+16*i,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i
v,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i

v,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i,7+16*(i-1),9+16*(i-1),9+16*i,7+16*i
v,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i,9+16*(i-1),12+16*(i-1),12+16*i,9+16*i
v,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i
*enddo


*do,i,2,2
v,1+16*(i-1),2+16*(i-1),2+16*i,1+16*i,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i
v,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i
v,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i,8+16*(i-1),7+16*(i-1),7+16*i,8+16*i

v,2+16*(i-1),15+16*(i-1),15+16*i,2+16*i,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i
v,15+16*(i-1),16+16*(i-1),16+16*i,15+16*i,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i
v,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i

v,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i,7+16*(i-1),9+16*(i-1),9+16*i,7+16*i
v,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i,9+16*(i-1),12+16*(i-1),12+16*i,9+16*i
v,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i
*enddo


*do,i,3,3
v,1+16*(i-1),2+16*(i-1),2+16*i,1+16*i,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i
v,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i
v,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i,8+16*(i-1),7+16*(i-1),7+16*i,8+16*i

v,2+16*(i-1),15+16*(i-1),15+16*i,2+16*i,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i
v,15+16*(i-1),16+16*(i-1),16+16*i,15+16*i,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i
v,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i

v,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i,7+16*(i-1),9+16*(i-1),9+16*i,7+16*i
v,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i,9+16*(i-1),12+16*(i-1),12+16*i,9+16*i
v,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i
*enddo

 

*do,i,4,4
v,1+16*(i-1),2+16*(i-1),2+16*i,1+16*i,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i
v,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i
v,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i,8+16*(i-1),7+16*(i-1),7+16*i,8+16*i

v,2+16*(i-1),15+16*(i-1),15+16*i,2+16*i,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i
v,15+16*(i-1),16+16*(i-1),16+16*i,15+16*i,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i
v,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i

v,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i,7+16*(i-1),9+16*(i-1),9+16*i,7+16*i
v,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i,9+16*(i-1),12+16*(i-1),12+16*i,9+16*i
v,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i
*enddo


*do,i,5,5
v,1+16*(i-1),2+16*(i-1),2+16*i,1+16*i,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i
v,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i
v,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i,8+16*(i-1),7+16*(i-1),7+16*i,8+16*i

v,2+16*(i-1),15+16*(i-1),15+16*i,2+16*i,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i
v,15+16*(i-1),16+16*(i-1),16+16*i,15+16*i,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i
v,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i

v,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i,7+16*(i-1),9+16*(i-1),9+16*i,7+16*i
v,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i,9+16*(i-1),12+16*(i-1),12+16*i,9+16*i
v,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i
*enddo

 

*do,i,6,6
v,1+16*(i-1),2+16*(i-1),2+16*i,1+16*i,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i
v,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i
v,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i,8+16*(i-1),7+16*(i-1),7+16*i,8+16*i

v,2+16*(i-1),15+16*(i-1),15+16*i,2+16*i,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i
v,15+16*(i-1),16+16*(i-1),16+16*i,15+16*i,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i
v,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i

v,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i,7+16*(i-1),9+16*(i-1),9+16*i,7+16*i
v,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i,9+16*(i-1),12+16*(i-1),12+16*i,9+16*i
v,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i
*enddo


*do,i,7,7
v,1+16*(i-1),2+16*(i-1),2+16*i,1+16*i,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i
v,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i
v,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i,8+16*(i-1),7+16*(i-1),7+16*i,8+16*i

v,2+16*(i-1),15+16*(i-1),15+16*i,2+16*i,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i
v,15+16*(i-1),16+16*(i-1),16+16*i,15+16*i,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i
v,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i

v,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i,7+16*(i-1),9+16*(i-1),9+16*i,7+16*i
v,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i,9+16*(i-1),12+16*(i-1),12+16*i,9+16*i
v,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i
*enddo


*do,i,8,8
v,1+16*(i-1),2+16*(i-1),2+16*i,1+16*i,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i
v,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i
v,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i,8+16*(i-1),7+16*(i-1),7+16*i,8+16*i

v,2+16*(i-1),15+16*(i-1),15+16*i,2+16*i,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i
v,15+16*(i-1),16+16*(i-1),16+16*i,15+16*i,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i
v,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i

v,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i,7+16*(i-1),9+16*(i-1),9+16*i,7+16*i
v,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i,9+16*(i-1),12+16*(i-1),12+16*i,9+16*i
v,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i
*enddo


*do,i,9,9
v,1+16*(i-1),2+16*(i-1),2+16*i,1+16*i,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i
v,4+16*(i-1),3+16*(i-1),3+16*i,4+16*i,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i
v,5+16*(i-1),6+16*(i-1),6+16*i,5+16*i,8+16*(i-1),7+16*(i-1),7+16*i,8+16*i

v,2+16*(i-1),15+16*(i-1),15+16*i,2+16*i,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i
v,15+16*(i-1),16+16*(i-1),16+16*i,15+16*i,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i
v,3+16*(i-1),14+16*(i-1),14+16*i,3+16*i,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i

v,6+16*(i-1),10+16*(i-1),10+16*i,6+16*i,7+16*(i-1),9+16*(i-1),9+16*i,7+16*i
v,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i,9+16*(i-1),12+16*(i-1),12+16*i,9+16*i
v,14+16*(i-1),13+16*(i-1),13+16*i,14+16*i,10+16*(i-1),11+16*(i-1),11+16*i,10+16*i
*enddo
allsel,all


!选择主钢筋所在处的线(共4条线)
VSEL,S,LOC,Z,1200,0!-1200
VSEL,R,LOC,X,30,120
VSEL,R,LOC,Y,40,220
ASLV,S
ASEL,R,LOC,Y,220
LSLA,S
LSEL,R,LOC,X,120
cm,LZ1,line
ALLSEL,ALL

VSEL,S,LOC,Z,1200,0!-1200
VSEL,R,LOC,X,30,120
VSEL,R,LOC,Y,40,220
ASLV,S
ASEL,R,LOC,Y,220
LSLA,S
LSEL,R,LOC,X,30
cm,LZ2,line
ALLSEL,ALL

VSEL,S,LOC,Z,1200,0!-1200
VSEL,R,LOC,X,30,120
VSEL,R,LOC,Y,40,220
ASLV,S
ASEL,R,LOC,Y,40
LSLA,S
LSEL,R,LOC,X,30
cm,LZ3,line
ALLSEL,ALL


VSEL,S,LOC,Z,1200,0!-1200
VSEL,R,LOC,X,30,120
VSEL,R,LOC,Y,40,220
ASLV,S
ASEL,R,LOC,Y,40
LSLA,S
LSEL,R,LOC,X,120
cm,LZ4,line
ALLSEL,ALL

cmgrp,lz,lZ1,lZ2,lZ3,lZ4


!选出箍筋所在的线

VSEL,S,LOC,Z,1200,0!-1200
VSEL,R,LOC,X,30,120
VSEL,R,LOC,Y,40,220
ASLV,S
ASEL,R,LOC,X,120
CM,A1,AREA
ALLSEL

VSEL,S,LOC,Z,1200,0!-1200
VSEL,R,LOC,X,30,120
VSEL,R,LOC,Y,40,220
ASLV,S
ASEL,R,LOC,X,30
CM,A2,AREA
ALLSEL

VSEL,S,LOC,Z,1200,0!-1200
VSEL,R,LOC,X,30,120
VSEL,R,LOC,Y,40,220
ASLV,S
ASEL,R,LOC,Y,40
CM,A3,AREA
ALLSEL

VSEL,S,LOC,Z,1200,0!-1200
VSEL,R,LOC,X,30,120
VSEL,R,LOC,Y,40,220
ASLV,S
ASEL,R,LOC,Y,220
CM,A4,AREA
ALLSEL

CMSEL,S,A1,AREA
CMSEL,A,A2,AREA
CMSEL,A,A3,AREA
CMSEL,A,A4,AREA
LSLA,S

CM,LG,LINE

!选出所有钢筋
cmgrp,lzong,lz,lg


!创建纤维布
VSEL,S,LOC,Z,1200,0!-1200
VSEL,R,LOC,X,0,150
VSEL,R,LOC,Y,0,40
ASLV,S
ASEL,R,LOC,Y,0
cm,fa,area
ALLSEL


!定义钢筋参数
rd11=12.0 !主钢筋直径
rd12=8.0!架立钢筋直径
rd13=6.5 !箍筋直径

!定义单元类型 
et,1,solid65
et,2,link8
ET,3,shell63


!定义混凝土参数
MPTEMP,,,,,,,,              
MPTEMP,1,0                  
                            
MPDATA,EX,1,,3.0E4      
MPDATA,PRXY,1,,0.3         
                            
!TB,MELA,1,1,5               
!TBTEMP,0                    
!TBPT,,0.0005,15    
!TBPT,,0.001,24    
!TBPT,,0.002,30    
!TBPT,,0.0025,29  
!TBPT,,0.0038,22  
!TBPT,,0.002347,30.87036536

              
MPTEMP,1,0                  
MPDATA,ALPX,1,,1e-5         
MPTEMP,1,0                  
MPDATA,DENS,1,,2.5e-3        
!TB,CONC,1,1,9,              
!TBTEMP,0                    
!TBDATA,,0.25,0.95,3.64,-1,, 
!TBDATA,,,,,,,


!定义钢筋实常数
mp,ex,2,2.1e5
mp,prxy,2,0.16

r,11,3.14*0.25*rd11*rd11
r,12,3.14*0.25*rd12*rd12
r,13,3.14*0.25*rd13*rd13


!定义纤维布参数
mp,ex,3,7.24e4
mp,ey,3,1.25e5
mp,prxy,3,0.26
mp,gxy,3,1.52e1

R,3,1.3,0, 


!对各条线的划分单元
*DO,I,0,9
KSEL,S,KP,,1+I*16,2+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,15+I*16,16+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO


*DO,I,0,9
KSEL,S,KP,,3+I*16,4+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,13+I*16,14+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,5+I*16,6+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,10+I*16,11+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,7+I*16,8+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,9+I*16,12+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,1+I*16,4+I*16
LSLK,S,1
lESIZE,all,20
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,2+I*16,3+I*16
LSLK,S,1
lESIZE,all,20
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,14+I*16,15+I*16
LSLK,S,1
lESIZE,all,20
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,13+I*16,16+I*16
LSLK,S,1
lESIZE,all,20
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,5+I*16,8+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,6+I*16,7+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,9+I*16,10+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,11+I*16,12+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,4+I*16,5+I*16
LSLK,S,1
lESIZE,all,20
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,3+I*16,6+I*16
LSLK,S,1
lESIZE,all,20
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,10+I*16,14+I*16
LSLK,S,1
lESIZE,all,20
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,11+I*16,13+I*16
LSLK,S,1
lESIZE,all,20
ALLSEL
*ENDDO


*DO,I,0,9
KSEL,S,KP,,2+I*16,15+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,3+I*16,14+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,6+I*16,10+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,9
KSEL,S,KP,,7+I*16,9+I*16
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO

*DO,I,0,15
*DO,J,1,8
KSEL,S,KP,,1+I+16*(J-1),17+I+16*(J-1)
LSLK,S,1
lESIZE,all,15
ALLSEL
*ENDDO
ALLSEL
*ENDDO

*DO,I,16,16
*DO,J,1,8
KSEL,S,KP,,1+I+16*(J-1),17+I+16*(J-1)
LSLK,S,1
lESIZE,all,25
ALLSEL
*ENDDO
ALLSEL
*ENDDO

 

MSHAPE,0,3d
MSHKEY,1

!建立混凝土单元
type,1
mat,1

vsel,all
vmesh,all

!建立钢筋单元
type,2
mat,2
real,11
allsel,all
cmsel,s,lz1,line
cmsel,a,lz2,line
!LCLEAR,ALL
lmesh,all
allsel

type,2
mat,2
real,12
allsel,all
cmsel,s,lz3,line
cmsel,a,lz4,line
!LCLEAR,ALL
lmesh,all
allsel

type,2
mat,2
real,13
allsel,all
cmsel,s,lg,line
!LCLEAR,ALL
lmesh,all
allsel
!建立纤维布单元
type,3
mat,3
real,3
allsel,all
cmsel,s,fa
amesh,fa
allsel,all


!创建梁的另一半的对称点和单元
NSYM,Z,20000,ALL
ESYM,,20000,ALL

NUMMRG,NODE, , , ,LOW       
NUMMRG,KP, , , ,LOW 

!对梁两侧支点进行约束
NSEL,S,LOC,Z,1110,1170
NSEL,R,LOC,Y,0
D,ALL,UY

NSEL,S,LOC,Z,-1110,-1170
NSEL,R,LOC,Y,0
D,ALL,UY


!对梁施加对称条件
NSEL,S,LOC,Z,0
DSYM,SYMM,Z
ALLSEL
NSEL,S,LOC,X,75
DSYM,SYMM,X
ALLSEL
/SOLU
ANTYPE,0             
allsel,all
!--------------定义非线性分析选项--------------------
cnvtol,f,,0.05,2 !定义收敛条件,使用缺省的value
!nsubst,10 !定义子步数
outres,all,all !输出每一子步的结果
autots,1 !打开自动时间步控制
lnsrch,1 !打开线性搜索
ncnv,2 !如果不收敛时结束而不退出
!neqit,50 !每一子步中方程的迭代次数限值
pred,on !打开预测器
time,1 !定义时间
ACEL,,9.8
solve !求解