我做了一个SRC的梁的有限元分析

小弟做了一个SRC梁的分析，设置100步，但在95步不收敛了，并自动退出ANSYS，同时产生一个ansabort.db的文件，错误文件中显示的是：
*** FATAL *** CP = 12872.016 TIME= 18:45:38
The supplied nonlinear material properties produced an incorrect
material tangent matrix. Please check your TB input carefully. A
very small yield stress may cause this failure.

************************************************************************
The above error is non-recoverable by ANSYS
ANSYS run terminated by the indicated error
Current data base saved if possible.
************************************************************************
以下是我的材料属性的定义：
finish$/clear$/prep7
!定义单元与材料性质
et,1,solid65 !定义混凝土单元
KEYOPT,1,7,1 !考虑混凝土应力松弛，有助于计算收敛
r,1 ! 定义混凝土的实常数
! 定义混凝土的材料属性
mp,ex,1,13585$mp,prxy,1,0.2$mp,dens,1,2.5e-6 !定义混凝土C30材料的弹性模量、泊松比和密度（kg/mm^3）
fc=14.3$ft=1.43 !混凝土单轴抗压和抗拉强度设计值
tb,concr,1$tbdata,,0.35,0.9,ft,-1 !定义混凝土破坏准则，关闭压碎
tb,miso,1,1,11 !定义混凝土应力应变曲线，用miso模型
tbpt,,0,0
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.002,fc$tbpt,,0.0033,fc*0.85

et,2,solid45 !定义型钢单元
r,2$r,3$r,4 !分别为型钢上、腹、下翼缘实常数，便于分别看上、下翼缘和腹板的结果

mp,ex,2,2.06e5$mp,prxy,2,0.3$mp,dens,2,7.8e-6 !型钢材料的弹性模量和泊松比
tb,bkin,2$tbdata,,300,0 !型钢（Q345-B）的应力应变关系，用BKIN模型

et,3,link8 !定义钢筋单元
mp,ex,3,2e5$mp,prxy,3,0.3$mp,dens,3,7.8e-6 !钢筋材料的弹性模量和泊松比
tb,bkin,3$tbdata,,360,0 !钢筋的应力应变关系，用BKIN模型
pi=acos(-1)$r1=32$r2=12$r3=16 !定义参数π、纵筋直径r1、箍筋直径r2、腰筋直径r3
r,5,0.25*pi*r1*r1 !纵筋的实常数
r,6,0.25*pi*r2*r2 !箍筋的实常数
r,7 !刚性垫板实常数（试算之前材料数据先与型钢一样的）
r,8,0.25*pi*r3*r3 !腰筋的实常数
mp,ex,4,2e6$mp,prxy,4,0.3 !刚性垫板材料的弹性模量和泊松比