CN111985019B - 带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法 - Google Patents

带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法 Download PDF

Info

Publication number
CN111985019B
CN111985019B CN202010245719.7A CN202010245719A CN111985019B CN 111985019 B CN111985019 B CN 111985019B CN 202010245719 A CN202010245719 A CN 202010245719A CN 111985019 B CN111985019 B CN 111985019B
Authority
CN
China
Prior art keywords
coefficient
wind
conventional single
tower
load
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010245719.7A
Other languages
English (en)
Other versions
CN111985019A (zh
Inventor
赵爽
晏致涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University of Science and Technology
Original Assignee
Chongqing University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing University of Science and Technology filed Critical Chongqing University of Science and Technology
Priority to CN202010245719.7A priority Critical patent/CN111985019B/zh
Publication of CN111985019A publication Critical patent/CN111985019A/zh
Application granted granted Critical
Publication of CN111985019B publication Critical patent/CN111985019B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/08Construction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Business, Economics & Management (AREA)
  • General Engineering & Computer Science (AREA)
  • Structural Engineering (AREA)
  • General Health & Medical Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Civil Engineering (AREA)
  • Architecture (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • Mathematical Optimization (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Wind Motors (AREA)
  • Buildings Adapted To Withstand Abnormal External Influences (AREA)

Abstract

本发明公开了一种带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法,建立了常规单塔设计风荷载的计算模型,通过计算风振系数进而求取常规单塔设计风荷载。采用统计的影响线方法计算等效背景风荷载,确定了中间变量γB的拟合公式,推导了风振系数考虑外形变化影响的修正系数θηB。通过采用γB、θηB和相关的修正系数,推导了带悬挑横担的常规单塔的风振系数。在实际设计当中,通过风振系数的计算结果、风振系数与风荷载相关的计算公式,即可得到常规单塔的风荷载。计算步骤简单,分析横隔面、横担和钢管混凝土对输电塔风振系数的各贡献量的影响。

Description

带悬挑横担常规单塔基于有效荷载法的风振系数简化计算 方法
技术领域
本发明涉及常规单塔设计技术领域,具体的说是一种带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法。
背景技术
中国规范沿用苏联的抗风设计方法,该方法采用振型分解法计算等效背景风荷载,对于悬臂构筑物而言,以1阶模态的参振为主,振型分解法的近似计算引起的塔顶位移误差不大。与惯性力法对比,由于有效荷载法采用准静态的方法计算背景风荷载,其确定的背景风荷载是准确的理论解。采用有效荷载法计算时,根据响应分量确定背景风荷载和共振风荷载的权重系数,能够达到计算响应与实际响应等效的作用。有效荷载法的物理意义清晰,计算输电塔的内力结果比采用惯性力法得到的结果更准确。澳大利亚和新西兰的杆塔规范、标准基于该方法设计风荷载。然而这些规范、标准中的输电塔设计风荷载都没有考虑横担、梯度风高度和钢管混凝土的影响。准确的输电塔设计风荷载不仅能够保证输电线路的安全运行,还能够降低工程造价,起到经济性作用。
众所周知的,求取风荷载的关键相当于求取单塔的风振系数,则针对现有技术的缺陷,对与常规单塔,有必要将横担、梯度风高度和钢管混凝土三者的影响考虑到计算风振系数的方法中,从而进一步的对单塔的风荷载进行设计,提高现有技术中常规单塔的设计效果。
发明内容
针对上述问题,本发明提供了一种带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法,基于有效荷载法确定单塔的设计风荷载,综合考虑到横担、梯度风高度和钢管混凝土对风荷载的影响。
为达到上述目的,本发明采用的具体技术方案如下:
一种带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法,其关键技术在于具体步骤为:
S1:确定常规单塔的计算参数,并确定常规单塔所在地面粗糙度类别,设定10m高度处的平均分速
Figure GDA0003593319720000021
常规单塔的总高度H;跟开b1;横担个数nc;横担平均外伸长度
Figure GDA0003593319720000022
S2:通过水平均布荷载作用下结构的挠曲线获得步骤S1中的常规单塔0°风向角的1阶侧弯振型φ1(z);
S3:根据荷载规范引入背景分量因子Bz(z),进而计算输电塔的脉动风荷载在水平方向的相关系数ρx和脉动风荷载在竖直方向的相关系数ρz
根据荷载规范引入的共振分量因子R;
根据荷载规范确定地面粗糙度指数α、峰值因子gs、10m高度处的湍流度 I10、瑞流度高度变化系数
Figure GDA0003593319720000023
风压随高度变化系数μz
分别计算常规单塔剩余塔身、横担、横隔面处的风振系数考虑对荷载响应相关系数ρfr′影响的修正系数θηB
S4:根据步骤S1中常规单塔所在地面粗糙度类别,获取背景分量因子的中间变量γ的拟合系数kγ、aγ
根据常规单塔的宽度、深度得到风振系数考虑整体外形变化的修正系数θv
根据常规单塔的横担平均外伸长度、总高度、横担个数得到风振系数考虑附加面积的修正系数θa与风振系数考虑附加质量的修正系数θm的乘积θl
根据常规单塔考虑所在地面粗糙度类别和带悬挑横担的影响因子,确定常规单塔荷载响应相关系数ρfr′的中间变量γB
根据常规单塔考虑带悬挑横担的影响因子,确定常规单塔剩余塔身、横担、横隔面处的新共振分量因子
Figure GDA0003593319720000034
考虑脉动风荷载空间相关性的修正系数θηR
S5:简化并得到常规单塔剩余塔身的风振系数考虑局部外形变化的修正系数的计算公式并计算得到对应的剩余塔身的风振系数考虑局部外形变化的修正系数的θb(z);
计算得到横担的风振系数考虑局部外形变化的修正系数θb(zI)和横隔面的风振系数考虑局部外形变化的修正系数θb(zJ);
S6:考虑外形规律变化的影响因素,当塔高H确定时,计算新背景分量因子
Figure GDA0003593319720000031
考虑横担、横隔面的附加质量和附加面积的影响,计算新共振分量因子
Figure GDA0003593319720000032
S7:计算带悬挑横担常规单塔风振系数β(z)。
进一步的技术方案,步骤S2中常规单塔的脉动风荷载在水平方向的相关系数ρx和脉动风荷载在竖直方向的相关系数ρz的计算公式为:
Figure GDA0003593319720000033
b为外轮廓宽度。
再进一步的技术方案为:步骤S3中:所述共振分量因子R的计算公式为:
Figure GDA0003593319720000041
其中,
Figure GDA0003593319720000042
n为脉动风速的频率;
所述风振系数考虑对荷载响应相关系数ρfr′影响的修正系数θηB的计算公式为:
Figure GDA0003593319720000043
再进一步的技术方案为:步骤S4中荷载响应相关系数ρfr′的中间变量γB的计算公式为:
Figure GDA0003593319720000044
Figure GDA0003593319720000045
其中,荷载响应相关系数的中间变量γB的计算参数取值范围见表1;
表1 常规单塔γB的计算参数取值详情表
Figure GDA0003593319720000046
再进一步的技术方案为:步骤S5所述的常规单塔剩余塔身的风振系数考虑局部外形变化的修正系数的计算公式为:
Figure GDA0003593319720000051
所述横担的风振系数考虑局部外形变化的修正系数θb(zI)计算公式为:
Figure GDA0003593319720000052
所述横隔面的风振系数考虑局部外形变化的修正系数θb(zJ)计算公式为:
Figure GDA0003593319720000053
再进一步的技术方案为:步骤S6中所述的新背景分量因子
Figure GDA0003593319720000054
的计算公式为:
Figure GDA0003593319720000055
所述考虑横担、横隔面的附加质量和附加面积的影响,计算新共振分量因子
Figure GDA0003593319720000056
其中,剩余塔身、横担、横隔面处的新共振分量因子
Figure GDA0003593319720000059
的表达是分别为:
Figure GDA0003593319720000057
再进一步的技术方案为:步骤S7带悬挑横担常规单塔风振系数的计算公式为;
Figure GDA0003593319720000058
本发明的有益效果:建立了常规单塔设计风荷载的计算模型,通过计算风振系数进而求取常规单塔设计风荷载。采用统计的影响线方法计算等效背景风荷载,确定了中间变量γB的拟合公式,推导了风振系数考虑外形变化影响的修正系数θηB。通过采用γB、θηB和相关的修正系数,推导了带悬挑横担的常规单塔的风振系数。在实际设计当中,通过风振系数的计算结果、风振系数与风荷载相关的计算公式,即可得到常规单塔的风荷载。
附图说明
图1是本发明常规单塔结构示意图;
图2是本发明常规单塔
Figure GDA0003593319720000061
沿高度的分布示意图;
图3是常规单塔在不同粗糙度类别下γB(z,H)的拟合值和数值解对比示意图;
图4是常规单塔的计算参数对比示意图;
图5是基于有效荷载法的常规单塔的风振系数分布示意图;
图6是常规塔的时域位移最大值与基于有效荷载法计算的位移最大值对比示意图;
图7是三个工况常规单塔的风振系数贡献量示意图;
图8是本发明方法流程图;
图9为风向角定义图。
具体实施方式
下面结合附图对本发明的具体实施方式以及工作原理作进一步详细说明。
从图1可以看出,常规单塔分为横隔面、横担和剩余塔身3部分,如1所示。图中b1为跟开;b2为横担单肢的外伸长度;H为总高;Hn为呼高;Hg为梯度风高度;H′为混凝土灌溉高度;ΔH1为横担近塔端厚度;ΔH2为相邻横担的中心距;M1(x1,z1)和M2(x2,z2)为空间任意两点。横隔面、横担、剩余塔身的质量和挡风面积沿高度的分布规律不同,在计算过程中需要区别对待。
在本实施例中,常规单塔所在地面粗糙度类别为B类地貌,设定10m高度处的平均分速
Figure GDA0003593319720000062
为速30m/s。整塔阻力系数取2.6。被风面遮挡系数按照DL/T 5154计算,取0.853。通过ANSYS确定基频为1.522Hz。
在本实施例中,针对常规单塔,采用有效荷载法的方式,对风振系数进行求解,进而再求取常规单塔的设计风荷载。
首先,不变外形的常规单塔风振系数的设计公式具体步骤为:
输电塔的等效静力风荷载计算图如图1所示,有效荷载法确定的等效静力风荷载由3部分组成,分别是平均风荷载、等效背景风荷载和等效共振风荷载。
计算风振系数时将阻力系数视为常数,平均风荷载为:
Figure GDA0003593319720000071
等效背景风荷载为:
Figure GDA0003593319720000072
Figure GDA0003593319720000073
式中,gB为背景峰值因子,与荷载规范保持一致,取gB=gs=2.5;ρfr′(z)为荷载响应相关系数;σf(z)为脉动风荷载的标准差。
第1阶振型的等效共振风荷载为:
Figure GDA0003593319720000074
式中,gR为共振峰值因子,同样取gR=gs=2.5;
Figure GDA0003593319720000075
为由共振风荷载激发的1阶模态广义坐标标准差。计算时与荷载规范保持一致,计算
Figure GDA0003593319720000076
时不考虑气动阻尼的影响。
由有效荷载法确定的等效静力风荷载为:
Figure GDA0003593319720000077
WB
Figure GDA0003593319720000078
分别通过响应分量的SRSS法确定。背景响应和共振响应的差别不会太大,为方便计算,美国输电线结构荷载设计指南ASCE No.74近似采用SRSS 方法考虑背景风荷载和共振风荷载的组合。
则上式可以改写成;
Figure GDA0003593319720000079
风振系数表达式为:
Figure GDA00035933197200000710
对于外形沿高度不变的输电塔,
采用Davenport提出的与高度无关的风速谱和Shiotani提出的与频率无关的相干函数,荷载响应相关系数的表达式:
Figure GDA0003593319720000081
式中,Ir(z)为输电塔响应的影响线,
Figure GDA0003593319720000082
Δ为单位荷载作用在塔顶时引起的响应,
Figure GDA0003593319720000083
为Ir(z)沿高度的分布系数。
等效背景风荷载对风振系数的贡献:
Figure GDA0003593319720000084
1阶模态共振响应的广义模态坐标标准差为:
Figure GDA0003593319720000085
该公式仅考虑输电塔1阶频率处的风速谱。
等效共振风荷载对风振系数的贡献:
Figure GDA0003593319720000086
下文重点叙述等效背景风荷载的求解思路。同样地,先从最简单分布的输电塔着手,逐步完善计算模型,进而推导出常规单塔的风振系数,进一步计算设计风荷载。
为了与荷载规范中风振系数计算公式的变量保持一致,在这里引入与荷载规范中物理意义一致的背景分量因子Bz。当塔高H确定时,Bz(z)的表达式为:
Figure GDA0003593319720000091
Figure GDA0003593319720000092
对于外形沿高度不变的输电塔,b(z)和δ(z)均为常数,上式简化为:
ρf′r′(z)=γB(z,H)ρx.
脉动风荷载在水平方向的相关系数ρx采用公式
Figure GDA0003593319720000093
计算,中间变量采用指数函数对其进行拟合,表达式为:
Figure GDA0003593319720000094
其中,中间变量γ的拟合系数kγ、aγ通过表2得到,具体为:
表2 kγ和aγ的取值
地面粗糙度类别 A B C D
k<sub>γ</sub> 1.276 0.910 0.404 0.155
a<sub>γ</sub> 0.186 0.218 0.292 0.376
则,CR的表达式简化为:
Figure GDA0003593319720000095
当塔高H确定时,β(z)表达式用以上变量简化表示为:
Figure GDA0003593319720000101
Figure GDA0003593319720000102
Figure GDA0003593319720000103
为不至于混淆,上式公式中
Figure GDA0003593319720000104
称为新背景分量因子,
Figure GDA0003593319720000105
称为新共振分量因子。
进一步的,常规单塔锥形变化单塔风振系数的设计公式的具体内容为:
以输电塔外形沿高度不变时的计算公式为基础,通过增加风振系数的修正系数来考虑外形规律变化的影响。为了不致混淆起见,同样地对外形沿高度变化的输电塔加以“a”脚码。当塔高H确定时,新背景分量因子为:
Figure GDA0003593319720000106
Figure GDA0003593319720000107
其中,
Figure GDA0003593319720000108
新共振分量因子为:
Figure GDA0003593319720000109
脉动风荷载空间相关性的修正系数θηR与以下θη的表达式一致:
Figure GDA0003593319720000111
剩余塔身的θηR≈1。对于剩余塔身而言,bs(z)沿高度的变化系数
Figure GDA0003593319720000112
和m(z)沿高度的变化系数μm(z)分别按公式
Figure GDA0003593319720000113
Figure GDA0003593319720000114
进行确认。此时的e=1。根据表3可以得到θv,其计算公式为:
Figure GDA0003593319720000115
Figure GDA0003593319720000116
和μm(z)的简化分布确定的θb(z)如下:
Figure GDA0003593319720000117
表3 常规单塔的宽度深度均沿高度作同一规律变化时θv的值
b<sub>s</sub>(H)/b<sub>s</sub>(0) 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1
θ<sub>v</sub> 1.00 1.10 1.20 1.32 1.50 1.75 2.08 2.53 3.30 5.60
则对于常规单塔锥形变化单塔,将公式
Figure GDA0003593319720000118
Figure GDA0003593319720000119
带入
Figure GDA00035933197200001110
即可得到锥形变化的常规单塔的风振系数设计公式。
最后,针对带悬挑横担的常规单塔风振系数的设计公式进行讨论:
为推导等效背景风荷载的设计表达式,需要弄清输电塔响应影响线的分布关系。带悬挑横担的常规单塔Iu(z)的表达式为:
Figure GDA0003593319720000121
式中,
Figure GDA0003593319720000122
为Iu(z)沿高度的分布系数。
在对高度进行积分计算时,Δ与积分变量无关,可以作为常数放在积分表达式外。并且,在计算等效背景风荷载时,分子与分母的Δ相互约去。如此,不用考虑Δ。
结合图1可以看出,带悬挑横担的常规单塔
Figure GDA0003593319720000123
的简化分布与实际分布的对比曲线详见图2;对比结果表示两者的吻合度好。
Figure GDA0003593319720000124
的简化分布公式
Figure GDA0003593319720000125
代入公式:
Figure GDA0003593319720000126
由于γB(z,H)的表达式涉及复杂的2重积分运算,为方便设计使用,采用非线性曲面拟合γB(z,H)。
Figure GDA0003593319720000127
Figure GDA0003593319720000128
以上参数在不同地面粗糙度类别下的取值如表1所示。
表1 常规单塔γB的计算参数取值表
Figure GDA0003593319720000129
Figure GDA0003593319720000131
不同粗糙度类别下,对比分别采用拟合公式
Figure GDA0003593319720000132
和数值计算求解γB(z,H)的结果,对比示意图如图3所示,对比结果表明两者的吻合度好。
风振系数考虑外形对ρfr′影响的修正系数θηB的表达式复杂、影响因素多,通过对输电塔的θηB进行计算分析,结果表明:θηB主要受轮廓宽度变化的影响,在剩余塔身和横隔面处,θηB略小于1。根据公式
Figure GDA0003593319720000133
θηB的取值对β的影响要小于对
Figure GDA0003593319720000134
的影响,从而可近似取θηB=1;在横担位置,轮廓宽度突然变大,θηB比剩余塔身和横隔面处的取值更小,需要计算考虑。θηB近似表示为:
Figure GDA0003593319720000135
对比γB(z,H)、θηB的实际分布与由公式
Figure GDA0003593319720000136
Figure GDA0003593319720000137
确定的简化分布,具体详见图4。两者对比的吻合度好,表明采用γB(z,H)、θηB的简化分布确定带悬挑横担的常规单塔的风荷载设计公式是可靠的。
新共振分量因子为:
Figure GDA0003593319720000141
对于带悬挑横担的常规单塔,横担两端的尺寸显著大于剩余塔身的外轮廓宽度,需要计入横担对脉动风荷载空间相关性的影响,此时θηR不为1。此时θηR的取值见表4所示:
表4 常规单塔的θηR取值范围表
Figure GDA0003593319720000142
修正系数θl为考虑横担、横隔面的附加质量和附加面积的影响,且θl=θaθm。其中,θa风振系数考虑附加面积的修正系数,θm风振系数考虑附加质量的修正系数。
修正系数θl的取值详见表5,根据常规单塔的横担平均外伸长度、总高度、横担个数得到。
表5 常规单塔的θl取值
Figure GDA0003593319720000143
Figure GDA0003593319720000151
其中,带悬挑横担的常规输电塔的横隔面在呼高以上和呼高以下按两种方式布置。呼高以上,横担近塔端的上下边缘各有1个横隔面。呼高以下,塔身坡度变化断面有1横隔面,假定该横隔面与上端相邻横隔面的间距为
Figure GDA0003593319720000152
塔腿处有1横隔面,距离地面高度为跟开b1;在上下两横隔面之间有nd个横隔面,
Figure GDA0003593319720000153
round表示四舍五入取整数,每个横隔面按等间距布置。以上确定了横隔面的在塔身的分布关系,还需要确定横隔面和横担的挡风面积、质量沿高度的分布关系。同样地,对多个带悬挑横担的常规输电塔进行统计分析,获得关系:
Figure GDA0003593319720000154
M(zJ)=1.334m(zJ)。
横担挡风面积沿高度的变化系数:
Figure GDA0003593319720000155
横隔面挡风面积沿高度的变化系数:
Figure GDA0003593319720000156
横担质量沿高度的变化系数:
Figure GDA0003593319720000157
横隔面质量沿高度的变化系数:
Figure GDA0003593319720000158
将θa和θm计算修改参数的独立影响参数均为
Figure GDA0003593319720000159
H、nc、bs(H)/bs(0)。为了制表方便,令θl=θaθm
其中,
Figure GDA00035933197200001510
Figure GDA00035933197200001511
尽管实际分布与简化分布之间有差异,但是通过公式θa和θm计算θl时,由于
Figure GDA0003593319720000161
在分子,μM(zI)、μM(zJ)在分母,误差部分抵消掉。在一定范围内,采用简化分布计算θl时引起的误差可以忽略,可采用简化分布进行近似计算,该结论在计算中得到了证明。
局部修正系数θb(z)、θb(zI)、θb(zJ)的计算公式为:
剩余塔身的风振系数考虑局部外形变化的修正系数的计算公式为:
Figure GDA0003593319720000162
计算得到横担的风振系数考虑局部外形变化的修正系数θb(zI)以及和横隔面的风振系数考虑局部外形变化的修正系数θb(zJ);
横担的风振系数考虑局部外形变化的修正系数θb(zI)计算公式为:
Figure GDA0003593319720000163
横隔面的风振系数考虑局部外形变化的修正系数θb(zJ)计算公式为:
Figure GDA0003593319720000164
为保证书写的简洁,将剩余塔身、横隔面和横担的新共振分量因子合并书写,三者z高度处的
Figure GDA0003593319720000165
可按下式计算:
Figure GDA0003593319720000166
将γB(z,H)和θηB代入公式
Figure GDA0003593319720000167
计算带悬挑横担的常规单塔的
Figure GDA0003593319720000168
Figure GDA0003593319720000169
和上式的
Figure GDA00035933197200001610
代入公式
Figure GDA00035933197200001611
即可得到带悬挑横担的常规单塔的风振系数设计公式。
综上,在本实施例中,一种带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法,结合图8可以看出,具体步骤为:
S1:确定常规单塔的计算参数,并确定常规单塔所在地面粗糙度类别,设定10m高度处的平均分速
Figure GDA00035933197200001612
常规单塔的总高度H;跟开b1;横担个数nc;横担平均外伸长度
Figure GDA0003593319720000171
S2:通过水平均布荷载作用下结构的挠曲线获得步骤S1中的常规单塔0°风向角的1阶侧弯振型φ1(z);具体公式为:
Figure GDA0003593319720000172
S3:根据荷载规范引入背景分量因子Bz(z),进而计算输电塔的脉动风荷载在水平方向的相关系数ρx和脉动风荷载在竖直方向的相关系数ρz
步骤S3常规单塔的脉动风荷载在水平方向的相关系数ρx和脉动风荷载在竖直方向的相关系数ρz的计算公式为:
Figure GDA0003593319720000173
b为外轮廓宽度。
根据荷载规范引入的共振分量因子R;
所述共振分量因子R的计算公式为:
Figure GDA0003593319720000174
其中,
Figure GDA0003593319720000175
n为脉动风速的频率;
根据荷载规范确定地面粗糙度指数α、峰值因子gs、10m高度处的湍流度 I10、瑞流度高度变化系数
Figure GDA0003593319720000176
风压随高度变化系数μz
分别计算常规单塔剩余塔身、横担、横隔面处的风振系数考虑对荷载响应相关系数ρfr′影响的修正系数θηB
风振系数考虑对荷载响应相关系数ρfr′影响的修正系数θηB的计算公式为:
Figure GDA0003593319720000177
S4:根据步骤S1中常规单塔所在地面粗糙度类别,获取背景分量因子的中间变量γ的拟合系数kγ、aγ
根据常规单塔的宽度、深度得到风振系数考虑整体外形变化的修正系数θv
根据常规单塔的横担平均外伸长度、总高度、横担个数得到风振系数考虑附加面积的修正系数θa与风振系数考虑附加质量的修正系数θm的乘积θl
根据常规单塔考虑所在地面粗糙度类别和带悬挑横担的影响因子,确定常规单塔荷载响应相关系数ρfr′的中间变量γB
步骤S4中荷载响应相关系数ρfr′的中间变量γB的计算公式为:
Figure GDA0003593319720000181
Figure GDA0003593319720000182
其中,荷载响应相关系数的中间变量γB的计算参数取值范围见表1;
根据常规单塔考虑带悬挑横担的影响因子,确定常规单塔剩余塔身、横担、横隔面处的新共振分量因子
Figure GDA0003593319720000183
考虑脉动风荷载空间相关性的修正系数θηR
S5:简化并得到常规单塔剩余塔身的风振系数考虑局部外形变化的修正系数的计算公式并计算得到对应的剩余塔身的风振系数考虑局部外形变化的修正系数的θb(z);
计算得到横担的风振系数考虑局部外形变化的修正系数θb(zI)和横隔面的风振系数考虑局部外形变化的修正系数θb(zJ);
步骤S5所述的常规单塔剩余塔身的风振系数考虑局部外形变化的修正系数的计算公式为:
Figure GDA0003593319720000191
所述横担的风振系数考虑局部外形变化的修正系数θb(zI)计算公式为:
Figure GDA0003593319720000192
所述横隔面的风振系数考虑局部外形变化的修正系数θb(zJ)计算公式为:
Figure GDA0003593319720000193
S6:考虑外形规律变化的影响因素,计算当塔高H确定时,新背景分量因子
Figure GDA0003593319720000194
步骤S6中所述的新背景分量因子
Figure GDA0003593319720000195
的计算公式为:
Figure GDA0003593319720000196
所述考虑横担、横隔面的附加质量和附加面积的影响,计算新共振分量因子
Figure GDA0003593319720000197
其中,剩余塔身、横担、横隔面处的新共振分量因子
Figure GDA0003593319720000198
的表达是分别为:
Figure GDA0003593319720000199
考虑横担、横隔面的附加质量和附加面积的影响,计算新共振分量因子
Figure GDA00035933197200001910
S7:计算带悬挑横担常规单塔风振系数β(z)。
步骤S7带悬挑横担常规单塔风振系数的计算公式为;
Figure GDA00035933197200001911
在本实施例中,通过公式
Figure GDA00035933197200001912
计算x1=60.871。查表1确定γB(z,H)的计算参数,通过公式
Figure GDA00035933197200001913
计算γB(z,H)。通过公式
Figure GDA0003593319720000201
计算θηB。将以上参数代入公式
Figure GDA0003593319720000202
计算得到
Figure GDA0003593319720000203
按钢结构取ζ1=0.01,通过公式
Figure GDA0003593319720000204
计算R=1.839。查表2取值kγ=0.91、aγ=0.218,塔高H=91.1m通过公式
Figure GDA0003593319720000205
计算γR(H)=2.433。跟开b1=17.28m,通过公式
Figure GDA0003593319720000206
计算ρz=0.730、ρx=0.945。塔底主材外径0.61m,塔底斜材外径0.219m,由两个主材和两个斜材构成的总迎风外径1.658m。塔顶主材外径 0.273m,塔顶斜材外径0.159m,由两个主材和两个斜材构成的总迎风外径 0.864m。由主材与斜材确定的bs(H)/bs(0)=0.521,查表3确定θv=1.697。横担个数nc=3,横担平均外伸长度14.822m,
Figure GDA0003593319720000207
查表4确定θηR=0.952,查表5确定θl=0.890。
分别按公式
Figure GDA0003593319720000208
计算θb(z)。然后将以上参数代入公式
Figure GDA00035933197200002012
计算
Figure GDA00035933197200002010
再代入公式
Figure GDA00035933197200002011
计算风振系数。β(z)的分布如图5所示,剩余塔身、横隔面和横担的风振系数随高度的增大而增大,同高度处三者的风振系数大体上为横隔面的最大,剩余塔身的次之,横担的最小。由于完整考虑背景风荷载,有效荷载风振系数要比惯性力风振系数大,随着高度的增大,差异逐渐减小,塔顶处的风振系数接近。相比较而言,两种方法计算剩余塔身风振系数时的差异最明显,计算横隔面时的差异最小。
对建立好的输电塔有限元模型进行模拟风场下的风振时域计算,来流风为 0°风向角。获取输电塔不同高度处节点的顺风向最大位移,并与本章推导的设计风荷载作用下输电塔的节点位移进行比较。
在本发明中,结合图9可以定义,当来流风平行于横担轴向时风向角θ=0°,当来流风平行于导线走向时风向角θ=90°。其中,x向表示横担轴向,y向表示顺线向。
风振系数随风向角的变化不大,并且风向角对塔身风振系数和横担风振系数的影响是相反的,对整塔而言该影响可以抵消。电力相关标准中仅考虑0o风向角下输电塔的风振系数。因此,可以忽略风向角对风振系数的影响,其它风向角下的等效静力风荷载通过DL/T 5154中的风荷载分配系数确定。
其中,电力相关标准包括:GB 50545-2010.110kV~750kV架空输电线路设计规范[S].北京:中国计划出版社,2010;GB 50665-2011.1000kV架空输电线路设计规范[S].北京:中国计划出版社,2011;DL/T 5154-2012.架空输电线路杆塔结构设计技术规定[S].北京:中国计划出版社,2012;DL/T 5504-2015.特高压架空输电线路大跨越设计技术规定[S].北京:中国计划出版社,2015。
如图6所示。图6(a)中,输电塔的风振位移由下到上非线性增大,时域和设计风荷载的计算结果吻合度好。图6(b)中,以时域结果作为准确值,由设计风荷载计算风振位移的误差由下到上逐渐增大,底部误差的绝对值最大。输电塔的底部位移小,尽管底部的相对误差大,但是绝对误差并不大,并且输电塔的抗风设计由顶部位移的大小控制。因此,对于带悬挑横担的常规单塔,由本章推导的设计风荷载计算的结果能够满足工程需要。常规单塔的影响参数对3种工况下图1的带悬挑横担的常规单塔进行计算分析。工况1为只考虑剩余塔身;工况2在工况1的基础上增加横隔面;工况3在工况2的基础上增加横担。3种工况的输电塔风振系数贡献量如图7所示。图7(a)中,增加横隔面、横担后对CB(z) 无影响。图7(b)中,增加横隔面后剩余塔身的CR(z)略微减小,增加横担后剩余塔身和横隔面的CR(z)均显著减小。横隔面对输电塔剩余塔身的CR(z)的影响小,横担对输电塔剩余塔身和横隔面的CR(z)的影响大。
应当指出的是,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的普通技术人员在本发明的实质范围内所做出的变化、改性、添加或替换,也应属于本发明的保护范围。

Claims (4)

1.一种带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法,其特征在于具体步骤为:
S1:确定常规单塔的计算参数,并确定常规单塔所在地面粗糙度类别,设定10m高度处的平均分速
Figure FDA0003557759870000011
常规单塔的总高度H;跟开b1;横担个数nc;横担平均外伸长度
Figure FDA0003557759870000012
S2:通过水平均布荷载作用下结构的挠曲线获得步骤S1中的常规单塔0°风向角的1阶侧弯振型φ1(z),具体公式为:
Figure FDA0003557759870000013
S3:根据荷载规范引入背景分量因子Bz(z),进而计算输电塔的脉动风荷载在水平方向的相关系数ρx和脉动风荷载在竖直方向的相关系数ρz
根据荷载规范引入的共振分量因子R;
根据荷载规范确定地面粗糙度指数α、峰值因子gs、10m高度处的湍流度I10、瑞流度高度变化系数
Figure FDA0003557759870000014
风压随高度变化系数μz
分别计算常规单塔剩余塔身、横担、横隔面处的风振系数考虑对荷载响应相关系数ρfr'影响的修正系数θηB
S4:根据步骤S1中常规单塔所在地面粗糙度类别,获取背景分量因子的中间变量γ的拟合系数kγ、aγ
根据常规单塔的宽度、深度得到风振系数考虑整体外形变化的修正系数θv
根据常规单塔的横担平均外伸长度、总高度、横担个数得到风振系数考虑附加面积的修正系数θa与风振系数考虑附加质量的修正系数θm的乘积θl
根据常规单塔考虑所在地面粗糙度类别和带悬挑横担的影响因子,确定常规单塔荷载响应相关系数ρfr'的中间变量γB
根据常规单塔考虑带悬挑横担的影响因子,确定常规单塔剩余塔身、横担、横隔面处的新共振分量因子
Figure FDA0003557759870000021
考虑脉动风荷载空间相关性的修正系数θηR
S5:简化并得到常规单塔剩余塔身的风振系数考虑局部外形变化的修正系数的计算公式并计算得到对应的剩余塔身的风振系数考虑局部外形变化的修正系数的θb(z);
计算得到横担的风振系数考虑局部外形变化的修正系数θb(zI)和横隔面的风振系数考虑局部外形变化的修正系数θb(zJ);
S6:考虑外形规律变化的影响因素,计算当塔高H确定时,新背景分量因子
Figure FDA0003557759870000022
考虑横担、横隔面的附加质量和附加面积的影响,计算新共振分量因子
Figure FDA0003557759870000023
S7:计算带悬挑横担常规单塔风振系数β(z);
步骤S3中常规单塔的脉动风荷载在水平方向的相关系数ρx和脉动风荷载在竖直方向的相关系数ρz的计算公式为:
Figure FDA0003557759870000024
b为外轮廓宽度,中间变量采用指数函数对其进行拟合,表达式为:
Figure FDA0003557759870000025
步骤S3中:所述共振分量因子R的计算公式为:
Figure FDA0003557759870000026
其中,
Figure FDA0003557759870000027
Figure FDA0003557759870000028
n为脉动风速的频率;
所述风振系数考虑对荷载响应相关系数ρfr'影响的修正系数θηB的计算公式为:
Figure FDA0003557759870000031
步骤S4中荷载响应相关系数ρfr'的中间变量γB的计算公式为:
Figure FDA0003557759870000032
Figure FDA0003557759870000033
其中,荷载响应相关系数的中间变量γB的计算参数取值范围见表1;
表1 常规单塔γB的计算参数取值详情表
Figure FDA0003557759870000034
2.根据权利要求1所述的带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法,其特征在于步骤S5所述的常规单塔剩余塔身的风振系数考虑局部外形变化的修正系数的计算公式为:
Figure FDA0003557759870000041
所述横担的风振系数考虑局部外形变化的修正系数θb(zI)计算公式为:
Figure FDA0003557759870000042
所述横隔面的风振系数考虑局部外形变化的修正系数θb(zJ)计算公式为:
Figure FDA0003557759870000043
3.根据权利要求1所述的带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法,其特征在于步骤S6中所述的新背景分量因子
Figure FDA0003557759870000044
的计算公式为:
Figure FDA0003557759870000045
所述考虑横担、横隔面的附加质量和附加面积的影响,计算新共振分量因子
Figure FDA0003557759870000046
其中,剩余塔身、横担、横隔面处的新共振分量因子
Figure FDA0003557759870000047
的表达是分别为:
Figure FDA0003557759870000048
4.根据权利要求1所述的带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法,其特征在于步骤S7带悬挑横担常规单塔风振系数的计算公式为;
Figure FDA0003557759870000049
CN202010245719.7A 2020-03-31 2020-03-31 带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法 Active CN111985019B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010245719.7A CN111985019B (zh) 2020-03-31 2020-03-31 带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010245719.7A CN111985019B (zh) 2020-03-31 2020-03-31 带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法

Publications (2)

Publication Number Publication Date
CN111985019A CN111985019A (zh) 2020-11-24
CN111985019B true CN111985019B (zh) 2022-05-20

Family

ID=73442195

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010245719.7A Active CN111985019B (zh) 2020-03-31 2020-03-31 带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法

Country Status (1)

Country Link
CN (1) CN111985019B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116150842B (zh) * 2022-12-30 2023-09-22 重庆科技学院 基于iwl法计算弯扭柱螺旋梁景观塔设计风荷载的方法
CN117848280B (zh) * 2024-03-07 2024-05-24 中建五局第三建设有限公司 一种异型大跨度悬挑钢结构的在线健康监测方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101692566A (zh) * 2009-09-09 2010-04-07 中国电力科学研究院 一种输电线路高塔风振控制方法
CN106227951A (zh) * 2016-07-26 2016-12-14 南京航空航天大学 适用于超规范特大型冷却塔的二维风振系数计算方法
CN107291985A (zh) * 2017-05-23 2017-10-24 南京航空航天大学 一种冷却塔施工全过程风振系数取值方法
CN108959742A (zh) * 2018-06-20 2018-12-07 重庆科技学院 大跨越输电塔线体系气动弹性模型设计方法
CN110378050A (zh) * 2019-07-25 2019-10-25 重庆科技学院 一种输电塔顺风向气动阻尼比计算方法
CN110851905A (zh) * 2019-11-14 2020-02-28 福州大学 一种峡谷微地形条件下的输电铁塔风振系数计算方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102609579B (zh) * 2012-01-31 2014-04-23 浙江省电力设计院 一种输电塔三维抗风设计方法
US9644606B2 (en) * 2012-06-29 2017-05-09 General Electric Company Systems and methods to reduce tower oscillations in a wind turbine
CN109101726B (zh) * 2018-08-10 2022-09-09 国网浙江省电力有限公司经济技术研究院 一种基于风载总计的输电塔整体脉动风荷载谱确定方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101692566A (zh) * 2009-09-09 2010-04-07 中国电力科学研究院 一种输电线路高塔风振控制方法
CN106227951A (zh) * 2016-07-26 2016-12-14 南京航空航天大学 适用于超规范特大型冷却塔的二维风振系数计算方法
CN107291985A (zh) * 2017-05-23 2017-10-24 南京航空航天大学 一种冷却塔施工全过程风振系数取值方法
CN108959742A (zh) * 2018-06-20 2018-12-07 重庆科技学院 大跨越输电塔线体系气动弹性模型设计方法
CN110378050A (zh) * 2019-07-25 2019-10-25 重庆科技学院 一种输电塔顺风向气动阻尼比计算方法
CN110851905A (zh) * 2019-11-14 2020-02-28 福州大学 一种峡谷微地形条件下的输电铁塔风振系数计算方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
"1000kV苏通大跨越输电塔线体系气弹模型的风洞试验研究";赵爽 等;《中国电机工程学》;20180905;第38卷(第17期);第5257-5264页 *
"Dynamic Analysis of Wind Power Turbine"s Tower under the Combined Action of Winds and Waves";CHENG You-liang等;《International Journal of Plant Engineering and Management》;20170915;第22卷(第03期);第140-149页 *
"Research on Wind-induced Vibration Hybrid Control of Hefei TV Tower";Zhiqiang Zhang 等;《2011 International Conference on Electric Technology andCivil Engineering (ICETCE)》;20110527;第5387-5391页 *
"双呼高输电塔的风振响应及风振系数研究";刘俊卿 等;《西安建筑科技大学学报(自然科学版)》;20200228;第52卷(第01期);第9-14页 *
"基于不同方法的超高杆塔结构风振系数计算对比";游溢 等;《科技创新与应用》;20170518(第2017年第14期);第41-42页 *
"基于风洞试验的苏通大跨越输电塔风振系数研究";赵爽 等;《建筑结构学报》;20190807;第40卷(第11期);第35-44页 *
"输电杆塔风荷载风振系数的数值计算与规范比较";汪大海 等;《建筑结构》;20180710;第48卷(第13期);第49-56页 *

Also Published As

Publication number Publication date
CN111985019A (zh) 2020-11-24

Similar Documents

Publication Publication Date Title
CN111985019B (zh) 带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法
CN111506946B (zh) 基于有效荷载法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法
CN111651806B (zh) 基于惯性力法的常规输电塔设计风载荷计算方法
CN111651805B (zh) 考虑线形与线长影响因子的悬垂绝缘子串最大风偏角及风振系数计算方法
CN112395752B (zh) 基于螺栓节点刚度的输电铁塔位移计算系统及方法
CN109614659B (zh) 一种悬垂绝缘子串风偏荷载调整系数确定方法
CN111651804B (zh) 考虑塔线耦合影响的常规输电线路设计风载荷计算方法
CN109800442B (zh) 航空发动机叶片模型重构方法
CN111985018A (zh) 基于惯性力法和塔线分离法考虑塔线耦合影响的超高大跨越塔、线设计风载荷的计算方法
CN111651807B (zh) 基于有效荷载法的超高单塔风振系数的简化计算方法
CN112949217A (zh) 风电机组尾流分区域建模及快速分析方法
CN111651808B (zh) 考虑梯度风和混凝土质量影响的超高输电塔风振系数简化计算方法
CN106484927B (zh) 一种导线舞动失稳分析方法
CN111143941A (zh) 一种复合材料加筋壁板轴压承载能力计算方法
CN108038292A (zh) 一种基于双代理技术的高效自适应采样方法
CN111506944B (zh) 考虑塔线耦合效应的杆塔风荷载脉动折减系数计算方法
CN108894120B (zh) 拱上建筑拆装过程中拱圈变形安全判断方法
CN116226986A (zh) 一种拉索无应力长度计算方法
CN109522648B (zh) 一种考虑运动气动力的尾流下弹性支撑圆柱驰振分析方法
CN111753435A (zh) 基于分段悬链线和索力连续算法的缆索吊装系统计算方法
CN111506945B (zh) 基于塔线耦合影响因子的输电塔共振响应的等效阻尼系数计算方法
CN113686287B (zh) 基于悬垂杆塔绝缘子串偏斜角的覆冰厚度测量方法及系统
CN110968935B (zh) 一种任意高差覆冰输电线路脱冰最大跳跃高度的算法
CN109190931B (zh) 一种基于三维空间坐标的硬跳线迭代计算的方法
CN110990962B (zh) 辅助变流柜用风机的智能优化方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant