, 温航1, 李正良1,21. 重庆大学 土木工程学院, 重庆 400045;
2. 重庆大学 山地城镇建设与新技术教育部重点实验室, 重庆 400045
收稿日期：2015-12-30
基金项目：国家自然科学基金资助项目 (51308568);中央高校基本科研业务费专项资金资助项目 (CDJZR12200061，CDJZR12200017)。
作者简介：肖正直 (1979-)，男，重庆开县人，重庆大学副教授;
李正良 (1963-)，男，江苏江阴人，重庆大学教授，博士生导师。
摘要：针对±800 kV悬索拉线塔立柱的静力稳定性问题, 首先基于广义变分原理推导斜拉线的等效弹性模量公式, 然后根据应变能相等得到了立柱的等效截面属性, 再根据悬索拉线塔立柱的受力特点, 将其简化为下端铰接、上端弹性支撑的等截面压杆, 最后基于压杆稳定理论得到拉线塔立柱的计算长度系数μ与初始预拉力和风速的关系.研究结果表明：背风侧立柱在垂直线路方向以及迎风侧立柱的静力计算长度系数可取值为1；风荷载及初始预拉力对背风侧立柱沿线路方向的静力计算长度系数有较大的影响，当初始预拉力较小时，其静力计算长度系数随风速的增大而迅速地增大，但当初始预拉力超过某值之后，背风侧立柱的静力计算长度系数受风速的影响较小.
关键词：悬索拉线塔静力稳定计算长度系数广义变分原理压杆稳定理论
Research on Static Stability of ±800kV Suspension Cable Guyed Tower
XIAO Zheng-zhi1,2
, WEN Hang1, LI Zheng-liang1,21. School of Civil Engineering, Chongqing University, Chongqing 400045, China;
2. Key Laboratory of New Technology for Construction of Cities in Mountain Area, Ministry of Education, Chongqing University, Chongqing 400045, China
Corresponding author: XIAO Zheng-zhi, E-mail: xiaozhengzhi@cqu.edu.cn
Abstract: To solve the problem of static stability of ±800 kV suspension cable guyed tower, equivalent elastic modulus formula of suspension cable is firstly derived based on the generalized variational principle. Then according to the strain energy, the equivalent column attributes of equivalent cross section are obtained. According to the force characteristics of suspension cable guyed tower column, it is simplified to the bar with uniform cross section which is hinged at the bottom and is supported elastically at the top. Finally, based on compressive bar stability theory, the relationship of calculation length coefficient of guyed tower column is got, the initial pretension and wind speed. The study results show that the static calculated length coefficient of column in windward side and column in leeward side along vertical transmission line direction is one, wind load and initial pretension have great impact on static calculated length coefficient of column in leeward side along transmission line, which rapidly increases with an increase in wind speed when initial pretension is lesser, but wind speed has a minor effect on it when initial pretension exceeds a certain value.
Key Words: suspension cable guyed towerstatic stabilitycalculated length coefficientgeneralized variational principlecompressive bar stability theory
悬索拉线塔的侧向刚度主要由拉线提供，拉线在力学上主要表现为高度柔性和极强的几何非线性，因此吸引了许多研究者对其进行研究.拉线初始预张力是拉线塔的重要参数，不少国外学者[1-4]对拉线初始预拉力大小对拉线塔静动力响应的影响进行了研究，但很少有国内外学者对拉线塔或拉线的简化分析进行研究.Ernst[5]提出通过修正材料的弹性模量来考虑索垂度引起的几何非线性, 但此方法的精度随着拉线的松弛和倾角的增大而降低.文献[6]采用一致切线刚度推导拉线的几何刚度，在Ernst公式的基础上通过引入一个修正因子来考虑由垂度引起的几何非线性，但该修正因子是通过拟合获得的经验公式.刘树堂[7]基于拉线的状态方程推导拉线的等效弹性模量.洪显诚[8]从能量出发，基于广义变分原理导出普遍适用的斜拉索等效弹性模量公式.
悬索拉线塔立柱为格构式构件，其分析模型主要有空间桁架模型和等效梁柱模型；Kahla[9]使用单位荷载法推导格构式构件的等效截面属性；Kollbrunner等[10]基于能量等效推导了桅杆的等效截面属性；由于上述等效方法均较为繁琐，H.Meshmesha等[11]提出一种较为简化的等效方法，该方法是基于应变能等效推导桅杆的等效截面属性，并从静力、模态、动力等多个角度验证这种等效方法的精度.
本文先基于广义变分原理推导拉线的等效弹性模量公式，在推导的过程中参考了文献[8]提供的方法，但在该方法的基础上做了适当改进；然后根据应变能相等得到了立柱的等效截面属性，再根据悬索拉线塔立柱的受力特点，将其简化为下端铰接、上端弹性支撑的压杆，最后基于压杆稳定理论推导了拉线塔立柱的计算长度系数与相对刚度之间的关系表达式，在此基础上得到拉线塔立柱的计算长度系数与初始预拉力和风速的关系曲线.
1 拉线等效弹性模量公式的推导基本假定:①只考虑斜拉线承受拉力, 即拉索是柔性的; ②大位移小应变, 即考虑斜拉线的几何非线性, 但材料处于弹性工作阶段.
图 1a为一平面斜拉线, 容重γ, 面积为A, 沿其弦向的均布荷载为q, 其中X, Z方向的分量分别为qx和qz; 微段变形前后的拉力分别为T0和T; 设u(x, z), w(x, z) 分别为微段沿X, Z方向的位移, 变形后C点移动到C′点; 用O, C端表示斜拉线的端部, 其边界条件为
 | (1) |
 | 图 1 斜拉线分板模型Fig.1 Analysis model of inclined cable (a)—平面斜拉索；(b)—斜拉索微段. |
