一种基于KL展开的分布随机动载荷识别方法与流程

技术领域：

本发明涉及一种基于kl展开的分布随机动载荷识别方法，属于结构动力学反问题技术领域。

背景技术：
：

工程结构上的动载荷信息是结构设计和安全评估的依据。许多情况下，动载荷难以通过直接测量获得，常直接测量结构上的动响应，通过间接识别的手段获取结构上的动载荷信息。

传统的动载荷识别方法是利用单次实测的结构动响应数据识别引起该次动响应的激励信息，是确定性动载荷识别方法。现有的确定性动载荷识别方法被用于获取工程结构上的集中动载荷，移动载荷以及分布式动载荷等信息。值得注意的是，分布式动载荷识别问题相当于识别无穷多个集中动载荷，难度更大，一般需要将分布式动载荷识别问题降维求解。

实际工程结构上作用的动载荷，如建筑物上的风载荷，海洋平台承受的海浪载荷以及飞行器表面的气动载荷等，不仅分布于结构上，而且还具有随机性。随机动载荷施加于结构时，动响应也将随之呈现“随机性”；因此，单次实测的结构动响应只能是结构随机动响应信息的其中一个样本，基于某个响应样本利用确定性动载荷识别方法获得的动载荷信息也只能部分反映该随机动载荷激励；另外，单次测量中包含的动响应误差在确定性动载荷识别中也被作为“真实响应”的一部分，引起载荷识别结果的偏差。针对此类分布式随机动载荷的识别问题，传统的确定性分布动载荷识别方法和集中随机动载荷的识别方法均无法适用，需要发展一种针对分布式随机动载荷识别的新方法。

技术实现要素：

本发明的目的是提供一种基于kl展开的分布随机动载荷识别方法，解决在时域内利用实测结构动响应样本识别结构上随机动载荷随空间分布的时变统计特征问题，为服役于分布随机动载荷环境下的工程结构设计与安全评估提供一种动载荷间接获取手段。

上述的目的通过以下技术方案实现：

一种基于kl展开的分布随机动载荷识别方法，该方法包括如下步骤：

s1.开展模态试验，获取结构的模态参数，包括固有频率和模态振型；

s2.将结构随机振动响应利用模态振型展开，获取结构在模态空间的动响应；

s3.利用kl展开由模态空间内的随机动响应求解模态空间内随机动载荷；

s4.求解结构上随机动载荷的随空间分布的时变统计特征。

所述的基于kl展开的分布随机动载荷识别方法，步骤s2中所述的将结构随机振动响应利用模态振型展开，获取结构在模态空间的动响应，具体步骤为：

s21：利用多次重复测量方式获取随机振动响应的样本集合；

s22：针对单个样本，即单次实测结构振动响应，利用模态振型展开获得结构振动在模态空间内响应，具体步骤如下：

在结构表面上(x1,y1)，(x1,y2)，…，(xn,yn)等位置处第r次测量获取位移响应样本向量wr表示为：

wr＝{wr(x1,y1,t)wr(x1,y2,t)…wr(xn,yn,t)}^t,r＝1,…,n(1)，

其中wr(xi,yj,t)表示第r次测量获得的结构上空间位置(xi,yj)处t时刻的动位移值，n为测量的次数，利用模态振型函数计算模态空间内第r次测量对应的模态位移向量：

其中qi,r(t)为第r次测量获得的结构位移在第i阶模态空间中的模态位移，表示第i阶模态振型函数在(xj,yk)处的值，右上角+号表示广义逆。

所述的基于kl展开的分布随机动载荷识别方法，步骤s3中所述的利用kl展开由模态空间内的随机动响应求解模态空间内随机动载荷，具体包括以下步骤：

s31：利用模态空间内随机动响应的样本集合qi,r(t)，求解第i阶模态空间内随机位移响应的协方差矩阵γqi；

s32：对协方差矩阵进行特征值分解γqi，计算其特征值λij和特性向量ηij(t)，进一步获取第i阶模态空间内随机位移响应的第j个kl向量zij(t)，可以表示为：

由此，第i阶模态空间内随机位移响应qi(t,θ)的kl展开可以表示为：

其中ξj(θ)为随机变量，θ表示随机维度，当j＝0时ξ0(θ)＝1；

s33：根据下式，由随机动位移的kl向量zij(t)反演随机动载荷对应向量uij(t)；

其中ωi，ζi和mi分别为第i阶固有频率，模态阻尼比和模态质量，和分别是zij(t)对时间t的一阶和两阶导数；

s34：由随机动载荷对应向量uij(t)求解模态空间内随机动载荷fi(t,θ)，

所述的基于kl展开的分布随机动载荷识别方法，步骤s4中所述的求解结构上随机动载荷的随空间分布的时变统计特征，包括以下步骤：

分布随机动载荷f(x,y,t,θ)的表达式为：

结构上随机动载荷的随空间分布的时变统计特征，包括均值μf(x,t)和方差varf(x,t)分别为：

有益效果：

本发明与现有技术相比，具有以下优点：

1、现有的随机动载荷识别技术一般只能由实测结构动响应样本识别结构上随机集中动载荷，目前已经出现的分布随机动载荷识别方法大多无法适用于非平稳随机动载荷的识别；本发明中提供的基于kl展开的分布随机动载荷时域识别技术能够利用有限测点处的实测结构动响应样本识别随机动载荷随空间分布的统计特征，具有一定的先进性；

2、利用kl展开在模态空间内由随机动响应反演随机动载荷，比基于随机样本的蒙特卡洛法具有更高的计算效率，具有一定的优势。

