首页 / 专利库 / 电脑编程 / 算法 / 一种基于分组稀疏性的动态载荷识别方法

一种基于分组稀疏性的动态载荷识别方法

阅读:551发布:2021-06-08

专利汇可以提供一种基于分组稀疏性的动态载荷识别方法专利检索,专利查询,专利分析的服务。并且本 发明 提供一种基于分组稀疏性的动态 载荷 识别方法。本发明的方法首先将动载荷识别的三维问题转化为二维问题,然后选择合适的惩罚函数,并对结构上的点按 位置 进行分组。最后使用G-FISTA方法优化求解,得到最优的动态载荷向量。本发明所述的方法,对于不便于直接安装 力 传感器 测量动载荷的情形,是非常有帮助的。只需测量结构上几个点的振动 加速 度,在已知结构频响函数的前提下,即可间接获得结构上的动载荷。,下面是一种基于分组稀疏性的动态载荷识别方法专利的具体信息内容。

1.一种基于分组稀疏性的动态载荷识别方法,其特征在于,该方法步骤如下:
步骤(1).将动态载荷识别由三维问题转化为二维问题;
由于矩阵必须以二维形式进行计算,将三维矩阵转换成二维形式,分别将X,H,F转换成步骤(2).选择合适的惩罚函数g;
步骤(3).对结构上的点按位置进行分组;
对同一组内的元素,采用连续性限制与l2范数惩罚;对不同组的元素,采用稀疏性限制与l1范数惩罚;
步骤(4).采用优化方法求解;
采用G-FISTA算法对损失函数进行优化,解决lp正则化问题;最后,使用贝叶斯信息准则方法,从所有的解 中找到最优的解
2.根据权利要求1所述的一种基于分组稀疏性的动态载荷识别方法,其特征在于,步骤(1)具体方法如下:
要求解的初始问题为式(1),在感兴趣的频段内,一共有Nf个频谱线,分别将Ff,连接成一个列向量得到新的矩阵 将所有的Hf矩阵放置在对
线上,形成一个矩阵 Hf是对应于频率f的频响函数矩阵, 中包含大量的0元素;重排列后的二维结构如式(2)所示;
X=HF               (1)
其中,X是n×1×Nf的加速度向量,n是加速度传感器的个数;
F是k×1×Nf的向量,k是力作用位置的个数;H(f)是n×k×Nf的频响函数(FRF)模型;
通过上述处理,三维问题转化为二维问题,公式如下,
通过最小二乘估计方法求解二维问题,公式如下;
由于 的行数少于列数,上式所描述的问题是病态的,所以采用惩罚最小二乘求解,公式如下;
式中, 是惩罚函数,λ是正则化参数。
3.根据权利要求2所述的一种基于分组稀疏性的动态载荷识别方法,其特征在于,步骤(2)具体方法如下:
在病态问题求解中,额外的信息可以补偿已知方程个数的缺少;根据二维问题的物理特点,选择的惩罚函数需要满足以下限制:
1)稀疏性限制:由于外部力通常只作用于结构上较少的位置,因此力向量在大部分点处为0,待估计力对应于各个位置是稀疏解;因此,采用l1范数惩罚;
2)连续性限制:当结构受到局部的外部力激励(锤子或激振器),力在频率轴上应该是连续的;在一个脉冲激励的情形下,力频谱是宽频段的;为保证力向量在频率轴上的连续性,采用l2范数惩罚;在结构上某个位置,力在频率轴的取值是一致的,对每个载荷作用位置,估计出的力在整个频率轴上是全部为0或全部非0;
3)损失函数是凸函数:为了确保对应一个λ,只有一个优化解,损失函数应该是凸函数;
当lp的两项都是凸函数,综合的损失函数也是凸函数;当p≥1时,lp范数是凸函数。
4.根据权利要求3所述的一种基于分组稀疏性的动态载荷识别方法,其特征在于,步骤(3)具体方法如下:
采用满足稀疏性、连续性和凸函数限制的惩罚函数;将不同频率的力按照作用点进行分组,分组被定义成如下形式:
Gi={Fi(f1),...,Fi(fNf)},i∈{1,...,k}           (6)
对同一组内的元素,使用连续性限制,使用l2范数惩罚;对不同组的元素,使用稀疏性限制,使用l1范数惩罚;如式(7)所示,新的惩罚函数是一个综合l2与l1范数惩罚;新的惩罚最小二乘问题定义成式(8),式(8)是式(5)的特例;
5.根据权利要求4所述的一种基于分组稀疏性的动态载荷识别方法,其特征在于,步骤(4)具体方法如下:
问题:对于一个给定的λ,
-2 8
输入 λ∈{λ1,...,λr}(预定义的集合,且10 ≤λ≤10);
输出
所述的G-FISTA方法具体如下:
①初始化,
Gk;
的最大特征值;
α(0)=1;
i=0;
②i=i+1;
Trλ=λ/L;
检查收敛准则:是,转到③;否,转到②;
③模型选择:
式中,nr是残差向量Rλ中零元素的个数; q是未知Nf×k的总数;
G-FISTA方法基于频响函数 和测量的加速度向量 估计力向量
对于每个λ,算法会分别估计出一个力向量;算法的迭代部分为最小化优化(Majorize–(i)
Minimization)算法的一个特例;参数 和α 是中间变量,用于算法收敛;在每次迭代中,通过阈值Rλ,进行算法收敛,估计力的值;算法中的pL(u)函数是投影算子;当获得的解 相对于上一步的获得的解 的相对误差小于tol时,算法终止,公式如下;
通过贝叶斯信息准则方法,从所有的解 中,找到最优的解

