一种基于POD插值的非平稳风场高效模拟方法

摘要

本发明公开了一种基于POD插值的非平稳风场高效模拟方法，涉及脉动风速时程模拟技术领域，其技术方案要点：S1、基于给定的风速演化功率谱密度函数S(ω,t)，定义时间‑频率插值点S2、基于给定的自演化功率谱密度函数和相干函数，计算时频插值点处的演化功率谱密度矩阵S3、分解时频插值点处的演化功率谱密度矩阵获得相应的分解后谱矩阵S4、通过本征正交分解(POD)解耦S5、通过插值主坐标获得通过插值特征向量获得S6、基于谱表示法基本模拟公式，采用快速傅里叶变换(FFT)生成模拟样本xj(t)。本发明能够解决非平稳风场模拟中时间和频率点数多导致需要进行Cholesky分解的次数多，计算量大，需要的电脑储存空间大的问题，可满足模拟的精度且能提高模拟效率。

说明书

技术领域

本发明涉及脉动风速时程模拟技术领域，更具体地说，它涉及一种基于POD插值的非平稳风场高效模拟方法。

背景技术

在对大跨度桥梁、高层建筑等结构的风振响应分析中，由于需要考虑结构非线性和气动非线性，传统的频域分析方法已不能满足要求，需要采用非线性响应时程分析法。在非线性时程分析中，模拟结构所受的脉动风场是一个重要的环节。

在众多模拟方法中，谱表示法的形式简单且具有严密的理论推导，是风场模拟最常用的方法之一。为了提高多点风场的模拟效率，该方法常与快速傅里叶变换(FFT)技术相结合。然而对于非平稳脉动风场的模拟，为了引入FFT，需要先对时频耦合的函数进行解耦，目前已有包括本征正交分解(POD)、小波分解、非负矩阵分解等处理手段，但是解耦计算量会随着模拟点数的增多而急剧增大。此外，若所模拟风场的相干函数具有时变特性，那么还需要大量的谱矩阵分解计算，这使得非平稳风场的模拟效率进一步大打折扣。

因此需要提出一种更为高效的非平稳风场谱表示模拟方法，以解决上述问题。

发明内容

本发明的目的是为了解决上述技术问题，提供一种基于POD插值的非平稳风场高效模拟方法，本发明的方法能够解决非平稳风场模拟中时间和频率点数多导致需要进行Cholesky分解的次数多，计算量大，需要的电脑储存空间大的问题，本发明的方法用极少量的时间-频率插值点代替原始时间-频率点，不仅能够满足模拟的精度也能提高模拟效率。

本发明的上述技术目的是通过以下技术方案得以实现的：一种基于POD插值的非平稳风场高效模拟方法，具体包括以下步骤：

S1、基于给定的风速演化功率谱密度函数S(ω,t)，定义时间-频率插值点

S2、基于给定的自演化功率谱密度函数和相干函数，计算时频插值点处的演化功率谱密度矩阵

S3、分解时频插值点处的演化功率谱密度矩阵

S4、通过本征正交分解(POD)解耦

S5、通过插值主坐标

S6、基于谱表示法基本模拟公式，采用快速傅里叶变换(FFT)生成模拟样本x

进一步地，步骤S1的具体方法为：

(1)沿时间方向的插值点采用均匀分布，将区间[0,T

其中，T

(2)采用一种可调四次多项式来描述非均匀分布的频率插值点，即非均匀频率插值点坐标为：

其中，ω

进一步地，步骤S2的具体方法为：时频插值点处的演化功率谱密度矩阵

其中，S

进一步地，步骤S3的具体方法为：

在每一个时频插值点处对演化功率谱密度矩阵进行Cholesky分解可得，分解公式如下：

此过程总共需要

进一步地，步骤S4的具体方法为：对于

1)将所有时频插值点处的

并将该常数矩阵的每一列看作一个列向量，使用POD找到一组最优标准正交基

其中，

2)然后，列向量在特征向量上的投影，即主坐标计算为：

其中，a

其中，N

3)

其中，

进一步地，步骤S5的具体方法为：

1)基于插值点处的

2)然后可近似重构出在原始时间和频率坐标处的H

进一步地，步骤S6的具体方法为：

1)基于谱表示方法，任意随机过程的模拟公式为：

其中，t是时间坐标；ω

2)将解耦表达式代入上式(12)中并进行数学等价变化得：

3)然后，采用FFT技术对模拟公式(12)中的三角函数求和项进行计算，快速地生成多点风场样本。

综上所述，本发明具有以下有益效果：

1、本发明的方法能够解决非平稳风场模拟中时间和频率点数多导致需要进行Cholesky分解的次数多，计算量大，需要的电脑储存空间大的问题；

2、采用本发明的方法模拟非平稳随机风场，通过少量的插值点来代替原始时间-频率点，不需要传统方法中的在每个时间-频率点处进行Cholesky分解，而且仅对时频插值点处的函数进行解耦，从而能够大大减少了计算量，提高了模拟效率；

