2. 重庆大学 山地城镇建设与新技术教育部重点实验室, 重庆 400045
2. Key Laboratory of New Technology for Construction of Cities in Mountain Area, Ministry of Education, Chongqing University, Chongqing 400045, China
悬索拉线塔的侧向刚度主要由拉线提供,拉线在力学上主要表现为高度柔性和极强的几何非线性,因此吸引了许多研究者对其进行研究.拉线初始预张力是拉线塔的重要参数,不少国外学者[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) |
选取斜拉线的微段如图 1b所示, A, B为微段变形前的端点, A(s), B(s+ds) 为微段变形后的端点, s为坐标参数; 设变形前A点的坐标为 (x, z),B点的坐标为 (x+dx, z+dz);
斜拉线处于弹性工作阶段, 由胡克定律有
(2) |
其中:ε为变形前后微段应变增量, 其表达式为
(3) |
变形前后微段应变能的增量为
(4) |
斜拉索变形前后外荷载px和pz及支座反力做功引起势能的改变量为
(5) |
其中:TCx, TCz分别为C端沿弦向与垂直弦向的支反力; 则微段变形前后的总势能泛函为
(6) |
其中, 应变ε增量受到式 (3) 的制约, 即式 (6) 取极值时属于条件极值, 此处用斜拉索变形后的拉力作为拉格朗日乘子, 即斜拉线的广义总势能泛函为
(7) |
对式 (7) 求变分并结合分部积分, 同时略去高阶微量得
(8) |
式中,H为微段变形后拉力T的水平分量.
对式 (8) 中的第三项进行积分得
(9) |
对式 (8) 中的第四项进行积分并代入边界条件得
(10) |
由式 (10) 可知w(x) 是H(0) 的函数,即要得到w(x),则必须先求得H(0).下面利用斜拉索的广义最小余能泛函来求H(0), 表达式为
(11) |
对式 (11) 求变分并略去高阶微量得
(12) |
对式 (11) 求导代入式 (12) 并积分得
(13) |
由于支座C端垂直于弦向的位移对斜拉线拉力变化的影响不大, 所以只考虑斜拉线沿其弦向的位移; 由力的分解与合成可知荷载qx和qz的表达式为
(14) |
将式 (14) 代入式 (13) 并积分得
(15) |
令等效弹性模量等于沿弦向应力改变量与应变改变量的比值, 即
(16) |
将式 (15) 代入式 (16) 并整理得
(17) |
式中:E为拉线的弹性模量; L为拉线弦向的长度; α为拉弦弦向与水平方向的夹角; A为拉线的截面面积, 其余参数的表达式为
(18) |
根据悬索拉线塔立柱的受力特点,将其简化为下端为铰接, 上端为弹性支座的截面压杆,立柱的等效截面属性根据文献[11]提供的方法得到,其简化模型见图 2a,EI为立柱的等效抗弯刚度.
在风荷载作用下,迎风侧的拉线将会被张紧,其垂度引起的几何非线性可以忽略,此时拉线给立柱提供非常大的抗侧刚度,立柱上端的支承可以视为铰接.背风侧的拉线将会出现松弛现象,但双柱悬索拉线塔的双柱通过悬索及绝缘子协同工作,因此背风侧立柱在垂直线路方向的支承也可简化为铰接.背风侧立柱在沿线路方向的抗侧刚度主要由松弛的拉线来提供,因此将上端边界视为弹性支座,其刚度值K的表达式需结合图 3进行推导.
设PA长度为L1, 等效弹性模量为EeqPA; 设PB长度为L2, 等效弹性模量为EeqPB; 假设P点在垂直OP的外力F1作用下产生Δ的位移, 则PA, PB内力改变量的合力在垂直OP方向的分量等于外力F; 即
(19) |
则背风侧立柱在沿线路方向的刚度K为
(20) |
式中:θ1, θ2分别为拉线AP, BP与水平方向的夹角.
3 立柱的计算长度系数μ与相对刚度KL3/EI之间的关系背风侧立柱在沿线路方向的分析模型如图 2b所示, 设拉线塔立柱的屈曲模态为y(x); 对屈曲后的任意截面列力矩平衡方程, 则
(21) |
式 (21) 对x求两次导得
(22) |
式中,
(23) |
由立柱的边界条件确定A, B, C, D;
当x=0时, y(x)=0;y″(x)=0, 即
(24) |
当x=L时, y″(x)=0;
(25) |
若A,C同时为零, 则方程只有零解, 则拉线塔立柱不会失稳; 所以A, C不可能同时为零, 即
(26) |
解式 (26) 得:
当P1≤P2, 即
则拉线塔立柱计算长度系数μ:
(27) |
拉线塔立柱的计算长度系数μ与相对刚度KL3/EI的关系如图 4所示.由图 4可知, 当相对刚度接近10时, 拉线塔立柱的计算长度系数μ等于1, 即立柱上端的边界条件接近铰接; 此时的临界刚度Kcr为
(28) |
本文算例分析依托于郑州-哈密±800 kV直流悬索拉线塔, 如图 5所示.悬索拉线塔立柱属于格构式构件, 如图 6所示, 截面尺寸为2 m×2 m, 由不同型号的角钢组成格构式塔架.文献[11]先根据应变能相等将格构式结构的缀条等效为薄板, 然后求出等效截面属性; 本文按照这种等效方法求得A, B区段的等效截面属性为
由于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, 图 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 |