基于谐波合成法的脉动风场模拟

基于谐波合成法的脉动风场模拟摘要在建筑结构的现代设计中，抗风设计在建筑安全设计逐渐占有的着越来越重要的地位。

但受困于现有实测数据的局限性，风模拟的重要性日益得到凸显。

本文首先简述了对风的认识、风荷载研究发展历程和风荷载对建筑的影响。

而后概述平均风、脉动风、旋涡等风场基本元素的特性，并且回顾近年来风模拟研究进展，介绍利用谐波合成法进行风的数值模拟的想法。

最后通过对模拟生成数据与现有实测数据的对比，验证了谐波合成法脉动风模拟的可行性。

关键词：风荷载；谐波合成法；脉动风1、绪论1.1、对风的基本认识空气的流动是风形成的本质原因，这是现象是由于太阳热辐射而产生的自然现象。

太阳的热辐射到达地表，使其温度上升高，地表附近的空气受热,因密度变小而上升。

向上的热空气因渐渐变冷密度增大而下落。

由于地表一直接受着太阳的照射，所以温度比高空空气高，它加热的空气又会继续上升，这种原因导致的空气循环就是风形成的本质。

虽然其流动的方向看似确定，但空气移动的速率和方向却是有着随时间、空间的变化而会产生随机性变化的特点。

风的这种特性是可以利用风速剖面、脉动风的湍流强度、湍流积分尺度、湍流功率谱密度函数等随机概率特性进行叙述解释的[1]。

1.2、脉动风模拟研究意义随着经济的飞速发展，生活水平不断提高，多种建筑形式也在不断进入我们的生活。

对于建筑结构提出了更高的要求。

从前的低矮建筑样式在防灾减灾方面主要考虑抗震性方面的因素。

而随着建筑结构向着更高、更大、更复杂的方面发展，对高耸结构、高层建筑结构、大跨度空间结构、高压输电塔线体系等结构的安全性研究就变得日益重要。

从前对于低矮建筑结构的抗灾要求已远远不能满足现今复杂建筑结构的安全需求。

其中风荷载是当今现代高层、复杂的建筑结构体系设计中必需要考虑到的一类非常重要的荷载。

事实上，我们为能计算出正确可靠的结构抗风性能，应该求出结构的反应概率密度分布关于时间变化的过程，并且在实际的大气边界层湍流中，脉动风速不仅随时间改变，而且还具有与空间位置的相关性。

但是现在的现实情况是没有足够的现场实测的风速数据，整合难度大并且参考价值比较的小,因为现在的实测数据大部分为点分布，没有区域性的资料，而且由于地理环境等因素的改变，数据也需要不断更新，当一个建筑要开工建设之前通常只有很短的时间留给风速资料的收集，而且这样的数据一般也不会考虑风速与空间的相关关系。

所以，由于实测资料的限制，使用实测风速资料的局限性很大，进而，利用计算机进行人工对风场的模拟是一个十分可行的方法。

2、风载荷2.1大气边界层简介需要首先说明的一点是，在实际工程设计中研究的风一般是在大气边界层之内，处于大气层的最底层，空气经过地面时，地面上各种障碍物，如植物、建筑、沟壑、山体等物体会使流动空气受到阻力而使其流动的速度减慢，因为其靠近大地表面受地面障碍物摩擦力影响最大，这种阻碍减速的作用借由大气中形成的层流而往上依次传播，并且伴随着高度的逐渐升高，它的作用会慢慢弱化，当超过某一高度后，这种作用的影响就基本不存在了，这个离地高度就被视作做大气边界层的厚度，它会随着天气状况、地理条件、地表租糙成度而变化，一般厚度处在三百米至一千米之间。

一般人类活动范围都在这个大气区间内，所以在工程建筑方面考虑风荷载的作用时，大气边界层的流体力学规律已能满足工程需要。

2.2风载荷的成分分类风载荷是由于空气流动而对建筑产生的压力，这种由空气流动产生的荷载也被称为风的动压力。