3、本发明的方法可为大跨度结构，如桥梁、屋盖等抗风设计的时程分析提供输入的风荷载时程。

附图说明

图1是本发明实施例中的流程图；

图2是本发明实施例中模拟风速的相关函数验证曲线图；

图3是本发明实施例中模拟风速的演化功率谱密度函数验证曲线图；

图4是本发明实施例中不同模拟点数对应的耗时曲线图。

具体实施方式

以下结合附图1-4对本发明作进一步详细说明。

实施例：本发明的技术方案提供一种基于POD插值的非平稳风场高效模拟方法，如图1所示，具体包括以下步骤：

S1、基于给定的风速演化功率谱密度函数S(ω,t)，定义时间-频率插值点

(1)由于目标演化功率谱随时间变化缓慢，则沿时间方向的插值点采用均匀分布，将区间[0,T

其中，T

(2)考虑到目标演化功率谱在低频处随频率变化剧烈，而在高频处变化缓慢，则采用一种可调四次多项式来描述非均匀分布的频率插值点，即非均匀频率插值点坐标为：

其中，ω

S2、基于给定的自演化功率谱密度函数和相干函数，计算时频插值点处的演化功率谱密度矩阵

其中，S

S3、分解时频插值点处的演化功率谱密度矩阵

在每一个时频插值点处对演化功率谱密度矩阵进行Cholesky分解可得，分解公式如下：

此过程总共需要

S4、通过本征正交分解(POD)解耦

1)将所有时频插值点处的

并将该常数矩阵的每一列看作一个列向量，使用POD找到一组最优标准正交基

其中，

2)然后，列向量在特征向量上的投影，即主坐标计算为：

其中，a

其中，N

3)

其中，

S5、通过插值主坐标

1)基于插值点处的

2)然后可近似重构出在原始时间和频率坐标处的H

在上述步骤S5的步骤1)和步骤2的计算中，POD仅仅被执行于插值点上，而插值技术也只对少量时间函数和频率函数使用。因此，所提方法比现有的时频解耦方法更高效。

S6、基于谱表示法基本模拟公式，采用快速傅里叶变换(FFT)生成模拟样本x

1)基于谱表示方法，任意随机过程的模拟公式为：

其中，t是时间坐标；ω

2)将解耦表达式代入上式(12)中并进行数学等价变化得：

3)该模拟公式(12)中的三角函数求和项即可采用FFT技术进行计算，从而快速地生成多点风场样本。

以下以模拟一个水平分布的非平稳风场为例对本发明的方法进行进一步说明：

(1)模拟点为沿着960m水平直线均匀分布的100个点，它们的高度均为60m。

(2)假定各点的自演化功率谱密度函数为：

u

其中，z

U(z

其中，U(z

考虑随时间变化的相干函数，取为：

其中，x

(3)基于上述定义的风速演化功率谱密度函数，确定时间-频率插值点

其中，ω

(4)根据上述定义的自演化功率谱与相干函数，即可以确定互演化功率谱密度函数，从而得到时频插值点处的演化功率谱矩阵

(5)在每一个时频插值点处对演化功率谱密度矩阵

(6)对于

其中，

(7)基于插值点处的

(8)随机生成在[0,2π]范围内均匀分布的随机相位角φ

其中，N为频率步数，取1024，Δω＝ω

通过上述模拟方法，模拟出风场中100个点的脉动风速时程。通过改变随机相位角，模拟2000组这样的风场，并对其自演化功率谱、自相关函数和互相关函数与定义的相应目标值进行比对。图2和图3分别给出了模拟风速的相关函数与演化功率谱密度函数的验证结果，由图可知，模拟值与目标值保持一致，从而可以说明模拟方法的合理性。

进一步，改变风场中模拟点的数目分别为：

50，100，150，200，250，300，350，400，450，500，统计每一种工况下模拟一个风场样本所需要的耗时。图4展示了相应的统计结果，可以看到模拟方法的计算效率相当高，对于500个模拟点的风场模拟也仅仅耗时200秒左右。

在本发明的上述实施例中，本发明的方法能够解决非平稳风场模拟中时间和频率点数多导致需要进行Cholesky分解的次数多，计算量大，需要的电脑储存空间大的问题；此外，采用本发明的方法模拟非平稳随机风场，不需要传统方法中的在每个时间-频率点处进行Cholesky分解，仅用少量的插值点来代替原始时间-频率点，从而能够大大减少了计算量，提高了模拟效率；同时，本发明的方法可为大跨度结构，如桥梁、屋盖等抗风设计的时程分析提供输入的风荷载时程。

本具体实施例仅仅是对本发明的解释，其并不是对本发明的限制，本领域技术人员在阅读完本说明书后可以根据需要对本实施例做出没有创造性贡献的修改，但只要在本发明的权利要求范围内都受到专利法的保护。