机械与运载工程学院显式算法与隐式算法比较显式算法平衡方程.ppt

Apr 2007 显式非线性动力学分析 显式算法与隐式算法的区别 ABAQUS/Explicit 适用问题类型 显式算法与隐式算法比较 显式算法与隐式算法比较 ABAQUS/Explicit 自动时间增量与稳定性 阻尼 ABAQUS/Explicit分析接触 问题描述 * * 机械与运载工程学院 机械与运载工程学院 * 机械与运载工程学院 * 机械与运载工程学院 机械与运载工程学院 显式算法的特点 显式方法需要小的时间增量 时间增量只与模型的最高自然频率相关，与载荷类型和载荷持续时间无关。 通常模拟需取10,000到1,000,000个增量，但是每个增量步内的计算成本相对较小。 隐式方法对于时间增量没有内在的限制 增量大小一般由精度和收敛性控制。 与显式方法相比，隐式模拟中的增量步成数量级的减少。但因必须在每个增量步内求解系统总体方程组，隐式方法在每个增量步中的成本比显式方法大很多。 隐式算法的特点 高速动力学事件 捕获瞬间动力响应，需精确跟踪应力波（爆炸） 利用显式动力学方法建立接触比隐式方法容易得多 显式算法特别适合分析受冲击载荷并在系统内发生复杂相 互接触的结构的瞬态动态响应。（手机跌落） 复杂接触问题 复杂后屈曲问题 结构刚度随载荷施加发生剧烈变化，通常后屈曲响应中还包含接触（薄壁筒的压溃） 高度非线性准静态过程 材料退化和失效 显式算法 平衡方程 隐式算法 平衡方程 求解动力平衡方程组，得加速度、速度、位移 对非线性问题，位移增量的修正值需要解如下瞬时方程 最大时间步长由数值稳定性确定 基于系统最高阶固频的稳定性限制 稳定性限制 基于单元最高阶固频的稳定性限制（保守估计） 膨胀波通过由单元特征长度定义的距离所需时间 提高稳定极限的策略 质量缩放 针对单元尺寸很小或形状极差的单元，通过增大这些单元的质量（密度），可显著增加稳定极限，而对整体响应影响可以忽略。 ABAQUS设置：缩放因子。 网格技术 尽可能采用均匀的网格 材料对稳定极限的影响 弹塑性材料屈服后刚度降低，波速减小，稳定极限增加 质量比例阻尼——减弱低频振荡 刚度比例阻尼——减弱高频振荡 材料阻尼（Rayleigh阻尼） 限制数值振荡或为系统施加物理阻尼 从属表面网格的细化，避免主面过渡侵入 网格细化 需注意的几个问题 避免初始过盈 非常小的过盈量，都可能导致节点加速度极大，导致分析初始阶段导致网格严重扭曲。 电路板跌落模拟(Explicit) 学习要点： （1）显式非线性动态分析；（2）显式接触定义；（3）装配技术 Compression yield stress ratio ： 泡沫材料屈服应力-塑性应变数据 （Foam Hardening） Hydrostatic yield stress ratio ： Deshpande,??V.?S., ?Fleck N.?A., Isotropic Constitutive Model for Metallic Foams, Journal of the Mechanics and Physics of Solids 2000, 48:1253–1276. 名义应力：变形前单位面积上的力；名义应变：变形前单位长度的长度变化。 变性前长度l0，当前长度l。体积不可压缩l0A0=lA. 名义应力：变形前单位面积上的力；名义应变：变形前单位长度的长度变化。 变性前长度l0，当前长度l。体积不可压缩l0A0=lA. 名义应力：变形前单位面积上的力；名义应变：变形前单位长度的长度变化。 变性前长度l0，当前长度l。体积不可压缩l0A0=lA. 名义应力：变形前单位面积上的力；名义应变：变形前单位长度的长度变化。 变性前长度l0，当前长度l。体积不可压缩l0A0=lA.