一种实验模态频率的丢失辨识及预估方法技术

本发明专利技术涉及一种实验模态频率的丢失辨识及预估方法，属于机械振动技术领域。本发明专利技术的目的是提供一种实验模态频率的丢失辨识及预估方法，依据计算模态频率与实验模态频率具有强的正相关性，将实验模态频率与相应的计算模态频率进行相关性分析，进而辨识实验模态频率丢失的阶数，同时预估相应丢失的实验模态频率。该方法克服了多点激励多点拾取中设备成本高，调试效率低的问题，也克服了计算模态振型与实验模态振型振型相关性分析中的计算量大，计算时间长，辨识效率低的问题。

【技术实现步骤摘要】

本专利技术涉及一种实验模态频率的丢失辨识及预估方法，属于机械振动

技术介绍
通过实验测量获得对象的实验模态频率的过程中，如果激励点或拾振点位于某一阶模态振型的节点处，亦或是实验测量平台系统的噪音对某一阶模态频率具有明显的干扰效应，都会引起实验模态频率阶数丢失的情况发生。目前，实验模态频率阶数丢失的解决方案主要有两种，一种是在对象上采用多点激励并多点拾取的措施，该方法能够避免实验模态频率阶数的丢失，但该方法适用于大中型结构，且需要较大的激励能量的情况，同时该方法中的各激励器需反复调试，致使调试效率低，设备成本高，另一种是将计算模态振型与实验模态振型进行振型相关性分析，该方法能够辨识实验模态频率丢失的阶数，但该方法中的计算量大，计算时间长，致使该方法对实验模态频率丢失阶数的辨识效率低。
技术实现思路
本专利技术正是为了克服上述现有技术的不足而提供的一种实验模态频率的丢失辨识及预估方法本专利技术的目的是提供一种实验模态频率的丢失辨识及预估方法，依据计算模态频率与实验模态频率具有强的正相关性，将实验模态频率与相应的计算模态频率进行相关性分析，进而辨识实验模态频率丢失的阶数，同时预估相应丢失的实验模态频率。该方法克服了多点激励多点拾取中设备成本高，调试效率低的问题，也克服了计算模态振型与实验模态振型振型相关性分析中的计算量大，计算时间长，辨识效率低的问题。本专利技术是通过如下技术方案来实现的：一种实验模态频率的丢失辨识及预估方法，包括如下步骤：步骤1在有限元软件中，导入加载对象的有限元模型；步骤2在有限元软件中，设置有限元模型的杨氏模量Ey；步骤3在有限元软件中，计算对象的各阶模态频率w1y，w2y……wny；步骤4通过实验平台，测出对象的各阶模态频率w1c，w2c……wnc；步骤5辨识丢失的实验模态频率阶数；步骤6预估丢失的实验模态频率。其中辨识丢失的实验模态频率阶数的具体步骤包括：5.1计算k，其中k＝w1c/w1y；5.2令w1＝w1c，wj为对象的各阶模态频率，j为对象各阶模态频率的下标号，wic为测出对象的各阶模态频率，i为测出对象各阶模态频率的下标号，wpy为计算对象的各阶模态频率，p为计算对象各阶模态频率的下标号，其中j，i，p均为整数，2≤j，2≤i≤n，2≤p≤z，令j＝2，i＝2，p＝2；5.3计算m，Δm，b，Δb，其中，5.4判断Δm，Δb的大小关系若Δm≤Δb，则未丢失该阶实验模态频率，即wj＝wic，判断i是否等于n，若i≠n，则i加1，j加1，p加1，否则输出对象的各阶模态频率wj，其中n为通过实验平台测出对象的模态频率最大阶数，同时舍弃剩余未使用的计算对象的各阶模态频率；若Δm>Δb，则为丢失该阶实验模态频率，即wj＝0，判断i是否等于n，若i≠n，则j加1，p加1，否则输出对象的各阶模态频率wj，其中n为通过实验平台测出对象的模态频率最大阶数，同时舍弃剩余未使用的计算对象的各阶模态频率；5.5重复步骤5.3-5.4。其中预估丢失的实验模态频率的具体步骤包括：6.1依据输出的wj＝0，读取丢失的相应实验模态频率wt，wt+1，wt+2……wt+r以及相邻的实验模态频率wt-1，wt+r+1，其中，t，t+1，t+2……t+r为丢失的实验模态频率的下标，t-1，t+r+1为相邻的实验模态频率的下标；6.2对丢失的实验模态频率进行赋值，即wt＝uw(t)y，wt+1＝uw(t+1)y，wt+2＝uw(t+2)y……wt+r＝uw(t+r)y，其中u为待求系数，等式右端的t，t+1，t+2……t+r为相应的计算对象各阶模态频率的下标；6.3构造关于待求系数u的函数f(u)，即并求解f(u)最小值时的u；6.4将所求的系数u代入到步骤302中的等式中，即可预估丢失的实验模态频率wt，wt+1，wt+2……wt+r。本专利技术具有以下有益效果：1、能够减少实验设备的调试量，提升实验设备的调试效率；2、对实验设备要求低，能够降低实验测量成本；3、能够减少计算量和计算时间，提高实验模态频率丢失阶数的辨识效率；4、能够预估相应丢失的实验模态频率。附图说明图1为实验模态频率的丢失辨识及预估方法示意图；图2为辨识丢失的实验模态频率阶数流程图；图3为预估丢失的实验模态频率流程图。具体实施方式下面结合附图和实施例，对本专利技术作进一步说明，但本专利技术的内容并不限于所述范围。实施例1：如图1所示，一种实验模态频率的丢失辨识及预估方法，包括如下具体步骤：步骤1在有限元软件中，导入加载对象的有限元模型；步骤2在有限元软件中，设置有限元模型的杨氏模量Ey；步骤3在有限元软件中，计算对象的各阶模态频率w1y，w2y…，wpy，…wzy，其中z为通过有限元软件计算出对象的模态频率最大阶数；步骤4通过实验平台，测出对象的各阶模态频率w1c，w2c…，wic，…wnc，其中n为通过实验平台测出对象的模态频率最大阶数，且n≤z；步骤5辨识丢失的实验模态频率阶数如图2所示，5.1计算k，其中k＝w1c/w1y；5.2令w1＝w1c，wj为对象的各阶模态频率，j为对象各阶模态频率的下标号，wic为测出对象的各阶模态频率，i为测出对象各阶模态频率的下标号，wpy为计算对象的各阶模态频率，p为计算对象各阶模态频率的下标号，其中j，i，p均为整数，2≤j，2≤i≤n，2≤p≤z，令j＝2，i＝2，p＝2；5.3计算m，Δm，b，Δb，其中，5.4判断Δm，Δb的大小关系若Δm≤Δb，则未丢失该阶实验模态频率，即wj＝wic，判断i是否等于n，若i≠n，则i加1，j加1，p加1，否则输出对象的各阶模态频率wj，其中n为通过实验平台测出对象的模态频率最大阶数，同时舍弃剩余未使用的计算对象的各阶模态频率；若Δm>Δb，则为丢失该阶实验模态频率，即wj＝0，判断i是否等于n，若i≠n，则j加1，p加1，否则输出对象的各阶模态频率wj，其中n为通过实验平台测出对象的模态频率最大阶数，同时舍弃剩余未使用的计算对象的各阶模态频率；5.5重复步骤5.3-5.4。步骤6预估丢失的实验模态频率如图3所示，6.1依据输出的wj＝0，读取丢失的相应实验模态频率wt，wt+1，wt+2……wt+r以及相邻的实验模态频率wt-1，wt+r+1，其中，t，t+1，t+2……t+r为丢失的实验模态频率的下标，t-1，t+r+1为相邻的实验模态频率的下标；6.2对丢失的实验模态频率进行赋值，即wt＝uw(t)y，wt+1＝uw(t+1)y，wt+2＝uw(t+2)y……wt+r＝uw(t+r)y，其中u为待求系数，等式右端的t，t+1，t+2……t+r为相应的计算对象各阶模态频率的下标；6.3构造关于待求系数u的函数f(u)，即并求解f(u)最小值时的u；6.4将所求的系数u代入到步骤302中的等式中，即可预估丢失的实验模态频率wt，wt+1，wt+2……wt+r。实施例2：步骤1在有限元软件中，导入加载对象的有限元模型。步骤2在有限元软件中，设置有限元模型的杨氏模量Ey＝1010Pa。步骤3在有限元软件中，计算对象的各阶模态频率w1y＝10.2Hz，w2y＝15.4Hz，w3y＝22.6Hz，w4y＝32.1Hz，w5y＝53.6Hz，w6y＝66.本文档来自技高网...

