CN110851905A - 一种峡谷微地形条件下的输电铁塔风振系数计算方法 - Google Patents

一种峡谷微地形条件下的输电铁塔风振系数计算方法 Download PDF

Info

Publication number
CN110851905A
CN110851905A CN201911110495.2A CN201911110495A CN110851905A CN 110851905 A CN110851905 A CN 110851905A CN 201911110495 A CN201911110495 A CN 201911110495A CN 110851905 A CN110851905 A CN 110851905A
Authority
CN
China
Prior art keywords
wind
iron tower
tower
calculating
terrain
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.)
Granted
Application number
CN201911110495.2A
Other languages
English (en)
Other versions
CN110851905B (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.)
Fuzhou University
Original Assignee
Fuzhou University
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 Fuzhou University filed Critical Fuzhou University
Priority to CN201911110495.2A priority Critical patent/CN110851905B/zh
Publication of CN110851905A publication Critical patent/CN110851905A/zh
Application granted granted Critical
Publication of CN110851905B publication Critical patent/CN110851905B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Wind Motors (AREA)

Abstract

本发明涉及一种峡谷微地形条件下的输电铁塔风振系数计算方法,包括以下步骤:步骤S1:利用Google地图提取待测峡谷一定尺度的微地形,并建立网格化3Dmax地形模型;步骤S2:进行流速仿真,确定风场CFD仿真风剖面幂指数;步骤S3:构建杆梁铁塔模型,并进行模态分析,得到铁塔前三阶振型及其模态频率;步骤S4:进行脉动风速时程模拟,得到风速时程数据;步骤S5:进行铁塔风荷载计算,计算得到风荷载时程数据;步骤S6:对进行瞬态动力分析,得到铁塔各塔段位移响应均方差;步骤S7:根据得到的各塔段位移响应均方差,计算得到铁塔风振系数,并将得到的铁塔风振系数用于铁塔抗风设计。本发明实现了输电工程沿线峡谷微地形处铁塔风振系数的高精度计算。

Description

