CN114611365A - 基于脉动风影响下输电杆塔结构的动力学分析方法及系统 - Google Patents

基于脉动风影响下输电杆塔结构的动力学分析方法及系统 Download PDF

Info

Publication number
CN114611365A
CN114611365A CN202210288857.2A CN202210288857A CN114611365A CN 114611365 A CN114611365 A CN 114611365A CN 202210288857 A CN202210288857 A CN 202210288857A CN 114611365 A CN114611365 A CN 114611365A
Authority
CN
China
Prior art keywords
transmission tower
wind
pulsating
tower
wind speed
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.)
Pending
Application number
CN202210288857.2A
Other languages
English (en)
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.)
Electric Power Research Institute of State Grid Zhejiang Electric Power Co Ltd
Original Assignee
Electric Power Research Institute of State Grid Zhejiang Electric Power Co Ltd
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 Electric Power Research Institute of State Grid Zhejiang Electric Power Co Ltd filed Critical Electric Power Research Institute of State Grid Zhejiang Electric Power Co Ltd
Priority to CN202210288857.2A priority Critical patent/CN114611365A/zh
Publication of CN114611365A publication Critical patent/CN114611365A/zh
Pending legal-status Critical Current

Links

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
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/04—Power grid distribution networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Wind Motors (AREA)

Abstract

本发明公开了一种基于脉动风影响下输电杆塔结构的动力学分析方法及系统。本发明动力学分析方法采用的技术方案包括：根据某转角输电杆塔的工程实际数据建立输电杆塔有限元模型；针对输电杆塔的结构特点，对输电杆塔有限元模型进行模态分析；结合快速傅里叶变换的谐波叠加法，通过仿真软件模拟输电杆塔的空间脉动风速；将动态风载荷沿0°风向角作用在塔顶节点上，对输电杆塔的动态响应方程进行时域分析求解，得到不同风向角下塔顶的位移响应。本发明解决了现有输电杆塔在自然灾害以及极端恶劣天气工况下的稳定性和安全性问题，有助于提高输电线路的安全性和稳定性，有助于提高输电线路的抗风性能，保障电网的安全运行。

Description

