当前位置:文档之家› 脉动风场模拟技术的研究与进展

脉动风场模拟技术的研究与进展

第28卷第6期地震工程与工程振动V01.28No.62008年12月JOURNALOFEARTHQUAKEENGINEERINGANDENGINEERINGVIBRATIONDec.2008

文章编号:1000—1300(2008)06—0206一07

脉动风场模拟技术的研究与进展

张希黔1’2,葛勇1,严春风3,晏致涛1

(1.重庆大学土木工程学院,重庆400030;2.中国建筑第i工程局,湖北武汉430070;3.成都市规划管理局,I匹t)ll成都610031)

摘要:通过对脉动风场模拟技术近2()年来圈内外研究成果的总结,阐述了脉动风场模拟研究的发

展历程。酋‘先简要介绍了近地风场的特性和脉动风场模拟技术的研究现状,然后从脉动风场模拟技

术的研究方法出发,分别介绍了谐波合成法、线性滤波法和小波方法在风场模拟技术中的研究发展

现状和所面临的问题,最后对风场模拟技术的未来趋势进行r展单。

关键词:风场模拟;谐波合成法;线性滤波法;小波方法

中图分类号:279.77“文献标志码:A

Advancesinresearchofsimulationtechnologyoffluctuationwindloading

ZHANGXiqianl一,GEYon91,YANChunfen93,YANZhita01

(1.CollegeofCivilEngineering,ChongqingUniversi!Y,Chongqing400030,China;2.ChinaConstructionThirdEngineeringDivision.CORP.LTD,Wuhan43(X)70,China;3.ChengduCityPlanningConservancy,Chengdu610031,China)

Abstract:Accordingtothefluctuatingwindfieldsimulationtechnologyathomeandabroadinthepast20years,asummaryofresearchresultsexpatiatesthedevelopmentprocessofthefluctuationwindfiledsimulationtechnology.Firstofallwegiveabriefingonthecharacteristicsofboundarylayerwindandthestatusquoofthefluctuatingwindfieldsimulationtechnology.Andthen,basedontheresearchofthesimulationtechnologyoffluctuationwindfield,weintroducethecurrentsituationandproblemsofharmonicwavesuperimposingmethod,linearfilteringmethodandwaveletmethod.Finally,weprospectthefutureoftheresearch.

Keywords:windloadingsimulation;harmonicwavesuperimposingmethod;linearfilteringmethod;waveletmethod引言

对于高耸结构、高层建筑结构、大跨度空间结构、大跨桥梁结构及高压输电塔线体系等,风荷载是结构设计时所必需考虑的一类重要的随机动力荷载¨1。由于经典随机振动理论一般是在二阶统计意义七考察随机过程,因而导致在结构非线性风振响应分析与动力可靠度分析中的一系列困难。事实上,为了正确地计算结构的抗风可靠性,应该设法求解结构反应的概率密度分布及其随时间的变化过程。在实际的大气边界层湍流中,脉动风速不仅是时间的函数,而且也随空间位置(x,Y,z)而变化,是一个单变量四维(1V一4D)随机场口J。由于目前条件的限制,现场实测的脉动风速时程信号非常有限,并且项目建成前所取得的短期风速记录资料未能考虑风速的空间相关性。用实测的风速时程进行风振分析有一定的局限性,利用计算机进行人工风速时程的模拟能很好地解决这个问题。

收稿日期:2008—08—06;修订日期:2008—09一】9

基金项目:国家自然科学基金项目(50708118);重庆市自然科学基金项目(2007BB0133)

作者简介:张希黔(1939一),男,教授,主要从事结构振动、土木水利施工技术研究.E-mail:geyongl202@sina.com

地震工程与T程振动第28卷

以丁):型挚;盯。

p掣(r)=堡学

(6)式中,肛。啦,和盯,、盯,分别为随机变量茏、Y的平均值和均方根值。

在频域中脉动风的相关性一般用相干函数(coherencefunction)来表示。风洞试验和实测表明,相干函数是一条指数衰减曲线一]。文献[8]中,Davenport给出了纵向脉动风速的竖向和横向相关函数的表达式为r

∞【Ci(Yl—Y2)2+c:(彳1一z2)21Ⅳ2、c。h(Ⅷ)=expI2盯

丢[u(z。)+u(Z2)]‘7)^L一一I,一、,J\二

