第9章 LS-DYNA运行时间控制和重启动.ppt

第九章 运行时间控制 重启动 运行时间控制 对于简单的问题（主要是理论性的），时间控制相对不重要，对于要花费大量时间的问题，时间步控制成为很重要的因素 ? 显式和隐式积分 ? 时间步计算 ? 运行时间 ? 控制时间步和运行 显式和隐式积分 一个简单的例子：弹簧质量系统 Newton第二定律 ? m a(t) + k x(t) = F(t) ? v(t) = dx/dt ? 数值积分离散微分方程为一步一步的求解过程. 在时间 已知： 和 得到: 向前推进：一旦 的值得到，计算 的值等等 显式积分—中心差分法 运动方程 ? 离散 ? 推进步 ? 稳定性 显式积分—概要 隐式积分 – 动力平衡 ?将来的加速度依赖于过去，现在和将来的位置 把 代进动力平衡方程求解： 与 无关 无条件稳定 隐式积分 – 动力平衡 显式和隐式比较 材料波速 MEDIUM 钢 铝 钛 树脂玻璃 水 空气 METERS / SECOND 5240 5328 5220 2598 1478 331 LS-DYNA 计算 梁单元 壳单元 体单元 材料波速 特征长度 体单元 ?体积/最大边的面积 壳单元 ? 面积/最大边（缺省） ? 面积/对角线 ? 面积/最小边 实体壳 ?体积/最大边的面积 梁单元 ? 长度 时间步计算 离散弹簧 ?与距离无关 ?如果所有的单元是刚体或离散的，用户则应该使用 *CONTROL_TIMESTEP设置时间步 时间步缩放因子 ?上面是基于线性分析，对于非线性，我们设置一个因子： 0.9 (缺省) 0.67 (高速) 运行时间主要依赖于： ? 材料特性 ? 网格尺寸 ? 单元数量 ? 接触 ? 计算机速度 CPU 时间估计 ? 时间步 =最小 ?循环数 = 终止时间/ ? CPU 时间 = (# 循环)(# 单元)(每一区域循环时间） ? 修改是必要的，为了降低时间步 ? 修改是必要的，对于接触的个数和尺寸 运行时间 控制运行时间 基于初始时间步尺寸百分比的侵蚀单元 避免蜕化单元（三角形、四面体等）: 非常小，高刚度，低密度 对于有较高精度要求，增加网格数量 监控 D3HSP 文件中 100 个最小时间步单元 缩放壳单元密度，维持一个最小时间步 基于初始时间步尺寸的一个百分比的终止求解 基于能量改变百分比的终止求解 基于质量改变百分比的终止求解 *CONTROL_TIMESTEP, *CONTROL_TERMINATION and *DATABASE_GLSTAT 停止和重启动能力是非常重要的 ? 重设置输出间隔 ?删除接触表面 ?删除单元和 parts ? 改变边界条件 ? 变形体和刚体转换开关 ? 控制时间步和求解时间 ? 改变阻尼选项 重启动