基于脉动风影响下输电杆塔结构的动力学分析方法及系统
技术领域
本发明属于输电杆塔结构分析领域，具体地说是一种基于脉动风影响下输电杆塔结构的动力学分析方法及系统。
背景技术
输变电塔线体系是电力系统的重要组成部分，是电力系统安全可靠运行的基石，将电力安全稳定地传输到各行各业，保证社会稳定有序地开展作业。作为电力系统的大动脉，开展其稳定性及安全性分析尤其重要。输变电塔线系统工作于外部曝露环境，常年受日光、风雨等天气影响，特别在自然灾害以及极端恶劣天气的工况下，对电力系统的稳定和安全将会带来影响，严重的将导致输电线路杆塔振动、线路舞动、线路器件损坏甚至失稳倒塌，不仅导致电力系统运行瘫痪，也会给社会带来巨大的经济损失。全面了解掌握输电杆塔在受风载荷、地震载荷等各种载荷综合作用下的动力特性，是提高输电线路安全性与稳定性的重中之重。
近年来，随着有限元分析的高速发展，其建模仿真的高效、准确等特点受到许多行业的青睐，通过建立输电线路杆塔结构准确的有限元模型，研究脉动风载荷下的输电线路杆塔的动力学响应，对于输电杆塔结构的风振响应设计以及电网的安全稳定性具有重要的工程意义。目前主流的研究方案主要综合了包含有限元、随机振动、结构分析理论等多类学科理论，主要针对输电线路动力计算模型、动力试验以及动力特性分析等方面研究输电杆塔线体系。而在众多灾害中，由风引起的振动灾害最为严重，甚者会导致倒塔和断线，因此，输电杆塔抗风研究是输电线路安全工作的关键。
发明内容
针对上述现有技术的不足，本发明提供一种基于脉动风影响下输电杆塔结构的动力学分析方法，以解决现有输电杆塔在自然灾害以及极端恶劣天气工况下的稳定性和安全性问题，提高输电线路的安全性和稳定性，提高输电线路的抗风性能，以保障电网的安全运行。
为此，本发明采用如下的技术方案：一种基于脉动风影响下输电杆塔结构的动力学分析方法，其包括步骤：
a)根据某转角输电杆塔的工程实际数据建立输电杆塔有限元模型；
b)针对输电杆塔的结构特点，对输电杆塔有限元模型进行模态分析；
c)结合快速傅里叶变换的谐波叠加法，通过仿真软件模拟输电杆塔的空间脉动风速；
d)将动态风载荷沿0°风向角作用在塔顶节点上，对输电杆塔的动态响应方程进行时域分析求解，得到不同风向角下塔顶的位移响应。
进一步地，步骤a)中，转角输电杆塔的工程实际数据包括塔高、呼高、档距和塔基尺寸，角钢与角钢之间用螺栓和连接板连接，忽略螺栓和连接板的几何尺寸效应，将角钢构件设置为刚接，采用Ansys有限元软件中Beaml88单元对输电杆塔的角钢构件进行建模。
进一步地，步骤b)中，所述的模态分析包括模态求解、扩展模态和结果观察，具体过程如下：
对输电杆塔有限元模型进行模态分析处理，模态分析的实质是计算输电杆塔结构振动特征方程的特征值和特征向量，根据结构动力学基本方程：
Figure BDA0003559348020000021
上式中：M为输电杆塔的整体刚度矩阵；C为输电杆塔的整体阻尼矩阵；K为输电杆塔的整体刚度矩阵；
Figure BDA0003559348020000022
为结构加速度向量；
Figure BDA0003559348020000023
为结构速度向量；X为结构位移向量；F(t)为加载在结构上的载荷；t为时间；
由于工程结构中可不考虑阻尼对结构频率和模态的影响，忽略结构阻尼，在无外荷载的作用下输电杆塔结构自由振动动力学方程为：
Figure BDA0003559348020000024
当输电杆塔结构在固有频率下发生振动时，其输电杆塔结构速度和加速度向量表示为：
Figure BDA0003559348020000025
Figure BDA0003559348020000026
上式中，
Figure BDA0003559348020000027
为结构节点位移的节点振幅矩阵；ω为固有频率；
Figure BDA0003559348020000028
为结构发生振动的相位；
将上式(3)与(4)带入式(2)得：
Figure BDA0003559348020000029
按自由振动理论，r阶自由度系统的自由振动方程式有r个固有频率ωi，i＝1，2，3，…，n，并且由频率行列式决定，即：
|K-ω2M|＝0， (6)
当输电杆塔有限元模型的自由度为y时，式(6)具有y个特征值，记λj为第j个特征值，其平方根ωj是结构的第j阶固有频率，并且
Figure BDA0003559348020000031
是相应的第j阶特征向量，特征向量即为输电杆塔结构的模态振型，是输电杆塔结构在第j阶模态振型下的变形状态。
进一步地，步骤c)中，模拟输电杆塔空间脉动风速的方法为：
将风荷载视为静力风/平均风与动力风/脉动风的共同作用，在任意时刻高度Z处的风速v(z，t)表示为平均风速与脉动风速之和：
Figure BDA0003559348020000032
式中，
Figure BDA0003559348020000033
为Z高度处的平均风速，vf(z，t)为Z高度处的脉动风速；
平均风速与高度之间的变化规律称为风剖面，风剖面用指数律模拟：
Figure BDA0003559348020000034
式中，
Figure BDA0003559348020000035
为任一点平均风速，z为任一点的高度；
Figure BDA0003559348020000036
为10m处基准平均风速；α为地面粗糙度系数；
脉动风是一个均值为零的高斯随机过程，功率谱密度是对其均方值的量度，脉动风自功率谱的一种表达式Davenport谱为：
Figure BDA0003559348020000037
式中，
Figure BDA0003559348020000038
k为与地面粗糙度有关的系数，γ为脉动风频率，运用基于快速傅里叶变换技术的谐波叠加法：
Figure BDA0003559348020000039
p＝0，1，2，…2N，j＝0，1，2，…n；
其中，q是p/2N的余数，q＝0，1，2，…2n-1；
hjm(qΔt)由下式给出：
Figure BDA00035593480200000310
其中：
Figure BDA00035593480200000311
其中，m为风谱在频域内划分的分部数；Δt为时间增量；Δω为频率增量；S(ω)为互功率谱密度矩阵；Gjm为快速傅里叶变换中分解矩阵的元素；
Figure BDA0003559348020000041
为两个不同作用点之间的相位角；l为步长；
计算得到脉动风速时程图和自功率谱对比图；其中，自功率谱对比图反映了通过脉动风速时程图形成的功率谱与达文波特功率谱之间拟合效果的优劣；
输电杆塔每段单个节点的动态风荷载计算公式为：
ωs＝ωkμsA0，
其中，ωk为对应高度的脉动风压；μs为转角输电杆塔的风荷载体型系数，A0为迎风投影面积；
脉动风压ωk通过脉动风速由公式ωk＝0.5ρv2转换得到，ρ表示迎风面的方向系数；v表示脉动风速；
风荷载体型系数μs由公式：μs＝1.3(1+η)计算得到，其中η为背风面的风载荷降低系数，最后模拟出塔顶的脉动风载荷时程曲线图。
进一步地，步骤d)中，对输电杆塔的动态响应方程进行时域分析求解方法为：
输电杆塔在稳态风的作用下，动态响应方程为：
Figure BDA0003559348020000042
其中，M为输电杆塔的整体质量矩阵，C为输电杆塔的整体阻尼矩阵，K为输电杆塔的整体刚度矩阵；
Figure BDA0003559348020000043
和u分别为节点的加速度、速度以及位移向量，F为输电杆塔线体系所受的动态风荷载；
将动态风载荷F沿0°风向角作用在塔顶上，对输电杆塔的动态响应方程进行时域分析求解。
本发明采用的另一种技术方案为：基于脉动风影响下输电杆塔结构的动力学分析系统，其包括：
有限元模型建立单元，根据某转角输电杆塔的工程实际数据建立输电杆塔有限元模型；
模态分析单元，针对输电杆塔的结构特点，对输电杆塔有限元模型进行模态分析；
空间脉动风速模拟单元，结合快速傅里叶变换的谐波叠加法，通过仿真软件模拟输电杆塔的空间脉动风速；
时域分析求解单元，将动态风载荷沿0°风向角作用在塔顶节点上，对输电杆塔的动态响应方程进行时域分析求解，得到不同风向角下塔顶的位移响应。
本发明具有的有益效果如下：
1)本发明使用工程实际数据建立输电杆塔的有限元模型，有效解决了模型实际性与计算准确性问题。
2)本发明使用基于快速傅里叶变换的谐波叠加法，仿真软件模拟输电杆塔的空间脉动风速，通过FFT技术(快速傅里叶变换技术)改进了其运算效率低的问题。
3)本发明解决了现有输电杆塔在自然灾害以及极端恶劣天气工况下的稳定性和安全性问题，提高了输电线路的安全性和稳定性，提高了输电线路的抗风性能，保障电网的安全运行。
附图说明
图1为本发明动力学分析方法的流程示意图；
图2为本发明输电杆塔的有限元模型图(图2a为立体状态的模型图，图2b为正视图的模型图，图2c为侧视图的模型图，图2d为俯视图的模型图)；
图3为本发明应用例中塔顶处脉动风速时程图(图3a)和自功率谱对比图(图3b)；
图4为本发明应用例中脉动风荷载时程变化曲线图；
图5为本发明应用例中塔顶处位移时程曲线图(风向角0°)。
具体实施方式
下面结合附图对本发明的技术方案作进一步描述。
实施例1
本实施例为一种基于脉动风影响下输电杆塔结构的动力学分析方法，如图1所示，其步骤如下：
a)根据某转角输电杆塔的工程实际数据建立输电杆塔有限元模型，包括简化模型、输电塔架建模和网格划分；
b)针对输电杆塔的结构特点，对输电杆塔有限元模型进行模态分析；
c)结合快速傅里叶变换的谐波叠加法，通过仿真软件模拟输电杆塔的空间脉动风速；
d)将动态风载荷沿0°风向角作用在塔顶节点上，对输电杆塔的动态响应方程进行时域分析求解，得到不同风向角下塔顶的位移响应。
步骤a)中，转角塔的工程实际数据包括塔高、呼高、档距和塔基尺寸，角钢与角钢之间用螺栓和连接板连接，忽略螺栓和连接板的几何尺寸效应，将角钢构件设置为刚接，采用Ansys有限元软件中Beam188单元对输电杆塔的角钢构件进行建模，有限元模型图如图2所示。
步骤b)中，所述的模态分析包括模态求解、扩展模态和结果观察，具体过程如下：
对输电杆塔有限元模型进行模态分析处理，模态分析的实质是计算输电杆塔结构振动特征方程的特征值和特征向量，根据结构动力学基本方程：
Figure BDA0003559348020000061
上式中：M为输电杆塔的整体刚度矩阵；C为输电杆塔的整体阻尼矩阵；K为输电杆塔的整体刚度矩阵；
Figure BDA0003559348020000062
为结构加速度向量；
Figure BDA0003559348020000063
为结构速度向量；X为结构位移向量；F(t)为加载在结构上的载荷；t为时间；
由于工程结构中可不考虑阻尼对结构频率和模态的影响，忽略结构阻尼，在无外荷载的作用下输电杆塔结构自由振动动力学方程为：
Figure BDA0003559348020000064
当输电杆塔结构在固有频率下发生振动时，其输电杆塔结构速度和加速度向量表示为：
Figure BDA0003559348020000065
Figure BDA0003559348020000066
上式中，
Figure BDA0003559348020000067
为结构节点位移的节点振幅矩阵；ω为固有频率；
Figure BDA0003559348020000068
为结构发生振动的相位；
将上式(3)与(4)带入式(2)得：
Figure BDA0003559348020000069
按自由振动理论，r阶自由度系统的自由振动方程式有r个固有频率ωi，i＝1，2，3，…，n，并且由频率行列式决定，即：
|K-ω2M|＝0， (6)
当输电杆塔有限元模型的自由度为y时，式(6)具有y个特征值，记λj为第j个特征值，其平方根ωj是结构的第j阶固有频率，并且
Figure BDA00035593480200000610
是相应的第j阶特征向量，特征向量即为输电杆塔结构的模态振型，是输电杆塔结构在第j阶模态振型下的变形状态。
步骤c)中，模拟输电杆塔空间脉动风速的方法为：
将风荷载视为静力风/平均风与动力风/脉动风的共同作用，在任意时刻高度Z处的风速v(z，t)表示为平均风速与脉动风速之和：
Figure BDA00035593480200000611
式中，
Figure BDA00035593480200000612
为Z高度处的平均风速，vf(z，t)为Z高度处的脉动风速；
平均风速与高度之间的变化规律称为风剖面，风剖面用指数律模拟：
Figure BDA00035593480200000613
式中，
Figure BDA0003559348020000071
为任一点平均风速，z为任一点的高度；
Figure BDA0003559348020000072
为10m处基准平均风速；α为地面粗糙度系数；
脉动风是一个均值为零的高斯随机过程，功率谱密度是对其均方值的量度，脉动风自功率谱的一种表达式Davenport谱为：
Figure BDA0003559348020000073
式中，
Figure BDA0003559348020000074
k为与地面粗糙度有关的系数，γ为脉动风频率，运用基于快速傅里叶变换技术的谐波叠加法：
Figure BDA0003559348020000075
p＝0，1，2，…2N，j＝0，1，2，…n；
其中，q是p/2N的余数，q＝0，1，2，…2n-1；
hjm(qΔt)由下式给出：
Figure BDA0003559348020000076
其中：
Figure BDA0003559348020000077
其中，m为风谱在频域内划分的分部数；Δt为时间增量；Δω为频率增量；S(ω)为互功率谱密度矩阵；Gjm为快速傅里叶变换中分解矩阵的元素；
Figure BDA0003559348020000078
为两个不同作用点之间的相位角；l为步长；
计算得到脉动风速时程图和自功率谱对比图；其中，自功率谱对比图反映了通过脉动风速时程图形成的功率谱与达文波特功率谱之间拟合效果的优劣；
输电杆塔每段单个节点的动态风荷载计算公式为：
ωs＝ωkμsA0，
其中，ωk为对应高度的脉动风压；μs为转角输电杆塔的风荷载体型系数，A0为迎风投影面积；
脉动风压ωk通过脉动风速由公式ωk＝0.5ρv2转换得到，ρ表示迎风面的方向系数；v表示脉动风速；
风荷载体型系数μs由公式：μs＝1.3(1+η)计算得到，其中η为背风面的风载荷降低系数，最后模拟出塔顶的脉动风载荷时程曲线图。
步骤d)中，对输电杆塔的动态响应方程进行时域分析求解方法为：
输电杆塔在稳态风的作用下，动态响应方程为：
Figure BDA0003559348020000081
其中，M为输电杆塔的整体质量矩阵，C为输电杆塔的整体阻尼矩阵，K为输电杆塔的整体刚度矩阵；
Figure BDA0003559348020000082
和u分别为节点的加速度、速度以及位移向量，F为输电杆塔线体系所受的动态风荷载；
将动态风载荷F沿0°风向角作用在塔顶上，对输电杆塔的动态响应方程进行时域分析求解。
实施例2
本实施例为一种基于脉动风影响下输电杆塔结构的动力学分析系统，其由有限元模型建立单元、模态分析单元、空间脉动风速模拟单元和时域分析求解单元组成。
有限元模型建立单元，根据某转角输电杆塔的工程实际数据建立输电杆塔有限元模型。
所述的有限元模型建立单元中，转角塔的工程实际数据包括塔高、呼高、档距和塔基尺寸，角钢与角钢之间用螺栓和连接板连接，忽略螺栓和连接板的几何尺寸效应，将角钢构件设置为刚接，采用Ansys有限元软件中Beaml88单元对输电杆塔的角钢构件进行建模。
模态分析单元，针对输电杆塔的结构特点，对输电杆塔有限元模型进行模态分析。
所述的模态分析单元中，所述的模态分析包括模态求解、扩展模态和结果观察，具体过程如下：
对输电杆塔有限元模型进行模态分析处理，模态分析的实质是计算输电杆塔结构振动特征方程的特征值和特征向量，根据结构动力学基本方程：
Figure BDA0003559348020000083
上式中：M为输电杆塔的整体质量矩阵；C为输电杆塔的整体阻尼矩阵；K为输电杆塔的整体刚度矩阵；
Figure BDA0003559348020000084
为结构加速度向量；
Figure BDA0003559348020000085
为结构速度向量；X为结构位移向量；F(t)为加载在结构上的载荷；t为时间；
由于工程结构中可不考虑阻尼对结构频率和模态的影响，忽略结构阻尼，在无外荷载的作用下输电杆塔结构自由振动动力学方程为：
Figure BDA0003559348020000091
当输电杆塔结构在固有频率下发生振动时，其输电杆塔结构速度和加速度向量表示为：
Figure BDA0003559348020000092
Figure BDA0003559348020000093
上式中，
Figure BDA0003559348020000094
为结构节点位移的节点振幅矩阵；ω为固有频率；
Figure BDA0003559348020000095
为结构发生振动的相位；
将上式(3)与(4)带入式(2)得：
Figure BDA0003559348020000096
按自由振动理论，r阶自由度系统的自由振动方程式有r个固有频率ωi，i＝1，2，3，…，n，并且由频率行列式决定，即：
|K-ω2M|＝0， (6)
当输电杆塔有限元模型的自由度为y时，式(6)具有y个特征值，记λj为第j个特征值，其平方根ωj是结构的第j阶固有频率，并且
Figure BDA0003559348020000097
是相应的第j阶特征向量，特征向量即为输电杆塔结构的模态振型，是输电杆塔结构在第j阶模态振型下的变形状态。
空间脉动风速模拟单元，结合快速傅里叶变换的谐波叠加法，通过仿真软件模拟输电杆塔的空间脉动风速。
所述的空间脉动风速模拟单元中，模拟输电杆塔空间脉动风速的方法为：
将风荷载视为静力风/平均风与动力风/脉动风的共同作用，在任意时刻高度Z处的风速v(z，t)表示为平均风速与脉动风速之和：
Figure BDA0003559348020000098
式中，
Figure BDA0003559348020000099
为Z高度处的平均风速，vf(z，t)为Z高度处的脉动风速；
平均风速与高度之间的变化规律称为风剖面，风剖面用指数律模拟：
Figure BDA00035593480200000910
式中，
Figure BDA00035593480200000911
为任一点平均风速，z为任一点的高度；
Figure BDA00035593480200000912
为10m处基准平均风速；α为地面粗糙度系数；
脉动风是一个均值为零的高斯随机过程，功率谱密度是对其均方值的量度，脉动风自功率谱的一种表达式Davenport谱为：
Figure BDA00035593480200000913
式中，
Figure BDA0003559348020000101
k为与地面粗糙度有关的系数，γ为脉动风频率，运用基于快速傅里叶变换技术的谐波叠加法：
Figure BDA0003559348020000102
p＝0，1，2，…2N，j＝0，1，2，…n；
其中，q是p/2N的余数，q＝0，1，2，…2n-1；
hjm(qΔt)由下式给出：
Figure BDA0003559348020000103
其中：
Figure BDA0003559348020000104
其中，m为风谱在频域内划分的分部数；Δt为时间增量；Δω为频率增量；S(ω)为互功率谱密度矩阵；Gjm为快速傅里叶变换中分解矩阵的元素；
Figure BDA0003559348020000105
为两个不同作用点之间的相位角；l为步长；
计算得到脉动风速时程图和自功率谱对比图；其中，自功率谱对比图反映了通过脉动风速时程图形成的功率谱与达文波特功率谱之间拟合效果的优劣。
输电杆塔的动态风荷载计算公式为：
ωs＝ωkμsA0，
其中，ωk为对应高度的脉动风压；μs为转角输电杆塔的风荷载体型系数，A0为迎风投影面积；
脉动风压ωk通过脉动风速由公式ωk＝0.5ρv2转换得到，ρ表示迎风面的方向系数；v表示脉动风速；
风荷载体型系数μs由公式：μs＝1.3(1+η)计算得到，其中，η为背风面的风载荷降低系数，最后模拟出塔顶的脉动风载荷时程曲线图。
时域分析求解单元，将动态风载荷沿0°风向角作用在塔顶节点上，对输电杆塔的动态响应方程进行时域分析求解，得到不同风向角下塔顶的位移响应。
所述的时域分析求解单元中，对输电杆塔的动态响应方程进行时域分析求解方法为：
输电杆塔在稳态风的作用下，动态响应方程为：
Figure BDA0003559348020000111
其中，M为输电杆塔的整体质量矩阵，C为输电杆塔的整体阻尼矩阵，K为输电杆塔的整体刚度矩阵；
Figure BDA0003559348020000112
和u分别为节点的加速度、速度以及位移向量，F为输电杆塔线体系所受的动态风荷载；
将动态风载荷F沿0°风向角作用在塔顶上，对输电杆塔的动态响应方程进行时域分析求解。
应用例
采用本发明的动力学分析方法进行应用，其步骤如下：
①对某500kV线路干字型转角塔的工程实际数据建立输电杆塔的有限元模型，塔高65m，呼高为40m，档距为220m，塔基尺寸为14.12m×14.12m，采用Ansys有限元软件中Beaml88单元对输电杆塔的角钢构件进行建模。
②针对建立好的高压铁塔模型，在ANSYS中对不同的材料施加属性，进入ANSYS求解器，设置分析类型为Modal，对输电杆塔采用Block Lanczos法进行模态分析，对杆塔塔腿施加固定约束，即约束其所有的自由度，得到输电杆塔的前6阶自振频率和振型。
③在Matlab中采用谐波合成法编写程序模拟输电杆塔的空间脉动风速，脉动风速谱选用Davenport谱，并且考虑到横向与竖向的空间相关性。脉动风速时间间隔为0.02s，共模拟200s。经过快速傅里叶变换以及谐波叠加法模拟得到塔顶的脉动风荷载时程曲线图以及自功率谱对比图，如图3所示。
④获得脉动风速后，因为作用在杆塔上的是风荷载，所以要将风速转换为风压，并且需要按照结构的外形特征考虑迎风面积以及风荷载系数。通过风荷载系数将动态风荷载按照等效静荷载来处理，并且要利用风压高度变化系数来描述风荷载随高度的变化。在模拟脉动风速时，已经考虑到风速的波动以及风速随高度的变化。
输电杆塔每段单个节点的动态风荷载计算公式为：
ωs＝ωkμsA0
其中，ωk为对应高度的风压；μs为转角输电杆塔的风荷载体型系数，A0为迎风投影面积。
脉动风压ωk可以通过脉动风速由公式ωk＝0.5ρv2转换得到。ρ表示迎风面的方向系数；v表示脉动风速；
风荷载体型系数μs可以由公式：μs＝1.3(1+η)计算得到，其中η为背风面的风载荷降低系数。通过计算公式对输电杆塔所承受的动态风荷载进行分析求解，如图4所示。
⑤输电杆塔在稳态风的作用下，动态响应方程为：
Figure BDA0003559348020000121
其中，M为输电杆塔的整体质量矩阵，C为输电杆塔的整体阻尼矩阵，K为输电杆塔的整体刚度矩阵；
Figure BDA0003559348020000122
和u分别为节点的加速度、速度以及位移向量，F为输电杆塔线体系所受的动态风荷载。
将动态风荷载F沿0°风向角作用在塔顶节点上。对转角输电杆塔的动态响应方程进行时域分析求解，得到位移时程曲线，如图5所示。
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包括在本发明的保护范围之内。

