请帮助看一下命令流，怎么考虑了材料非线性后就无法运行了？谢谢！

shimly521

finish
/clear,start
/filname,pipeline-5,1
/title,buckling of pipeline
/com,structural
/unit,si
/config,nres,1100000
/PREP7
!定义管参数
!define material parameters
l1=5
l2=10
l3=40
f1=1
f2=1
f3=1
n6=l1/f1+1
n7=n6+(l2-l1)/f2
n8=n7+(l3-l2)/f3
Radius_inner=0.142875
Radius_outer=0.161925
pi=3.1415926
ET,1,PIPE20,
mptemp,,,,,,,,
mptemp,1,0
MPdata,EX,1,1,2.05e11,
MPdata,nuXY,1,1,0.3,
mpdata,dens,1,1,7850,
mpdata,alpx,1,1,1.17e-5,
!应力应变关系
tb,miso,1,1,100,
tbtemp,0
tbpt,,0,0
tbpt,,60e-3/205,60e6
tbpt,,0.1588261135e-2,325.516595e6
tbpt,,0.1608809036e-2,329.5297242e6
tbpt,,0.1613983911e-2,330.5330447e6
tbpt,,0.1629612629e-2,333.5431034e6
tbpt,,0.1640127757e-2,335.5498936e6
tbpt,,0.165072915e-2,337.5567551e6
tbpt,,0.166142583e-2,339.5636915e6
tbpt,,0.1672227704e-2,341.5707067e6
tbpt,,0.1683145637e-2,343.5778051e6
tbpt,,0.1694191543e-2,345.5849915e6
tbpt,,0.171102925e-2,348.5959479e6
tbpt,,0.1728233991e-2,351.6071346e6
tbpt,,0.1745862353e-2,354.6185751e6
tbpt,,0.1757881529e-2,356.6263562e6
tbpt,,0.1770138949e-2,358.634271e6
tbpt,,0.1782657983e-2,360.6423292e6
tbpt,,0.1795464143e-2,362.6505419e6
tbpt,,0.1808585266e-2,364.6589207e6
tbpt,,0.1822051707e-2,366.6674788e6
tbpt,,0.1835896559e-2,368.6762305e6
tbpt,,0.1850155877e-2,370.6851913e6
tbpt,,0.1864868924e-2,372.6943785e6
tbpt,,0.1880078441e-2,374.7038107e6
tbpt,,0.1895830931e-2,376.7135086e6
tbpt,,0.1912176969e-2,378.7234944e6
tbpt,,0.1929171531e-2,380.7337928e6
tbpt,,0.1946874352e-2,382.7444304e6
tbpt,,0.1974899999e-2,385.7610878e6
tbpt,,0.2004909372e-2,388.7786852e6
tbpt,,0.215082966e-2,400.8612577e6
tbpt,,0.2195648414e-2,403.8858184e6
tbpt,,0.2244720339e-2,406.9123801e6
tbpt,,0.2280071796e-2,408.9313306e6
tbpt,,0.2337511228e-2,411.9618408e6
tbpt,,0.2400854839e-2,414.995148e6
tbpt,,0.2470864264e-2,418.0316244e6
tbpt,,0.2548392693e-2,421.0716899e6
tbpt,,0.2634395056e-2,424.1158182e6
tbpt,,0.272993926e-2,427.164543e6
tbpt,,0.2836218464e-2,430.2184649e6
tbpt,,0.2954565169e-2,433.2782596e6
tbpt,,0.3086465209e-2,436.3446865e6
tbpt,,0.3133724822e-2,437.3684471e6
tbpt,,0.3233575148e-2,439.4185983e6
tbpt,,0.3341001076e-2,441.4724989e6
tbpt,,0.345659974e-2,443.5304607e6
tbpt,,0.3581012154e-2,445.5928196e6
tbpt,,0.3714926183e-2,447.6599384e6
tbpt,,0.3859079709e-2,449.7322079e6
tbpt,,0.4014263978e-2,451.8100494e6
tbpt,,0.4181327151e-2,453.8939167e6
tbpt,,0.4361178078e-2,455.9842986e6
tbpt,,0.4554790276e-2,458.0817217e6
tbpt,,0.4987541511e-2,462.3e6
tbpt,,0.5626068e-2,467.6234947e6
tbpt,,0.6392815864e-2,473.014248e6
tbpt,,0.8414692614e-2,484.0560939e6
tbpt,,0.113089197e-1,495.5728226e6
tbpt,,0.2124190040e-1,520.9492489e6
tbpt,,0.407157826e-1,552.0246986e6
tbpt,,0.77507046e-1,594.3244154e6
tbpt,,0.1438569586,658.1897867e6
tbpt,,0.2565384051,762.5445569e6
tbpt,,0.4338004523,941.297657e6

！定义实常数
r,1,2*Radius_outer,Radius_outer-Radius_inner

pi=3.1415926
！定义节点

allsel
n,n6,l1,0,0
n,n7,l2,0,0
n,n8,l3,0,0
fill,n6,n7,
fill,n7,n8,
!缺陷
*do,j,0,n6-2
*afun,rad
x4=j*f1
pi=3.1415926
n,j+1,f1*j,,0.2*(1+cos(pi*x4/l1))/2
*enddo

*do,i,1,n8-1,1
type,1
real,1
mat,1
e,i,i+1,
*enddo
nlgeom,on
/solu
antype,0,
allsel
！定义边界条件
allsel
nsel,s,node,,1
dsym,symm,x
allsel
nsel,s,node,,1,n8-1
d,all,uy
allsel
allsel
d,n8,all
!施加内压
allsel
esel,s,element,,1,n8-1
sfe,all,1,pres,,10400000
allsel
eqsl,sparse
nropt,full,,on
nsubst,5
allsel
solve,
allsel

!加温度
allsel
tref,0
allsel
tunif,300
!设置求解
allsel
solc,off
sstif,off
arclen,on,1,
outres,all,all      
arctrm,u,6,1,uz          !设置收敛准则     
nsubst,25
allsel
solve

/post26
pltime,0,1
nsol,2,1,u,Z,uy_2
prod,3,1,,,,,,30,
/axlab,x,Dis-x(m)
/axlab,y,Pres(℃)
xvar,2
plvar,3
prvar,2
/ui,copy,save,bmp
/post1
set,last
prnsol,u,sum
save,
fini

skylarkcrane

这样的问题还是自己好好钻研一下吧，你看懂了，你会进步很多。
刚看了一点，命令不熟悉，惭愧。

rodge

无法运行，应该是无法收敛是不是？
在引入材料的非线性的时候，收敛这个问题是需要考虑的，但引起不收敛的情况有很多种，需要自己慢慢摸索

shimly521

不是不收敛，是根本不计算