选取斜拉线的微段如图 1b所示, A, B为微段变形前的端点, A(s), B(s+ds) 为微段变形后的端点, s为坐标参数; 设变形前A点的坐标为 (x, z)，B点的坐标为 (x+dx, z+dz);
斜拉线处于弹性工作阶段, 由胡克定律有
 | (2) |
 | (3) |
 | (4) |
 | (5) |
 | (6) |
 | (7) |
 | (8) |
对式 (8) 中的第三项进行积分得
 | (9) |
 | (10) |
 | (11) |
 | (12) |
 | (13) |
 | (14) |
 | (15) |
 | (16) |
 | (17) |
 | (18) |
图 2(Fig. 2)
 | 图 2 拉线塔立柱Fig.2 Stand column on the guyed tower (a)—简化模型；(b)—分析模型. |
在风荷载作用下，迎风侧的拉线将会被张紧，其垂度引起的几何非线性可以忽略，此时拉线给立柱提供非常大的抗侧刚度，立柱上端的支承可以视为铰接.背风侧的拉线将会出现松弛现象，但双柱悬索拉线塔的双柱通过悬索及绝缘子协同工作，因此背风侧立柱在垂直线路方向的支承也可简化为铰接.背风侧立柱在沿线路方向的抗侧刚度主要由松弛的拉线来提供，因此将上端边界视为弹性支座，其刚度值K的表达式需结合图 3进行推导.
图 3(Fig. 3)
 | 图 3 背风侧拉线塔立柱的简化模型Fig.3 Simplified model of stand column on the guyed towerin leeward side |
设PA长度为L1, 等效弹性模量为EeqPA; 设PB长度为L2, 等效弹性模量为EeqPB; 假设P点在垂直OP的外力F1作用下产生Δ的位移, 则PA, PB内力改变量的合力在垂直OP方向的分量等于外力F; 即
 | (19) |
 | (20) |
3 立柱的计算长度系数μ与相对刚度KL3/EI之间的关系背风侧立柱在沿线路方向的分析模型如图 2b所示, 设拉线塔立柱的屈曲模态为y(x); 对屈曲后的任意截面列力矩平衡方程, 则
 | (21) |
 | (22) |
 | (23) |
当x=0时, y(x)=0;y″(x)=0, 即
 | (24) |
 | (25) |
 | (26) |
当P1≤P2, 即
则拉线塔立柱计算长度系数μ:
 | (27) |
 | (28) |
 | 图 4 计算长度系数与相对刚度的关系曲线Fig.4 Relationship curve between calculation length coefficient and relative rigidity |
4 算例分析本文算例分析依托于郑州-哈密±800 kV直流悬索拉线塔, 如图 5所示.悬索拉线塔立柱属于格构式构件, 如图 6所示, 截面尺寸为2 m×2 m, 由不同型号的角钢组成格构式塔架.文献[11]先根据应变能相等将格构式结构的缀条等效为薄板, 然后求出等效截面属性; 本文按照这种等效方法求得A, B区段的等效截面属性为
 |
 | 图 5 悬索拉线塔示意图Fig.5 Diagram of suspension cable tower |
图 6(Fig. 6)
 | 图 6 悬索拉线塔立柱Fig.6 Stand column of suspension guyed tower |