Claims (10)

1.基于脉动风影响下输电杆塔结构的动力学分析方法，其特征在于，包括步骤：
a)根据某转角输电杆塔的工程实际数据建立输电杆塔有限元模型；
b)针对输电杆塔的结构特点，对输电杆塔有限元模型进行模态分析；
c)结合快速傅里叶变换的谐波叠加法，通过仿真软件模拟输电杆塔的空间脉动风速；
d)将动态风载荷沿0°风向角作用在塔顶节点上，对输电杆塔的动态响应方程进行时域分析求解，得到不同风向角下塔顶的位移响应。
2.根据权利要求1所述的基于脉动风影响下输电杆塔结构的动力学分析方法，其特征在于，步骤a)中，转角输电杆塔的工程实际数据包括塔高、呼高、档距和塔基尺寸，角钢与角钢之间用螺栓和连接板连接，忽略螺栓和连接板的几何尺寸效应，将角钢构件设置为刚接，采用Ansys有限元软件中Beam188单元对输电杆塔的角钢构件进行建模。
3.根据权利要求1所述的基于脉动风影响下输电杆塔结构的动力学分析方法，其特征在于，步骤b)中，所述的模态分析包括模态求解、扩展模态和结果观察，具体过程如下：
对输电杆塔有限元模型进行模态分析处理，模态分析的实质是计算输电杆塔结构振动特征方程的特征值和特征向量，根据结构动力学基本方程：
Figure FDA0003559348010000011
上式中：M为输电杆塔的整体质量矩阵；C为输电杆塔的整体阻尼矩阵；K为输电杆塔的整体刚度矩阵；
Figure FDA0003559348010000012
为结构加速度向量；
Figure FDA0003559348010000013
为结构速度向量；X为结构位移向量；F(t)为加载在结构上的载荷；t为时间；
忽略结构阻尼，在无外荷载的作用下输电杆塔结构自由振动动力学方程为：
Figure FDA0003559348010000014
当输电杆塔结构在固有频率下发生振动时，其输电杆塔结构速度和加速度向量表示为：
Figure FDA0003559348010000015
Figure FDA0003559348010000016
上式中，
Figure FDA0003559348010000017
为结构节点位移的节点振幅矩阵；ω为固有频率；
Figure FDA0003559348010000018
为结构发生振动的相位；
将上式(3)与(4)带入式(2)得：
Figure FDA0003559348010000019
按自由振动理论，r阶自由度系统的自由振动方程式有r个固有频率ωi，i＝1,2,3,…,n，并且由频率行列式决定，即：
|K-ω2M|＝0， (6)
当输电杆塔有限元模型的自由度为y时，式(6)具有y个特征值，记λj为第j个特征值，其平方根ωj是结构的第j阶固有频率，并且
Figure FDA0003559348010000021
是相应的第j阶特征向量，特征向量即为输电杆塔结构的模态振型，是输电杆塔结构在第j阶模态振型下的变形状态。
4.根据权利要求1所述的基于脉动风影响下输电杆塔结构的动力学分析方法，其特征在于，步骤c)中，模拟输电杆塔空间脉动风速的方法为：
将风荷载视为静力风/平均风与动力风/脉动风的共同作用，在任意时刻高度Z处的风速v(z，t)表示为平均风速与脉动风速之和：
Figure FDA0003559348010000022
式中，
Figure FDA0003559348010000023
为Z高度处的平均风速，vf(z，t)为Z高度处的脉动风速；
平均风速与高度之间的变化规律称为风剖面，风剖面用指数律模拟：
Figure FDA0003559348010000024
式中，
Figure FDA0003559348010000025
为任一点平均风速，z为任一点的高度；
Figure FDA0003559348010000026
为10m处基准平均风速；α为地面粗糙度系数；
脉动风是一个均值为零的高斯随机过程，功率谱密度是对其均方值的量度，脉动风自功率谱的一种表达式Davenport谱为：
Figure FDA0003559348010000027
式中，
Figure FDA0003559348010000028
k为与地面粗糙度有关的系数，γ为脉动风频率，运用基于快速傅里叶变换技术的谐波叠加法：
Figure FDA0003559348010000029
p＝0，1，2，…2N，j＝0，1，2，…n；
其中，q是p/2N的余数，q＝0，1，2，…2n-1；
hjm(qΔt)由下式给出：
Figure FDA00035593480100000210
其中：
Figure FDA0003559348010000031
其中，m为风谱在频域内划分的分部数；Δt为时间增量；Δω为频率增量；S(ω)为互功率谱密度矩阵；Gjm为快速傅里叶变换中分解矩阵的元素；
Figure FDA0003559348010000032
为两个不同作用点之间的相位角；l为步长；
计算得到脉动风速时程图和自功率谱对比图；其中，自功率谱对比图反映通过脉动风速时程图形成的功率谱与达文波特功率谱之间拟合效果的优劣；
输电杆塔每段单个节点的动态风荷载计算公式为：
ωs＝ωkμsA0，
其中，ωk为对应高度的脉动风压；μs为转角输电杆塔的风荷载体型系数，A0为迎风投影面积；
脉动风压ωk通过脉动风速由公式ωk＝0.5ρv2转换得到，ρ表示迎风面的方向系数；v表示脉动风速；
风荷载体型系数μs由公式：μs＝1.3(1+η)计算得到，其中η为背风面的风载荷降低系数，最后模拟出塔顶的脉动风载荷时程曲线图。
5.根据权利要求1所述的基于脉动风影响下输电杆塔结构的动力学分析方法，其特征在于，步骤d)中，对输电杆塔的动态响应方程进行时域分析求解方法为：
输电杆塔在稳态风的作用下，动态响应方程为：
Figure FDA0003559348010000033
其中，M为输电杆塔的整体质量矩阵，C为输电杆塔的整体阻尼矩阵，K为输电杆塔的整体刚度矩阵；
Figure FDA0003559348010000034
和u分别为节点的加速度、速度以及位移向量，F为输电杆塔线体系所受的动态风荷载；
将动态风载荷F沿0°风向角作用在塔顶上，对输电杆塔的动态响应方程进行时域分析求解。
6.基于脉动风影响下输电杆塔结构的动力学分析系统，其特征在于，包括：
有限元模型建立单元，根据某转角输电杆塔的工程实际数据建立输电杆塔有限元模型；
模态分析单元，针对输电杆塔的结构特点，对输电杆塔有限元模型进行模态分析；
空间脉动风速模拟单元，结合快速傅里叶变换的谐波叠加法，通过仿真软件模拟输电杆塔的空间脉动风速；
时域分析求解单元，将动态风载荷沿0°风向角作用在塔顶节点上，对输电杆塔的动态响应方程进行时域分析求解，得到不同风向角下塔顶的位移响应。
7.根据权利要求6所述的基于脉动风影响下输电杆塔结构的动力学分析系统，其特征在于，所述的有限元模型建立单元中，转角输电杆塔的工程实际数据包括塔高、呼高、档距和塔基尺寸，角钢与角钢之间用螺栓和连接板连接，忽略螺栓和连接板的几何尺寸效应，将角钢构件设置为刚接，采用Ansys有限元软件中Beam188单元对输电杆塔的角钢构件进行建模。
8.根据权利要求6所述的基于脉动风影响下输电杆塔结构的动力学分析系统，其特征在于，所述的模态分析单元中，所述的模态分析包括模态求解、扩展模态和结果观察，具体过程如下：
对输电杆塔有限元模型进行模态分析处理，模态分析的实质是计算输电杆塔结构振动特征方程的特征值和特征向量，根据结构动力学基本方程：
Figure FDA0003559348010000041
上式中：M为输电杆塔的整体质量矩阵；C为输电杆塔的整体阻尼矩阵；K为输电杆塔的整体刚度矩阵；
Figure FDA0003559348010000042
为结构加速度向量；
Figure FDA0003559348010000043
为结构速度向量；X为结构位移向量；F(t)为加载在结构上的载荷；t为时间；
忽略结构阻尼，在无外荷载的作用下输电杆塔结构自由振动动力学方程为：
Figure FDA0003559348010000044
当输电杆塔结构在固有频率下发生振动时，其输电杆塔结构速度和加速度向量表示为：
Figure FDA0003559348010000045
Figure FDA0003559348010000046
上式中，
Figure FDA0003559348010000047
为结构节点位移的节点振幅矩阵；ω为固有频率；
Figure FDA0003559348010000048
为结构发生振动的相位；
将上式(3)与(4)带入式(2)得：
Figure FDA0003559348010000049
按自由振动理论，r阶自由度系统的自由振动方程式有r个固有频率ωi，i＝1，2，3，…，n，并且由频率行列式决定，即：
|K-ω2M|＝0， (6)
当输电杆塔有限元模型的自由度为y时，式(6)具有y个特征值，记λj为第j个特征值，其平方根ωj是结构的第j阶固有频率，并且
Figure FDA00035593480100000410
是相应的第j阶特征向量，特征向量即为输电杆塔结构的模态振型，是输电杆塔结构在第j阶模态振型下的变形状态。
9.根据权利要求6所述的基于脉动风影响下输电杆塔结构的动力学分析系统，其特征在于，所述的空间脉动风速模拟单元中，模拟输电杆塔空间脉动风速的方法为：
将风荷载视为静力风/平均风与动力风/脉动风的共同作用，在任意时刻高度Z处的风速v(z,t)表示为平均风速与脉动风速之和：
Figure FDA0003559348010000051
式中，
Figure FDA0003559348010000052
为Z高度处的平均风速，vf(z,t)为Z高度处的脉动风速；
平均风速与高度之间的变化规律称为风剖面，风剖面用指数律模拟：
Figure FDA0003559348010000053
式中，
Figure FDA0003559348010000054
为任一点平均风速，z为任一点的高度；
Figure FDA0003559348010000055
为10m处基准平均风速；α为地面粗糙度系数；
脉动风是一个均值为零的高斯随机过程，功率谱密度是对其均方值的量度，脉动风自功率谱的一种表达式Davenport谱为：
Figure FDA0003559348010000056
式中，
Figure FDA0003559348010000057
k为与地面粗糙度有关的系数，γ为脉动风频率，运用基于快速傅里叶变换技术的谐波叠加法：
Figure FDA0003559348010000058
p＝O，1，2，…2N，j＝0，1，2，…n；
其中，q是p/2N的余数，q＝0，1，2，…2n-1；
hjm(qΔt)由下式给出：
Figure FDA0003559348010000059
其中：
Figure FDA00035593480100000510
其中，m为风谱在频域内划分的分部数；Δt为时间增量；Δω为频率增量；S(ω)为互功率谱密度矩阵；Gjm为快速傅里叶变换中分解矩阵的元素；
Figure FDA0003559348010000061
为两个不同作用点之间的相位角；l为步长。
计算得到脉动风速时程图和自功率谱对比图；其中，自功率谱对比图反映了通过脉动风速时程图形成的功率谱与达文波特功率谱之间拟合效果的优劣；
输电杆塔每段单个节点的动态风荷载计算公式为：
ωs＝ωkμsA0，
其中，ωk为对应高度的脉动风压；μs为转角输电杆塔的风荷载体型系数，A0为迎风投影面积；
脉动风压ωk通过脉动风速由公式ωk＝0.5ρv2转换得到，ρ表示迎风面的方向系数；v表示脉动风速；
风荷载体型系数μs由公式：μs＝1.3(1+η)计算得到，其中η为背风面的风载荷降低系数，最后模拟出塔顶的脉动风载荷时程曲线图。
10.根据权利要求6所述的基于脉动风影响下输电杆塔结构的动力学分析系统，其特征在于，所述的时域分析求解单元中，对输电杆塔的动态响应方程进行时域分析求解方法为：
输电杆塔在稳态风的作用下，动态响应方程为：
Figure FDA0003559348010000062
其中，M为输电杆塔的整体质量矩阵，C为输电杆塔的整体阻尼矩阵，K为输电杆塔的整体刚度矩阵；
Figure FDA0003559348010000063
和u分别为节点的加速度、速度以及位移向量，F为输电杆塔线体系所受的动态风荷载；
将动态风载荷F沿0°风向角作用在塔顶上，对输电杆塔的动态响应方程进行时域分析求解。
CN202210288857.2A 2022-03-22 2022-03-22 基于脉动风影响下输电杆塔结构的动力学分析方法及系统 Pending CN114611365A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210288857.2A CN114611365A (zh) 2022-03-22 2022-03-22 基于脉动风影响下输电杆塔结构的动力学分析方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210288857.2A CN114611365A (zh) 2022-03-22 2022-03-22 基于脉动风影响下输电杆塔结构的动力学分析方法及系统

