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

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


背景技术：

2.在对大跨度桥梁、高层建筑等结构的风振响应分析中，由于需要考虑结构非线性和气动非线性，传统的频域分析方法已不能满足要求，需要采用非线性响应时程分析法。在非线性时程分析中，模拟结构所受的脉动风场是一个重要的环节。
3.在众多模拟方法中，谱表示法的形式简单且具有严密的理论推导，是风场模拟最常用的方法之一。为了提高多点风场的模拟效率，该方法常与快速傅里叶变换(fft)技术相结合。然而对于非平稳脉动风场的模拟，为了引入fft，需要先对时频耦合的函数进行解耦，目前已有包括本征正交分解(pod)、小波分解、非负矩阵分解等处理手段，但是解耦计算量会随着模拟点数的增多而急剧增大。此外，若所模拟风场的相干函数具有时变特性，那么还需要大量的谱矩阵分解计算，这使得非平稳风场的模拟效率进一步大打折扣。
4.因此需要提出一种更为高效的非平稳风场谱表示模拟方法，以解决上述问题。


技术实现要素：

5.本发明的目的是为了解决上述技术问题，提供一种基于pod插值的非平稳风场高效模拟方法，本发明的方法能够解决非平稳风场模拟中时间和频率点数多导致需要进行cholesky分解的次数多，计算量大，需要的电脑储存空间大的问题，本发明的方法用极少量的时间
‑
频率插值点代替原始时间
‑
频率点，不仅能够满足模拟的精度也能提高模拟效率。
6.本发明的上述技术目的是通过以下技术方案得以实现的：一种基于pod插值的非平稳风场高效模拟方法，具体包括以下步骤：
7.s1、基于给定的风速演化功率谱密度函数s(ω,t)，定义时间
‑
频率插值点
8.s2、基于给定的自演化功率谱密度函数和相干函数，计算时频插值点处的演化功率谱密度矩阵
9.s3、分解时频插值点处的演化功率谱密度矩阵获得相应的分解后谱矩阵
10.s4、通过本征正交分解(pod)解耦
11.s5、通过插值主坐标获得通过插值特征向量获得
12.s6、基于谱表示法基本模拟公式，采用快速傅里叶变换(fft)生成模拟样本x
j
(t)。
13.进一步地，步骤s1的具体方法为：
14.(1)沿时间方向的插值点采用均匀分布，将区间[0,t0]均匀分割，确定时间插值点
坐标为：
[0015][0016]
其中，t0为演化功率谱密度总时长；为时间插值点数且m是目标演化功率谱及模拟时程的时间点数；
[0017]
(2)采用一种可调四次多项式来描述非均匀分布的频率插值点，即非均匀频率插值点坐标为：
[0018][0019]
其中，ω1和ω
u
分别是第一个频率点和最后一个频率点；κ是调节非均匀分布的参数，当κ＝0.25时，表示均匀分布；是频率插值点数且n是目标演化功率谱的频率点数。
[0020]
进一步地，步骤s2的具体方法为：时频插值点处的演化功率谱密度矩阵由下式(3)求得：
[0021][0022]
其中，s
jj
(ω,t)是第j点的自演化功率谱密度函数；s
jk
(ω,t)是第j点和第k点的互演化功率谱密度函数；γ
jk
(ω,t)是第j点和第k点的相干函数。
[0023]
进一步地，步骤s3的具体方法为：
[0024]
在每一个时频插值点处对演化功率谱密度矩阵进行cholesky分解可得，分解公式如下：
[0025][0026]
此过程总共需要次cholesky分解操作，求得的分解后谱矩阵是下三角矩阵，可表示为：
[0027][0028]
进一步地，步骤s4的具体方法为：对于的每一个下三角元素，即分别进行如下的pod解耦：
[0029]
1)将所有时频插值点处的取值表示为一个确定的常数矩阵：
[0030][0031]
并将该常数矩阵的每一列看作一个列向量，使用pod找到一组最优标准正交基使得列向量的投影最大化，且这组正交基通过特征向量分解来确定：
[0032][0033]
其中，r为时间平均的频率相关矩阵；φ
q
是第q个特征向量；λ
q
是第q个特征值；
[0034]
2)然后，列向量在特征向量上的投影，即主坐标计算为：
[0035][0036]
其中，a
q
是第q个投影向量；然后，将特征值按降序排列，选取特征值更小的少量模态，近似重构l为：
[0037][0038]
其中，n
r
为所选低阶项的个数；
[0039]
3)近似表示为几个低阶特征向量与对应主坐标的乘积之和，将常数矩阵表达的重构式还原为函数形式，即：
[0040][0041]
其中，为主坐标；为特征向量。
[0042]
进一步地，步骤s5的具体方法为：
[0043]
1)基于插值点处的和采用插值技术，近似求得其他时频坐标处的和
[0044]
2)然后可近似重构出在原始时间和频率坐标处的h
jk
(ω,t)，即
[0045][0046]
进一步地，步骤s6的具体方法为：
[0047]
1)基于谱表示方法，任意随机过程的模拟公式为：
[0048][0049]
其中，t是时间坐标；ω
l
是第l个频率点；δω是频率间隔；e是指数函数；i是虚数单位；φ
kl
是一组随机相位角，对于具体的样本为确定量；re表示获取复数的实部；
[0050]
2)将解耦表达式代入上式(12)中并进行数学等价变化得：
[0051][0052]
3)然后，采用fft技术对模拟公式(12)中的三角函数求和项进行计算，快速地生成多点风场样本。
[0053]
综上所述，本发明具有以下有益效果：
[0054]
1、本发明的方法能够解决非平稳风场模拟中时间和频率点数多导致需要进行cholesky分解的次数多，计算量大，需要的电脑储存空间大的问题；
[0055]
2、采用本发明的方法模拟非平稳随机风场，通过少量的插值点来代替原始时间
‑
频率点，不需要传统方法中的在每个时间
‑
频率点处进行cholesky分解，而且仅对时频插值点处的函数进行解耦，从而能够大大减少了计算量，提高了模拟效率；
[0056]
3、本发明的方法可为大跨度结构，如桥梁、屋盖等抗风设计的时程分析提供输入的风荷载时程。
附图说明
[0057]
图1是本发明实施例中的流程图；
[0058]
图2是本发明实施例中模拟风速的相关函数验证曲线图；
[0059]
图3是本发明实施例中模拟风速的演化功率谱密度函数验证曲线图；
[0060]
图4是本发明实施例中不同模拟点数对应的耗时曲线图。
具体实施方式
[0061]
以下结合附图1
‑
4对本发明作进一步详细说明。
[0062]
实施例：本发明的技术方案提供一种基于pod插值的非平稳风场高效模拟方法，如图1所示，具体包括以下步骤：
[0063]
s1、基于给定的风速演化功率谱密度函数s(ω,t)，定义时间
‑
频率插值点具体方法为：
[0064]
(1)由于目标演化功率谱随时间变化缓慢，则沿时间方向的插值点采用均匀分布，将区间[0,t0]均匀分割，可以容易地确定时间插值点坐标为：
[0065][0066]
其中，t0为演化功率谱密度总时长；为时间插值点数且m是目标演化功率谱及模拟时程的时间点数；
[0067]
(2)考虑到目标演化功率谱在低频处随频率变化剧烈，而在高频处变化缓慢，则采
用一种可调四次多项式来描述非均匀分布的频率插值点，即非均匀频率插值点坐标为：
[0068][0069]
其中，ω1和ω
u
分别是第一个频率点和最后一个频率点；κ是调节非均匀分布的参数，当κ＝0.25时，表示均匀分布；是频率插值点数且n是目标演化功率谱的频率点数。
[0070]
s2、基于给定的自演化功率谱密度函数和相干函数，计算时频插值点处的演化功率谱密度矩阵具体方法为：时频插值点处的演化功率谱密度矩阵由下式(3)求得：
[0071][0072][0073]
其中，s
jj
(ω,t)是第j点的自演化功率谱密度函数；s
jk
(ω,t)是第j点和第k点的互演化功率谱密度函数；γ
jk
(ω,t)是第j点和第k点的相干函数。
[0074]
s3、分解时频插值点处的演化功率谱密度矩阵获得相应的分解后谱矩阵具体方法为：
[0075]
在每一个时频插值点处对演化功率谱密度矩阵进行cholesky分解可得，分解公式如下：
[0076][0077]
此过程总共需要次cholesky分解操作，求得的分解后谱矩阵是下三角矩阵，可表示为：
[0078][0079]
s4、通过本征正交分解(pod)解耦具体方法为：对于的每一个下三角元素，即分别进行如下的pod解耦：
[0080]
1)将所有时频插值点处的取值表示为一个确定的常数矩阵：
[0081][0082]
并将该常数矩阵的每一列看作一个列向量，使用pod找到一组最优标准正交基使得列向量的投影最大化，且这组正交基通过特征向量分解来确定：
[0083][0084]
其中，r为时间平均的频率相关矩阵；φ
q
是第q个特征向量；λ
q
是第q个特征值；
[0085]
2)然后，列向量在特征向量上的投影，即主坐标计算为：
[0086][0087]
其中，a
q
是第q个投影向量；然后，将特征值按降序排列，选取特征值更小的少量模态，近似重构l为：
[0088][0089]
其中，n
r
为所选低阶项的个数；
[0090]
3)近似表示为几个低阶特征向量与对应主坐标的乘积之和，将常数矩阵表达的重构式还原为函数形式，即：
[0091][0092]
其中，为主坐标；为特征向量。
[0093]
s5、通过插值主坐标获得通过插值特征向量获得s5的具体方法为：
[0094]
1)基于插值点处的和采用插值技术，近似求得其他时频坐标处的和
[0095]
2)然后可近似重构出在原始时间和频率坐标处的h
jk
(ω,t)，即
[0096][0097]
在上述步骤s5的步骤1)和步骤2的计算中，pod仅仅被执行于插值点上，而插值技术也只对少量时间函数和频率函数使用。因此，所提方法比现有的时频解耦方法更高效。
[0098]
s6、基于谱表示法基本模拟公式，采用快速傅里叶变换(fft)生成模拟样本x
j
(t)，
具体方法为：
[0099]
1)基于谱表示方法，任意随机过程的模拟公式为：
[0100][0101]
其中，t是时间坐标；ω
l
是第l个频率点；δω是频率间隔；e是指数函数；i是虚数单位；φ
kl
是一组随机相位角，对于具体的样本为确定量；re表示获取复数的实部；
[0102]
2)将解耦表达式代入上式(12)中并进行数学等价变化得：
[0103][0104]
3)该模拟公式(12)中的三角函数求和项即可采用fft技术进行计算，从而快速地生成多点风场样本。
[0105]
以下以模拟一个水平分布的非平稳风场为例对本发明的方法进行进一步说明：
[0106]
(1)模拟点为沿着960m水平直线均匀分布的100个点，它们的高度均为60m。
[0107]
(2)假定各点的自演化功率谱密度函数为：
[0108][0109]
u
*
(z
j
,t)＝ku(z
j
,t)/ln(z
j
/z0)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(15)
[0110][0111]
其中，z
j
是第j点的高度，u
*
(z
j
,t)是剪切波速，k＝0.4，z0＝0.01m为地面粗糙长度。u(z
j
,t)是在高度z
j
处的平均风速，被定义为：
[0112]
u(z
j
,t)＝d(t)u(z
j
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(17)
[0113][0114]
其中，u(z
j
)为平均风速取u(z
j
)＝40m/s，d(t)为时间调制函数，参数取值为t
max
＝β0/λ,t
max
＝600,β0＝4。
[0115]
考虑随时间变化的相干函数，取为：
[0116][0117]
其中，x
j
是第j点的水平坐标，c
x
＝7。
[0118]
(3)基于上述定义的风速演化功率谱密度函数，确定时间
‑
频率插值点各点坐标为
[0119][0120]
[0121]
其中，ω
u
取πrad/s，ω1取0.0031rad/s，t0取2048s，取20，取15，κ取0.8。
[0122]
(4)根据上述定义的自演化功率谱与相干函数，即可以确定互演化功率谱密度函数，从而得到时频插值点处的演化功率谱矩阵即：
[0123][0124][0125]
(5)在每一个时频插值点处对演化功率谱密度矩阵进行cholesky分解，获得相应的分解后谱矩阵
[0126][0127][0128]
(6)对于的每一个下三角元素即的每一个下三角元素即分别进行pod解耦可得：
[0129][0130]
其中，为主坐标；为特征向量；n
r
取4。
[0131]
(7)基于插值点处的和采用插值技术，可近似求得其他时频坐标处的和
[0132]
(8)随机生成在[0,2π]范围内均匀分布的随机相位角φ
jk
，带入上述结果到以下模拟公式并借助fft生成各点的风速时程：
[0133][0134]
其中，n为频率步数，取1024，δω＝ω
u
/n为频率增量，ω
l
＝lδω为频率点坐标。
[0135]
通过上述模拟方法，模拟出风场中100个点的脉动风速时程。通过改变随机相位角，模拟2000组这样的风场，并对其自演化功率谱、自相关函数和互相关函数与定义的相应目标值进行比对。图2和图3分别给出了模拟风速的相关函数与演化功率谱密度函数的验证结果，由图可知，模拟值与目标值保持一致，从而可以说明模拟方法的合理性。
[0136]
进一步，改变风场中模拟点的数目分别为：
[0137]
50，100，150，200，250，300，350，400，450，500，统计每一种工况下模拟一个风场样
本所需要的耗时。图4展示了相应的统计结果，可以看到模拟方法的计算效率相当高，对于500个模拟点的风场模拟也仅仅耗时200秒左右。
[0138]
在本发明的上述实施例中，本发明的方法能够解决非平稳风场模拟中时间和频率点数多导致需要进行cholesky分解的次数多，计算量大，需要的电脑储存空间大的问题；此外，采用本发明的方法模拟非平稳随机风场，不需要传统方法中的在每个时间
‑
频率点处进行cholesky分解，仅用少量的插值点来代替原始时间
‑
频率点，从而能够大大减少了计算量，提高了模拟效率；同时，本发明的方法可为大跨度结构，如桥梁、屋盖等抗风设计的时程分析提供输入的风荷载时程。
[0139]
本具体实施例仅仅是对本发明的解释，其并不是对本发明的限制，本领域技术人员在阅读完本说明书后可以根据需要对本实施例做出没有创造性贡献的修改，但只要在本发明的权利要求范围内都受到专利法的保护。