新手求助：关于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以上的值就可以算了
这是为什么呢？
硬件问题吗？
请高手指教啊~~~！！！ 我也出现了这种问题，不知道是什么原因，阿有高手回答一下啊
页:
[1]