基于ANSYS／LS—DYNA的护栏冲击模拟分析精度研究.pdfVIP

振 动 与 冲 击 第27卷第4期 JOURNAL OF VIBRATION AND SHOCK 基于ANSYS／LS—DYNA的护栏冲击模拟分析精度研究 张 鹏， 周德源 (同济大学结构工程与防灾研究所，上海 200092) 摘 要：借助ANSYS／LS—DYNA程序从网格密度、壳单元厚度方向的积分点个数、单元公式、接触类型、沙漏控制 等五个方面，对碰撞模拟计算的精度和效率进行了探讨，得出主要结论为：①加密网格、增加壳单元沿厚度方向积分点的 个数可以有效评估和提高计算精度，时间成本以N／l、积分点个数的幂函数方式增加；②B—T公式和H—L公式计算结 果相差不大，B—W—C公式的计算结果用于指导设计是偏于安全的，计算效率按由低到高的顺序排列为H—L、B—T、B — W—C；③单面接触(ASSC)计算精度相对较高，可以作为接触类型计算结果的参考标准；通过合理设定主、从接触面， 采用单向接触类型，可以提高计算效率，并提出了主从接触面的设定原则，结合理论与实践经验给出了相应的解释。 关键词：护栏；冲击；分析精度；计算效率；有限单元法；ANSYS；LS—DYNA 中图分类号：0242．21。U417．1+2 文献标识码：A 目前汽车一护栏碰撞安全性能的研究主要有两种 几种评估手段着重探讨以下两方面的问题： 方法：实车碰撞试验和数值模拟¨l2 J。尽管实车碰撞试 1)在数值分析的范围内，评估和提高计算精度的 验被公认为综合评价汽车碰撞安全性能的最基本、最 方法； 有效的方法，然而实车碰撞试验费用高昂，由于碰撞过 2)模拟精度和计算效率的平衡问题，即如何在确 程涉及的因素(碰撞速度、碰撞角度等)较多，为了对碰 保模拟精度的条件下尽可能提高计算效率。 撞过程进行全面研究，需要进行多次试验，研究往往要 计算所用计算机主要性能指标：①CPU—Intel 受到经济条件的限制。相比之下，数值模拟在时间和 Pentium 4，主频3．0GHz；②内存1GB。 经济方面的优势就非常明显。 1 有限元模型描述 碰撞是一个瞬态的大位移和大变形过程，涉及几 何、材料、接触状态等多重非线性，大型复杂结构碰撞 本文建立的3跨波形梁护栏一冲击块有限元模型 过程的模拟计算往往需要耗费很长的机时，如何确保 如图1所示，跨度为4m。护栏波形梁和立柱的几何尺 模拟精度(特别是无试验数据比较的情况下的模拟精 寸参照规范 规定，其截面厚度分别为3mm、5mm，详 度)和提高计算效率成为一个值得研究的问题。LS— 细尺寸见图2，冲击块 、Y、 向的几何尺寸为0．4m× DYNA程序中，提供了大量的可控变量(如接触类型、 0．2m×0．4m，质量约为2．5t，与一小型轿车质量相当。 单元公式)，ANSYS的参数化语言APDL可以方便地对 波形梁、立柱和冲击块材料均采用Q235结构钢，其主 这些变量进行设定，这些都为工程技术人员评估和提 要力学性能参数见表 1，真实应力应变曲线如图3 ‘ 高计算精度与计算效率提供了便利条件。 所示。 鉴于已有的有限元理论，有限元计算结果精度可 表1 材料力学性能参数 以通过以下四种方法进行评估： 1)不同的网格形式和不同网格密度的计算结果 比较 ， ； 2)同类单元采用不同单元计算公式的计算结果 比较； 碰撞过程材料应变变化速率较大，这将对弹塑性 3)非线性问题梁、壳单元厚度方向取不同积分点