Ansys瞬态分析的收敛问题

用 ansys做瞬态分析，计算不收敛。!C型钢阻尼元件Ansys模拟
FINISH
/CLE
/PREP7
!几何参数，C型钢所在曲线为椭圆
A_in=700$B_in=370 !内椭圆的长轴半径以及短轴半径
A_out=780$B_out=426 !外椭圆的长轴半径以及短轴半径
A0=(A_in+A_out)/2$B0=(B_in+B_out)/2 !轴线椭圆的长轴半径以及短轴半径
D_in=36 !螺栓孔的直径
R_in=D_in/2 !螺栓孔的半径
R_out=36 !C型钢两端铰支处的外半径
Rd=50 !倒角
H=40 !C型钢的厚度
D1=273.8 !下螺栓孔中心距C型钢所在椭圆中的水平距离
D2=109.7 !两螺栓孔中心之间的水平距离
!实体建模
k,1,,B_out
k,2,,B_in
k,3,,-B_in
k,4,,-B_out
k,5,,B0
k,6,,-B0
CSWPLA,11,1,B_out/A_out,1,
L,1,4
CSWPLA,12,1,B_in/A_in,1,
L,2,3
CSWPLA,13,1,B0/A0,1,
L,5,6
CSYS,0
LSYMM,X,all
LDELE,1,3
NUMMRG,ALL
k,7,-(D1+D2)
k,8,-(D1+D2),B_out
l,7,8
LOVLAP,1,6
LDELE,2
LDELE,7
LDELE,8
k,10,-D1
k,11,-D1,-B_out
l,10,11
LOVLAP,1,3
LDELE,2
LDELE,7
LDELE,8
KDELE,5,8,1
KDELE,10,11
NUMMRG,ALL
NUMCMP,ALL
CIRCLE,5,R_out
CIRCLE,6,R_out
LOVLAP,1,5
LDELE,12
LOVLAP,2,6
LDELE,1
LOVLAP,5,9
LDELE,1
LOVLAP,13,10
LDELE,1
LDELE,3
A,16,15,18,17
AL,4,7,16,12,15,14
AL,2,8,11,9,5,6
LFILLT,14,18,Rd
LFILLT,16,17,Rd
LFILLT,9,18,Rd
LFILLT,2,17,Rd
AL,10,13,19
AL,20,21,22
AL,23,24,25
AL,26,27,28
AADD,ALL
CIRCLE,5,R_in
CIRCLE,6,R_in
AL,1,3,5,6
AL,12,13,15,19
ASBA,8,1
ASBA,3,2
lcomb,10,14
lcomb,20,16
lcomb,2,26
lcomb,9,23
NUMMRG,ALL
NUMCMP,ALL
CSYS,1
AGEN,2,all,,,,4.455,,,,1
KGEN,2,5,6,,,4.455,,,,1
CSYS,0
KDELE,1,4
NUMMRG,ALL
NUMCMP,ALL
!网格划分
ET,1,PLANE82
KEYOPT,1,3,3
R,1,H
kwpave,16
wprota,,,90
wprota,,-30
asbw,1
kwpave,17
wprota,,40
asbw,3
lcomb,22,16
LESIZE,21,,,80
LESIZE,16,,,80
LESIZE,17,,,10
LESIZE,20,,,10
MSHAPE,0,2D
MSHKEY,1
AMESH,1
lcomb,19,10
LESIZE,10,,,10
LESIZE,4,,,10
LESIZE,7,,,10
LESIZE,15,,,4
LESIZE,1,,,10
LESIZE,6,,,10
LESIZE,5,,,10
LESIZE,3,,,10
MSHAPE,0,2D
MSHKEY,0
AMESH,2
LESIZE,9,,,8
LESIZE,11,,,8
LESIZE,8,,,8
LESIZE,23,,,4
LESIZE,12,,,10
LESIZE,13,,,10
LESIZE,14,,,10
LESIZE,18,,,10
MSHAPE,0,2D
MSHKEY,0
AMESH,4
NUMMRG,ALL
NUMCMP,ALL
NPLO
ET,2,MPC184
KEYOPT,2,1,1
TYPE,2
NKPT,8000,1
NKPT,8001,2
*DO,I,2649,2728,1
E,8000,I
*ENDDO
*DO,I,3228,3307,1
E,8001,I
*ENDDO
!材料属性，Q345B钢材
MP,EX,1,2.1e5
MP,PRXY,1,0.3
tb,bkin,1,1,2,1
tbdata,1,340,6100,,,,
TBPLOT
finish
save
/solu
antype,4
trnopt,full
D,8000,UX,0
D,8000,UY,0
D,8001,UX,0
DT=0.1
*DO,I,1,40
D,8001,UY,100*SIN(2*3.14*DT*I)
TIME,DT*I
NSUBST,1,,,1
KBC,0
OUTRES,ALL,ALL
SOLVE
*ENDDO
TIME,40*DT+2
NSUBST,20,,,1
KBC,0
SOLVE
FINISH
/POST26
NSOL,2,8001,U,Y
RFORCE,3,8001,F,Y
XVAR,2
PLVAR,3