结构动力学的有限元法.ppt

有限单元法 动力分析 固有特性分析 响应分析 固 有 频 率 振 型 位 移 响 应 速 度 响 应 加 速 度 响 应 动 应 变 动 应 力 固有特性：是一组模态参数构成，它由结构本身（质量与刚度 分布）决定，而与外部载荷无关，但决定了结构对动载荷的响应； 响应分析：是计算结构对给定动载荷的各种响应特性。 第一节 动态分析有限元法的特点 一、载荷特点 结构所受的载荷是随时间变化的动载荷。 这是与静力分析的一个根 本区别。 二、位移特点 1、节点位移{q}不仅是坐标的函数，而且也是时间的函数。仍以节 点位移{q}作为基本未知量。 2、节点具有速度 加速度。 3、利用节点位移插值表示单元内任一点的位移 一般仍采用与静力分析相同的形函数，[N]。当单元数量较多时，上述 插值可以得到较好的插值精度。 4、在线弹性条件下，单元内的应变和应力与节点位移的关系仍为 但这时的位移、应变和应力都是某一时刻的瞬时值，它们都是随时间t 变化的函数。 5、由于节点具有速度和加速度，结构将受到阻尼和惯性力的作用。 根据达朗伯原理，引入惯性力和阻尼力之后结构仍处于平衡状态，因 此动态分析中仍可采用虚位移原理来建立单元特性方程，然后再集成。 整个结构的平衡方程为 式又称运动方程，它不再是静力问题那样的线性方程，而是一个二阶 常微分方程组。 求解过程复杂，建立有限元模型时要特别注意控制模型规模。 第二节 动态分析有限元法的一般步骤 一、结构离散 该步骤与静力分析完全相同，只是应该分析内容不同，对网格形式的要求 有可能不一样。 静力分析：要求在应力集中部位加密网格； 动态分析：由于固有频率和振型主要与结构的质量和刚度分布有关，要求 整个结构采用尽可能均匀的网格形式。 二、单元分析 单元分析的任务仍是建立单元特性矩阵，形成单元特性方程。 动态分析中，单元特性矩阵：刚度矩阵、质量矩阵和阻尼矩阵。 动态分析中，仍采用虚位移原理建立单元特性矩阵。 在动载荷作用下，对于任一瞬时，设单元节点发生虚位移 ，则单元 内也产生相应的虚位移 和虚应变 。单元内产生的虚应变能为: 单元除受动载荷外，还有加速度和速度引起的惯性力 和阻尼力 ，其中ρ为材料密度，v是线性阻尼系数。外力所做的虚功为： 式中，{Pv}、{Ps}、{Pc}分别为作用于单元上的动态体力、动态面力和动态 集中力；V为单元面积；A为单元面积。 由于 且形函数仅为坐标x、y、z的函数，与时间无关，因此有 根据虚位移原理，有 代入经整理，可得单元运动方程为 式中 分别称为单元的刚度矩阵、质量矩阵和阻尼矩阵，它们就是决定单元动态性能的 特性矩阵。 称为单元节点动载荷列阵，它是作用在单元上的体力、面力和集中力向单元节 点移置的结果。 在动态分析和静力分析中，单元的刚度矩阵是相同的，外部载荷的移置原理也一样。 在动态分析中，单元的质量矩阵通常采用以下两种形式。 1、一致质量矩阵 按 形成的单元质量矩阵称为一致质量矩阵，因为 它采用了和刚度一致的形函数。这种质量矩阵取决于单元的类型和形函 数的形式。 2、集中质量矩阵 集中质量矩阵将单元的分布质量按等效原则分配在各个节点上，等效原则 就是要求不改变原单元的质量中心，这样形成的质量矩阵称为集中质量矩 阵。集中质量矩阵是一个对角阵， 集中质量矩阵：是一个对角阵，因而可简化动态计算，减小存储容量。利 用这种矩阵计算出的结构固有频率偏低。不过有限元模型本身比实际结构 偏刚，两者相互补偿，计算出的固有频率反而更接近真实值。 一致质量矩阵：由于分布较合理，因此可以求得更精确的振型，另外，整 个模型的质量分布还受网格划分形式的影响。 三、总体矩阵集成 总体矩阵集成的任务是将各单元特性矩阵装配成整个结构的特性矩阵， 从而建立整体平衡方程，即 式中，{q}为所以节点位移分量组成的n阶列阵，n为结构总自由度数； （i为节点数），称为节点载荷列阵；[K]、[M]、[C] 分别为结构的刚度矩阵、质量矩阵和阻尼矩阵。 其中[K]与静力分析中的总刚度矩阵完全相同，矩阵[M]、[C]也采用与 [K]相同的集成方式，即 矩阵[K]、[M]和[C]均为n阶对称阵。 四、固有特性分析 结构的固有特性由结构本身决定，与外部载荷无关，它由一组模态参数 定量描述。包括：固有频率、模态振型、模态质量、模态刚度和模