一种基于自相关函数的结构风振响应高效频域估计方法技术

一种基于自相关函数的结构风振响应高效频域估计方法，本发明专利技术涉及基于自相关函数的结构风振响应高效频域估计方法。本发明专利技术的目的是为了解决现有风振响应分析方法计算效率低以及计算误差大的缺点。具体过程为：一：基于pi(t)计算σpi、Rpi(τ)和Cohpij(ω)；二：得到ωmi；三：得到该结构的相干指数kc，计算i、j两点的风荷载相干系数；四：提取[M]、[K]、[C]、[R]，[I]；五：根据一、二、三和四计算解析积分得到[Σkr]；六：形成模态响应协方差矩阵；七：计算位移响应协方差矩阵，进而可求得其任意响应的协方差矩阵。本发明专利技术用于结构风阵响应领域。

【技术实现步骤摘要】

本专利技术涉及基于自相关函数的风振响应高效频域估计方法。
技术介绍
风振响应分析是大型复杂结构(例如大跨度屋盖结构、超高层建筑等)抗风设计的重要环节之一。大型复杂结构的振动特性难以用简单的振型表示，要考虑高阶振型的贡献以及振型之间的耦合，因此，脉动风荷载作用下的结构动力响应极为复杂。传统的时域风振响应分析，需要输入大量风荷载时程信息，采用Newmark-β积分法，计算量较大，计算时所需存储空间较大，导致计算效率低。而传统的频域风振响应分析，同样需要输入大量信息，包括风荷载的协方差矩阵，采用数值积分法，计算耗时较长，导致计算效率低。此外，传统频域计算方法难以考虑振型耦合效应，计算误差大。
技术实现思路
本专利技术的目的是为了解决现有风振响应分析方法计算效率低以及计算误差大的缺点，而提出一种基于自相关函数的结构风振响应高效频域估计方法。一种基于自相关函数的结构风振响应高效频域估计方法具体过程为：步骤一：基于风洞试验测得的结构表面脉动风压时程数据pi(t)，计算各风荷载加载点的风压时程数据的平均值风压时程数据的标准差σpi、风压时程数据的自相关函数Rpi(τ)和风压时程数据的相干函数Cohpij(ω)；其中，i、j为风荷载加载点总数Q中任意两个点，1≤i≤Q，1≤j≤Q，Q为风荷载加载点总数，为正整数；步骤二：采用指数函数拟合自相关函数，得到i点的风荷载加载点的无量纲风压谱曲线的峰值频率ωmi；由自功率谱的定义可知Spi(ω)σpi2=1π∫-∞∞Rpi(τ)exp(-iωτ)dτ=2π·ωmiωmi2+ω2]]>式中，Spi(ω)为各加载点的风压自功率谱密度函数；ω为频率；τ为时差；则无量纲风压谱表示为可得出，式中，Si(ω)为无量纲风压谱；步骤三：采用指数函数拟合相干函数，得到该结构的相干指数kc，计算i、j两点的风荷载相干系数；步骤四：提取结构的质量矩阵[M]、刚度矩阵[K]、阻尼矩阵[C]、风荷载加载点与结构自由度的转换矩阵[R]，响应的影响面矩阵[I]；对结构进行模态分析，求解特征方程得到结构的各阶自振频率ωnk，k＝1,2,…,N，k为自振频率所对应的模态阶数，N为结构自由度数，为正整数；取结构振型矩阵其中为第k、r阶振型，k,r＝1,2,…W；W为计算所选取的模态阶数，W为正整数，10≤W≤N；k、r为自振频率ωnk、ωnr所对应的模态阶数；使得计算结构k阶振型的广义阻尼比式中，T为转置矩阵；步骤五：根据步骤一、步骤二、步骤三和步骤四计算解析积分；得到[Σkr]；步骤六：计算模态响应协方差形成模态响应协方差矩阵；步骤七：计算位移响应协方差矩阵[Σx]＝[Ψ]T[Σy][Ψ]，进而可求得其任意响应的协方差矩阵。本专利技术的有益效果为：本方法设计基于自相关函数的大型复杂结构风振响应高效频域估计方法，解决传统风振响应分析方法在大型复杂结构分析中计算效率低以及计算精度较差的缺点。该方法首先将风洞试验得到的脉动风荷载处理为供计算用的统计量和频谱参数信息，大大简化了风荷载的输入。然后，结合振型分解法，对每个振型，在每个自由度上对风荷载频谱进行简化的解析积分，得到模态响应协方差。最后，通过对耦合振型的组合得到风振响应的协方差。该方法在保证模态耦合计算精度的同时大大提高了计算效率。本专利技术在保证计算精度的前提下，大大提升了大型复杂结构风振计算的效率，降低了计算所占用的存储空间。充分利用了解析积分的计算原理，将数值积分转化为解析积分的代数运算，避免了数值积分带来的误差，同时也大大提高了计算效率。在实施过程中，步骤明确，可操作性较强，通过实例分析发现，本专利技术较传统算法效率大为提升，较为实用。采用该专利技术对某扇形平面大跨度悬挑网架结构进行风振响应分析，与传统分析方法相比，最大误差不超过1％，远优于传统的简化频域算法误差(16.7％)，效率较传统算法提高200倍。附图说明图1为本专利技术的流程图；图2a为取-6.42，σpi取2.48的结构表面某点风荷载时程按步骤一统计建模的结果示意图，纵坐标为为归一化的风荷载时程，横坐标为t，为时间(秒)，pi(t)为风荷载时程，pi(t)为平均风荷载，σpi为风荷载标准差；图2b为Rpi(τ)＝exp(-ωmi|τ|)，ωmi＝2.40rad/s的结构表面某点风荷载时程按步骤一统计建模的结果示意图，纵坐标为Rpi(τ)，Rpi(τ)为风荷载的自相关函数，ωmi为风荷载频率，横坐标为τ，τ为时差(秒)；图2c为的结构表面某点风荷载时程按步骤一统计建模的结果示意图，纵坐标为为无量纲风荷载功率谱，横坐标为ω，ω为频率(弧度/秒)，Spi(ω)为风荷载功率谱密度函数；图3a为扇形平面悬挑网架结构的有限元模型示意图；图3b为扇形平面悬挑网架结构按步骤二进行模态分析给出的一阶振型(ωn1＝7.4rad/s)结果示意图；图3c为扇形平面悬挑网架结构按步骤二进行模态分析给出的二阶振型(ωn2＝7.5rad/s)结果示意图；图3d为扇形平面悬挑网架结构按步骤二进行模态分析给出的二阶振型(ωn3＝9.9rad/s)结果示意图；图4a为步骤三计算的模态响应均方根结果示意图，RMSmodaldisp为模态响应均方根值，modalorder为模态阶数，frequency为结构自振频率，条形图为模态响应标准差，圆圈为结构各阶模态的频率；图4b为步骤四计算的模态响应协方差矩阵结果示意图；图5是本专利技术的算法与传统算法在精度和效率上的比较。具体实施方式具体实施方式一：结合图1说明本实施方式，本实施方式的一种基于自相关函数的结构风振响应高效频域估计方法具体过程为：步骤一：基于风洞试验测得的结构(例如大跨度屋盖结构(如两端支承钢结构屋盖跨度大于36m，悬挑长度大于10m的结构，例如多数普通体育场、机场候机楼、会展中心展览馆等的屋盖结构)、超高层建筑(如高度超过100m的高层结构)，)表面脉动风压时程数据pi(t)，计算各风荷载加载点的风压时程数据的平均值风压时程数据的标准差σpi、风压时程数据的自相关函数Rpi(τ)和风压时程数据的相干函数Cohpij(ω)；其中，i、j为风荷载加载点总数Q中任意两个点，1≤i≤Q，1≤j≤Q，Q为风荷载加载点总数，为正整数；步骤二：采用指数函数拟合自相关函数，得到i点的风荷载加载点的无量纲风压谱曲线的峰值频率ωmi；由自功率谱的定义可知Spi(ω)σpi2=1π∫-∞∞Rpi(τ)exp(-iωτ)dτ=2π·ωmiωmi2+ω2]]>式中，Spi(ω)为各加载点的风压自功率谱密度函数；ω为频率；τ为时差；则无量纲风压谱表示为可得出，式中，Si(ω)为无量纲风压谱；步骤三：采用指数函数拟合相干函数，得到该结构的相干指数kc，计算i、j两点的风荷载相干系数；步骤四：提取结构的质量矩阵[M]、刚度矩阵[K]、阻尼矩阵[C]、风荷载加载点与结构自由度的本文档来自技高网...