一种峡谷微地形条件下的输电铁塔风振系数计算方法
技术领域
本发明涉及输电铁塔抗风设计领域,具体涉及一种峡谷微地形条件下的输电铁塔风振系数计算方法。
背景技术
现有的风工程研究方法有理论分析、现场监测、风洞实验和数值模拟四大类。有人通过自主开发的新型加速度传感器,对强风作用下能盘线输电塔进行了实时在线监测,所得数据接近实际值,但增加了工程造价,不利于推广。有人通过风洞试验揭示了塔线耦合的机理,较为完整分析了整体系统的风振响应,所得结果与事故记录较为一致,但实验条件较为严格,无法通用于实际工程。有人采用有限元建模的方法,通过数值模拟得到不同高度层风荷载并施加到铁塔模型上,分析得到塔身和横担处的位移时程、轴力时程和风振系数,并与相关规范值进行比较,能够满足工程需要。一般来说,输电铁塔的抗风设计还要考虑线路经过复杂微地形时的风场变化问题,微地形影响风场分布,如越山风和峽谷风等。有人给出了输电线路塔架位置周围的一般微地形风场特征,得到了塔架位置风速加速度随高度和风向角的变化规律,但其所建微地形模型较为简单,所得数据精度不够。我国为多山国家,与平地地形不同,山区微地形对近地风场必然造成很大的影响,且山地环境中一般没有气象监测,难以获得风场特性,就无法准确分析输电铁塔的风振响应特性,增加了铁塔风振系数计算的难度,这对铁塔在微地形风场下的安全性提出了新的要求。
发明内容
有鉴于此,本发明的目的在于提供一种峡谷微地形条件下的输电铁塔风振系数计算方法以解决现在国内输电工程中难以确定峡谷地形处铁塔风振系数的问题,以提升铁塔在峡谷微地形风场下的安全性能。
为实现上述目的,本发明采用如下技术方案:
一种峡谷微地形条件下的输电铁塔风振系数计算方法,包括以下步骤:
步骤S1:利用Google地图提取待测峡谷一定尺度的微地形,并建立网格化3Dmax地形模型;
步骤S2:将得到的3Dmax地形模型导入CFD环境中进行流速仿真,确定风场CFD仿真风剖面幂指数;
步骤S3:构建杆梁铁塔模型,并进行模态分析,得到铁塔前三阶振型及其模态频率;
步骤S4:根据杆梁铁塔模型,进行脉动风速时程模拟,得到风速时程数据;
步骤S5:根据风速时程数据,进行铁塔风荷载计算,计算得到风荷载时程数据;
步骤S6:根据杆梁铁塔模型合风荷载时程数据,对进行瞬态动力分析,得到铁塔各塔段位移响应均方差;
步骤S7:根据得到的各塔段位移响应均方差,计算得到铁塔风振系数,并将得到的铁塔风振系数用于铁塔抗风设计。
进一步的,所述步骤S2具体为:
步骤S21:将得到的3Dmax地形模型导入CFD,中,并以该地形为底面,构造一个立方体,并划分网格;
步骤S22:根据近地气象站实测资料统计出特定年限10min最大平均风速,作为仿真10m高度入口风速基准风速,并确定风场CFD仿真风剖面幂指数。
进一步的,所述步骤S4中将铁塔沿高度简化为分段加载模型,采用线性滤波法中的自回归模型数值方法模拟脉动风速时程,具体包括如下:
步骤S41:采用davenport谱为校验谱:
Figure BDA0002272570640000031
其中,K为地面粗糙度系数;n为频率;
Figure BDA0002272570640000033
为离地10m高度处的基准风速;
步骤S42:计算脉动风协方差矩阵R:
其中,为空间两点脉动风互功率谱;
步骤S43:计算自回归系数矩阵ψk:
Figure BDA0002272570640000036
其中,ψk为M*M阶矩阵,
Figure BDA0002272570640000037
为pM*M阶矩阵,p为模型阶数;
步骤S44:求解给定方差的随机过程N(t):
Figure BDA0002272570640000038
其中,Δt为时间步长,由上述步骤得出的R,R0,结合上式可以求出N(t);
步骤S45:求解4维脉动风时程;平均风与脉动风叠加获得总风速。
进一步的,所述步骤S5具体为:
步骤S51:在已知风速vi条件下,由
Figure BDA0002272570640000041
可得到铁塔第i塔段的风荷载时程
式中,vi(t)为模拟得到的铁塔第i塔段中心风速时程;μsi、Ai分别为第i塔段风荷载体型系数和挡风面积
步骤S52:通过模拟出的风速时程数据,得到风荷载时程数据。
进一步的,所述铁塔瞬态动力分析具体为:根据得到的风荷载时程数据,在ANSYS软件环境中加载到铁塔有限元模型相应节点上,设置铁塔材料阻尼比、荷载步和荷载子步,采用完全法对铁塔进行瞬态动力分析,统计铁塔各塔段位移响应均方差。
进一步的,所述铁塔风振系数计算具体为:
在ANSYS环境中对铁塔进行风振响应时程分析,根据动力时程计算结果,采用随机振动理论对铁塔进行相应计算,
铁塔z高度处风振系数βi为:
Figure BDA0002272570640000042
式中,i为铁塔分段,i=1,2,…,n,Mi为铁塔分段质量;ω0为基本风压,ω0=v2/1600,v为设计风速;σli、μzi分别为铁塔z高度处塔段位移均方差、风压变化系数;g为峰值因子。
本发明与现有技术相比具有以下有益效果:
本发明实现了输电工程沿线峡谷微地形处铁塔风振系数的高精度计算,通过常规地形和微地形条件下不同风振系数计算值的对比,可以更为直观的显现出微地形对铁塔风振系数的影响,对电力铁塔的抗风设计提供了精确的数据。
附图说明
图1是本发明方法流程图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
请参照图1,本发明提供一种峡谷微地形条件下的输电铁塔风振系数计算方法利用Google地图提取一定尺度的微地形,建立网格化3Dmax地形模型,导入CFD仿真软件中,以收集近地气象监测风场数据为入口风速,可以获得经过峡谷微地形后到达铁塔位置的风场特性,进而获得总风速样本;然后在ANSYS软件环境中建立铁塔杆梁混合模型,通过风速样本计算铁塔风荷载,进而针对铁塔进行风振响应分析,最终通过随机振动理论获得准确的风振系数。具体包括如下步骤:
步骤S1:利用Google地图提取待测峡谷一定尺度的微地形,并建立网格化3Dmax地形模型;选取待测峡谷,气象监测点一般在近地面分布设定,因此需要提取一定尺度的该峡谷地形,利用近地可获取的气象风场数据,为峡谷输入风速,进行风场经过峡谷微地形后到达铁塔位置的变化风场研究。从google地图中提取铁塔所在方圆峡谷地形,导入3Dmax中显示具体模型并进行网格划分,最终把3Dmax图导入计算流体仿真软件(CFD)环境中进行流速仿真
步骤S2:将得到的3Dmax地形模型导入CFD环境中进行流速仿真,以该地形为底面,构造一个立方体,并划分网格。根据近地气象站历年实测资料统计出10min最大平均风速,为仿真10m高度入口风速基准风速。并确定风场CFD仿真风剖面幂指数为α=0.12。
步骤S3:采用AutoCAD铁塔组装图建立铁塔有限元杆梁混合模型,横隔梁及塔身主材采用梁单元,其余杆件采用杆单元。在ANSYS中对铁塔混合模型进行模态分析,得到铁塔前三阶振型及其模态频率。
步骤S4:将铁塔沿高度简化为分段加载模型,采用线性滤波法中的自回归模型数值方法模拟脉动风速时程,主要分为以下几个阶段:
采用davenport谱为校验谱:
Figure BDA0002272570640000061
计算脉动风协方差矩阵R:
Figure BDA0002272570640000062
计算自回归系数矩阵ψk
求解给定方差的随机过程
N(t):
求解4维脉动风时程;平均风与脉动风叠加获得总风速。
步骤S5:在已知风速vi条件下,由可得到铁塔第i塔段的风荷载时程。式中,vi(t)为模拟得到的铁塔第i塔段中心风速时程;μsi、Ai分别为第i塔段风荷载体型系数和挡风面积,通过模拟出的风速时程数据即可得到风荷载时程数据。
步骤S6:在已知铁塔各塔段风荷载时程后,在ANSYS软件环境中加载到铁塔有限元模型相应节点上,正确设置铁塔材料阻尼比,设置荷载步和荷载子步,采用完全法对铁塔进行瞬态动力分析,统计铁塔各塔段位移响应均方差。
步骤S7:根据得到的各塔段位移响应均方差,计算得到铁塔风振系数,并将得到的铁塔风振系数用于铁塔抗风设计。在ANSYS环境中对铁塔进行风振响应时程分析,根据动力时程计算结果,采用随机振动理论对铁塔进行相应计算。铁塔z高度处风振系数βi(i为铁塔分段,i=1,2,…,n)为:
Figure BDA0002272570640000071
式中,Mi为铁塔分段质量;ω0为基本风压,ω0=v2/1600,v为设计风速;σli、μzi分别为铁塔z高度处塔段位移均方差、风压变化系数;g为峰值因子。
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。

