风速时程数值模拟研究(重要).docVIP

文章编号:1008-3723(200601-0001-03 风速时程数值模拟研究 白泉,朱浮声,康玉梅 (东北大学资源与土木工程学院,辽宁沈阳110004 摘　要:随着结构的大型化、复杂化,对于风荷载更加敏感,风荷载常常起主要甚至决定作用。时域分析法是对结构进行抗风分析的较好方法,而时域分析首先要有时程样本,利用已知的频域信息通过计算机数值模拟重现时程样本是比较好的解决办法。本文利用谐波叠加法模拟了风速时程样本。通过算例表明,该方法模拟精度较高。另外对部分参数的合理选取进行了讨论,并给出了取值的建议。 关键词:风速时程;数值模拟;谐波叠加法 中图分类号:T K81　　　　文献标识码:A 0　引言 近年来,随着生产生活的需要,大型结构日益受到人们的青睐。大跨度桥梁跨度已从20世纪初的500米到世纪末已逼近2000米;大跨度空间结构跨度也已超过了200米;高层建筑方面,新近建成的台北101大楼达到了508米。而且随着新技术、新材料、新工艺、新型式、新设计方法的应用,工程结构将日趋大型化、复杂化,这也必然使结构更轻更柔,对于风荷载将更加敏感。对于这些大跨、高层、高耸结构,风荷载往往起主要甚至决定作用,因此对于这些结构在风荷载作用下的响应研究越来越受到科研工作者的重视[1][2][3][4]。 结构的抗风分析方法一般分为频域法和时域法[3][5][6][7]。频域法是建立在线性假设前提下的,因此只能对结构进行线性或线性化分析,尽管概念清晰、简便,在工程中得到广泛应用,但它不能给出结构响应的相关函数、瞬态反应,一般也不能分析强非线性结构。因此,要对结构进行较精确的非线性分析只能借助于时域法。进行时域分析必须要利用风速时程样本,而风速时程样本的获取只能通过实际强风记录、风洞实验或数值模拟来得到。由于实际记录到的强风作用过程应用于工程实际还不能普遍实现,因此风洞实验和风速时程数值模拟是解决问题的有效方法。在没有风洞实验数据的情况下,利用已知的频域信息通过计算机数值模拟重现时程样本是比较好的解决办法。时程样本模拟的好坏,即所模拟的风与自然风在统计特性、时间相关性和空间相关性等方面是否一致或接近,对于时程分析的结果具有较大的影响,因此模拟出尽可能接近和满足自然风特性的时程样本模拟具有重要意义。 1　风的基本特性 大量实测记录表明,顺风向的风速可看作为由两部分组成,即周期在10分钟以上的长周期部分和周期在几秒钟至几十秒以内的短周期部分。通常长周期部分远离一般结构物的自振周期,其作用属静力性质,而短周期部分则与结构 收稿日期:2005-06-21 作者简介:白泉(1973-,男,吉林梨树人,东北大学资源与土木程学院博士研究生.的自振周期较为接近,因而其作用具有动力性质,根据风荷载的这一特点,在工程结构应用中通常将作用在结构物上的风荷载视为平均风(静力风和脉动风两部分的共同作用[8]。 大气边界层内平均风速的变化可以通过平均风剖面来描述。土木工程领域风速廓线常用指数函数来表达。Dav2 enport根据多次观测资料整理出不同场地条件下的风剖面,提出平均风速沿高度变化的规律可用指数函数予以描述,即: v(z v b =(z z b α(1 式中,z b,v b—标准参考高度和标准参考高度处的平均风速; z,v(z—任一高度和任一高度处的平均风速;α—地面粗糙度指数。规范中,取z b=10m。 而脉动风具有随机性,随时间和空间随机地变化,在数学上属随机过程范畴,通常将其假定为具有零均值的平稳高斯随机过程,风速时程的模拟主要针对脉动风速而言。脉动风速的特性可用功率谱和相关函数予以描述。脉动风速的功率谱主要反映脉动风中各种频率成分对应的能量分布规律。相关函数则反映各点脉动风速之间在时间或空间的相互影响关系。功率谱与相关函数之间可通过维纳-辛钦公式进行变换。风速谱按方向可划分为:水平阵风功率谱、竖向阵风功率谱和横向阵风功率谱。目前应用考虑较多的是水平阵风功率谱。按是否考虑湍流积分尺度随高度的变化有两大类:一类是不考虑湍流积分尺度随高度的变化,如Davenport谱,Harris谱等,另一类是考虑随高度的变化,如Simiu谱、K aimal谱等。 我国规范采用的Davenport谱具有如下的形式: S v(n=4k v210 x2 n(1+x24/3 (2 式中,S v(n为脉动风速功率谱;n为脉动风频率(H z;x= 1200n V10 ;V10为10米高度处的平均风速(m/s;k-表面粗糙度系数。 2　风速时程模拟的算法 目前,对于风速时程的模拟方法主要有二大类即谐波迭加法和线性滤波法。其中谐波迭加法是基于三角级数求和 的频谱表示法,采用以离散谱逼近目标随机过程的随机模型的一种离散数值模拟方法,该算法简单直观,数学基础严密,适用于任意指定谱特征的平稳高斯