CN111651804A - 基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法 - Google Patents

基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法 Download PDF

Info

Publication number
CN111651804A
CN111651804A CN202010245662.0A CN202010245662A CN111651804A CN 111651804 A CN111651804 A CN 111651804A CN 202010245662 A CN202010245662 A CN 202010245662A CN 111651804 A CN111651804 A CN 111651804A
Authority
CN
China
Prior art keywords
tower
wind
coefficient
line
power transmission
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
CN202010245662.0A
Other languages
English (en)
Other versions
CN111651804B (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.)
Chongqing University of Science and Technology
Original Assignee
Chongqing University of Science and Technology
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 Chongqing University of Science and Technology filed Critical Chongqing University of Science and Technology
Priority to CN202010245662.0A priority Critical patent/CN111651804B/zh
Publication of CN111651804A publication Critical patent/CN111651804A/zh
Application granted granted Critical
Publication of CN111651804B publication Critical patent/CN111651804B/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/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • 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
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/08Construction
    • 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
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Economics (AREA)
  • Computer Hardware Design (AREA)
  • Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • General Health & Medical Sciences (AREA)
  • General Business, Economics & Management (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Tourism & Hospitality (AREA)
  • Strategic Management (AREA)
  • Primary Health Care (AREA)
  • Marketing (AREA)
  • Mathematical Analysis (AREA)
  • Water Supply & Treatment (AREA)
  • Public Health (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Wind Motors (AREA)

Abstract

本发明公开了一种基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法,步骤为:搭建常规输电塔的塔线体系,获取塔线体系的物理参数;基于塔线耦合影响因子,求取杆塔等效阻尼系数、塔线体系常规输电塔的风振系数、塔线体系悬垂绝缘子串最大风偏角的风振系数;塔线体系风荷载脉动折减系数;对塔线体系常规输电塔的风振系数、风偏角的风振系数进行修正计算,得到塔线体系常规输电塔的修正风振系数和塔线体系输电线的修正风振系数;在等效振动惯性力作用下计算塔线体系中常规输电塔的设计风荷载;基于塔线分离法计算输电线设计风荷载。有益效果:单塔设计精度高,可靠。

Description

基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电 塔、线设计风载荷的计算方法
技术领域
本发明涉及技术领域,具体的说是一种基于惯性力法和塔线分离法考虑塔 线耦合影响的常规输电塔、线设计风载荷的计算方法。
背景技术
常规输电塔,是与超高输电塔相比较,其塔高在梯度风高度范围内,未采 用钢管混凝土复合材料的输电塔。相对而言,常规输电塔在输电线路中更加普 遍地被采用。
采用杆塔设计规范计算的风荷载简明、方便、省时,现阶段该方法仍被设 计人员广泛采用。通过规范计算的风荷载应该具有能够使输电塔的风振响应和 实际最大风振响应等效的作用。采用准确的效静力风荷载进行杆塔设计是保证 输电线路正常运行的前提。现有电力相关标准中:例如文献(1)GB 50545-2010. 110kV~750kV架空输电线路设计规范[S].北京:中国计划出版社,2010;(2) GB 50665-2011.1000kV架空输电线路设计规范[S].北京:中国计划出版社, 2011;(3)DL/T 5154-2012.架空输电线路杆塔结构设计技术规定[S].北京:中 国计划出版社,2012和(4)DL/T 5504-2015.特高压架空输电线路大跨越设计技 术规定[S].北京:中国计划出版社,2015中给出了60m以下的单塔风振系数的取 值,并推荐了60m以上时采用荷载规范计算风振系数。荷载规范的风振系数适 用于外形和质量有规律变化的密实建筑物。输电塔为格构式构筑物,横担和横 隔面的质量和挡风面积对风振系数的影响大。另外,采用钢管混凝土是荷载规 范计算风振系数时没有考虑到的。采用随机振动理论计算输电塔的等效静力风 荷载时,表达式涉及复杂的多重积分,并且输电塔的外形和质量分布不规律, 很难用一个统一的表达式概括。并且,导线发生风振时的气动阻尼随平均风速 的增加而增大,风振响应的共振分量因气动阻尼而大幅降低,计算中可忽略不 计。并且常规输电塔组成的单塔体系还需要考虑塔线耦合效应产生的,这对最 终计算出高精度的单塔体系风荷载具有重要的意义。
发明内容
针对上述问题,本发明提供了一种基于惯性力法和塔线分离法考虑塔线耦 合影响的常规输电塔、线设计风载荷的计算方法,来提高塔线体系风载荷计算 精度。为达到上述目的,本发明采用的具体技术方案如下:
一种基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计 风载荷的计算方法,具体步骤为:S1:搭建常规输电塔的塔线体系,并获取塔 线体系的常规输电塔、输电线、绝缘子串的物理参数;上述数据至少包括常规 输电塔所在地面粗糙度类别、设定10m高度处的平均分速
Figure BDA0002433914740000021
常规输电塔的 总高度H、跟开b1、横担个数nc、横担平均外伸长度
Figure BDA0002433914740000022
以及输电塔与导线、绝 缘子串的布置方案;还包括导线线性和导线线长等。S2:基于塔线耦合影响因 子,根据塔线体系中杆塔等效阻尼系数δe;S3:将步骤S2得到的塔线体系中杆 塔等效阻尼系数δe来替换阻尼系数ζ1,求取塔线体系常规输电塔的风振系数 β(z);考虑线形与线长影响因子,计算塔线体系悬垂绝缘子串最大风偏角的风振 系数β;S4:考虑塔线耦合效应,求取塔线体系风荷载脉动折减系数εc;S5:根 据步骤S4得到的塔线体系风荷载脉动折减系数,对步骤S3中的塔线体系常规 输电塔的风振系数、塔线体系输电线的风振系数进行修正计算,得到塔线体系 常规输电塔的修正风振系数β*(z)和塔线体系输电线的修正风振系数β*;S6:基 于塔线分离法,在等效振动惯性力作用下计算塔线体系中常规输电塔的设计风 荷载fESWL(z)和输电线的设计风荷载WX
再进一步的技术方案,步骤S2的具体步骤为:
S21:根据步骤S1的常规输电塔的塔线体系,得到常规输电塔塔线体系计 算模型图;所述塔线体系计算模型中的杆塔为密实结构,塔身为正方形的变截 面,由下至上尺寸变小,横担为等截面;所述塔线体系计算模型中的导线两端 等高,与固定铰支座连接;所述塔线体系计算模型中的杆塔高度为H,横担悬 臂长度为lca,绝缘子长度为lin,导线跨度为L。导线挂点无高差。S22:设定输 电线和绝缘子串振动的假设条件,得到的常规输电塔塔线体系中输电线和绝缘 子串的振型图以及迎风面、被风面输电线和绝缘子串的广义质量、广义刚度和 广义阻尼;并将常规输电塔塔线体系中输电线和绝缘子串组合形成索结构体系;所述导线和绝缘子串振动的假设条件为:迎风面和背风面的导线在风荷载下的 振动同步;绝缘子串的频率和阻尼比由导线控制,与导线的一致。所述迎风面、 被风面导线的广义质量计算公式为:
Figure BDA0002433914740000031
所述迎风面、被风面导线的广义刚度计算公式为:
Figure BDA0002433914740000032
所述迎风面、 被风面导线的广义阻尼计算公式为:
Figure BDA0002433914740000033
mc为单根导线单位线长的质 量;单根导线振型
Figure BDA0002433914740000034
γg为导线的自重比载;σ0为导线的水平初 应力;Γ为导线的线长,
Figure BDA0002433914740000035
ζc=ζscac;ζsc为导线结构阻尼比;ζac为导线 启动阻尼比;Nc为分裂导线的个数;Tw为平均风状态下单根导线的水平张力;ζc为导线阻尼比;所述迎风面、被风面绝缘子串的广义质量计算公式为:
Figure BDA0002433914740000041
所述迎风面、被风面绝缘子串的广义刚度计算公式为:
Figure BDA0002433914740000042
所述迎风面、被风面绝缘子串的广义阻尼计算公式为:
Figure BDA0002433914740000043
其中,min为绝缘子串单位高度质量;Din为绝缘子串迎风外径; 绝缘子串振型
Figure BDA0002433914740000044
H-lin≤z≤H;ζin为绝缘子串阻尼比;
Figure BDA0002433914740000045
lin为绝缘 子长度;其中,索结构体系对应的广义质量、广义刚度和广义阻尼的计算公式 为:
Figure BDA0002433914740000046
本发明中,下标ci表示索结构。
S23:将步骤S22得到的数据构建索结构体系结合杆塔结构组成塔线耦合简 化计算模型;S24:基于杆塔结构组成塔线耦合简化计算模型,求取常规输电塔 塔线体系下杆塔顺风向位移的共振分量的均方值和单塔时杆塔顺风向位移共振 分量的均方值;从而得到二者的比例式;所述塔线体系下杆塔顺风向位移的共 振分量的均方值的计算公式为:
Figure BDA0002433914740000047
Figure BDA0002433914740000049
Figure BDA00024339147400000410
Figure BDA00024339147400000411
杆塔1阶模态的振型φt(z)=(z/H)2,0≤z≤H;
Figure BDA0002433914740000051
为杆塔的广义质量,
Figure BDA0002433914740000052
Mca为横担的质量,mt(z)为随高 度变化的杆塔单位高度质量;coh(z1,z2)为z1和z2高度处两点的脉动风速的相干函 数;Sf(nt)为归一化风速谱,nt为杆塔脉动风速的频率;
Figure BDA0002433914740000053
σv'为 脉动风速的标准差;
Figure BDA00024339147400000511
为索结构与杆塔的广义质量比值,
Figure BDA0002433914740000054
λn为索结构与杆塔的频率比值;λn=nci/nt;导线悬挂于杆塔的顶部,
Figure BDA0002433914740000055
ζt为 总阻尼比;ζt=ζstat;ζst为杆塔结构阻尼比;ωt为杆塔无阻尼振动的圆频率;δci为 索结构总阻尼比,近视取导线阻尼比,δci≈ζc,ζc=ζscac
ρa为空气密度,μs(z)为风压随高度变化系数;bs(z)为随高度变化的迎风面宽 度;
Figure BDA0002433914740000056
随高度变化的平均风速,σv'为脉动风速的标准差;
ζat为杆塔气动阻尼比;
Figure BDA0002433914740000057
As,ca为横 担的挡风面积;所述单塔时杆塔顺风向位移共振分量的均方值为:
Figure BDA0002433914740000058
所述塔线体系下杆塔顺风向位移的共振分量和所述单塔时杆塔顺风向位移 共振分量的比例式为:
Figure BDA00024339147400000510
S25:基于步骤S24得到的计算公式,推导常规输电塔悬挂输电线后塔线等 效阻尼系数的计算公式,并计算塔线体系中杆塔等效阻尼系数。计算公式的步 骤为:悬挂导线后杆塔的等效阻尼比为:
Figure BDA0002433914740000061
其中,ρ与ζe的关系为:
Figure BDA0002433914740000062
对于输电塔线体系而言,索结构为柔性 体系,卓越频率远小于杆塔的频率;则忽略λn的高阶项;杆塔的阻尼比约为0.01, 索结构的阻尼比小于1,则忽略
Figure BDA0002433914740000063
项;导线悬挂于杆塔的顶部,
Figure BDA0002433914740000064
故悬 挂导线后塔线等效阻尼系数的计算公式为:
Figure BDA0002433914740000069
再进一步的技术方案,步骤S3求取塔线体系常规输电塔的风振系数β(z)的 步骤为:S311:根据步骤S1中常规输电塔的物理参数,确定常规输电塔所在地 面粗糙度类别,设定10m高度处的平均分速
Figure BDA0002433914740000065
常规输电塔的总高度H;跟 开b1;横担个数nc;横担平均外伸长度
Figure BDA0002433914740000066
自立式输电塔分为横隔面、横担和 剩余塔身3部分;横隔面、横担、剩余塔身的质量和挡风面积沿高度的分布规 律不同,在计算过程中需要区别对待。S312:构建常规输电塔的风荷载的计算 模型,通过水平均布荷载作用下结构的挠曲线获得常规输电塔0°风向角的1阶 侧弯振型φ1(z);
Figure BDA0002433914740000067
z为实际高度值。对于该弯曲振型,有如下积分关系:
Figure BDA0002433914740000068
根据荷载规范引入背景分量因子 Bz(z),进而计算输电塔的脉动风荷载在水平方向的相关系数ρx和脉动风荷载在 竖直方向的相关系数ρz
根据荷载规范引入的共振分量因子R;确定地面粗糙度指数α;峰值因子gs;10m高度处的湍流度I10
Figure BDA0002433914740000071
b为外轮廓宽度;共振分量因子R 的计算公式为:
Figure BDA0002433914740000072
其中,ξ1=δe
Figure BDA0002433914740000073
n 为脉动风速的频率;当塔高H确定时,Bz(z)的表达式为:
Figure BDA0002433914740000075
S314:根据地面粗糙度类别获取背景分量因子的中间变量γ的拟合系数kγ和 aγ;根据常规输电塔的宽度、深度得到风振系数考虑整体外形变化的修正系数 θv;根据常规输电塔的横担平均外伸长度、总高度、横担个数得到风振系数考 虑附加面积的修正系数θa与风振系数考虑附加质量的修正系数θm的乘积θl
S315:简化并得到剩余塔身的风振系数考虑局部外形变化的修正系数的计 算公式并计算得到对应的剩余塔身的风振系数考虑局部外形变化的修正系数的 θb(z);计算得到横担的风振系数考虑局部外形变化的修正系数θb(zI)和横隔面 的风振系数考虑局部外形变化的修正系数θb(zJ);
其中,
Figure BDA0002433914740000076
S316:根据塔身的实际高度值z,结合步骤S315对应得到的风振系数考虑 局部外形变化的修正系数计算公式计算所述背景分量因子Bz(z),具体计算公式 为:
Figure BDA0002433914740000077
其中,
Figure BDA0002433914740000081
θv是风振系数考虑整体外形变化的修正系数;θη是风 振系数考虑脉动风空间相关性的修正系数;
其中,
Figure BDA0002433914740000082
Figure BDA0002433914740000083
S317:计算风振系数β(z);
Figure BDA0002433914740000084
其中,gs为峰值因子。
再进一步的技术方案,塔线体系悬垂绝缘子串最大风偏角的风振系数步骤 为β;S321:根据步骤S1中常规输电塔塔线体系的输电线、绝缘子串的物理参 数,以重力和平均风荷载作用下作为导线和悬垂绝缘子串计算的初始条件,通 过LRC方法确定悬垂绝缘子串风偏角的计算模型;所述导线物理参数至少包括 导线型号、导线计算截面积、导线弹性模量、线密度、导线外径;所述输电塔 上的绝缘子串物理参数至少包括绝缘子串长度、绝缘子串弹性模量、绝缘子串 质量、绝缘子串挡风面积。S322:计算常规输电塔塔线体系中输电塔之间导线 单位面积的等效静力风荷载;所述输电塔之间导线单位面积的等效静力风荷载pESWL的计算公式为:
Figure BDA0002433914740000085
式中,(:,i)表示矩阵的第i列元素;
Figure BDA0002433914740000086
为 等效背景风压;
Figure BDA0002433914740000087
为平均风荷载;导线在风荷载作用下的振动方程矩阵表达式为:
Figure BDA0002433914740000091
式中,
Figure BDA0002433914740000092
Y′分别为脉动风荷载作用下导线节点顺风向的加速度、速度 和位移;
Figure BDA0002433914740000093
为平均风荷载作用下导线节点顺风向的位移
M为质量矩阵;C为阻尼矩阵;K刚度矩阵;Ls为节点从属面积矩阵;
导线在脉动风荷载作用下的振动方程矩阵表达式为:
Figure BDA0002433914740000094
根据上述内容可以得到等效静力风荷载计算悬垂绝缘子串的最大风偏角; 等效静力风荷载计算悬垂绝缘子串的最大风偏角的计算公式为:
Figure BDA0002433914740000095
式中,
Figure BDA0002433914740000096
为脉动风荷载作用下B点的顺风向峰值位移
Figure BDA0002433914740000097
IAB为A、B 两点间的绝缘子串长度;
Figure BDA0002433914740000098
为平均风荷载作用下B点的顺风向位移,
Figure BDA00024339147400000911
为平均风偏角;具体计算公式为:
Figure BDA00024339147400000912
Figure BDA00024339147400000913
Gv分别为目标点处悬垂绝缘子串的平均风荷载和竖向重力荷载;Wh、Wv分别为目标点处导线传递给悬垂绝缘子串的平均风荷载和竖向荷载。
S323:计算塔线体系悬垂绝缘子串最大风偏角的风振系数β;
Figure BDA00024339147400000914
c表示对计算域内的元素 进行求和;Γc为计算域内导线的线长;
Figure BDA00024339147400000915
为平均风荷载;
Figure BDA00024339147400000916
为等效背景风压。
另外的,为了设计塔线体系最好求取到等效静力风荷载计算悬垂绝缘子串 的最大风偏角
Figure BDA00024339147400000917
具体的:
Figure BDA00024339147400000918
式中,
Figure BDA00024339147400000919
为脉动风荷载作用下B点的顺 风向峰值位移
Figure BDA00024339147400000920
IAB为A、B两点间的绝缘子串长度;
Figure BDA00024339147400000921
Figure BDA00024339147400000922
为平均风荷载作用下B点的顺风向位移,
Figure BDA00024339147400000923
Figure BDA00024339147400000924
为平均风偏角;具 体计算公式为:
Figure BDA00024339147400000925
Gv分别为目标点处悬垂绝缘子串的平 均风荷载和竖向重力荷载;Wh、Wv分别为目标点处导线传递给悬垂绝缘子串的平 均风荷载和竖向荷载。
再进一步的技术方案为:塔线体系风荷载脉动折减系数εc的计算步骤为:
S41:构建常规输电塔塔线体系计算模型,并得到塔线体系计算模型图;
S42:根据常规输电塔,建立杆塔响应与杆塔风振系数的关系,得到塔高H 处建立杆塔荷载引起塔顶位移的均方根值σut(H)与杆塔风振系数β(H)的关系式;
Figure RE-GDA00026102828300000911
其中,ω0为基本风压;μz(H)为风压随常规单塔高 度变化系数;μs(H)为杆塔随高度阻力系数;bs(H)随高度变化的迎风面宽度;gs为峰值因子;ω1为顺风向1阶模态的自振圆频率;m(H)为随高度变化的单位高 度质量;当导线悬挂于杆塔顶部时,建立导线与导线风振系数的关系,得到导 线荷载引起塔顶位移的均方根值σuc(H)的计算公式;
Figure RE-GDA00026102828300000912
其中,Np为导线的相数;μsc为导线阻力系数;μz(Hc)为风压随导线高度变化系 数Nc为分裂导线的个数;Dc为子导线/地线的计算外径;Lp为水平档距;H为塔 高高度;Et为弹性模量。σut(z)与杆塔风振系数β(z)的关系为:
Figure RE-GDA00026102828300000913
在塔高H处建立σut(H)与β(H)的的关系:
Figure RE-GDA00026102828300000914
S43:根据步骤S41得到的内容,采用SRSS的方法确定塔线体系下杆塔的 峰值响应计算公式;
S44:基于步骤S43的塔线体系下杆塔的峰值响应计算公式,采用塔线分离 方法,引入杆塔风荷载脉动折减系数,并得到所述杆塔的峰值响应计算公式的 等价峰值响应计算公式:
Figure BDA0002433914740000112
表示 杆塔荷载引起的峰值响应,
Figure BDA0002433914740000114
表示输电线荷载引起的峰值响应。
S45:以塔顶位移响应为目标,对步骤S44得到的所述杆塔的峰值响应计 算公式的等价峰值响应计算公式进一步更新,得到带未知导线荷载引起塔顶位 移的均方根值和未知杆塔荷载引起塔顶位移的均方根值的脉动折减系数更新计 算公式;
S46:将步骤S42中的计算得到的导线荷载引起塔顶位移的均方根值和杆塔 荷载引起塔顶位移的均方根值带入步骤S45得到脉动折减系数的更新计算公式 中,得到脉动折减系数的最终计算公式,并计算杆塔风荷载脉动折减系数;
Figure BDA0002433914740000115
其中,
Figure BDA0002433914740000116
再进一步的技术方案,所述塔线体系常规输电塔的修正风振系数β*(z)和所 述塔线体系输电线的修正风振系数β*的计算公式为:
Figure BDA0002433914740000117
再进一步的技术方案,在等效振动惯性力作用下计算塔线体系中常规输电 塔的设计风荷载fESWL(z)与所述塔线体系常规输电塔的修正风振系数β*(z)的关 系式为:
Figure BDA0002433914740000118
其中,ξ1=ξe
Figure BDA0002433914740000119
Figure BDA00024339147400001110
Figure BDA0002433914740000121
m(z)=m(0)μm(z);
Figure BDA0002433914740000122
Sf(n)为归一化风速谱,
Figure BDA0002433914740000123
Iz(z)为z高度处的脉动风湍流密度;
Figure BDA0002433914740000124
I10为10m高度处的脉动风湍流密度;x′1为公式
Figure BDA0002433914740000125
中,n=n1时 的取值,n1为输电塔的1阶模态频率;u1和ηxz1是与风场湍流特性和空间相关性 等有关的系数,分别称为综合影响系数和空间相关性折减系数。
再进一步的技术方案,基于塔线分离法计算输电线设计风荷载WX的计算公 式为:
Figure BDA0002433914740000126
其中,β=α'βc;式中,α′为取值小 于1的风压不均匀系数;μsc为导线阻力系数;βc为风荷载调整系数,计算风偏角时 取1;Dc为子导线/地线的计算外径;Lp为杆塔的水平档距;Bl为覆冰时风荷载的 增大系数;ω0为基本风压;μz为风压随高度变化系数;Bl为覆冰时风荷载的增 大系数;Nc为分裂导线的个数;θ为风向角。
本发明的有益效果:采用等效振动惯性力法,并对阻尼系数进行精细化计 算以及考虑考虑塔线耦合效应产生的情况,对塔线体系中常规输电塔设计风荷 载进行计算。通过塔线分离法,并且考虑了阻尼系数进行精细化计算以及考虑 考虑塔线耦合效应产生的情况对塔线体系输电线设计风荷载进行计算,使最终 设计出的塔线体系更为接近实际,设计精度高。
附图说明
图1是塔线体系计算模型图;
图2是导线和绝缘子串的振型图;
图3是塔线耦合简化计算模型图;
图4是常规输电塔计算图;
图5是悬垂绝缘子串风偏角计算模型示意图;
图6是本发明计算流程图;
图7是塔线体系中杆塔等效阻尼系数计算流程图;
图8是塔线体系常规输电塔的风振系数计算流程图;
图9是塔线体系悬垂绝缘子串最大风偏角的风振系数计算流程图;
图10是塔线体系风荷载脉动折减系数计算流程图;
图11为风向角的定义图。
具体实施方式
下面结合附图对本发明的具体实施方式以及工作原理作进一步详细说明。
一种基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计 风载荷的计算方法,结合图6可以看出,具体步骤为:S1:搭建常规输电塔的 塔线体系,并获取塔线体系的常规输电塔、输电线、绝缘子串的物理参数;
结合图1可以看出为常规输电塔的塔线体系。在本实施例中,物理参数为: H=90m,
Figure BDA0002433914740000131
As,ca=0.372m2,Mca=5514.522kg,
Figure BDA0002433914740000132
μs(z)=1,lin=10m,L=700m,Nc=6,mc=1.917kg/m,Dc=0.032m,μsc=1,T0=70500.446N, min=123.808kg/m,Din=0.36m,α=0.15,
Figure BDA0002433914740000133
以上参数中T0为导线的水平初张力。 来流风向为90°,杆塔不考虑被风面的风荷载,横担的被风面风荷载与迎风面风 荷载相同,不考虑遮挡效应。杆塔和导线的体形系数取1。
S2:基于塔线耦合影响因子,根据塔线体系中杆塔等效阻尼系数δe;具体 的,结合图7可以看出,步骤S2的具体步骤为:S21:根据步骤S1的常规输电 塔的塔线体系,得到常规输电塔塔线体系计算模型图,详见图2。所述塔线体 系计算模型中的杆塔为密实结构,塔身为正方形的变截面,由下至上尺寸变小, 横担为等截面;所述塔线体系计算模型中的导线两端等高,与固定铰支座连接; 所述塔线体系计算模型中的杆塔高度为H,横担悬臂长度为lca,绝缘子长度为 lin,导线跨度为L;导线挂点无高差。S22:设定输电线和绝缘子串振动的假设 条件,得到的常规输电塔塔线体系中输电线和绝缘子串的振型图以及迎风面、 被风面输电线和绝缘子串的广义质量、广义刚度和广义阻尼;并将常规输电塔 塔线体系中输电线和绝缘子串组合形成索结构体系;迎风面、被风面导线的广 义质量计算公式为:
Figure BDA0002433914740000141
所述迎风面、被风面导线的广 义刚度计算公式为:
Figure BDA0002433914740000142
所述迎风面、被风面导线的广义阻尼计算公式 为:
Figure BDA0002433914740000143
mc为单根导线单位线长的质量;单根导线振型
Figure BDA0002433914740000144
γg为导线的自重比载;σ0为导线的水平初应力;Γ为导线的线长,
Figure BDA0002433914740000145
ζc=ζscac;ζsc为导线结构阻尼比;ζac为导线启动阻尼比;Nc为分裂导线的个数; Tw为平均风状态下单根导线的水平张力;ζc为导线阻尼比;
所述迎风面、被风面绝缘子串的广义质量计算公式为:
Figure BDA0002433914740000146
所述迎风面、被风面绝缘子串的广义刚度计算公式为:
Figure BDA0002433914740000147
所述迎风面、 被风面绝缘子串的广义阻尼计算公式为:
Figure BDA0002433914740000148
其中,min为绝缘子串单 位高度质量;Din为绝缘子串迎风外径;绝缘子串振型
Figure BDA0002433914740000149
H-lin≤z≤H;ζin为绝 缘子串阻尼比;
Figure BDA00024339147400001410
lin为绝缘子长度;其中,索结构体系对应的广义质量、 广义刚度和广义阻尼的计算公式为:
Figure BDA0002433914740000151
S23:将步骤S22得到的数据构建索结构体系结合杆塔结构组成塔线耦合简 化计算模型,在本实施例中,该模型详见图3;S24:基于杆塔结构组成塔线耦 合简化计算模型,求取常规输电塔塔线体系下杆塔顺风向位移的共振分量的均 方值和单塔时杆塔顺风向位移共振分量的均方值;从而得到二者的比例式;所 述塔线体系下杆塔顺风向位移的共振分量和所述单塔时杆塔顺风向位移共振分 量的比例式为:
Figure BDA0002433914740000152
所述单塔时杆塔顺风向 位移共振分量的均方值为:
Figure BDA0002433914740000153
所述塔线体系下杆塔顺风向位移的共振分量的均方值的计算公式为:
Figure BDA0002433914740000154
S25:基于步骤S24得到的计算公式,推导常规输电塔悬挂输电线后塔线 等效阻尼系数的计算公式,并计算塔线体系中杆塔等效阻尼系数。
塔线体系中杆塔等效阻尼系数计算公式为:
Figure BDA0002433914740000155
对于输电塔线体系而言,索结构为柔性体系,卓越频率远小于杆塔的频率。 因此,公式(1)可以忽略λn的高阶项。此外,杆塔的阻尼比约为0.01,索结构 的阻尼比小于1,可以忽略
Figure BDA0002433914740000156
项。对于图3的计算模型,导线悬挂于杆塔的顶 部,
Figure BDA0002433914740000161
此时,公式(1)可以简化为:
Figure BDA0002433914740000165
S3:结合图8可以看出,将步骤S2得到的塔线体系中杆塔等效阻尼系数δe来替换阻尼系数ζ1,求取塔线体系常规输电塔的风振系数β(z);在本实施例中, 结合图4可以看出,为常规输电塔计算图;具体步骤:S311:根据步骤S1中常 规输电塔的物理参数,确定常规输电塔所在地面粗糙度类别,设定10m高度处 的平均分速
Figure BDA0002433914740000162
常规输电塔的总高度H;跟开b1;横担个数nc;横担平均外伸 长度
Figure BDA0002433914740000163
S312:构建常规输电塔的风荷载的计算模型,通过水平均布荷载作用下结 构的挠曲线获得常规输电塔0°风向角的1阶侧弯振型φ1(z);
Figure BDA0002433914740000164
z为实际高度值;
在本发明中,结合图10可以定义,当来流风平行于横担轴向时风向角 θ=0°,当来流风平行于导线走向时风向角θ=90°。其中,x向表示横担轴向, y向表示顺线向。
风振系数随风向角的变化不大,并且风向角对塔身风振系数和横担风振系 数的影响是相反的,对整塔而言该影响可以抵消。电力相关标准中仅考虑0o风 向角下输电塔的风振系数。因此,可以忽略风向角对风振系数的影响,其它风 向角下的等效静力风荷载通过DL/T 5154中的风荷载分配系数确定。
其中,电力相关标准包括:GB 50545-2010.110kV~750kV架空输电线路设 计规范[S].北京:中国计划出版社,2010;GB 50665-2011.1000kV架空输电 线路设计规范[S].北京:中国计划出版社,2011;DL/T 5154-2012.架空输电线路 杆塔结构设计技术规定[S].北京:中国计划出版社,2012;DL/T 5504-2015.特 高压架空输电线路大跨越设计技术规定[S].北京:中国计划出版社,2015。
S313:根据荷载规范引入背景分量因子Bz(z),进而计算输电塔的脉动风荷 载在水平方向的相关系数ρx和脉动风荷载在竖直方向的相关系数ρz
S314:根据地面粗糙度类别获取背景分量因子的中间变量γ的拟合系数kγ和 aγ;详见表1;根据常规输电塔的宽度、深度得到风振系数考虑整体外形变化的 修正系数θv;取值详见表2;根据常规输电塔的横担平均外伸长度、总高度、横 担个数得到风振系数考虑附加面积的修正系数θa与风振系数考虑附加质量的修 正系数θm的乘积θl;取值详见表4.
表1 kγ和aγ的取值
地面粗糙度类别 A B C D
k<sub>γ</sub> 1.276 0.910 0.404 0.155
a<sub>γ</sub> 0.186 0.218 0.292 0.376
Figure BDA0002433914740000171
表2常规输电塔的宽度深度均沿高度作同一规律变化时θv的值
b<sub>s</sub>(H)/b<sub>s</sub>(0) 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1
θ<sub>v</sub> 1.00 1.10 1.20 1.32 1.50 1.75 2.08 2.53 3.30 5.60
根据常规输电塔的横担平均外伸长度、总高度、横担个数得到风振系数考虑附加面积的 修正系数θa与风振系数考虑附加质量的修正系数θm的乘积θl
Figure BDA0002433914740000172
Figure BDA0002433914740000173
表4常规输电塔的θl取值
Figure BDA0002433914740000181
S315:简化并得到剩余塔身的风振系数考虑局部外形变化的修正系数的计 算公式并计算得到对应的剩余塔身的风振系数考虑局部外形变化的修正系数的 θb(z);计算得到横担的风振系数考虑局部外形变化的修正系数θb(zI)和横隔面的 风振系数考虑局部外形变化的修正系数θb(zJ);
其中,
Figure BDA0002433914740000182
S316:根据塔身的实际高度值z,结合步骤S315对应得到的风振系数考虑 局部外形变化的修正系数计算公式计算所述背景分量因子Bz(z),具体计算公式 为:
Figure BDA0002433914740000183
其中,
Figure BDA0002433914740000184
θv是风振系数考虑整体外形变化的修正系数,取值详 见表2;θη是风振系数考虑脉动风空间相关性的修正系数;θη是风振系数考虑脉 动风空间相关性的修正系数,取值详见表3.
表3常规输电塔的θn取值
Figure BDA0002433914740000185
其中,
Figure BDA0002433914740000191
Figure BDA0002433914740000192
S317:计算风振系数β(z);其中,风振系数表达式为:
Figure BDA0002433914740000193
在本实施例,计算塔线体系输电线的风振系数β计算公式为;
Figure BDA0002433914740000194
Figure BDA0002433914740000195
c表示对计算域内的元素进行求和;Γc为计算域内导线的线长;
Figure BDA0002433914740000196
为平均风荷载;
Figure BDA0002433914740000197
为等效背景风压。
本实施例中,结合图9可以看出,计算塔线体系悬垂绝缘子串最大风偏角 的风振系数β的步骤为:S321:根据步骤S1中常规输电塔塔线体系的输电线、 绝缘子串的物理参数,以重力和平均风荷载作用下作为导线和悬垂绝缘子串计 算的初始条件,通过LRC方法确定悬垂绝缘子串风偏角的计算模型;具体见图 5.S322:计算常规输电塔塔线体系中输电塔之间导线单位面积的等效静力风荷 载;输电塔之间导线单位面积的等效静力风荷载pESWL的计算公式为:
Figure BDA0002433914740000198
式中,(:,i)表示矩阵的第i列元素;
Figure BDA0002433914740000199
为等效背景风压;
Figure BDA00024339147400001910
为平均风荷载;导线在风荷载作用下的振动方程矩阵表达式为:
Figure BDA0002433914740000201
式中,
Figure BDA0002433914740000202
Y′分别为脉动风荷载作用下导线节点顺风向的加速度、速度 和位移;
Figure BDA0002433914740000203
为平均风荷载作用下导线节点顺风向的位移;M为质量矩阵;C为 阻尼矩阵;K刚度矩阵;Ls为节点从属面积矩阵;导线在脉动风荷载作用下的 振动方程矩阵表达式为:
Figure BDA0002433914740000204
323:计算悬垂绝缘子串的风振 系数;
Figure BDA0002433914740000205
c表示对计算域内的元素进行求和;Γc为计算域内导线的线长;
Figure BDA0002433914740000206
为平均 风荷载;
Figure BDA0002433914740000207
为等效背景风压。本实施例中,DL/T 5154的导/地线的水平风荷载标 准值表达式为:
Figure BDA0002433914740000208
其中,β=α'βc
式中,α′为取值小于1的风压不均匀系数;μsc为阻力系数;βc为风荷载调整系 数,计算风偏角时取1;Dc为子导线/地线的计算外径;Lp为杆塔的水平档距;Bl为覆冰时风荷载的增大系数。
α′βc和荷载规范中β的物理意义一致,考虑了脉动的风动力效应。通过考虑风 压不均匀性的平均风荷载乘以βc,以此确定导/地线的等效静力风荷载。因此, α′βc=β。根据物理意义,
Figure BDA0002433914740000209
采用LRC计算的β不是常数,为方便设计使用, 根据pESWL的分布特性,采用平均化方式处理,计算一致β。pESWL在目标点位置凸出, 在远离目标点位置逼近于
Figure BDA00024339147400002010
为非均匀分布。为此,设定计算域,将目标点的等 效静力风荷载在计算域内进行平均化处理。当目标点与邻近杆塔导线挂点的高 差为0时,选取目标点水平档距为其计算域。有高差时,目标点位置处的等效 静力风荷载显得更加凸出,从而选取目标点左右1/4跨为其计算域。
S4:考虑塔线耦合效应,求取塔线体系风荷载脉动折减系数εc,结合图10 可以看出,具体为:S41:构建常规输电塔塔线体系计算模型,并得到塔线体系 计算模型图详见图1;所述塔线体系计算模型中的杆塔为密实结构,塔身为正方 形的变截面,由下至上尺寸变小,横担为等截面;所述塔线体系计算模型中的 导线两端等高,与固定铰支座连接;所述塔线体系计算模型中的杆塔高度为H, 横担悬臂长度为lca,绝缘子长度为lin,导线跨度为L。导线挂点无高差;
S42:根据常规输电塔,建立杆塔响应与杆塔风振系数的关系,得到塔高H 处建立杆塔荷载引起塔顶位移的均方根值σut(H)与杆塔风振系数β(H)的关系 式;当导线悬挂于杆塔顶部时,建立导线与导线风振系数的关系,得到导线荷 载引起塔顶位移的均方根值σuc(H)的计算公式;
Figure BDA0002433914740000211
其中,ω0为基本风压;μz(H)为风压随 常规单塔高度变化系数;μs(H)为杆塔随高度阻力系数;bs(H)随高度变化的迎 风面宽度;gs为峰值因子;ω1为顺风向1阶模态的自振圆频率;m(H)为随高度 变化的单位高度质量;所述导线荷载引起塔顶位移的均方根值σuc(H)的计算公式 为:
Figure BDA0002433914740000212
其中,Np为导线的相数;μsc为导线阻力系数;μz(Hc)为风压随导线高度变 化系数Nc为分裂导线的个数;Dc为子导线/地线的计算外径;Lp为水平档距;H 为塔高高度;Et为弹性模量。S43:根据步骤S41得到的内容,采用SRSS的方 法确定塔线体系下杆塔的峰值响应计算公式;
Figure BDA0002433914740000213
其中,所述
Figure BDA0002433914740000214
为由杆塔平均风荷载引起的杆塔响应;
Figure BDA0002433914740000215
为由导线平均风荷载引起 的杆塔响应;
Figure BDA0002433914740000216
为塔线体系平均风荷载引起的杆塔响应σr为塔线体系下杆塔响应 的标准差;gs为峰值因子;σrt为由杆塔脉动风荷载引起的塔体均方根响应;σrc为 由导线脉动风荷载引起的塔体均方根响应。S44:基于步骤S43的塔线体系下杆 塔的峰值响应计算公式,采用塔线分离方法,引入杆塔风荷载脉动折减系数, 并得到所述杆塔的峰值响应计算公式的等价峰值响应计算公式;
Figure BDA0002433914740000221
Figure BDA0002433914740000222
表示杆塔荷载引起的峰值响应,
Figure BDA0002433914740000224
表示 输电线荷载引起的峰值响应。S45:以塔顶位移响应为目标,对步骤S44得到的 所述杆塔的峰值响应计算公式的等价峰值响应计算公式进一步更新,得到带未 知导线荷载引起塔顶位移的均方根值和未知杆塔荷载引起塔顶位移的均方根值 的脉动折减系数更新计算公式;
S46:将步骤S42中的计算得到的导线荷载引起塔顶位移的均方根值和杆塔 荷载引起塔顶位移的均方根值带入步骤S45得到脉动折减系数的更新计算公式 中,得到脉动折减系数的最终计算公式,并计算杆塔风荷载脉动折减系数;计 算杆塔风荷载脉动折减系数的公式为:
Figure BDA0002433914740000225
其中,
Figure BDA0002433914740000226
在本实施例中,所述塔线体系常规输电塔的修正风振系数β*(z)和所述塔线 体系输电线的修正风振系数β*的计算公式为:
Figure BDA0002433914740000227
在等效振动惯性力作用下计算得到塔线体系中常规输电塔的设计风荷载 fESWL(z)与所述塔线体系常规输电塔的修正风振系数β*(z)的关系式存在:
Figure BDA0002433914740000228
其中,ξ1=ξe
Figure BDA0002433914740000229
Figure BDA00024339147400002210
Figure BDA0002433914740000231
m(z)=m(0)μm(z);
Figure BDA0002433914740000232
Sf(n)为归一化风速谱,
Figure BDA0002433914740000233
Iz(z)为z高度处的脉动风湍流密度;
Figure BDA0002433914740000234
I10为10m高度处的脉动风湍流密度;x′1为公式
Figure BDA0002433914740000235
中,n=n1时的取值,n1为输电塔的1阶模态频率;
u1和ηxz1是与风场湍流特性和空间相关性等有关的系数,分别称为综合影响 系数和空间相关性折减系数。
在本实施例中,基于塔线分离法计算输电线设计风荷载WX的计算公式为:
Figure BDA0002433914740000236
其中,β=α'βc;α′为取值小于1的风压不均匀系数;μsc为导线阻力系数;βc为 风荷载调整系数,计算风偏角时取1;Dc为子导线/地线的计算外径;Lp为杆塔 的水平档距;Bl为覆冰时风荷载的增大系数;ω0为基本风压;μz为风压随高度变 化系数;Bl为覆冰时风荷载的增大系数;Nc为分裂导线的个数;θ为风向角。
综上可以基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线 设计风载荷,对常规输电塔进行设计。
应当指出的是,上述说明并非是对本发明的限制,本发明也并不仅限于上 述举例,本技术领域的普通技术人员在本发明的实质范围内所做出的变化、改 性、添加或替换,也应属于本发明的保护范围。

Claims (8)

1.一种基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法,其特征在于:具体步骤为:
S1:搭建常规输电塔的塔线体系,并获取塔线体系的常规输电塔、输电线、绝缘子串的物理参数;
S2:基于塔线耦合影响因子,根据塔线体系中杆塔等效阻尼系数ζe
S3:将步骤S2得到的塔线体系中杆塔等效阻尼系数ζe来替换阻尼系数ζ1,求取塔线体系常规输电塔的风振系数β(z);
考虑线形与线长影响因子,计算塔线体系悬垂绝缘子串最大风偏角的风振系数β;
其中,步骤S3求取塔线体系常规输电塔的风振系数β(z)的步骤为:
S311:根据步骤S1中常规输电塔的物理参数,确定常规输电塔所在地面粗糙度类别,设定10m高度处的平均分速
Figure FDA0002433914730000011
常规输电塔的总高度H;跟开b1;横担个数nc;横担平均外伸长度
Figure FDA0002433914730000012
S312:构建常规输电塔的风荷载的计算模型,通过水平均布荷载作用下结构的挠曲线获得常规输电塔0°风向角的1阶侧弯振型φ1(z);
Figure FDA0002433914730000013
z为实际高度值
S313:根据荷载规范引入背景分量因子Bz(z),进而计算输电塔的脉动风荷载在水平方向的相关系数ρx和脉动风荷载在竖直方向的相关系数ρz
根据荷载规范引入的共振分量因子R;确定地面粗糙度指数α;峰值因子gs;10m高度处的湍流度I10
Figure FDA0002433914730000021
b为外轮廓宽度;共振分量因子R的计算公式为:
Figure FDA0002433914730000022
其中,ξ1=ζe
Figure FDA0002433914730000023
n为脉动风速的频率;
S314:根据地面粗糙度类别获取背景分量因子的中间变量γ的拟合系数kγ和aγ;根据常规输电塔的宽度、深度得到风振系数考虑整体外形变化的修正系数θv;根据常规输电塔的横担平均外伸长度、总高度、横担个数得到风振系数考虑附加面积的修正系数θa与风振系数考虑附加质量的修正系数θm的乘积θl
S315:简化并得到剩余塔身的风振系数考虑局部外形变化的修正系数的计算公式并计算得到对应的剩余塔身的风振系数考虑局部外形变化的修正系数的θb(z);
计算得到横担的风振系数考虑局部外形变化的修正系数θb(zI)和横隔面的风振系数考虑局部外形变化的修正系数θb(zJ);
其中,
Figure FDA0002433914730000024
S316:根据塔身的实际高度值z,结合步骤S315对应得到的风振系数考虑局部外形变化的修正系数计算公式计算所述背景分量因子Bz(z),具体计算公式为:
Figure FDA0002433914730000025
其中,
Figure FDA0002433914730000026
θν是风振系数考虑整体外形变化的修正系数;θη是风振系数考虑脉动风空间相关性的修正系数;
其中,
Figure FDA0002433914730000031
Figure FDA0002433914730000032
S317:计算风振系数β(z);
Figure FDA0002433914730000033
其中,gs为峰值因子;
S4:考虑塔线耦合效应,求取塔线体系风荷载脉动折减系数εc
S5:根据步骤S4得到的塔线体系风荷载脉动折减系数,对步骤S3中的塔线体系常规输电塔的风振系数、风偏角的风振系数β进行修正计算,得到塔线体系常规输电塔的修正风振系数β*(z)和塔线体系输电线的修正风振系数β*
S6:基于塔线分离法,在等效振动惯性力作用下计算塔线体系中常规输电塔的设计风荷载fESWL(z)和输电线的设计风荷载WX
2.根据权利要求1所述的基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法,其特征在于:步骤S2的具体步骤为:
S21:根据步骤S1的常规输电塔的塔线体系,得到常规输电塔塔线体系计算模型图;
所述塔线体系计算模型中的杆塔为密实结构,塔身为正方形的变截面,由下至上尺寸变小,横担为等截面;所述塔线体系计算模型中的导线两端等高,与固定铰支座连接;所述塔线体系计算模型中的杆塔高度为H,横担悬臂长度为lca,绝缘子长度为lin,导线跨度为L;导线挂点无高差
S22:设定输电线和绝缘子串振动的假设条件,得到的常规输电塔塔线体系中输电线和绝缘子串的振型图以及迎风面、被风面输电线和绝缘子串的广义质量、广义刚度和广义阻尼;并将常规输电塔塔线体系中输电线和绝缘子串组合形成索结构体系;
所述迎风面、被风面导线的广义质量计算公式为:
Figure FDA0002433914730000041
所述迎风面、被风面导线的广义刚度计算公式为:
Figure FDA0002433914730000042
所述迎风面、被风面导线的广义阻尼计算公式为:
Figure FDA0002433914730000043
mc为单根导线单位线长的质量;单根导线振型
Figure FDA0002433914730000044
γg为导线的自重比载;σ0为导线的水平初应力;Γ为导线的线长,
Figure FDA0002433914730000045
ζc=ζscac;ζsc为导线结构阻尼比;ζac为导线启动阻尼比;Nc为分裂导线的个数;Tw为平均风状态下单根导线的水平张力;ζc为导线阻尼比;
所述迎风面、被风面绝缘子串的广义质量计算公式为:
Figure FDA0002433914730000046
所述迎风面、被风面绝缘子串的广义刚度计算公式为:
Figure FDA0002433914730000047
所述迎风面、被风面绝缘子串的广义阻尼计算公式为:
Figure FDA0002433914730000048
其中,min为绝缘子串单位高度质量;Din为绝缘子串迎风外径;绝缘子串振型
Figure FDA0002433914730000051
H-lin≤z≤H;ζin为绝缘子串阻尼比;
Figure FDA0002433914730000052
lin为绝缘子长度;其中,索结构体系对应的广义质量、广义刚度和广义阻尼的计算公式为:
Figure FDA0002433914730000053
S23:将步骤S22得到的数据构建索结构体系结合杆塔结构组成塔线耦合简化计算模型;
S24:基于杆塔结构组成塔线耦合简化计算模型,求取常规输电塔塔线体系下杆塔顺风向位移的共振分量的均方值和单塔时杆塔顺风向位移共振分量的均方值;从而得到二者的比例式;
S25:基于步骤S24得到的计算公式,推导常规输电塔悬挂输电线后塔线等效阻尼系数的计算公式,并计算塔线体系中杆塔等效阻尼系数。
3.根据权利要求2所述的基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法,其特征在于:步骤S24中所述塔线体系下杆塔顺风向位移的共振分量的均方值的计算公式为:
Figure FDA0002433914730000054
其中,
Figure FDA0002433914730000055
Figure FDA0002433914730000056
λn=nci/nt
Figure FDA0002433914730000061
Figure FDA0002433914730000062
Figure FDA0002433914730000063
Figure FDA0002433914730000064
杆塔1阶模态的振型φt(z)=(z/H)2,0≤z≤H;
Figure FDA0002433914730000065
为杆塔的广义质量,
Figure FDA0002433914730000066
Mca为横担的质量,mt(z)为随高度变化的杆塔单位高度质量;
coh(z1,z2)为z1和z2高度处两点的脉动风速的相干函数;
Sf(nt)为归一化风速谱,nt为杆塔脉动风速的频率;
Figure FDA0002433914730000067
σv'为脉动风速的标准差;
Figure FDA0002433914730000068
为索结构与杆塔的广义质量比值,
Figure FDA0002433914730000069
λn为索结构与杆塔的频率比值;λn=nci/nt;导线悬挂于杆塔的顶部,
Figure FDA00024339147300000610
ζt为总阻尼比;ζt=ζstat;ζst为杆塔结构阻尼比;ωt为杆塔无阻尼振动的圆频率;ζci为索结构总阻尼比,近视取导线阻尼比,ζci≈ζc,ζc=ζscac
ρa为空气密度,μs(z)为风压随高度变化系数;bs(z)为随高度变化的迎风面宽度;
Figure FDA00024339147300000611
随高度变化的平均风速,σv'为脉动风速的标准差;
ζat为杆塔气动阻尼比;
Figure FDA00024339147300000612
As,ca为横担的挡风面积;
所述单塔时杆塔顺风向位移共振分量的均方值为:
Figure FDA0002433914730000071
其中,
Figure FDA0002433914730000072
所述塔线体系下杆塔顺风向位移的共振分量和所述单塔时杆塔顺风向位移共振分量的比例式为:
Figure FDA0002433914730000073
步骤S25中推导常规输电塔悬挂导线后塔线等效阻尼系数的计算公式的步骤为:
悬挂导线后杆塔的等效阻尼比为:
Figure FDA0002433914730000074
其中,ρ与ζe的关系为:
Figure FDA0002433914730000075
对于输电塔线体系而言,索结构为柔性体系,卓越频率远小于杆塔的频率;则忽略λn的高阶项;杆塔的阻尼比约为0.01,索结构的阻尼比小于1,则忽略
Figure FDA0002433914730000076
项;
导线悬挂于杆塔的顶部,
Figure FDA0002433914730000077
故悬挂导线后塔线等效阻尼系数的计算公式为:
Figure FDA0002433914730000078
4.根据权利要求1或3所述的基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法,其特征在于:塔线体系悬垂绝缘子串最大风偏角的风振系数β的步骤为:
S321:根据步骤S1中常规输电塔塔线体系的输电线、绝缘子串的物理参数,以重力和平均风荷载作用下作为导线和悬垂绝缘子串计算的初始条件,通过LRC方法确定悬垂绝缘子串风偏角的计算模型;
S322:计算常规输电塔塔线体系中输电塔之间导线单位面积的等效静力风荷载;
所述输电塔之间导线单位面积的等效静力风荷载pESWL的计算公式为:
Figure FDA0002433914730000081
式中,(:,i)表示矩阵的第i列元素;
Figure FDA0002433914730000082
为等效背景风压;
Figure FDA0002433914730000083
为平均风荷载;
导线在风荷载作用下的振动方程矩阵表达式为:
Figure FDA0002433914730000084
式中,
Figure FDA0002433914730000085
Y′分别为脉动风荷载作用下导线节点顺风向的加速度、速度和位移;
Figure FDA0002433914730000086
为平均风荷载作用下导线节点顺风向的位移;
M为质量矩阵;C为阻尼矩阵;K刚度矩阵;Ls为节点从属面积矩阵;
导线在脉动风荷载作用下的振动方程矩阵表达式为:
Figure FDA0002433914730000087
S323:计算悬垂绝缘子串的风振系数;
Figure FDA0002433914730000088
Figure FDA0002433914730000089
Figure FDA00024339147300000810
c表示对计算域内的元素进行求和;Γc为计算域内导线的线长;
Figure FDA00024339147300000811
为平均风荷载;
Figure FDA00024339147300000812
为等效背景风压。
5.根据权利要求1所述的基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法,其特征在于:塔线体系风荷载脉动折减系数εc的计算步骤为:
S41:构建常规输电塔塔线体系计算模型,并得到塔线体系计算模型图;
S42:根据常规输电塔,建立杆塔响应与杆塔风振系数的关系,得到塔高H处建立杆塔荷载引起塔顶位移的均方根值σut(H)与杆塔风振系数β(H)的关系式;
当导线悬挂于杆塔顶部时,建立导线与导线风振系数的关系,得到导线荷载引起塔顶位移的均方根值σuc(H)的计算公式;
S43:根据步骤S41得到的内容,采用SRSS的方法确定塔线体系下杆塔的峰值响应计算公式;
S44:基于步骤S43的塔线体系下杆塔的峰值响应计算公式,采用塔线分离方法,引入杆塔风荷载脉动折减系数,并得到所述杆塔的峰值响应计算公式的等价峰值响应计算公式;
S45:以塔顶位移响应为目标,对步骤S44得到的所述杆塔的峰值响应计算公式的等价峰值响应计算公式进一步更新,得到带未知导线荷载引起塔顶位移的均方根值和未知杆塔荷载引起塔顶位移的均方根值的脉动折减系数更新计算公式;
S46:将步骤S42中的计算得到的导线荷载引起塔顶位移的均方根值和杆塔荷载引起塔顶位移的均方根值带入步骤S45得到脉动折减系数的更新计算公式中,得到脉动折减系数的最终计算公式,并计算杆塔风荷载脉动折减系数;
Figure FDA0002433914730000091
其中,
Figure FDA0002433914730000101
6.根据权利要求1所述的基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法,其特征在于:
所述塔线体系常规输电塔的修正风振系数β*(z)和所述塔线体系输电线的修正风振系数β*的计算公式为:
Figure FDA0002433914730000102
7.根据权利要求1所述的基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法,其特征在于:在等效振动惯性力作用下计算得到塔线体系中常规输电塔的设计风荷载fESWL(z)与所述塔线体系常规输电塔的修正风振系数β*(z)的关系式存在:
所述常规输电塔的塔线体系中常规输电塔设计风荷载fESWL(z)与所述塔线体系常规输电塔的修正风振系数β*(z)的关系式为:
Figure FDA0002433914730000103
其中,ξ1=ξe
Figure FDA0002433914730000104
Figure FDA0002433914730000105
Figure FDA0002433914730000106
Figure FDA0002433914730000111
m(z)=m(0)μm(z);
Figure FDA0002433914730000112
Sf(n)为归一化风速谱,
Figure FDA0002433914730000113
Iz(z)为z高度处的脉动风湍流密度;
Figure FDA0002433914730000114
I10为10m高度处的脉动风湍流密度;x′1为公式
Figure FDA0002433914730000115
中,n=n1时的取值,n1为输电塔的1阶模态频率;
u1和ηxz1是与风场湍流特性和空间相关性等有关的系数,分别称为综合影响系数和空间相关性折减系数。
8.根据权利要求1所述的基于惯性力法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法,其特征在于:基于塔线分离法计算输电线设计风荷载WX的计算公式为:
Figure FDA0002433914730000116
其中,β=α'βc;α′为取值小于1的风压不均匀系数;μsc为导线阻力系数;βc为风荷载调整系数,计算风偏角时取1;Dc为子导线/地线的计算外径;Lp为杆塔的水平档距;Bl为覆冰时风荷载的增大系数;ω0为基本风压;μz为风压随高度变化系数;Bl为覆冰时风荷载的增大系数;Nc为分裂导线的个数;θ为风向角。
CN202010245662.0A 2020-03-31 2020-03-31 考虑塔线耦合影响的常规输电线路设计风载荷计算方法 Active CN111651804B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010245662.0A CN111651804B (zh) 2020-03-31 2020-03-31 考虑塔线耦合影响的常规输电线路设计风载荷计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010245662.0A CN111651804B (zh) 2020-03-31 2020-03-31 考虑塔线耦合影响的常规输电线路设计风载荷计算方法

Publications (2)

Publication Number Publication Date
CN111651804A true CN111651804A (zh) 2020-09-11
CN111651804B CN111651804B (zh) 2022-08-26

Family

ID=72346435

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010245662.0A Active CN111651804B (zh) 2020-03-31 2020-03-31 考虑塔线耦合影响的常规输电线路设计风载荷计算方法

Country Status (1)

Country Link
CN (1) CN111651804B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113032882A (zh) * 2021-03-31 2021-06-25 广东电网有限责任公司 风振响应模型生成方法、设备、装置及介质
CN115048694A (zh) * 2022-06-09 2022-09-13 国网山东省电力公司临沂供电公司 输电塔体系振型聚类方法、装置以及计算机设备
CN116108725A (zh) * 2023-04-10 2023-05-12 国网浙江省电力有限公司温州供电公司 输电线路铁塔中预应力拉索的参数优化方法、装置及介质

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101692566A (zh) * 2009-09-09 2010-04-07 中国电力科学研究院 一种输电线路高塔风振控制方法
CN105740548A (zh) * 2016-02-01 2016-07-06 西安交通大学 一种随机风载荷下输电线路风振计算方法
CN106326584A (zh) * 2016-08-29 2017-01-11 国网江西省电力公司电力科学研究院 一种输电线路风‑冰耦合荷载断线效应的分析方法
JP2017053275A (ja) * 2015-09-10 2017-03-16 株式会社日立製作所 風力発電システムまたは風力発電システムの制御方法
CN108073758A (zh) * 2017-10-17 2018-05-25 广东电网有限责任公司电力科学研究院 一种输电塔线风致振动响应仿真方法及装置
CN108121893A (zh) * 2017-12-18 2018-06-05 广西电网有限责任公司电力科学研究院 一种基于塔线耦合的配电杆塔抗风性能评估方法
CN108710763A (zh) * 2018-05-22 2018-10-26 国网江西省电力有限公司经济技术研究院 220kV输电线路塔线耦联体系覆冰仿真方法
CN108959742A (zh) * 2018-06-20 2018-12-07 重庆科技学院 大跨越输电塔线体系气动弹性模型设计方法
CN110298101A (zh) * 2019-06-24 2019-10-01 国网浙江省电力有限公司电力科学研究院 一种耦合局地风场的输电线路风致响应有限元分析方法
CN110378050A (zh) * 2019-07-25 2019-10-25 重庆科技学院 一种输电塔顺风向气动阻尼比计算方法
CN110851905A (zh) * 2019-11-14 2020-02-28 福州大学 一种峡谷微地形条件下的输电铁塔风振系数计算方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101692566A (zh) * 2009-09-09 2010-04-07 中国电力科学研究院 一种输电线路高塔风振控制方法
JP2017053275A (ja) * 2015-09-10 2017-03-16 株式会社日立製作所 風力発電システムまたは風力発電システムの制御方法
CN105740548A (zh) * 2016-02-01 2016-07-06 西安交通大学 一种随机风载荷下输电线路风振计算方法
CN106326584A (zh) * 2016-08-29 2017-01-11 国网江西省电力公司电力科学研究院 一种输电线路风‑冰耦合荷载断线效应的分析方法
CN108073758A (zh) * 2017-10-17 2018-05-25 广东电网有限责任公司电力科学研究院 一种输电塔线风致振动响应仿真方法及装置
CN108121893A (zh) * 2017-12-18 2018-06-05 广西电网有限责任公司电力科学研究院 一种基于塔线耦合的配电杆塔抗风性能评估方法
CN108710763A (zh) * 2018-05-22 2018-10-26 国网江西省电力有限公司经济技术研究院 220kV输电线路塔线耦联体系覆冰仿真方法
CN108959742A (zh) * 2018-06-20 2018-12-07 重庆科技学院 大跨越输电塔线体系气动弹性模型设计方法
CN110298101A (zh) * 2019-06-24 2019-10-01 国网浙江省电力有限公司电力科学研究院 一种耦合局地风场的输电线路风致响应有限元分析方法
CN110378050A (zh) * 2019-07-25 2019-10-25 重庆科技学院 一种输电塔顺风向气动阻尼比计算方法
CN110851905A (zh) * 2019-11-14 2020-02-28 福州大学 一种峡谷微地形条件下的输电铁塔风振系数计算方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
D.M. HUANG等: "Covariance proper transformation-based pseudo excitation algorithm and simplified SRSS method for the response of high-rise building subject to wind-induced multi-excitation", 《ENGINEERING STRUCTURES》 *
仲崇硕: "地表变形作用下输电塔线体系的风振响应和承载性能研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
刘俊卿 等: "双呼高输电塔的风振响应及风振系数研究", 《西安建筑科技大学学报(自然科学版)》 *
赵爽 等: "基于风洞试验的苏通大跨越输电塔风振系数研究", 《建筑结构学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113032882A (zh) * 2021-03-31 2021-06-25 广东电网有限责任公司 风振响应模型生成方法、设备、装置及介质
CN115048694A (zh) * 2022-06-09 2022-09-13 国网山东省电力公司临沂供电公司 输电塔体系振型聚类方法、装置以及计算机设备
CN115048694B (zh) * 2022-06-09 2023-05-16 国网山东省电力公司临沂供电公司 输电塔体系振型聚类方法、装置以及计算机设备
CN116108725A (zh) * 2023-04-10 2023-05-12 国网浙江省电力有限公司温州供电公司 输电线路铁塔中预应力拉索的参数优化方法、装置及介质
CN116108725B (zh) * 2023-04-10 2023-07-14 国网浙江省电力有限公司温州供电公司 输电线路铁塔中预应力拉索的参数优化方法、装置及介质

Also Published As

Publication number Publication date
CN111651804B (zh) 2022-08-26

Similar Documents

Publication Publication Date Title
CN111506946B (zh) 基于有效荷载法和塔线分离法考虑塔线耦合影响的常规输电塔、线设计风载荷的计算方法
CN111651804B (zh) 考虑塔线耦合影响的常规输电线路设计风载荷计算方法
CN111985018B (zh) 基于惯性力法和塔线分离法考虑塔线耦合影响的超高大跨越塔、线风载荷的计算方法
CN108959742B (zh) 大跨越输电塔线体系气动弹性模型设计方法
CN112287424B (zh) 基于有效荷载法和塔线分离法考虑塔线耦合影响的超高大跨越塔、线设计风载荷的计算方法
CN111651805B (zh) 考虑线形与线长影响因子的悬垂绝缘子串最大风偏角及风振系数计算方法
CN106096105A (zh) 输电线路风偏瞬态响应计算方法
CN111506944B (zh) 考虑塔线耦合效应的杆塔风荷载脉动折减系数计算方法
CN109902351A (zh) 一种覆冰导线动态风偏的简化计算方法
CN111783201A (zh) 一种三跨自锚式悬索桥动力特性的快速分析方法
CN111651806B (zh) 基于惯性力法的常规输电塔设计风载荷计算方法
CN109271751A (zh) 一种悬垂绝缘子串的最大动态风偏角确定方法
CN111783198A (zh) 一种基于双梁模型的悬索桥精细化动力分析方法
CN111651907B (zh) 一种复杂索网体系的建模分析方法
CN107143463B (zh) 一种控制风力机叶片挥舞的空间拉索系统及优化方法
CN117271949B (zh) 考虑弹性边界和大垂度影响的悬索自振分析方法及系统
CN111783199A (zh) 一种多段式索缆结构动力特性的精细化快速求解方法
CN104636603B (zh) 一种特高压单柱拉线塔扭振频率的两自由度计算方法
CN111651920B (zh) 塔线体系下大跨越输电塔风致机械能的传递分析方法
CN111985019B (zh) 带悬挑横担常规单塔基于有效荷载法的风振系数简化计算方法
CN111506945B (zh) 基于塔线耦合影响因子的输电塔共振响应的等效阻尼系数计算方法
CN111241731B (zh) 一种依据水平张力的输电导线找形方法
CN114912324A (zh) 一种基于分布式mtmdi的大跨度桥梁多阶涡激振动控制方法
CN114329926A (zh) 一种计算架空输电线路覆冰不平衡张力的方法
CN111651808A (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