CN112115616B - 一种输电塔的动力学特性分析方法及装置 - Google Patents

一种输电塔的动力学特性分析方法及装置 Download PDF

Info

Publication number
CN112115616B
CN112115616B CN202010996555.1A CN202010996555A CN112115616B CN 112115616 B CN112115616 B CN 112115616B CN 202010996555 A CN202010996555 A CN 202010996555A CN 112115616 B CN112115616 B CN 112115616B
Authority
CN
China
Prior art keywords
transmission tower
power transmission
equivalent
stiffness
segmented
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
CN202010996555.1A
Other languages
English (en)
Other versions
CN112115616A (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.)
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Original Assignee
Electric Power Research Institute of Guangdong Power Grid 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 Guangdong Power Grid Co Ltd filed Critical Electric Power Research Institute of Guangdong Power Grid Co Ltd
Priority to CN202010996555.1A priority Critical patent/CN112115616B/zh
Publication of CN112115616A publication Critical patent/CN112115616A/zh
Application granted granted Critical
Publication of CN112115616B publication Critical patent/CN112115616B/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/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
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • 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

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

一种输电塔的动力学特性分析方法及装置
技术领域
本发明涉及输电规划技术领域,尤其涉及一种输电塔的动力学特性分析方法及装置。
背景技术
鉴于输电塔是一种对风载荷非常敏感的结构,容易发生震动疲劳损伤甚至在极端条件下倒塌,因此需要对输电塔的动力学特性进行分析,确保输电塔及塔线体系在安全范围内。
现阶段输电塔的动力学特性分析方法,大部分情况下通过串联多质点系力学模型进行分析。串联多质点系力学模型是输电塔体系风振响应计算分析中使用最广泛的一种简化模型,其原理是将输电塔简化为由多个集中质量串联起来的多自由度体系。但串联多质点系力学模型在建立时一般只考虑一阶横向弯曲振型,且必须先通过有限元分析得到简化模型的柔度矩阵,再求逆得到简化模型的参数,才能根据简化模型分析输电塔的动力学特性,其复杂的构建过程导致整个输电塔的动力学特性的分析过程用时相对较长,分析效率并不能满足现阶段的研究要求。
发明内容
本发明提供了一种输电塔的动力学特性分析方法及装置,用于快速分析输电塔的动力学特性,确保输电塔及塔线体系在安全范围内。
本发明提供的一种输电塔的动力学特性分析方法,包括:
获取输电塔信息,所述输电塔信息包括:输电塔结构、受力特点、截面尺寸和节点坐标;
根据所述输电塔结构、受力特点对输电塔进行分段,并用等效构件替代分段后的输电塔中的每个分段,得到输电塔简化示意图;
基于所述截面尺寸和节点坐标,结合预先设定的截面属性计算公式,计算等效构件的截面属性;
组合所述输电塔简化示意图和所述截面属性,生成输电塔简化模型;
根据所述输电塔简化模型和每个等效构件的风载荷属性,分析所述输电塔的动力学特性。
可选地,所述根据所述输电塔结构、受力特点对输电塔进行分段,并用等效构件替代分段后的输电塔中的每个分段,得到输电塔简化示意图的步骤,包括:
若相邻的等效构件不能直接连接,则用刚性构件连接所述相邻的等效构件。
可选地,所述截面属性包括:轴向刚度、抗弯刚度、扭转刚度、剪切刚度、等效密度;所述基于所述截面尺寸和分段后的输电塔的节点坐标和拓扑结构,结合预先设定的截面属性计算公式,计算等效构件的截面属性的步骤包括:
根据轴向刚度计算公式计算所述轴向刚度,所述轴向刚度计算公式为:
Figure BDA0002692714650000021
式中,ki,N是第i个分段等效构件的轴向刚度,li是第i个分段的长度,ni是第i个分段中构件的数量,Ei,j、Ai,j是第i个分段中第j根构件的弹性模量、横截面面积,θi,j是第j根构件与第i个分段等效构件轴向的倾角;
根据抗弯刚度计算公式计算所述抗弯刚度,所述抗弯刚度计算公式为:
Figure BDA0002692714650000022
Figure BDA0002692714650000023
式中,ki,M是第i个分段等效构件的抗弯刚度,ki,N是第i个分段等效构件的轴向刚度,ri是第i个分段等效构件横截面的等效回转半径,hi,1和hi,2分别是第i个分段两端横向构件的长度,li是第i个分段的长度;
根据剪切刚度计算公式计算所述剪切刚度,所述剪切刚度计算公式为:
Figure BDA0002692714650000024
式中,ki,S是第i个分段等效构件的剪切刚度,Ei,j、Ai,j是第i个分段中第j根构件的弹性模量、横截面面积,
Figure BDA0002692714650000031
是第j根构件与第i个分段等效构件剪切方向的夹角;
根据扭转刚度计算公式计算所述扭转刚度,所述扭转刚度计算公式为:
Figure BDA0002692714650000032
式中,ki,T是第i个分段等效构件的扭转刚度,ki,S是第i个分段等效构件的剪切刚度,ri是第i个分段等效构件横截面的等效回转半径,li是第i个分段的长度;
根据等效密度计算公式计算所述等效密度,所述等效密度计算公式为:
Figure BDA0002692714650000033
式中:ρi是第i个分段等效构件的密度,ρi,j、Ai,j、li,j分别是第i个分段中第j根构件的密度、横截面面积和长度,ni是第i个分段中构件的数量,mk是第i个分段中第k个节点板或螺栓等结构的质量。
可选地,所述计算等效构件的截面属性的步骤,还包括:
若所述等效构件相邻的横向构件的长度不相同,则分别计算所述所述等效构件的抗弯刚度和剪切刚度。
可选地,所述基于所述输电塔简化模型,计算所述等效构件的风载荷属性,分析所述输电塔的动力学特性的步骤包括:
从预先设定风载荷属性参数表中确定所有所述等效构件的风载荷属性;
结合所有所述等效构件的风载荷属性和所述输电塔简化模型,分析所述输电塔的动力学特性。
为了解决上述问题,本申请实施例还公开了一种输电塔的动力学特性分析装置,包括:
获取模块,用于获取输电塔信息,所述输电塔信息包括:输电塔结构、受力特点、截面尺寸和节点坐标;
分段模块,用于根据所述输电塔结构、受力特点对输电塔进行分段,并用等效构件替代分段后的输电塔中的每个分段,得到输电塔简化示意图;
计算模块,用于基于所述截面尺寸和节点坐标,结合预先设定的截面属性计算公式,计算等效构件的截面属性;
生成模块,用于组合所述输电塔简化示意图和所述截面属性,生成输电塔简化模型;
分析模块,用于根据所述输电塔简化模型和每个等效构件的风载荷属性,分析所述输电塔的动力学特性。
所述分段模块包括:
连接子模块,用于若相邻的等效构件不能直接连接,则用刚性构件连接所述相邻的等效构件。
可选地,所述截面属性包括:轴向刚度、抗弯刚度、扭转刚度、剪切刚度、等效密度;所述计算模块包括:
轴向刚度计算子模块,用于根据轴向刚度计算公式计算所述轴向刚度,所述轴向刚度计算公式为:
Figure BDA0002692714650000041
式中,ki,N是第i个分段等效构件的轴向刚度,li是第i个分段的长度,ni是第i个分段中构件的数量,Ei,j、Ai,j是第i个分段中第j根构件的弹性模量、横截面面积,θi,j是第j根构件与第i个分段等效构件轴向的倾角;
抗弯刚度计算子模块,用于根据抗弯刚度计算公式计算所述抗弯刚度,所述抗弯刚度计算公式为:
Figure BDA0002692714650000042
Figure BDA0002692714650000043
式中,ki,M是第i个分段等效构件的抗弯刚度,ki,N是第i个分段等效构件的轴向刚度,ri是第i个分段等效构件横截面的等效回转半径,hi,1和hi,2分别是第i个分段两端横向构件的长度,li是第i个分段的长度;
剪切刚度计算子模块,用于根据剪切刚度计算公式计算所述剪切刚度,所述剪切刚度计算公式为:
Figure BDA0002692714650000051
式中,ki,S是第i个分段等效构件的剪切刚度,Ei,j、Ai,j是第i个分段中第j根构件的弹性模量、横截面面积,
Figure BDA0002692714650000052
是第j根构件与第i个分段等效构件剪切方向的夹角;
扭转刚度计算子模块,用于根据扭转刚度计算公式计算所述扭转刚度,所述扭转刚度计算公式为:
Figure BDA0002692714650000053
式中,ki,T是第i个分段等效构件的扭转刚度,ki,S是第i个分段等效构件的剪切刚度,ri是第i个分段等效构件横截面的等效回转半径,li是第i个分段的长度;
根据等效密度计算公式计算所述等效密度,所述等效密度计算公式为:
Figure BDA0002692714650000054
式中:ρi是第i个分段等效构件的密度,ρi,j、Ai,j、li,j分别是第i个分段中第j根构件的密度、横截面面积和长度,ni是第i个分段中构件的数量,mk是第i个分段中第k个节点板或螺栓等结构的质量。
可选地,所述计算模块还用于:若所述等效构件相邻的横向构件的长度不相同,则分别计算所述等效构件的抗弯刚度和剪切刚度。
可选地,所述分析模块包括:
风载荷属性确定子模块,用于从预先设定风载荷属性参数表中确定所有所述等效构件的风载荷属性;
动力学特性分析子模块,用于结合所有所述等效构件的风载荷属性和所述输电塔简化模型,分析所述输电塔的动力学特性。
从以上技术方案可以看出,本发明具有以下优点:
本申请实施例通过获取输电塔信息,所述输电塔信息包括:输电塔结构、受力特点、截面尺寸和节点坐标;根据所述输电塔结构、受力特点对输电塔进行分段,并用等效构件替代分段后的输电塔中的每个分段,得到输电塔简化示意图;基于所述截面尺寸和节点坐标,结合预先设定的截面属性计算公式,计算等效构件的截面属性;组合所述输电塔简化示意图和所述截面属性,生成输电塔简化模型;根据所述输电塔简化模型和每个等效构件的风载荷属性,分析所述输电塔的动力学特性。利用建立的简化模型,结合等效构件的风载荷属性,从而快速分析输电塔的动力学特性,解决传统的分析输电塔的动力学特性方法分析效率较低的问题,进而确保输电塔及塔线体系在安全范围内。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。
图1为本发明的一种输电塔的动力学特性分析方法实施例一的步骤流程图;
图2为本发明的一种输电塔的动力学特性分析方法实施例二的步骤流程图;
图3为本发明的一种输电塔的动力学特性分析方法实施例二的输电塔完整模型单线图;
图4为本发明的一种输电塔的动力学特性分析方法实施例二的输电塔简化示意图;
图5为本发明的一种输电塔的动力学特性分析方法实施例二的刚性构件连接图;
图6为本发明的一种输电塔的动力学特性分析方法实施例二的输电塔其中一个原输电塔模型的模型分段几何示意图;
图7为本发明的一种输电塔的动力学特性分析方法实施例二的输电塔横切截面图;
图8为本发明的一种简化模型的前五阶模态;
图9为本发明的简化模型与原输电塔模型的顺线向塔顶位移时程对比示意图;
图10为本发明的简化模型的高处与原输电塔模型高处的位移均方根值对比示意图;
图11为本发明的一种输电塔的动力学特性分析装置实施例的结构框图。
具体实施方式
本发明实施例提供了一种输电塔的动力学特性分析方法及装置,用于快速分析输电塔的动力学特性,确保输电塔及塔线体系在安全范围内。
为使得本发明的发明目的、特征、优点能够更加的明显和易懂,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,下面所描述的实施例仅仅是本发明一部分实施例,而非全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
实施例一,请参阅图1,图1为本发明的一种输电塔的动力学特性分析方法实施例一的步骤流程图,具体可以包括如下步骤:
步骤S101,获取输电塔信息,所述输电塔信息包括:输电塔结构、受力特点、截面尺寸和节点坐标;
步骤S102,根据所述输电塔结构、受力特点对输电塔进行分段,并用等效构件替代分段后的输电塔中的每个分段,得到输电塔简化示意图;
步骤S103,基于所述截面尺寸和节点坐标,结合预先设定的截面属性计算公式,计算等效构件的截面属性;
步骤S104,组合所述输电塔简化示意图和所述截面属性,生成输电塔简化模型;
步骤S105,根据所述输电塔简化模型和每个等效构件的风载荷属性,分析所述输电塔的动力学特性。
在本申请实施例中,本申请实施例通过获取输电塔信息,所述输电塔信息包括:输电塔结构、受力特点、截面尺寸和节点坐标;根据所述输电塔结构、受力特点对输电塔进行分段,并用等效构件替代分段后的输电塔中的每个分段,得到输电塔简化示意图;基于所述截面尺寸和节点坐标,结合预先设定的截面属性计算公式,计算等效构件的截面属性;组合所述输电塔简化示意图和所述截面属性,生成输电塔简化模型;根据所述输电塔简化模型和每个等效构件的风载荷属性,分析所述输电塔的动力学特性。利用建立的简化模型,结合等效构件的风载荷属性,从而快速分析输电塔的动力学特性,解决传统的分析输电塔的动力学特性方法分析效率较低的问题,进而确保输电塔及塔线体系在安全范围内。
实施例二,请参阅图2,为本发明的一种输电塔的动力学特性分析方法实施例二的步骤流程图,具体包括:
步骤S201,获取输电塔信息,所述输电塔信息包括:输电塔结构、受力特点、截面尺寸和节点坐标;
在本发明实施例中,从输电塔的设计图纸当中获得输电塔的信息,其中,输电塔信息包括:输电塔结构、受力特点、截面尺寸和节点坐标等信息。
需要说明的是,现有的输电塔绝大部分情况下是通过角钢一类的构件连接而成的,而各个构件之间的连接点即为构件节点,输电塔结构是指输电塔的各个构件节点的连接关系;受力特点具体是指承受轴力的方向,与输电塔主材的方向基本一致,例如:塔身竖直方向是主要承受轴力的方向,干字塔横担一般水平方向主要承受轴力的方向,猫头塔分支处的斜向方向主要承受轴力的方向;截面尺寸是指输电塔各个构件的分段长度、横截面面积、横向构件的长度等;节点坐标是根据输电塔的设计图纸上记载的设计尺寸计算得出。
步骤S202,根据所述输电塔结构、受力特点对输电塔进行分段,并用等效构件替代分段后的输电塔中的每个分段,得到输电塔简化示意图;
请参阅图3,图3为本发明的一种输电塔的动力学特性分析方法实施例二的输电塔完整模型单线图,在本发明实施例中,根据输电塔的结构、受力特点从输电塔完整模型单线图的端部和主材节点板位置进行自然分段,每个分段用等效构件替代,结合原输电塔模型的节点坐标,得到分段后的得到输电塔的简化示意图,并且简化示意图中的每个等效构件都具备风载荷属性,请参阅图4,图4为本发明的一种输电塔的动力学特性分析方法实施例二的输电塔简化示意图。
在具体实现中,使用梁单元作为等效构件简化分段,且梁单元轴向与原输电塔模型分段主要承受轴力方向保持一致。
步骤S203,若相邻的等效构件不能直接连接,则用刚性构件连接所述相邻的等效构件;
在本发明实施例中,将步骤202中分段后的输电塔中的每个分段两个端部中心点作为等效构件两个节点,相邻两个分段可能存在因等效构件节点坐标缩聚导致没有直接连接的情况,此时用刚性构件进行连接。
请参阅图5,图5为本发明实施例提供的一种输电塔的动力学特性分析方法实施例二的刚性构件连接图,其中A为塔身分段等效构件节点,B为塔身与横担的共同节点,C点为横担分段等效构件节点,1为刚性构件。在具体实现中,在利用等效构件替代原输电塔模型分段时,由于简化模型中一个原输电塔模型分段用一个梁单元模拟,伴随坐标缩聚的原因,横担分段构件节点C与塔身分段等效构件节点A并不在同一位置,为了简化考虑塔身分段与横担的连接关系,将上述两个节点需要使用刚性构件1进行连接,一般地,刚性构件1的截面积、惯性矩等参数与其相邻等效构件的相同,但弹性模量远高于等效构件的弹性模量,用以提高连接的刚度。
步骤S204,基于所述截面尺寸和分段后的输电塔的节点坐标和拓扑结构,结合预先设定的截面属性计算公式,计算等效构件的截面属性,所述截面属性包括:轴向刚度、抗弯刚度、扭转刚度、剪切刚度、等效密度;
请参阅图6,图6为本发明实施例提供的一种输电塔的动力学特性分析方法实施例二的输电塔其中一个原输电塔模型的模型分段几何示意图,其显示有:第i个分段两端横向构件的长度hi,1和hi,2、第i个分段的长度弹性模量li、第j根构件与第i个分段等效构件轴向的倾角θi,j、第j根构件与第i个分段等效构件剪切方向的夹角
Figure BDA0002692714650000091
以及第i个分段中第j根构件的长度li,j
在一个优选实施例中,所述基于所述截面尺寸和分段后的输电塔的节点坐标和拓扑结构,结合预先设定的截面属性计算公式,计算等效构件的截面属性的步骤包括:
根据轴向刚度计算公式计算所述轴向刚度,所述轴向刚度计算公式为:
Figure BDA0002692714650000101
式中,ki,N是第i个分段等效构件的轴向刚度,li是第i个分段的长度,ni是第i个分段中构件的数量,Ei,j、Ai,j是第i个分段中第j根构件的弹性模量、横截面面积,θi,j是第j根构件与第i个分段等效构件轴向的倾角;
根据抗弯刚度计算公式计算所述抗弯刚度,所述抗弯刚度计算公式为:
Figure BDA0002692714650000102
Figure BDA0002692714650000103
式中,ki,M是第i个分段等效构件的抗弯刚度,ki,N是第i个分段等效构件的轴向刚度,ri是第i个分段等效构件横截面的等效回转半径,hi,1和hi,2分别是第i个分段两端横向构件的长度,li是第i个分段的长度;
根据剪切刚度计算公式计算所述剪切刚度,所述剪切刚度计算公式为:
Figure BDA0002692714650000104
式中,ki,S是第i个分段等效构件的剪切刚度,Ei,j、Ai,j是第i个分段中第j根构件的弹性模量、横截面面积,
Figure BDA0002692714650000105
是第j根构件与第i个分段等效构件剪切方向的夹角;
根据扭转刚度计算公式计算所述扭转刚度,所述扭转刚度计算公式为:
Figure BDA0002692714650000106
式中,ki,T是第i个分段等效构件的扭转刚度,ki,S是第i个分段等效构件的剪切刚度,ri是第i个分段等效构件横截面的等效回转半径,li是第i个分段的长度;
根据等效密度计算公式计算所述等效密度,所述等效密度计算公式为:
Figure BDA0002692714650000111
式中:ρi是第i个分段等效构件的密度,ρi,j、Ai,j、li,j分别是第i个分段中第j根构件的密度、横截面面积和长度,ni是第i个分段中构件的数量,mk是第i个分段中第k个节点板或螺栓等结构的质量。
步骤S205,若所述等效构件相邻的横向构件的长度不相同,则分别计算所述等效构件的抗弯刚度和剪切刚度;
请参阅图7,图7为本发明实施例提供的一种输电塔的动力学特性分析方法实施例二的输电塔横切截面图,其包含有横向构件相邻两段的长度hx和hy,由于抗弯刚度和剪切刚度的大小由横向构件的长度决定,在使用抗弯刚度计算公式和剪切刚度计算公式计算等效构件的抗弯刚度和剪切刚度时,对于相邻的hx和hy不相等时,两个对称轴方向的抗弯刚度和剪切刚度不同,需要分开计算。
步骤S206,组合所述分段的输电塔的节点坐标、所述拓扑结构和所述截面属性,生成输电塔简化模型;
在本发明实施例中,将步骤S204得到的截面属性赋予用等效构件替代的分段后的输电塔模型,得到输电塔的简化模型。
步骤S207,从预先设定风载荷属性参数表中确定所有所述等效构件的风载荷属性;
在本发明实施例中,预先设定有风载荷属性参数表,如下所示:
Figure BDA0002692714650000112
Figure BDA0002692714650000121
在具体实现中,根据原输电塔模型分段的属性,结合预先设定的风载荷属性参数表,确定原输电塔模型分段对应的等效构件的风载荷属性。
步骤S208,结合所有所述等效构件的风载荷属性和所述输电塔简化模型,分析所述输电塔的动力学特性。
在具体实现中,将步骤S207确定的风载荷属性和步骤S206的简化模型数据导入至有限元分析软件,如ANSYS和ABAQUS等,得到输电塔的动力学特性。其中,有限元分析软件可以根据可以实际情况进行选择,本发明在此不再限制。
在本申请实施例所提供的一种输电塔的动力学特性分析方法,通过获取输电塔信息;并根据所述输电塔结构、受力特点对输电塔进行分段,并用等效构件替代分段后的输电塔中的每个分段,得到输电塔简化示意图;基于所述截面尺寸和节点坐标,结合预先设定的截面属性计算公式,计算等效构件的截面属性;组合所述输电塔简化示意图和所述截面属性,生成输电塔简化模型;根据所述输电塔简化模型和每个等效构件的风载荷属性,分析所述输电塔的动力学特性。利用建立的简化模型,结合等效构件的风载荷属性,从而快速分析输电塔的动力学特性以及加装阻尼器等装置时的减震效果,进而确保输电塔及塔线体系在安全范围内,与现有的动力学评测方法相比,可以更好反映结构的扭转模态和高阶弯曲模态,因而具有更高的计算精度,同时在与现有技术的分析结果比较接近的情况下,计算效率更高。
为方便本领域技术人员对本发明的有益效果的理解,下面对本发明分析得到的动力学特性与采用现有技术分析得到的动力学特性进行举例说明。
示例一,使用本发明实施例提供的一种输电塔的动力学特性分析方法,分析前五阶固有频率,包括:(1)建立简化模型;(2)结合风载荷属性,运用ANSYS有限元分析软件对该简化模型进行了模态分析,得到图8,图8为本发明的一种简化模型的前五阶模态;(3)基于该前五阶模态得到前五阶固有频率。同时使用原输电塔模型(即完整模型)得到的前五阶固有频率。两种方法得到的固有频率比较结果如下表所示:
Figure BDA0002692714650000131
可以看出,使用本发明实施例提供的一种输电塔的动力学特性分析方法得出的结果,与使用原输电塔模型得出的结果相比,最大误差在2.21%,即相比于原输电塔模型得出的结果相差并不大。
示例二,对塔底进行固定约束,在塔顶施加大小为10kN的顺线向集中力从而模拟实际情况,然后使用本发明实施例提供的一种输电塔的动力学特性分析方法和现有的分析方法分析得到的结构位移响应,分析结果请参阅图9和图10,图9为本发明的简化模型与原输电塔模型的顺线向塔顶位移时程对比示意图,其中①为简化模型计算得到的顺线向塔顶位移时程,②为原输电塔模型计算得到的顺线向塔顶位移时程,可以看出,基于原输电塔模型计算的常规分析方法得出的点,基本上排布在简化模型的线上;图10为本发明的简化模型的高处与原输电塔模型高处的位移均方根值对比示意图,①为简化模型计算得到的均方根变化曲线,②为原输电塔模型计算得到的均方根变化点集合,可以看出,基于完整模型计算得到的多个均方根点,几乎全部分布在简化模型计算得到的均方根值变化曲线上。概括而言,本发明提供的一种输电塔的动力学特性分析方法能很好地反映输电塔的动力学特性,具有较高的精度。
请参阅图11,示出了一种输电塔的动力学特性分析装置实施例的结构框图,包括如下模块:
获取模块101,用于获取输电塔信息,所述输电塔信息包括:输电塔结构、受力特点、截面尺寸和节点坐标;
分段模块102,用于根据所述输电塔结构、受力特点对输电塔进行分段,并用等效构件替代分段后的输电塔中的每个分段,得到输电塔简化示意图;
在本发明的一个优选实施例中,所述分段模块102包括:
连接子模块,用于若相邻的等效构件不能直接连接,则用刚性构件连接所述相邻的等效构件。
计算模块103,用于基于所述截面尺寸和节点坐标,结合预先设定的截面属性计算公式,计算等效构件的截面属性;
在本发明的一个优选实施例中,所述计算模块103,还用于若所述等效构件相邻的横向构件的长度不相同,则分别计算所述等效构件的抗弯刚度和剪切刚度。
在本发明的一个优选实施例中,所述截面属性包括:轴向刚度、抗弯刚度、扭转刚度、剪切刚度、等效密度;所述计算模块103包括:
轴向刚度计算子模块,用于根据轴向刚度计算公式计算所述轴向刚度,所述轴向刚度计算公式为:
Figure BDA0002692714650000141
式中,ki,N是第i个分段等效构件的轴向刚度,li是第i个分段的长度,ni是第i个分段中构件的数量,Ei,j、Ai,j是第i个分段中第j根构件的弹性模量、横截面面积,θi,j是第j根构件与第i个分段等效构件轴向的倾角;
抗弯刚度计算子模块,用于根据抗弯刚度计算公式计算所述抗弯刚度,所述抗弯刚度计算公式为:
Figure BDA0002692714650000142
Figure BDA0002692714650000143
式中,ki,M是第i个分段等效构件的抗弯刚度,ki,N是第i个分段等效构件的轴向刚度,ri是第i个分段等效构件横截面的等效回转半径,hi,1和hi,2分别是第i个分段两端横向构件的长度,li是第i个分段的长度;
剪切刚度计算子模块,用于根据剪切刚度计算公式计算所述剪切刚度,所述剪切刚度计算公式为:
Figure BDA0002692714650000151
式中,ki,S是第i个分段等效构件的剪切刚度,Ei,j、Ai,j是第i个分段中第j根构件的弹性模量、横截面面积,
Figure BDA0002692714650000152
是第j根构件与第i个分段等效构件剪切方向的夹角;
扭转刚度计算子模块,用于根据扭转刚度计算公式计算所述扭转刚度,所述扭转刚度计算公式为:
Figure BDA0002692714650000153
式中,ki,T是第i个分段等效构件的扭转刚度,ki,S是第i个分段等效构件的剪切刚度,ri是第i个分段等效构件横截面的等效回转半径,li是第i个分段的长度;
根据等效密度计算公式计算所述等效密度,所述等效密度计算公式为:
Figure BDA0002692714650000154
式中:ρi是第i个分段等效构件的密度,ρi,j、Ai,j、li,j分别是第i个分段中第j根构件的密度、横截面面积和长度,ni是第i个分段中构件的数量,mk是第i个分段中第k个节点板或螺栓等结构的质量。
生成模块104,用于组合所述输电塔简化示意图和所述截面属性,生成输电塔简化模型;
分析模块105,用于根据所述输电塔简化模型和每个等效构件的风载荷属性,分析所述输电塔的动力学特性。
在本发明的一个优选实施例中,所述分析模块105包括:
风载荷属性确定子模块,用于从预先设定风载荷属性参数表中确定所有所述等效构件的风载荷属性;
动力学特性分析子模块,用于结合所有所述等效构件的风载荷属性和所述输电塔简化模型,分析所述输电塔的动力学特性。
对于装置实施例而言,由于其与方法实施例基本相似,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的装置和模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在本申请所提供的几个实施例中,应该理解到,所揭露的装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-OnlyMemory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (8)

