索单元找形后是处于悬链线状态吗？
最近使用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 !如满足迭代条件，退出循环
页:
[1]