说明书全文

一种基于分组稀疏性的动态载荷识别方法

技术领域

[0001] 本发明涉及结构振动技术领域,具体涉及一种基于分组稀疏性的动态载荷识别方法。

背景技术

[0002] 当机械结构受到动态载荷作用时,会产生振动现象。机械振动的产生会破坏系统的安全运行状态,并使整个系统的性能变差。对动态载荷大小和位置的获取是十分必要的。一方面,获取动态载荷后,可以用来指导结构的设计。对结构进行改进以此降低动态载荷。
另一方面,动态载荷对结构的强度和寿命有重要影响,可以利用动态载荷更加准确的计算其可靠性信息。
[0003] 为了获得结构所受动态载荷的大小,传统的方法是在指定位置上安装传感器。但是有的结构在工作状态下,不便于安装力传感器。因此需要设计一种新的方法,通过测量结构上某个部位的振动位移大小,间接的求出结构受到的动态载荷。
[0004] 已有的动态载荷求解方法中,只能求解载荷作用点个数少于加速度测试点的情况,即超静定方程的求解。而对于载荷作用点个数少于加速度测试点的情况,即病态方程的求解,还没有人提出过解决方法。

发明内容

[0005] 针对在结构受到动态载荷作用,不能直接测量出动态载荷大小时,如何利用振动加速度和结构频响函数求出动态载荷的大小的问题,本发明提供了一种基于分组稀疏性的动态载荷识别方法。
[0006] 本发明的方法是先将动载荷识别的三维问题转化为二维问题,然后选择合适的惩罚函数,并对结构上的点按位置进行分组。最后使用G-FISTA方法优化求解,得到最优的动态载荷向量
[0007] 一种基于分组稀疏性的动态载荷识别方法,步骤如下:
[0008] 步骤(1).将动态载荷识别由三维问题转化为二维问题。
[0009] 由于矩阵必须以二维形式进行计算,将三维矩阵转换成二维形式,分别将X,H,F转换成
[0010] 步骤(2).选择合适的惩罚函数g。
[0011] 步骤(3).对结构上的点按位置进行分组;
[0012] 对同一组内的元素,采用连续性限制与l2范数惩罚。对不同组的元素,采用稀疏性限制与l1范数惩罚。
[0013] 步骤(4).采用优化方法求解;
[0014] 采用G-FISTA算法对损失函数进行优化,解决lp正则化问题。最后,使用贝叶斯信息准则方法,从所有的解 中找到最优的解
[0015] 步骤(1)具体方法如下:
[0016] 要求解的初始问题为式(1),在感兴趣的频段内,一共有Nf个频谱线,分别将Ff,Xf连接成一个列向量得到新的矩阵 将所有的Hf矩阵放置在对线上,形成一个矩阵 Hf是对应于频率f的频响函数矩阵, 中包含大量的0元素。重排列后的二维结构如式(2)所示。
[0017] X=HF  (1)
[0018] 其中,X是n×1×Nf的加速度向量,n是加速度传感器的个数。
[0019] F是k×1×Nf的力向量,k是力作用位置的个数。H(f)是n×k×Nf的频响函数(FRF)模型。
[0020]
[0021] 通过上述处理,三维问题转化为二维问题,公式如下,
[0022]
[0023] 通过最小二乘估计方法求解二维问题,公式如下。
[0024]
[0025] 由于 的行数少于列数,上式所描述的问题是病态的,所以采用惩罚最小二乘求解,公式如下。
[0026]
[0027] 式中, 是惩罚函数,λ是正则化参数。
[0028] 步骤(2)具体方法如下:
[0029] 在病态问题求解中,额外的信息可以补偿已知方程个数的缺少。根据二维问题的物理特点,选择的惩罚函数需要满足以下限制:
[0030] 1)稀疏性限制:由于外部力通常只作用于结构上较少的位置,因此力向量在大部分点处为0,待估计力对应于各个位置是稀疏解。因此,采用l1范数惩罚。
[0031] 2)连续性限制:当结构受到局部的外部力激励(锤子或激振器),力在频率轴上应该是连续的。在一个脉冲激励的情形下,力频谱是宽频段的。为保证力向量在频率轴上的连续性,采用l2范数惩罚。在结构上某个位置,力在频率轴的取值是一致的,对每个载荷作用位置,估计出的力在整个频率轴上是全部为0或全部非0。
[0032] 3)损失函数是凸函数:为了确保对应一个λ,只有一个优化解,损失函数应该是凸函数。当lp的两项都是凸函数,综合的损失函数也是凸函数。当p≥1时,lp范数是凸函数。
[0033] 步骤(3)具体方法如下:
[0034] 采用满足稀疏性、连续性和凸函数限制的惩罚函数。将不同频率的力按照作用点进行分组,分组被定义成如下形式:
[0035] Gi={Fi(f1),...,Fi(fNf)},i∈{1,...,k}  (6)
[0036] 对同一组内的元素,使用连续性限制,使用l2范数惩罚。对不同组的元素,使用稀疏性限制,使用l1范数惩罚。如式(7)所示,新的惩罚函数是一个综合l2与l1范数惩罚。新的惩罚最小二乘问题定义成式(8),式(8)是式(5)的特例。
[0037]
[0038]
[0039] 步骤(4)具体方法如下:
[0040] 问题:对于一个给定的λ,
[0041]
[0042] 输入 λ∈{λ1,...,λr}(预定义的集合,且10-2≤λ≤108);
[0043] 输出
[0044] 所述的G-FISTA方法具体如下:
[0045] ①初始化,
[0046]
[0047] Gk;
[0048] 的最大特征值;
[0049] α(0)=1;
[0050]
[0051] i=0;
[0052] ②i=i+1;
[0053]
[0054]
[0055] Trλ=λ/L;
[0056]
[0057]
[0058]
[0059] 检查收敛准则:是,转到③;否,转到②。
[0060] ③模型选择:
[0061]
[0062] 式中,nr是残差向量Rλ中零元素的个数。 q是未知Nf×k的总数。
[0063]
[0064] G-FISTA方法基于频响函数 和测量的加速度向量 估计力向量
[0065] 对于每个λ,算法会分别估计出一个力向量。算法的迭代部分为最小化优化(Majorize–Minimization)算法的一个特例。参数 和α(i)是中间变量,用于算法收敛。在每次迭代中,通过阈值Rλ,进行算法收敛,估计力的值。算法中的pL(u)函数是投影算子。当获得的解 相对于上一步的获得的解 的相对误差小于tol时,算法终止,公式如下;
[0066]
[0067] 通过贝叶斯信息准则方法,从所有的解 中,找到最优的解
[0068] 本发明的有益效果如下:
[0069] 本发明所述的方法,对于不便于直接安装力传感器测量动载荷的情形,是非常有帮助的。只需测量结构上几个点的振动加速度,在已知结构频响函数的前提下,即可间接获得结构上的动载荷。附图说明
[0070] 图1为本发明的动载荷识别问题的三维结构;
[0071] 图2为本发明的G-FISTA方法的流程图

