第9章动力学有限元.pptVIP

计算结构的固有频率和振型是结构动力学分析的主要内容，也是分析结构动力响应和其它动力特性问题的基础。由于一般结构阻尼对结构的固有频率和振型影响极小，所以，求结构的固有频率和振型时，直接用无阻尼的自由振动方程求解。即因任意弹性体的自由振动都可分解为一系列的简谐振动的迭加：即结构上各节点位移为为节点位移振幅向量（即振型），与时间t无关的位移幅值；为与该振型对应的频率。9.3结构无阻尼自由振动

第31页，共74页，星期日，2025年，2月5日1、固有频率和振型计算将节点位移代入动力方程，化简得广义特征值问题：由于结构自由振动时，各个节点的振幅不可能全为零，则称为结构的特征方程，即求结构的固有频率和振型归结为特征值问题。设计结构的自由度为n，则特征方程为的n次代数方程，其n个根称为特征值，记为它们的平方根称为系统的固有频率，即将这些固有频率从小到大依次排列为最低的频率称为基频，它是所有频率中最重要的一个。第32页，共74页，星期日，2025年，2月5日这个过程称之为正规化利用正规化，可得2、特征向量对应每个固有频率，可有方程由此求得一组节点振幅不全为0的向量称为特征向量，也称为振型或模态向量。由于上述方程为齐次方程，显然解不唯一，也就是说：振型的形状是唯一的，但其振幅不是唯一的；或一个特征值可对应有多个特征向量，但一个特征向量只对应一个特征值。实际中，常选特征向量使第33页，共74页，星期日，2025年，2月5日则对应所有的特征值问题:3、特征向量的性质正交性：任意两个特征值对应的特征向量关于质量矩阵或刚度矩阵正交。即设则有若将所有的特征值对应的特征向量组装成特征向量矩阵，即第34页，共74页，星期日，2025年，2月5日考虑到正规化:可进一步记为：可简记为矩阵形式：第35页，共74页，星期日，2025年，2月5日1、幂迭代法特点：用于计算最大（主）特征值十分有效。这里[D]称为动力矩阵，也即一个变换矩阵，它可将任一特征向量变换为一常数与其自身的乘积.9.4特征值问题的解法结构固有频率和振型的计算归结为求的特征值和特征向量。由于有限元法将结构离散为n个自由度，n一般相当大，故n次特征方程的直接求解十分困难，常求其近似解，常用的求解方法有幂迭代法、逆迭代法、子空间迭代法等。第36页，共74页，星期日，2025年，2月5日由于任两个特征值对应的特征向量是正交的，则n个特征向量可组成特征向量空间中的一个特征向量基，其特征向量空间中的任一特征向量可表示为基向量的线性组合。即存在任一向量：设这个向量被[D]变换后形成一新的特征向量为：类推，可得：第37页，共74页，星期日，2025年，2月5日由于所有的特征值排列为：即存在考虑到问题为齐次方程，特征向量前的系数可以略去，则上式在p趋近无穷时，其第一项就趋近实际计算，只需迭代有限次即可得精确解。第38页，共74页，星期日，2025年，2月5日幂法迭代格式1、选初始特征向量，如单位向量2、构造新特征向量，并归一化3、计算特征值近似值4、计算相邻两次迭代的特征值误差，检查是否收敛若需计算二阶、三阶等特征值，则需构造新的动力矩阵第39页，共74页，星期日，2025年，2月5日2、逆迭代法逆迭代法也称为反幂法，类似于幂法，特征值问题改写为：其具体迭代格式为：1）选初始向量如单位向量2）计算中间向量3）求解线性方程组4）归一化5）计算特征值近似值6）计算相邻两次迭代的特征值误差，检查是否收敛第40页，共74页，星期日，2025年，2月5日9.5动力响应的计算对于受迫振动，基本方程为求解此方程通常有两种数值方法：振型迭加法和逐次积分法1、振型迭加法振型迭加法的基本思想是利用结构固有振型的正交性，把结构的复杂振动分解为一组相互独立的单自由度振动（即解耦），从而求得结构的位移响应。设结构无阻尼自由振动的各阶固有频率和相应的固有振型为：则结构任意时刻的受迫振动产生的位移可认为是n个固有振型为基的线性组合，即为组合系数，是时间t的函数，也称为振形坐标第41页，共74页，星期日，2025年，2月5日广义质量阵广义阻尼阵广义刚度阵广义激振力上式可记为这里代入动力学方程：左乘第42页，共74页，星期日，2025年，2月5日据正交性可知，这些广义矩阵均为对角矩阵，即表示方程各个变量之间是没有耦合项的，从而动力方程