动态载荷识别 模型修正与结构动力修改.pptVIP

* 研究有阻尼情形下能将复模态简单转化到实模态空间进行拓扑重分析的方法，从而建立实、复模态空间都适用的研究思路和方法体系；提出Benchmark算例，以便建立统一的重分析评价标准，编制拓扑重分析的程序或者软件包，嵌套到结构拓扑优化设计软件中，从而使拓扑重分析方法可以真正解决实际问题。 * J.J. He, J.S. Jiang, B. Xu. Modal reanalysis methods for structural large topological modifications with added degrees of freedom and non-classical damping. Finite Elements in Analysis and Design [J]. 2007 (44):75– 85 (SCI, EI收录) 2. 何建军，姜节胜，康兴无, 结构拓扑大修改的动力学重分析方法[J]. 振动工程学报， 2007，20(4):407-411 (EI收录) 3. 何建军, 姜节胜.结构拓扑修改动力学重分析的单步摄动逆迭代法[J]. 西北工业大学学报，2006, 24(3):313-316 (EI收录) 4. 何建军, 姜节胜.基于降阶模型的非经典阻尼结构拓扑修改的复模态重分析 方法.振动与冲击，2009, 28(7): 50-54(EI收录) 发表的动力学拓扑重分析的文章 结构动力学优化设计 结构动力特性修改的方法有很多种，目前已发展起来的有矩阵摄动法，加权欧式范数法， 传递函数法，及灵敏度法等等。本节着重讨论结构动力修改反问题的灵敏度方法。 灵敏度方法是建立在结构特征灵敏度分析的基础上，运用多元函数的泰勒展开式来确定结 构特性参数的改变量。 设特征值λr与特征向量Ψr均为结构参数mij，kij和cij的多元函数 式中mij，kij和cij分别为质量矩阵M，刚度矩阵K和阻尼矩阵C中第i行，第j列元素。将上式 展开成泰勒级数。 在实际计算时，当结构参数的修改量较小时，常略去二阶修正项。此时特征值向量为 式中的导数均对原特征值 。 根据一阶灵敏度公式，在结构参数修改量Δmij 、 Δkij 、 Δcij 确定后，即可由上式 求出特征值的修正量，从而求得修改后结构的特征值。 同样，对特征向量亦可得类似上式的公式，在略去修正量后，可得 式中的导数均对 取值。 有了上面两式，就可以在已知灵敏度的前提下，根据结构参数的改变量 Δmij 、Δkij 、Δcij求修改后结构的特征值及特征向量的改变量Δλr及ΔΨr ，求 Δmij 、Δkij 、Δcij，这是结构动力修改的逆问题。 对各阶特征值λr（r=1，2，…，N）都可以写出类似的公式，对各阶特征向量中每个元素 φir（i=1,2 … ，N；r=1,2， … N）。都可以写成类似公式。这些方程组便可求解结构 动力修改的正问题及逆问题。 当结构参数的改变量较大时，可采取分步计算。 【例】如图所示，一个4自由度，无阻尼系统，各质量及刚度如下： 其质量矩阵及刚度矩阵分别为 其各阶固有频率及振型分别为 1） 修改固有频率 要求系统的第一阶固有频率由原来的13.294rad/s改变到15.00rad/s，而将第一个质量m11作为修改对象，求出其修改量Δ m11。 由式 式中， 为原系统的第一阶固有频率及振型向量中的第一个元素。 已知 ，为了提高计算精度，将分成五步计算， 每次取 。分步计算时按如下公式迭代 假设在每一步计算时，振型不变。质量 总修正量为 计算结果的每次修改为 总质量修改量 经过修改后总质量变为 对修改后的系统重新计算其第一阶固有频率为15.005rad∕s，可见误差是极小的。 2）修改振型 由给出的原始振型数据可见，该系统的第二阶振型的节点位于第一点与第二点点之间。要求将节点位置移到第2点，仍以m11作为修改对象。在此情况下，修改的目标使（α2）2=0，即为第二阶振型向量中的第二个元素为零 修改公式为 设 当s≠r，i=j时， 当s=r，i=j时， 故 故 因此，修改后第一个