1.一种输电塔的动力学特性分析方法,其特征在于,包括:
获取输电塔信息,所述输电塔信息包括:拓扑结构、受力特点、截面尺寸和节点坐标;
根据所述拓扑结构、受力特点对输电塔进行分段,并用等效构件替代分段后的输电塔中的每个分段,得到输电塔简化示意图;
基于所述截面尺寸、分段后的输电塔的节点坐标和所述拓扑结构,结合预先设定的截面属性计算公式,计算等效构件的截面属性;
组合所述输电塔简化示意图和所述截面属性,生成输电塔简化模型;
根据所述输电塔简化模型和每个等效构件的风载荷属性,分析所述输电塔的动力学特性;
所述截面属性包括:轴向刚度、抗弯刚度、扭转刚度、剪切刚度、等效密度;所述基于所述截面尺寸、分段后的输电塔的节点坐标和所述拓扑结构,结合预先设定的截面属性计算公式,计算等效构件的截面属性的步骤包括:
根据轴向刚度计算公式计算所述轴向刚度,所述轴向刚度计算公式为:
Figure FDA0003153747430000011
式中,ki,N是第i个分段等效构件的轴向刚度,li是第i个分段的长度,ni是第i个分段中构件的数量,Ei,j、Ai,j是第i个分段中第j根构件的弹性模量、横截面面积,θi,j是第j根构件与第i个分段等效构件轴向的倾角;
根据抗弯刚度计算公式计算所述抗弯刚度,所述抗弯刚度计算公式为:
Figure FDA0003153747430000012
Figure FDA0003153747430000013
式中,ki,M是第i个分段等效构件的抗弯刚度,ki,N是第i个分段等效构件的轴向刚度,ri是第i个分段等效构件横截面的等效回转半径,hi,1和hi,2分别是第i个分段两端横向构件的长度,li是第i个分段的长度;
根据剪切刚度计算公式计算所述剪切刚度,所述剪切刚度计算公式为:
Figure FDA0003153747430000021
式中,ki,S是第i个分段等效构件的剪切刚度,Ei,j、Ai,j是第i个分段中第j根构件的弹性模量、横截面面积,
Figure FDA0003153747430000022
是第j根构件与第i个分段等效构件剪切方向的夹角,ni是第i个分段中构件的数量;
根据扭转刚度计算公式计算所述扭转刚度,所述扭转刚度计算公式为:
Figure FDA0003153747430000023
式中,ki,T是第i个分段等效构件的扭转刚度,ki,S是第i个分段等效构件的剪切刚度,ri是第i个分段等效构件横截面的等效回转半径,li是第i个分段的长度;
根据等效密度计算公式计算所述等效密度,所述等效密度计算公式为:
Figure FDA0003153747430000024
式中:ρi是第i个分段等效构件的等效密度,ρi,j、Ai,j、li,j分别是第i个分段中第j根构件的密度、横截面面积和长度,ni是第i个分段中构件的数量,mk是第i个分段中第k个节点板或螺栓结构的质量。
2.根据权利要求1所述的输电塔的动力学特性分析方法,其特征在于,所述根据所述拓扑结构、受力特点对输电塔进行分段,并用等效构件替代分段后的输电塔中的每个分段,得到输电塔简化示意图的步骤,包括:
若相邻的等效构件不能直接连接,则用刚性构件连接所述相邻的等效构件。
3.根据权利要求1所述的输电塔的动力学特性分析方法,其特征在于,所述计算等效构件的截面属性的步骤,还包括:
若所述等效构件相邻的横向构件的长度不相同,则分别计算所述等效构件的抗弯刚度和剪切刚度。
4.根据权利要求1所述的输电塔的动力学特性分析方法,其特征在于,所述根据所述输电塔简化模型和每个等效构件的风载荷属性,分析所述输电塔的动力学特性的步骤包括:
从预先设定风载荷属性参数表中确定所述每个等效构件的风载荷属性;
结合所述每个等效构件的风载荷属性和所述输电塔简化模型,分析所述输电塔的动力学特性。
5.一种输电塔的动力学特性分析装置,其特征在于,包括:
获取模块,用于获取输电塔信息,所述输电塔信息包括:拓扑结构、受力特点、截面尺寸和节点坐标;
分段模块,用于根据所述拓扑结构、受力特点对输电塔进行分段,并用等效构件替代分段后的输电塔中的每个分段,得到输电塔简化示意图;
计算模块,用于基于所述截面尺寸、分段后的输电塔的节点坐标和所述拓扑结构,结合预先设定的截面属性计算公式,计算等效构件的截面属性;
生成模块,用于组合所述输电塔简化示意图和所述截面属性,生成输电塔简化模型;
分析模块,用于根据所述输电塔简化模型和每个等效构件的风载荷属性,分析所述输电塔的动力学特性;
所述截面属性包括:轴向刚度、抗弯刚度、扭转刚度、剪切刚度、等效密度;所述计算模块包括:
轴向刚度计算子模块,用于根据轴向刚度计算公式计算所述轴向刚度,所述轴向刚度计算公式为:
Figure FDA0003153747430000031
式中,ki,N是第i个分段等效构件的轴向刚度,li是第i个分段的长度,ni是第i个分段中构件的数量,Ei,j、Ai,j是第i个分段中第j根构件的弹性模量、横截面面积,θi,j是第j根构件与第i个分段等效构件轴向的倾角;
抗弯刚度计算子模块,用于根据抗弯刚度计算公式计算所述抗弯刚度,所述抗弯刚度计算公式为:
Figure FDA0003153747430000041
Figure FDA0003153747430000042
式中,ki,M是第i个分段等效构件的抗弯刚度,ki,N是第i个分段等效构件的轴向刚度,ri是第i个分段等效构件横截面的等效回转半径,hi,1和hi,2分别是第i个分段两端横向构件的长度,li是第i个分段的长度;
剪切刚度计算子模块,用于根据剪切刚度计算公式计算所述剪切刚度,所述剪切刚度计算公式为:
Figure FDA0003153747430000043
式中,ki,S是第i个分段等效构件的剪切刚度,Ei,j、Ai,j是第i个分段中第j根构件的弹性模量、横截面面积,
Figure FDA0003153747430000044
是第j根构件与第i个分段等效构件剪切方向的夹角,ni是第i个分段中构件的数量;
扭转刚度计算子模块,用于根据扭转刚度计算公式计算所述扭转刚度,所述扭转刚度计算公式为:
Figure FDA0003153747430000045
式中,ki,T是第i个分段等效构件的扭转刚度,ki,S是第i个分段等效构件的剪切刚度,ri是第i个分段等效构件横截面的等效回转半径,li是第i个分段的长度;
根据等效密度计算公式计算所述等效密度,所述等效密度计算公式为:
Figure FDA0003153747430000046
式中:ρi是第i个分段等效构件的等效密度,ρi,j、Ai,j、li,j分别是第i个分段中第j根构件的密度、横截面面积和长度,ni是第i个分段中构件的数量,mk是第i个分段中第k个节点板或螺栓结构的质量。
6.根据权利要求5所述的输电塔的动力学特性分析装置,其特征在于,所述分段模块包括:
连接子模块,用于若相邻的等效构件不能直接连接,则用刚性构件连接所述相邻的等效构件。
7.根据权利要求5所述的输电塔的动力学特性分析装置,其特征在于,
所述计算模块,还用于若所述等效构件相邻的横向构件的长度不相同,则分别计算所述等效构件的抗弯刚度和剪切刚度。
8.根据权利要求5所述的输电塔的动力学特性分析装置,其特征在于,所述分析模块包括:
风载荷属性确定子模块,用于从预先设定风载荷属性参数表中确定所有所述等效构件的风载荷属性;
动力学特性分析子模块,用于结合所有所述等效构件的风载荷属性和所述输电塔简化模型,分析所述输电塔的动力学特性。
CN202010996555.1A 2020-09-21 2020-09-21 一种输电塔的动力学特性分析方法及装置 Active CN112115616B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010996555.1A CN112115616B (zh) 2020-09-21 2020-09-21 一种输电塔的动力学特性分析方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010996555.1A CN112115616B (zh) 2020-09-21 2020-09-21 一种输电塔的动力学特性分析方法及装置

