3-第二章 风振系数计算

第2章 风振系数计算

2.1 引言

在随机脉动风压的作用下，高耸结构会产生随机振动，除了顺风向的风振响应外，结构还会产生横风向的风振响应。但在通常情况下，对于非圆截面，顺风向风振响应占据主要地位，对于一般的塔架结构，可以忽略横风向共振的作用[13]。因此，本章主要研究输电塔结构在随机风荷载作用下的顺风向风振系数的计算。

作用于结构物上的脉动风荷载对结构产生的动力响应与结构物本身的动力特性有关。当结构物刚性很强时，由脉动风所引起的结构物风振惯性力并不明显，可以略去，但需要考虑由脉动风所引起的瞬时阵风荷载；当结构物刚性较弱即为柔性结构时，除静力风荷载(z)外，还应计及风振惯性力的大小，即风振动力荷载。如果风振动力荷载用d(z,t)表示，则柔性结构物的总风荷载W(z,t)表达如下[4]：

W(z,t)(z)d(z,t) （2－1）

工程计算中，常采用集中风荷载的表达式，则式（2－1）改写为

P(z)Pc(z)Pd(z) （2－2a）

或

P d （2－2b） PiPci式中，P(z)，Pi —— 顺风向z高度处第i点的总风荷载（kN）；

Pc(z)，Pci —— 顺风向z高度处i点总静力风荷载（kN）； Pd(z)，Pdi——顺风向z高度处i点风振动力荷载（kN），其中

Pd(z)d(z)Az，()id()zA或PidziAz(Ai)。在这里，

为z高度（第i点）处相关的迎风面竖向投影面积（m2）。

本章下面将讨论风振动力荷载的计算原理和表达式，以及可在实际输电塔设计中应用的风振系数的计算方法。

2.2 顺风向风振系数的计算方法

2.2.1结构风振随机振动理论[4][10][7]

风荷载是输电塔结构的各类荷载中起主要作用的荷载，由静、动两部分风荷载组成，动力风荷载即脉动风是一种随机动力干扰，引起结构的振动。为了便于工程的实际应用，我国的《建筑结构荷载规范》引入了风振系数作为等效静态放大系数，将风荷载的静力作用与动力作用一并考虑在内。

在风作用下，结构可在各个方向产生振动。通常结构抗风验算可在结构两个主轴进行。当验算的主轴方向与风向一致时，结构发生顺风向响应。由于风可分解为平均风和脉动风，前者变化缓慢，周期很长，可作为静态作用来处理，而后者变化很快，周期很短，引起结构激烈振动。由于顺风向脉动风的作用是随机的，引起了结构随机振动[10]。

基于结构随机振动理论，可导出风振时各类风响应。工程中将风敏感的高耸结构等结构设为竖向一维悬臂结构，它们的风振响应以第一振型为主，常以先求风振系数再求各类风响应较为方便。

本文按照张相庭教授提出的由基本理论推出的通用计算公式计算输电高塔的风振系数，现行荷载规范中的计算公式是以结构物质量和刚度线性分布为前提对该通式的简化，并不适用于输电塔这类外形变化不规则及附有集中质量的结构[10][7]。

顺风向风振分析应按结构随机振动理论进行。对于任意一个n维自由度结构，采用矩阵表示的运动方程为

My(t)Cy(t)Ky(t)Pf(t) （2－3）

式中：M，C，K —— 分别是结构的质量、阻尼和刚度矩阵；

y(t)，y(t)，y(t) —— 分别为结构节点的的加速度、速度和位移向量；

Pf(t) —— 脉动风力向量。

令：

yi(t)jiqj(t) （2－4）

j1n式中：ji —— 振型j点i的振型系数；

qj(t) —— 振型j的广义坐标。

在取瑞雷阻尼符合正交性条件下，将式（2－4）代入式（2－3）可得振型广义坐标表示的形如单自由度的运动方程如下：

