空间分布动态载荷的等效识别方法技术

本发明专利技术公开一种结构所受空间分布动态载荷的等效识别方法。首先，利用本征正交分解将时间与空间耦合的分布动态载荷分解成一系列时间与空间独立的子载荷；其次，对结构响应进行测量，通过模态变换获取模态响应，并结合正则化方法对模态载荷进行准确稳定识别；再次，对识别的模态载荷进行盲源分离，获取有限个等效集中动态载荷的时间历程，并计算其混叠系数；最后，考虑混叠系数为载荷空间分布函数与结构模态振型的内积，对其进行稀疏分解，识别等效集中载荷的空间位置。本发明专利技术为难以测量的空间分布动态载荷提供了一种有效的等效识别方法，所识别的等效集中载荷施加于结构时，能产生与原空间分布动态载荷相同的输出响应。

An equivalent recognition method for dynamic load of spatial distribution

The invention discloses an equivalent recognition method for the dynamic load of the spatial distribution of the structure. First of all, by using the proper orthogonal decomposition, distributed dynamic load time and space coupling is decomposed into a series of time and space independent sub load; secondly, to measure the response of structure modal response by modal transformation, and combined with the regularization method for accurate and stable identification of modal loads; thirdly, the mode of load identification blind source separation, time limited access to a centralized equivalent dynamic load, and calculate the mixing coefficient; finally, considering the mixing coefficient of load inner product space distribution function and structure mode, on the sparse decomposition, space position recognition equivalent concentrated load. The invention provides an effective equivalent identification method for the difficult to measure spatial distribution dynamic load. When the equivalent concentrated load is applied to the structure, it can produce the same output response as the original space distribution dynamic load.