Claims (6)

1.一种峡谷微地形条件下的输电铁塔风振系数计算方法,其特征在于,包括以下步骤:
步骤S1:利用Google地图提取待测峡谷一定尺度的微地形,并建立网格化3Dmax地形模型;
步骤S2:将得到的3Dmax地形模型导入CFD环境中进行流速仿真,确定风场CFD仿真风剖面幂指数;
步骤S3:构建杆梁铁塔模型,并进行模态分析,得到铁塔前三阶振型及其模态频率;
步骤S4:根据杆梁铁塔模型,进行脉动风速时程模拟,得到风速时程数据;
步骤S5:根据风速时程数据,进行铁塔风荷载计算,计算得到风荷载时程数据;
步骤S6:根据杆梁铁塔模型合风荷载时程数据,对进行瞬态动力分析,得到铁塔各塔段位移响应均方差;
步骤S7:根据得到的各塔段位移响应均方差,计算得到铁塔风振系数,并将得到的铁塔风振系数用于铁塔抗风设计。
2.根据权利要求1所述的一种峡谷微地形条件下的输电铁塔风振系数计算方法,其特征在于,所述步骤S2具体为:
步骤S21:将得到的3Dmax地形模型导入CFD,中,并以该地形为底面,构造一个立方体,并划分网格;
步骤S22:根据近地气象站实测资料统计出特定年限10min最大平均风速,作为仿真10m高度入口风速基准风速,并确定风场CFD仿真风剖面幂指数。
3.根据权利要求1所述的一种峡谷微地形条件下的输电铁塔风振系数计算方法,其特征在于,所述步骤S4中将铁塔沿高度简化为分段加载模型,采用线性滤波法中的自回归模型数值方法模拟脉动风速时程,具体包括如下:
步骤S41:采用davenport谱为校验谱:
Figure FDA0002272570630000021
其中,K为地面粗糙度系数;
Figure FDA0002272570630000022
n为频率;为离地10m高度处的基准风速;
步骤S42:计算脉动风协方差矩阵R:
Figure FDA0002272570630000024
其中,为空间两点脉动风互功率谱;
步骤S43:计算自回归系数矩阵ψk:
Figure FDA0002272570630000026
其中,ψk为M*M阶矩阵,
Figure FDA0002272570630000027
为pM*M阶矩阵,p为模型阶数;
步骤S44:求解给定方差的随机过程N(t):
Figure FDA0002272570630000028
其中,Δt为时间步长,由上述步骤得出的R,R0,结合上式可以求出N(t);
步骤S45:求解4维脉动风时程;平均风与脉动风叠加获得总风速。
4.根据权利要求1所述的一种峡谷微地形条件下的输电铁塔风振系数计算方法,其特征在于,所述步骤S5具体为:
步骤S51:在已知风速vi条件下,由
Figure FDA0002272570630000031
可得到铁塔第i塔段的风荷载时程
式中,vi(t)为模拟得到的铁塔第i塔段中心风速时程;μsi、Ai分别为第i塔段风荷载体型系数和挡风面积
步骤S52:通过模拟出的风速时程数据,得到风荷载时程数据。
5.根据权利要求1所述的一种峡谷微地形条件下的输电铁塔风振系数计算方法,其特征在于,所述铁塔瞬态动力分析具体为:根据得到的风荷载时程数据,在ANSYS软件环境中加载到铁塔有限元模型相应节点上,设置铁塔材料阻尼比、荷载步和荷载子步,采用完全法对铁塔进行瞬态动力分析,统计铁塔各塔段位移响应均方差。
6.根据权利要求1所述的一种峡谷微地形条件下的输电铁塔风振系数计算方法,其特征在于,所述铁塔风振系数计算具体为:
在ANSYS环境中对铁塔进行风振响应时程分析,根据动力时程计算结果,采用随机振动理论对铁塔进行相应计算,
铁塔z高度处风振系数βi为:
Figure FDA0002272570630000032
式中,i为铁塔分段,i=1,2,…,n,Mi为铁塔分段质量;ω0为基本风压,ω0=v2/1600,v为设计风速;σli、μzi分别为铁塔z高度处塔段位移均方差、风压变化系数;g为峰值因子。
CN201911110495.2A 2019-11-14 2019-11-14 一种峡谷微地形条件下的输电铁塔风振系数计算方法 Active CN110851905B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911110495.2A CN110851905B (zh) 2019-11-14 2019-11-14 一种峡谷微地形条件下的输电铁塔风振系数计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911110495.2A CN110851905B (zh) 2019-11-14 2019-11-14 一种峡谷微地形条件下的输电铁塔风振系数计算方法

