abaqus第十讲：准静态分析课件.ppt

第十讲 ABAQUS/Explicit中的准静态分析;概述;简介;显式动力学方法还可以模拟准静态问题，比如金属成型过程，但是需要特殊的考虑： 如果以自然时间周期计算，用显式动力学方法求解准静态问题是不切实际的。一般需要上百万的时间增量。;为节省计算时间，可以在模拟过程中人为的增加轧制过程的速度。 在增加轧制速度之后，静平衡问题演化为动平衡问题。惯性力的影响将会增加。 准静态分析的目标就是：在惯性力的影响较小的前提之下，尽量缩短计算的时间周期。;在模拟过程中，人为的增加准静态成型过程的速度是必要的，它可以提高解的经济性。 但是，在不使结果退化的前提下，究竟可以把速度提高多少呢？ 比如，金属成型过程中，典型的工具速度大约为1 m/s的数量级。 这个速度与金属中的典型波速相比是非常小的（钢中的波速为5000 m/s）。 一般推荐的载荷速率为材料中波速的1%。;推荐的方法： 以不同的速率多次模拟（比如，工具的速度为100, 50, 5 m/s）。 因为以低的载荷速率进行分析的时间比较长，所以以从高的载荷速度到底的载荷速度进行分析。 检查结果（变形形状、应力、应变、能量），分析不同载荷速率对结果的影响。 在显式板金成型模拟过程中，过大的工具速度将抑制起皱现象，并激起非真实的局部拉伸。 在屈曲成型过程中，过大的工具速度将引起“喷注”效应—水动力学响应（下页中有图形）。;工具速度对变形形状的影响;例子：板金 右图为汽车门标准门梁的干涉测试的简化模型。 圆梁在每个端点固定，接触刚体圆柱后变形。 测试为准静态的。;速度 25 m/s: 好的全局结果;为什么25 m/sec的速度是合适的？ 一阶频率 ( f ) 大约为 250 Hz。 相应的周期为 t=0.004 秒。 在此周期内，刚体圆柱被推向梁 d = 0.1 m。 这样，估计的速度 v 为 v = d/t = 0.1/0.004 = 25 m/sec。 金属波速为5000 m/sec，所以碰撞速度25 m/sec为波速的 0.5% 。 碰撞速度应该小于材料波速的1%。 在分析步内，以光滑缓坡的方式，把碰撞速度从零增加到所施加的碰撞速度可以得到更加精确的解。;SMOOTH STEP幅值曲线 通过逐步施加载荷的方式，可以提高准静态解的精度： 工具中的常速度条件导致金属毛坯遭受突然的冲击载荷。 这将导致应力波穿过毛坯，产生不希望发生的结果。 以斜坡的方式，使工具速度从零逐步增加，将减小这些不利的效应。 出于同样的原因，在分析结束移除工具的过程中，也以斜坡的方式将工具速度逐步降为零。;SMOOTH STEP幅值定义两个幅值之间以5阶多项式过渡。比如，在过渡开始和结束时一阶和二阶时间导数为零。 在使用SMOOTH STEP定义位移时间历程时，每个指定的幅值处的速度和加速度为零。 ;能量平衡;考虑单轴拉伸试样的拉伸测试 如果实际测试是准静态的，拉伸试件的外力功等于试件的内能。;准静态测试的能量历程显示在右图中： 惯性力是可以忽略的。 测试试件的材料速度是很小的。 动能是可以忽略的。 当测试的速度增加以后： 试件的响应偏离静态、趋于动态。 因此，材料速度和动能更加明显。;因此，能量检查为ABAQUS/Explicit金属成型过程的结果是否反应了准静态解提供了另外的评估方法。 在主要的成型过程中，变形材料的动能不可以超过内能的一小部分。 这个小部分一般为1–5%。 因为毛坯在产生显著变形之前将被移动，因此在成型过程的早期，一般很难达到这个值。 使用光滑幅值曲线将改进早期响应。 不关心工具的动能。 从模型的全部动能中减掉工具的动能，或者限制变形组件的能量输出。;例子：圆柱杯的深冲压 右图为有限元模型的1/4模型。 在所有的接触面中定义摩擦： 冲头和毛坯： m = 0.25。 砧和毛坯： m = 0.125。 毛坯夹具和毛坯： m = 0。 通过为毛坯夹具施加22.87 KN的向下力模拟深冲压过程，并将冲头向下移动36 mm。;尝试三种不同的冲压速度： 3 m/s 30 m/s 150 m/s 下表中总结了每个圆柱杯深冲压的计算费用：;最终构型中毛坯厚度云图 过大的冲压速度导致结果与实际的物理现象不符。 尽管计算费用相差10倍，以30 m/s和3 m/s的速度冲压的结果非常接近。;比较内能和动能 当冲压速度为150 m/s时，毛坯中的动能与内能相比占很大的比例。 当冲压速度为3 m/s和30 m/s 时，在成型过程中，动能与内能相比只占很小的一部分。;质量缩放;显式动力学过程稳定极限的估计公式为 其中 Le 是最小的特征单元长度， cd 是材料的膨胀波速。 泊松比为零线弹性材料的膨胀波速为 其中 E 为杨氏模量， ? 为材料密度。;如果以f 2 的方式人为的增加材料密度：