Publications (1)

Publication Number Publication Date
CN114611365A true CN114611365A (zh) 2022-06-10

Family

ID=81864370

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210288857.2A Pending CN114611365A (zh) 2022-03-22 2022-03-22 基于脉动风影响下输电杆塔结构的动力学分析方法及系统

Country Status (1)

Country Link
CN (1) CN114611365A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115062432A (zh) * 2022-06-29 2022-09-16 国网甘肃省电力公司电力科学研究院 一种基于数据驱动的输电杆塔倾斜预警方法及系统
CN115144173A (zh) * 2022-09-05 2022-10-04 济南百顿机械设备有限公司 一种基于智能压电阻尼器的输电塔监测方法及设备
CN116108725A (zh) * 2023-04-10 2023-05-12 国网浙江省电力有限公司温州供电公司 输电线路铁塔中预应力拉索的参数优化方法、装置及介质

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115062432A (zh) * 2022-06-29 2022-09-16 国网甘肃省电力公司电力科学研究院 一种基于数据驱动的输电杆塔倾斜预警方法及系统
CN115062432B (zh) * 2022-06-29 2023-09-08 国网甘肃省电力公司电力科学研究院 一种基于数据驱动的输电杆塔倾斜预警方法及系统
CN115144173A (zh) * 2022-09-05 2022-10-04 济南百顿机械设备有限公司 一种基于智能压电阻尼器的输电塔监测方法及设备
CN115144173B (zh) * 2022-09-05 2022-11-25 济南百顿机械设备有限公司 一种基于智能压电阻尼器的输电塔监测方法及设备
CN116108725A (zh) * 2023-04-10 2023-05-12 国网浙江省电力有限公司温州供电公司 输电线路铁塔中预应力拉索的参数优化方法、装置及介质

