9页矩阵特征值反问题ch3–7.docVIP

3.7 线性(谱)约束下矩阵束最佳逼近问题 3.7.1 问题的提法 个自由度无阻尼线性振动系统的动力学方程为 ， (3.7.1) 其中通常分别称为质量矩阵和刚度矩阵，分别称为系统的位移和对系统的作用力。结构动力分析首先考虑系统的固有动力特性，即在无外力（）作用下系统的固有振动特性。为此，需要研究如下齐次微分方程组解的特性 ， (3.7.2) 假定矩阵与时间无关，则常系数微分方程组(3.7.2)的解（即主振动的一般形式）为 ， (3.7.3) 其中称为系统的固有频率，称为系统的振型。将(3.7.3)代入(3.7.2)，即得 ， (3.7.4) 记，则(3.7.4)即为如下广义特征值问题 ， (3.7.5) 结构振动的正问题就是由结构的质量矩阵和刚度矩阵求结构的特征值和相应的特征向量（振型），然后利用求得的特征值和特征向量将微分方程组(3.7.1)解耦，最终求得系统的解（响应）。 用有限元方法可得到系统的质量矩阵和刚度矩阵，通常称和分别为分析质量矩阵和分析刚度矩阵。有限元方法对结构进行离散化必须对结构部件之间的连接条件、边界条件等作力学上的简化。因此，有限元方法所得到的分析质量矩阵和分析刚度矩阵与实际结构的质量矩阵和刚度矩阵往往存在一定的误差。通常用有限元方法建立的分析质量矩阵和分析刚度矩阵计算结构的低阶频率和振型（工程力学家认为：用有限元方法求得系统的三分之一低阶频率是可信的），但如果直接用有限元方法所得到的分析质量矩阵和分析刚度矩阵求系统的响应则误差较大。 最近三十年来，模态参数（频率和振型）测试技术得到了很大发展，复杂结构模态参数的实测已不再是一件耗资巨大、十分困难的事情。因此，利用结构的实测模态参数（低阶部分频率和相应的振型）对有限元方法所得的分析质量矩阵和分析刚度矩阵进行修正，使得修正模型的动态特性至少低阶部分与实测结果相一致，再利用修正模型求系统的响应。三十余年的工程实践表明，这样做不仅可行，而且确实有效。 如果系统的质量矩阵和刚度矩阵已知，且非奇异，由(3.7.4)求得结构的动态特性，即系统的固有频率和相应的振型，则 (3.7.6) 记 ， (3.7.7) 和分别称为频率矩阵和振型矩阵，则 ， (3.7.8) 如果质量矩阵对称正定，刚度矩阵对称，则可求频率矩阵和振型矩阵使得 ， (3.7.9) 或 。 (3.7.10) 因此，结构动力学模型的质量矩阵、刚度矩阵、频率矩阵和振型矩阵理论上应满足关系式(3.7.8)－(3.7.10)。通常称(3.7.8)为特征方程，(3.7.9)或(3.7.10)为正交关系。 由振动试验可测得结构系统的低阶个固有频率和相应的振型，记 ， (3.7.11) 如果和是结构动力模型的质量矩阵和刚度矩阵，则和至少应满足下列三个方程之一： ， (3.7.12) ， (3.7.13) 。 (3.7.14) 结构有限元模型修正的提法是：基于振动试验的实测数据(3.7.11)是可靠的，对结构有限元模型作必要的修正，使修正模型的前阶模态数据与实测模态数据一致。因此，结构有限元模型修正可概述为如下广义特征值反问题。 问题GI 给定矩阵，其中是块对角矩阵，其对角块是1阶或2阶实矩阵，，求，使得和满足矩阵方程(3.7.12)－(3.7.14)的一部分。 问题GII 给定，求，使得 ， (3.7.15) 其中是问题GI的解集合。 问题GI与问题GII称为谱约束下矩阵束最佳逼近问题。问题GI与问题GII本质上是一个约束最优化问题，工程技术人员通常用Lagrange乘子法求解结构有限元模型修正问题，读者可参阅文献［21-23, 85-93］。本节将用代数方法研究结构有限元模型修正中的数学问题。 3.7.2 谱约束下实矩阵束最佳逼近问题 考虑如下线性约束下矩阵束最佳逼近问题［94］： 问题3.7.1 给定，求使得 。 (3.7.16) 问