想要计算输电导线在各种气象条件下的线形和内力情况。
第一步根据在年平均气温条件（荷载为自重荷载，气温15度）下导线张力为定值，计算出导线的初始形态。这一步已经计算出来了，跨中最大弧垂36.006m，结果合理。
第二步要计算大风条件（荷载为Y方向自重和Z方向风载，气温-5度）下的线形和内力。但这一步需要考虑导线降温20度。我试着用tunif命令，但使用tunif,20和tunif,-20，计算结果居然是一样的。BFE命令也试过，但是结果也不合理，不知道该如何考虑降温的影响。
附上命令流，请大家帮助分析分析。

/clear
/filname,daoxian
/prep7
!导线参数
LH=1043           !档距
h=0               !高差
area=729.5e-6     !截面积m2
g1=3.275*9.80665  !单重N/m
em=9.7e10         !弹模Pa
h0=166.39*729.5   !最大使用张力
err0=0.001        !控制误差
enum=1000         !单元数量
istran=1.0E-7

!创建单元
et,1,link10
r,1,area,istran
mp,ex,1,em
mp,prxy,1,0.3
mp,dens,1,g1/area    !换算密度，N/m3
mp,ALPX,1,,16.2E-6   !热膨胀系数
!mp,reft,1,15        !参考温度

!创建模型，施加自重荷载
k,1
k,2,LH,-h
l,1,2
lesize,all,,,enum
lmesh,all
d,node(0,0,0),all
d,node(LH,-h,0),all
!d,ALl,UX
node1=nelem(enum/2,1)       !获得中间单元的两个节点号
node2=nelem(enum/2,2)
acel,,1     !施加值为1的加速度（自重），由于前面密度采用了N/m3的单位
finish
!求解、更新、判别收敛条件
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)         !cos方向角
nfor=nfor*abs(cosref)       !水平方向张力
err1=abs(nfor-h0)/h0        !与最大使用张力的误差
finish
/prep7
*if,nfor,lt,h0,exit           !张力过小时，跳出循环
*if,err1,gt,0.1,then            !误差大于0.1
upgeom,0.5,last,last,daoxian,rst  !将位移缩小为0.5倍加到有限元几何体上
*else
upgeom,0.03,last,last,daoxian,rst    !将位移缩小为0.03倍加到有限元几何体上
*endif
finish
*if,err1,lt,err0,exit       !误差小于err0（本例为0.001），则跳出循环
*enddo

/prep7
upgeom,1,last,last,daoxian,rst  

!求解大风荷载作用下的内力和变形
/prep7
!TOFFST,273     !定义绝对温度偏差
ALLSEL
!TREF,0        !定义参考温度
!TUNIF,-20      !降温20度
/solu
esel,all
!BFE,all,TEMP,,-20,-20,-20,-20
!acel,,1.0,48.321E-3/(g1/area*1.0E-6),        !Y方向自重，Z方向风载
nlgeom,on
nsubst,20
outres,all,all
solve
finish