六、-动力学问题的有限元法

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第六单元 动力学问题的有限元法,第一节 变形体动力学问题概述,变形体动力学问题在工程和科学问题中特别普遍。该类问题由随时间变化的载荷或边界条件产生。,这类动力学问题涉及的对象包括各种机械零部件、工程构造、弹性介质。,依据问题的特点和载荷及受力体的动态特性，一般意义上的变形体动力学问题按如下三个途径处理。,指边界条件和/或体力变化缓慢，或者物体内加速度分布均匀等类型的问题。这类变形体问题的平衡微分方程中无视了惯性项，但载荷是时间的函数。在某时刻t，承受动静法将整体惯性力转化为体力，或者无视惯性力。对应此刻载荷的静力学解作为t时刻的解。工程上可取随时间变化载荷的最大值的静力学解作为问题的准静态解。,尽管这种静态状况在实际上并不存在，但作为一种根本力学模型，在工程实践上具有重要意义。很多实际问题可近似归入准静态问题，而满足工程上的精度要求。,准静态问题,通过这种近似处理，可以避开大量的动力学模型解算，而在有限的计算机资源下，可把实际问题的模型在准静态假设前提下考虑得更细致、更有用。在很多状况下，由此带来的对实际状况的靠近将大大抵消由于准静态假设产生的误差。,至于哪些问题可作准静态来处理，需要综合考虑分析目的与精度要求，构件的尺度和动态特性固有振动周期，载荷的特性上升前沿和作用时间，计算机资源状况等。,构造动力学问题,该领域争论以下问题：弹性构造系统的自由振动特性频率和振型分析；瞬态响应分析；频率响应分析；响应谱分析等。,就构造的瞬态响应分析而言，典型的有构造在冲击载荷下的响应问题。构造动力学中这类问题的特点是，载荷作用前沿时间与构件的自振基频周期相近，远大于应力波在构件中的传播时间。或者构件上长时间作用随时间猛烈变化的载荷。,构造动力学问题在工程中具有普遍性。,弹塑性动力学问题,这是连续介质变形体动力学问题的另一个重要领域。涉及很多科学和工程领域，如高速碰撞，爆炸冲击，人工地震勘探，无损探伤等。,这类问题的争论要深入到介质中的弹塑性波的传播过程以及考虑波动效应前提下介质中应力应变的响应。,这类问题中载荷的特点是构件上载荷作用前沿时间远少于应力波在构件中的传播时间。该状态通常由构件高速碰撞或爆炸载荷产生。,对于上述后两类问题，描述质点平衡和运动的微分方程一样，包含惯性力项和阻尼力项。其数值求解方法主要是有限元法。,其次节 动力学问题的有限元方程,在连续介质的动力学问题中，描述力学参量的坐标是四维：3个空间坐标和一个时间坐标。进展有限元法求解时，只对空间区域进展离散化，得到离散多自由度系统的动力学模型。,其有限元法步骤与静力学问题一样。只是在单元上对随时间变化的节点位移进展插值，得到单元内随时间变化的假设位移场：,为建立有限元动力学响应控制方程，利用达朗倍尔原理，在每个时刻,t,，,将连续介质中质点加上惯性力 和阻尼力 ，则系统的动力学问题转化为等效静力学问题。对等效系统应用虚功原理：,将前面位移空间离散表达式和单元的几何方程、物理方程代入上式虚功方程，并考虑到变分的任意性，得到离散系统掌握方程构造有限元动力学方程：,方程中的系数矩阵分别为：系统质量矩阵，阻尼矩阵，整体刚度矩阵。右端项为整体节点载荷向量。,上述矩阵由相应的单元矩阵组集而成：,其中：,单元质量矩阵,单元刚度矩阵,单元阻尼矩阵,单元等效节点力向量,假设无视阻尼，则构造动力学方程简化为：,上式动力学方程的右端项为零时就得到构造自由振动方程。,从动力学方程导出过程可以看出，动力学问题的有限元分析中，由于平衡方程中消失了惯性力和阻尼力，从而引入了质量矩阵和阻尼矩阵，运动方程是耦合的二阶常微分方程组，而不是代数方程组。该方程又称为有限元半离散方程，由于对空间是有限元离散的，对时间是连续的。,当求解该微分方程组，得出节点位移响应后，其它计算步骤与静力分析一样。,有限元动力学方程的求解虽然可以承受常规的常微分方程组解法，但由于实际问题有限元模型的阶数往往很高，用常规方法不经济，通常承受一些对有限元方程有效的解法，主要分为两类：直接积分法和振型叠加法。,第三节 质量矩阵和阻尼矩阵,1、协调质量矩阵和集中质量矩阵,该矩阵称为协调质量矩阵或全都质量矩阵。由于它和刚度矩阵依据同样的原理、过程和插值函数导出，还表示质量在单元上呈某种分布。,此外，有限元中还常常承受集中质量矩阵，它是一个对角矩阵，由假定单元质量集中在节点上得到。,上节导出的单元质量矩阵为：,对于3节点三角形单元，按上述公式计算得到的全都质量矩阵为：,该单元的集中质量矩阵为：,实际应用中，两种质量矩阵都有应用，得到的计算结果相差不多。承受集中质量矩阵可以使计算得到简化，提高计算效率，由此得到的自振频率常低于准确解。,在波传播问题和高速瞬态非线性分析中，通常承受显式动力学求解方法协作使用线性位移单元和集中质量阵。,2、阻尼矩阵,单元阻尼矩阵：,称为协调阻尼矩阵。这种阻尼是由阻尼力正比于质点运动速度得到的，属于粘性阻尼。明显，这种阻尼阵与质量矩阵成正比。,对构造而言，阻尼并非粘性的，而主要是由于材料内部摩擦效应引起的能量耗散，但这种耗散机理尚未完全清晰，更难以用数学模型表达，故通常假设这种状况的阻尼力正比于应变速率，从而可导出比例于单元刚度矩阵的单元阻尼阵，大多数情形下足够准确。,上述两种阻尼矩阵称为比例阻尼或振型阻尼。其比例系数一般依靠于频率，很难准确确定。,一个通行的方法是将构造的阻尼矩阵简化为构造刚度阵和构造质量阵的线性组合：,其中，是不依靠于频率的常数。这种振型阻尼称为Rayleigh阻尼。当=0时，较高阶振型受到的阻尼较大；当=0 时，较低阶振型受到的阻尼大。,由于系统的固有振型对于构造质量矩阵和构造刚度矩阵具有正交性，因此，系统振型对上述Rayleigh阻尼矩阵也是正交的。所以这类阻尼矩阵又称为振型阻尼。,承受振型阻尼矩阵后，可以利用系统振型对动力学方程进展变换，得到解耦的方程组，使每个方程可以独立求解，给计算带来便利。,第四节 构造自振频率和振型,其中 是n阶向量，表示有限元离散结构所有自由度的振幅，是该向量振动的频率。将上式代入自由振动方程得到：,争论构造自由振动特性。设阻尼和外力均为零，则构造自由振动有限元运动方程为：,设各自由度作简谐运动：,该方程描述的问题称为广义特征值问题。,其中特征值：代表系统的n个固有频率，并有：,特征向量：代表系统的,n,个固有振型，或称为主振型。其幅度是不确定的，但可以用下列方法对其正则化：,求解该问题可以得到n对特征解特征对,这样的固有振型又称为正则振型，以后默认所谓固有振型即是指这种正则振型。,简洁证明，固有振型具有对M和K的正交性：,定义：,它们分别称为固有振型矩阵和固有频率矩阵,利用固有振型矩阵和固有频率矩阵，构造固有振型的正交性质可以表示成：,原来的特征值问题可以表示成：,固有频率和固有振型是一个构造自由振动的根本特性，也是构造动态特性的根本要素。,求解构造自由振动的广义特征值问题，由于系统自由度很多，而争论系统动态响应和动态特性时，往往只需要少数低阶特征值和特征向量。因此在有限元分析中进展了很多针对上述特点的效率较高的算法。其中应用最广泛的有Lanczos法、子空间迭代法、逆迭代法等。,第五节 瞬态响应分析,瞬态响应分析是计算动力强迫响应分析的最一般方法。其目的是计算构造受随时间变化鼓励作用下的行为。瞬态鼓励定义在时间域中，每个瞬时的大小。鼓励可以是作用力和强迫运动。,依据构造和载荷的性质，可以用两种不同的数值方法进展瞬态响应分析：直接积分法和振型叠加法。前者对全耦合的有限元离散运动方程直接进展积分；后者利用主振型对运动方程进展变换和解耦，构造的响应依据相应于各振型的响应累加而成。,1、直接积分法,第一，将求解时间域0tT内任何时刻t都满足运动方程的要求降低为在相隔t的离散时间点上满足运动方程。,直接积分法的两个前提：,其次，在离散时间点之间的t区域，对位移，速度，加速度进展假设。相当于对运动微分方程组在时间域进展离散化，并逐点求解。,直接积分法概述：,直接积分法的时间离散化方程有显式和隐式两类。在显式方程中，由,t,时刻的运动方程求,t+t,时刻的位移；而隐式方法是从与,t+t,时刻运动方程关联的表达式中求,t+t,时刻的位移。,显式方法要求很小的时间步长，但每步求解所需计算量较小；而隐式方法允许较大的时间步长，但每一步求解方程的消耗较大。,大多数显式方法是条件稳定的：当时间步长大于构造最小周期的肯定比例时，计算得到的位移和速度将发散或得到不正确的结果；,隐式方法往往是无条件稳定的，步长取决于精度，而不是稳定性方面的考虑。,典型的显式方法是所谓的“中心差分法”，其根本思想如下。,t+t,时刻的位移解 从,t,时刻的运动方程建立：,中心差分法,将某时刻的加速度和速度用中心差分表示：,将加速度和速度的差分格式代入上式，得到：,上式就是求离散时间点上位移解的递推公式。但该算法有起步问题见P449。,中心差分法特点如下：,1是显式算法，并且当质量阵和阻尼阵都是对角阵时，利用该递推公式求解运动方程时不需要进展矩阵求逆，这个特点在非线性问题中将更有意义。,是有限元系统的最小固有振动周期，通常用最小尺寸单元的最小固有振动周期代替。因此，有限元网格中最小单元尺寸将决定中心差分法时间步长的选择。有限元网格划分时要考虑到这个因素，避免个别单元尺寸太小。,2）是条件稳定算法。时间步长必须小于某个临界值：,3中心差分法适合用于考虑波传播效应的线性、非线性响应分析。但是对于构造动力学问题中的瞬态响应分析，不适合承受中心差分法，由于这类问题，重要的是较低频的响应成分，允许承受较大的时间步长。通常承受无条件稳定的隐式算法。,应用最广泛的一种,隐式算法是,Newmark,方法。,在积分区间上承受如下的速度、位移假设：,Newmark,方法,通过,t+t,时刻的运动方程来决定,t+t,时刻的位移解 ，即：,从上面三个方程联立可推出从,t,时刻的运动参量计算,t+t,时刻位移的公式：,由于从上式求解,t+t,时刻位移时需要对非对角的等效,刚度阵求逆，因此称为隐式算法。,当算法中的参数满足肯定条件时，该算法是无条件稳定,的。此时，步长的选择取决于解的精度，可以依据对结,构响应有主要奉献的假设干根本振型的周期来确定。通常,可取为所要考虑的根本振型周期中最小周期的二特别之,一。,对构造动力学问题，所关心的较低阶振型的周期比全系统的最小周期大得多，也就是无条件稳定的隐 式算法可以承受比有条件稳定的显式算法大得多的时间步长，而承受较大时间步长还可以滤掉不准确的高阶响应成分。,2、振型叠加法,振型叠加法是计算构造瞬态响应的另一种数值方法。该方法利用构造固有振型对动力学方程组进展变换，缩减未知量规模，并对运动方程组进展解耦，大幅度提高数值求解的效率。,由于在构造动力学分析过程中，计算固有频率、固有振型是评价构造动态特性之必需，故振型法响应分析是常规模态分析的自然延长。,此变换的意义：,变换把结构的瞬态位移响应从以有限元网格节点位移为基向量n维空间转换到以固有振型为基向量的n维空间。这里 看成广义位移基向量，x,i,是广义位移分量。数学上看，是离散系统位移在两个不同向量空间之间的变换。,振型叠加法的根本思想如下。,首先引入变换：,其中：,并且两边左乘 ，并考虑到 的正交性，则得到新向量空间内的运动方程：,将上述变换代入有限元运动方程,假设方程中的阻尼矩阵是振型阻尼阵，利用主振型的正交性可得：,成为对角方阵。,则,其中i定义为第i阶振型的阻尼比。在此状况下，上面变换后的方程就成为n个相互独立的二阶常微分方程：,在求出每个振型坐标上的位移重量xi后，再按前面的变换关系得到节点自由度上的位移响应：,从上面的过程可以看出，振型叠加法似乎需要以求解n阶广义特征值问题为代价来获得非耦合的n个单自由度系统运动方程，以提高方程求解效率。,当模型的规模很大时，求n阶广义特征值问题的计算消耗极高。事实上，对于实际问题，通常只需要利用前面较少数