/式中,∞为圆频率;r空间两点的距离;y。、Y:、z。、彳:分别为空间两点的横向坐标和竖向坐标,两点的连线与平均风速的方向垂直;V(z。)、U(z:)分别为高度z。、z:处的平均风速;C,、C;分别为横向和竖向相关的指数衰减系数,其取值可以由试验确定,文献[8]根据风洞试验结果建议C,=16、C;=10。

Panofsky[91给出了纵向脉动风速的纵向柑干函数的表达式为Coh(删)=exp(一告掣1(8)、二"II-u\z,/

式中,C。为顺风向的指数衰减系数,在陆地I:取6.0,在海面上取3.0,茗。,z:为窄问两点的纵向^|^标。2脉动风场模拟技术的研究现状大气风脉动可近似认为是各态历经、零均值、平稳高斯过程。MonteCarlo方法是模拟随机过程的一种直接有效的方法,是包括一系列随机参数模拟以及用确定性方法估计其后的系统响应的一种统计样本试验【l0I。经过几十年各国学者的努力,目前无论是稳态或非稳态,均匀或非均匀,一维或多维,单变量或多变量,高斯或非高斯随机过程都提出了一系列的模拟方法,但归纳起来,主要有三类。一类是利用i角函数叠加的谐波合成法(WAWS);第二类是基于数字滤波技术的线性滤波法如自回归算法(AR)、移动平均算法(MA)及自同归移动平均算法(ARMA)等等;第三类是利用小波在时域和频域上同时具有良好局部化特性,采用离散小波逆变换重构风速时程。

谐波合成法足一种利用谱分解和三角级数叠加来模拟随机过程样本的传统方法。1954年,Rice。¨o首先提出了谐波合成法的基本思想,当时只能模拟一维单变量的平稳高斯随机过程。Borgman_l副和Shinozu—ka【【卜蝤1解决了多维、多变量甚至非平稳随机过程的模拟问题。在此之后,他们提出的力.法得到了广泛应用,并在应用中不断得到改进。Yang]16I在计算中综合快速傅里叶变换技术,大大提高了该类方法的计算效率。YamazakimI还提出了用于模拟非高斯随机过程的迭代谐波合成方法。Deotatis叫8。应川FFT技术结合双索频率的概念模拟了各态历经的多变餐平稳随机过程——沿不同高度分布的多点脉动风速,模拟了各态历经的多变量平稳高斯随机过程。目前,由于模拟时问较长,风场的模拟大多是简化成一个“一维m个变量”的高斯平稳随机过程.厂(t),将目标谱密度近似考虑为实数形式。屈应辉¨9i、李永乐∞0j、韩大建』21I、白海峰‘223等将三维脉动风场简化为多个一维脉动风速场。曹映泓旧引应用特殊情况下谱矩阵的显式分解,大大提高了该方法的计算效率,成功模拟了大跨度桥梁沿主梁的脉动风速。丁泉顺。24等通过对谱分解矩阵引入插值近似,减少谱分解的次数,从而提高该谐波合成法的计算效率。罗俊杰125。-26’采用三次均匀B样条插值方法来拟合分解谱密度函数曲线,然后用递归优化算法进行矩阵的Cholesky分解,加快了脉动风场合成的速度。李黎m,28I、李杰㈣]、李元齐[3叫等将特征正交分解的强形式,即基于功率谱矩阵的随机过程特征正交分解(SeT)引入到谱表示法中,代替原有的对功率谱矩阵的Cholesky分解,极大地提高了计算效率。

第二类方法是利用数字滤波技术的线性滤波法。Gerch“1首先将线性滤波技术用于生成随机时间序列等工程问题,该方法以其方便快捷开辟了随机模拟的新思路。此后许多学者对线性滤波方法进行了不断改进,发展到了可以用AR,MA和ARMA模型模拟多维、多变量的随机过程。在这方面,Spanos[32J、Sama,ras[331、Kareem[341等都做了重要的工作。Reed[35|、Iazzunni[=16]、Huang[”1及曾宪武‘川等都在脉动风场模拟中应用了这些方法。其中,多变量随机过程的自回归模型AR(P)在线性滤波法中运用最多。胡雪莲等∞引利

第6期张希黔等:脉动风场模拟技术的研究与进展

小波分析是一个范围可变的窗口方法,可以用长的时间间隔来获得更精确的低频信息,用短的时间间隔来获得高频信息,其主要优点是提供局部分析与细化的能力,可以揭示信号的高阶导数不连续性和自相似性。正是这种特性,使小波被广泛用于信号分析中。基于Meyer小波基,Yamada[3纠等提出了一种有效的数值方法,利用离散的小波变换分析实测的脉动风,发现自然脉动风具有一定的间歇性和局部相似性。Kitaga—wa[加1根据文[39]的研究,同样基于Meyer小波基,利用小波逆变换的方法模拟随机风场的脉动风,这样模拟出来的随机风场的脉动风具有像自然脉动风类似的间歇性和局部相似性。陈艾荣H【’训等根据小波的多分辨率分析性质,首先提出了利用小波重构算法,选用db3小波基,模拟随机风场的脉动风。韩艳【433基于Meyer小波基,根据文献[39]的数值方法,利用小波变换分析洞庭湖大桥处自然脉动风的特征,根据实际脉动风的特征,利用Log—Poisson湍流模型㈣。拍1,得到一系列不同尺度的小波系数,再利用离散的小波逆变换得到随机风场的脉动风。徐闻H7’4刮等采用ARMA技术产生每一尺度小波系数,然后依据风场的空间相关特性修正随机场样本,最终合成了具有指定相关特性的空间风场。马昌恒I49J推导了含有时间信息的自相关函数估计和功率谱密度函数估计的表达式,利用小波函数消失矩建立了脉动风速时程模拟的小波分析模型,并取得了很好的模拟效果。

谐波合成法和线性滤波器方法在这一领域应用最广。然而,这些方法都不适用于非均一网格的随机场模拟。而且,这些过程很难加强随机场样本的局部求解。而小波变换正是具有非均一网格并能检测时间历程局部相似性的一种恰当方法。

3脉动风场的模拟方法

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

一维m变量的平稳高斯随机过程{Z(t)}(歹=1,2,…,m),其互相关函数矩阵和互谱密度矩阵分别为足(7-)和S(r)。互相关函数矩阵的元素与互谱密度矩阵的元素存在维纳一辛钦关系式。根据Shinozuka【5叫的理论,随机过程{.f(t)}的样本可由下式来模拟:

Z(t)=2 ̄/△∞,∑,,乏IHjk(tO?)Icos[∞zt一0jk(∞z)+咖肼](9)其中,巩(∞,)为S(OJ;)的Cholesky分解矩阵H(O.J。)中的元素,即s(OJ。)=H(∞,)日。’(∞:);日7’(山,)为日(∞,)的共轭转置矩阵。Ⅳ为采样频率点数,理论上应有¨∞,实际上Ⅳ为一充分大的正整数即可保证模

拟的精度,为了能在计算中使用F盯技术,一般取N=2“,OL为正整数。△∞=∞。/N为频率增量,OJ。。为上界截止频率,即当∞>∞。。时S(∞。)=o;咖。。为均匀分布于[o,2盯)区间的随机相位角;tOf为频率

f∞z=(1一—砉-)△∞,z=1,…,Ⅳ

k¨乏,u∞式中,酗,为均匀分布于(一△∞7/2,△甜7/2)的随机频率,△∞’<<△∞。鼬。的引入是为了消除模拟随机过程的周期性,但它以牺牲效率为代价,因此最好不用随机频率措施。但是如果不考虑跏,,则上式模拟随机过程的周期为4盯/Ato,要获得较长的模拟周期只能减小△∞,但上界截止频率不可能太小,而过大的J7\r值将大大增加计算的耗费,为了获得较长的模拟周期,Shinozuka"¨和Deodatis”21建议∞舡采用双索引频率

(11)

∞村=(1-1)△to+旦m—△∞,凡=1,…,m

则式(9)可表示为

iN

Z(£)’=2/△∞主,,二1日砖(0J“)Icos[∞^ft一0j,(∞“)+9“](12)式中oi,(∞f)为巩(咄)的幅角。上式模拟的随机过程的周期为21Tm/Aw。

3.2基于线性滤波法的脉动风场模拟

m个点空间相关脉动风速时程列向量的AR模型可表示为:

v(t)=一.∑.蟹≮口(t—kat)+Ⅳ(£)(13)式中,P为AR模型阶数;At是模拟风速时程的时间步长;帆为AR模型自回归系数矩阵;Ⅳ(t)为独立随机

脉动风场模拟技术的研究与进展

作者:张希黔, 葛勇, 严春风, 晏致涛, ZHANG Xiqian, GE Yong, YAN Chunfeng, YAN Zhitao

作者单位:张希黔,ZHANG Xiqian(重庆大学,土木工程学院,重庆,400030;中国建筑第三工程局,湖北,武汉,430070), 葛勇,晏致涛,GE Yong,YAN Zhitao(重庆大学,土木工程学院,重庆,400030)

, 严春风,YAN Chunfeng(成都市规划管理局,四川,成都,610031)

刊名:

地震工程与工程振动

英文刊名:JOURNAL OF EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION

年,卷(期):2008,28(6)

被引用次数:5次

参考文献(51条)

1.She Z;Leveque E Universal scaling laws in fully developed turbulence 1994

2.Dubrulle B Intermittency in fully developed turbulence:log-poisson statistics and generalized scale covariance 1994

3.Anselment F;Gangne Y;Hopfinger E J High-order velocity structure functions in turbulence shearflow 1984

4.李杰;刘章军随机脉动风场的正交展开方法[期刊论文]-土木工程学报 2008(02)

5.Harris R I The nature of wind.In the modern design of wind sensitive structures 1971

6.Davenport A G The spectrum of horizontal gustiness near the ground in high winds 1961

7.Yang W W;Chang T Y P;Chang C C An efficient wind field simulation technique for bridges[外文期刊] 1997(0)

8.Shinozuka M;Deotatis G Simulation of stochastic processes by spectral representation 1991

9.Shinozuka M;Jan C M Digital Simulation of random processes and its applications 1972(01)

10.Panofsky H A;Thomson D W;Sullivan D A Two-Point velocity statistics over lake Ontario 1974

11.Davenport A G The dependence of wind load upon meteorological parameters 1968

12.陈贤川大跨度屋盖结构风致响应和等效风荷载的理论研究及应用[学位论文] 2005

13.贺德馨风工程与工业空气动力学 2006

14.Kaimal J C Spectral characteristics of surface-layer turbulence 1972

15.Deodatis G Simulation of ergodic multivariate stochastic processed 1996(08)

16.Shinozuka M;Deodatis G Simulation of stochastic processes by spectral representation 1991(04)

17.马昌恒基于小波分析的风场模拟及大型储罐风致屈曲初步研究 2007

18.徐闻;叶继红;单建小波分析在随机风场模拟中的应用 2007

19.徐闻;叶继红;单建小波分析在空间随机风场模拟中的应用[期刊论文]-振动与冲击 2008(02)

20.韩艳;陈政清利用小波逆变换模拟随机风场的脉动风[期刊论文]-振动工程学报 2007(01)

21.王毅;陈艾荣基于小波方法的随机脉动风模拟 2003

22.陈艾荣;王毅基于小波方法的随机脉动风模拟[期刊论文]-同济大学学报(自然科学版) 2005(04)

23.Kitagawa T;Nomura T A wavelet-based method to generate artificial wind fluctuation data 2003

24.Yamada M;Ohkitani K Orthonormal wavelet analysis of turbulence 1991

25.胡雪莲;李正良;晏致涛大跨度桥梁结构风荷载模拟研究[期刊论文]-重庆建筑大学学报 2005(03)

26.Huang Z;Chalabi Z S;Bedford U.K Use of time-series analysis to model and forecast wind speed[外文

期刊] 1995(2/3)

27.Iazzunni A;Spinell P Artificial wind generation and structural response 1987(12)

28.Reed A;Scanlan R H Time series analysis of cooling tower wind loading 1983(02)

29.Li Y;Kareem A ARMA representation of wind field 1995(01)

30.Samaras E;Shinozuka M;Tsurui A ARMA representation of random processes 1985(03)

31.Spanos P D;Zeldin B A Monte Carlo treatment of random fields:a broad perspective 1998(03)

32.Gerch W;Yonemoto J Synthesis of multi-variate random vibration systems:a two-stage least squares ARMA model approach 1977(04)

33.李元齐;沈祖炎本征正交分解法在曲面模型风场重构中的应用[期刊论文]-同济大学学报 2006(01)

34.李国豪桥梁结构稳定与振动 1997

35.胡亮;李黎;樊剑基于特征正交分解的桥梁风场模拟[期刊论文]-武汉理工大学学报(交通科学与工程版)

2008(01)

36.李黎;胡亮;樊剑具有桥塔风效应的桥梁风场数值模拟[期刊论文]-振动与冲击 2007(05)

37.罗俊杰;韩大建大跨度结构随机脉动风场的快速模拟方法[期刊论文]-工程力学 2008(03)

38.罗俊杰;韩大建谐波合成法模拟随机风场的优化算法[期刊论文]-华南理工大学学报 2007(07)

39.丁泉顺;陈艾荣;项海帆大跨度桥梁空间脉动风场的计算机模拟[期刊论文]-力学季刊 2006(02)

40.曹映泓;项海帆;周颖大跨度桥梁随机风场的模拟 1998(03)

41.白海峰;李宏男大跨越输电塔线体系随机脉动风场模拟研究[期刊论文]-工程力学 2007(07)

42.韩大建;邹小江;苏成大跨度桥梁考虑桥塔风效应的随机风场模拟[期刊论文]-工程力学 2003(06)

43.李永乐;周述华;强士中大跨度斜拉桥三维脉动风场模拟[期刊论文]-土木工程学报 2003(10)

44.屈应辉;徐旭大跨度悬索桥脉动风场的数值模拟[期刊论文]-工业建筑 2006(增刊)

45.Deodatis G Simulation of ergodic multivariate stochastic processes 1996(08)

46.Yamazaki F;Shinozuka M Digital generation of non-Gaussian stochastic fields[外文期刊] 1986(07)

47.Shinozuka M Simulation of multivariate and multidimensional random process 1971(01)

48.Borgman L E Ocean wave simulation for engineering design 1969(04)

49.Rice S O Mathematical analysis of random noise 1954

50.曾宪武;韩大建大跨度桥梁风场模拟方法对比研究[期刊论文]-地震工程与工程振动 2004(01)

51.李杰;刘章军随机脉动风场的正交展开方法[期刊论文]-土木工程学报 2008(02)

本文读者也读过(10条)

1.李永乐.廖海黎.李明水.卢伟风洞短试验段大气边界层被动模拟技术[会议论文]-

2.王志春.宋丽莉.刘荣.叶燕翔边界层风速随高度变化的曲线模型分析[会议论文]-2007

3.李正农.刘艳萍.王莺歌.苏万林近地边界层脉动风数值模拟[期刊论文]-湖南大学学报(自然科学版)

2010,37(4)

4.韩艳.陈政清利用小波逆变换模拟随机风场的脉动风[会议论文]-

5.张京穗.刘章军.ZHANG Jingsui.LIU Zhangjun谐波合成法在动力荷载模拟中的应用[期刊论文]-四川建筑科学研究2007,33(4)

6.赵超.黎景宇.ZHAO Chao.LI Jing-yu脉动风的数值模拟[期刊论文]-山西建筑2009,35(34)

8.杨鑫炎高耸结构抖振时域分析[学位论文]2005

9.甘凤林.杨振伟.代晓光.GAN Feng-lin.YANG Zhen-wei.DAI Xiao-guang基于改进谐波合成法模拟输电塔线风速[期刊论文]-电力建设2009,30(10)

10.王娟.苏英志AR法在空间结构风场模拟中的应用[会议论文]-2007

引证文献(5条)

1.贾玉琢.肖茂祥.王永杰500kV架空输电线路风偏数值模拟研究[期刊论文]-广东电力 2011(2)

2.陈俊儒.吕西林上海中心大厦脉动风荷载模拟研究[期刊论文]-力学季刊 2010(1)

3.曹树森.秦剑.柯坚.刘晓红接触网空间相关随机风场的模拟研究[期刊论文]-机车电传动 2010(2)

4.赵超.黎景宇脉动风的数值模拟[期刊论文]-山西建筑 2009(34)

5.汪之松.孙波.李正良.肖正直输电塔线体系风洞试验模型的变比例问题分析[期刊论文]-重庆大学学报(自然科学版) 2009(11)

本文链接:https://www.doczj.com/doc/1d16180321.html,/Periodical_dzgcygczd200806029.aspx

相关主题
文本预览
相关文档 最新文档