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

直接积分法 直接积分法直接对动力学方程中的速度和加速度进行展开，用差分表示其中的导数关系，从而将动力学微分方程转换为线性代数方程 直接积分有两大类，分别是中心差分表示的显式积分和以NEWMARK法为代表的隐式差分 显式直接积分 将分析的时间平均划分为n等分，利用函数泰勒展开，有 将上述差分格式代入动力学基本方程，有 整理后，可以得到 求解该线性代数方程，就可以得到该步的位移，并连续进入下一时间步的求解，最终得到所有结果 显式差分是条件稳定的差分格式，要使得计算误差控制在有效范围内，时间步 选取受到限制 一般而言，显式直接积分的时间步长较短，迭代次数较多 显式积分方程的系数矩阵是固定的，不随计算迭代步而变，因此，显式积分只需要对方程矩阵进行一次分解，以后各迭代步只需要做矩阵乘法就可以，而不需要每一迭代步都重头进行方程求解 综合上述特点，显式积分通常用在总的计算时间较短，但变形剧烈，应力及应变激烈波动，需要在极短的时间步分析结果的类型 冲击、碰撞、爆炸、侵彻等都需要应用显式直接积分 隐式直接积分 隐式差分格式比较多，著名的Newmark是其中的代表 Newmark法是线性加速度法的一种推广 Newmark法基本假设 将上述方程代入动力学基本方程，有 其中， 隐式差分是一个无条件稳定格式，因此，迭代时间步的选取比较自由，更多的受问题求解目的和精度的控制，而不是受到算法的限制 隐式差分的线性代数方程系数矩阵 是一个随时间迭代结果而变化的矩阵，因此，在求解过程中，每迭代步的系数矩阵都不相同 方程求解过程的每一迭代步都需要从头求解，计算耗费大量的时间，计算效率不高 综合上述，隐式积分时间步长较大，但迭代计算非常耗费时间，通常应用在分析时间较长而迭代步较少（且不至于引起分析结果不可靠）的问题，如机械机构的运动，机构受普通动载荷作用的响应等 面力－压强－应力的区别 * 主对角元置1法 总刚的处理 对角元充1 对应的行和列充0 主对角元乘大数法 已知边界条件为 M是很大的数，远大于 其它元素如M=1.0E+30 应力结果与后处理 在求解有限元基本方程得到解答后，其解是整个有限元系统的所有节点位移，需要进一步处理得到各种分析结果 与单元分析过程类似，从节点位移列向量分离出各个单元的节点位移，并根据形函数得到单元位移向量 从单元位移向量经几何矩阵可以得到单元应变矩阵，并经弹性矩阵得到单元应力 各个单元应力是按照单元内部的位移、应变计算的，周围单元应力相差甚大，需要按照一定规则进行平均，并用图形显示 应力计算及结果整理 求解整体结构平衡方程的结构节点位移 总刚用一维变带宽存储技术 常采用波阵法求解 根据节点位移求单元应力 3节点三角形单元为常应力单元 可以理解应力为单元中心处的应力 结果的整理 采用绕节点平均法求节点应力 插值法求边界节点应力 3点的应力可由2点和4点插值而得 1 2 3 4 1 3 2 4 目前国内常见的有限元系统 专业有限元分析系统 ANSYS ADINA NASTRAN ABQUS ALGOR SUPER SAP CAD软件挂带 I－DEAS软件中的有限元系统 PRO/E软件中的有限元系统 UG软件带的有限元系统 有限元软件评价要点 单元库 材料库 算法库 前后处理能力 数据接口和数据转换标准 与CAD软件UG、PRO/E、Parasolid、IGES 有限元软件之间的数据接口 IGES等数据转换标准 操作方便性 第8章 有限元分析中的若干问题 有限元模型建模准则 模型应满足平衡条件 单元应满足变形协调条件 满足边界条件 刚度等价原则 谨慎选择单元 单元网格划分 几何模型的精确性 等效节点载荷的准确性 结构的静力学、动力学等效 超单元划分单级化 边界条件处理准则 指定位移边界，直接在整体刚度矩阵形成后处理 弹性边界，如果支撑的刚度已知，则把支撑作为弹簧处理；如果支撑的刚度未知，则把一部分支撑材料也一起包括进来，共同分析 有限元能够解决静定和超静定问题，但对于机构类的不足约束问题，静力有限元不能处理，必须要补充施加足够的约束 对剖分结构，施加的附加约束必须合理，不能引起额外的约束力，即不改变结构的力学形态 减小解题规模的常用措施 对称性与反对称性的合理应用 对于对称的几何结构，如果载荷也是对称的，则为对称问题 如果载荷是反对称的，则为反对称问题 对称或反对称问题，只需要分析一半结构，这样，可以减小一半的计算工作量 当分割对称或反对称面时，需要在分割出来的面上施加对称或反对称约束 如果模型结构对称，载荷不对称，可以将载荷看作对称与反对称两部分的叠加，