基于非平稳风速的风场湍流特性研究

基于非平稳风速的风场湍流特性研究

风负荷是结构中应考虑的负荷，其动态特征主要由风场的水流特征决定。风场的湍流特性主要包括湍流强度、阵风因子、湍流积分尺度和脉动风速谱等。基于现场实测技术得到风速记录并加以分析,是目前最直接最有效的研究风场湍流特性的手段。近年来,国内已有不少文献[1-5]对建筑结构一定高度处(一般在顶部)风速的湍流特性进行了实测研究,为认识脉动风速湍流特性提供了宝贵资料;然而,我国规范以空旷平坦地貌下10m高度处的风场为参考标准,而已有的文献对该地貌同高度下的风场湍流特性的研究鲜有涉及。因此,开展10m高度空旷平坦地貌(标准地貌)的风场湍流特性研究具有重要的现实意义。在以往的风特性研究中,通常假定风荷载为平稳高斯随机过程;而大量现场实测风速记录本文根据某空旷平坦地貌下10m高度实测得到的1128h强风风速样本,采用EMD法建立非平稳风速模型,基于该模型计算各风速下湍流特征参数,包括湍流强度、阵风因子、湍流积分尺度和风速谱等,并与平稳风速模型的对应值进行对比,以研究风场湍流特性的合理表达。1非稳定型风速模型1.1平稳风速模型通常采用矢量分解法对实测风速样本进行分解,得到沿风速主向的纵向风速时程和垂直风速主向的横向风速时程,具体求解过程参见文献[10]。现有的风场特性分析假定实测风速记录为平稳过程,即基本时距内平均风速为恒定不变值,故实测风速可表示为:上述过程得到的风速模型称为平稳风速模型。计算该模型的脉动风速自相关函数,并通过快速傅里叶变换(FFT),可以估计风场的功率谱密度。然而,基于平稳过程假定导致平稳风速模型的脉动风速时程具有较强的非平稳性,且矢量分解法将偏离主方向的风速强行分解到该风向的垂向和平行向,加剧了该模型的脉动风速非平稳性。由于FFT只对平稳信号具有较好的处理能力,直接对该模型的脉动风速进行计算,将产生较大的误差。1.2平均风分量的确定大量强风实测记录表明,复杂地形等高湍流风场的风速平均值并非恒定的,而是随时间在不断变化,通过EMD法可以提取该平均风速变化历程。EMD法是由美国NASA的黄锷一般认为实测风速记录中风速周期远大于结构固有周期的成分可视为平均风分量,EMD法则将单调的趋势项作为平均风成分处理,采用这种筛选标准需要循环多次,比较耗机时。事实上,对各阶IMF的功率谱分析表明,低频段的IMF频率已经足够低,可将这部分IMF视为平均风成分。因此,一般将原信号筛选几次后,即可获得风速的平均风成分由EMD法筛选得到的趋势项是随时间不断变化的,即平均风分量是时间的函数。因此,实测风速可表示为:上述风速模型即为非平稳风速模型。由于平均风速随时间变化,故称之为时变平均风速;脉动风速可认为是统计参量随时间不变的零均值平稳过程。2风速时程的平稳性假设为研究我国规范规定的标准地貌(空旷平坦地貌)在标准参考高度(离地面10m)处的风场特性,设立了测风塔并采用三维风速仪进行连续观测。实测风场位于河南省开封市,场地周边平坦,四周为农田,按照《建筑结构荷载规范》有关地表粗糙度的分类标准,该场地为典型的B类风场。对2013年2月1日~3月19日期间,总计观测时间为1128h的风速进行了连续记录。将实测风速样本分割为1128个子样本,每个子样本的记录时长约为1h;对记录过程中数据缺失时段不予插值处理。图1(a)为某实测子样本的风速时程。采用矢量分解法计算每个子样本的纵、横风向风速时程,并应用轮次检验法对各风速时程进行处理,以检验信号的平稳性。处理前对数据缺失较多或风速失真严重的子样本予以剔除,得到970个有效序列。平稳性检验结果显示,当显著性水平为0.05时,仅436个子样本的顺风向风速时程具有平稳性,占总数的45%;而横风向风速时程样本中满足平稳性假定的子样本数为924个,占总数的95%左右。这说明,与台风等强风气候类似,本文研究的大部分实测风速子样本的顺风向时程具有非平稳性,风速平稳性假设与实际情形不符;而横风向风速时程则较好地符合了平稳性假定,可视其为一平稳的随机过程。下文主要研究实测风速的顺风向湍流特征参数。从实测风速记录中选出平均风速较大的52组风速子样本,选取原则为:1h风速记录较完整且平均风速大于6m/s。建立每个子样本的平稳风速模型和非平稳风速模型,分别得到两类模型的平均风速和脉动风速,如图1(b)~(d)所示。从图中可以看出,相比于平稳风速模型,非平稳风速模型的时变平均风速更好地反映了风速的变化趋势;该模型的脉动风速明显更具有平稳性,基本为零均值平稳过程。3.2非平稳风速模型湍流强度用于描述风场的脉动强度,是计算结构脉动风荷载不可缺少的参数,通常将其定义为基本时距内的脉动风速标准差与平均风速的比值式中,I(z)为高度z处的湍流强度;σ(z)、U(z)分别为同一高度处的风速均方根和平均值;基本时距取为T=10min。对前述52组风速子样本,计算两类模型下的湍流强度,计算结果如图3所示。从图中可以看出,非平稳风速模型的湍流强度比平稳风速模型的计算值略小,前者主要分布在0.16~0.20之间,后者大部分在0.18~0.22间;湍流强度在各平均风速段的变化不明显。3风场流特征分析3.1模型概率密度分布研究风速的概率密度分布可以检验信号是否具有高斯性,并为设计风荷载的确定提供依据。通常认为,脉动风的概率分布呈0均值的标准正态分布,其概率密度函数表达式为:式中,σ针对52组有效样本的脉动风速时程,计算两类风速模型的概率密度分布,并与标准正态函数进行对比。某样本的计算结果如图2所示。从图中可以看出,平稳风速模型的概率密度曲线宽而缓和,与实测风速分布偏离较大;非平稳风速模型的概率密度曲线窄而陡峭,且与实测结果吻合得很好。平稳风速模型的52组样本的平均偏度系数和平均峰度系数分别为0.28和6.23;非平稳风速模型的两个系数对应为0.21和2.96。这说明,采用非平稳风速模型描述的风场风速更符合高斯假定。3.3非平稳风速模型计算结果对比阵风因子也可用来表征风场的脉动强度。结构风工程中定义阵风因子为阵风持续期t从图中可以看出,非平稳风速模型计算得到的阵风因子明显比平稳风速模型的计算结果小;由非平稳风速模型计算得到的阵风因子为1.17,而由平稳风速模型计算得到的结果为1.54,两者差异显著。对图4各平均风速段的平均阵风因子计算表明,阵风因子随平均风速增大而略有减小。3.4纵向湍流积分尺度湍流积分尺度是气流中湍流涡旋平均尺度的量度式中,L符合Karman谱的风场其湍流积分尺度可根据自功率谱密度函数直接计算得到式中,S针对前述52组风速子样本,分别采用自相关函数直接积分法和Karman谱反推法计算两类模型下的湍流积分尺度,计算结果如图5所示。从图5可以看出,由自相关函数直接积分法计算得到的湍流积分尺度大部分集中在50m~200m,而由Karman谱反推法计算得到的湍流积分尺度较小,其结果大部分集中在0~100m;这可能是因为自相关函数直接积分法在计算中引入了Taylor假设,使两者结果有差异。采用自相关函数直接积分法计算时,平稳风速模型和非平稳风速模型计算得到的湍流积分尺度统计值分别为172m、139m,采用Karman反推法的计算值分别为117m、90m。两种方法计算的湍流积分尺度值均表明,非平稳风速模型的结果比平稳风速模型小,且分布更加集中。图5还表明,风场的纵向湍流积分尺度随平均风速而增大,由自相关函数直接积分法得到的湍流积分尺度随平均风速变化较缓慢;而由Karman谱反推法计算的结果则随平均风速增大而显著变大。说明当风场移动速度增大时,风场漩涡对行径范围的影响将扩大。为考察纵向湍流积分尺度随高度的变化关系,对国内学者的研究成果进行归纳,结果列于表1。从表1可以看出,随着竖向高度的增加,纵向湍流尺度逐渐增大;当竖向高度大于梯度风高度后,该值减小。3.5风速模型密度脉动风速功率谱反映风场能量在频域的分布情况,是描述风场脉动特性的重要指标。按照风速风向关系,有纵向、横向和竖向风速功率谱密度。常用的纵向风速功率谱有Davenport谱、VonKarman谱和Kaimal谱等,形式如式(8)。Davenport谱:VonKarman谱:Kaimal谱:式中,S研究表明,结构主要对处于低频段的风荷载作用较敏感,一般认为这一范围为0.001~0.4Hz;本次风场实测中,风速的采集频率为1Hz,根据采样定理可知,实测数据能够较精确地反映0.5Hz以下的实际风场能量。另外,由于风暴持续时间一般为1h左右,因此本文在计算风速谱时,时距选为T=1h。对实测的970个有效风速子样本,分别计算两类模型的风速自功率谱密度,并与Davenport谱、Kaimal谱和VonKarman谱三条经验谱曲线进行对比。以莫宁坐标(fz/u)为横坐标,归一化自谱(fS从图6可以看出,对于两类风速模型,在峰值频率右端,Davenport谱和VonKarman谱均能较好地描述风场能量随频率的分布关系,Kaimal谱则有所高估;而在峰值频率左端,VonKarman谱和Kaimal谱的功率谱估计均偏大,Davenport谱与实测曲线最接近。因此,Davenport谱能最准确地描述风场的能量分布。为定量评估两类模型下,各平均风速段的自谱与经验谱的偏差,采用式(9)的计算指标进行衡量式中,由式(9)计算得到的两类模型下各风速段的实测纵向风速功率谱与Davenport经验谱的偏差列于表2。从表中可以看出,高风速下Davenport谱与实测谱偏差较小,Davenport谱能够较精确地描述强风的风场能量分布。4平均风速、风速本文根据某空旷平坦地貌下10m高度处实测得到的1128h非平稳强风风速样本,建立平稳风速模型和非平稳风速模型,并计算两类模型下的各风场湍流特征参数,得到以下结论:(1)实测风速样本具有较强的非平稳性,非平稳风速模型考虑了平均风速的时变性,能更好地反映实际风场风速;(2)通过计算各湍流特征参数随平均风速的变化规律,表明,湍流强度随风速变化不明显,主要集中在0.16~0.20之间,规范规定同类场地下10m高度处的湍流强度为0.14,比实测湍流强度略小;阵风因子随平均风速减小,湍流积分尺度则随之增大;