【技术保护点】
一种实验模态频率的丢失辨识及预估方法，其特征在于：包括如下具体步骤：步骤1在有限元软件中，导入加载对象的有限元模型；步骤2在有限元软件中，设置有限元模型的杨氏模量Ey；步骤3在有限元软件中，计算对象的各阶模态频率w1y，w2y…，wpy，…wzy，其中z为通过有限元软件计算出对象的模态频率最大阶数；步骤4通过实验平台，测出对象的各阶模态频率w1c，w2c…，wic，…wnc，其中n为通过实验平台测出对象的模态频率最大阶数，且n≤z；步骤5辨识丢失的实验模态频率阶数5.1计算k，其中k＝w1c/w1y；5.2令w1＝w1c，wj为对象的各阶模态频率，j为对象各阶模态频率的下标号，wic为测出对象的各阶模态频率，i为测出对象各阶模态频率的下标号，wpy为计算对象的各阶模态频率，p为计算对象各阶模态频率的下标号，其中j，i，p均为整数，2≤j，2≤i≤n，2≤p≤z，令j＝2，i＝2，p＝2；5.3计算m，Δm，b，Δb，其中，5.4判断Δm，Δb的大小关系若Δm≤Δb，则未丢失该阶实验模态频率，即wj＝wic，判断i是否等于n，若i≠n，则i加1，j加1，p加1，否则输出对象的各阶模态频率wj，其中n为通过实验平台测出对象的模态频率最大阶数，同时舍弃剩余未使用的计算对象的各阶模态频率；若Δm>Δb，则为丢失该阶实验模态频率，即wj＝0，判断i是否等于n，若i≠n，则j加1，p加1，否则输出对象的各阶模态频率wj，其中n为通过实验平台测出对象的模态频率最大阶数，同时舍弃剩余未使用的计算对象的各阶模态频率；5.5重复步骤5.3‑5.4；步骤6预估丢失的实验模态频率6.1依据输出的wj＝0，读取丢失的相应实验模态频率wt，wt+1，wt+2……wt+r以及相邻的实验模态频率wt‑1，wt+r+1，其中，t，t+1，t+2……t+r为丢失的实验模态频率的下标，t‑1，t+r+1为相邻的实验模态频率的下标；6.2对丢失的实验模态频率进行赋值，即wt＝uw(t)y，wt+1＝uw(t+1)y，wt+2＝uw(t+2)y……wt+r＝uw(t+r)y，其中u为待求系数，等式右端的t，t+1，t+2……t+r为相应的计算对象各阶模态频率的下标；6.3构造关于待求系数u的函数f(u)，即并求解f(u)最小值时的u；6.4将所求的系数u代入到步骤302中的等式中，即可预估丢失的实验模态频率wt，wt+1，wt+2……wt+r。...