【技术实现步骤摘要】
空间分布动态载荷的等效识别方法
本专利技术属参数识别
，尤其涉及一种结构所受空间分布动态载荷的等效识别方法。
技术介绍
分布动态载荷广泛存在于实际工程中，如建筑物或飞行器所受的风载荷，堤坝或航船所受的浪载荷等。对于工程结构，分析研究其由分布动态载荷引起的振动问题，掌握振动机理，能够有效地防止或消除振动带来的不良影响，合理地利用振动有利的一面。然而，这需要以确定结构所受的分布动态载荷为前提。分布动态载荷既与时间历程相关，又在空间上呈现分布形式，由于受技术或经济条件的限制，其很难甚至根本无法通过直接测量对其确定。结构动态位移响应的测量相对容易和准确，因此利用有限个结构测量响应识别结构受到的分布动态载荷，将具有重要实际意义。
技术实现思路
针对分布动态载荷难以直接测量的困难，本专利技术提出一种基于结构动态响应的分布动态载荷等效识别方法，从而为结构振动控制、动力修改和动力学优化设计提供准确有效的外部载荷参量。本专利技术提供一种基于本征正交分解、模态分析、盲源分离及稀疏分解的空间分布动态载荷等效识别方法，包括如下步骤：1、一种空间分布动态载荷的等效识别方法，其特征在于，包括如下步骤：步骤1：将结构所受的未知分布动态载荷进行本征正交分解，转化为一系列时间与空间独立的子载荷；步骤2：建立结构有限元模型，并进行模态分析，获得结构的模态振型系数矩阵、模态质量、模态阻尼、模态刚度、无阻尼固有频率、阻尼固有频率及阻尼比；步骤3：利用传感器测量分布动态载荷作用下结构的动态响应，并进行模态变换获取相应的截断模态响应；步骤4：基于模态响应和模态格林核函数，利用正则化方法对有限阶次的模态载荷进行准确稳定识别；步骤5：将模态载荷等效为有限个集中动态载荷的混叠，利用盲源分离技术从识别的模态载荷中获取等效集中载荷的时间历程，并计算相应混叠系数；步骤6：考虑混叠系数为载荷空间分布与结构模态振型的内积，对混叠系数向量进行稀疏分解，识别等效集中载荷的空间作用位置。优选地，在所述步骤1中，待识别的空间分布动态载荷形式可以是时间和空间相互独立的，也可以时间和空间相互耦合的。当时间和空间相互耦合时，通过本征正交分解能将其转化为如式(1)所示的一系列时间与空间独立的子载荷。式中，F(x,t)表示作用于结构上的时间与空间耦合的分布动态载荷，通过本征正交分解可将其近似为L项时间和空间独立的子载荷之和，ψn(x)表示第n个子载荷的空间分布函数，sn(t)表示第n个子载荷的时间历程。优选地，在所述步骤3中，用于空间分布动态载荷识别的测量响应可以是位移、速度或加速度。优选地，在所述步骤4中，结构模态载荷反求的正向模型是基于格林核函数法建立的，并通过正则化提高识别结果的精度和稳定性。优选地，在所述步骤5中，各阶模态载荷的时间历程可等效为有限个集中载荷时间历程的混叠，其具有如下表达式，式中，pi(t)表示第i阶模态载荷，表示结构的第i阶模态振型，上标T表示向量或矩阵的转置，表示时间历程sn(t)的第i个混叠系数，其为第i阶模态振型向量与第n个空间分布函数向量的内积。优选地，在所述步骤6中，求解子载荷空间分布函数的式(7)为欠定方程组，故载荷空间分布函数的反求结果存在多解。将最稀疏解作为反求结果，即通过稀疏分解确定等效集中载荷的作用位置。式中，表示求解得到的对应第n个空间分布函数的混叠系数，H为截断模态数，N表示结构加载区域的离散点数，其通常远大于截断模态数H，(j＝1,2,…,N)表示结构第i阶模态下第j个离散点处的模态振型系数，ψnj表示第n个空间分布函数在第j个离散点处的等效值。本专利技术创新性地提出一种利用结构动态响应来实现分布动态载荷等效识别的方法。该方法将时间和空间耦合的分布动态载荷识别等效为有限个集中载荷时间历程的识别及其空间作用位置的确定，利用易于测量的结构动态响应准确稳定地识别出结构受到的等效集中载荷。根据上述技术方案，本专利技术的有益效果包括：(1)工程结构受到空间连续的分布动态载荷作用，受测试设备及经济技术条件的限制，常规的测试方法很难对分布动态载荷进行直接测量。本专利技术将分布动态载荷的识别等效为有限个集中载荷时间历程及作用位置的识别，对于实际工程中的结构设计、优化和评估具有十分重要的意义。(2)本专利技术在进行分布动态载荷时间历程识别时，充分利用了相对容易测量的结构动态响应信息及准确的结构模态信息，并结合了正则化和盲源分离技术，实现了集中载荷时间历程的稳定识别。(3)本专利技术在对分布载荷空间分布函数进行识别时，通过稀疏分解技术，求取了空间分布函数的最稀疏解，实现了集中载荷作用位置的确定。附图说明附图仅用于示出具体实施例的目的，而并不认为是对本专利技术的限制，在整个附图中，相同的参考符号表示相同的部件。图1是本专利技术中分布动态载荷等效识别方法的步骤流程图；图2是时间与空间耦合的分布动态载荷；图3是本征正交分解得到时间与空间独立的分布动态载荷；图4是受分布动态载荷作用的引擎盖内板；图5是识别的前6阶模态载荷；图6a是盲源分离得到的等效集中载荷时间历程；图6b是稀疏分解得到的等效集中载荷空间位置；图6c是识别的等效集中动态载荷；图7是等效载荷作用下结构动态响应及验证。具体实施方式下面结合附图1-7以等效识别某车引擎盖内板边缘受到的分布动态载荷为例，对本专利技术的具体实施方式进行详细说明。结构受到的分布动态载荷通常无法通过直接测量的手段来获取，而结构因分布动态载荷产生的响应却相对容易测量。因此，首先通过本征正交分解，将时间和空间耦合的分布动态载荷解耦成一系列时间与空间独立的子分布动态载荷；其次，利用测量响应反求结构的模态响应并在模态空间对有限阶模态载荷进行识别；再次，基于盲源分离技术获取有限个集中载荷的时间历程；最后通过稀疏分解确定集中载荷的空间作用位置。图1示出了本专利技术具体实施方式中分布动态载荷等效识别方法的步骤流程，具体实施步骤如下：步骤1：将结构所受的未知分布动态载荷进行本征正交分解，转化为一系列时间与空间独立的子载荷；考虑如图2所示的时间与空间耦合的分布动态载荷F(x,t)＝1000sin(20πxt)，0≤x≤1m，0≤t≤0.1s。通过本征正交分解能将其转化为如式(1)所示的一系列时间与空间独立的子载荷，式中，F(x,t)表示作用于结构上的时间与空间耦合的分布动态载荷，通过本征正交分解可将其近似为L项时间和空间独立的子载荷之和，ψn(x)表示第n个子载荷的空间分布函数，sn(t)表示第n个子载荷的时间历程。在该实施例中，待识别的分布动态载荷可等效为转化如图3所示的4个时间与空间独立的子载荷之和。步骤2：建立结构有限元模型，并进行模态分析，获得结构的模态振型系数矩阵、模态质量、模态阻尼、模态刚度、无阻尼固有频率、阻尼固有频率及阻尼比；如图4所示，该实施例中采用了受分布动态载荷作用的某汽车引擎盖内板模型，共包含10775个壳单元，10685个节点，均匀厚度为0.8mm，四个角被固定约束，材料的弹性模量为E＝210GPa,泊松比为μ＝0.3，密度为ρ＝7.65×103kg/m3。分布动态载荷作用在其前端边缘的中间区域。结构的模态振型系数矩阵、模态质量、模态阻尼、模态刚度、无阻尼固有频率、阻尼固有频率及阻尼比通过有限元方法分析得到。步骤3：利用传感器测量分布动态载荷作本文档来自技高网...

