单索的找形和加载。
这是沈士钊老师《悬索结构设计》书中例题2-1.设有承受均布荷载的抛物线索，已知：A=0.674cm2,E=17000KN/ cm2,即EA=11460KN;又l=8m,H0=10KN, q0=0.2kn/m, q=0.5kn/m。求索内水平张力H以及索在始态和终态时的跨中垂度。（答案：H=18.98 kn, f0=0.160m, f=0.211m）。

我将索分为1m一段时，加载计算不成功（索内力不变或不收敛：提示位移过大，刚体运动）；索分段长度为0.5，0.2m时，连找形后，更新几何的计算都不收敛，提示small equation solver
pivot terms。只有采用了直接加荷载的办法找形（通过节点力加载或重力加载，结果相同。0.5m，0.2m分段结果相同。这样做正确吗？）

! 索分为1m一段,支座位移法找形
!找形时该而采用较小的弹性模量和较大的应变
/prep7
ET,1,LINK10
!Keyopt,1,2,2
!Keyopt,1,3,0
R,1,0.674e-4,0.9,
MP,EX,1,164853280
MP,NUXY,1,0.3
!定义节点
*do,i,1,9
n,i,1*i-5,0,0
*enddo  
!定义单元
Type1,1
mat,1
real,1
*do,i,1,8
e,i,i+1
*enddo
!施加位移约束
*do,i,1,9
d,i,uz,0
*enddo
*do,i,-4,4
d,i+5,uy, (0.01*i*i-0.160)
*enddo  

d,1,ux, 0
d,9,ux, 0
!设定时间步
!TIME,1
AUTOTS,0
NSUBST,20, , ,1
KBC,0
!输出每个子步的结果
OUTPR,BASIC,ALL,
OUTRES,ALL,ALL,
!设定大变形，应力刚化
NLGEOM,1
SSTIF,ON
FINISH
!求截
/SOLU
!/STATUS,SOLU
SOLVE
FINISH
后处理
/POST1
!PRNSOL,DOF,
PLNSOL,U,y,0,1
!PRESOL,SMISC,1
PLESOL,SMISC,1
!N=10001~10026N
finish
/PREP7
!更新几何形状
UPCOORD,1,0
!弹性模量恢复为真值
R,1,0.674e-4,8.7275e-4
MP,EX,1,1.7e11
!删除位移约束
*do,i,2,8
ddele,i,uy
*enddo
d,1,uy, 0
d,9,uy, 0
!求截

!TIME,2
/SOLU
!/STATUS,SOLU
SOLVE
FINISH
/POST1
!观察施加荷载之前结构的位移
!可见位移基本为零,预张力基本不变,为10000N
PLNSOL,U,y,0,1
PLESOL,SMISC,1
FINISH
!写荷载工况文件01
!/PREP7
!LSWRITE,01,
/prep7
d,1,ux, 0
d,9,ux, 0
d,1,uy, 0
d,9,uy, 0
*do,i,1,9
d,i,uz,0
*enddo
!二次施加节点荷载
*do,i,2,8
f,i,fy,-300
*enddo
!设置第2荷载步
Finish
/solu
!TIME,3
AUTOTS,on
NSUBST,15, , ,1
KBC,0
!neqit，50
l!nsrch,on
Cnvtol,f,,0.05
!pred,on
sstif,on
nlgeom,on
nropt,full
/SOLU
solve
!eqslv,spar
!LSWRITE,02,
FINISH

!LSSOLVE,1,2,1,
!后处理
Finish
/POST26
!求解 ,不成功？！

先顶一下，别沉了以后有时间在研究。