Publications (2)

Publication Number Publication Date
CN110851905A true CN110851905A (zh) 2020-02-28
CN110851905B CN110851905B (zh) 2022-06-14

Family

ID=69601622

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911110495.2A Active CN110851905B (zh) 2019-11-14 2019-11-14 一种峡谷微地形条件下的输电铁塔风振系数计算方法

Country Status (1)

Country Link
CN (1) CN110851905B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111651804A (zh) * 2020-03-31 2020-09-11 重庆科技学院 基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法
CN111985019A (zh) * 2020-03-31 2020-11-24 重庆科技学院 带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法
CN112580231A (zh) * 2020-11-24 2021-03-30 广东电网有限责任公司阳江供电局 一种输电杆塔在强风下力学特性仿真方法
CN113358313A (zh) * 2021-06-08 2021-09-07 国网河南省电力公司电力科学研究院 一种输电铁塔螺栓松动试验方法
CN114048777A (zh) * 2021-11-11 2022-02-15 南方电网数字电网研究院有限公司 输电线路铁塔振动分析方法、装置和计算机设备
CN114741760A (zh) * 2022-04-15 2022-07-12 国网陕西省电力有限公司经济技术研究院 一种概率密度可调的风速场数值模拟方法及系统
CN118410744A (zh) * 2024-07-01 2024-07-30 国网浙江省电力有限公司经济技术研究院 一种输电线路导线风载荷监测方法、装置和设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010072060A1 (zh) * 2008-12-23 2010-07-01 徐州中煤钢结构建设有限公司 超大跨度双层网架穹顶施工工艺
CN105354377A (zh) * 2015-10-29 2016-02-24 中国电力科学研究院 一种确定输电杆塔的脉动风振荷载的方法
CN107291985A (zh) * 2017-05-23 2017-10-24 南京航空航天大学 一种冷却塔施工全过程风振系数取值方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010072060A1 (zh) * 2008-12-23 2010-07-01 徐州中煤钢结构建设有限公司 超大跨度双层网架穹顶施工工艺
CN105354377A (zh) * 2015-10-29 2016-02-24 中国电力科学研究院 一种确定输电杆塔的脉动风振荷载的方法
CN107291985A (zh) * 2017-05-23 2017-10-24 南京航空航天大学 一种冷却塔施工全过程风振系数取值方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HONGJIE ZHANG ET AL: "Impact Analysis of Wind Direction upon Wind Load on Tower Line in", 《IOP CONFERENCE SERIES: EARTH AND ENVIRONMENTAL SCIENCE》 *
楼文娟 等: "山地风场特性及其对输电线路风偏响应的影响", 《土木工程学报》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111651804A (zh) * 2020-03-31 2020-09-11 重庆科技学院 基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法
CN111985019A (zh) * 2020-03-31 2020-11-24 重庆科技学院 带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法
CN111985019B (zh) * 2020-03-31 2022-05-20 重庆科技学院 带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法
CN111651804B (zh) * 2020-03-31 2022-08-26 重庆科技学院 考虑塔线耦合影响的常规输电线路设计风载荷计算方法
CN112580231A (zh) * 2020-11-24 2021-03-30 广东电网有限责任公司阳江供电局 一种输电杆塔在强风下力学特性仿真方法
CN112580231B (zh) * 2020-11-24 2024-05-28 广东电网有限责任公司阳江供电局 一种输电杆塔在强风下力学特性仿真方法
CN113358313A (zh) * 2021-06-08 2021-09-07 国网河南省电力公司电力科学研究院 一种输电铁塔螺栓松动试验方法
CN114048777A (zh) * 2021-11-11 2022-02-15 南方电网数字电网研究院有限公司 输电线路铁塔振动分析方法、装置和计算机设备
CN114741760A (zh) * 2022-04-15 2022-07-12 国网陕西省电力有限公司经济技术研究院 一种概率密度可调的风速场数值模拟方法及系统
CN118410744A (zh) * 2024-07-01 2024-07-30 国网浙江省电力有限公司经济技术研究院 一种输电线路导线风载荷监测方法、装置和设备

