现代机械设计方法第2版教学课件作者谢里阳主编第2篇第5-9章节课件幻灯片.ppt

响应问题常见求解方法 振型叠加法 此法适用于外力只激起较少振型的问题（如地震） 或所需计算响应时间较长的情况 逐步积分法 适用于激发较多的振型（如冲击） 适用于计算响应时间短促的情况 振型叠加法 结构的模态包含各阶固有频率和对应的振型 振型关于结构的质量矩阵和刚度矩阵正交 振型构成n维空间的基向量，位移向量可由振型线性表示 其中，振型可以取全部叠加，也可以取前S个振型进行叠加，称为模态截断 高阶振型对系统影响越来越小 对于系统运动初始条件 和 ，可相应地用模态坐标表示 对有限元动力学基本方程 作模态坐标变换，并左乘 ，注意到振型的正交性，有 其中，阻尼系数 可以由比例阻尼和振型阻尼确定 求解解耦的动力学方程，并进行模态叠加，可以得到系统的动力学响应 直接积分法 直接积分法直接对动力学方程中的速度和加速度进行展开，用差分表示其中的导数关系，从而将动力学微分方程转换为线性代数方程 直接积分有两大类，分别是中心差分表示的显式积分和以NEWMARK法为代表的隐式差分 显式直接积分 将分析的时间平均划分为n等分，利用函数泰勒展开，有 将上述差分格式代入动力学基本方程，有 整理后，可以得到 求解该线性代数方程，就可以得到该步的位移，并连续进入下一时间步的求解，最终得到所有结果 显式差分是条件稳定的差分格式，要使得计算误差控制在有效范围内，时间步 选取受到限制 一般而言，显式直接积分的时间步长较短，迭代次数较多 显式积分方程的系数矩阵是固定的，不随计算迭代步而变，因此，显式积分只需要对方程矩阵进行一次分解，以后各迭代步只需要做矩阵乘法就可以，而不需要每一迭代步都重头进行方程求解 综合上述特点，显式积分通常用在总的计算时间较短，但变形剧烈，应力及应变激烈波动，需要在极短的时间步分析结果的类型 冲击、碰撞、爆炸、侵彻等都需要应用显式直接积分 隐式直接积分 隐式差分格式比较多，著名的Newmark是其中的代表 Newmark法是线性加速度法的一种推广 Newmark法基本假设 将上述方程代入动力学基本方程，有 其中， 隐式差分是一个无条件稳定格式，因此，迭代时间步的选取比较自由，更多的受问题求解目的和精度的控制，而不是受到算法的限制 隐式差分的线性代数方程系数矩阵 是一个随时间迭代结果而变化的矩阵，因此，在求解过程中，每迭代步的系数矩阵都不相同 方程求解过程的每一迭代步都需要从头求解，计算耗费大量的时间，计算效率不高 综合上述，隐式积分时间步长较大，但迭代计算非常耗费时间，通常应用在分析时间较长而迭代步较少（且不至于引起分析结果不可靠）的问题，如机械机构的运动，机构受普通动载荷作用的响应等 面力－压强－应力的区别 * 实际上，只要结点编号合理，这些稀疏的非零元素将集中分布在以主对角线为中心的一条带状区域内，具有带状分布的特点，如图所示。 结构刚度矩阵的带状分布图 结构刚度矩阵K具有以下性质： 对称性； 奇异性； 稀疏性； 非零元素呈带状分布。 边界约束条件处理 划行划列法 当给定某自由度位移为0时，在结构刚度矩阵中，可以将位移为0的自由度所对应的该行与该列元素从整体结构刚度矩阵中划掉，从而不出现该自由度的刚度信息；同时，在节点载荷列向量和位移列向量中，也划去该自由度对应的载荷信息。 经过划行划列处理的刚度矩阵、载荷列向量和位移列向量，都减少了一个自由度，平衡方程阶数也随之降低。 这种处理方法，简单明了，容易理解，但不容易编程实现。更主要的是，用划行划列方法处理约束，破坏了方程的节点-自由度结构，给后续分析带来一定困难。 主对角元置1法 总刚的处理 对角元充1 对应的行和列充0 主对角元乘大数法 已知边界条件为 M是很大的数，远大于 其它元素如M=1.0E+30 应力结果与后处理 在求解有限元基本方程得到解答后，其解是整个有限元系统的所有节点位移，需要进一步处理得到各种分析结果 与单元分析过程类似，从节点位移列向量分离出各个单元的节点位移，并根据形函数得到单元位移向量 从单元位移向量经几何矩阵可以得到单元应变矩阵，并经弹性矩阵得到单元应力 各个单元应力是按照单元内部的位移、应变计算的，周围单元应力相差甚大，需要按照一定规则进行平均，并用图形显示 应力计算及结果整理 求解整体结构平衡方程的结构节点位移 总刚用一维变带宽存储技术 常采用波阵法求解 根据节点位移求单元应力 3节点三角形单元为常应力单元 可以理解应力为单元中心处的应力 结果的整理 采用绕节点平均法求节点应力 插值法