摘要
根据输电线路张力放线施工中牵引走板与放线滑车的实际接触运动关系,构建了由钢丝绳、牵引走板、导线及双滑车等部件组成的多体运动系统,提出了符合实际情况的钢丝绳、导线、牵引走板等与双滑车的接触判断方法.基于有限质点法理论,结合部件间的空间位置关系,给出了单元节点的作用力表达式,并通过罚函数法给出了不同接触状态下牵引走板对滑车的接触力,形成了可分析牵引走板与双滑车结构接触碰撞的动力计算方法.根据张力放线施工实际情况,计算分析了牵引走板与双滑车结构接触过程中的滑车载荷、摆动角度的变化,并针对悬挂长度、走板厚度等施工因素,总结了牵引力、最大接触力、滑车载荷、摆动幅度等参数对系统安全的影响.本文提出的计算方法可指导工程中牵引走板与滑车的结构设计及施工安装,为输电线路放线施工提供技术参考.
牵引走板是输电线路工程张力放线施工中用于连接牵引钢丝绳和各子导线的专用工具,有一牵二、一牵四、一牵六等不同型式.走板由钢板焊接而成,由牵引钢丝绳带动,牵引导线通过输电线路铁塔上的各个放线滑车.当线路铁塔间的转角、高差较大时,需悬挂双放线滑车,增加滑车承载能力,减小导线对滑车的包络
随着特高压工程建设的开展,大截面导线普遍应用,牵引走板质量、体积也不断增加,在冲击力的作用下对滑车和导线会产生一定损伤,严重时会发生滑车破损而导致放线事故,造成较大的经济损失.因此牵引走板、放线滑车的设计与安装成为特高压工程放线施工中关注的重
有限质点法以向量力学和数值计算为基础,将结构离散为质点群,采用牛顿第二定律描述这些质点的运动,避免求解非线性方程组和整体刚度矩阵,适于计算发生刚体位移和几何大变形的结构,在处理结构大变形、接触、碰撞等非线性问题方面取得良好的应用效
本文在有限质点法的基础上提出牵引走板与双放线滑车系统的动力计算方法,分析牵引走板与双放线滑车在接触下的作用影响,为张力放线施工安全提供技术参考.
在牵引机的牵引力及张力机的张力作用下,钢丝绳、牵引走板及导线连续通过线路中各个铁塔横担下悬挂的放线滑车(组),直至放线过程结束.如

图1 牵引走板通过滑车
Fig.1 Traction plate passing through the pulley
将牵引钢丝绳、牵引走板、导线及双滑车、滑车悬臂、滑车连接杆看作由多个部件组成的多体运动系统.放线滑车悬挂在铁塔横担上,悬挂点J1、J2,并以J1为原点建立平面坐标系,如

图2 牵引走板与双滑车运动系统
Fig.2 An illustrative system of traction plate and double-pulley
在牵引力与张力的作用下,牵引绳、导线可视为柔性的线弹性单元,仅承受轴向拉力,运动过程中弹性模量及截面积保持不变.因此设置牵引绳单元AB,单元初始长度为sAB,弹性模量与截面积乘积为EAR.导线单元CD,单元初始长度为sCD,弹性模量与截面积乘积为EAC.牵引走板为钢制,在接触过程中无变形,视为矩形刚性单元,单元长度为sBC,厚度为bP,走板与钢丝绳、导线连接节点B、C位于单元轴心线,下部顶点为P、Q,如

图3 走板与滑车示意图
Fig.3 Schematic diagram of traction plate and pulley
悬挂单元O1J1、O2J2的初始长度为sOJ(在此认为2个悬挂单元长度相同),弹性模量与截面积乘积为EAJ.双滑车为完全相同的圆形刚性单元,中心为O1、O2(下文中以O1、O2表示2个滑车),半径为R.双滑车间隔单元O1O2,单元初始长度为sOO,弹性模量与截面积乘积为EAO.将单元质量均匀集中于质点B、C、D、O1、O2上,则各质点质量为:
mB=ρR·sAB/2+mP/2
mC=ρC·sCD/2+mP/2 | (1) |
mD=ρC·sCD/2
mO=ρJ·sOJ/2+mK
式中,ρR、ρC、ρJ分别为牵引绳、导线、悬挂单元的单位长度质量;mP、mK分别为牵引走板、滑车质量.
在实际中,牵引绳、导线只会与滑车的上部(高于中心O1、O2的部分)发生接触,因此为简化分析,从牵引绳、导线与滑车上部的接触点入手判断接触情况.节点A到滑车O1、O2上部的切点为A1、A2,节点B到滑车O1、O2上部的切点为B1、B2,滑车O1、O2的上部外切线切点为U1、U2,如