Publications (2)

Publication Number Publication Date
CN112115616A CN112115616A (zh) 2020-12-22
CN112115616B true CN112115616B (zh) 2021-09-24

Family

ID=73800305

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010996555.1A Active CN112115616B (zh) 2020-09-21 2020-09-21 一种输电塔的动力学特性分析方法及装置

Country Status (1)

Country Link
CN (1) CN112115616B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114065547B (zh) * 2021-11-26 2024-06-25 北京机电工程研究所 一种固支梁多自由度系统的振动响应谱计算方法
CN115048694B (zh) * 2022-06-09 2023-05-16 国网山东省电力公司临沂供电公司 输电塔体系振型聚类方法、装置以及计算机设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105426599A (zh) * 2015-11-10 2016-03-23 华北电力大学(保定) 一种简化的拉线塔主柱有限元模型的计算方法
CN108763827A (zh) * 2018-06-21 2018-11-06 广东电网有限责任公司电力科学研究院 一种输电塔有限元模型建立方法及装置
CN111651806A (zh) * 2020-03-31 2020-09-11 重庆科技学院 基于惯性力法的常规输电塔设计风载荷计算方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002041195A2 (en) * 2000-11-17 2002-05-23 Battelle Memorial Institute Structural stress analysis
CN103675197B (zh) * 2012-08-31 2015-09-02 上海宇航系统工程研究所 半刚性电池基板的结构动力学分析方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105426599A (zh) * 2015-11-10 2016-03-23 华北电力大学(保定) 一种简化的拉线塔主柱有限元模型的计算方法
CN108763827A (zh) * 2018-06-21 2018-11-06 广东电网有限责任公司电力科学研究院 一种输电塔有限元模型建立方法及装置
CN111651806A (zh) * 2020-03-31 2020-09-11 重庆科技学院 基于惯性力法的常规输电塔设计风载荷计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于改进多质点模型的输电塔_线体系动力特性评估;周颖 等;《特种结构》;20200228;第37卷(第1期);第87-93页 *

