基于插值技术谐波叠加法的风速场模拟.pdfVIP

第28卷第8期 振动与冲击 JOURNAI。OF VIBRATlON AND SHOCK 基于插值技术谐波叠加法的风速场模拟 李春祥，刘晨哲 (上海大学土木工程系，上海200072) 摘 要：风荷载的数值模拟在结构设计中起着非常莺要的作用。在土木工程脉动风速时程的各种数值模拟方法 中，谐波叠加法最为常用。而且，通过引入F丌和不同插值技术可以在不显著地影响模拟精度的情况下，大大地缩短模 拟计算所花费的时间。基于两次多项式、三次多项式、以及两种不同边界条件三次样条函数插值技术的谐波叠加法来模 拟一幢100 m高高层建筑上10个点的脉动风速时程，通过均方根误差(Root Mean Square Error，RMSE)和相对误差系数 (Error factor，E，)两项指标来评价不同插值技术与传统谐波叠加法(未引入插值技术)相比较的插值模拟计算精度，并且 记录各自所耗费的时间，综合比较模拟效率和精度以找出最为适合的一种插值技术。结果表明：基于第一种边界条件三 次样条函数插值技术的谐波叠加法为最佳的模拟方法。 关键词：脉动风速；数值模拟；谐波叠加法；土木工程；插值技术 中图分类号：TU973+．213 文献标识码：A 风荷载是各类建筑物、构筑物的重要设计荷载之 一。对大跨度桥梁、大跨空间结构、高层建筑、高耸塔 桅结构而言，风荷载常常起主要甚至决定性作用。为 了在风荷载作用下对这些结构进行时域分析，首先必 须获得一定时间段内的脉动风荷载时程或脉动风速时 程。确定风荷载的主要手段有风洞试验、现场实测、数 值模拟等方法。相比前两种方法，数值模拟方法更加 简便高效，可以作为前两种方法的有效补充和检验。 运用Monte—Carlo思想模拟多变量随机过程目前主要 有两类方法，即谐波叠加法与线性滤波器法。 谐波叠加法适用于任意指定谱特征的平稳高斯随 机过程，算法简单直观，数学基础严密，通过叠加一系 列具有随机相位角的三角函数来模拟多变量或多维随 机过程。Shinozuka?和Shinozuka和Jan旧1提出了使用 CAwS法和WAWS法来模拟平稳随机场的一般理论。 Li和Kareem¨1和DeodatisMl引入了F盯技术之后，大 大提高了谐波叠加法的效率，Ding等∞o将脉动风速互 谱密度矩阵简化为实对称矩阵，利用三次托格朗日多 项式插值方法减少了互谱密度矩阵Cholesky分解的次 数，进一步缩短了谐波叠加法模拟风速场所花费的时 间。当所需模拟风速场离散点的个数非常多时，运用 插值技术可以显著地提高模拟的效率，节省的时间是 惊人的。然而，插值技术众多，插值的精度和效率各 异。本文的主要目标：通过引入均方根误差(Root Mean Square Error，RMSE)和相对误差系数(Error fae． tor，耳)两项指标来评价四种插值技术的模拟精度和效 基金项目：国家自然科学基金资助项目(No 收稿日期：2008—08—07 第一作者李舂样男，博士，教授，博士生导师，1964年12月生 率，试图找出一种最适合谐波叠加法的插值技术。 1风速场基本特性 在三维风速场中，任一点风速可表示为： U=U(z)+“(髫，Y，z，t) (1a) 口=口(髫，Y，z，￡) (1b) 叫=彬(z，Y，z，￡) (1c) 式中：U(z)为顺风向平均风速；／／．(石，Y，彳，t)、移(戈，Y，z， ￡)、w(x，Y，z，t)分别为脉动风速在顺风向、横风向和垂 直方向上的投影，其大小随时问变化；￡表示时间；z，Y，石 分别为对应三个风向的坐标轴。 本文仅研究工程结构顺风向脉动风速的数值模 拟，即如下一维风速场的模拟。 U=U(z)+u(髫，Y，彳，￡) (2) 2脉动风速功率谱密度(PSD) 顺风向脉动风速功率谱密度(PSD)函数按照其研 究方法可分为两类。第一类风速谱是对强风观测记录 进行相关分析，获得相关曲线和相关函数，再通过傅里 叶变换求得风速谱的数学表达式。这类谱沿高度不 变。第二类风速谱将强风记录通过超低频滤波器，直 接测出风速的功率谱曲线，拟合出风速谱的数学表达 式。在风工程中，广泛使用的Davenport风速谱就属于 第一类谱： s(∞)=磊1×面熹等矛 (3) 式中：算：6—00—to；∞为脉动风的圆频率(rad／s)；石。o为 ’玎口10 10 nl处的平均风速(rn／s)；K为地面粗糙度系数。 本文在数值模拟中采用Davenport脉动风速功 万方数据 振动与冲击 2009年第28卷 率谱。 3数值模拟风速场的谐波叠加法 对于脉动风。Davenport提出相干函数： Co^雄(光，，Yj，zf；茹I，Y，z^；co)= expf一旦焦互互臣互2————2—1 ～2,rr v(z。) 。 (4) 式中：V(z。)=(V(彳j