新手求助：关于ANSYS非线性分析时出现的问题，急！！！

忘记歌词

命令流如下：
!Analysis on a simple supported reinforced concrete beam
!分离式模型，1/4模型分析，力加载，位移收敛准则，收敛误差设为3%，关闭压碎
!Units: N, mm, MPa (N/mm2)
finish
/clear
/config,nres,2000  !maximum number of results sets (substeps) allowed on the result file, set to 2000
!
/prep7
!1. Define element types, material properties, real constants-----------------------------
!---set element types---
et,1,solid65  !concrete
keyopt,1,1,0  !Include extra displacement shapes
keyopt,1,7,1  !Include tensile stress relaxation after cracking to help convergence
et,2,link8  !steel/rebars
!
!---set material properties---
!(a)concrete
mp,ex,1,16588
mp,prxy,1,0.2
!constitutive relationship of concrete
fc=14.3
ft=1.43
tb,concr,1
tbdata,,0.5,0.95,ft,-1,,  !βt, βc, ft, no crush,
tbdata,,,,0.6  !Tc=0.6
tb,miso,1,,19
tbpt,,0.0002,fc*0.232
tbpt,,0.0004,fc*0.432
tbpt,,0.0006,fc*0.602
tbpt,,0.0008,fc*0.741
tbpt,,0.0010,fc*0.850
tbpt,,0.0012,fc*0.930
tbpt,,0.0014,fc*0.979
tbpt,,0.00164,fc*1.000
tbpt,,0.0018,fc*0.988
tbpt,,0.0020,fc*0.949
tbpt,,0.0022,fc*0.894
tbpt,,0.0024,fc*0.834
tbpt,,0.0026,fc*0.773
tbpt,,0.0028,fc*0.715
tbpt,,0.0030,fc*0.662
tbpt,,0.0032,fc*0.613
tbpt,,0.0034,fc*0.570
tbpt,,0.0036,fc*0.531
tbpt,,0.003772,fc*0.5
tbplot,miso,1
!----Steel----
mp,ex,2,2.0e5
mp,prxy,2,0.3
tb,biso,2
tbdata,,300,0
tbplot,biso,2
!-------Real constants--------
pi=acos(-1)
r,1,pi*22*22/4  !rebar at the tensile region corner
r,2,pi*22*22/4/2  !rebar on middle of the tensile region（on the symmetric plane）
r,3,pi*10*10/4  !rebar at the compressional region corner; stirrup
r,4,pi*10*10/4/2  !stirrup on the symmetric plane
r,5,,,,,, !concrete
rmore,,,,,,
!
!2. Create geometry model--------------------------------------
/view,1,1,1,1
/replot
blc4,,,200/2,300,3000/2
*do,i,1,14
wpoffs,,,100
vsbw,all
*enddo
wpcsys,-1
wpoffs,,,50
vsbw,all
wpcsys,-1
wprota,,-90
wpoffs,,,30
vsbw,all
wpoffs,,,240
vsbw,all
wpcsys,-1
wpoffs,30
wprota,,,90
vsbw,all
wpcsys,-1
!3. Mesh the steel-------------------------
!create components and assign the label of element types/real constants/material properties
elemsize=50  !set element size
lsel,s,loc,x,30
lsel,r,loc,y,30
lplot
cm,As-t,line !As(corner)
latt,2,1,2
lesize,all,elemsize
!--------
lsel,s,loc,x,100
lsel,r,loc,y,30
lplot
cm,As-tm,line !As(middle, on the symmetry plane)
latt,2,2,2
lesize,all,elemsize
!--------
lsel,s,loc,x,30
lsel,r,loc,y,270
lplot
cm,As-c,line !As'
latt,2,3,2
lesize,all,elemsize
!---------
lsel,s,tan1,z
lsel,r,loc,y,30,270
lsel,r,loc,x,30,95
lsel,u,loc,z,50
lplot
cm,Asv,line  !stirrups (general)
latt,2,3,2
lesize,all,elemsize
!----------
lsel,s,loc,z,0
lsel,r,loc,y,30,270
lsel,r,loc,x,30,95
lplot
cm,Asv1,line  !stirrup(on the symmetry plane)
latt,2,4,2
lesize,all,elemsize
!
lsel,all
cmsel,s,As-t
cmsel,a,As-tm
cmsel,a,As-c
cmsel,a,Asv
cmsel,a,Asv1
cm,allrebars,line  !all rebars
!
lmesh,all
lsel,all
!
/eshape,1
/number,1
/pnum,real,1
eplot
!
!4. Mesh concrete-----------------------------------
allsel,all
vatt,1,5,1
mshkey,1
esize,elemsize
vmesh,all
allsel,all
nummrg,all
!5. Apply loads and restraints----------------------
!simple support restraints
lsel,s,loc,y,0
lsel,r,loc,z,1400
lplot
dl,all,,ux,0
dl,all,,uy,0
!
allsel,all
/replot
!
asel,s,loc,z,0
da,all,symm
!
asel,s,loc,x,100
da,all,symm
!
p=115000  !concentrating load
qe=p/200/100  !equivalent distributing load
!
asel,s,loc,z,0,50
asel,r,loc,y,300
sfa,all,1,pres,qe
allsel,all
/replot
finish
!
!6. Set solution options and solve
/solu
antype,0
nsubst,250
outres,all,all
autots,on
neqit,150
cnvtol,u,,0.03
gplot



以上命令流是老师给的，其他人都可以正常计算没问题，但是我每次算到一半ANSYS就自动退出了！
错误日志里的描述是：

*** FATAL ***                           CP =      33.328   TIME= 22:47:23
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.


马上要交作业了！！！
请各位大侠指教！！！

忘记歌词

是与切线模量的设置有关吗？

但是为什么别人计算都没问题呢？

忘记歌词

我把屈服应力改成310以上的值就可以算了

这是为什么呢？

硬件问题吗？

请高手指教啊~~~！！！

civilstar

我也出现了这种问题，不知道是什么原因，阿有高手回答一下啊