matlab谐波合成法模拟风速时程,基于谐波合成法的输电塔线体系风致响应分析

0引言塔线的耦合作用使动力风荷载作用下大跨、高柔的输电塔线体系表现出很强的非线性特征,导线与输电塔之间的风动力耦联作用不容忽略[1-2]。动力风致响应分析的正确性和精确度将直接影响到输电塔线结构的安全性[3],由于时域法可以考虑脉动风的时间相关性、空间相关性、结构非线性的影响,更精确地反映结构的耦合风振[4],所以输电塔线体系风致响应分析通常采用时域法。时域内进行风振响应分析须以脉动风速时程为基础[5],而我国风场的时程记录资料较少[6],只有采用数值模拟风速时程才能满足场地、风速谱特性、结构特点等实际工程的要求,因此采用数值模拟风速时程的方法得到广泛应用。目前国内外对风速时程进行数值模拟的主要方法有谐波合成(weightedamplitudewavesuperposition,WAWS)法和线性滤波法[7]。谐波合成法是以离散谱逼近目标随机过程的一种离散化数值模拟方法,适合于任意指定谱特征的平稳高斯随机过程,其特点是精度比较高、算法简单直观、数学基础严密,但计算工作量十分巨大。为此,学者们进行了相关研究,文献[8-9]提出了使用恒波叠加(constantamplitudewavesuperposition,CAWS)法和WAWS法模拟平稳随机场的理论,模拟多维或多变量均匀高斯随机过程;文献[10-11]成功地引入傅里叶变换技术以提高模拟计算效率。尽管如此,对于大跨越输电塔线体系,需要模拟的点数比较多,计算工作量仍然相当大。针对上述问题,本文以在建的1000kV汉江大跨越输电塔线体系为例,基于按时间抽取的快速傅里叶变换技术(decimation-in-timefastFouriertransform,DIT-FFT)改进的谐波合成法,利用MATLAB软件,对输电塔线体系进行了风致响应分析,得到了输电塔线体系的脉动风速时程。通过对模拟风速谱和目标风速谱的吻合程度、输电塔线体系不同部位的脉动风标准化互相关函数以及改进前后计算效率的检验,验证应用按时间抽取的快速傅里叶变换技术改进的谐波合成法仿真输电塔线体系风速时程的可行性。1脉动风1.1概念脉动风是由风的不规则性引起的,它的强度是随时间按随机规律变化的。它周期较短,作用性质完全是动力的[12-13],脉动风的基本特性主要包括功率谱密度和空间相关性。1.2脉动风功率谱密度在频域内功率谱密度函数的研究中,风工程界采用的是Davenport提出的沿高度不变的风速谱[14]221024310()4(1)1200vSnvkxfxxfv????=+???=(1)式中:Sv(n)为脉动风速功率谱;f为脉动风频率;k为地面粗糙度系数;v10为标准高度10m处的风速。1.3脉动风的空间相关性在自然风场中,由于湍流的不均匀性,不同高度计算点的风速和相位均不相同,会出现阵风首先作用于结构的较高点,经过一定时差后,作用于较低点的现象。即使在同一水平高度上的2点,也会因相位差的不同,使脉动风速时程存在相关性。顺风向风场脉动风速v为0均值随机过程,其统计特性可用互谱密度函数来描述。函数(为复数形式)表达式(对空间任意2点)为[15]sij(f)=sii(f)sjj(f)ch(f)exp(i?(f))(i,j=1,2,n)(2)式中:?(f)为相位角,其表达式为*****/8,0.2()51.25,0.20.250.25f?????+=??随机数,且*=?z/v;=2f;?z为高度变化量;v=(vi+vj)/2;随机数取[?,]上均匀分布的随机数;对于3维问题,相干函数ch(f)一般有如下形式[16]2ch(,,)2()()()e