基于有限元的风偏角计算的动力时程分析.docxVIP

基于有限元的风偏角计算的动力时程分析 0 线路风偏系统研究 通过分析近年来经常出现的电网运营故障的原因，我们可以看到，风致是一个相对频繁的类别。据统计, 2007—2010, 国家电网辖区高压、超高压输电线路因风偏分别发生跳闸157、93、79、174次, 涉及北京、山东、山西、河南、江苏、安徽、湖北、浙江等省市。风偏跳闸后重合闸成功率低, 影响供电可靠性;此外, 随着特高压输电技术被列入我国的“十二五”规划, 特高压线路在我国的国民经济中的地位越来越重要, 保障其安全运行具有重大的战略意义, 因此对风偏进行系统研究是十分有必要的。 关于线路风偏的计算, 我国现行电力行业标准中采用的是刚性直棒法, 该方法是一个简化计算公式, 其对真实的线路环境做了过多的简化处理, 得到的结果误差较大, 因此很多学者都借助有限元的方法对规范方法开展了大量优化工作。肖正直、娄文娟等研究得到了分裂导线的体型系数;卢文生、李宏男等对大跨越输电塔线体系的脉动风场进行了模拟, 得到了可用于塔线体系动力计算的脉动风时程;邵瑰玮等对输电导线悬挂点处的等效风荷载进行了仿真分析;李黎等运用ANSYS建立了输电塔-绝缘子串-导线的耦合模型, 并且施加脉动风速时程计算了该模型的风偏角;严波等运用ABAQUS分别建立了2档绝缘子串-导线耦合模型和绝缘子串多刚体模型, 并进行了脉动风时程分析, 在此基础上提出了考虑动力效应的风荷载调整系数。 以上研究都是关于超高压及以下电压等级的线路, 对于特高压线路而言, 还存在以下问题:①没有考虑导线与绝缘子串之间的耦合作用;②特高压线路所用导线分裂数多, 设计中只是按单导线的风荷载简单叠加, 没有考虑分裂导线之间屏蔽尾流的影响;③没有得到脉动风对风偏角的动力放大系数。针对这些问题, 本文基于有限元仿真, 考虑绝缘子串-导线耦合作用、屏蔽尾流影响以及脉动风动力放大效应, 对风偏角计算的刚性直棒法提出系统修正建议。 1 风荷载和重力荷载 刚性直棒法由静力平衡和力的三角形原理推导而来, 基本模型如图1所示。图1中, GH、GV分别为直棒形心处绝缘子串的风荷载和重力荷载, 按文献中规定计算;WH、WV分别为绝缘子串末端传来的导线风荷载和重力荷载;φ为风偏角。刚性直棒法按下式对风偏角进行估算: 式中:pH、pV分别为导线单位长度上的水平风荷载和垂向荷载, 根据荷载规范的规定计算;l1、l2分别为杆塔两侧的左、右档距;h1、h2分别为杆塔两侧导线挂点相对中间挂点的左、右高差;T0为导线张力。 2 刚性直棒法的修正建议 风偏计算的刚性直棒法中, 绝缘子串被当作独立的系统, 忽略了导线对绝缘子串的耦合作用, 由此而带来了一定的误差, 以下就绝缘子串和导线的耦合作用进行研究, 提出对刚性直棒法的修正建议。 2.1 参数参数 本文的研究中, 选取陶瓷绝缘子为54×XWP-300, 导线为8×JL/G1A-630/45, 具体参数如表1和表2所示。 对于输电线路的表征参数:档距和高差, 本文根据实际的特高压线路, 选取等档距800m和不同高差的组合, 确定了5种基本工况, 如表3所示。 2.2 大径管网初始张力的计算 采用ANSYS软件来建立串-线有限元模型。根据陶瓷绝缘子串片与片之间的连接具有一定转动能力, 以及文献的研究:相邻绝缘子片间的摩擦阻尼可以忽略, 以BEAM188单元来模拟绝缘子片, 每片绝缘子简化为1根圆杆, 圆杆之间为铰接连接, 圆杆直径取为圆杆长度的1/10, 由此计算出等效密度为ρeq=2.23×105kg/m3, 绝缘子串上端铰接与杆塔 (固定端) 。 对于导线, 其局部受力具有梁的特性, 故采用BEAM188单元对其进行模拟。单元直径取为导线直径, 长度取10m一个单元, ρeq=2.317×103kg/m3, 初始张力根据拉断力进行估算, 取为44MPa, 将初始张力代入悬链线方程求得导线的初始形态。 根据以上要点分别建立2档、4档和6档绝缘子串-导线耦合模型, 以工况3为例, 如图2、3所示。 2.3 刚性直棒法与4档耦合有限元模型的对比 图4为绝缘子串的风偏形态。取风速30m/s, 分别对2、4、6档串-线耦合模型施加平均风荷载进行静力分析, 同时按照刚性直棒法对5种工况的风偏进行计算, 结果对比如图5所示。图5中, M为工况编号。 由图5可见: 1) 4档和6档耦合有限元模型的风偏角相同, 大于2档模型, 相差约2.5°, 说明4档模型对串-线的耦合作用已进行充分考虑, 比2档模型更精确, 比6档模型计算量小, 因此本文以下的研究中均只考虑4档耦合有限元模型。 2) 当高差为0时, 刚性直棒法与4档模型的风偏角基本相同;当高差不为0时, 随着高差绝对值的增大, 两者的误差逐渐增大。而高差又与导线垂直荷载相关, 高差由