基本风压、地形地貌、大气的粗糙程度、距地面的距离，工程结构、建筑物外观设计等很多的原因都会影响到风荷载的强度。

风荷载由两种成分组成：平均风作用和脉动风作用，平均风是由周期在十分钟以上的长周期成分组成，其变化周期远远大于一般建筑物的固有振动周期，所以平均风的作用可近似使用静力方法进行计算。

但是脉动风就大大的不同了，它的变化十分的频繁，周期只有几秒到几十秒，这样问题就来了，脉动风变换的周期与建筑物的固有振动频率就比较的接近，这样的话脉动风就会引起结构的共振，因此在高度较高、结构较为复杂建筑的设计里，风荷载中脉动风的成分引起的振动是造成建筑结构破坏的最重要原因，所以在风荷载对建筑所造成的应力影响中，平均风所造成的静力作用和脉动风造成的动力作用，现在研究重点就在于对脉动风的动力作用的研究，从另一个角度讲，由于力学在静力方面的理论都已经很成熟，所以在现今抗风性能的研究主要对象也就是针对脉动风的研究。

2.3平均风的基本性质平均风的定义是在一定的时间段之内其所探测到的风向、风速进行平均取值，所以在求取平均风中可分为对平均风风向和平均风风速的求取，相对比较简单的计算方法是分别对所探测到的风向、风速求出它们的算术平均值，从严格意义上讲应是此时段内各时刻瞬时风的矢量合成后的平均量，以直角坐标中的三个风速分量的平均值或平均风矢量的方向和速度表示，平均风的风向和风速决定了大气传输的方向和速率。

在大气边界层厚度以下，平均风速与研究时使用的时距有关（我国普遍设定标准时距为十分钟），并且，它同样也会伴随高度的变化而发生变化，它的变化规律使用风剪切和风速廓线来进行描述，风速廓线可利用对数或指数律分布进行表达[2]。

2.3.1对数律从近些年的研究实验的成果可以看出，对数律分布律在离地一百米以下的大气底层时有着非常好的效果，大气学家认为在离地一百米以下的大气范围内的风速廓线描述中，可以近似忽略掉剪切应力的变化，此时使用普朗特对数律进行描述的效果很理想，所以在气象学遇到的实际问题中一般使用对数律进行描述[3]。

其数学表达式为：()⎪⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛=0*ln z z k V z V （2-1） 式中 ()z V ——距地面距离为z 处的平均水平风速；*V ——流动剪切速度或摩擦速度，它主要是用来量度气体流动时其内部的摩擦力，其可表示为()5.00*ρτ=V ，其中的0τ表示的是近地面空气的剪切应力，ρ表示空气密度；k ——卡门常数，普遍近似取为0.4；0z ——地面粗糙长度，不同地面情况的0z 可见表2—1；z ——样本距地高度；表 2—1 不同地表情况下的粗糙长度地形m z 0 沿海地区0.005—0.01 开阔场地0.03—0.10 建筑物较少的郊区0.20—0.40 建筑物较多的郊区0.80—1.20 大城市中心2.00—3.002.3.2指数律指数律分布可用来描述靠近地面处的风速随高度的变化规律，指数率分布描述方法的来源是由经验的总结形成。

