CN111651920B - 塔线体系下大跨越输电塔风致机械能的传递分析方法 - Google Patents
塔线体系下大跨越输电塔风致机械能的传递分析方法 Download PDFInfo
- Publication number
- CN111651920B CN111651920B CN202010488810.1A CN202010488810A CN111651920B CN 111651920 B CN111651920 B CN 111651920B CN 202010488810 A CN202010488810 A CN 202010488810A CN 111651920 B CN111651920 B CN 111651920B
- Authority
- CN
- China
- Prior art keywords
- wind
- tower
- power transmission
- average
- mechanical energy
- 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
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 256
- 238000004458 analytical method Methods 0.000 title description 4
- 238000012360 testing method Methods 0.000 claims abstract description 37
- 238000004364 calculation method Methods 0.000 claims abstract description 26
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000012546 transfer Methods 0.000 claims abstract description 19
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 8
- 238000005516 engineering process Methods 0.000 claims abstract description 6
- 238000006073 displacement reaction Methods 0.000 claims description 76
- 230000004044 response Effects 0.000 claims description 28
- 239000011159 matrix material Substances 0.000 claims description 15
- 230000009467 reduction Effects 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 4
- 239000004020 conductor Substances 0.000 claims description 3
- 230000010355 oscillation Effects 0.000 claims description 3
- 238000011160 research Methods 0.000 abstract description 19
- 238000013016 damping Methods 0.000 description 18
- 230000009471 action Effects 0.000 description 14
- 230000000694 effects Effects 0.000 description 14
- 238000013461 design Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 230000007246 mechanism Effects 0.000 description 8
- 238000005265 energy consumption Methods 0.000 description 5
- 238000005381 potential energy Methods 0.000 description 5
- 229910000831 Steel Inorganic materials 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 239000010959 steel Substances 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- RWSOTUBLDIXVET-UHFFFAOYSA-N Dihydrogen sulfide Chemical compound S RWSOTUBLDIXVET-UHFFFAOYSA-N 0.000 description 2
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000004870 electrical engineering Methods 0.000 description 2
- 229910000037 hydrogen sulfide Inorganic materials 0.000 description 2
- 239000012212 insulator Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000009022 nonlinear effect Effects 0.000 description 2
- 230000010349 pulsation Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 239000004411 aluminium Substances 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 229910001220 stainless steel Inorganic materials 0.000 description 1
- 239000010935 stainless steel Substances 0.000 description 1
- 239000005437 stratosphere Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 125000000341 threoninyl group Chemical group [H]OC([H])(C([H])([H])[H])C([H])(N([H])[H])C(*)=O 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/06—Wind turbines or wind farms
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种塔线体系下大跨越输电塔风致机械能的传递分析方法,具体步骤为:基于大跨越输电塔搭建塔线体系,并获取塔线体系物理参数;构建刚性模型,进行风洞试验,测得风阻力参数;构建大跨越气动弹性模型,进行风洞试验,测得动力参数;使用正交分解技术确定其参振模态;设定工况,并得出每一种工况下输电塔风致机械能最大值的计算模型;计算与工况对应的输电塔风致机械能最大值;分析影响因子输电线、风速和梯度风高度对风致机械能分量的影响以及子结构之间能量的传递关系。有益效果:采用大跨越模型确定的阻力系数和动力参数,通过机械能指标对其风振特性进行拓展性研究。
Description
技术领域
本发明涉及大跨越输电体系技术领域,具体的说是一种塔线体系下大跨越输电塔风致机械能的传递分析方法。
背景技术
由于我国能源分布的极不平衡,实现大规模、高效率、远距离输电的特高压交直流输电方案将作为我国经济社会健康可持续发展要求的国家电网发展模式。拟建的大跨越输电塔线路(下文简称大跨越)高度高、跨度大,受风荷载的影响显著。
由于刚性模型的设计和制作相对容易,早期开展了输电塔的气动力研究。通过研究密实度、风向角、风攻角等参数对阻力系数的影响,提出了阻力系数的计算公式。然而受模型比例和边界条件的影响,风洞试验获取的阻力系数通常要比实测结果大得多。输电塔线体系表现出复杂的风致耦合振动现象,对风振响应的研究引起了多数学者的关注。
初期,通过采用小比例的设计方案制作了输电塔线体系的气动弹性模型。然而这样会使得导线单位长度的质量和刚度过小,导致模型在结构基本缩尺律和构件层面的缩尺律等方面都存在一些问题。现场测试是理论研究的重要依据,通过在多山地区对全尺寸塔线体系进行了实测,结果表明导线影响输电塔的振动特性,导线的气动阻尼显著。然而,现场实测研究的风场特性很难全面掌握,同时也不利于参数化的研究分析。
文献Loredo-Souza A M,Davenport A G.A novel approach for wind tunnelmodelling of transmission lines[J].Journal of Wind Engineering and IndustrialAerodynamics,2001,89(11-12):1017-1029.提出一种变比例的输电线模型设计方案,从而有效地解决了塔线体系气动弹性模型设计的难题。在此基础上,诸多学者进行了特高压输电塔线路气动弹性模型的风洞试验研究。通过对1 000kV特高压输电线路进行数值模拟和风洞试验研究,揭示了塔线耦联体系的风荷载传递机制。
此外,采用设计好的塔线体系气动弹性模型,研究500kV高压输电线路的塔线耦联作用和1 000kV特高压输电线路的风致响应特性。在特高压拉线塔线体系方面,通过气动弹性模型的风洞试验研究表明拉线塔的风振机制有别于自立式杆塔的。
上述研究为特高压输电塔线体系的抗风提供了重要参考。与传统的特高压输电塔线路相比,大跨越结构复杂、需要考虑的影响因素多,使得其风振特性有别于以往的研究。例如:苏通大跨越。同时,风洞试验研究不能够对塔身风荷载和输电线风荷载进行分解探讨。
发明内容
针对上述问题,本发明提供了一种塔线体系下大跨越输电塔风致机械能的传递分析方法,采用大跨越模型(气动弹性模型和刚性模型)确定的阻力系数和动力参数,通过机械能指标对其风振特性进行拓展性研究。
为达到上述目的,本发明采用的具体技术方案如下:
一种塔线体系下大跨越输电塔风致机械能的传递分析方法,具体步骤为:
S1:基于大跨越输电塔搭建塔线体系,并获取塔线体系物理参数;
S2:基于塔线体系,设定模型缩小比例,构建刚性模型,基于刚性模型进行风洞试验,测得风阻力参数;
基于塔线体系,设定模型缩小比例,构建大跨越气动弹性模型,基于大跨越气动弹性模型进行风洞试验,测得动力参数;
使用正交分解技术(POD)确定其参振模态;
S3:以输电线大变形的平均风状态作为计算初始条件,基于随机振动理论和线性叠加原理,设定工况,并得出每一种工况下输电塔风致机械能最大值的计算模型;
S4:根据步骤S3中得到的输电塔风致机械能最大值的计算模型,根据工况,计算与工况对应的输电塔风致机械能最大值;
S5:分析影响因子输电线、风速和梯度风高度对风致机械能分量的影响以及子结构之间能量的传递关系。
通过上述设计,依据随机振动理论提出输电塔风致机械能最大值的计算模型。通过数值计算研究输电线、风速、梯度风高度对输电塔风致机械能的影响,以及塔身风荷载和输电线风荷载激发的输电塔机械能成分,进而揭示出超高大跨越输电塔风致机械能的传递机制。通过本发明技术方案,能够让抗风设计人员对此类结构的风振特性有更为深入的理解。
进一步的技术方案为:所述塔线体系物理参数至少包括输电塔布置方案、导线地线的跨度垂度参数、导线型号参数、地线型号参数。
在本发明中,还可能用到其他参数,具体详见文献:赵爽,晏致涛,李正良等.1000kV苏通大跨越输电塔线体系气弹模型的风洞试验研究[J].中国电机工程学报,2018,38(17):5257-5265;相关的行业专有名词或概念见文献:张殿生.电力工程高压送电线路设计手册(第2版)[M].中国:中国电力出版社,2003.
再进一步的技术方案,步骤S2中,测量所述风阻力参数时:刚性模型采用五分量应变天平测量气动力,并设定采样时间、采样频率、来流风向、试验风速转化为实际基本风速比例;其中,刚性模型按照设定的比例进行设计,并且确定输电塔构件的材质。
式中,α为地面粗糙度指数。
再进一步的技术方案,在步骤S2中,在进行动力参数测量时,采用气动弹性模型进行风洞试验研究,并且设计到两个问题:1)钢管混凝土材料的非线性影响;2)梯度风高度对风振响应的影响;对于上述两个问题以及具体的测量步骤,在文献“赵爽,晏致涛,李正良等.1000kV苏通大跨越输电塔线体系气弹模型的风洞试验研究[J].中国电机工程学报,2018,38(17):5257-5265;”中有详细提及。
使用正交分解技术(POD)确定其参振模态;具体方法详见文献“Yang J N,LeiYing,Pan Shuwen,et al.System identification of linear structures based onhilber-huang spectral analysis.Part I:normal modes[J].Earthquake Engineeringand Structural Dynamics,2003,32:1443-1467.”和文献“Ibrahim S R.Randomdecrement technique for modal identification of structures[J].Journal ofSpacecraft and Rockets,1997,14(11):696-698.”
再进一步的技术方案,所述输电塔风致机械能最大值的计算模型为:
其中,塔体的平均风荷载和输电线传递的平均风荷载在各自对应的平均位移下做的实功总和为:
塔体的平均风荷载在对应的平均位移下做的实功为:
输电线传递的平均风荷在对应的平均位移下做的实功为:
塔体的平均风荷载和输电线传递的平均风荷载在峰值背景位移下做的虚功总和为:
塔线体系下,输电塔的峰值背景位移包括输电塔峰值背景风荷载和输电线峰值背景风荷载两部分激发的输电塔位移。塔体的平均风荷载在峰值背景位移下做的虚功为:
输电线传递的平均风荷载在峰值背景位移下做的虚功为:
塔体的峰值背景风荷载和输电线传递的峰值背景风荷载在各自对应的峰值背景位移下做的实功总和为:
塔体的峰值背景风荷载在对应的峰值背景位移下做的实功总和为:
输电线传递的峰值背景风荷载在对应的峰值背景位移下做的实功总和为:
塔体的峰值背景风荷载和输电线传递的峰值背景风荷载在平均位移下做的虚功总和为:
塔线体系下,输电塔的平均位移包括输电塔平均风荷载和输电线平均风荷载两部分激发的输电塔位移。塔体的峰值背景风荷载在平均位移下做的虚功为:
输电线传递的峰值背景风荷载在平均位移下做的虚功为:
塔体的共振响应激发风致动能的最大值为:
式中,下标“t2”表示挂线时的输电塔;W1为塔体荷载做功;W2为传递荷载做功;和Rb,max分别为传递荷载的平均分量和峰值背景分量;IR为KIR=E的解,E为单位矩阵;分别通过对输电线水平档距内的平均风荷载和峰值背景风荷载求和计算和Rb,max;
本发明的有益效果:依据随机振动理论提出输电塔风致机械能最大值的计算模型。通过数值计算研究输电线、风速、梯度风高度对输电塔风致机械能的影响,以及塔身风荷载和输电线风荷载激发的输电塔机械能成分,进而揭示出超高大跨越输电塔风致机械能的传递机制。通过本文的研究,能够让抗风设计人员对此类结构的风振特性有更为深入的理解。
附图说明
图1是塔线体系中输电塔结构示意图;
图2是苏通大跨越的塔线体系结构示意图;
图3是刚性模型现场搭建示意图;
图4是刚性模型的风洞试验确定输电塔各段的阻力系数分布示意图;
图5是在各风速工况下对设计制作好的气动弹性模型进行风洞试验的单塔结构示意图;
图6是各风速工况下对设计制作好的气动弹性模型进行风洞试验的塔线体系结构示意图;
图7是输电塔风致机械能最大值的分布图;
图8是不同荷载激发的输电塔风致机械能最大值的分布示意图;
图9是不挂线的输电塔时,风速对输电塔风致机械能最大值的影响示意图;
图10是挂线的输电塔时,风速对输电塔风致机械能最大值的影响示意图;
图11是梯度风高度对不挂线时的输电塔风致机械能最大值的影响示意图;
图12是梯度风高度对挂线时的输电塔风致机械能最大值的影响示意图;
图13是梯度风高度对输电塔风塔体荷载致机械能最大值的影响示意图;
图14是梯度风高度对输电塔传递荷载风致机械能最大值的影响示意图;
图15是本发明方法流程图。
具体实施方式
下面结合附图对本发明的具体实施方式以及工作原理作进一步详细说明。
在本实施例中,以1 000kV苏通大跨越输电线路为例,进行大跨越输电塔风致机械能传递计算。并且基于在西南交通大学XNJD-3风洞实验室进行的大气边界层风洞试验,测得风阻力和动力参数。使用正交分解技术(POD)确定其参振模态。
苏通大跨越工程介绍:
苏通大跨越属于淮南-南京-上海(北环)走向的1 000kV双回路交流输变电线路试验示范工程的关键跨越。
一种塔线体系下大跨越输电塔风致机械能的传递分析方法,在本实施例中,结合图15可以看出,分析方法具体步骤为:
S1:基于大跨越输电塔搭建塔线体系,并获取塔线体系物理参数;
所述塔线体系物理参数至少包括输电塔布置方案、导线地线的跨度垂度参数、导线型号参数、地线型号参数。
在本实施例中,采用耐张塔-直线塔-直线塔-耐张塔布置,详见图1和图2;整个耐张段约5公里。导线与地线的跨度和垂度如表1所示。JLHAl/G6A-500/400钢芯铝绞线物理参数如表2所示。OPGW-350T地线物理参数如表3所示。其它设计参数见文献“赵爽,晏致涛,李正良等.1000kV苏通大跨越输电塔线体系气弹模型的风洞试验研究[J].中国电机工程学报,2018,38(17):5257-5265.”,相关的行业专有名词或概念见文献“张殿生.电力工程高压送电线路设计手册(第2版)[M].中国:中国电力出版社,2003.”。
表1导线与地线的跨度和垂度(m)
类别 | 边跨 | 中间跨 | 边跨 |
档距 | 1113 | 2600 | 1344 |
地线垂度 | 42.5 | 231.5 | 63.4 |
导线垂度 | 45 | 245.9 | 67.1 |
表2 JLHA1/G6A-500/400钢芯铝绞线物理参数
表3 OPGW-350T物理参数
在本实施例中,通过模拟现行荷载规范(详见中华人民共和国住房和城乡建设部.GB 50009-2012,建筑结构荷载规范[S].北京:中国建筑工业出版社,2012.)中B类地貌的风场特性,在西南交通大学的XNJD-3风洞实验室进行苏通大跨越模型的大气边界层风洞试验。分别采用刚性模型和气动弹性模型确定下文计算输电塔风致机械能需要的阻力系数和动力参数。
S2:基于塔线体系,设定模型缩小比例,构建刚性模型,基于刚性模型进行风洞试验,测得风阻力参数;
在本实施例中,刚性模型按照1∶120的几何相似比设计,输电塔构件为圆形截面,采用不锈钢管制作输电塔刚性模型,通过增加钢管厚度保证其刚度。制作好的刚性模型详见图3。
步骤S2中,测量所述风阻力参数时:刚性模型采用五分量应变天平测量气动力,并设定采样时间、采样频率、来流风向、试验风速转化为实际基本风速比例;
在本实施例中,刚性模型采用五分量应变天平测量气动力,采样时间10s,采样频率800Hz,试验风速转化为实际基本风速为40.156m/s,来流风向垂直于线路走向(顺风向)。
将输电塔的塔身和横担沿高度分段,通过扣除法确定各段的风阻力,其中第i段的阻力系数表达式为:
式中,α为地面粗糙度指数。通过刚性模型的风洞试验确定输电塔各段的阻力系数分布如图4。
基于塔线体系,设定模型缩小比例,构建大跨越气动弹性模型,基于大跨越气动弹性模型进行风洞试验,测得动力参数;使用正交分解技术(POD)确定其参振模态;
荷载规范的梯度风是一种理想模型,由于风洞上壁面对风场具有剪切效应,难以模拟梯度风高度以上的平流层,因此实际试验中不考虑梯度风高度。苏通大跨越气动弹性模型的风洞试验涉及到两个关键问题:1)钢管混凝土材料的非线性影响;2)梯度风高度对风振响应的影响。在各风速工况下对设计制作好的气动弹性模型进行风洞试验研究,详见图5和图6。
对试验数据联合运用经验模态分解和随机减量方法,得到模型顺风向的基频和对应的阻尼比参数详见表4。
表4输电塔不挂线与挂线时的基频与对应阻尼比参数
S3:以输电线大变形的平均风状态作为计算初始条件,基于随机振动理论和线性叠加原理,设定工况,并得出每一种工况下输电塔风致机械能最大值的计算模型;
S4:根据步骤S3中得到的输电塔风致机械能最大值的计算模型,根据工况,计算与工况对应的输电塔风致机械能最大值;
S5:分析影响因子输电线、风速和梯度风高度对风致机械能分量的影响以及子结构之间能量的传递关系。
输电塔的机械能由动能和弹性势能构成,由于风荷载和阻尼耗能的共同作用,机械能的大小会发生变化。以下将建立单塔和塔线体系下输电塔风致机械能最大值的计算模型。
一、风振方程
来流风荷载可以分解为平均风荷载和脉动风荷载两部分,输电塔在风荷载作用下的振动方程矩阵表达式为:
式中,M、C、K分别为质量矩阵、阻尼矩阵、刚度矩阵;L为节点从属面积上作用单位风压荷载时节点等效力组成的转换矩阵;和pd分别为节点处的平均风压荷载和脉动风压荷载;分别为脉动风荷载作用下节点顺风向的加速度、速度和位移;为平均风荷载作用下节点顺风向的位移。和pd的表示为:
式中,·*表示矩阵的对应元素相乘;v′为节点的脉动风速。pd作用下节点的位移响应包含背景响应和共振响应,可以采用随机振动理论求解。
二、不挂线时风致机械能最大值
风致机械能由其平均分量、背景分量和共振分量组成。平均分量和背景分量具有准静态特性,由外力做功转化,不引起阻尼耗能。共振分量等于动能和弹性势能之和,并且动能幅值和弹性势能幅值相同,两者的相位差为90°,会引起阻尼耗能。当背景位移达到最大时,共振位移很小,若此时共振位移为0,阻尼耗能可以忽略不计。因此,当塔顶背景位移达到最大,共振动能达到幅值,共振位移为0时,输电塔风致机械为最大值。由于输电塔的最大位移响应是背景位移和共振位移两者处于相对较大值时的组合,输电塔风致机械能最大值发生时,并不是风致位移响应的最不利情况。
对苏通大跨越气动弹性模型多个测点的风致位移响应使用正交分解技术(POD),结果表明1阶振型的位移能量贡献度高达97.1%,下文计算仅考虑1阶模态参振。
不挂线时输电塔风致机械能最大值的表达式为:
塔体的平均风荷载在对应的平均位移下做的实功为:
塔体的平均风荷载在峰值背景位移下做的虚功为:
塔体的峰值背景风荷载在对应的峰值背景位移下做的实功为:
塔体的峰值背景风荷载在平均位移下做的虚功为:
塔体的共振响应激发风致动能的最大值为:
式中,上标“R”表示实功;上标“I”表示虚功;下标“t1”表示不挂线时的输电塔;下标“t”表示挂线与不挂线时的输电塔;为平均风荷载对结构做功;Wb,max为背景风荷载对结构做功的最大值;Er,max为共振响应激发风致动能的最大值。对于线弹性结构,由功的互等定理确定XP为KXP=L的解。
采用LRC方法确定等效背景风荷载,峰值背景风荷载的表达式为:
Pb,max=gbρPr·*σP (12)
式中,gb为结构的背景峰值因子,按3.5取值;ρPr为荷载响应相关系数;σP为脉动风压的标准差。
ρPr的表达式为:
式中,上横线表示对时间求平均;./表示矩阵的对应元素相除;Ir为节点响应的影响线,本文取塔顶顺风向位移影响线。节点背景响应标准差的表达式为:
式中,diag( )表示提取矩阵的对角元素;SPP为节点处的脉动风压互谱矩阵。公式(13)中塔顶处的σb通过公式(15)对应节点确定。
sPP中采用Davenport风速谱和Shiotani相干函数。Davenport风速谱的表达式为:
Shiotani相干函数为:
coh(x1,z1,x2,z2)=cohx(x1,x2)cohz(z1,z2) (19)
式中,coh(x1,z1,x2,z2)为相干函数;cohx(x1,x2)为水平相干函数;cohz(z1,z2)为竖向相干函数。
σP的表达式为:
式中,σv′ 2为脉动风速的方差;n为脉动风速的频率。
对于具有钢管混凝土的超高输电塔,塔顶顺风向位移影响线可采用以下表达式:
式中,Δ为单位荷载作用在塔顶时产生的塔顶位移,与输电塔的结构特性有关;H为塔高。
仅考虑1阶模态参振时,公式(11)可以改写为以下表达式:
式中,ω1为1阶振型的圆频率;gr1为1阶振型的共振峰值因子,通过文献“Davenport A G.Gust loading factors[J].Journal of the Structural Division,1967,93(ST3):24-26.”提出的首次穿越理论确定;为1振型的广义质量;
共振响应下1振型广义坐标标准差的表达式为:
式中,ζ1为结构第1阶振型的阻尼比;φ1为1阶振型。
三、挂线时风致机械能最大值
塔线体系下,输电塔的风荷载来源主要由两部分构成,一部分是杆塔自身风荷载(塔体荷载),另一部分是输电线通过绝缘子串传递到杆塔的风荷载(传递荷载)。来流风荷载引起输电线的风致响应,在气动阻尼作用下,共振响应可以忽略,从而以准静态的背景响应为主,ASCE规范(详见ASCE NO.74.Guidelines for Electrical Transmission LineStructural Loading[S].Reston:American Society of Civil Engineering,2009:116-117.)同样不考虑输电线共振响应。输电线的风致响应表现为在平均响应的平衡位置处作小幅度、往返的振动。可以假定输电线在脉动风荷载作用下满足小位移,荷载-响应关系近似为线性。此时,以输电线的平均风和自重状态作为计算的初始条件,这样可以基于线性叠加原理,采用荷载响应相关法(LRC)求解输电线的背景响应。
需要强调的是,本实施例分析的机械能达到最大值时的状态是一种可能性状态,即是该可能性小,对此状态研究得出的结论仍不失一般性。如此,当塔顶位移背景分量和输电线在挂点处的位移背景分量同时达到最大时,输电塔风致机械能达到最大值。
挂线时输电塔风致机械能最大值的表达式为:
塔体的平均风荷载和输电线传递的平均风荷载在各自对应的平均位移下做的实功总和为:
塔体的平均风荷载在对应的平均位移下做的实功为:
输电线传递的平均风荷在对应的平均位移下做的实功为:
塔体的平均风荷载和输电线传递的平均风荷载在峰值背景位移下做的虚功总和为:
塔线体系下,输电塔的峰值背景位移包括输电塔峰值背景风荷载和输电线峰值背景风荷载两部分激发的输电塔位移。塔体的平均风荷载在峰值背景位移下做的虚功为:
塔体的峰值背景风荷载在对应的峰值背景位移下做的实功总和为:
输电线传递的峰值背景风荷载在对应的峰值背景位移下做的实功总和为:
塔体的峰值背景风荷载和输电线传递的峰值背景风荷载在平均位移下做的虚功总和为:
塔线体系下,输电塔的平均位移包括输电塔平均风荷载和输电线平均风荷载两部分激发的输电塔位移。塔体的峰值背景风荷载在平均位移下做的虚功为:
输电线传递的峰值背景风荷载在平均位移下做的虚功为:
塔体的共振响应激发风致动能的最大值为:
式中,下标“t2”表示挂线时的输电塔。W1为塔体荷载做功;W2为传递荷载做功;和Rb,max分别为传递荷载的平均分量和峰值背景分量;IR为KIR=E的解,E为单位矩阵。根据前面叙述的计算思路确定输电线的平均风荷载和峰值背景风荷载,计算输电线峰值背景风荷载时取输电线挂点处的顺风向位移影响线。
根据上述描述,可以求取苏通大跨越输电塔的风致机械能。采用风洞试验确定的阻力系数、输电塔自由振动频率、输电塔的阻尼比,并假定这些参数受梯度风高度的影响小。采用ANSYS有限元模型确定计算需要的挡风面积、质量、计算点坐标、输电塔振型、影响线。导线的结构阻尼取0.4%,气动阻尼由文献“Loredo-Souz A M,Davenport A G.Theeffects of high winds on transmission lines[J].Journal of Wind Engineeringand Industrial Aerodynamics,1998,74-76:987-994.”中的公式进行确定。风场特性由荷载规范确定
有限元计算结果与试验结果对比为:
采用ANSYS软件建立苏通大跨越的有限元模型,输电塔采用Beam188单元模拟,导线采用Link10单元模拟,绝缘子串采用Link8单元模拟。截面特性和材料特性均按照实际设计参数赋予。风速40.15m/s时,采用有限元模型计算平均风荷载作用下单塔和塔线体系的塔顶位移,并与风洞试验的对应统计值(通过相似比转化为原型值)进行对比,如表5所示。数据对比吻合度好,验证了有限元模型的有效性。因而,采用此有限元模型输出的参数作为后续数值计算是可行的。
表5风速40.15m/s时由不同方法确定的塔顶平均位移
类别 | 气动弹性模型/mm | 有限元模型/mm | 误差/% |
单塔 | 228.34 | 235.15 | 2.98 |
塔线体系 | 773.22 | 791.41 | 2.35 |
风致机械能传递机制
为了研究输电塔的风致机械能传递机制,下文分析了不同工况下输电塔风致机械能的平均分量、背景分量和共振分量的变化情况,阐述了该型结构的风致机械能传递规律。
1)输电线的影响
基本风速时,输电塔挂线与不挂线两种情况下,输电塔风致机械能最大值的分布如图7所示,各分量的比例如表6所示。
表6输电塔风致机械能最大值时各分量的比例(%)
类别 | 平均分量 | 背景分量 | 共振分量 |
不挂线时的输电塔 | 65.46 | 18.70 | 15.84 |
挂线时的输电塔 | 76.90 | 17.51 | 5.59 |
不挂线时,输电塔风致机械能最大值以平均分量为主,脉动分量比例显著,不能忽视。其中,脉动分量以背景分量为主,共振分量比例比背景分量的略小。挂线后,输电塔风致机械能最大值的准静态分量增大,共振分量减小。输电塔风致机械能最大值以平均分量为主,背景分量比例显著,共振分量比例较小。对比挂线前的输电塔风致机械能最大值,挂线后的平均分量增大2.106倍,脉动分量增大0.768倍,其中背景分量增大1.476倍,共振分量减小0.068倍。
图7中由塔体荷载和传递荷载激发输电塔风致机械能最大值的成分如图8所示,各分量的比例如表7所示。输电塔挂线后阻尼比增大,共振能量比挂线前减小。与图7的单塔对比,图8塔体荷载作用下,输电塔风致机械能最大值的平均分量增大0.398倍,背景分量增大1.048倍,共振分量减小0.068倍。塔体荷载作用下,输电塔风致机械能最大值以平均分量为主,脉动分量比例显著,不能忽视。其中,脉动分量以背景分量为主。传递荷载作用下输电塔风致机械能最大值以平均分量为主,背景分量的比例较小。传递荷载作用下输电塔风致机械能最大值的平均分量是塔体荷载作用下的1.221倍,背景分量是塔体荷载作用下的0.209倍。
塔线耦联对输电塔风致机械能传递的影响表现为:由于横担挂点位置高,传递荷载对输电塔风致机械能的贡献显著,向输电塔传递了大量的准静态机械能,该部分能量大小接近于由塔体荷载激发的,十分显著。此外,塔线的耦联振动会增大整个体系的阻尼比,抑制输电塔的共振机械能。对于超高大跨越输电线路,输电线的行为起到了控制作用,输电塔的风致机械能以准静态分量为主。
表7不同荷载激发输电塔风致机械能最大值时各分量的比例(%)
类别 | 平均分量 | 背景分量 | 共振分量 |
塔体荷载 | 63.30 | 26.48 | 10.22 |
传递荷载 | 93.32 | 6.68 | - |
2)风速的影响
挂线、不挂线两种情况下,风速对超高输电塔抖振机械能最大值的影响如图9、图10所示。为基本风速,与成正比,并且对于线弹性结构,与成正比,因此与成正比。σP与成正比,ρPr受的影响小,同样Pb,max与XPPb,max成正比,因此近似与成正比。增大时,ζ1增大,ω1减小,输电塔1阶模态频率处的Spp与的更高次方成正比,因此Er,max也与的更高次方成正比。这样,Er,max受的影响要比和大。
图9和图10中,风速对输电塔风致机械能传递的影响表现为:随着风速的增大,输电塔风致机械能最大值将会增大,但是各成分增大的速率不一致。其中,共振分量增大最快,平均分量增大最慢。随着风速的增大,输电塔将有可能发生动力破坏。当输电塔挂线后,整个体系的阻尼比增大,抑制了共振效应的增强,在一定风速范围内,输电塔的风致机械能仍然以准静态分量为主。
3)梯度风高度的影响
基本风速时,梯度风高度对输电塔风致机械能最大值的影响如图11所示,各分量的比例如表所示。与不考虑梯度风高度时的输电塔风致机械能最大值对比,考虑梯风高度后,不挂线时输电塔风致机械能最大值的平均分量减小0.126倍,背景分量减小0.538倍,共振分量减小0.575倍;挂线时输电塔风致机械能最大值的平均分量减小0.089倍,背景分量减小0.515倍,共振分量减小0.578倍;由塔体荷载激发输电塔风致机械能最大值的平均分量减小0.09倍,背景分量减小0.578倍,共振分量减小0.578倍;由传递荷载激发输电塔风致机械能最大值的平均分量减小0.088倍,背景分量减小0.215倍。
梯度风高度对输电塔风致机械能传递的影响表现为:考虑梯度风高度后,由于风场整体的湍流特性变弱,输电塔风致机械能最大值将会减小,但是各成分减小的速率不一致。其中,共振分量减小最快,平均分量减小最慢。由于输电线仅有小部分区域超出梯度风高度范围,梯度风高度对输电线风振的影响要比输电塔小,进而梯度风高度对挂线时输电塔的影响比不挂线时的小。考虑梯度风高度后,输电塔的风致机械能以准静态分量为主,共振分量偏小。
表8梯度风高度下输电塔风致机械能最大值时各分量的比例(%)
类别 | 平均分量 | 背景分量 | 共振分量 |
不挂线时的输电塔 | 78.83 | 11.90 | 9.27 |
挂线时的输电塔 | 86.60 | 10.49 | 2.91 |
塔体荷载 | 78.81 | 15.29 | 5.9 |
传递荷载 | 94.19 | 5.81 | - |
经上述说明后得到以下结论:采用风洞试验数据,结合本文提出的计算模型分析苏通大跨越输电塔风致机械能传递机制,结论如下:
1)输电塔的阻力系数沿高度变化,横担处的阻力系数比同高度塔身的小。随着风速的增大,输电塔频率减小,阻尼比增大。
2)当输电塔顶部的背景位移达到最大,并且由共振响应激发的弹性势能为0,动能达到最大,此时的阻尼耗能为0,输电塔的风致机械能达到最大状态。
3)超高大跨越输电塔的风致机械能传递机制为:传递荷载向输电塔传递了大量的准静态能量,该部分能量的大小接近于由塔体荷载激发的,十分显著。由于输电线的行为起主要作用,输电塔的风致机械能以准静态分量为主,共振分量偏小。风速增大时,输电塔风致机械能各成分增大的速率不一致,塔线耦联振动会抑制共振机械能的增长幅度。考虑梯度风高度后,输电塔风致机械能各成分减小的速率不一致,输电线的存在会削弱梯度风高度对能量的影响。
应当指出的是,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的普通技术人员在本发明的实质范围内所做出的变化、改性、添加或替换,也应属于本发明的保护范围。
Claims (3)
1.一种塔线体系下大跨越输电塔风致机械能的传递分析方法,其特征在于:具体步骤为:
S1:基于大跨越输电塔搭建塔线体系,并获取塔线体系物理参数;
S2:基于塔线体系,设定模型缩小比例,构建刚性模型,基于刚性模型进行风洞试验,测得风阻力参数;
基于塔线体系,设定模型缩小比例,构建大跨越气动弹性模型,基于大跨越气动弹性模型进行风洞试验,测得动力参数;
使用正交分解技术POD确定其参振模态;
S3:以输电线大变形的平均风状态作为计算初始条件,基于随机振动理论和线性叠加原理,设定工况,并得出每一种工况下输电塔风致机械能最大值的计算模型;
S4:根据步骤S3中得到的输电塔风致机械能最大值的计算模型,根据工况,计算与工况对应的输电塔风致机械能最大值;
S5:分析影响因子输电线、风速和梯度风高度对风致机械能分量的影响以及子结构之间能量的传递关系;
所述输电塔风致机械能最大值的计算模型为:
塔体的平均风荷载和输电线传递的平均风荷载在各自对应的平均位移下做的实功总和为:
塔体的平均风荷载在对应的平均位移下做的实功为:
输电线传递的平均风荷在对应的平均位移下做的实功为:
塔体的平均风荷载和输电线传递的平均风荷载在峰值背景位移下做的虚功总和为:
塔线体系下,输电塔的峰值背景位移包括输电塔峰值背景风荷载和输电线峰值背景风荷载两部分激发的输电塔位移;塔体的平均风荷载在峰值背景位移下做的虚功为:
输电线传递的平均风荷载在峰值背景位移下做的虚功为:
塔体的峰值背景风荷载和输电线传递的峰值背景风荷载在各自对应的峰值背景位移下做的实功总和为:
塔体的峰值背景风荷载在对应的峰值背景位移下做的实功总和为:
输电线传递的峰值背景风荷载在对应的峰值背景位移下做的实功总和为:
塔体的峰值背景风荷载和输电线传递的峰值背景风荷载在平均位移下做的虚功总和为:
塔线体系下,输电塔的平均位移包括输电塔平均风荷载和输电线平均风荷载两部分激发的输电塔位移;塔体的峰值背景风荷载在平均位移下做的虚功为:
输电线传递的峰值背景风荷载在平均位移下做的虚功为:
塔体的共振响应激发风致动能的最大值为:
式中,下标“t2”表示挂线时的输电塔;W1为塔体荷载做功;W2为传递荷载做功;和Rb,max分别为传递荷载的平均分量和峰值背景分量;IR为KIR=E的解,E为单位矩阵;分别通过对输电线水平档距内的平均风荷载和峰值背景风荷载求和计算和Rb,max;
2.根据权利要求1所述的塔线体系下大跨越输电塔风致机械能的传递分析方法,其特征在于:所述塔线体系物理参数至少包括输电塔布置方案、导线地线的跨度垂度参数、导线型号参数、地线型号参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010488810.1A CN111651920B (zh) | 2020-06-02 | 2020-06-02 | 塔线体系下大跨越输电塔风致机械能的传递分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010488810.1A CN111651920B (zh) | 2020-06-02 | 2020-06-02 | 塔线体系下大跨越输电塔风致机械能的传递分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111651920A CN111651920A (zh) | 2020-09-11 |
CN111651920B true CN111651920B (zh) | 2023-03-21 |
Family
ID=72351073
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010488810.1A Active CN111651920B (zh) | 2020-06-02 | 2020-06-02 | 塔线体系下大跨越输电塔风致机械能的传递分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111651920B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113639945B (zh) * | 2021-06-28 | 2024-02-09 | 上海宇航系统工程研究所 | 基于经验模态分解的航天器随机振动试验条件设计方法 |
CN113688774B (zh) * | 2021-09-03 | 2023-09-26 | 重庆大学 | 基于深度学习的高层建筑风致响应预测、训练方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102692309A (zh) * | 2012-05-10 | 2012-09-26 | 南方电网科学研究院有限责任公司 | 一种应用于台风风场的输电线风洞试验系统及方法 |
CN108073758A (zh) * | 2017-10-17 | 2018-05-25 | 广东电网有限责任公司电力科学研究院 | 一种输电塔线风致振动响应仿真方法及装置 |
CN108959742A (zh) * | 2018-06-20 | 2018-12-07 | 重庆科技学院 | 大跨越输电塔线体系气动弹性模型设计方法 |
CN110298101A (zh) * | 2019-06-24 | 2019-10-01 | 国网浙江省电力有限公司电力科学研究院 | 一种耦合局地风场的输电线路风致响应有限元分析方法 |
CN110378050A (zh) * | 2019-07-25 | 2019-10-25 | 重庆科技学院 | 一种输电塔顺风向气动阻尼比计算方法 |
-
2020
- 2020-06-02 CN CN202010488810.1A patent/CN111651920B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102692309A (zh) * | 2012-05-10 | 2012-09-26 | 南方电网科学研究院有限责任公司 | 一种应用于台风风场的输电线风洞试验系统及方法 |
CN108073758A (zh) * | 2017-10-17 | 2018-05-25 | 广东电网有限责任公司电力科学研究院 | 一种输电塔线风致振动响应仿真方法及装置 |
CN108959742A (zh) * | 2018-06-20 | 2018-12-07 | 重庆科技学院 | 大跨越输电塔线体系气动弹性模型设计方法 |
CN110298101A (zh) * | 2019-06-24 | 2019-10-01 | 国网浙江省电力有限公司电力科学研究院 | 一种耦合局地风场的输电线路风致响应有限元分析方法 |
CN110378050A (zh) * | 2019-07-25 | 2019-10-25 | 重庆科技学院 | 一种输电塔顺风向气动阻尼比计算方法 |
Non-Patent Citations (5)
Title |
---|
1000kV特高压输电塔线体系风荷载传递机制风洞试验研究;谢强 等;《中国电机工程学报》;20130105;第33卷(第01期);第109-116页 * |
1000kV苏通大跨越输电塔线体系气弹模型的风洞试验研究;赵爽 等;《中国电机工程学报》;20180905;第38卷(第17期);第5257-5265,5323页 * |
Parameter analysis of wind-induced vibration of electrical transmission line models;YOU Yi 等;《2nd International Conference on Civil, Transportation and Environmental Engineering》;20170731;第129-134页 * |
Wind-induced vibration of UHV transmission tower line system: Wind tunnel test on aero-elastic model;Qiang Xie 等;《Journal of Wind Engineering & Industrial Aerodynamics》;20171017;第171卷;第219-229页 * |
大跨越高压输电塔线体系风致振动的研究与进展;郭惠勇 等;《河海大学学报(自然科学版)》;20080125(第01期);第106-111页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111651920A (zh) | 2020-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | The effect of amplitude-dependent damping on wind-induced vibrations of a super tall building | |
CN108959742B (zh) | 大跨越输电塔线体系气动弹性模型设计方法 | |
Ke et al. | A new methodology for analysis of equivalent static wind loads on super-large cooling towers | |
CN111651920B (zh) | 塔线体系下大跨越输电塔风致机械能的传递分析方法 | |
CN101866395B (zh) | 一种导线舞动时输电线路杆塔挂点荷载计算方法 | |
CN105354377A (zh) | 一种确定输电杆塔的脉动风振荷载的方法 | |
CN106096257B (zh) | 一种非线性索单元分析方法及系统 | |
CN104408295B (zh) | 一种大跨桥梁下部结构风‑浪耦合作用荷载数值模拟方法 | |
Ma et al. | Investigations of aerodynamic effects on streamlined box girder using two-dimensional actively-controlled oncoming flow | |
CN111125612B (zh) | 一种架空导线覆冰舞动幅值的计算方法 | |
CN106096105A (zh) | 输电线路风偏瞬态响应计算方法 | |
CN101655588B (zh) | 一种光纤复合架空地线opgw光缆覆冰机械性能模拟方法 | |
CN104200102A (zh) | 一种基于cfd的光伏阵列风荷载预测方法 | |
CN114611365A (zh) | 基于脉动风影响下输电杆塔结构的动力学分析方法及系统 | |
CN109543258A (zh) | 一种基于模态叠加法的风电机组塔架应力状态计算方法 | |
CN108268729A (zh) | 输电导线的弹性模量频率灵敏度分析方法及系统 | |
Vemuri et al. | Wind Analysis of Tall Buildings Using Codal Provisions, Stochastic Approach and CFD–A Comparative Study | |
CN101261651A (zh) | 主梁断面气动自激力的全过程数值模拟方法 | |
Ma et al. | Case study of wind-induced performance and equivalent static wind loads of large-span openable truss structures | |
CN110363677A (zh) | 一种覆冰输电线路找形方法 | |
Li et al. | A reconstruction method for structural stress distribution of wind turbine tower using optimised mathematical model | |
Ke et al. | Wind-induced responses of super-large cooling towers | |
CN109340060A (zh) | 一种基于模态叠加法的风电机组塔架振动状态计算方法 | |
CN110909512B (zh) | 一种基于非定常方法的双分裂背风子导线气动荷载仿真方法 | |
Gu et al. | Parametric vibration of stay cables under axial narrow-band stochastic excitation |
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 |