搜索
您的当前位置:首页正文

窄基塔随机风致振动响应及风振系数研究

来源:榕意旅游网
第29卷第1期 2013年1月 科技通报 BULLET【N 0F SCIENCE AND TECHN0L0GY V01.29 No.1 Jan.2013 窄基塔随机风致振动响应及风振系数研究 潘 峰 ,陈稼苗 ,聂建波 ,叶 尹 ,沈建国 (1.浙江省电力设计院,杭卅I 310012;2.浙江省电力物资供应公司,杭州310003) 摘 要:小根开塔(窄基塔)是输电线路一种新塔型,属于风敏感结构,风与结构的相互作用十分复杂,风 荷载常常是设计的主要控制荷载。本文以某一220 kV同塔双回路直线窄基塔为工程背景,建立了详细 的有限元计算模型,获得了结构的自振动力特性。在考虑节点风荷载空间相关性的基础上,利用Kaimal 谱对结构的风荷载进行了准确的数值模拟。结合动力时程积分法,计算得到了结构的风致振动响应时 程。研究了结构位移平均值、位移均方根值和加速度均方根值的分布特点,同时分析了结构各层的荷载 风振系数和位移风振系数,并与我国规范进行了对比。通过研究,揭示了窄基塔的风致振动特性,结果 可作为窄基塔结构抗风设计的参考。 关键词:窄基塔;有限元模型;风荷栽模拟:时程分析;风致响应;风振系数 中图分类号:TU311 文献标识码:A 文章编号:1001—7119(2013)01—0086—07 Random Wind-Induced Dynamic Response and Vibration Coeficifent of Narrow Base Transmission Tower Pan Feng1,Chen Jiamiao ,Nie Jianbo ,Ye Yin ,Shen Jinguo a(1.Zhejiang Electirc Power Design Institute,Hangzhou 310012,China; 2.Zhejiang Provincial Electirc Equipment&Materials Supply Company,Hangzhou 3 10003,China) Abstract:Small root span tower(narrow base transmission tower)is a new tower type of the transmission line.With the complex interaction between wind and structure.narrow base transmission tower belongs to wind—sensitive structures and the wind load becomes the main problem of the structure design.Based on the engineering background of a double—circuit 220kV narrow base transmission tower,a precise finite element model Was established and the vibration characteristics of structure have been obtained.The time history of lfuctuating wind load was accurately simulated by considering the spatial correlation of Kaimal wind spectrum.Combination of dynamic time-integration method,the wind—induced vibration responses of the tower were obtained in the time domain.The characteristics of MeaR displacement、Rms displacement and Rms acceleration at different layer of transmission tower have been analyzed in a detailed.Meanwhile,the wind—induced vibration coefficient at diferent layer of transmission tower were also studied and compared,respectively.Through research,some characteristics of wind—induced vibration of narrow base transmission tower were revealed and the procedure and the results presented in this paper can beusedasreferenceforthe desin gofnarrowbasetransmissiontower. Key words:narrow base transmission tower;finite element model;wind load simulation;time domain analysis; wind——induced dynamic response;wind——induced vibration c ̄fficient; 收稿日期:2011-09—20 作者简介:潘峰(1980一),男,浙江绍兴人,博士,主要从事输电塔结构分析研究,E-mail:qiushif@163.con 第1期 潘峰等蹿基塔随机风致振动响应及风振系数研究 87 0 引言 量, × 。 随着我国经济飞速发展,根据建设“统一坚 强智能电网”的战略部署,我国东部地区加快了 电网建设力度。输电线路杆塔作为电网建设的重 要组成部分,其设计直接关系到电网建设的安 全、经济和电网运行的稳定。 。鼬 爻 × 苗 .d< 爻 ×  ̄0110.∞ × 餐 7蝴 受 )( ; 。d< 爻 >b。 但随着城市建设的高速发展,土地资源越来 萎:~ 950 越稀缺,人们节约、环保的观念也不断增强,尤其 对于输电线路在市区通过遇到了越来越大的困 难。而小根开塔(以下简称窄基塔)可以设立在城 市道路的绿化带中,可以较大程度上避免因线路 经过引起的纠纷,同时占地面积大大小于传统铁 塔,对节约城市的土地资源意义很大。 研究发现窄基塔具有柔性大、阻尼小等特 点,是一种风敏感性结构,风荷载常常是设计的 主要甚至是决定性的控制荷载。而我国规范对窄 基塔的抗风设计参数并不十分完善,这对设计和 研究人员提出了更高更新的要求。近些年,由于 大风引起的铁塔倒塌事故时有发生。据不完全统 计,2005年我国由于强风引起的输电塔倒塌就有 75基[1]。在输电塔结构抗风设计中,风振系数是结 构设计的关键数据。至今,对于在强风作用下输 电塔的风荷载及其风振响应尚无法从理论上得 到统一的模型。如何快速、可靠地获得输电塔结 构的动力特性及其在强风作用下的风致振动响 应,是目前输电塔结构研究的重要方向。 本文将以某一220 kV同塔双回路直线窄基 塔为工程背景,建立了详细的有限元计算模型, 研究了结构的动力特性。并在时域内对窄基塔结 构的风致响应和风振系数进行了详细的研究和 分析,得出了一些可供设计参考的结论。 1 结构有限元模型及动力特性研究 该输电跨越塔塔高48.75 m,呼高33 m,塔身 平面形状为正方形,基底根开约3.0 m,全塔结构 的单线图及三维图如图1所示。定义沿横担方向 为 向(横风向),垂直横担方向为Y向(/1 ̄风向), 高度方向为Z-向。全塔的主立杆采用高强钢管,铁 塔其它杆件采用角钢。这样就可以最大限度发挥 钢管、角钢的优势。 量 爻 × 萎 × : × i × ; 量 × × : -× _一;  5101. 5100 豢 × { l 图1窄基塔外形图与坐标系统 Fig.1 The configuration and coordinate system of narrow base transmission tower 在有限元计算中,采用两种计算模式:一种为 梁系计算模式(梁系模型),钢管采用空间梁单元 模拟,角钢采用空间杆单元模拟;另一种为全桁 架计算模式(杆系模型),所有杆件均采用空间杆 单元模拟。而输电塔的中间爬梯、节点板及其它 连接构件对结构刚度影响较小,在建模时仅考虑 其质量,以质量单元来模拟。 由于输电塔结构自振频率较为稀疏,因此重 点考虑输电塔的前六阶频率 采用ANSYS通用 有限元计算软件对窄基塔结构进行建模,钢材的 弹性模量取E=2.06×10”N/m ,密度取p=7850 kg/m。,泊松比取/z=O.3,利用分块Lanczos算法,获 得了窄基塔的前六阶频率,如表1所示,图2为 结构的模态示意图。从表1中可以得出,采用两 种不同的计算模式得到的窄基塔的前六阶频率 是相当接近的,说明结构建模计算是正确的。同 时可以发现,窄基塔的前三阶模态分别是x向一 阶弯曲、Y向一阶弯曲和一阶扭转振型,并且x向 与Y向弯曲振型的频率非常接近。 88 表1窄基塔前六阶频率 科技通报 第28卷 经基本得到解决。 Table 1 The first six natural frequencies of transmission tower 第一阶 第二阶 第三阶 第四阶 第五阶第六阶 图2输电塔前六阶模态 Fig.2 The ifsrt six modals of narrow base transmission tower 2窄基塔风荷载模拟 结构抗风分析法主要有频域法和时域法。而 对于像输电塔一类的重要结构,由于频域法需要 进行结构数学模型的简化并且无法考虑结构非 线形、时间相关性等因素,实际计算往往采用时 域法。风荷载的准确模拟是输电塔结构进行时域 风致响应分析的第一步,而风荷载模拟主要是结 构风场的准确模拟。 目前,脉动风速数值模拟方法主要有线性回 归滤波法[21和谐波合成法 。前者基于线性滤波 技术,主要有AR(autoregressive)法、MA(moving average)法和ARMA (autoregressive moving average)法;后者基于三角级数求和,主要有 CAWS(constant amplitude wave superposition)法和 WAWS(weighted maplitude wave superposition)法; AR法和ARMA法运算量相对较少,但模型精度 随各处风谱差异的增大而降低。CAWS法和 WAWS法运算精度较高,但为保证风场的各态历 经性要在每个频率上作大量运算,故计算效率较 低,但随着计算机技术的快速发展,这个问题已 本文采用谐波合成法中的经典方法 Shinozuka—Deodatis法。该法由Deodatis【4]于1996 提出,目前认为是谐波合成法中精度较高、稳定 性较强的一种方法,方法中引入FFT变换技术和 频率双索引(double—indexing ̄equency)的概念, 极大地提高了谐波合成法的模拟精度,可以模拟 各态历经的多变量平稳高斯随机过程。输电塔的 10 m高度处5O年重现期10 min平均风速定为 30 m/s,基本风压为0.5625 kN/mz,水平顺风向脉 动风速谱采用随高度变化的Kaimal谱,其表达式 为: s“(n)=2 丽 (1) 式中为 。标准高度10 m处的平均风速;rt为频 率;K为表面阻力系数,取值在0.003—0.03之间; v(z)为 高度处的平均风速。 脉动风的相干函数的表达形式为: Coh(r, )= r 2n expL—— ] (z)+ ( ) (2) 式中,c :分别为 ,y, 方向的衰减系数,取值 为:c =6,Cy=16,cz=lO。 计算的主要参数如表2所示。输电塔从底部 至顶部的每一节间定义为计算层,共有33层。利 用Shinozuka—Deodatis法模拟得到了所有层的顺 风向脉动风速时程。图3为采用Kaimal谱得到的 输电塔顶层的脉动风速时程曲线图。图4为相应 的风速功率谱密度图。从图4中可以得出, Shinozuka—Deodatis法模拟得到的脉动风速时程 的功率谱与目标谱非常吻合,说明模拟得到的结 果精度较高。 表2风速时程模拟时的主要参数 Table 2 The main parameters in wind speed simulation of narrow ba岛e transmissi0n tower 主要参数名称 参数值 参考高度处风J/Um/s 30.0 地形地貌类型 B类 卡门常数 0.4 地面粗糙长度Z0/m 0.05 摩擦风速u*/ngs 2.03 模拟时间步长/s 0.10 模拟风速时程时间长度/s 400 第1期 潘峰等蹿基塔随机风致振动响应及风振系数研究 89 嘲 于是,作用与输电塔结构某一层的风荷载可 由下式得到: P(z,t)=P(z)+p(z,t) ,。2。 ——亩—— —I{ — —— — — — o = ( (z)v。( ) ,( )A(z)V(z)u(z,£)(3) 二 式中,P( )、P( , )为结构某一层的平均风荷载和 脉动风荷载; ( )为结构的体型系数,本文按照 杆塔规范钢管取为0.85,角钢取为1.3;A(z)为风 荷载作用面积;p为空气密度; (z)为结构某一层 高度的平均风速,u(z, )为相应的脉动风速。 3窄基塔随机风致响应研究 利用动力时程积分法对窄基塔结构进行风致 Fr日q z) 动力响应分析,得到了结构各层的位移、速度和 图4输电塔顶层的风速功率谱密度图 Fig.4 Comparison between the target and simulated PSD of 加速度等信息。图5为采用Kaimal谱计算得到的 输电塔第80(塔顶)、300(上横担挂线点)、5100 (塔底)号节点的位移和加速度时程曲线图。图6 lfuctuating wind speed at top of tower ●  ’■ 下 '_ 图5窄基塔不同层的位移和加速度时程曲线图 Fig.5 Time history of displacement and acceleration at diferent layer of narrow base transmission tower ∥ ; I—。e—v期 崖均  l图6窄基塔各层位移平均值、位移均方根值和加速度均方根值的对比示意图 Fig.6 Variation of Mean displacement、Rms displacement and Rms acceleration at different layer of naITOW base tower 科技通报 第28卷 为输电塔各层位移平均值、位移均方根值和加速 度均方根值的对比示意图。 从图中可以发现,塔顶和横担挂线点处的振 动明显较强,且处在同一量级,塔顶处振动略大 于横担挂线点;沿塔身高度的位移和加速度平均 动力风荷载。目前,最常用的两种风振系数是荷 载风振系数 和位移风振系数 。 计算荷载风振系数时,需得到铁塔各节点的 加速度响应,公式可表示为: 31L=1+g (5) 值和均方根值的分布基本与前二阶振型相似。 图7为窄基塔第8o(塔顶)、3oo)上横担挂线 点)、950(边坡点)、5100(塔底)号节点的加速度功 式中,m 是节点对应的集中质量; : 为节点的加 率谱密度示意图。从图中可以发现,对于窄基塔 的上部,结构的振动以第一、二阶弯曲模态响应 为主,第三阶扭转模态响应的贡献相对较小;而 对于窄基塔的底部,结构的振动仍以第一、二阶 速度均方根值; 为节点的体型系数; 为节点 的风压高度变化系数;OJo=V 02/1600为基本风压; A 为节点对应的受风面积;g为峰值因子,本文取 3.Of ol,对应的保证率为98.65%。 位移风振系数 为结构最大位移反应与平 均位移反应之比值,本质上等同于Davenport提 出的阵风荷载因子法。由于位移响应已经包含了 弯曲模态响应占主导,同时部分第三阶扭转模态 响应被激起。因此,输电塔结构的风致响应需考 虑前三阶模态的贡献。 4窄基塔风振系数研究 风振系数 是指在一定时间范围内由平均风 风致响应的背景分量和共振分量,因此概念上比 较清晰。平均风载及等效静风载作用下构架结构 的静力平衡方程为: [K]{U}=JR],[K]{U+gtru}_[Pe] (6) 压和脉动风压共同作用的总效应(位移、内力、内 矩、荷载)与平均风压产生的效应之比,定义为: : :一式中, 为平均位移,(rⅡ为位移均方根,K为结构 总刚度矩阵。 于是,风振系数列阵{ }可表示为: P+Pd一二:1+ P (4) P P { }- i ‘ }:1+[ 3i } [K]{ } (7) [K]{一U} 式中, 为等效静风荷载,P为静力风荷载,Pd为 140。 12。0 I12o。 ‘目‘塔80号节点 向加进度功率昔1 1。o。 ——窜基蟮3∞峙节煮v向施速度功率昔 辜∞0o 量e0。 V ^ 80O 重 邑 6。。 6o。 惜 麒4o。 盏 星 。0 20。 瑙 异 2(30 。 J 窜基塔950号节点 ̄f ̄30速度功率谮 』 LL . ^ 16。 14o 毒 00 量 锝 臼。 菩∞ 譬4。 2。  lI』 L .  一1J 、^^ | J。 J . 圈7窄基塔不同层的加速度功率谱密度图 Fig.7 The power spectral density of acceleration at different layer of nalTow base transmission tower 铝第1期 对于节点i,可得: 潘峰等席基塔随机风致振动响应及风振系数研究 91 式中, 为脉动增大系数;s.为风压脉动和风压高 ∑k ̄go" 十 :度变化等的影响系数; 为振型、结构外形的影 (8) 1+ 响系数。 ∑k + 户 ≠ 按照我国《建筑结构荷载规范》(GB 50009— 2001)同和《建筑结构荷载规范))(GBJ9—87)t81的定 义,结构在名高度处的风振系数可按下式计算: :式中,为i节点的相邻节点, 、 为总刚度矩阵 . 第i行的对角元和非对角元。定义比例参数: 1+鱼 (12) j-=lj#i ∑k ̄go" , 2=∑k j=t/ei 1+哗 g=1+ (9) (9) 式中, 为脉动增大系数; 为脉动影响系数; 为 振型系数; 为风压高度变化系数。 而《架空送电线路杆塔结构设计技术规定》 (DLT5154—2002)[91中只给出全塔高度小于6O m 的风振系数,当高度大于60 131时,风振系数按照 则第 节点的位移风振系数 可以表示为: o'./’ O ̄2+kii _, ~-, 写成标准形式为: =l+ =l+g (10) 式中,Y 、 分别为节点i的位移响应的平均值和 均方根值;叼为修正系数;其它定义与 的相同。 般而言,各节点位移平均值比均方根来得大, 且往往各节点的均方根与平均值成同一比例,此 时OL 一 ,因此修正系数叩一般取为1.0。 按照《高耸结构设计规范))(GB50135—2006) ̄ 中的定义,自立式高耸结构在 高度处的风振系 一数可按下式确定: =l+ge】 2 表3窄基塔沿高度的风振系数比较 (11) 《建筑结构荷载规范》计算,由上到下逐渐减小,并 且按高度加权平均值不低于1.60。 表3列出了输电塔结构沿高度荷载风振系数 和位移风振系数 的分布,并与《高耸结构设 计规范》取值、《建筑结构荷载规范》取值进行了 比较。图8为相应的风振系数分布图。从中可以 发现,计算得到的窄基塔结构的荷载风振系数在 沿高度上的数值分布不均匀,底部小,上部大;而 位移风振系数在在沿高度上的数值基本一致,如 果全塔风振系数取同一数值时建议采用位移风 振系数作为设计参考较为合适。采用《高耸结构 设计规范》和《建筑结构荷载规范》计算得到的风 Table 3 Comparison of wind-induced vibration coefficient at different height of na ̄ow base transmission tower 高度,m 荷载风 位移风 高耸设 荷载规 荷载规 振系数 振系数 计规范 范2006 范1987 图8窄基塔沿高度的风振系数分布示意图 Fig.8 Distribution of wind-induced vibration coefficient at different height of na/Tow base transmission tower 科技通报 第28卷 振系数规律是从底部至塔身数值逐渐增大,总的 趋势也是底部小,上部大。值得注意的是按照 如果全塔风振系数取同一数值时建议采用位移 风振系数作为设计参考较为合适。 (4)综合时域法和我国规范计算得到风振系 数,设计时风振系数塔身顶部取值1.96,塔身上 部取1.85,塔身中部取1.7O,塔底部取1.60;地线 1987版《建筑结构荷载规范》计算得到的风振系 数略小于2006版,原因是2006版计算时采用的 结构一阶频率取自有限元分析结果为1.057 Hz, 而1987版计算时结构一阶频率按照规范中的经 验公式得到,约为0.6 Hz,故风振系数略小。 最终窄基塔设计时,塔身的风振系数按照图 架和上、中、下横担同样可按照上述方法分析,分 别取值为2.30、2.20、2.10和2.00。 参考文献: 【l】张勇.输电线路风灾防御的现状与对策IJ】.华东电力. 2006,34(3):210-213. 【2】 Iannuzzi A,Spinelli P.Artiifcial wind generation and 8中的粗实线选取,而地线架和上、中、下横担同 样可按照上述方法分析,分别取值为2.30、2.20、 2.10和2.00。 5 结论 本文以某一220 kV同塔双回路小根开杆塔 structural response.Joumal of Stmcturla Engineering【J】. ASCE 1987,113(12):2382—2398. [3】 Shinozuka M,Jan CM.Digital simulation of random processes and its applications[J].Journal of Sound and 为工程背景,研究了结构的动力特性,并在时域 内对结构的风致响应和风振系数进行了详细的 研究,可以得出以下结论: Vibration 1972,25(1):1 1 1-128. [4]Deodatis G.Simulation of ergodic multivariate stochastic processes[J1.Journal of Engineering Mechanics,ASCE 1996;122(8):778—787. (1)窄基塔结构有限元分析采用梁系模式和 [5]张相庭,结构风压和风振计算【M】_上海:同济大学出版 社出版,1985. 杆系模式计算得到的结构自振频率非常接近;模 态振型分别是 向弯曲、Y向和扭转振型,且 向 与Y向弯曲振型的频率非常接近。 【6】GBJ 135—90,高耸结构设计规范【s】.北京:中国建筑工业 出版社,1991. [7】GB50009—2001,建筑结构荷载规范[s】.北京:中国建筑工 业出版社.2002. (2)窄基塔的上部,结构的振动以弯曲模态 响应为主,扭转模态响应的贡献几乎没有;而对 于窄基塔的底部,结构的振动以前两阶弯曲模态 占主导,扭转模态响应为次要。因此,输电塔结构 的风致响应需考虑前三阶模态的贡献。 (3)计算得到的输电塔结构的荷载风振系数 在沿高度上的数值分布不均匀,底部小,上部大; 而位移风振系数在在沿高度上的数值基本一致, [8】GBJ 9-87,建筑结构荷载规范[s】.北京:中国建筑工业出 版社。1987. [9】DLT5154—2002,架空送电线路杆塔结构设计技术规定 [s].北京:中国电力出版社,2002. 【10]张建胜,武岳,沈世钊.结构风振极值分析中的峰值因 子取值探讨[J】.铁道科学与工程学报,2007,4(1):28— 32. (上接第85页) [3】陈祥桥,涂佳黄,周岱等.某大型机场航站楼结构风荷载 Industirla Aerodynamics.201 1(99):573—583. 体型系数的数值模拟[J】,上海交通大学学报,2011,45 (4):517—522. 【6】燕辉,楼文娟,唐锦春,等.萧山国际商务中心风载风洞 试验研究【J].空气动力学报,2004,22(2):225—230. [7]侯家健,韩小雷,谢壮宁.复杂连体双塔高层建筑的风荷 载特性『J1.华南理工大学学报(自然科学版),2008,36 (3):121—126. I4】Wonsul Kim,Yukio Tamura,Akihito Yoshida.Interference effects on local peak pressure between two buildings[J】, Journal of Wind Engineering and Industiral Aerodynamics. 201 1(99):584-600. 【5]K M Lam,J G Zhao,M Y H Leung.Wind—induced load— ing and dynamic responses of a row of tall buildings under [8】PATANKAR S V.Numerical heat transfer and fluid flow 【M].New York Mcgraw2H ill,1980. 【9】黄本才.结构抗风分析原理及应用[M】.上海:同济大学出 版社,2001. stongr interference【J],Journal of Wind Engineeirng and 

因篇幅问题不能全部显示,请点此查看更多更全内容

Top