因为使用指数律分布计算廓线相较使用对数律分布比较简单方便且两种方法的描述水平差别不是很大，所以在土木工程等工程领域的计算上一般使用的指数律分布来进行描述风速廓线[4]，描述风速廓线的指数律分布可表示成：()αss z z z V z V =)()( （2—2） 式中 ： V —— 平均风速；z —— 所要研究风速的高度；s z —— 标准参考高度；α —— 地面粗糙度指数 （详见表2—2）；上述式子中基于地面粗糙度指数α在一定高度（设为c z ）一下是常数，并且高度c z （详见表2—2）仅仅是式中指数α的函数的假定。

表2—2 Davenport 建议[3]的α值和c z 值地形类别α值 ()m z c 值 开阔地形0.16 275 郊区地形0.28 400 大城市中心 0.40 520在我国的工程规范中使用的指数风剖面中，标准参考高度一般取为m z c 10=。

2.4脉动风的基本性质生活中，人们虽然能够感觉出风大概的风向，使用简单的仪器也可以测量出风的速度。

但真实的风并不像我们感觉的那么稳定，而正是由于它的这种不稳定性才引起脉动风的产生，而脉动风的强度是与时间相关的随机变化的，并且地面的粗糙程度也是影响脉动风的原因之一，脉动风与地面粗糙度的关系从历史资料可以知道。

在地面粗糙成度较大的地区，脉动风的幅值比较大且频率也比较高。

由于脉动风极其不稳定的特性，对脉动风大多采用数理统计的方法来进行描述，在描述的脉动风时最常使用的统计性质有：湍流强度、湍流尺度、湍流功率谱密度，还有相干函数等等，这些统计特性可被用来描述大气运动中的湍流结构。

2.4.1湍流强度湍流强度是用来形容风速脉动的相对强度的量，它是由脉动风速和平均风速来定义的，其中前者的均方根值与后者之比就被定义为湍流强度。

脉动风可分为三个方向的分量，而在各个风量上都有自己的湍流强度，一般来讲顺风向的湍流强度u I 要大于横向湍流强度v I ，而最小的是竖风向的ωI 。

对于强度最大的顺风向的湍流强度，每个国家都有自己的经验公式，而我国的经验公式如下：()α-⎪⎭⎫ ⎝⎛=1010z I z I u （2—3） 式中： α——地面粗糙度指数（详见表2—2）；10I ——十米处的名义湍流水平，国家对四种地貌进行了规范取值，分别为0.12、0.14、0.23、0.39；此外，我国还规定在没有现场实测数据的时候可取：u v I I 88.0= u I I 55.0=ω2.4.2湍流积分尺度大气边界层中的湍流运动可以将其看成是由大小不同的旋涡组成的，这些尺度不同的旋涡有着不同的特性，它们在大气湍流的运动中也发挥着不同的作用，这些旋涡的平均尺度就被称为湍流的积分尺度，在一般的实际生活和工程运用中，通常情况只需要对顺风向的平均湍流尺度x u L 进行了解，并且由于湍流的冻结假设，即湍流在以u 的平均风速、在顺风向移动时，其内部湍流形态在运动过程中不会发生显著变化，即可列出x u L 的表达式：()⎰∞=02ττσd R UL u u x u （2—4）式中 ()τu R ——脉动风速1u 的自相关函数x ——顺风向脉动速度的方向2u σ——顺风向脉动风速分量的均方值当空间中的两点之间距离小于平均尺度，就可认为这两点经常同时处在一个旋涡之中，这样的两个点，他们的脉动速度就是相关的，涡旋作用增强；若这两个点间的距离大于平均尺度，则他们经常不在同一个涡旋之中，这样的话，它们的速度不相关，且涡旋作用也会相应的消减。

根据长期大数量的观测实验结果显示出，大气边界层中的湍流积分尺度的增减与地面粗糙度、离地高度有关系，地面粗糙度越大湍流积分尺度则越小；离地高度越高湍流积分尺度越大。

2.4.3脉动风速功率谱大气湍流运动是由前文所提到的很多大小不同旋涡组成的，其运动能量也是有这些大大小小的旋涡提供，而脉动风速的功率谱就是用来描述脉动风能量的分布情况，它可以体现出脉动风中不同频率的部分对总能量的贡献，旋涡的尺度大小与它的出现频率是呈现反比关系，旋涡的尺度大，脉动频率就会较低；反之旋涡尺度较小，脉动频率就会较高。

脉动风速谱可以按照不同方向分为顺风向和竖风向的脉动风速功率谱，本文仅考虑顺风向情况，经过多年的研究实验，学者先后提出了很多功率谱，如达文波特（Davenport ）谱、哈里斯（Harris ）谱、卡曼(Kaimal)谱、西蒙(Simiu)谱等等。