等效风荷载计算方法

等效静力风荷载的物理意义从风洞试验获取屋面风荷载气动力信息，到得到结构的风振响应整个过程来看，计算过程中涉 及到风洞试验和随机振动分析等复杂过程，不易为工程设计人员所掌握，因此迫切需要研究简便的 建筑结构抗风设计方法。等效静力风荷载理论就是在这一背景下提出的。其基本思想是将脉动风的 动力效应以其等效的静力形式表达出来，从而将复杂的动力分析问题转化为易于被设计人员所接受 的静力分析问题。等效静力风荷载是联系风工程研究和结构设计的纽带3，是结构抗风设计理论的 核心内容，近年来一直是结构风工程师研究的热点之一。等效静力风荷载的物理意义可以用单自由度体系的简谐振动来说明世108。图1.3气动力作用下的单自由度体系对如图1.3的单自由度体系，在气动力P(t)作用下的振动方程为: mx + cx + kx = P (t)考虑粘滞阻尼系统，则振动方程可简化为：(1.4.2)无+ 2&(2兀 f * +(2兀 f )2 x = P 00式中f =顽为该系统的自振频率，& =二为振动系统的临界阻尼比。0 2兀2展m假设气动力为频率为f的简谐荷载，即P(t)= Feg，那么其稳态响应为： 0/)F k(1.4.3)X U )=7r-07r ei2 冗 ft1 (f f + 八(2& f f )进一步化简有:J1-(f./f0)2 2 +(2& fx (t )= Aei(2n ft -v)2& f f，V = arctan一7p，A 为振幅，1 -(f fj(1.4.4)为气动力和位移响应之间的相位角。现在假设该系统在某静力F作用下产生幅值为A的静力响应，那么该静力应该为:F = kA =0(1.4.5)1 -( f/fo2 +(2& ff如果不考虑相位关系，静力F与简谐气动力P(t)将产生一致的幅值响应，则这两种荷载之间 存在一种“等效”的关系，那么F可以称为P(t)的“等效静力风荷载”。从上面这个简单的实例可以很清楚的体会到，所谓等效静力风荷载是指这样一种静力荷载，当 把它作用于结构上时，其在结构上产生的静力响应(不仅指代位移响应，也包括内力响应等)与外 加气动力荷载产生的动力响应最大幅值是完全相等的。本文中，将动力响应的最大幅值称为峰值响 应，或目标响应。等效静力风荷载理论的提出和发展等效静力风荷载(Equivalent static wind loading, ESWL)理论研究始于高层、高耸结构。1967 年，A.G. Davenport率先引入随机振动理论，建立了结构抖振响应分析的理论框架，并借助阵风荷 载因子(Gust Loading Factor, GLF)这一概念将复杂的动力分析问题转化为易于被设计者接受的静 力分析问题，从而开创了等效静风荷载理论研究的先河109。其后，先后有很多学者进行过等效静 力风荷载的探讨，并且提出了多种计算方法，但大多是针对高层结构而提出的一系列改进措施108, 110-119。上世纪九十年代，Kasperski(1992)在研究低矮房屋的等效静力风荷载时，重新审视了阵风荷 载因子法的不足，提出了适用于刚性屋面的荷载一响应相关(Load-Response-Correlation, LRC)法 120, 121 ，用于计算其背景等效静力风荷载。LRC法的提出和发展，使得等效静力风荷载的物理概念 更加清晰。随后，LRC法被广泛的应用于大跨度屋盖结构等效静力风荷载的计算71,96。LRC法的优 点是，它利用荷载和响应之间的相关系数来确定等效静力风荷载，这使得求得的等效静力风荷载是 实际可能发生的。在LRC法的基础上，Holmes等人(1996,1999)建议采用LRC法和等效风振惯性力相结合的办 法来表示等效静力风荷载，并且给出了平均风荷载、背景风荷载以及代表多阶共振分量的惯性风荷 载一起组合的等效静力风荷载形式122(或称为三分量组合形式)。之后，不断有学者对三分量法 提出改进和完善2 6, 7, 45, 97-105, 107。到目前为止，已经出现多种静力等效方法，下面详细介绍几种主要的方法。1.4.2.1阵风荷载因子(GLF )法Davenport ( 1967)引入“阵风荷载因子”(Gust Loading Factor, GLF)来考虑脉动风荷载对结 构响应的放大109，这种简单可行的方法得到发展并运用到实际工程中，成为制定高层建筑风荷载 规范的主要依据。阵风荷载因子法定义峰值响应与平均响应之比一一“阵风荷载因子” G来表征结构对脉动荷载 的放大作用。作用在结构上以某个响应等效的静力等效风荷载可用下式计算(1.4.6)p (z)= G(Z)p(Z)式中，p (z)为平均风荷载，阵风荷载因子G (z )由下式确定：c() r(z)(1.4.7)G (=汩其中r(z)表示峰值响应，r (z)为平均响应。r(z)可以表示为:(1.4.8)得到(1.4.9)r(z)=r(z)+ 舛(z)r其中幺为峰值因子，b r (z)为计算得到的某个响应的均方根值。将(1.4.8)代入(1.4.7)，b (z)G (z ) = 1 + g j()利用阵风荷载因子法来表示静力等效风荷载简单方便，因而在近年来的大跨度屋盖结构抗风研 究中应用也很广泛。目前对封闭平屋盖等效静力风荷载的研究一般都采用了阵风荷载因子法。例如 Marukawa(1993)针对来流紊流度、屋盖的几何特性和梁的结构特性为阵风荷载因子提供了经验公 式123。Ueda(1994)采用同步测压技术研究了梁柱框架结构平屋盖的风振响应124，特别研究了来 流紊流对风荷载的影响，提供了比文献123更详尽的阵风荷载因子表达形式。Uematsu根据封闭平坦 矩形屋盖的结构形式，把平坦矩形屋盖分为主次梁体系屋盖和空间整体体系屋盖两大类，前者由互 相平行的主梁作为承重结构，主梁之间通过次梁连接，结构振型为主梁在竖向的振动，第一阶振型 可以用一维的正弦曲线描述；而后者为空间网架，在风荷载作用下屋面发生类似弹性板的竖向振 动，振型可以用两个正弦曲线的乘积形式描述。Uematsu(1997)对不同跨高比的第一类平屋盖在 不同流场中进行了刚性模型试验125，用第一阶模态力计算了主梁的动力反应，发现靠近屋盖边缘 的主梁最大风振反应发生在风向垂直于梁轴线的情况；而位于屋盖中央的主梁其最大风振反应发生 在来流平行于梁轴线的情况。根据这个规律对第一阶模态力推导的梁阵风荷载因子公式进行了简 化，提出了适合工程运用的经验公式，其中考虑了紊流度、结构跨高比、主梁位置等因素。 Uematsu(1996，1997)还研究了第二类平坦矩形屋盖126,127，研究方法与第一类矩形平屋盖基本相 同。由于其振动形式与第一类矩形平屋盖不同，所以最不利的工况为来流垂直于屋盖边缘的情况。 对阵风荷载因子的研究表明，当折减频率比较小时，阵风响应因子受结构跨高比的影响较大，并且 此时的等效风荷载比按准定常方法得到的风荷载要大很多。Uematsu(1999)采用类似平坦矩形屋 盖的方法进一步研究了圆形平屋盖的风振响应83。文中用考虑第一阶模态的阵风荷载因子经验公式 (包含了高跨比及来流紊流的影响)计算了几个圆形平屋盖的位移及弯矩，发现计算结果与时程分 析结果吻合得很好。Uematsu的方法优点在于计算简便、快捷，但仅考虑了一阶模态的贡献，忽略 了高阶振型的影响。阵风荷载因子法同样被用于结构外形相对复杂的大跨度屋盖结构128。尽管阵风荷载因子法使用很简单，但有很大的局限性。从式(1.4.6)可知，该方法给出的静力 等效风荷载是与平均风荷载同分布的。由于大跨度屋盖结构各响应的阵风响应因子常常差别很大， 就可能导致某响应对应静力等效风荷载作用下的该响应大小，并不是所有静力等效风荷载作用下的 最大响应，这样易导致设计人员的误解。另外，如果结构的平均响应(荷载)为零时，GLF法给出 的阵风荷载因子将会出现无穷大(零)的情况。1.4.2.2惯性风荷载(IWL)法实际上，保证控制点响应等效的静风荷载分布形式存在无穷多个，Davenport提出的GLF法及其 改进方法都是假定等效静力风荷载的分布形式同平均风荷载，并没有体现响应出现极值时结构真实 的最不利荷载分布。惯性风荷载(IWL)法129-134从结构动力方程出发研究等效静力风荷载的分 布，认为脉动风对应的等效静力风荷载可以用结构的惯性力表示，其分布形式是真实的最不利荷载 分布。其主要思想是：如果结构第，阶振型七(z )在结构上的模态坐标标准差为b尸则相应于该 振型的惯性力为m (z)也2。出(z ) 135。下面证明在惯性力m(z)s 2。(z)作用下结构产生的响应为b 8 (z)。 j j j j j在此惯性力下的广义力为(因振型对质量的正交性，其它阶振型的广义力均为零)(1.4.10)j L 8. (z )m (z )s 2b 8 (z )dz而在此广义力作用下的广义模态坐标为，(1.4.11)由此可以证明惯性力m(z)sj2b j8j (z)作用下结构产生的相应为b惯性风荷载法实际上也是一种阵风荷载因子法，只不过其阵风荷载因子由惯性力来表示。由于 中国建筑结构荷载规范GBJ中采用此方法，因而惯性风荷载法习惯上也称为GBJ法。在中国建筑结 构荷载规范中，对于主要为第一阶振型起作用的结构(对于多阶模态作用的结构可用相同的方法计 算阵风荷载因子)，阵风荷载因子(中国规范GBJ9-87称风振系数)为：()m (z )gb8(z)(1.4.12)心=以搭(J其中s 1为第一阶自振圆频率。显然，GBJ法给出的阵风荷载因子与结构的质量分布和动力特性有关，其静力等效风荷载与平 均风荷载的分布是不同的，GBJ法赋予了静力等效风荷载明确的物理意义。但GBJ法也有不足，虽 然它给出的共振等效风荷载和响应与实际值是相同的，但背景等效风荷载和其它响应则与实际情况 不同，另外GBJ法无法处理多模态的耦合情况，因而不适用于大跨度屋盖结构。类似于GLF法，如 果结构的平均荷载为零时，GBJ法给出的风振系数也将会出现无穷大的情况。阵风荷载因子法和惯性风荷载法都用阵风荷载因子来反映总等效风荷载和平均风荷载之间的关 系；不同之处在于对阵风荷载因子的计算，前者认为阵风荷载因子等于动力响应与平均响应的比 值，而后者则将风振惯性力与平均风荷载的比值作为阵风因子来反映风荷载的脉动放大作用。以上 根据“阵风荷载因子”思想提出的静力等效风荷载方法写入了许多国家的高层建筑结构抗风规范。 使用阵风荷载因子法虽然简单方便，但直接把研究高层结构的方法搬到大跨度屋盖结构显然不合 适，因为大跨度屋盖结构相对高层结构而言，不论荷载还是响应特性都要复杂很多。1.4.2.3荷载响应相关(LRC)法Kasperski(1992)年提出的荷载-响应相关法，即LRC法120,121，是在研究低矮建筑风洞试验95, 136-138基础上发展起来的一种计算静力等效风荷载的方法。LRC法利用准静力的方法计算背景响应， 能够得到背景风荷载的等效分布形式，它的出现使等效背景风荷载的计算有了坚实的理论基础。 Kasperski(1992)认为，即使对于非高斯过程的荷载，LRC法仍然能给出具有很好近似程度的等效 荷载分布121。Holmes(1992)将LRC法与正交分解法结合在一起来表示等效的背景风荷载叩。从结构动力学可以知道，结构在低频部分的响应可以认为仅在弹性恢复力来抵抗外力的(俞载 道,1987)140。根据这一原理，LRC法考虑了结构上脉动风荷载之间的相关性，用准静力的方法 计算出结构表而等效背景风荷载。用准静态方法得到t时刻结构上某点的瞬态背景响应为(1.4.13)其中，p (z,t)为作用在结构上的脉动风荷载，I (z )为对应于响应r的影响线。相应的平均响应r(t)=j lp(z,t)I (zDz0rr = jlp (z)l (z)dz 0r 根据(1.4.13)，可得到响应的标准差b:r, Bb 2 =j l j lp (z , t) p (z , t )I (z ) I (z )dzdz r ,B0 012 r 1 r 21 2式中下标“B”表示背景响应。作为上式计算的一个中间步骤，响应和荷载的协方差可以表示为：b 2 =jl p (z , t) p (z , t)I (z )dz = p b b (z ) r, p012 r 11r, p r, B p 2(1.4.14)(1.4.15)(1.4.16)Prp定义为z0位置的背景响应与z2位置处荷载之间的相关系数，那么可以把(1.4.15)改写为:b 2 r ,B0 r把(1416)代入(1417)得：I (z 用孚(z,t)p(z ,t、)I (z 匾01dZ 2(1.4.17)=lI (z)Tpb (z2)dz2(1.4.18)P (z)= g p b (z)eBB r, p p(1.4.18)即为:b= jlI (z)P (z)dzB r, B0 reB(1.4.19)式中为为背景响应的峰值因子。那么(1.4.20)如果定义:可见，peB (z )为对应于背景响应值gb r B的等效静力风荷载。从上述的推导可知，LRC法通过荷载与响应之间的相关性分析，过滤掉了对所考察响应没有贡 献或贡献很小的脉动荷载，从而体现了对响应有效的脉动荷载分布。对于响应和荷载的均值为零的 情况，LRC法均可以给出解答。利用LRC法，某响应对应静力等效风荷载作用下的该响应的大小， 一定是所有静力等效风荷载作用下的最大响应。由于Prp小于1，所以LRC法定义的等效荷载分布 在峰值压力分布(p = p + gf )的包络线内，由此可见LRC法是一个易于理解的等效静荷载概 念。另外，LRC法采用了准静态方法，实际上已经包含了构成背景分量的所有模态的贡献。但由于 LRC法只能计算背景响应，因此将它应用于大跨度屋盖结构时，只适用于强刚性屋盖的等效荷载分 布的计算，而对于弱刚性或柔性屋盖则不适用，因为这时共振响应占了相当大的比重。1.4.2.4三分量法周印、顾明(1998)根据Davenport提出的把顺风向响应处理为平均、背景和共振分量的思想提出 了用这三个分量的组合来表示高层建筑静力等效风荷载108,113 的方法，其中背景等效风荷载用LRC 法表示；由于共振时结构等效风荷载等于惯性力，因此用代表共振分量的等效风振惯性力表示共振 等效风荷载。和以上方法类似，对于背景、共振响应均不能忽略的风致结构振动问题，Holmes (2002)采用了将LRC法和等效风振惯性力相结合的办法来表示静力等效风荷载，并给出了平均 风荷载、背景风荷载以及代表共振分量的惯性风荷载一起组合的等效风荷载形式122。下面简述其 计算的基本过程。对于背景响应等效静力风荷载，上小节已经给出了详细的推导；而对于共振响应等效静力风荷 载，与1.4.2.2小节的惯性力略有不同，区别在于这里仅考虑各阶模态的共振响应。对于某结构，其 第i阶模态的位移均方响应可以表示为45：b 2 =s (n)dn = !J H (n)2 S (n)dn(1.4.21)qi o qiK 2 0 iQi式中Sqj (n)为第i阶模态的位移响应谱，可以表示为：S . (n)= K H，(n)|2 SQ (n)(1.4.22)i其中k.为第i阶模态的广义刚度；H (n)为第i阶模态的频响函数，Ih. (n)|2则称为第i阶模态 的位移导纳；SQi(n)为第.阶模态力自谱。对于一般的荷载，其输入力谱的卓越频率均比结构的 自振频率小很多，因而对于位移导纳|H. (n)|2，其值在第.阶自振频率七附近有较大的峰值，而在其他频率上的值则很小，因此可以将位移导纳分为两个部分：在其自振频率ni附近，这部分认为 是由于结构的共振引起，称为共振区域；另一部分是除共振区域外的其他区域，该部分与结构的共 振无关，称为背景区域。通过推导45 可以得到：1./ 、 冗 n / 、b 2 三 j S (n)dn +S (n )(1.4.23)qi K 2 0 Qi4&Qi iii上式的前面部分为第i阶模态位移响应的背景分量，后面部分为第i阶模态响应的共振分量。如用 qi,R表示第，阶模态响应的共振分量，那么它可以表示为:(1.4.24)(1.4.25)b2 =-L-当S (n )qi,RK 2 4& Qi iii因此第i阶模态的共振响应等效静力风荷载可以表示为：P = gm (z ) 2巾(z )beR,iR,ii i qi, R上式考虑了共振响应峰值因子gR.。根据三分量法的思想，总的峰值响应可表示为r = r + ；g2b2 +Z g2 a2(1.4.26)T B r ,B.R,i qi ,R上式其实就是常用的SRSS法，它基于小阻尼及各模态之间的交叉项可忽略的假定。相应的静力等效风荷载(包括平均、背景和共振分量)为p (z )= p (z)+WPB (z )+Z WP (z)(1.4.27)i式中P (z)为LRC法得到的背景响应等效静力风荷载，见(1.4.19)式；P (z)为第i阶共振eBeR,i响应的等效静力风荷载；WB、WRi分别为PeB (z )和tRi (z )的权值系数，由下式确定122W = % r，B(1.4.28)B(g 2b 2 +乙 g 2 b 2 )1/2W = r R(1.4.29)R,i(g 2b 2 +乙 g 2 b 2 )1/2i上述方法虽然考虑了背景和多个模态共振响应，但必须假定参振模态之间能够很好的分离。而 在工程实践中，常见的大跨度屋盖结构不仅要包含多振型的贡献，而且应该考虑不同振型响应之间 的耦合影响。这使得Holmes提出的方法在处理大跨度屋盖结构的静力等效风荷载问题上遇到了障 碍。