Similar Documents

Publication Publication Date Title
CN114611365A (zh) 基于脉动风影响下输电杆塔结构的动力学分析方法及系统
Zhao et al. Studies on application of scissor-jack braced viscous damper system in wind turbines under seismic and wind loads
Chen et al. Parametric study on the along-wind response of the CAARC building to downbursts in the time domain
Selvam Computation of pressures on Texas Tech University building using large eddy simulation
Huang et al. Experimental and computational simulation for wind effects on the Zhoushan transmission towers
CN105354377A (zh) 一种确定输电杆塔的脉动风振荷载的方法
CN106934185A (zh) 一种弹性介质的流固耦合多尺度流动模拟方法
Liang et al. Constructing simplified models for dynamic analysis of monopile-supported offshore wind turbines
Zhang et al. Modal parameters of a transmission tower considering the coupling effects between the tower and lines
Li et al. Frequency domain dynamic analyses of freestanding bridge pylon under wind and waves using a copula model
CN109543258A (zh) 一种基于模态叠加法的风电机组塔架应力状态计算方法
CN112507414B (zh) 评估下击暴流下输电铁塔安全性的方法、系统及存储介质
Fan et al. Collapse analysis of wind turbine tower under the coupled effects of wind and near‐field earthquake
Karmakar et al. Conditional simulation of non-Gaussian wind velocity profiles: Application to buffeting response of Vincent Thomas suspension bridge
Qu et al. Comparative study of short-term extreme responses and fatigue damages of a floating wind turbine using two different blade models
Liu et al. Proposal of a novel GPU-accelerated lifetime optimization method for onshore wind turbine dampers under real wind distribution
Vemuri et al. Wind Analysis of Tall Buildings Using Codal Provisions, Stochastic Approach and CFD–A Comparative Study
Gu et al. Equivalent static wind loads for stability design of large span roof structures
Frandsen Computational fluid structure interaction applied to long-span bridge design
CN109340060A (zh) 一种基于模态叠加法的风电机组塔架振动状态计算方法
Duthinh et al. The use of wind tunnel measurements in building design
Cheng et al. Wind resource assessment at mountainous wind farm: Fusion of RANS and vertical multi-point on-site measured wind field data
Mendes et al. Horizontal and vertical axis wind turbines on existing jacket platforms: Part 1–A comparative study
Hu et al. Vortex-induced force model and reliability analysis of a steel box girder with projecting slab in a cable-stayed bridge
Hua et al. Turbulence integral scale corrections to experimental results of aeroelastic models with large geometric scales: Application to gust loading factor of a transmission line tower

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