侵彻数值仿真中网格密度对单元失效应变影响分析

Analysis on influence of grid density on element failure strain in penetration numerical simulation

  • 摘要:
      目的  旨在解决数值模拟平头弹对金属板的侵彻问题中,网格尺寸对单元失效应变取值和弹体剩余速度的影响。
      方法  应用有限元软件LS-DYNA对Q235钢材料单轴拉伸实验进行数值仿真,结合实验中试样断裂时的伸长量得到该网格密度下的单元失效应变,绘制出失效应变随网格密度变化的曲线并进行动态修正;然后,进行平头弹侵彻Q235钢板的数值仿真,并对靶板进行不同尺寸的网格划分。其中Q235钢材料的单元失效应变按修正曲线取值,将弹体的剩余速度与实验结果对比,分析数值仿真中网格密度对金属板抗侵彻问题仿真结果的影响。
      结果  研究表明,数值仿真中选取的单元失效应变应随网格密度的增大而增加,在金属板的抗侵彻问题中随着网格密度的增大,对弹体剩余速度预报的仿真结果逐渐收敛并与实验结果相近。当网格尺寸为0.5 mm时,研究速度段内数值仿真与实验拟合曲线的平均相对误差最小,为5.13%;仿真与实验的误差在低速段更大,且低速度段弹体剩余速度对网格密度更加敏感。
      结论  相关计算方法及研究结果对弹体侵彻问题中网格密度和材料失效应变的取值及变化规律具有一定参考价值。

     

    Abstract:
      Objectives  In order to solve the problem of numerical simulation of flat-nosed projectile penetration into metal plates, the influence of mesh size on element failure strain value and residual velocity of projectiles was studied.
      Methods  The finite element software LS-DYNA was used to simulate the process of uniaxial tensile test of Q235 steel sample, and the failure strain of the element under the grid density is obtained by the elongation of the tensile sample during fracture. In the meantime, the correction curve of the failure strain with the grid density was plotted and dynamically corrected. Then, the numerical simulation of flat-nosed projectile penetrating Q235 steel plate was carried out with the target plate meshed with different sizes. The failure strain of Q235 steel material was selected according to the correction curve. Finally, the residual velocity of the projectile is compared with the experimental results to analyze the influence of mesh size on the simulation results of the penetration resistance problem of the metal plate in the numerical simulation.
      Results  The results show that the element failure strain selected in the numerical simulation should increase with the increase in grid density, and in the case of the metal plate anti-penetration problem, it should increase with the increase of the grid density. In the problem of penetration resistance of metal plates, the simulation results of residual velocity prediction gradually converge with the experimental results. With the increase of mesh density, when the grid size is 0.5 mm, the average relative error of the numerical simulation and the test fitting curve in the velocity section is 5.13%, and the error between the numerical simulation and the test is larger in the low-speed section. Moreover, the residual velocity of the projectile body is more sensitive to mesh density in the low velocity range.
      Conclusions  The related calculation methods and research results have a certain reference value for the selection of mesh size and material failure strain in the projectile penetration problem.

     

/

返回文章
返回