网壳 模态分析 出错（怀疑是模型问题）求助！！

进行模态分析的时候，算三次都显示
不同的单元，计算长度和实际长度不符，但再次SOLVE的时候就可以算出结果。
不知道是不是模型除了问题，这样的话后面的动力计算也没法算下去。
小白，求大神指教！！
下面是模型的命令流，是用师兄的修改的，不知道是不是哪里有错误
!建立鞍型网壳的有限元模型
!***************************
/PREP7
!(1)设置模型几何常数
B_span=4 !长度
SPAN=B_span*2**(1/2) !跨度
DIVSPAN=4 !对角线等分4份
ratio_HS=1/5
!高跨比1/5
AHIGH=0.5*SPAN*ratio_HS !拱高
SHIGH=0.5*SPAN*ratio_HS !跨高
!设定双曲抛物面方程中的P、Q参数
CONSP=SPAN**2/(4*SHIGH)
CONSQ=SPAN**2/(4*AHIGH)
!(2)定义关键点
PI=3.1415926735
*DO,I,1,DIVSPAN+1,1
*DO,J,1,DIVSPAN+1,1
X_COORD=-SPAN*SIN(PI/4)/2+SPAN*SIN(PI/4)/DIVSPAN*(J-1)
Y_COORD=-SPAN*SIN(PI/4)/2+SPAN*SIN(PI/4)/DIVSPAN*(I-1)
Z_COORD=(X_COORD+Y_COORD)**2/(2*CONSP)-(X_COORD-Y_COORD)**2/(2*CONSQ)
K,(DIVSPAN+1)*(I-1)+J,X_COORD,Y_COORD,Z_COORD
*ENDDO
*ENDDO
!建杆单元
!连接内部横向线
*DO,I,2,DIVSPAN
*DO,J,1,DIVSPAN
L,(DIVSPAN+1)*(I-1)+J,(DIVSPAN+1)*(I-1)+J+1
*ENDDO
*ENDDO
!连接内部纵向线
*DO,I,2,DIVSPAN
*DO,J,1,DIVSPAN
L,I+(DIVSPAN+1)*(J-1),I+(DIVSPAN+1)*J
*ENDDO
*ENDDO
!连接斜向线
*DO,I,1,DIVSPAN,
*DO,J,1,I,1
L,(DIVSPAN+1)-I+(DIVSPAN+1+1)*(J-1),(DIVSPAN+1)-I+(DIVSPAN+1+1)*J
*ENDDO
*ENDDO
*DO,I,1,(DIVSPAN-1),1
*DO,J,1,DIVSPAN-I
L,(DIVSPAN+1)*(I-1)+(DIVSPAN+1+1)*J,(DIVSPAN+1)*(I-1)+(DIVSPAN+1+1)*(J+1)
*ENDDO
*ENDDO
!画边框线
*DO,I,1,DIVSPAN
L,I,I+1
*ENDDO
*DO,I,1,DIVSPAN
L,(DIVSPAN+1)*I,(DIVSPAN+1)*(I+1)
*ENDDO
*DO,I,1,DIVSPAN
L,(DIVSPAN+1)**2-(I-1),(DIVSPAN+1)**2-I
*ENDDO
*DO,I,1,DIVSPAN
L,DIVSPAN*(DIVSPAN+1)+1-(DIVSPAN+1)*(I-1),DIVSPAN*(DIVSPAN+1)+1-(DIVSPAN+1)*I
*ENDDO
!(3)定义单元类型及是常数
!定义杆件属性
!定义常数π
PI=3.1415926735
!定义边梁BEAM4几何参数
B=0.25
H=0.25
T=0.008
BM_AREA=B*H-(B-T)*(H-T)
IYY=B*H**3/12-(B-T)*(H-T)**3/12
IZZ=B**3*H/12-(B-T)**3*(H-T)/12
!内部网壳PIPE20几何参数Φ16×1.5
OUTDIA1=25*1.0E-3
THCK1=2.5*1.0E-3
OUTDIA2=25*1.0E-3
THCK2=2.5*1.0E-3
!设置KP单元类型MASS21
ET,1,MASS21
KEYOPT,1,1,0
KEYOPT,1,2,0
KEYOPT,1,3,2
!设置LINE单元类型PIPE20
ET,2,PIPE20
!设置边梁单元类型BEAM4:B×H=0.25×0.25×0.008
ET,3,BEAM4
!设置实常数
!静力计算的内力组合
P=175 !已知的屋面荷载标准值(kg/m2)
PG=(P*9.81+0.5*500)/9.81 !重力荷载代表值
!设置节点的实常数
P0=1/3*PG*0.5*(B_span/4)**2 !单块三角形所承担的屋面荷载
P1=P0 !单块一块三角形屋面的节点所承受荷载
P2=2*(P0) !单块两块三角形屋面的节点所承受荷载
P3=3*(P0) !单块三块三角形屋面的节点所承受荷载
P6=6*(P0) !单块六块三角形屋面的节点所承受荷载
*DO,I,1,1 !最内侧横向边框节点实常数设置
*DO,J,1,DIVSPAN+1
*IF,J,EQ,1,THEN
R,5*(I-1)+J,P2
*ELSEIF,J,GE,2,AND,J,LE,DIVSPAN
R,5*(I-1)+J,P3
*ELSEIF,J,EQ,DIVSPAN+1
R,5*(I-1)+J,P1
*ENDIF
*ENDDO
*ENDDO
*DO,I,2,DIVSPAN !中间部分横向边框节点实常数设置
*DO,J,1,DIVSPAN+1
*IF,J,GE,2,AND,J,LE,DIVSPAN,THEN
R,5*(I-1)+J,P6
*ELSE
R,5*(I-1)+J,P3
*ENDIF
*ENDDO
*ENDDO
*DO,I,DIVSPAN+1,DIVSPAN+1 !最外侧横向边框节点实常数设置
*DO,J,1,DIVSPAN+1
*IF,J,EQ,1,THEN
R,5*(I-1)+J,P1
*ELSEIF,J,GE,2,AND,J,LE,DIVSPAN
R,5*(I-1)+J,P3
*ELSEIF,J,EQ,DIVSPAN+1
R,5*(I-1)+J,P2
*ENDIF
*ENDDO
*ENDDO
R,1+(DIVSPAN+1)**2,OUTDIA1,THCK1, !设置内部直纹杆件实常数
R,2+(DIVSPAN+1)**2,OUTDIA2,THCK2, !设置内部斜杆实常数
R,3+(DIVSPAN+1)**2,BM_AREA,IZZ,IYY,B,H !设置边梁实常数
!定义梁截面类型和几何数据
SECTYPE,1,BEAM,HREC !定义箱形截面
SECDATA,0.25,0.25,0.008,0.008,0.008,0.008 !输入箱形截面几何参数
!设置材料属性
MPTEMP,1,0
MPDATA,EX,1,,2.06E11
MPDATA,PRXY,1,,0.3
MPDATA,DENS,1,,7850
TB,BKIN,1,1,2,1
TBTEMP,0
TBDATA,,235E6,0,,,,
!(5)几何模型的单元属性与网格划分
!关键点单元属性赋值
*DO,I,1,(DIVSPAN+1)**2,1
KSEL,S,KP,,I
KATT,1,I,1,0
KSEL,ALL
*ENDDO
!关键点网格划分
KESIZE,ALL,0
KMESH,ALL
GPLOT
!赋予内部直纹杆件单元类型、数值及网格划分
!其中(DIVSPAN-1)*DIVSPAN为纵线(或横线)总数
LSEL,S,LINE,,1,(DIVSPAN-1)*DIVSPAN*2,,
LATT,1,1+(DIVSPAN+1)**2,2,0
LESIZE,ALL,,,3,,,,,1
LMESH,ALL
LSEL,ALL
GPLOT
!赋予内部斜杆件单元类型、数值及网格划分
!其中(DIVSPAN-1)*DIVSPAN为纵线(或横线)总数，DIVSPAN**2为斜向线总数
LSEL,S,LINE,,(DIVSPAN-1)*DIVSPAN*2+1,(DIVSPAN-1)*DIVSPAN*2+DIVSPAN**2,,
LATT,1,2+(DIVSPAN+1)**2,2,0
LESIZE,ALL,,,3,,,,,1
LMESH,ALL
LSEL,ALL
GPLOT
!赋予边梁单元类型、数值及网格划分
LSEL,S,LINE,,(DIVSPAN-1)*DIVSPAN*2+DIVSPAN**2+1,(DIVSPAN-1)*DIVSPAN*2+DIVSPAN**2+DIVSPAN*4,,
LATT,1,3+(DIVSPAN+1)**2,3,0
LESIZE,ALL,,,3,,,,,1
LMESH,ALL
LSEL,ALL
GPLOT
/REPLOT
!施加约束
DK,1,ALL,0
DK,5,ALL,0
DK,21,ALL,0
DK,25,ALL,0
!(6)保存模型文件
SAVE,MODE,DB,,ALL