（广西博阳电力勘察设计有限公司）
摘要：对覆冰情况下的高压输电线路进行了风险评估研究，建立了有限元计算模型。高压输电线路覆冰的种类、形成机理及影响覆冰的因素.。分析输电线路冰害事故的成因和线路覆冰的主要危害。本文首先建立了输电线路塔-线系统力学分析模型，然后就覆冰塔-线系统安全裕度进行了分析，最后以实际线路为例计算了其杆塔在不同风速和覆冰厚度情况下的单元应力值，对杆塔安全裕度进行了分析并提出结构改进措施。
关键词：高压输电线路；覆冰；塔-线系统；安全裕度；风险评估
0引言
覆冰是自然界常见的现象之一，而输电线路的导线、地线和绝缘子覆冰都将不同程度对输电线路的安全稳定运行造成不利影响。2008年我国部分省份出现的罕见冰冻灾害，对电网造成了大面积的损毁，诸多科研工作者和科研机构开始对覆冰情况下的输电线路设计标准进行了修订，电网运行单位也有针对性地开展了地区差异化的电网规划工作。
目前，对输电塔-线系统覆冰受力分析计算及风险评估的文献相对较少。文献建立了恶劣天气情况下的电网稳定性模型，该模型主要用于研究在暴风雪和大风等恶劣天气情况下，瑞典部分区域输电线路和铁塔的失效特性；提出了在风荷的塔线系统解析数值模型，应用该模型计算得到的数据用于分析塔线系统的动态特性和失效机理；建立了多跨塔-线系统的二维有限元数值计算模型，分析了脱冰和断线等情况下，输电线路的铁塔和导线动力响应方面的问题。国内方面，在总结中国输电线路典型覆冰事故的基础上，对中国输电线路覆冰事故特点、事故原因及防治技术措施进行了总结分析；通过对一例架空地线因覆冰滑动的事故，分析了在微地形、微气象区导地线覆冰后产生滑动的原因及危害，同时提出了解决方案；采用梁单元模拟覆冰输电线，研究了风速、档距长度对舞动幅值及导线舞动过程中张力的影响，为导线舞动的预防及减小舞动的幅值提供了依据。
1高压架空输电线路力学分析模型
按照有限元力学分析建模要求，高压架空输电线路数值计算模型的建立需要选取正确的杆塔、导、地线单元精确模拟各部件结构特性，设定贴近实际的材料特性，建立精细的塔-线耦合仿真模型，进而通过力学有限元计算来分析高压架空输电线路塔-线结构系统受力情况。
1.1塔-线结构系统计算模型
有限元模型建模所应用的力学单元如图1所示，其中：LINK10截面为四边形，用以模拟导地线的悬链线模型；铁塔使用Q345和Q235两种型号角钢，钢截面为“L”形，采用BEAM188模拟；覆冰时，利用梁单元BEAM188的ctube模型在导(地)线外包围一层厚度为15mm的冰。
根据构建塔-线系统的材料型号和特性选定相应单元，建立用于数值计算的塔-线耦合系统有限元计算模型。在塔建模过程中，为提高工作效率，考虑塔的平面对称特性，首先建立1/4塔模型，随即通过两次复制得到塔的整体模型。
1.2荷载的施加与计算
将导线离散成单元形式进行数值计算，对导线风荷载采用单位荷载进行施加。导地线水平档距为lH时，其风荷载为
1.3塔-线系统初始形态确定
为在重力求解后能够真实模拟线路实际运行状
态，同时要考虑塔线间的耦合作用，减小其它工况下(覆冰荷载、风荷载等)有限元计算的误差，必须对塔-线系统的初始形态进行确定。
塔-线结构系统找形采用如图2所示迭代修正法，该方法以重力求解后的导地线应力与其年平均运行应力之差、初始弧垂加上位移之和与设计弧垂之差为判据，通过迭代法来修正导地线的初始应力和初始应变，从而使塔-线系统在重力求解后的状态与线路正常运行状态一致。
迭代修正法采用线性静力求解，不仅避免了刚度矩阵的奇异性，还具有适用于曲线索单元和直线杆单元、求解控制简单、实现容易等特点。
2案例研究
2.1塔-线系统建模
针对我国西南地区某实际线路10mm冰区(易发生覆冰事故微气候段)的两个耐张段内41基杆塔，建立1:1三维有限元模型。按前文所述方法对逐基杆塔进行循环加载计算，对导、地线及杆塔进行加载，并对导、地线端部及杆塔脚部施加约束，得到杆塔在不同风速和覆冰情况下的单元应力值。图3所示为冰厚10mm、风速12m/s时一塔两线系统应力分布计算结果示意图。
2.2应力计算分析
以某杆塔为例分析其在不同风速冰厚临界情况下的应力分布情况，随着荷载的增大应力值逐步增大且在20-60、380-510、750-870及1150-1250四个塔头位置的单元应力均出现较大值，这四个位置对应于杆塔的塔头部分，说明钢构应力较大值不是单独出现，在一定范围内钢构的应力值均较大。
杆塔失效并非一两个钢材超过其屈服强度就失效，而是大范围的钢材均超过屈服强度时杆塔才会失效。通过临界冰厚风速情况下部分单元应力变化可知，35和409号单元在冰厚为28mm，风速为14m/s时，其应力值最大，此部分最为薄弱；1239-1241号单元在冰厚为12mm，风速为22m/s时应力值最大，此部分最为薄弱。
杆塔最大屈服强度随覆冰和风速变化情况，覆冰厚度在0~12mm之间杆塔的轴向应力基本保持不变即此时杆塔的稳定性较好；12~16mm之间随着覆冰厚度的增加轴向应力上升较快且上升的曲率较大，在覆冰厚度为16mm时杆塔的轴向应力达到最大。虽然风速在20m/s和22m/s时最大轴向应力超过了钢材的屈服强度，但绝大部分钢材还在屈服强度内，因此塔-线系统还没有失效；16~18mm时随着覆冰厚度的增加杆塔的轴向应力有所下降，即覆冰的影响对塔线系统的稳定有促进作用；随着覆冰厚度的增加，轴向应力继续增大。覆冰厚度相同时，随着风速的增大杆塔最大应力增大，且风速越大效果越明显。
2.3位移计算分析
杆塔合位移随覆冰和风速变化情况，在风速一定时，覆冰厚度为0~14mm时杆塔上的节点X向、Y向、Z向和合位移增加的较缓慢即节点位移却会改变较大，在覆冰厚度为16mm时杆塔上的节点X向、Y向、Z向和合位移达到最大，覆冰厚度为16~18mm时杆塔上的节点X向、Y向和合位移有减小的趋势，即此时覆冰对整个系统的稳定性起到促进的作用，最后再增加覆冰厚度杆塔上的节点X向、Y向、Z向和合位移均增加。此规律与应力的变化规律类似，且风速越大杆塔X向、Y向、Z向和合位移增大的趋势越明显。
2.4安全裕度计算
根据杆塔的临界情况下的风速与覆冰厚度，将获得的数据作为函数关系的因变量，通过回归分析，对仿真结果数据进行拟合，绘出失效评定曲线，得到精确的安全裕度表达式，为单塔裕度计算提供基础。
经过对多项式、指数形式等多种表达式反复计算验证，对杆塔失效评定曲线对比得到最精确的单塔裕度表达式为
杆塔安全裕度曲线，由于该塔在线路中的档距较大，因此曲线在起始位置的值不大，随着冰厚的增加，曲线呈下降趋势，由于力矩的作用冰厚逐渐达到20mm时呈上升的趋势，当冰厚继续增大时曲线迅速下降。
对于杆塔受力来说各变量是独立变量，以风速与覆冰厚度两个独立变量求解安全裕度，将拟合的失效函数作为反映各变量分布的失效临界函数F=f(v,T)，假定[VisTi]取某一组值为[10,10]，则可得杆塔安全裕度范围P。
2.5杆塔结构改进
通过有限元数值计算结果表明，杆塔薄弱处的位置大部分位于塔头部分，其原因在于相邻的杆塔高差档距较大，承受着两端线路水平档距和高差不同产生的不平衡张力，由此导致塔头部分受到较大的拉力，当应力较大时产生大的变形。
位于塔头之上薄弱钢材的屈服强度为Q235，通过提高钢材屈服强度，采用Q345钢材进行结构改进。标红的钢材为采用高强度钢材的部位，采用与之前相同的加载求解方式，得出此种情况下耐张段薄弱位置受力情况，将其进行改进后，该塔的临界条件有了明显的优化。
为杆塔应力比值最大的十个单元改进前后应力比值对比，证明提高局部钢材强度能够改变杆塔薄弱位置分布，使得较高强度的钢材承受更大应力，从而更充分发挥高强度材料的作用。
将薄弱处的钢材换成高强度的钢材后，在相同风速条件下，杆塔所能承受的极限冰厚值都有所提高，提高的范围在2~6mm不等，改善效果较为明显。
3结论
(1)高压架空输电线路塔-线系统有限元计算模型结合基于迭代方法的塔-线系统找形分析方法，实现了多档距线路在冰载荷及风载荷作用下杆塔及线路的受力分布计算。
(2)将安全预度概念引入杆塔风险分析中，通过安全预度函数以及安全预度曲线直观给出了杆塔在给定条件下的风险等级。
(3)以我国西南地区某500kV输电线路为研究对象，计算不同冰厚及风速情况下的杆塔应力值，通过回归分析方法，得出单塔失效安全裕度曲线数据。通过分析杆塔薄弱点位置受力情况，提出的杆塔改进措施切实可行，具有实际工程指导意义。
参考文献
[1]蒋兴良,常恒,胡琴,等.输电线路综合荷载等值覆冰厚度预测与试验研究[J].中国电机工程学报,2013,33(10):177-183.
[2]刘赟,俞集辉,程鹏.基于电磁-热耦合场的架空输电线路载流量分析与计算[J].电力系统保护与控制,2015,43(9):28-34.
[3]朱斌,潘玲玲,邹扬,等.考虑融冰因素的输电线路覆冰故障概率计算[J].电力系统保护与控制,2015,43(10):79-84.
[4]彭迎,阮江军.模拟电荷法计算特高压架空线路3维工频电场[J].高电压技术,2006,32(12):69-74.
[5]王少华.输电线路典型覆冰事故及防治技术分析[J].高压电器,2010,46(10):85-89.
[6]刘云,钱振东,夏开全,等.鼓型塔输电线路绝缘子破坏非线性动响应分析[J].振动工程学报,2009,22(1):6-13.