具体实施方式

[0072] 以下结合附图,对本发明进行进一步的详细说明。
[0073] 本发明的动态载荷识别方法,步骤如下:
[0074] 步骤(1).将动态载荷识别由三维问题转化为二维问题。
[0075] 由于矩阵必须以二维形式进行计算,将三维矩阵转换成二维形式,分别将X,H,F转换成
[0076] 要求解的初始问题为式(1),在感兴趣的频段内,一共有Nf个频谱线,由于不能直接计算三维矩阵 ,因此该三维结构需要被重排列为二维结 构。分 别将连接成一个列向量得到新的矩阵 将所有的Hf矩阵放置在
f
对角线上,形成一个块矩阵 H是对应于频率f的频响函数矩阵, 中包含大量的0元素。
重排列后的二维结构如式(2)所示。图1为本发明的动载荷识别问题的三维结构。
[0077] X=HF  (1)
[0078] 式中,X是n×1×Nf的加速度向量,n是加速度传感器的个数。
[0079] F是k×1×Nf的力向量,k是力作用位置的个数。H(f)是n×k×Nf的频响函数(FRF)模型。
[0080]
[0081] 通过上述处理,三维问题转化为二维问题,公式如下,
[0082]
[0083] 通过最小二乘估计方法求解二维问题,公式如下。
[0084]
[0085] 由于 的行数少于列数,上式所描述的问题是病态的,所以采用惩罚最小二乘求解,公式如下。
[0086]
[0087] 式中, 是惩罚函数,λ是正则化参数。
[0088] 步骤(2).选择合适的惩罚函数g。
[0089] 在病态问题求解中,额外的信息可以补偿已知方程个数的缺少。根据二维问题的物理特点,选择的惩罚函数需要满足以下限制:
[0090] 1)稀疏性限制:由于外部力通常只作用于结构上较少的位置,因此力向量在大部分点处为0,待估计力对应于各个位置是稀疏解。因此,采用l1范数惩罚。
[0091] 2)连续性限制:当结构受到局部的外部力激励(锤子或激振器),力在频率轴上应该是连续的。在一个脉冲激励的情形下,力频谱是宽频段的。为保证力向量在频率轴上的连续性,采用l2范数惩罚。在结构上某个位置,力在频率轴的取值是一致的,对每个载荷作用位置,估计出的力在整个频率轴上是全部为0或全部非0。
[0092] 3)损失函数是凸函数:为了确保对应一个λ,只有一个优化解,损失函数应该是凸函数。当lp的两项都是凸函数,综合的损失函数也是凸函数。当p≥1时,lp范数是凸函数。
[0093] 步骤(3).对结构上的点按位置进行分组;
[0094] 对同一组内的元素,采用连续性限制与l2范数惩罚。对不同组的元素,采用稀疏性限制与l1范数惩罚。
[0095] 采用满足稀疏性、连续性和凸函数限制的惩罚函数。将不同频率的力按照作用点进行分组,分组被定义成如下形式:
[0096] Gi={Fi(f1),...,Fi(fNf)},i∈{1,...,k}  (6)
[0097] 对同一组内的元素,使用连续性限制,使用l2范数惩罚。对不同组的元素,使用稀疏性限制,使用l1范数惩罚。如式(7)所示,新的惩罚函数是一个综合l2与l1范数惩罚。新的惩罚最小二乘问题定义成式(8),式(8)是式(5)的特例。
[0098]
[0099]
[0100] 步骤(4).采用优化方法求解;
[0101] 采用G-FISTA算法对损失函数进行优化,解决lp正则化问题。最后,使用贝叶斯信息准则方法,从所有的解 中找到最优的解
[0102] 式(5)所示的正则化问题一般采用任意稀疏优化方法求解。因为式(8)中的惩罚函数比较特殊,在选择优化方法时有限制。式(8)中描述的问题包含光滑部分(最小二乘)和非光滑部分(综合l2与l1范数惩罚)。对这种损失函数进行优化的方法中,Lasso和FISTA是最常用的。Lasso只能解决l1正则化问题,FISTA可以解决lp正则化问题。但原始的FISTA方法,是对于实参数解决惩罚最小二乘问题,并且变量之间没有分组,因此需要对其进行两个方面的修改,以适应式(8)的求解。修改后算法可以针对分组变量求解,并能解决l1&2惩罚最小二乘。
[0103] 问题:对于一个给定的λ,
[0104]
[0105] 输入 λ∈{λ1,...,λr}(预定义的集合,且10-2≤λ≤108);
[0106] 输出
[0107] 如图2所示,G-FISTA方法具体如下:
[0108] ①初始化,
[0109]
[0110] Gk;
[0111] 的最大特征值;
[0112] α(0)=1;
[0113]
[0114] i=0;
[0115] ②i=i+1;
[0116]
[0117]
[0118] Trλ=λ/L;
[0119]
[0120]
[0121]
[0122] 检查收敛准则:是,转到③;否,转到②。
[0123] ③模型选择:
[0124]
[0125] 式中,nr是残差向量Rλ中零元素的个数。 q是未知Nf×k的总数。
[0126]
[0127] G-FISTA方法基于频响函数 和测量的加速度向量 估计力向量
[0128] 对于每个λ,算法会分别估计出一个力向量,算法的迭代部分为最小化优化(Majorize–Minimization)算法的一个特例。参数 和α(i)是中间变量,用于算法收敛。在每次迭代中,通过阈值Rλ,进行算法收敛,估计力的值。算法中的pL(u)函数是投影算子。当获得的解 相对于上一步的获得的解 的相对误差小于tol时,算法终止,公式如下;
[0129]
[0130] 通过贝叶斯信息准则方法,从所有的解 中,找到最优的解 本发明提出的一种求解动载荷的方法,对于不便于直接安装力传感器测量动载荷的情形,可以使用该方法间接获得动载荷。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