由于A, B区段的抗弯刚度差距不大, 按照区段的高度取加权平均值为3.112 75×109 N·m2.
为了研究拉线的初始预张力和风速对双柱悬索拉线塔立柱静力稳定性的影响, 分别取拉线极限承载力的5%~40%, 增量取5%;风速为30, 33, 35, 40, 42 m/s进行研究.
双柱悬索拉线塔单根拉线的容许拉应力取780 N/mm2, 单根拉线的计算截面面积为464.95 mm2, 拉线的弹性模量为185 kN/mm2; 每单位长度的质量为3.87 kg; 预拉力对背风侧立柱计算长度的影响如图 7所示, 风荷载对背风侧立柱计算长度系数的影响如图 8所示.
图 7(Fig. 7)
 | 图 7 计算长度系数与初始预拉力的关系曲线Fig.7 Relationship curve between calculation length coefficient and initial pretension |
图 8(Fig. 8)
 | 图 8 计算长度系数与初始预拉力的关系曲线Fig.8 Relationship curve between calculation length coefficient and wind speed |
由图 7, 图 8可知, 当风速一定时, 背风侧立柱的计算长度系数随拉线初始预拉力的增大而减小, 但当拉线初始预拉力超过某数值时, 计算长度系数始终为1;当拉线初始预拉力一定时, 风速越大, 背风侧立柱的计算长度系数越大.
5 结论1) 迎风侧拉线在风荷载作用下将被拉紧, 为立柱提供较大的侧向刚度; 由拉线塔立柱临界刚度计算表达式可知, 迎风侧立柱在自重、预拉力及风荷载共同作用下的计算长度系数为1.
2) 背风侧拉线在风荷载作用下将会松弛, 由于双柱悬索拉线的双柱在悬索作用下协同工作, 背风侧立柱在垂直线路方向的计算长度系数为1;但在沿线路方向的计算长度系数受到风荷载大小及预拉力大小的影响.
3) 当风速一定时, 背风侧立柱在垂直线路方向的计算长度系数随初始预拉力的增大而减小; 当预拉力一定时, 背风侧立柱在垂直线路方向的计算长度系数随风速的增大而增大.
参考文献
| [1] | Sparling B F.The dynamic behavior of guys and guyed masts in turbulent winds[D]. Ontario:University of Western Ontario, 1995.http://www.nrcresearchpress.com/doi/full/10.1139/l10-130 |
| [2] | Wahba M M Y, Monforton G. Dynamic response of guyed masts[J].Engineering Structures, 1998, 20(12): 1097–1101.DOI:10.1016/S0141-0296(97)00206-X |
| [3] | Ballaben J, Guzman M, Rosales M. Parametric studies of guyed towers under wind and seismic loads[J].Asociación Argentina de Mecánica Computacional, 2011, 30: 1019–1032. |
| [4] | Carrasco-Luzardo A, Parnas V E, Martin-Rodriguez P. Guy tension influence on the structural behavior of a guyed mast[J].International Association for Shell and Spatial Structures, 2012, 53(2): 111–116. |
| [5] | Ernst H J. Der E-modul von seilen unter berucksiehtigung des durchanges[J].Battingenieur, 1965, 40(2): 52–55. |
| [6] | Kiureghian A D, Sackman J L. Tangent geometric stiffness of inclined cables under self-weight[J].ASCE, 2005, 131: 941–945.DOI:10.1061/(ASCE)0733-9445(2005)131:6(941) |
| [7] | 刘树堂. 建立斜拉线单元刚度矩阵的等效线性方法[J].广州大学学报, 2006, 5(5): 77–81. ( Liu Shu-tang. Equivalent linearization method for setting up stayed cable element stiffness matrices[J].Journal of Guangzhou University, 2006, 5(5): 77–81.) |
| [8] | 洪显诚. 精确的斜拉索等效弹性模量公式的推导[C]//全国桥梁结构学术大会. 上海: 同济大学出版社, 1992: 834-838. ( Hong Xian-cheng.Deduction of exact formula of the inclined cable's equivalent elasticity modulus[C]//The National Bridge Structure Academic Conference.Shanghai:Tongji University Press, 1992:834-838.) |
| [9] | Kahla N B. Equivalent beam-column analysis of guyed towers[J].Computers & Strucutre, 1995, 55(4): 631–645. |
| [10] | Kollbrunner C F, Basler K. Torsion in structures[M]. Berlin: Springer-Verlag, 2013: 10-94. |
| [11] | Meshmesha H, Sennah K, Kennedy J B. Simple method for static and dynamic analysis of guyed towers[J].Strucutral Engineering and Mechanics, 2006, 23(6): 635–649.DOI:10.12989/sem.2006.23.6.635 |