*Pfjf(t)qj(t)2jjq(t)q(t)Fj(t)2jM*j (j1,2,n, （2－5）

式中：j，j —— 分别为振型j的固有频率和阻尼；

Fj(t) —— 脉动风动力作用的振型j广义力；

* —— 将时间分量分离后，脉动风作为静力时振型j的广义力，常简Pfj*fiAiji； 称为脉动风的振型j广义力，Pfjifi —— 结构上点i的脉动面力；

Ai —— 点i的迎风面积；

f(t) —— 脉动风的时间函数；

*2M*j —— 振型j的广义质量；MjMiji，Mi为点i的团集质量。

i由于在风力输入时，脉动风的时间函数f(t)包含有随机性，因而需要根据随机振动理论来求解式（2－5）。此时风力输入为统计值，常以输入功率谱密度

Sf()为代表，为圆频率，由于脉动风具有空间相关性，因而不同点i和i之

间风压应考虑空间相关性系数ii()。输出亦为统计值，常以输出位移的功率谱密度Sy()为代表。由此，根据随机振动理论由式（2－5）可求出点i位移响应的根方差yi：

yinSyi()d2ji*2jj1n2yji2jiSFF()Hj(i)d j1jjn2 =(Mj1)Sf()Hj(i)2iifiAifiAiii()jijid （2－6）

式中：SFjFj() —— 振型j脉动风动力作用的自功率谱密度；

Hj(i) —— 振型j频率响应函数。

此时由脉动风引起的振型j点i的风振力即等效惯性力为：

PjiMi2jyji （2－7）

由式（2－7），就可以求出结构振型j点i的任何响应ri，如弯矩、剪力等的根方差rji。

在结构频率比较稀疏可以略去不同振型之间的相互影响时，结构总响应（包括位移）的根方差可由下式求出：

rij1n2rji （2－8）

由于根方差仅代表均方振幅，所以必须考虑一定的保证系数或称为峰值因子（简称峰因子）（Peak factor）才可以求出在一定保证率下的响应幅值，峰因子用表示[4]。

如果总响应最大值出现的概率和各振型响应最大值出现的概率都相同，则

ri(rji)2j1n2rji （2－9） j1n对于高耸结构等第一振型起主要作用的结构，式（2－9）可简化为：

yi1iM1*Sf()H1(i)2iifiAifiAiii()1i1id （2－10）

此时，采用风振系数进行结构抗风分析较为方便，结构点i的风荷载为：

wiiwsiisiziw0 （2－11）

式中：w0 —— 基本风压；

si —— 点i的体型系数；

wsi —— 平均风压；

i —— 点i的风振系数；

zi —— 点i的风压高度变化系数，与地面粗糙程度及高度变化有关。

z0.32zz35 （2－12）

HT2式中：、HT —— 分别为任意地貌的粗糙程度系数及梯度风高度；

于是，可以得到风振系数i的计算式：

iPsiPdfiPsi*11Mi12yisizi0Ai

2 = 1Mi121iMsizi0AiSf()H1(i)iifiAifiAiii()1i1id

（2－13）

式（2－13）即是风振系数的基本表达式。我国许多现行规范中的公式均为该式的变形和简化，本文也是根据该式推出程序中应用的表达式。

2.2.2 输电塔设计中可应用的风振系数表达式

根据风振系数基本计算公式（2－13），即可进行风振系数的分析计算。但要用该公式作风振系数的计算就需要求出脉动风动力作用的功率谱密度Sf()与频率响应函数Hj(i)，这使得该公式很难在实际工程设计中直接使用。因此，为了应用更为方便，一些更重视实际工程应用的论著以及各国规范提出了多种更简单的表达形式。

本文中采用以脉动增大系数来反映脉动风主要动力特性的表达形式。将式（2－13）按我国荷载规范可改写成

Mi1i12siziAii1Sf()H1(i)d2

fiAifiAiii()1i1id2Sf()H1(i)ii02M*1Sf()H1(i)d2

11u1r1i （2－14） 式中：1，u1，r1i —— 分别称为第一振型脉动增大系数、脉动影响系数和脉

动补充系数。

下面分别介绍这三个系数。 1．脉动增大系数1

在式（2－14）中，第一振型脉动增大系数1为：

112Sf()H1(i)d （2－15）

22上式中，传递函数的模H1(i)为 H1(i)21ff121f1f122412 （2－16）

根据随机振动理论，式（2－15）可以积分成三部分，但起主要作用的是前面两项，即常称的背景部分和共振部分。我国规范对风的频谱采用世界各国规范通常采用的Davenport实测拟合的脉动风功率谱：

Sf4KV102x2f1x24 （2－17）

3式中： —— 表面阻力系数，与地表特征有关；

V10 —— 10m高度处的平均风速（m/s）；

x1200f，f为频率(Hz)； V10从而可以得到：

1f1x1211Sf(f1)1 （2－18） 2243(21)61(1x1)式（2－15）到式（2－18）中，f为频率（2f）。当ff1（f1为结构第一自振频率）时，结构产生共振，相应有关部分也就称为共振部分；1为

2v10130结构第一振型阻尼比；f1，x1，w0。

2T11600w0T1式（2－18）表明，第一振型脉动增大系数1此时仅与1及w0T12有关，在我国《建筑结构荷载规范》（GB50009－2001）中，已将1按式（2－18）制成表格，供直接查用，如表2.1所示。

表2.1 脉动增大系数

w0T12(kNs2m2) 0.01（一般钢结构） 0.01 0.02 0.04 0.06 0.08 0.10 1.47 1.57 1.69 1.77 1.83 1.88 1.26 1.32 1.39 1.44 1.47 1.50 0.20 2.04 1.61 1.28 0.40 2.24 1.73 1.34 0.60 2.36 1.81 1.38  0.02（有填充墙的房屋钢结构） 0.05（混凝土及砌体结构） 1.11 1.14 1.17 1.19 1.21 1.23 w0T12(kNs2m2) 0.01（一般钢结构） 0.80 1.00 2.00 4.00 6.00 8.00 10.00 20.00 30.00 2.46 2.53 2.80 3.09 3.28 3.42 1.88 1.93 2.10 2.30 2.43 2.52 3.54 2.60 1.82 3.91 2.85 1.96 4.14 3.01 2.06  0.02（有填充墙的房屋钢结构） 0.05（混凝土及砌体结构） 1.42 1.44 1.54 1.65 1.72 1.77 注： 1. w0为该地区10m高的实际风压；

2. 计算w0T1时，对地面粗糙度B类地区可直接代入基本风压，而对A类、C类

和D类地区应该按当地的基本风压分别乘以1.38、0.62和0.32后代入。

22．脉动影响系数u1

在式（2－14）中，第一振型脉动影响系数u1为

u1Sf()H1(i)M1*2fiAifiAiii()1i1id20iiSf()H1(i)d2 （2－19）

其中，脉动风压fi可用脉动系数fi来表示，表达式为：

fifi （2－20）

sizi0建议采用： fz0.5351.80.16z （2－21） 10根据强风的观察资料表明，一次阵风作用在结构迎风面上各点处的风速和风向并不是完全同步的，有时甚至是几乎无关的。为此，对于构筑物上受到脉动风的作用是要考虑其空间相关性的。像输电塔这类高耸结构，高度远大于宽度，一般认为结构前后的脉动风压相同。本文讨论输电塔结构的顺风向风振，因此主要考虑风荷载的竖向相关性。根据文献[4]，脉动风压的竖向相关性可表达为：

zz （2－22） Lzzz,zexp式中：zz,z —— 脉动风压垂直方向的相干系数；

zz —— 垂直二点间距离。

Lz大部分落在40～60之间，根据试验资料的统计分析可知，故建议取Lz=50

米。

用考虑风压空间相关性的折算系数z1来表示输电塔受到风荷载的作用时脉动风压的竖向相关性，它只与结构物的总高度H有关。z1已制作成计算用表，如表2.2所示。

表2.2 迎风面风压空间相关性的折算系数z1

总高H（m） 10 0.98 20 0.97 30 0.95 40 0.93 50 0.92 60 0.90 70 0.89 80 0.87 z1 总高H（m） 90 0.86 100 0.85 150 0.79 200 0.74 250 0.70 300 0.67 350 0.64 400 0.61 z1 因而，对于连续分布质量体系或者集团质量体系，可得

u1fisiziAi1ifisiziAiii1ifisiziAi1iiii1M1*nMi12ii1nz1 （2－23）

3．脉动补充系数r1

在式（2－14）中，第一振型脉动补充系数r1为：

r1iMi1i （2－24）

siziAi式中各个系数的含义同前。

综合以上三式（2－18），（2－23）、（2－24），可以得到：

ni11u1r1i11w0T12fisiziAi1ii1Mi12ii1nz1HMi1i （2－25）

siziAi由式（2－25）可以知道，只要取得结构第一自振周期T1、基本风压w0、集团质量Mi、集团挡风面积Ai、第一振型系数1i、体型系数si、高度系数zi的值，就可以分别求出脉动增大系数1、脉动影响系数u1、脉动补充系数r1，进

而求得风振系数zi的值。

2.3 我国现行规范中风振系数的计算方法

2.3.1 《建筑结构荷载规范》中风振系数的表达式[2]

z1Bvz （2－26） z式中： —— 脉动增大系数，见表2.1；

—— 脉动影响系数，当结构迎风面和侧风面的宽度沿高度按直线或接近

直线变化，而质量沿高度按连续规律变化时，脉动影响系数应再乘以修正系数B、v，若结构的外形、质量沿高度比较均匀则不需要乘；

—— 构筑物迎风面在z高度处的宽度Bz与底部宽度B0的比值BzB0；

v —— 脉动影响系数按构筑物迎风面顶部的宽度BH与底部宽度B0的

比值BHB0的修正的系数；

z—— 振型系数，荷载规范提供了截面沿高度规律变化的高耸结构第一

振型系数的近似值；

z—— 风压高度变化系数；

在式（2－26）中，系数、v、z都是以质量和外形沿高度方向均匀变化的高耸结构为模型简化式（2－14）而得到的。

然而在高耸结构中，在像输电塔结构这样的沿高度方向质量和外形不规则变化的情况下，由于其质量和外形很难找到统一的规律性，因此由规范的表格或规定其实无法包括这类结构的风荷载计算。

2.3.2 《架空送电线路杆塔结构设计技术规定》中风振系数的规定[1]

对杆塔本身，当全高不超过 60m时，按照表2.3规定选用，全高采用一个系数；当杆塔全高超过60m时，应按《建筑结构荷载规范》[2]采用由下到上逐段增大的数值，但其加权平均值对自立式铁塔不应小于1.6。

表2.3 自立式杆塔风荷载调整系数z

杆塔全高H（m） 20 1.0 30 1.25 40 1.35 50 1.5 60 1.6 z 注：1. 中间值按插入法计算；

2. 对自立式铁塔，表中数值适用于高度与根开之比为4－6。

2.3.3 《架空送电线路大跨越设计技术规定》中风振系数的规定[17]

对于总高度小于50m的塔按《杆塔技术规定》[1]取用；对于总高度在50m

以上的塔按《荷载规范》[2]计算。

2.3.4《高耸结构设计规范》中风振系数的规定[3]

我国《高耸结构设计规范》中的风振系数计算公式实际上为《建筑结构荷载规范》（GB5009－2001）中的公式，即式（2－25）进行参数分析和约化的结果，必要时某些参数用近似式代入，则计算可以得到很大的简化。

高耸结构的迎风向宽度很小，如果结构物迎风面和侧风面的宽度沿高度按直线或接近直线变化，而质量沿高度按连续规律变化时，可以只选择高度作为变数。将式（2－25）各参数所依赖的变数写在该参数后面得到[4]：

zlx(H),,eHlx(0)1z11(1,w0T12)(,H)z(,z)1l(H)lx(z)xB lx(0)lx(0)l(H)lx(z)xB lx(0)lx(0)11(1,w0T12)(,H)z(,H)zlx(H),,eHl(0)xzH2zl(H)l(z)111,w0T121,H2,x,e,,x

lx(0)Hlx(0)1112 （2－26）

式中： —— 脉动增大系数；

1 —— 风压脉动和风压高度等的影响系数，它只涉及和H两个参数，

规范已经制作为表格可以根据参数来查取；

2—— 振型、结构外形的影响系数，规范中也已经制作为表格供查用。 式（2－26）即为《高耸结构设计规范》中规定的自立式高耸结构在z高度处的风振系数z的表达式。

与式（2－25）类似，式（2－26）中的系数1、2也同样是以质量和外形沿高度方向均匀变化的高耸结构为模型简化式（2－14）而得到的，因此同样不适合用于输电塔结构的风振系数计算。

2.4 小结

除了自重和导地线荷载之外，风荷载是输电塔结构的最重要的荷载，在随机脉动风压的作用下，塔架结构产生随机振动。本章主要研究输电塔结构在随机风荷载作用下的顺风向风振系数的计算式，得到的结论如下：

（1） 基于结构随机振动理论推导出风振系数的基本表达式可以化为能用

于集团质量体系的形式，只要知道结构第一自振周期T1、基本风压

w0、集团质量Mi、集团挡风面积Ai、第一振型系数1i、体型系数

si、高度系数zi的值，就可以求出风振系数zi的值。

（2） 我国《建筑结构荷载规范》和《高耸结构设计规范》中风振系数的

计算公式都是以质量和外形沿高度方向均匀变化的高耸结构为模型简化而得到的，但输电塔结构沿高度方向布置有数个横担结构，导致其质量和外形不规则变化，在横担处有突变。因此，理论上荷载规范和高耸规范的计算公式并不适用于输电塔结构。

（3） 根据随机振动理论，将输电塔结构简化为一个集团质量体系，由风

振系数的基本表达式（2－25）就可以较为准确地计算出输电塔结构的风振系数zi的值。