Also Published As

Publication number Publication date
CN110851905B (zh) 2022-06-14

Similar Documents

Publication Publication Date Title
CN110851905B (zh) 一种峡谷微地形条件下的输电铁塔风振系数计算方法
CN111651885B (zh) 一种智慧型海绵城市洪涝预报方法
CN102289570B (zh) 基于降雨-径流-洪水演进计算的洪水预报方法
CN116596321B (zh) 基于神经网络的建筑风险监测与评估方法和系统
Greenway An analytical approach to wind velocity gust factors
CN109543356A (zh) 考虑空间非平稳性的海洋内部温盐结构遥感反演方法
CN111144656A (zh) 一种基于gis的灾害评价分析方法
CN116434475A (zh) 一种基于降雨强度与地表位移变化率的滑坡预警方法
Zhao et al. Radial pressure profile of typhoon field near ground surface observed by distributed meteorologic stations
CN116823172B (zh) 基于模型优化的工程造价评估方法及系统
CN115471679A (zh) 一种天然河道水位流量同步同化方法及智能系统
CN112100922A (zh) 一种基于wrf和cnn卷积神经网络的风资源预测方法
Tarbell et al. An example of the initialization of the divergent wind component in a mesoscale numerical weather prediction model
CN116738553B (zh) 抗震结构的不确定参数敏感性分析方法
CN117094167A (zh) 滑坡全过程表面三维形变北斗雷达面场景数据仿真方法
CN113486295A (zh) 基于傅里叶级数的臭氧总量变化预测方法
CN112883339B (zh) 一种地震有感范围确定方法及系统
CN112697386B (zh) 一种风载荷作用下山区输电塔稳定性监控方法与系统
CN106338254B (zh) 基于3d激光扫描的地下工程施工快速监测预测系统及方法
CN109933921B (zh) 一种滚石灾害风险评估方法、装置、系统及存储介质
CN208887684U (zh) 一种具有mems传感器的应力及应变测试仪
Anning et al. Effects of the physical process ensemble technique on simulation of summer precipitation over China
CN206627073U (zh) 基于3d激光扫描的地下工程施工快速监测预测装置
CN116595898B (zh) 平原河道桥群阻水叠加影响量化分析方法及系统
CN114861566B (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