索单元找形后是处于悬链线状态吗？

fflyer

最近使用LINK 10 索单元，对两端固定的悬挂导线进行找形。
我想请教下，在ANSYS中，导线找形是基于什么理论的？
因为我找到的初始位形跟悬链线理论计算的结果差距比较大，望大家指点迷津！
使用的命令流如下：

FINISH
/CLEAR
/FILNAME,EX820
/PREP7
! 1.定义儿何参数和荷载参数等，定义单元类型和材料性质
L0=1055  
AREA=6.336E-4
EM=10.094E10                              !定义几何参数、面积、弹性模量
Q0=26.999                           !定义索单位重量（N/m）和集中荷载（N）
H0=64280.2                               !定义自重作用下的水平张力(已知)
ERR0=1/1000
ENUM=100
ISTRAN=1.0051E-3                         !定义迭代条件，单元数目、初应变
ET,1,LINK10
R,1,AREA,ISTRAN                            !定义单元和实常数
MP,EX,1,EM
MP,PRXY，1，0.3                                !定义材料性质
MP,DENS,1,Q0/AREA                       !采用换算密度，且为N/m3单位
! 2.在弦线位置创建模型，施加约束和自重荷载
K,1
K,2,L0
L,1,2                                       !创建几何模型
LESIZE,ALL,,,ENUM
LMESH,ALL                                !生成有限元模型
D,NODE(0,0,0),ALL
D,NODE(L0,0,0),ALL                       !施加约束
NODE1=NELEM(ENUM/2,1)                   !获得中间单元的两个节点号
NODE2=NELEM(ENUM/2,2)                   !以备后面使用
ACEL,,1.0
FINISH                                     !施加值为1的加速度(自重)
! 求解、进人后处理获得索内力、更新有限元模型、判别收敛条件是否满足
PASS1=1
*DOWHILE,PASS1                  
/SOLU
ANTYPE,0
NLGEOM,ON
SSTIF,ON                       !定义静态求解、打开大变形与应力刚度选顶
NSUBST,20
OUTRES,ALL,ALL
SOLVE
FINISH                                     !定义子步数、输出结果、求解
/POST1
SET,LAST,LAST                   !进入后处理，选择最后荷载步的最后子步
*GET,NFOR,ELEM,ENUM/2,SMISC,1      !获得跨中单元的索力并计算其余弦
COSREF=(NX(NODE2)-NX(NODE1))/DISTND(NODE1,NODE2)
NFOR=NFOR*ABS(COSREF)       !计算跨中单元的水平张力(可用其他单元)
ERR1=ABS(NFOR-H0)/H0
FINISH                                   !计算当前索水平张力误差
/PREP7                                   !进人前处理，更新有限元模型
*IF,ERR1,LT,0.05,THEN                     !如果误差小于5%时
UPGEOM,0.1,LAST,LAST,EX820,RST
*ELSE                                    !模型更新系数采用较小数值
UPGEOM,1,LAST,LAST,EX820,RST
*ENDIF                                !否则模型更新系数采用较大数值
*IF,ERR1,LT,ERR0,EXIT
*ENDDO                                   !如满足迭代条件，退出循环