ABAQUS显式动力学仿真方法及相关实例.docx

ABAQUS显式动力学仿真方法及相关实例——友荣一、动力学显式有限元方法显式时间积分Explicit应用中心差分法进行运动方程时间积分，由一个增量步的动力学条件计算下一个增量步的动力学条件。求逆，显式算法总是采用一个对角的或者集中的质量矩阵，不必同时求解联立方程。任何节点的加速度完全取决于节点质量和作用在节点上的合力，计算成本非常低。中心差分方法：假定加速度为常数，应用这个速度的变化值加上前一个增量步中点的速度来确定当前增量步中点的速度：同理，速度对时间的积分加上在增量步开始时的位移可以确定增量步结束时的位移：即：增量步开始时提供了满足动力学平衡的而加速度。之后，在时间上“显式地”前推速度和位移。所谓显式，即增量步结束时的状态仅依赖于该增量步开始时的位移、速度、加速度。为保证精确，时间增量必须相当小，在增量步中加速度几乎为常数。由于时间增量很小，典型分析需要成千上万的增量步。由于不用联立方程组，计算成本主要消耗在单元的计算上，以此确定作用在节点上的单元内力，计算成本很低。主要步骤：①动力学平衡方程： ②对时间积分： ③单元计算：根据应变速率，计算应变增量 根据本构关系计算应力 集成节点内力，再把时间变为，继续计算；显式和隐式时间积分程序，都是求解动力平衡方程的节点加速度，并应用同样的单元计算以获得单元内力。两个程序最大的不同在于求解节点加速度的方式。隐式采取直接求解一组线性方程组，在完全Newton迭代求解方法的基础上，使用自动增量步。隐式算法无条件稳定，时间增量比显式大。而对于非线性问题，隐式增量步需要经过几次迭代才能获得满足给定容差的解答，代价极大。隐式分析中，每一次迭代都需要求解大型的线性方程组，占用极大的计算资源、磁盘空间、内存。显式方法特别适用于求解高速动力学问题，它需要许多小的时间增量来获得高精度的解答。如果持续时间非常短，则可能得到高效率的解答。显式可以很容易模拟接触条件和其他一些极度不连续的情况，并且能够一个节点和一个节点地求解不必迭代。通过调整节点加速度来平衡在接触时的外力和内力。显式方法最显著的特点是不需要迭代。而稳定性在显式方法中尤为重要。稳定性限制了explicit求解器所能采用的最大步长，若时间增量大于这个最大的时间步长，即超出稳定性限制（stability Limite），可能造成数值不稳定，解答不收敛。通常不能精确确定稳定性限制，常采用保守的估计值。因为稳定性限制对可靠性和精确性有很大的影响，所以必须一致性和保守地确定这个值。为了提高效率，explicit尽量选择接近但不超过稳定性限制的时间增量。稳定性限制的定义：无阻尼：有阻尼：其中是系统内最高频率，为最高频率模态的临界阻尼部分。而实际最高频率是基于一组复合的相互作用因素，不大肯那个计算出确切值。代替方法是估算每个模型中每个个体单元的最高频率，取最大。计算方法：其中，——单元长度 ——材料波速材料波速可以通过以下公式求得：其中，——杨氏模量 ——材料密度可以说，稳定性极限是膨胀波通过由单元特征长度（取最小）定义的距离所需要的时间。在显式分析中所采用的时间增量必须小于稳定极限。对稳定性极限有影响的因素：质量缩放（密度），材料刚度变化（线性时，波速恒定，非线性时，材料屈服则刚度变小，波速减小），网格（单元长度）等等。二、仿真小钢球撞击圆盘1、问题简介半径0.05M的小钢球撞在半径0.25M、厚度0.05M的圆钢板。钢球以垂直钢板500M/S、平行钢板1000M/S的速度冲击，冲击时间4e-5S。2、步骤建模；设置材料：钢球杨氏模量2.068e11,泊松比0.3；钢板杨氏模量2.078e11,泊松比0.3，塑性1.22e9,0和2.44e9,1；更改关键字（keyword），添加数据块（shear failure，1.0）；装配，分割钢板；建立显式分析步，时间为4e-5；建立相互作用，使钢球为刚体（不重点考虑钢球）；建立钢球与钢板接触区域节点的接触关系；建立边界条件，钢板约束所有自由度，在钢板横切面施加对称约束，钢球施加速度；画网格，钢板取显式分析线性三维缩减积分单元C3D8R，钢球取显式线性三维应力单元C3D4R； 提交模拟； 后处理：1系统内能、动能曲线： 内能曲线 动能曲线内能与动能曲线出现问题：1、Conflicts is generated by Abaqus/CAE to highlight changes made via the Keywords Editor that conflict with changes made subsequently via the native GUI.解决：把含conflict的错误词条在keyword中删掉，它由CAE自动产生；2、The?independent?variables?must?be