【技术保护点】
一种基于自相关函数的结构风振响应高效频域估计方法，其特征在于：一种基于自相关函数的结构风振响应高效频域估计方法具体过程为：步骤一：基于风洞试验测得的结构表面脉动风压时程数据pi(t)，计算各风荷载加载点的风压时程数据的平均值风压时程数据的标准差σpi、风压时程数据的自相关函数Rpi(τ)和风压时程数据的相干函数Cohpij(ω)；其中，i、j为风荷载加载点总数Q中任意两个点，1≤i≤Q，1≤j≤Q，Q为风荷载加载点总数，为正整数；步骤二：采用指数函数拟合自相关函数，得到i点的风荷载加载点的无量纲风压谱曲线的峰值频率ωmi；步骤三：采用指数函数拟合相干函数，得到该结构的相干指数kc，计算i、j两点的风荷载相干系数；步骤四：提取结构的质量矩阵[M]、刚度矩阵[K]、阻尼矩阵[C]、风荷载加载点与结构自由度的转换矩阵[R]，响应的影响面矩阵[I]；对结构进行模态分析，求解特征方程得到结构的各阶自振频率ωnk，k＝1,2,…,N，k为自振频率所对应的模态阶数，N为结构自由度数，为正整数；取结构振型矩阵其中为第k、r阶振型，k,r＝1,2,…W；W为计算所选取的模态阶数，W为正整数，10≤W≤N；k、r为自振频率ωnk、ωnr所对应的模态阶数；使得计算结构k阶振型的广义阻尼比式中，T为转置矩阵；步骤五：根据步骤一、步骤二、步骤三和步骤四计算解析积分，得到[Σkr]；步骤六：计算模态响应协方差形成模态响应协方差矩阵[Σy]；步骤七：计算位移响应协方差矩阵[Σx]＝[Ψ]T[Σy][Ψ]，进而可求得其任意响应的协方差矩阵[Σs]。...

