六动力学问题的有限元法.pptxVIP

第六单元 动力学问题的有限元法;指边界条件和/或体力变化缓慢，或者物体内加速度分布均匀等类型的问题。这类变形体问题的平衡微分方程中忽略了惯性项，但载荷是时间的函数。在某时刻t，采用动静法将整体惯性力转化为体力，或者忽略惯性力。对应此刻载荷的静力学解作为t时刻的解。工程上可取随时间变化载荷的最大值的静力学解作为问题的准静态解。;通过这种近似处理，可以避免大量的动力学模型解算，而在有限的计算机资源下，可把实际问题的模型在准静态假设前提下考虑得更细致、更实用。在许多情况下，由此带来的对实际情况的逼近将大大抵消由于准静态假设产生的误差。;结构动力学问题;弹塑性动力学问题;对于上述后两类问题，描述质点平衡和运动的微分方程相同，包含惯性力项和阻尼力项。其数值求解方法主要是有限元法。;第二节 动力学问题的有限元方程;为建立有限元动力学响应控制方程，利用达朗倍尔原理，在每个时刻 t，将连续介质中质点加上惯性力 和阻尼力 ，则系统的动力学问题转化为等效静力学问题。对等效系统应用虚功原理：;上述矩阵由相应的单元矩阵组集而成：;如果忽略阻尼，则结构动力学方程简化为： ;当求解该微分方程组，得出节点位移响应后，其它计算步骤与静力分析相同。;第三节 质量矩阵和阻尼矩阵 ;对于3节点三角形单元，按上述公式计算得到的一致质量矩阵为：;实际应用中，两种质量矩阵都有应用，得到的计算结果相差不多。采用集中质量矩阵可以使计算得到简化，提高计算效率，由此得到的自振频率常低于精确解。;2、阻尼矩阵;上述两种阻尼矩阵称为比例阻尼或振型阻尼。其比例系数一般依赖于频率，很难精确确定。;由于系统的固有振型对于结构质量矩阵和结构刚度矩阵具有正交性，因此，系统振型对上述Rayleigh阻尼矩阵也是正交的。所以这类阻尼矩阵又称为振型阻尼。;第四节 结构自振频率和振型;;;;求解结构自由振动的广义特征值问题，由于系统自由度很多，而研究系统动态响应和动态特性时，往往只需要少数低阶特征值和特征向量。因此在有限元分析中发展了许多针对上述特点的效率较高的算??。其中应用最广泛的有Lanczos法、子空间迭代法、逆迭代法等。;第五节 瞬态响应分析;1、直接积分法;直接积分法概述：;大多数显式方法是条件稳定的：当时间步长大于结构最小周期的一定比例时，计算得到的位移和速度将发散或得到不正确的结果；;t+Δt时刻的位移解 从t时刻的运动方程建立： ;上式就是求离散时间点上位移解的递推公式。但该算法有起步问题（见P449）。;;在积分区间上采用如下的速度、位移假设： ;; 当算法中的参数满足一定条件时，该算法是无条件稳定 的。此时，步长的选择取决于解的精度，可以根据对结 构响应有主要贡献的若干基本振型的周期来确定。通常 可取为所要考虑的基本振型周期中最小周期的二十分之 一。;2、振型叠加法;此变换的意义： 变换把结构的瞬态位移响应从以有限元网格节点位移为基向量n维空间转换到以固有振型为基向量的n维空间。这里 看成广义位移基向量，xi 是广义位移分量。数学上看，是离散系统位移在两个不同向量空间之间的变换。;并且两边左乘 ，并考虑到 的正交性，则得到新向量空间内的运动方程：;其中ξi定义为第i阶振型的阻尼比。在此情况下，上面变换后的方程就成为n个相互独立的二阶常微分方程：;从上面的过程可以看出，振型叠加法似乎需要以求解n阶广义特征值问题为代价来获得非耦合的n个单自由度系统运动方程，以提高方程求解效率。 当模型的规模很大时，求n阶广义特征值问题的计算耗费极高。事实上，对于实际问题，通常只需要利用前面较少数目的振型进行坐标变换就能得到足够满意的解。也就是说只需要考虑前面部分振型对解的贡献。高阶振型对结构响应的影响通常较小。 此外，对于非线性问题通常不能采用振型叠加法，因为结构刚度矩阵随时间变化，系统特征解也随时间变化。