一個彈性拱在自重下的變形計算，很簡單的模型，在ansys14.5上用弧長法正常計算，換到2021R2上就算不了了，2019版也試了，也是不收斂。不知道是為什么，請大佬指點。附模型命令流如下：

fini

/clear

/prep7

! 拱軸參數

L  = 7500

H = 1500

L_num= 60

len = L/L_num

MP, DENS, 1, 2.6e-6

mp,ex,1,34500

mp,prxy,1,0.2

MP, DENS, 2, 7.9e-6

mp,ex,2,206000

mp,prxy,2,0.3

! 定義截面--混凝土

SECTYPE, 1, BEAM, CSOLID, Concrete,

SECDATA, 56, 18, 10

! 定義截面--鋼管

SECTYPE, 2, BEAM, CTUBE, Tube,

SECDATA, 56, 60.5, 36

ET, 1, BEAM188

keyopt,1,7,2

keyopt,1,2,1

! 建立拱軸節點

*do, i, 1, L_num+1

x = (i-1)*len

y = -4*H/L/L*x*x+4*H/L*x

n, i, x, y

n, 10000+i, x, -1000, 0

*enddo

! 建立混凝土拱單元

TYPE, 1

mat, 1

secnum, 1

*do, i, 1, L_num

e, i, i+1, 10000+i

*enddo

! 建立鋼管單元

TYPE, 1

mat, 2

secnum, 2

*do, i, 1, L_num

e, i, i+1, 10000+i

*enddo

! 施加約束

d,1,all

d,L_num+1,all

ALLSEL

! 求解

/SOLU

antype, 0             ! 求解類型

time, 1

eqslv,sparse

nsubst,1000, 100000, 10  ! 荷載子步

neqit,500

nlgeom,on

pstres,on

arclen,on

allsel

acel,,10,

outres,all,all

solve

fini

/post26

numvar,200

RFORCE,2,1,F,y, Fz_2

nsol,3,11,u,y,Uy

nsol,4,21,u,y,Uy

nsol,5,31,u,y,Uy

xvar,2

PLVAR,3,4,5

/gropts,view,1