【技术特征摘要】
1.一种基于自相关函数的结构风振响应高效频域估计方法，其特征在于：一种基于自相关函数的结构风振响应高效频域估计方法具体过程为：步骤一：基于风洞试验测得的结构表面脉动风压时程数据pi(t)，计算各风荷载加载点的风压时程数据的平均值风压时程数据的标准差σpi、风压时程数据的自相关函数Rpi(τ)和风压时程数据的相干函数Cohpij(ω)；其中，i、j为风荷载加载点总数Q中任意两个点，1≤i≤Q，1≤j≤Q，Q为风荷载加载点总数，为正整数；步骤二：采用指数函数拟合自相关函数，得到i点的风荷载加载点的无量纲风压谱曲线的峰值频率ωmi；步骤三：采用指数函数拟合相干函数，得到该结构的相干指数kc，计算i、j两点的风荷载相干系数；步骤四：提取结构的质量矩阵[M]、刚度矩阵[K]、阻尼矩阵[C]、风荷载加载点与结构自由度的转换矩阵[R]，响应的影响面矩阵[I]；对结构进行模态分析，求解特征方程得到结构的各阶自振频率ωnk，k＝1,2,…,N，k为自振频率所对应的模态阶数，N为结构自由度数，为正整数；取结构振型矩阵其中为第k、r阶振型，k,r＝1,2,…W；W为计算所选取的模态阶数，W为正整数，10≤W≤N；k、r为自振频率ωnk、ωnr所对应的模态阶数；使得计算结构k阶振型的广义阻尼比式中，T为转置矩阵；步骤五：根据步骤一、步骤二、步骤三和步骤四计算解析积分，得到[Σkr]；步骤六：计算模态响应协方差形成模态响应协方差矩阵[Σy]；步骤七：计算位移响应协方差矩阵[Σx]＝[Ψ]T[Σy][Ψ]，进而可求得其任意响应的协方差矩阵[Σs]。2.根据权利要求1所述一种基于自相关函数的结构风振响应高效频域估计方法，其特征在于：所述步骤一中基于风洞试验测得的结构表面脉动风压时程数据pi(t)，计算各风荷载加载点的风压时程数据的平均值风压时程数据的标准差σpi、风压时程数据的自相关函数Rpi(τ)和风压时程数据的相干函数Cohpij(ω)；具体过程为：由MATLAB中的mean函数求解；σpi由MATLAB中的std函数求解；Rpi(τ)由MATLAB中的xcorr函数求解；Cohpij(ω)由MATLAB中的mscohere函数求解。3.根据权利要求2所述一种基于自相关函数的结构风振响应高效频域估计方法，其特征在于：所述步骤二中采用指数函数拟合自相关函数，得到i点的风荷载加载点的无量纲风压谱曲线的峰值频率ωmi；具体过程为：采用指数函数拟合自相关函数Rpi(τ)＝exp(-ωmi|τ|)，得到i点的风荷载加载点的无量纲风压谱曲线的峰值频率ωmi。4.根据权利要求3所述一种基于自相关函数的结构风振响应高效频域估计方法，其特征在于：所述步骤三中采用指数函数拟合相干函数，得到该结构的相干指数kc，计算i、j两点的风荷载相干系数；具体过程为：采用指数函数拟合相干函数得到该结构的相干指数kc，计算i、j两点的风荷载相干系数式中，Dij表示两风荷载加载点i、j间距离，U表示参考风速，cij为i、j两点的风荷载相干系数，为非负实数。5.根据权利要求4所述一种基于自相关函数的结构风振响应高效频域估计方法，其特征在于：所述步骤五中根据步骤一、步骤二、步骤三和步骤四计算解析积分，得到[Σkr]；具体过程为：根据步骤一、步骤二、步骤三和步骤四计算解析积分；σijkr=2πσpiσpjωmiωmj∫0∞[(ω2-ωnk2)(ω2-ωnr2)-4ξnkξnrωnkωnrω2]·(ω2+ωmiωmj)[(ω2-ωnk2)2+4ξnk2ξnk2ω2]·[(ω2-ωnr2)2+4ξnr2ωnr2ω2]·(ωmi2+ω2)·(ωmj2+ω2)·[1+(cijω)2]dω=2πσpiσpjωmiωmj∫0∞Θ(ω)Λ(--1ω)·Λ(-1ω)dω=σpiσpjωmiωmj·I1I0]]>式中，σijkr为k、r阶耦合模态在i、j两点间的耦合风振响应均方根；σpi为i...

【专利技术属性】
技术研发人员：孙瑛，苏宁，武岳，沈世钊，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江;23