高层结构脉动风场空间相关性数值模拟研究.pdf

高层结构脉动风场的空间相关性数值模拟研究 汪大洋 （广州大学土木工程学院，广东广州，510006 ） E-mail ：wadaya2006@ 摘 要：本文针对高层结构风振时域分析的随机风场和沿高度变化的Kaimal 脉动风速谱的 特点，采用线性自回归模型模拟考虑空间相关特性的高层建筑顺风向脉动风速时程，基于离 散傅立叶变换（DFT ）和快速傅立叶变换（FFT ）技术改进了自相关Toeplitz 矩阵的计算方 法，通过两个高层结构算例验证了所提思路的正确性和有效性，并研究了脉动风速时程沿高 度方向的风振特性、变化规律及其空间相关性效应。 关键词：高层结构 脉动风场模拟 AR 模型 傅立叶变换技术 空间相关性 中图分类号：TU352.1 O328 1. 引言 高层建筑高宽比较大，阻尼小，水平方向的抗侧刚度也小，属于柔性结构体系。高层 结构承担的水平荷载除地震作用外，主要为风荷载，虽然一般来讲风力作用没有地震强烈， 但结构所承受的风载作用极其频繁。因此，对于高层结构而言，风荷载是一种重要的设计荷 [1,2] 载，甚至起决定性作用 。 目前结构风效应分析主要有频域分析法和时域分析法两种，而随着计算机技术的不断 提高，基于随机模拟的结构响应时域分析法越来越受到工程设计和研究人员的重视。结构响 应时域分析的首要问题在于确定荷载时程，而解决这一问题的途径是依据风荷载的统计特 性，将其时程变化模拟出来。随机模拟的方法主要有谐波叠加法和线性滤波器法。谐波合成 法[3-5]及改进的谐波合成法是用一系列具有随机频率的余弦函数序列来模拟随机过程。当所 需模拟的维数比较大时，随机频率的生成相当耗费机时，因为它要在每个频率上进行大量运 算。线性滤波器法[6-11]则是将均值为零的白噪声随机系列通过滤波器，使其输出为具有指定 谱特征的随机过程，占用内存少，计算快捷。 此外，由于风力在空间上的分布变化及时间上的强度变化，风对结构的作用显得尤其 复杂多样，影响结构风效应的因素很多，分布变化及时间上的强度变化，风对结构的作用显 得尤其复杂多样，影响结构风效应的因素很多，不同地区各类地面条件的差异，结构本身条 件的变化，同一建筑物不同区域内形状、表明粗糙度等的不同都会导致结构上风荷载的改变。 因此，本文考虑了高层建筑风速谱随高度变化的特点，采用线性滤波法中的自回归模型模拟 高层建筑顺风向脉动风响应时程，利用离散傅立叶变换和快速傅立叶变换改进了自相关 Toeplitz 矩阵的计算方法，基于 Matlab 语言编制分析程序验证所提思路的正确性和有效性， 研究脉动风速时程沿高度方向的风振特性、变化规律及其空间相关性。 2. 脉动风场的数值模拟技术 风速观测记录表明瞬时风速包含两种成分：周期在 10 分钟以上的平均风和周期在几秒 钟的脉动风[12] 。在满足工程计算精度要求的前提下，可对风速时程作以下假定：(1) 任意一 点处平均风速不随时间改变；(2) 脉动风速时程是零均值平稳随机过程；(3) 风速时程间具 有空间相关性。 本课题得到建设部科技项目(项目编号：06-K6-11)和广州市属高校科技计划项目(项目编号：2046) 的资助 - 1 - 风速时程本质上是随机时间系列，可用自回归模型（以下简称 AR 模型）模拟。M 个点 空间相关脉动风速时程V (X , Y, Z, t)列向量的 AR 模型可表示为[13] ： p V X, Y, Z, t =− Ψ V X, Y, Z, t −k∆t + N t （1） ( ) ∑ k (