【技术特征摘要】
1.一种实验模态频率的丢失辨识及预估方法，其特征在于：包括如下具体步骤：步骤1在有限元软件中，导入加载对象的有限元模型；步骤2在有限元软件中，设置有限元模型的杨氏模量Ey；步骤3在有限元软件中，计算对象的各阶模态频率w1y，w2y…，wpy，…wzy，其中z为通过有限元软件计算出对象的模态频率最大阶数；步骤4通过实验平台，测出对象的各阶模态频率w1c，w2c…，wic，…wnc，其中n为通过实验平台测出对象的模态频率最大阶数，且n≤z；步骤5辨识丢失的实验模态频率阶数5.1计算k，其中k＝w1c/w1y；5.2令w1＝w1c，wj为对象的各阶模态频率，j为对象各阶模态频率的下标号，wic为测出对象的各阶模态频率，i为测出对象各阶模态频率的下标号，wpy为计算对象的各阶模态频率，p为计算对象各阶模态频率的下标号，其中j，i，p均为整数，2≤j，2≤i≤n，2≤p≤z，令j＝2，i＝2，p＝2；5.3计算m，Δm，b，Δb，其中，5.4判断Δm，Δb的大小关系若Δm≤Δb，则未丢失该阶实验模态频率，即wj＝wic，判断i是否等于n，若i≠n，则i加1，j加1，p加1，否则输出对象的各阶模态频率wj，其中n为通过...

【专利技术属性】
技术研发人员：刘孝保，刘洋，刘璐，
申请(专利权)人：昆明理工大学，
类型：发明
国别省市：云南;53