Also Published As

Publication number Publication date
CN112115616A (zh) 2020-12-22

Similar Documents

Publication Publication Date Title
Liu et al. Aeroelastic analysis of a floating offshore wind turbine in platform‐induced surge motion using a fully coupled CFD‐MBD method
CN107346357B (zh) 一种基于整体耦合模型的海上风机疲劳分析系统
CN112115616B (zh) 一种输电塔的动力学特性分析方法及装置
Rezaei et al. Modal-based damage identification for the nonlinear model of modern wind turbine blade
Liang et al. Constructing simplified models for dynamic analysis of monopile-supported offshore wind turbines
CN103984813A (zh) 一种离心压缩机裂纹叶轮结构的振动建模与分析方法
CN106570204A (zh) 一种基于cpu+gpu异构并行计算的透平机械叶片静强度特性分析方法
CN105808829A (zh) 一种基于cpu+gpu异构并行计算的透平机械叶片固有频率特性分析方法
CN107529644B (zh) 一种电力系统静态电压稳定域边界线性近似方法
CN109558621B (zh) 一种结构损伤识别方法及系统
CN102682175A (zh) 基于屈曲模态组合的网格结构施工误差可靠性分析方法
Stieng et al. Load case reduction for offshore wind turbine support structure fatigue assessment by importance sampling with two‐stage filtering
CN111539138A (zh) 基于阶跃函数的结构动力学峰值时域响应灵敏度求解方法
CN109101752B (zh) 一种复杂水工建筑物局部结构自振频率计算方法
CN104123471A (zh) 重型施工机械重载结构安全评估方法
CN115049317B (zh) 一种风力资源评估工具的选择方法、装置及电子设备
CN113836625A (zh) 一种基于能力谱法的输电塔抗震性能等级划分方法及系统
CN116882104A (zh) 机床薄弱环节的识别方法、装置、电子设备及存储介质
CN106295159B (zh) 一种基于自相关函数的结构风振响应高效频域估计方法
CN107423509B (zh) 一种输电线塔体系气弹模型构建方法及装置
CN114117847A (zh) 一种轮毂圆锥滚子轴承疲劳寿命预测方法
CN109726454B (zh) 管路系统的流固耦合建模方法及装置
CN104778325B (zh) 基于表面单元的面载荷处理方法及装置
CN113901552A (zh) 一种堰体模型的建立方法及装置
Araújo et al. Modal identification of a wind turbine.

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