质量时变系统动载荷识别-修改稿.docVIP

质量时变系统的动载荷识别1） 彭 凡*，王 樑* , 肖 健+ *（湖南大学机械与运载工程学院，长沙 100190） +（北京强度与环境研究所，北京 100076) 摘要：研究质量时变系统的动态载荷识别技术。引入广义模态概念，在初瞬时的模态空间中利用模态叠加法得到部分解耦的动力学方程，结合考虑质量变化效应的虚拟变形法和Duhamel积分，构建反演动态力的两种分析途径。一种是基于加速度反演的连锁递推格式，另一种是基于位移和加速度反演的正则化求解格式。分析两个动载荷作用下的变质量多自由度系统，利用Euler中点辛差分方法求得系统的响应，将其叠加随机噪声以模拟测点响应，然后，应用两种模态力反求格式反演动载荷。结果对比表明，连锁递推格式计算效率高，但对加速度测量噪声敏感；正则化求解格式的稳健性较高。 关键词：变质量系统，载荷识别，模态力，递归计算，正则化 中图分类号： O326, O347.1 文献标识码：A 引 言 许多工程领域存在质量等振动参数随时间发生变化的动力学问题 [1-5]。质量时变系统的动力学响应分析方法主要包括两类，一是直接积分法，包括Newmark法和各类差分算法 [6]；另一种是基于短时时不变假设的模态叠加法[7,8]。与动力响应分析相对应的是质量时变振动系统的载荷反求，如考虑质量快速减少的运载工具的脉动推力反求，这在工程上也是非常有意义的课题。针对时不变结构的动态载荷识别，已经发展了较为丰富的方法与技术[9]。然而，关于时变结构的动态载荷识别方法的研究还较少，张青霞和段忠东等[10]将车桥系统视为质量时变系统，移动荷载的反求转化为移动质量的识别，基于虚拟变形法识别移动移动质量。本文研究质量时变系统载荷识别方法，在初始模态空间中，结合模态叠加和虚拟变形技术，构造模态力反求格式，经模态转换得到相应的物理量。分析变质量多自由度系统，首先由差分法计算受激振动响应，将其叠加噪声以模拟测量信号，然后，反求模态力得到动载荷，比较两种反求格式的有效性及优劣。 1 初始模态空间中的模态力计算格式 1.1基于初始模态求解时变质量系统的动力学方程 考虑自由度数为n的变质量系统，其动力学方程表示为[7] (1) 式中，M为时变的质量矩阵，F(t)是外激励力，是质量矩阵的时间变化率，ur为变化的质量相对原质点的速度矢量。式(1)右边第二项是系统质量变化获得的作用力，如火箭推力。不失一般性，不考虑该项，假设质量减少或并入的相对速度是零。由此，式(1)写为 (2) 式中，M0为初始质量矩阵，△M(t)为与时间相关的增量质量矩阵，即 (3) 考虑初始振动参数构成的自由振动系统，有 (4) 采用比例阻尼模型，C =αM0+βK。将基于式(4)得到的振型空间称为初始模态空间，取前面ns阶模态，相应的模态矩阵表示为Φ，其维数为n×ns。对式(3)所描述的振动系统，进行如下变量转换 (5) 式中，q为初始模态空间中的ns个模态坐标组成的列向量。式(5)代入式(2)，两边同时左乘ΦT，得到 (6) 式中，ωi和ξ i (i=1,2,…,ns)分别为初始瞬时的频率和模态阻尼比，diag(…)表示对角矩阵，右端各项分别为 (7) 式(7)第一式中的Q即为初始模态空间中的模态力。 式(6)右边的质量改变项若为零，对应的结构为初始结构。如果已知初始结构的响应，根据虚拟变形法， 结构模型变化后的响应可通过引入的虚拟力求出[11]。由此，在零初始条件下，可由初始系统的Duhamel积分得到式(6)的解为 (8) 式中，H(t)是单位脉冲响应函数构成的对角矩阵，H(t) = Diag(hi)，其中的对角项为脉冲响应函数 (9) 式中，。式(8)可具体展开为 (10) 下面构建由已知响应求模态力的两种反演途径。 1.2 模态力反演的连锁递推格式 引入以下参向量 (11) 由式(8)导出模态加速度为 (12) 式中，I是n阶单位矩阵， 针对式(12)构造等步长递推格式，时间步长为△t，第l步(l=0，1，2，…)所对应的时间tl=l△t，采用梯形求积法离散，可得如下递推式 l ≥1， (13a) (13b) 式中，，上标括号中的l与l-1分别表示第l和l-1递推步。结合式(12)和(11)，将当前时间步的模态力表示为 (14) 式中，， 1.3 模态力反演的正则化求解 在式(10)中，令 (15) 将式(10)的求积区间划分成长为△t的m等分，