附图说明

图1为本发明方法的逻辑流程框图。

图2为分布随机载荷作用下简支梁示意图。

图3(a)为梁跨中处随机动载荷均值识别结果。

图3(b)为梁跨中处随机动载荷方差识别结果。

图4为梁上随机动载荷空间分布识别结果。

具体实施方式

下面通过实施例的方式，对本发明技术方案进行详细说明，但实施例仅是本发明的其中一种实施方式，应当指出：对于本技术领域的技术人员来说，在不脱离本发明原理的前提下，还可以以改变结构和载荷形式等方式做出若干改进和等同替换，这些对本发明权利要求进行改进和等同替换后的技术方案，均落入本发明的保护范围。

实施例：对如图2所示一维简支梁上作用分布随机动载荷工况，利用本发明的技术由实测结构随机动响应样本识别随机动载荷随空间分布的时变统计特征。实施例中简支梁长l＝40m，横截面积a＝4.8m²，截面惯性矩i＝2.5498m⁴，结构的阻尼采用瑞利阻尼，各阶模态阻尼比ξi＝0.02，材料的弹性模量e＝5×10¹⁰n/m²，密度ρ＝2.5×10³kg/m³。待识别的梯形分布随机动载荷分布函数为：

分布式随机动载荷的随机性动载荷组分f(t,θ)分为确定性动载荷和随机性动载荷两个部分。

f(t,θ)的确定性动载荷部分：

fd(t)＝20000[1+0.1sin(2πt)]n(2)

f(t,θ)的随机性动载荷部分假定为零均值非平稳高斯随机过程，功率谱函数s(ω,t)为：

s(ω,t)＝cfpd(t)φ(ω)(3)

其中：cf表示随机水平，取cf＝0.2；φ(ω)表示零均值非平稳高斯随机过程的功率谱密度函数，有φ(ω)＝(1/2π)(2/ω²+1)。

具体包括以下步骤：

s1：获取结构的前五阶固有频率分别为3.9hz,15.6hz,35.1hz,62.5hz和97.6hz，同时获得各阶固有频率所对应的模态振型；

s2：将结构随机振动响应利用模态振型展开，利用多次测量获取位移响应信号，求解结构在模态空间的动响应，包括以下步骤：

s21：在梁结构上均匀布置19个测点，重复测量各测点处的动态位移信号，获取在分布随机激励下随机振动响应的样本集合；

s22：针对单个样本，即单次实测结构振动响应，利用模态振型展开获得结构振动在模态空间内响应。具体步骤如下：

在梁结构上(x1,x2,…xn)位置处第r次测量获取位移响应样本向量wr表示为：

wr＝{wr(x1,t)wr(x2,t)…wr(xn,t)}^t,r＝1,…,n(4)，

其中wr(xj,t)表示第r次测量获得的结构位移在xj处t时刻的值，n为测量的次数，即为样本总数。利用模态振型函数计算模态空间内第r次测量对应的模态位移向量：

其中qi,r(t)为第r次测量获得的结构位移在第i阶模态空间中的模态位移，表示第i阶模态振型函数在xj处的值，右上角+号表示广义逆。此时，测点数n＝19，模态数m＝5。

s3.利用kl展开由模态空间内的随机动响应求解模态空间内随机动载荷，包括以下步骤：

s31：利用模态空间内随机动响应的样本集合qi,r(t)，求解第i阶模态空间内随机位移响应的协方差矩阵γqi；

s32：对协方差矩阵进行特征值分解γqi，计算其特征值λij和特性向量ηij(t)。进一步获取第i阶模态空间内随机位移响应的第j个kl向量zij(t)，可以表示为：

kl展开中kl向量的截断原则为保证截断后的k-l向量组分在最小均方误差准则下接近原向量，具体根据下式：

其中γ＝0.99。

由此，第i阶模态空间内随机位移响应qi(t,θ)的kl展开可以表示为：

其中ξj(θ)为随机变量，θ表示随机维度，当j＝0时ξ0(θ)＝1。

s33：根据下式，由随机动位移的kl向量zij(t)反演随机动载荷对应向量uij(t)；

其中ωi，ζi和mi分别为第i阶固有频率，模态阻尼比和模态质量，和分别是zij(t)对时间t的一阶和两阶导数。简支梁的模态质量mi可由下式计算：

s34：由随机动载荷对应向量uij(t)求解模态空间内随机动载荷fi(t,θ)，

s4：求解结构上随机动载荷的随空间分布的时变统计特征，包括以下步骤：

分布随机动载荷f(x,t,θ)的识别结果可以表示为：

结构上随机动载荷的随空间分布的时变统计特征，包括均值μf(x,t)和方差varf(x,t)可以分别由下面公式得到：

图3(a)中给出了识别获得的梁跨中处随机动载荷均值随时间变化规律与真实值的对比；图3(b)中给出了识别获得的梁跨中处随机动载荷方差随时间变化规律与真实值的对比；图4中给出了识别获得的各时刻梁上随机动载荷的空间分布与真实分布的对比结果。由此可知，本发明中的识别方法能够利用有限测点处响应样本准确识别随机动载荷随空间的分布以及随时间变化的统计特征，适用于非平稳随机动载荷的情况；同时，跟基于样本的蒙特卡洛法相比，当实测响应样本数量较多时，在计算效率上有明显的优势，例如，在本实施例中，实测样本数等于5000时，在同等识别精度下，基于kl展开的识别方法所用计算时间仅为基于蒙特卡洛法的15％，大幅提高了识别效率。综上所述，本发明提出的方法具有一定的先进性。