【技术保护点】
一种空间分布动态载荷的等效识别方法，其特征在于其特征在于，包括如下步骤：步骤1：将结构所受的未知分布动态载荷进行本征正交分解，转化为一系列时间与空间独立的子载荷；步骤2：建立结构有限元模型，并进行模态分析，获得结构的模态振型系数矩阵、模态质量、模态阻尼、模态刚度、无阻尼固有频率、阻尼固有频率及阻尼比；步骤3：利用传感器测量分布动态载荷作用下结构的动态响应，并进行模态变换获取相应的截断模态响应；步骤4：基于模态响应和模态格林核函数，利用正则化方法对有限阶次的模态载荷进行准确稳定识别；步骤5：将模态载荷等效为有限个集中动态载荷的混叠，利用盲源分离技术从识别的模态载荷中获取等效集中载荷的时间历程，并计算相应混叠系数；步骤6：考虑混叠系数为载荷空间分布与结构模态振型的内积，对混叠系数向量进行稀疏分解，识别等效集中载荷的空间作用位置。

【技术特征摘要】
1.一种空间分布动态载荷的等效识别方法，其特征在于其特征在于，包括如下步骤：步骤1：将结构所受的未知分布动态载荷进行本征正交分解，转化为一系列时间与空间独立的子载荷；步骤2：建立结构有限元模型，并进行模态分析，获得结构的模态振型系数矩阵、模态质量、模态阻尼、模态刚度、无阻尼固有频率、阻尼固有频率及阻尼比；步骤3：利用传感器测量分布动态载荷作用下结构的动态响应，并进行模态变换获取相应的截断模态响应；步骤4：基于模态响应和模态格林核函数，利用正则化方法对有限阶次的模态载荷进行准确稳定识别；步骤5：将模态载荷等效为有限个集中动态载荷的混叠，利用盲源分离技术从识别的模态载荷中获取等效集中载荷的时间历程，并计算相应混叠系数；步骤6：考虑混叠系数为载荷空间分布与结构模态振型的内积，对混叠系数向量进行稀疏分解，识别等效集中载荷的空间作用位置。2.根据权利要求1所述的方法，其特征在于，在所述步骤1中，待识别的空间分布动态载荷形式可以是时间和空间相互独立的，也可以时间和空间相互耦合的；当时间和空间相互耦合时，通过本征正交分解能将其转化为如式(1)所示的一系列时间与空间独立的子载荷：式中，F(x,t)表示作用于结构上的时间与空间耦合的分布动态载荷，通过本...

【专利技术属性】
技术研发人员：刘杰，李锟，姜潮，韩旭，
申请(专利权)人：湖南大学，
类型：发明
国别省市：湖南,43