(a) AB与双滑车同时接触
(a) AB is contacted with two pulleys simultaneously

(b) AB与滑车O1接触
(b) AB is contacted with pulley O1

(c) AB与滑车O2接触
(c) AB is contacted with pulley O2
图4 AB与双滑车接触状态
Fig.4 AB is contacted with two pulleys
(1)当AB与滑车O1、O2同时接触,如
(2) |
其中,β1、β2为AB对滑车O1、O2的包络角,,;lO1A、lU1U2、lO2B、lA1U1、lU2B2分别为O1到A、U1到U2、O2到B、A1到U1、U2到B2的距离.
(3)当AB仅与滑车O2接触,如
(4) |
其中,;lAO2、lA2B2分别为A1到O2、A2到B2的距离.
(4)其他情况下AB与滑车O1、O2不接触.单元AB的长度SAB为A、B的距离lAB.
导线CD与滑车的接触判断可按上述方法实现,仅需将A、B替换为C、D即可.
牵引走板BC在运动体系中为刚性体,与滑车接触过程中不发生弯曲变形,牵引走板下部边界PQ与滑车O1、O2接触.因此可以通过直线PQ与滑车的关系判断接触情况.如

图5 牵引走板与滑车接触
Fig.5 Traction plate is contacted with the pulley
(1)当lO1P≤R,且lO1Q≥,滑车O1与P点发生点接触.
(2)当R≤lO1P<,且R≤lO1Q<,滑车O1与PQ发生区域接触,接触区域的中点为M.
(3)当lO1P≥,且lO1Q≤R,滑车O1与Q点发生点接触.
其中,lO1P、lO1Q分别为O1到P、O1到Q的距离.其他情形下走板与滑车O1无接触.走板与滑车O2的接触判断与O1相同,仅需将O1替换为O2即可,不再赘述.
在此采用有限质点法分析牵引走板与放线滑车多体系统的运动.系统中各单元的变形通过具有集中质量的质点位移表示,质点在单元内力和节点外力的共同作用下遵循牛顿第二定律进行运动.质点运动方程可表示为
(5) |
式中,m为质点质量;为质点加速度矢量;F为质点所受合力矢量.采用中心差分法来求解质点的运动方程.令h为时间增量步长,节点运动时间tn = nh,n ≥ 0;
(6) |
式中,,.
牵引绳AB与滑车发生柔性接触,线弹性单元AB内力可表示为
(7) |
单元AB对质点A、B、O1、O2的作用力fAB-A、fAB-B、fAB-O1、fAB-O2表达式可根据牵引绳与滑车的不同接触状态(
为简化叙述,以e表示单位向量,下标表示向量方向,即eAA1表示从A到A1的单位向量.其他向量采用同样的表示方法,不再一一说明.同样可以获得单元CD对C、D、O1、O2的作用力fCD-C、fCD-D、fCD-O1、fCD-O2.
走板与滑车O1、O2先后发生碰撞接触,其法向接触力等于侵入量与法向刚度系数的乘
fBC-O1= k(R-lO1P) ePO1 | (8) |
其中k为法向刚度系数(可取值1

图6 牵引走板与滑车接触
Fig.6 Contact force between traction plate and pulley
fBC-O1= k(R-lO1M) eMO1 | (9) |
其中lO1M为O1到直线PQ的距离,可根据三角形O1PQ的面积关系得到:
lO1M | (10) |
式中.同样,当Q点与滑车O1发生点接触时,侵入量为R-lO1Q,则单元BC对O1的作用力为
fBC-O1= k(R-lO1Q) eQO1 | (11) |
当不发生接触时,fBC-O1=0.走板为刚性体,与滑车的接触也对走板产生反作用力.由于走板厚度较小,认为滑车对走板的作用力直接作用于B、C.根据最大侵入点与端点B、C的相对位置,将接触力按比例分配到走板端点上,则滑车O1对质点B、C的作用力可表示
fO1B =-(1-λ) fBC-O1 | (12) |
fO1C =-λ fBC-O1 | (13) |
式中λ=.同样可得到牵引走板与滑车O2的接触力fBC-O2及O2对B、C的作用力fO2B、fO2C.走板为刚性体,在任意时刻应满足长度约束条件:
(14) |
已知tn时刻的位移,将tn+1的B、C位移表达
fBC-B = f2 eBC | (15) |
fBC-C = f2 eCB | (16) |
悬挂单元O1J1、O2J2和间隔单元O1O2均为线弹性单元,单元内力分别可根据单元变形得到:
f3=
f4= | (17) |
f5=
因此,单元O1J1对质点O1的作用力、O2J2对O2的作用力、O1O2对质点O1、O2的作用力分别为:
fO1J1-O1 = f3 eO1J1 | (18) |
fO2J2-O2 = f4 eO2J2 | (19) |
fO1O2-O1 = f5 eO2O1 = -fO1O2-O2 | (20) |
已知当前时刻
根据特高压工程中张力放线施工的典型技术参数,开展牵引走板通过双放线滑车的运动过程计算,分析施工中重点因素对滑车载荷、摆动角等参数的影响.
典型工况1:走板通过滑车时运行速度1.0m/s、牵引角度20°(牵引方向与水平的夹角)、张力角度-10°(张力方向与水平的夹角)、张力25kN、走板厚度0.1m、悬挂长度(sOJ)1.0m.计算中,参数ζ=0.2,h=1


图7 牵引走板通过双滑车全过程
Fig.7 The whole process of traction plate passing through the double-pulley
施工中主要关注滑车所受最大载荷,因此针对滑车所受载荷与滑车摆动角度进行分析,计算结果如

(a) 滑车1与滑车2所受载荷随时间变化曲线
(a) Time history of loads on pulley 1 and pulley 2

(b) 摆动角随时间变化曲线
(b) Time history of swing angle
图8 走板通过双滑车过程分析(典型工况1)
Fig.8 Process analysis of traction plate passing through double-pulley(typical condition 1)
通过计算发现牵引走板厚度和滑车悬挂长度对滑车载荷有较大影响,因此设置典型工况2:走板厚度0.3m、悬挂长度4m,其他参数与典型工况1相同.滑车载荷、滑车摆动角度结果如

(a) 滑车1与滑车2所受载荷随时间变化曲线
(a) Time history of loads on pulley 1 and pulley

(b) 摆动角随时间变化曲线
(b) Time history of swing angle
图9 走板通过双滑车过程分析(典型工况2)
Fig.9 Process analysis of traction plate passing through double-pulley(typical condition 2)
滑车的摆动角度与典型工况1近似,在中期因走板与滑车1、2的反复脱离、接触而有波动.从典型工况1、2结果中可知,随着放线滑车悬挂长度增大和走板厚度增加,走板通过双滑车的运行状态发生了复杂变化,产生了与滑车1和滑车2反复接触的变化过程,滑车载荷、摆动角度也有大幅增大.
特高压输电线路工程中,受实际条件影响,施工安装中滑车的悬挂长度有较大变化;另外牵引走板通过滑车时因空间位置不同而呈现不同厚度变化.因此分别计算悬挂长度3.0m-4.5m,走板厚度在0.1m、0.2m情况下,走板牵引力、走板与滑车的最大接触力、悬臂载荷、摆动幅度(向牵引侧摆动为正,向张力侧摆动为负)等参数的变化.计算结果如

图10 牵引力、接触力随悬挂长度的变化曲线
Fig.10 Curve of traction force and contact force with varying suspension length

图11 滑车1、2载荷随悬挂长度的变化曲线
Fig.11 Curve of loads of pulley 1 and 2 with varying suspension length

图12 滑车摆动角最大幅值随悬挂长度的变化曲线
Fig.12 Curve of maximum amplitude of pulley swing angle with varying suspension length
本文建立了输电线路张力放线施工中牵引走板与双滑车结构的多体运动系统,提出了运动系统中钢丝绳、导线、牵引走板与双滑车间的接触判别方法,通过罚函数法得到牵引走板与滑车的接触力,结合有限质点法形成了牵引走板与双滑车的碰撞接触计算方法,为输电线路工程施工提供了技术支持,可用于大高差、大转角等地形条件下牵引走板对双滑车通过性的精细计算.通过模型简化,可直接退化得到牵引走板与单滑车的接触计算方法.
针对输电线路施工中的典型工况,进行了牵引走板与滑车的接触过程分析,总结了运动过程中接触状态、滑车载荷、滑车摆动角度的变化规律.针对工程中关注的悬挂长度、走板厚度等因素,进行多个工况的走板过双滑车过程计算,分析了牵引力、最大接触力、滑车载荷、摆动幅度等主要影响参数的变化.计算分析结果为牵引走板、放线滑车的设计与施工安装提供了参考依据,为张力放线施工安全提供保障.
参 考 文 献
孟丽娟.输电线路张力架线施工技术探讨.电源技术应用,2013, 12:11~13 [百度学术]
Meng L J.Discussion on tension stringing construction technology of transmission line. Application of Power Supply Technology, 2013,12:11~13(in Chinese) [百度学术]
徐守琦,金连勇,孟昭清.二牵六放线施工工艺的实施.电力建设, 2010,31(8):38~42 [百度学术]
Xu S Q,Jin L Y,Meng Z Q. The implementation of two-drag-six construction technology. Electric Power Construction, 2010, 31(8):38~42(in Chinese) [百度学术]
季伟.对输电线路牵引绳走板施工可行性的探讨.华东科技(学术版), 2012, 12(1):315 [百度学术]
Ji W. Discussion on construction feasibility of traction rope walking plate of transmission line. East China Science & Technology, 2012, 12(1):315(in Chinese) [百度学术]
徐守琦,刘凯,易明阳.特高压双牵引走板的研制.电力建设, 2010, 31(3):102~105 [百度学术]
Xu S Q, Liu K, Yi M Y.Development of UHV double-towing-walking-plate.Electric Power Construction,2010, 31(3):102~105(in Chinese) [百度学术]
汤良友,鲁飞.±800kV特高压直流输电线路张力架线滑车悬挂施工工艺.电力建设,2011,32(4):121~124(Tang L Y,Lu F. Construction techniques of tension stringing pulley hanging for the ±800 kV UHVDC transmission line. Electric Power Construction,2011. 32(4):121~124(in Chinese)) [百度学术]
罗尧治,郑延丰,杨超等.结构复杂行为分析的有限质点法研究综述.工程力学,2014,31(8):1~7 [百度学术]
Luo Y Z,Zheng Y F,Yang C, et al. Review of the finite particle method for complex behaviors of structures. Engineering Mechanics, 2014, 31(8):1~7(in Chinese) [百度学术]
赵阳,彭涛,王震.基于向量式有限元的索杆张力结构施工成形分析.土木工程学报,2013,46(5):13~21 [百度学术]
Zhao Y,Peng T,Wang Z. Construction analysis of cable-strut tensile structures based on vector form intrinsic finite element.China Civil Engineering Journal, 2013,46(5):13~21(in Chinese) [百度学术]
Wu T Y. Dynamic nonlinear analysis of shell structures using a vector form intrinsic finite element.Engineering Structures,2013, 56(6):2028~2040 [百度学术]
梁立孚,郭庆勇.非保守非线性刚-热-弹耦合动力学的Lagrange方程.动力学与控制学报,2019,17(5):487~496 [百度学术]
Liang L F,Guo Q Y. The Lagrange equation of non conservative non-linear rigid-thermos-elastic coupling dynamics. Journal of Dynamics and Control,2019,17(5):487~496(in Chinese) [百度学术]
华卫江,章定国.柔性机器人系统碰撞动力学建模.机械工程学报,2008, 43(12): 222~228 [百度学术]
Hua W J, Zhang D G. Modeling of impact dynamics of flexible robots. Chinese Journal of Mechanical Engineering,2008, 43(12):222~228(in Chinese) [百度学术]
董富祥,洪嘉振.平面柔性多体系统正碰撞动力学建模理论研究,计算力学学报,2010,27(6):1042~1048 [百度学术]
Dong F X, Hong J Z. Study on the modeling theory of the normal impact dynamics for the planar flexible multibody system.Chinese Journal of Computational Mechanics, 2010,27(6):1042~1048(in Chinese) [百度学术]
Zheng J,An X,Huang M.GPU-based parallel algorithm for particle contact detection and its application in self-compacting concrete flow simulations. Computers & Structures, 2012, 112-113(DEC.):193~204 [百度学术]
俞锋,尹雄,罗尧治,等.考虑接触点摩擦的索滑移行为分析.工程力学,2017,34(8):42~50 [百度学术]
Yu F,Yin X,Luo Y Z,et al. Cable sliding analysis considering frictional effect. Engineering Mechanics,2017, 34(8): 42~50(in Chinese) [百度学术]
张鹏飞.结构破坏行为的数值模拟计算方法研究[博士学位论文]. 杭州:浙江大学, 2016(Zhang P F.Research on numerical methods for structural failure simulation[Ph.D Thesis].Hangzhou:Zhejiang Univesity, 2016(in Chinese)) [百度学术]
秦剑,乔良,张映晖,等.多档货运索道动力计算方法及结构冲击影响研究.动力学与控制学报,2020,18(2):59~68 [百度学术]
Qin J,Qiao L, Zhang Y H,et al. Calculation method for dynamics of multi⁃span freight cableway and influence of structure impact. Journal of Dynamics and Control,2020,18(2):59~68(in Chinese) [百度学术]
李君章.输电线路张力架线实用技术与计算.北京:中国电力出版社, 2018:43~64(Li J Z. Practical technology and calculation of tension stringing of transmission line.Beijing:China Electric Power Press,2018:43~64(in Chinese)) [百度学术]