覆冰分裂导线舞动的ANSYS有限元分析法.docVIP

覆冰分裂导线舞动的ANSYS有限元分析法 　　摘 要：为分析覆冰分裂导线舞动时的应力变化情况，研究导线舞动和指导连接金具设计，根据现有比较全面的三自由度导线舞动模型，利用ANSYS有限元软件建立了山西地区典型的紧凑型输电线路覆冰导线的舞动模型。利用该有限元模型对6分裂覆冰导线的舞动情况进行了动态仿真，输出了导线舞动时变化图像，可以直观的看到舞动时导线和间隔棒内部应力分布，为研究导线舞动及设计连接金具受力提供有力支持。 　　关键词：舞动 ANSYS 有限元 分裂导线 间隔棒 　　中图分类号：TM72 文献标识码：A 文章编号：1672-3791（2014）11（b）-0063-04 　　对于输电线路导线舞动的研究，国外从20世纪30年代就已开始，具有代表性的研究成果是美国Den Hartog的垂直舞动机理和加拿大O.Nigel的扭动舞动机理[1]，此外还有惯性耦合机理[2]和稳定性机理[3]等。近几十年来，输电线路舞动的计算机数值模拟分析方面取得了巨大成就，P.Yu和 Desai等人提出了导线舞动的三自由度模型[4～5]，Q.Zhang等人建立了混合模型用以研究分裂导线的舞动[6]，国内王丽新，杨文兵等人利用梁单元，模拟覆冰单导线，探讨了风速对舞动幅值的影响[7]。 　　尽管前人在这一领域已做了大量工作，但仍有一些问题没有解决，比如覆冰分裂导线中间隔棒的精细3维模型的建立、舞动中应力变化的分析等。由于随着紧凑型输电线路在各种气候地带架设，研究分裂导线间隔棒随导线舞动时的应力变化情况对指导设计间隔棒具有重要意义。 　　该文以P.Yu[4]提出的数学模型为基础，利用有限元软件ANSYS计算了山西常见紧凑型输电线路的不同风速下的舞动情况和间隔棒受力变化情况，这种应力分析模型为设计合理的间隔棒连接金具[8]应对断裂提供了理论基础。 　　1 覆冰分裂导线的有限元建模 　　考虑到导线覆冰时水平、垂直与扭转的三自由度共同作用情况，导线未覆冰时的受力情况如图1（a）所示，以及覆冰时的情况如图1（b）所示，图中的气动力包括阻力FD，升力FL和扭矩FM，θ为覆冰模型与风构成的角度。 　　这里使用ANSYS软件选用beam188三维梁单元建立输电线的有限元模型，导线两端用梁单元beam4单元来模拟相邻档对导线产生的影响，分裂导线的间隔棒采用采用beam4单元，这是一种可用于承受拉、压、弯、扭的单轴受力单元，可用于计算应力刚化及大变形的问题，实现3D弹性受力分析。假定导线在空中的几何形态视为悬链线形态，用与悬链线方程相接近的抛物线方程模拟架空导线的形态，对导线施加初应变、重力加速度和冰载荷，通过静力求解进行迭代，直到导线重力作用下的几何形态与理论值误差满足工程要求为止。 　　初始应力采用命令INISTATE直接在beam188单元上施加覆冰时导线拉力值。 　　输电线路的非线性动力学分析求解的基本方程是： 　　（1） 　　其中，M，K和C分别为单元质量矩阵、刚度矩阵和阻尼矩阵，F为单元载荷向量，在本程序中包括静载和动载，动载主要是空气动力载荷。根据流体诱发振动理论，长为L的覆冰导线在速度为U的水平风作用下，输电导线受到的气动力包括阻力FD，升力FL和扭矩FM动力，可分别表示为： 　　（2） 　　（3） 　　（4） 　　其中ρ为气流密度，U为吹向导线的水平风速，D为输电导线的直径；CL，CD，CM分别为升力、阻力和扭矩系数，与导线截面、覆冰形状和厚度、运动状态以及攻角有关。要得到3个气动系数，必须首先确定攻角，其计算公式是： 　　（5） 　　其中是初始角度，是由导线垂直振动引起的攻角变化是由导线扭振运动引起的攻角变化。 　　实际作用在导线上的水平力及垂直力为： 　　， 　　（6） 　　2 覆冰分裂导线舞动的ANSYS有限元计算 　　基于上述的三自由度模型，选取山西地区的一条紧凑型输电线路物理数据来设计导线舞动的ANSYS计算模型。 　　（1）山西地区紧凑型500kV导线规格及物理性能数据：这里选取LGJ-300/40钢芯铝绞线（六分裂，等边形布置，分裂间距选用500mm）；导线计算截面338.99 mm2，计算拉断力92.22 kN，计算材料材料选择为DENS=2518.3、EX=7.3E10、PRXY=0.2；导线跨度为450 m。 　　（2）气动参数选取列表如下。 　　表2选取参考文献《架空输电线舞动的计算机仿真》[9]中的计算内容。 　　其中空气密度通常取如下情况，即20°时，取1.205 kg/m3， 　　（3）导线初始情况设定：导线初始张力（预应力）为年均21935 N/根；导线弧垂按照抛物线架设情况选取16 m；导线年平均覆冰情况为10mm。 　　3 ANSYS有限元计算结果