基于ANSYS／LS—DYNA的非线性碰撞问题仿真分析.pdfVIP

装《备制造技术)2009~年第8期 基于ANSYS／LS—DYNA的非线性碰撞问题仿真分析 佘淑华 ’。陈新连z (1．柳州市科技情报所，广西柳州 545006；2．柳州市科技局，广西柳州545006) 摘要：4r绍显式动力学算法的基本理论，通过大型商用显式动力学有限元软件ANSYS／LS-DYNA对非线性碰撞问题进行了仿真研究。 关键词：显式动力学；ANSYS／LS—DYNA；非线性；碰撞 中图分类号：TB115 文献标识码：A 文章编号：1672—545X(2009)08--0039-02 传统的有限元法，在算法上属于静态隐式算法。求解中每 乘得到。 个时间步都要重新计算大型非线性刚度矩阵，并进行反复迭 如果 为对角阵，则它的转置为三角阵，且矩阵方程可看 代以获得收敛解。需要较长的计算时间和较大的存贮空间。近 作是每一个 自由度上的独立方程组： 年来 ，显式动力学有限元法 ，在碰撞冲击领域已得到广泛的应 = ， (4) 用 ，并显示了在处理大规模接触问题时的优势。不同于静态隐 在时域内应用中心差分法得 式算法 ，显式动力学有限元法无需建立刚度矩阵和求逆运算， 而是采用中心差分法显式求解有限元方程，并通过单点高斯 “ } “一}+ (△ ‘}。卜△ ‘÷) (5) + + 积分和集中质量，提高了求解速度，具有节省计算时间和存储 d叶t= +上 4 (6) 空间以及易于处理问题的优点1]【。 这里假定加速度在整个时间步长内恒定。 1 显式动力学有限元法基本理论 如果已经求得时间步长n时的节点位置和加速度以及时 间步长 一 )时的节点速度 ，则时间步长(n+1)时的位移 d 可 显式动力学有限元基本方程为圈： 假设当前时间步长为第rt．步 ，有如下运动方程 ： 以由式(5)和式(6)解出。因此在整个时域范围内，可由上述积分 递推公式 ，求得各个离散时间点处的位移、速度和加速度。 Mu+ = (1) 中心差分法是一种显式算法 ，但中心差分法是条件稳定 式中， — — 的，即用它求解具体问题时，时间步长必须小于由该问题求解 质量矩阵； c——阻尼矩阵； 方程性质所决定的某个临界值△ ，否则算法将是不稳定的。 K——刚度矩阵； 对显式算法来说，当时间步小于临界时间步时稳定： — — 外界作用力矢量； A￡≤△￡ ： f71 m “ u——时间步长n时的加速度； 式中，甜 为最大 自然角频率。