基于Ansys的高压输电塔风振特性分析_耿贞伟.pdf

基于 的高压输电塔风振特性分析耿贞伟，苏文伟，于凤荣，钱晶（云南电网有限责任公司信息中心，昆明；昆明理工大学 冶金与能源工程学院，昆明）摘要：在电力大范围应用的背景下，研究输电线路的振动特性对电网整体的安全性有着极其重要的意义。文中推导了三线两塔耦合系统的力学简化模型，利用 软件建立了典型的 交流双回路耐张塔模型，在此基础上利用 o 风速谱对电力塔模型节点的风载荷进行了模拟，阐明了分析的有效性。利用理论方法和有限元方法分别计算了塔线耦合体系的动力特性，结果表明理论计算与有限元结果保持一致，证明了理论的有效性。此外通过有限元的分析还能看出耦合系统出现了一部分由塔或线诱发的新振型，且同阶下耦合的振动频率要小于单塔的振动频率。关键词：输电线路；耐张塔；力学分析；有限元分析；风速响应谱：中图分类号：；文献标识码：文章编号：（），o，（Y ，g 5，c g g gg，g cc cg，g 5，）：o o o o ，o o o o o o o o o ，o o oo o ，ooo o o o ，o o o o o o o o o，o o o o o o o o o o o ，o o o o，o o o o o o o，o o o o o o o：o，o o，o o基金项目：国家自然科学基金资助项目（）引 言随着西电东输等大型电力工程的稳步推进，电力能源在全国范围内均得到了更大范围的应用。特高压输电塔是电力传输系统中的最重要的部件之一，其维修维护难度大，成本高，因此分析其可靠性和安全性至关重要。过去常用的低电压线路已难以满足日益增长的经济需求，相比之下特高压输电线路能够创造更大的经济效益。对这类架空输电线路，尤其是特高压输电线路来说，由于输电量和工作电压均有大幅提升，电气间隙和多分裂的导线截面积都会增加。这一特点不仅为塔身及塔臂带来了更多负荷，同时也导致在塔线耦合这个体系中，塔与线二者振动的相互影响已无法忽略。统计表明，输电塔坍塌事故主要是由风振导致的，李宏男等研究了电力塔在多种激励下的动力学特性，并综述了电力塔的安全评估体系。文献提出了电力塔的多质点简化模型，在多质点模型中，将导线与塔身等效为利用刚性杆件连接的多个质点模型，并利第 卷 第 期电测与仪表 年 月 日 ，用该理论模型计算了典型塔的振动特性。文献利用多波屈曲单元对电塔进行理论分析，有效预测了输电塔结构的抗风极限载荷。文献对 特高压长悬臂输电塔线耦合体系的风震特性进行了分析，并得出了悬臂长度对稳定性影响的定性分析结果。目前的研究集中在风振响应或塔本身的振动特性上，因此本文中对风振响应和耐张塔的振动特性进行联合分析具有一定的工程价值与意义。本文综合考虑了现有的电力塔建模方法，利用 有限元软件建立梁与杆的混合模型。由于要建立塔线耦合的模型，进一步使用了悬链法对导线与地线进行几何找形，同时合理地处理塔线耦合的边界条件。最终建立了两塔三线的耦合有限元模型，并对两塔三线的简化模型进行了振动特性的理论分析。在得到振动特性的基础上，又利用有限元动力时程分析方法，进一步分析了塔线耦合体系的风振响应。两塔三线耦合体系振动特性理论对一个多塔、线耦合的体系，通常选用一塔两线的分析方法，这样既可以简化模型，又能得到相对精确的结果。但相比于一塔两线模型，两塔三线模型不仅能反映塔线耦合特性与线线耦合特性，还能反映塔塔耦合特性，因此本文选择两塔三线模型进行研究。两塔三线模型如图 所示，其中是塔的抗弯刚度，是导线的抗弯刚度。由于电塔的跨度较大，边界条件变得模糊，因此导线的边界约束对整体的影响可以忽略，简化模型中没有考虑绝缘子等连接部件。zHE1I1E1I1E1I1lA1lA1lA1x1x2x3x4图 两塔三线耦合简化模型 oo o o对于图 所示的系统，采用模态综合分析法分析塔振动与线振动相耦合的动力特性，塔的固定截面主模态与约束模态可分别表示为：（）（）（）（）（）（）式中 表示塔身总高；代表沿塔身高度。三根长为 的导线的固定截面主模态和约束模态可表示为：（）()（）（）（）（）()()()()()()()()()()()()（）()（）()（）（）（）式中、分别为图 中标注的沿导线长度。以式（）为基础可以得到塔的横向位移：（，）（）（）（）（），（，）（），（）（）（）（）同理也可以得到三根导线的面外位移：（，）（）（）（）（）（，）（）（）（）（）（，）（）（）（）（）（，）（）（）（）（）（）式中（）是广义坐标。子部件之间相互的位移关系满足：（，）（，）（，）（，）（，），()（）因此可得力协调条件：（，），()，()（，），()，()（）可以进一步得到：（）（）（）（）（）（）（）（）（）（）（）（）（）第 卷 第 期电测与仪表 年 月 日 ，式中 ()；代表阻尼系数。本方法中我们将导线等效为受张力的梁，抗弯刚度可等效为：k()()()（）式中 是载荷的角标；是导线长度；是导线直径；是导线横截面积。通过式（），可以将阻尼系数化简为：（）综合式（），塔线耦合运动系统的张量表达式为：（）通过式（）中的坐标关系，可知 中实际只含两个变量，因此又有 ，那么 与 间的关系可以描述为：（）其中的转置矩阵 为：（）因此式（）转换为：（）和 分别为质量阵和刚度阵：（）其中元素分别为：o()()cgc ()o()o()cc（）o()()（）式中 为模态的阶数。塔线耦合系统的振动圆频率 k 可以利用下式计算：k（）得到 的特征向量后，将特征向量代入式（）可以得到 的模态向量，综合式（）式（），可以得到总耦合体系的振动模态。单塔的风速响应谱分析为了进一步验证所建立的两塔三线模型的合理性，首先对单塔模型进行了风速响应谱分析。单塔有限元模型用于分析模态及塔身受力状况的线框模型如图 所示，许多复杂结构简化为了桁架形式。塔高，呼高 ，结构尺寸由下至上线性变化，根开为 。结构主要使用的材料是角钢，材料型号选用，两种。为了合理计算电力塔结构，主材等重要部件使用 梁单元，其余辅材使用 杆单元，在 中赋予角钢截面后的有限元模型如图 所示。图 单塔线框模型 o o 单塔风速响应谱分析根据现有的大量对风的记录，脉动风可作为一个高斯平稳过程来考察。本文中使用常用的 o 风速谱对脉动风速进行模拟：（）k ()（）式中（）表示脉动风速功率谱；表示脉动风频率；k 表示地面粗糙度系数；表示标准高度 处的平均风速。第 卷 第 期电测与仪表 年 月 日 ，图 单塔有限元模型 o o o由于 o 风速谱与高度无关，所以引入一个随机函数来反映不同高度下的风速脉动，此函数的均值为，任意两点间的互相关函数定义为：，()()（）式中，分别为两点的高度。这里引入高度是为了方便在有限元软件中对不同位点施加载荷。模拟总时长为 ，时间增量 ，风速 ，地面粗糙度系数 ，通过上述分析可得到各节点处的脉动风速，图 以塔顶节点处风速即高度 风速为例展示了计算结果。同时也可得到风速功率谱密度，以塔顶节点处的风速功率谱密度为例，如图 所示，模拟出的 o 风速谱与目标模拟谱基本一致，自功率谱的能量分布与实际情况在趋势上一致，数值上接近，说明模拟出的 o风速谱是真实可用的。3020100-10-20-30风速/（m/s）050010001500时间/s图 风速时程曲线 o o 为了避免使用 带来的复杂计算流程，本文将塔顶的风速转换为风对塔身不同位点的等效载荷，这里使用标准风压公式：g（）式中 表示空气的重度；空气g；表示风速。结合脉动风压、塔的形体系数、风投影面宽，可以进一步得到各节点处的风载荷。10010-5功率谱密度频率/Hz10-2100图 风速功率谱密度 o 以图 风速时程曲线为例，进一步将风速转换为塔顶的风载荷，结果如图 所示。3020100050010001500时间/s位移/mm图 风载荷时程曲线 o o 单塔的风致响应分析在 软件中定义水平风载荷，将上面得到的载荷施加到各节点上，以分析塔身在耦合体系中的节点位移变化规律。载荷时程曲线如图 所示，以塔顶位节点移情况为例，可得到各节点处位移，如图 所示。200150100500050010001500时间/s位移/mm图 塔顶节点位移曲线 o o 第 卷 第 期电测与仪表 年 月 日 ，根据图、图 的结果，数值符合预期，证明单塔模型的建立是合理有效的，可以进行下一步耦合模型的分析与建立。塔线耦合有限元模型由于导线属于柔性体，在重力作用下形成自然下垂的形状，建立耦合模型就需要寻找该状态下导线的精确几何外形。针对这类问题，我们常用悬链线法进行计算，图 给出了典型的悬链线问题。，分别为固定端，是下垂的最低点，以 为原点建立标准，坐标系，定义，坐标分别为（，），（，），两点间的水平距离为。yBAOx图 悬链线图示 悬链法的找形方程为：c ()（）式中 为悬链线系数。该系数可以利用下式计算：()（）可以积分得出悬链线长度：()（）塔线耦合模型中还需要考虑线与塔之间如何连接，本文中将绝缘子视为刚体，绝缘子与导线间使用关节（o）连接，即该节点拥有三个转动自由度，如图 所示。导线的单元定义为梁单元，固定边界定义为塔底四角与架空导线的两端。为了使得仿真结果与实际情况更加贴合，建立了三线两塔双回路模型，塔距取典型值 ，导线与地线型号分别为 ，最终的有限元模型如图 所示。图 绝缘子模型 o o图 三线两塔有限元模型 oo o 塔的风致响应和振动特性分析 单塔振动特性分析单塔前十阶振型的有限元计算结果如表 所示，其中一阶振型和二阶振型分别是垂直于导线和平行于导线的两个方向，这两个振型的频率在数值上十分接近。表 单塔振动特性 o o o振动阶次振动频率 振动方式 横向弯曲 纵向弯曲 扭转 横向弯曲 局部模态 横向弯曲 纵向弯曲 扭转 局部模态 扭转 塔线耦合体系振动特性分析由于导线属于细长的柔性体，固有频率很低，且各个模态的振型与频率差别很小，因此频率的分布相当集中。在本算例中，我们共计算了前 阶模态，并选取了部分塔线耦合的典型模态进行了说明。根据有限元的计算结果，塔线耦合体系中绝大多数振型都是导线的振型，且相比于单塔的情况，耦合系第 卷 第 期电测与仪表 年 月 日 ，统的基频要更低，这也是因为导线的影响。例如前 阶模态频率主要分布在 ，附近，这是由于建模时考虑到实际条件，因此每根导线并不完全一致，尺寸在悬链线找形方程的附近有一定的浮动，因此每一根导线的振动频率都有细微差别。这四个频率分别是导线的水平一阶、二阶模态，以及垂直方向上的一阶、二阶模态。表 括号中的反代表两座塔振动方向相反，结合单塔的振动特性，可以直观地看出 与 阶是在线的影响下，两塔具有的共同耦合模态，频率略低于单塔的弯曲固有频率（约 ）。阶是个较为特殊的振型，首次出现了竖直振动的方式，如图 所示。这在单塔的模态分析中是没有出现的，因此这个模态是因为导线在竖直方向上振动而诱发的特殊模态。表 塔线耦合振动特性 o o o o振动阶次振动频率 振动方式 横向弯曲 纵向弯曲 纵向弯曲 纵向弯曲 纵向弯曲（反）扭转 扭转 扭转 扭转（反）纵向弯曲纵向弯曲（反）扭转（反）竖向（反）竖向横向弯曲 局部 扭转（反）横向弯曲（反）局部 扭转为了验证理论的正确性，依据第二节中的近似理论计算了一阶塔与线的，阶耦合理论值，如表 所示。表 塔线耦合体系振动频率理论值 o o o o o o 塔振动阶次与线 阶耦合与线 阶耦合与线 阶耦合与线 阶耦合与线 阶耦合 表 展示了近似计算理论反映的塔线耦合体系的趋势。低阶耦合状态下，系统的频率高于单塔频率，随着耦合模态阶数升高，频率逐渐趋于一个最小值，该最小值与有限元得到的塔线耦合体系的频率很接近（），对应耦合状态下的，阶振型，证明了理论切实反映了塔线耦合体系具有的物理特性。图 阶模态振型 o o 结束语推导了三线两塔耦合系统的简化动力学方程，并利用 软件建立了典型 高压耐张塔的线框模型，利用 有限元软件对耦合系统的振动特性和风致响应特性做了分析，最终得出以下结论：（）由于导线在本模型中是跨度极大（相对于导线尺寸）的柔性体，因此导线的振动特性对边界条件不敏