[发明专利]一种基于小波多分辨率分析的动载荷识别方法

说明书

技术领域

本发明涉及载荷识别领域，特别涉及一种基于小波多分辨率分析的动载荷识别方法。

背景技术

工程结构上的动载荷是结构的动力学设计、减振隔振设计、强度校核、动力学优化等的重要原始参数，然而由于技术或经济条件的限制，很多情况下结构所受的动态载荷很难直接进行测量。载荷识别就是在这样的背景下产生的，其可以在许多无法直接量测载荷的情况下通过响应等信息反求载荷。

传统载荷识别方法按信号所处域可分为频域识别和时域识别。载荷识别是一类非常基础的问题，研究中还存在一定问题，如频域下的载荷识别要求测量数据的样本具有一定的长度，一般只适用于稳态或平稳随机载荷的识别，对冲击、突变等瞬态载荷的识别不理想，当频率接近共振区时会出现数值不稳定问题；时域下的载荷识别包含递推，对边值条件和初始条件敏感，初始的误差会随着递推推导不断传递造成误差累积，有些方法需要知道一种类型以上的响应，给使用带来困难；小波分析是一种新颖的时频分析方法，其通过伸缩平移运算对信号逐步进行多尺度细化，最终达到高频处时间细分，低频处频率细分，能自动适应时频信号分析的要求，从而可聚焦到信号的任意细节。因此，一些研究将小波分析用于载荷识别来获得更理想的识别结果。基于小波分析的载荷识别方法是在小波的时频域下识别载荷，利用小波如尺度特性、平移特性、紧支性、移位正交性、多分辨率分析等优良性质，可识别包括冲击载荷在内的任意确定性载荷，但由于发展时间短，方法的理论基础并不完善，部分推导的表达式针对特定问题，局限性大，一些参数没有定性或定量的选择依据，一些方法没有利用二维分析特性。

发明内容

本发明的目的在于，提出一种基于小波多分辨率分析的动载荷识别方法，用以解决目前载荷识别研究中存在的问题。

为实现上述目的，本发明提出的技术方案是，一种基于小波多分辨率分析的动载荷识别方法，其特征是所述方法包括下列步骤：

步骤1：求解识别参数，积分上下限J和I、分解层级j；

步骤2：基于时域卷积模型，利用小波基函数对载荷进行小波重构，得到小波响应函数其中，t表示时间，h(t)是单位脉冲响应函数，φ(2^jt-k)是小波基函数，k为平移因子，j为分解层级(尺度因子)；

步骤3：对响应和小波响应函数进行小波变换，得到小波域下的系统响应和小波域下的小波响应函数其中y(t)是系统的响应；

步骤4：反求载荷在小波域下的权系数P＝G^-1u，其中u和G分别为小波域下的系统响应u(i)和小波域下的小波响应函数g_ik的矩阵形式；

步骤5：由权系数计算载荷其中P(k)是权系数矩阵P的每一个元素，并完成识别。

本发明实现了一种基于小波多分辨率分析的动载荷识别方法，其有益效果如下：

1给出一种基于小波多分辨率分析的动载荷识别方法，方法可识别非平稳载荷，对冲击、突变等快变瞬态载荷识别能力优秀；

2给出的方法对多输入多输出系统中同时包含快慢变化的多路载荷间的相互干扰不敏感，可以区别并辨识多路载荷；

3提出了多个识别参数的定量/定性确定方法，可用于基于小波分析的动载荷识别方法的参数确定。

附图说明

图1本发明方法的流程图。

图2是实验示意图。

图3是悬臂梁节点布置示意图。

图4是原载荷和识别载荷的时域图。

具体实施方式

下面结合附图，对优选实施例作详细说明。应该强调的是，下述说明仅仅是示例性的，而不是为了限制本发明的范围及其应用。

本发明解决问题的思路是基于小波和其多分辨率分析性质完成动载荷识别：首先求解识别参数；其次，计算小波响应函数；再次，计算小波域下的系统响应和小波域下的小波响应函数；最后，计算权系数，反求载荷，完成识别。

图1是本发明完成的一种基于小波多分辨率分析的动载荷识别方法的流程图，该方法具体步骤如下：

步骤1：量测频响函数(Frequency Response Functions，FRF)。

实验示意图如图2所示，首先使用力锤激励法量测悬臂梁系统的频响函数信息，悬臂梁节点布置如图3所示，节点2、3布置传感器，节点1、4不放置激振器，使用力锤敲击节点1，敲击5次取平均值，得到FRF21和FRF31，同理可得FRF24和FRF34，其中数字如21表示第1个节点激励下第2个节点的频响函数。

步骤2：计算单位脉冲响应函数(Impulse Response Functions，IRF)。