CN112507414B - 评估下击暴流下输电铁塔安全性的方法、系统及存储介质 - Google Patents
评估下击暴流下输电铁塔安全性的方法、系统及存储介质 Download PDFInfo
- Publication number
- CN112507414B CN112507414B CN202011210259.0A CN202011210259A CN112507414B CN 112507414 B CN112507414 B CN 112507414B CN 202011210259 A CN202011210259 A CN 202011210259A CN 112507414 B CN112507414 B CN 112507414B
- Authority
- CN
- China
- Prior art keywords
- wind
- power transmission
- load
- tower
- downburst
- 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 204
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000003860 storage Methods 0.000 title claims abstract description 10
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims abstract description 192
- 229910052742 iron Inorganic materials 0.000 claims abstract description 96
- 238000004088 simulation Methods 0.000 claims abstract description 89
- 238000004458 analytical method Methods 0.000 claims abstract description 71
- 238000013461 design Methods 0.000 claims abstract description 17
- 230000003068 static effect Effects 0.000 claims abstract description 13
- 230000001052 transient effect Effects 0.000 claims abstract description 10
- 238000006073 displacement reaction Methods 0.000 claims description 53
- 230000009471 action Effects 0.000 claims description 29
- 230000008569 process Effects 0.000 claims description 21
- 230000006870 function Effects 0.000 claims description 20
- 239000000463 material Substances 0.000 claims description 19
- 230000004044 response Effects 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000011156 evaluation Methods 0.000 claims description 14
- 238000009826 distribution Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 210000001503 joint Anatomy 0.000 claims description 4
- 238000004519 manufacturing process Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 16
- 238000004590 computer program Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 4
- 229910000831 Steel Inorganic materials 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- 239000012212 insulator Substances 0.000 description 2
- 239000000725 suspension Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000009439 industrial construction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000011160 research Methods 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/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- 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
- G06Q—INFORMATION 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/82—Energy audits or management systems therefor
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Strategic Management (AREA)
- Development Economics (AREA)
- Mathematical Analysis (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Water Supply & Treatment (AREA)
- Operations Research (AREA)
- Structural Engineering (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Civil Engineering (AREA)
- Architecture (AREA)
- Primary Health Care (AREA)
- Game Theory and Decision Science (AREA)
- Computational Mathematics (AREA)
- Quality & Reliability (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种评估下击暴流下输电铁塔安全性的方法、系统及存储介质,方法包括以下步骤:根据输电铁塔设计信息,划分输电铁塔塔段,设置风荷载模拟点,建立输电铁塔三维有限元数值模型,在模型文件中预设载荷信息文件接口;确定下击暴流风场结构,完成对风场结构各组成成分的模拟,遍历风荷载模拟点,生成各风荷载模拟点处的风速以及风荷载;将载荷导入输电铁塔三维数值模型中,对铁塔进行瞬态动力分析,根据分析结果计算输电铁塔结构整体风振系数;根据天气预报数据,生成下击暴流参数,计算各点时变平均风,用结构整体风振系数将平均风放大,得到等效静力风荷载,对铁塔进行准静态分析,利用分析结果评估输电铁塔结构安全性。
Description
技术领域
本发明属于输电技术领域,具体涉及一种输电铁塔结构的安全评估方法。具体是下击暴流天气下输电铁塔力学响应分析以及快速评估安全状态的方法、系统及存储介质。
背景技术
近年来,我国的科学技术和国民经济迅猛发展,国家对工业建设投入了大量人力物力,而电力工业更是其中的重要环节,倍受国家关注。输电铁塔结构是由不同规格角钢连接而成的空间体系,不同输电铁塔的高度也不尽相同,从几十米到几百米不等。输电铁塔体系结构复杂,对它的研究也较一般的建筑结构难度更大,且现今的输电线路规模增大,研究输电铁塔结构的安全性,特别是在极端天气下的健康状况是非常有必要的。
下击暴流作为极端天气中的一种,能产生破坏力极强的瞬时大风,对于输电铁塔这种高耸柔性结构,极易受其影响,会严重影响输电线路的安全运行。输电线路作为国家的大型生命线工程,一旦发生倒塌破坏事故,会导致国家经济受损,人民生命财产安全受到威胁。而工程中对下击暴流风场结构的确定方法不够具体,且缺乏较为简洁快速的铁塔抗下击暴流的设计的安全评估方法。因此,有必要发明一种有效的方法能够从理论上确定下击暴流的风场结构,从技术上快速解决输电铁塔结构的力学响应分析以及安全状态评估的问题,从而提高铁塔抗强对流天气的快速应急管理能力。
发明内容
发明目的:本发明目的在于能够提供一种评估下击暴流下输电铁塔安全性的方法,利用该方法能够快速分析在下击暴流这种强对流天气下输电铁塔结构的风振响应,实现对输电铁塔结构的快速安全分析,满足实际的工程应用。
技术方案:一种评估下击暴流下输电铁塔安全性的方法,包括以下步骤:
1)根据输电铁塔设计信息,设置风荷载模拟点,建立输电铁塔三维有限元数值模型;
2)确定下击暴流风场结构,完成对风场结构各组成成分的模拟,遍历风荷载模拟点,生成各风荷载模拟点处的风速以及风荷载;
3)将载荷信息文件对接至输电铁塔三维有限元数值模型中,对铁塔进行瞬态动力分析,根据分析结果计算输电铁塔结构整体风振系数;
4)根据天气预报数据,生成下击暴流参数,并对铁塔进行准静态分析;
5)利用准静态分析结果对输电铁塔结构进行快速安全评估。
进一步地,步骤1)包括以下具体步骤:
11)根据实际输电线路中铁塔的设计信息,将输电铁塔划分为不同的塔段以及不同的分块区域,把每个分块区域简化为一个模拟点进行风荷载模拟,以模拟点的风速代替所在分块区域所有节点的风速,所述输电铁塔设计信息包括输电铁塔构件布局、结构尺寸、材质、牌号、数量等信息;
12)根据实际输电铁塔的材料设计信息,利用有限元软件对输电铁塔进行三维有限元数值建模,采用梁单元对输电铁塔杆件进行模拟,单元之间的连接方式为固接;对梁单元进行设置,分别设置铁塔杆件的参数,包括杨氏模量、泊松比、密度、截面形状等;
13)将建立完成的输电铁塔三维有限元数值模型保存成为模型文件,并在文件中预留接口与载荷信息文件对接。
进一步地,步骤2)包括以下具体步骤:
式中:V(z)为竖向风剖面平均风速;f(t)为最大值为1的时间函数,描述输电铁塔处竖向风速剖面随时间的变化;Vc(t)是输电铁塔处合成水平风速;Vt为风暴中心移动速度,Vr(r,t)为径向风剖面风速,θ为输电线路与径向风剖面风速的夹角,α是合成水平风速与输电线路的夹角,r为输电铁塔至风暴中心距离;
b)竖向风剖面V(z)为:
式中,z为风载荷模拟点处的高度,Vmax为下击暴流的最大水平风速,Zmax为最大水平风速的产生高度;
径向风剖面Vr(r,t)为:
式中:Vr,max最大径向风速,rmax径向风速最大出距雷暴中心的距离,r观测点距雷暴中心的距离,Rr雷暴径向特征尺寸,T为时间常数;
c)脉动风速可表达为:
d)将各模拟点的时变平均风与脉动风相结合,计算出各点的下击暴流合成风速,并通过下式确定输电铁塔上各模拟点的下击暴流风荷载:
式中,Fx(t)为模拟点x方向的风荷载,N;Fy(t)模拟点y方向的风荷载,N;V(t)为下击暴流风速,m/s;ρ为空气密度,取1.226kg/m3;Cx、Cy分别为x、y方向阻力系数,角钢取1.3,导线、地线、悬索、绝缘子以及拉线取1.0;Ax和Ay分别为x、y方向输电铁塔各节点所辖的迎风面积。
进一步地,步骤3)包括以下具体步骤:
A)将各模拟点的风荷载信息存储于载荷信息文件之中,把载荷信息文件对接于模型文件接口,完成动力荷载对输电铁塔三维有限元数值模型的施加,并采用瞬态动力分析法计算输电铁塔的风振响应;
式中:和分别为时变平均风荷载作用下的节点位移最大值和合成风荷载作用下节点位移最大值,βdi为节点i处位移风振系数;为平均风作用下该节点的位移平均值,m;μ为峰值保证因子,取为3.0;σdi脉动风作用下该节点的位移均方根值,m。
进一步地,步骤4)包括以下具体步骤:
41)根据天气预报数据,编制程序生成下击暴流参数,将天气数据转换成荷载数据,把天气预报中的最大风速Vmax作为风速数据导入步骤2)b)中得到下击暴流风场结构竖向风剖面,根据风场结构模型及参数,生成下击暴流时变平均风(见步骤2)b))及风荷载(见步骤2)d));
42)根据步骤3)B)计算得到结构整体风振系数,再将所述结构整体风振系数与步骤2)d)计算得到的下击暴流风荷载相乘得出等效静力风荷载,将载荷信息文件与模型文件接口对接,并利用有限元软件对铁塔进行准静态分析。
进一步,步骤5)包括以下具体步骤:
根据准静态分析结果,对输电铁塔塔顶位移进行校核;遍历整个下击暴流过程,计算最不利时刻关键主材杆件的应力值,选择应力最大的杆件,计算该杆件在整个下击暴流过程的应力,完成对输电铁塔的安全评估。
一种评估下击暴流下输电铁塔安全性的系统,包括以下功能模块:
模型建立模块:根据输电铁塔设计信息,划分输电铁塔塔段,设置风荷载模拟点,建立输电铁塔三维有限元数值模型,在模型文件中预设载荷信息文件接口;
风荷载生成模块:确定下击暴流风场结构,完成对风场结构各组成成分的模拟,遍历风荷载模拟点,生成各风荷载模拟点处的风速以及风荷载;
风振系数计算模块:将载荷信息文件对接至输电铁塔三维有限元数值模型中,对铁塔进行强动力分析,根据分析结果计算输电铁塔结构整体风振系数;
准静态分析模块:根据天气预报数据,生成下击暴流参数,并对铁塔进行准静态分析;
安全评估模块:利用准静态分析结果对输电铁塔结构进行快速安全评估。
一种评估下击暴流下输电铁塔安全性的系统的存储介质,存储以下功能模块:
模型建立模块:根据输电铁塔设计信息,划分输电铁塔塔段,设置风荷载模拟点,建立输电铁塔三维有限元数值模型,在模型文件中预设载荷信息文件接口;
风荷载生成模块:确定下击暴流风场结构,完成对风场结构各组成成分的模拟,遍历风荷载模拟点,生成各风荷载模拟点处的风速以及风荷载;
风振系数计算模块:将载荷信息文件对接至输电铁塔三维有限元数值模型中,对铁塔进行强动力分析,根据分析结果计算输电铁塔结构整体风振系数;
准静态分析模块:根据天气预报数据,生成下击暴流参数,并对铁塔进行准静态分析;
安全评估模块:利用准静态分析结果对输电铁塔结构进行安全评估。
本发明所达到的有益效果:本发明的方法根据输电铁塔设计信息,划分输电铁塔塔段,设置风荷载模拟点,建立输电铁塔三维有限元数值模型,在模型文件中预设载荷信息文件接口;确定下击暴流风场结构,完成对风场结构各组成成分的模拟,遍历风荷载模拟点,生成各风荷载模拟点处的风速以及风荷载;将载荷导入输电铁塔三维数值模型中,对铁塔进行强动力分析,根据分析结果计算输电铁塔结构整体风振系数;根据天气预报数据,生成下击暴流参数,计算各点时变平均风,用结构整体风振系数将平均风放大,得到等效静力风荷载,对铁塔进行准静态分析,利用分析结果评估输电铁塔结构安全性。本发明的方法可用于强对流天气预测路线覆盖范围的大量铁塔快速计算,提高铁塔抗强对流天气的快速应急管理能力。
附图说明
图1为本发明一实施例中输电铁塔单线模型的示意图;
图2为输电铁塔塔段示意图;
图3为输电铁塔风荷载模拟点示意图;
图4为输电铁塔有限元数值模型的示意图;
图5为输电铁塔不同高度处时变平均风的风速时程曲线图;
图6为输电铁塔塔顶位移时程图;
图7为输电铁塔最不利时刻关键位置主材应力云图;
图8为输电铁塔最不利主材杆件应力时程图;
图9为实施例1输电铁塔主材应力时程响应对比图;
图10为实施例1输电铁塔塔顶位移时程响应对比图;
图11为实施例2输电铁塔杆塔结构示意图;
图12为实施例2输电铁塔主材应力时程响应对比图;
图13为实施例2输电铁塔塔顶位移时程响应对比图;
图14为本发明的评估下击暴流下输电铁塔安全性方法的流程图。
具体实施方式
实施例1
以下结合某输电线路中的一个典型输电铁塔对本发明方案作进一步说明,所结合的输电铁塔(见图1)杆塔类型为直线塔,塔高57米,呼高40m,评估该输电铁塔在下击暴流下的安全状态,包括以下步骤1)—步骤5):
1)根据输电铁塔设计信息,设置风荷载模拟点,建立输电铁塔三维有限元数值模型。
所述输电铁塔设计信息中包括且不限于输电铁塔的材料汇总图。输电铁塔的材料汇总图中包括构件结构尺寸、材料牌号内容,根据该内容,能够将输电铁塔划分为不同的塔段和不同的分块区域。
11)根据实际输电线路中铁塔的材料汇总图,将输电铁塔划分为不同的塔段以及不同的分块区域,把每个分块区域简化为一个模拟点进行风荷载的模拟,以风荷载模拟点风速代替所在分块区域所有节点的风速。
12)根据实际输电线路中铁塔的材料汇总图,利用有限元软件对输电铁塔进行三维有限元数值建模,采用梁单元对铁塔杆件进行模拟,单元之间的连接方式为固接。对梁单元进行设置,分别设置杆件的杨氏模量、泊松比、密度、截面形状等。
13)将建立完成的输电铁塔三维有限元数值模型保存成为模型文件,并在文件中预留接口与载荷信息文件对接。
如图3-4所示,本实施例中,将输电铁塔以横隔层为分界划分成12个塔段(见图2),以及26个不同的分块区域,把每个分块区域简化为一个模拟点进行风速模拟。利用有限元软件建立输电铁塔三维有限元数值模型,约束模型底部所有自由度,将模型命令流存储于模型文件之中,在文件中预留与载荷信息文件对接的接口。
2)确定下击暴流风场结构,完成对风场结构各组成成分的模拟,遍历风荷载模拟点,生成各风荷载模拟点处的风速以及风荷载。
22)竖向风剖面为:
式中,z为风载荷模拟点处高度,Vmax为下击暴流的最大水平风速,Zmax为最大水平风速的产生高度。
径向风剖面表达为:
式中:Vr,max最大径向风速,rmax径向风速最大出距雷暴中心的距离,r观测点(构筑物)距雷暴中心的距离,Rr雷暴径向特征尺寸,T为时间常数。
23)脉动风速可表达为:
24)将各模拟点的时变平均风与脉动风相结合,计算出各点的下击暴流合成风速,并通过下式确定输电铁塔上各模拟点的下击暴流风荷载:
式中,ρ为空气密度,取1.226kg/m3;Cx、Cy分别为x、y方向阻力系数,角钢取1.3,导线、地线、悬索、绝缘子以及拉线取1.0;Ax和Ay分别为x、y方向输电铁塔各节点所辖的迎风面积。
在本实施例中,假设y方向为顺导、地线方向,x向为垂直于导、地线方向,先模拟各风速模拟点位置处的时变平均风、脉动风以及下击暴流合成风速,计算各模拟点所辖的迎风面积,再根据风荷载计算式计算出每个风速模拟点的风荷载。本案例给出输电铁塔不同高度处时变平均风的风速时程曲线(见图5),根据脉动表达式计算出脉动风后,将平均风速与脉动风速相加,即可得到下击暴流合成风速。表1给出了时变平均风风速的模拟参数。
表1:时变平均风风速的模拟参数
3)将载荷信息文件对接至输电铁塔三维数值模型中,对铁塔进行强动力分析,根据分析结果计算输电铁塔结构整体风振系数。
31)将各模拟点的风荷载信息存储于载荷信息文件之中,把载荷信息文件对接于模型文件接口,完成动力荷载对输电铁塔三维有限元数值模型的施加,并采用瞬态动力分析法计算输电铁塔的风振响应。
式中:和分别为时变平均风荷载作用下的节点位移最大值和合成风荷载作用下节点位移最大值,βdi为节点i处位移风振系数;为平均风作用下该节点的位移平均值,m;μ为峰值保证因子,取为3.0;σdi脉动风作用下该节点的位移均方根值,m。
在本实施例中,将载荷信息存储于文件之中,并与模型文件的接口对接,对铁塔模型施加风荷载。先计算输电铁塔在时变平均风作用下的响应,得到各节点在平均风作用下的位移平均值,再计算输电铁塔在脉动风作用下的响应,得到各节点在脉动风作用下的位移均方根值。根据节点i的和σdi即可计算出该节点的位移风振系数βdi,再根据各节点的位移风振系数计算出输电铁塔结构的整体风振系数。最后提取出输电铁塔在下击暴流合成风速作用下的风振响应结果,以备用。
4)根据天气预报数据,生成下击暴流参数,计算各点时变平均风,用结构整体风振系数将平均风放大,得到等效静力风荷载,并对铁塔进行准静态分析。
41)根据天气预报数据生成下击暴流参数,将天气数据转换成荷载数据。根据风场结构模型及参数,生成下击暴流时变平均风。
42)利用结构整体风振系数将平均风放大,计算等效静力风荷载,将载荷信息文件与模型文件接口对接,对铁塔进行准静态分析。
在本案例中,根据预先获得的天气预报数据,通过程序转换成荷载数据,将载荷信息文件与模型文件对接,完成对输电铁塔等效静力风荷载的施加,计算铁塔的风振响应,并提取出该工况下的风振响应结果,以备用。
5)利用准静态分析结果对输电铁塔结构的安全评估。
I)根据准静态分析结果,对输电铁塔塔顶位移进行校核;遍历整个下击暴流过程,计算得出最不利时刻关键主材杆件的应力值,选择响应最大的杆件,计算该杆件在整个下击暴流过程的应力,完成对输电铁塔的快速安全评估。
在本实施例中,根据有限元软件对输电铁塔结构的准静态分析结果,计算输电铁塔的塔顶位移时程(见图6);遍历整个下击暴流过程,计算最不利时刻关键主材杆件的应力值(见图7),选择响应最大的杆件,计算该杆件在整个下击暴流过程的应力(见图8),根据规范完成输电铁塔安全状态的评估。此外,比较了在下击暴流作用和等效静风荷载作用下的主材应力与塔顶位移时程图对比(见图9和图10),验证了该方法的准确性。
一种评估下击暴流下输电铁塔安全性的系统,包括以下功能模块:
模型建立模块:根据输电铁塔设计信息,设置风荷载模拟点,建立输电铁塔三维有限元数值模型;
风荷载生成模块:确定下击暴流风场结构,完成对风场结构各组成成分的模拟,遍历风荷载模拟点,生成各风荷载模拟点处的风速以及风荷载;
风振系数计算模块:将载荷信息文件对接至输电铁塔三维有限元数值模型中,对铁塔进行强动力分析,根据分析结果计算输电铁塔结构整体风振系数;
准静态分析模块:根据天气预报数据,生成下击暴流参数,计算各点时变平均风,用结构整体风振系数将平均风放大,得到等效静力风荷载,并对铁塔进行准静态分析;
安全评估模块:利用准静态分析结果对输电铁塔结构进行安全评估。
一种评估下击暴流下输电铁塔安全性的系统的存储介质,存储以下功能模块:
模型建立模块:根据输电铁塔设计信息,设置风荷载模拟点,建立输电铁塔三维有限元数值模型;
风荷载生成模块:确定下击暴流风场结构,完成对风场结构各组成成分的模拟,遍历风荷载模拟点,生成各风荷载模拟点处的风速以及风荷载;
风振系数计算模块:将载荷信息文件对接至输电铁塔三维有限元数值模型中,对铁塔进行强动力分析,根据分析结果计算输电铁塔结构整体风振系数;
准静态分析模块:根据天气预报数据,生成下击暴流参数,计算各点时变平均风,用结构整体风振系数将平均风放大,得到等效静力风荷载,并对铁塔进行准静态分析;
快速安全评估模块:利用准静态分析结果对输电铁塔结构进行快速安全评估。
实施例2
针对另外一条输电线路中的一个典型输电铁塔进行分析,所结合的输电铁塔(见图11)杆塔类型为直线塔,塔高38.7米,呼高28m。根据有限元软件对输电铁塔结构的准静态分析结果,比较了在下击暴流作用和等效静风荷载作用下的最不利主材应力与塔顶位移时程图对比(见图12和图13),同样验证了该方法的准确性。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护之内。
Claims (4)
1.一种评估下击暴流下输电铁塔安全性的方法,其特征在于,包括以下步骤:
步骤1)根据输电铁塔设计信息,设置风荷载模拟点,建立输电铁塔三维有限元数值模型;
步骤2)确定下击暴流风场结构,完成对风场结构各组成成分的模拟,遍历风荷载模拟点,生成各风荷载模拟点处的风速以及风荷载,包括以下具体步骤:
式中:V(z)为竖向风剖面平均风速;f(t)为最大值为1的时间函数,描述输电铁塔处竖向风速剖面随时间的变化;Vc(t)是输电铁塔处合成水平风速;Vt为风暴中心移动速度,Vr(r,t)为径向风剖面风速,θ为输电线路与径向风剖面风速的夹角,α是合成水平风速与输电线路的夹角,r为输电铁塔至风暴中心距离;
b)竖向风剖面表述为:
式中,z为风载荷模拟点处的高度,Vmax为下击暴流的最大水平风速,Zmax为最大水平风速的产生高度;
径向风剖面模型Vr(r,t)为:
式中:Vr,max为最大径向风速,rmax为径向风速最大出距雷暴中心的距离,r为观测点距雷暴中心的距离,Rr为雷暴径向特征尺寸,T为时间常数;
c)脉动风速表达为:
d)将各模拟点的时变平均风与脉动风相结合,计算出各点的下击暴流合成风速,并通过下式确定输电铁塔上各模拟点的下击暴流风荷载:
式中,Fx(t)为模拟点x方向的风荷载;Fy(t)模拟点y方向的风荷载;V(t)为下击暴流风速;ρ为空气密度;Cx、Cy分别为x、y方向阻力系数;Ax和Ay分别为x、y方向输电铁塔各节点所辖的迎风面积;
步骤3)将载荷信息文件对接至输电铁塔三维有限元数值模型中,对输电铁塔进行瞬态动力分析,根据分析结果计算输电铁塔结构整体风振系数,包括以下具体步骤:
A)将各模拟点的风荷载信息存储于载荷信息文件之中,把载荷信息文件对接于模型文件接口,完成动力荷载对输电铁塔三维有限元数值模型的施加,并采用瞬态动力分析法计算输电铁塔的风振响应;
式中:和分别为时变平均风荷载作用下的节点位移最大值和合成风荷载作用下节点位移最大值,βdi为节点i处位移风振系数;为平均风作用下该节点的位移平均值;μ为峰值保证因子;σdi脉动风作用下该节点的位移均方根值;
步骤4)根据天气预报数据,生成下击暴流参数,并对铁塔进行准静态分析,包括以下具体步骤:
41)根据天气预报数据生成下击暴流参数,将天气数据转换成荷载数据,把天气预报中的最大风速Vmax作为风速数据导入步骤2)中的步骤b)中得到下击暴流风场结构竖向风剖面,根据风场结构模型及参数,生成下击暴流时变平均风及风荷载;
42)根据步骤3)中的步骤B)计算得到结构整体风振系数,再将所述结构整体风振系数与步骤2)中的步骤d)计算得到的下击暴流风荷载相乘得出等效静力风荷载,将载荷信息文件与模型文件接口对接,并利用有限元软件对铁塔进行准静态分析;
步骤5)利用准静态分析结果对输电铁塔结构进行安全评估,包括以下过程:
根据准静态分析结果,对输电铁塔塔顶位移进行校核;遍历整个下击暴流过程,计算输电铁塔主材杆件的应力值,选择应力最大的杆件,计算该杆件在整个下击暴流过程的应力,完成对输电铁塔的安全评估。
2.根据权利要求1所述的评估下击暴流下输电铁塔安全性的方法,其特征在于:在步骤1)中,具体包括以下步骤:
11)根据实际输电铁塔的设计信息,将输电铁塔划分为多个塔段以及多个分块区域,把每个分块区域简化为一个模拟点进行风荷载模拟,以风荷载模拟点的风速代替所在分块区域所有节点的风速;
12)根据实际输电铁塔的设计信息,利用有限元软件对输电铁塔进行三维有限元数值建模,采用梁单元对输电铁塔杆件进行模拟,所述梁单元之间的连接方式为固接;对梁单元进行设置,分别设置所述输电铁塔杆件的参数;
13)将建立完成的所述输电铁塔三维有限元数值模型保存成为模型文件,并在文件中预设载荷信息文件接口以与载荷信息文件对接。
3.一种评估下击暴流下输电铁塔安全性的系统,其特征在于,包括以下功能模块:
模型建立模块:根据输电铁塔设计信息,设置风荷载模拟点,建立输电铁塔三维有限元数值模型;
风荷载生成模块:确定下击暴流风场结构,完成对风场结构各组成成分的模拟,遍历风荷载模拟点,生成各风荷载模拟点处的风速以及风荷载,包括以下具体步骤:
式中:V(z)为竖向风剖面平均风速;f(t)为最大值为1的时间函数,描述输电铁塔处竖向风速剖面随时间的变化;Vc(t)是输电铁塔处合成水平风速;Vt为风暴中心移动速度,Vr(r,t)为径向风剖面风速,θ为输电线路与径向风剖面风速的夹角,α是合成水平风速与输电线路的夹角,r为输电铁塔至风暴中心距离;
b)竖向风剖面表述为:
式中,z为风载荷模拟点处的高度,Vmax为下击暴流的最大水平风速,Zmax为最大水平风速的产生高度;
径向风剖面模型Vr(r,t)为:
式中:Vr,max为最大径向风速,rmax为径向风速最大出距雷暴中心的距离,r为观测点距雷暴中心的距离,Rr为雷暴径向特征尺寸,T为时间常数;
c)脉动风速表达为:
d)将各模拟点的时变平均风与脉动风相结合,计算出各点的下击暴流合成风速,并通过下式确定输电铁塔上各模拟点的下击暴流风荷载:
式中,Fx(t)为模拟点x方向的风荷载;Fy(t)模拟点y方向的风荷载;V(t)为下击暴流风速;ρ为空气密度;Cx、Cy分别为x、y方向阻力系数;Ax和Ay分别为x、y方向输电铁塔各节点所辖的迎风面积;
风振系数计算模块:将载荷信息文件对接至输电铁塔三维有限元数值模型中,对输电铁塔进行瞬态动力分析,根据分析结果计算输电铁塔结构整体风振系数,包括以下具体步骤:
A)将各模拟点的风荷载信息存储于载荷信息文件之中,把载荷信息文件对接于模型文件接口,完成动力荷载对输电铁塔三维有限元数值模型的施加,并采用瞬态动力分析法计算输电铁塔的风振响应;
式中:和分别为时变平均风荷载作用下的节点位移最大值和合成风荷载作用下节点位移最大值,βdi为节点i处位移风振系数;为平均风作用下该节点的位移平均值;μ为峰值保证因子;σdi脉动风作用下该节点的位移均方根值;
准静态分析模块:根据天气预报数据,生成下击暴流参数,并对铁塔进行准静态分析,包括以下步骤:
41)根据天气预报数据生成下击暴流参数,将天气数据转换成荷载数据,把天气预报中的最大风速Vmax作为风速数据导入风荷载生成模块中的步骤b)中得到下击暴流风场结构竖向风剖面,根据风场结构模型及参数,生成下击暴流时变平均风及风荷载;
42)根据风振系数计算模块中的步骤B)计算得到结构整体风振系数,再将所述结构整体风振系数与风荷载生成模块中的步骤d)计算得到的下击暴流风荷载相乘得出等效静力风荷载,将载荷信息文件与模型文件接口对接,并利用有限元软件对铁塔进行准静态分析;
安全评估模块:利用准静态分析结果对输电铁塔结构进行安全评估,包括以下过程:
根据准静态分析结果,对输电铁塔塔顶位移进行校核;遍历整个下击暴流过程,计算输电铁塔主材杆件的应力值,选择应力最大的杆件,计算该杆件在整个下击暴流过程的应力,完成对输电铁塔的安全评估。
4.一种评估下击暴流下输电铁塔安全性的存储介质,存储以下功能模块:
模型建立模块:根据输电铁塔设计信息,设置风荷载模拟点,建立输电铁塔三维有限元数值模型;
风荷载生成模块:确定下击暴流风场结构,完成对风场结构各组成成分的模拟,遍历风荷载模拟点,生成各风荷载模拟点处的风速以及风荷载,包括以下具体步骤:
式中:V(z)为竖向风剖面平均风速;f(t)为最大值为1的时间函数,描述输电铁塔处竖向风速剖面随时间的变化;Vc(t)是输电铁塔处合成水平风速;Vt为风暴中心移动速度,Vr(r,t)为径向风剖面风速,θ为输电线路与径向风剖面风速的夹角,α是合成水平风速与输电线路的夹角,r为输电铁塔至风暴中心距离;
b)竖向风剖面表述为:
式中,z为风载荷模拟点处的高度,Vmax为下击暴流的最大水平风速,Zmax为最大水平风速的产生高度;
径向风剖面模型Vr(r,t)为:
式中:Vr,max为最大径向风速,rmax为径向风速最大出距雷暴中心的距离,r为观测点距雷暴中心的距离,Rr为雷暴径向特征尺寸,T为时间常数;
c)脉动风速表达为:
d)将各模拟点的时变平均风与脉动风相结合,计算出各点的下击暴流合成风速,并通过下式确定输电铁塔上各模拟点的下击暴流风荷载:
式中,Fx(t)为模拟点x方向的风荷载;Fy(t)模拟点y方向的风荷载;V(t)为下击暴流风速;ρ为空气密度;Cx、Cy分别为x、y方向阻力系数;Ax和Ay分别为x、y方向输电铁塔各节点所辖的迎风面积;
风振系数计算模块:将载荷信息文件对接至输电铁塔三维有限元数值模型中,对输电铁塔进行瞬态动力分析,根据分析结果计算输电铁塔结构整体风振系数,包括以下具体步骤:
A)将各模拟点的风荷载信息存储于载荷信息文件之中,把载荷信息文件对接于模型文件接口,完成动力荷载对输电铁塔三维有限元数值模型的施加,并采用瞬态动力分析法计算输电铁塔的风振响应;
式中:和分别为时变平均风荷载作用下的节点位移最大值和合成风荷载作用下节点位移最大值,βdi为节点i处位移风振系数;为平均风作用下该节点的位移平均值;μ为峰值保证因子;σdi脉动风作用下该节点的位移均方根值;
准静态分析模块:根据天气预报数据,生成下击暴流参数,并对铁塔进行准静态分析,包括以下步骤:
41)根据天气预报数据生成下击暴流参数,将天气数据转换成荷载数据,把天气预报中的最大风速Vmax作为风速数据导入风荷载生成模块中的步骤b)中得到下击暴流风场结构竖向风剖面,根据风场结构模型及参数,生成下击暴流时变平均风及风荷载;
42)根据风振系数计算模块中的步骤B)计算得到结构整体风振系数,再将所述结构整体风振系数与风荷载生成模块中的步骤d)计算得到的下击暴流风荷载相乘得出等效静力风荷载,将载荷信息文件与模型文件接口对接,并利用有限元软件对铁塔进行准静态分析;
安全评估模块:利用准静态分析结果对输电铁塔结构进行安全评估,包括以下过程:根据准静态分析结果,对输电铁塔塔顶位移进行校核;遍历整个下击暴流过程,计算输电铁塔主材杆件的应力值,选择应力最大的杆件,计算该杆件在整个下击暴流过程的应力,完成对输电铁塔的安全评估。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011210259.0A CN112507414B (zh) | 2020-11-03 | 2020-11-03 | 评估下击暴流下输电铁塔安全性的方法、系统及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011210259.0A CN112507414B (zh) | 2020-11-03 | 2020-11-03 | 评估下击暴流下输电铁塔安全性的方法、系统及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112507414A CN112507414A (zh) | 2021-03-16 |
CN112507414B true CN112507414B (zh) | 2022-09-06 |
Family
ID=74955251
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011210259.0A Active CN112507414B (zh) | 2020-11-03 | 2020-11-03 | 评估下击暴流下输电铁塔安全性的方法、系统及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112507414B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113361212B (zh) * | 2021-06-02 | 2022-05-03 | 重庆科技学院 | 下击暴流作用下输电塔顺风向响应频域分析方法 |
CN114662350A (zh) * | 2022-05-26 | 2022-06-24 | 国网山东省电力公司泰安供电公司 | 风沙耦合载荷模拟方法、系统及终端 |
CN117077234B (zh) * | 2023-10-18 | 2024-01-30 | 中国铁塔股份有限公司 | 一种用于确定格构塔的有限元模型的方法与设备 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104933628A (zh) * | 2014-03-18 | 2015-09-23 | 国家电网公司 | 一种下击暴流作用下输电铁塔风荷载确定方法 |
-
2020
- 2020-11-03 CN CN202011210259.0A patent/CN112507414B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104933628A (zh) * | 2014-03-18 | 2015-09-23 | 国家电网公司 | 一种下击暴流作用下输电铁塔风荷载确定方法 |
Non-Patent Citations (1)
Title |
---|
"下击暴流作用下双柱悬索拉线塔—线体系风振响应研究";邵鹏飞;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20190415;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112507414A (zh) | 2021-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112507414B (zh) | 评估下击暴流下输电铁塔安全性的方法、系统及存储介质 | |
Amirzadeh et al. | A computational framework for the analysis of rain-induced erosion in wind turbine blades, part I: Stochastic rain texture model and drop impact simulations | |
CN105260953A (zh) | 一种综合风电场理论出力评估仿真系统和评估方法 | |
CN102663251A (zh) | 基于计算流体力学模型的风电场功率物理预测方法 | |
CN104699936A (zh) | 基于cfd短期风速预测风电场的扇区管理方法 | |
CN104894968A (zh) | 一种可抗强/台风的大跨度桥梁用桥塔 | |
CN114611365A (zh) | 基于脉动风影响下输电杆塔结构的动力学分析方法及系统 | |
CN110847210A (zh) | 一种海上风机基础结构以及轻量化设计方法 | |
CN115659756A (zh) | 一种输电杆塔防风性能分析方法 | |
CN112036061A (zh) | 全方位长短腿塔线体系有限元建模及动力响应分析方法 | |
CN111639450A (zh) | 基于支持向量机的塔式起重机损伤识别方法 | |
Vemuri et al. | Wind Analysis of Tall Buildings Using Codal Provisions, Stochastic Approach and CFD–A Comparative Study | |
Xu et al. | Collapse analysis of a wind turbine tower with initial-imperfection subjected to near-field ground motions | |
CN116433091A (zh) | 一种基于环境模拟的钢结构质量检测系统 | |
Li1a et al. | A review of the transmission tower-line system performance under typhoon in wind tunnel test | |
CN115062432A (zh) | 一种基于数据驱动的输电杆塔倾斜预警方法及系统 | |
Santiago-Hernández et al. | Predicting topographic effect multipliers in complex terrain with shallow neural networks | |
Cheng et al. | Full-scale/model test comparisons to validate the traditional ABL wind tunnel simulation technique: a literature review | |
Htet | Wind Load Design of a Triangular Shaped Building Using Finite Element Analysis | |
Li et al. | Pressure distributions on prism-shaped buildings in experimentally simulated downburst | |
Sheikh Aleslami et al. | Effect of Wind Load on Milad Tower and Adjacent Buildings Case study: Shed Adjacent to Tower | |
CN107143469A (zh) | 基于gps定位的风机移动监测方法、装置及系统 | |
CN211312581U (zh) | 一种海上风机基础结构 | |
CN117034711B (zh) | 一种基于风力机叶片有限元模型的分布力加载方法 | |
Jamińska-Gadomska et al. | Wind flow around a church–Case study |
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 |