有限元第五章-有限元动力学基本原理.ppt

四、机械结构动力响应的计算 1.振型叠加法 于是，振动方程解耦为 依次用求解常微分方程的方法可以得到解： 再代回 可以得到问题的解。 四、机械结构动力响应的计算 1.振型叠加法 通常情况下，高阶振型对动力响应影响较小，因此，只取最低的3～5阶（≤10）阶振型就可以得到满意精度的动力响应。 2.直接积分法 对于动力学方程 假设已知其初始条件 如果将求解时间T等分成nT个时间区间 ，并能通过前若干个时刻的解来确定下一时刻的解，就可获得问题的解，直接积分法可以解决这一问题。 四、机械结构动力响应的计算 2.直接积分法 直接积分中对动力学方程是逐步地进行数值积分的，“直接”的意思是指：进行数值积分前没有进行把方程变为另一种形式的变换。 直接积分有两个假设，一个是动力学方程的解只在相隔 的一些离散时间区间上满足方程，而不要求在任意时刻都满足方程；另一个是假设位移速度和加速度在每一个时间区间 内按一定的规律变化。 直接积分法有中心差分法、Houbolt法、Wilson- 法、Newmark和龙格-库塔法。 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 对于加速度 上式的误差为 高阶小量 对于速度 方程在t时刻为 将速度和加速度的差分格式代入方程，可以得到 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 显然，求解 ，需要2个初始条件 和 我们有初始条件 ，根据 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 可以解得 注意，有 于是，可以求得方程的解。归纳起来后，中心差分法的计算步骤为2大步，9小步。 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 第一步：初始计算 1.形成 2.计算初始值 3.选取时间步长 ,并计算积分常数 4.计算 5.形成有效质量阵: 6.做三角分解: 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 第二步：对每一时间步长 7.计算在时刻t的有效载荷 8.求解在时刻 的位移 9.如果需要，计算时刻 的速度和加速度 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 说明： 1. 的取值 是结构最高频率对应的周期，是最小周期。 2.求解方程是显式差分，对于对角线质量阵和无阻尼振动求解很方便，可以直接求解，不需进行三角分解。 四、机械结构动力响应的计算 2.直接积分法 ●中心差分法 例题：一个二自由度振动系统，其动力方程为： 已知该系统的自由振动周期为 ，试用中心差分法求解步长为 和 时，方程的解。 解：取12个步长的系统响应，假设 计算 * 第五章 有限元动力学分析基本原理 一、单元质量矩阵的计算 1.单元一致质量矩阵 2.单元集中质量矩阵 3.常用单元的一致质量矩阵 二、单元阻尼矩阵 1.速度阻尼矩阵 2.应变阻尼矩阵 三、机械结构的固有频率和振型 1.无阻尼自由振动方程 2.矩阵迭代法 3.其他方法 四、机械结构的动力响应计算 1.振型叠加法 2.直接积分法 在前面的介绍中，我们均假设作用在弹性体（或结构）上的载荷与时间无关，与此相应的，位移、应力及应变等也都和时间无关，即前面介绍的全部内容皆称结构静力学有限元方法。但工程实际中还存在着另外一类载荷与时间有关的动载荷作用于结构或弹性体，此时，相应的位移、应力、应变等都与时间有关，而且必须考虑惯性力和加速度等因素，这类分析或问题，成为动力学分析。 对于质点—弹簧系统的振动，大家比较熟悉，例如一个自由度为n的质点—弹簧振系，其动平衡方程为 第五章 有限元动力学分析基本原理 上式中每一项的含义不同 对于单元体而言，可以得到类似的上述方程 第五章 有限元动力学分析基本原理 单元质量矩阵根据其形成过程分为一致质量阵和集中质量阵，各有自身的优点和缺点。 1.一致质量矩阵 一、单元质量矩阵的计算 在离散后的结构中，取出一个单元，根据达朗贝尔原理，单位体积上作用的惯性力为： 惯性力是分布力，按分布力向节点等效的原则和实施过程，有： 一、单元质量矩阵的计算 1.一致质量矩阵 于是，令 一、单元质量矩阵的计算 1.一致质量矩阵 的计算式是通式，并因为计算质量矩阵和刚度矩阵使用的形状函数一致，因此被称为一致质量阵。 2.集中质量矩阵 在工程实际中，为了求解方便，有人