一种用于动态拉伸装置新型试件装卡方式.docVIP

文章编号： 一种用于动态拉伸试验装置的新型试件装卡方式* 朱耀，庞宝君，盖秉政，管公顺 （哈尔滨工业大学航天学院，黑龙江 哈尔滨 150001） 摘要：针对杆杆型动态拉伸试验系统，设计了一种新型的楔楔 * 收稿日期：200X-XX-XX；修订日期：200X-XX-XX 通讯作者：朱耀（1983—），男，博士研究生，研究方向为材料动态力学性能，zhuyao@hit.edu.cn 。 (a) Shape of the fastener and connection (b) Meshing of the specimen and the fastener 图2 试件与杆的连接 Fig.2 Connection of specimen and the bars 试件与杆的连接如图2所示。通过在入射杆与透射杆的一端开一个楔v = 6m/s v = 20m/s 图4 直接动态拉伸波形 Fig.4 The immediate tensile waveform 直接式杆杆拉伸系统的几何参数为子弹长l0为200mm，杆1长l1为1000mm，杆2长l2为1000mm，应变片1和应变片2的贴片位置都在两杆的中间，x1及x2均为400mm。模拟得到的波形如图4所示，其中图4(a)为子弹速度为6m/s时的波形，透射波信号放大了6倍，图4(b)为子弹速度为20m/s时的波形，透射波信号放大了10倍。对比图4(a)和图4(b)可以发现，子弹速度越大，透射波的振荡就越大。 1.3 材料模型及其他有限元参数 在本文的数值模拟中，试件材料模型采用双线性各向同性硬化模型，具体形式为 （1） 其中σy为动态流动应力，σ0为屈服应力， Ep为切线模量，εP为塑性应变。 子弹、杆1和杆2的材料模型为线弹性。具体材料参数见表1和表2。 表1 试件材料参数 Table 1 Material parameters of specimen 材料 E / GPa μ ρ / kg/m3 σ0 / MPa Ep / MPa 失效应变 2024Al 71 0.3 2780 370 1000 0.75 表2 杆的材料参数 Table 2 Material parameters of bars 材料 E / GPa μ ρ / kg/m3 45号钢 210 0.3 7800 本文采用三维拉格朗日六面体solid164实体单元，通过ANSYS/LS-DYNA建立拉伸系统的1/4平面对称三维模型。子弹、杆1、杆2的径向单元尺寸为1mm，轴向单元尺寸为5mm，在卡口附近单元轴向尺寸为1mm，试件实验段单元尺寸为0.25mm，试件及卡口处的网格划分见图3。所有网格全部采用Mapped方式划分，对于1/4圆柱体来说，这种方式比Sweep方式划分的网格质量更好。整个直接式杆杆动态拉伸系统的1/4平面对称三维模型共包含45200个节点与36796个六面体单元。所有接触都采用自动面面接触，设置动摩擦与静摩擦系数均为0.1，接触刚度（界面的法函数因子）取为0.5，其余都是默认设置。计算时间为0.4ms，时间步长程序自动选择。每个算例计算时间大约为1个小时。 2 波形计算结果 与经典SHB实验计算公式一样，在动态拉伸实验中可以利用如下三个方程计算试件的拉伸应力-应变曲线[11]。 (2) (3) (4) 其中，s是试件中的平均应变率，εs是试件中的平均应变，σs为试件中的平均应力；εr和εt 分别为反射波和透射波对应的应变，由贴在入射杆和透射杆上的应变片可以直接读出；C0和A0分别是杆的声速与横截面积；L0和As分别是试件的初始长度与横截面积。 根据数值计算得到的应变波形信号，运用方程(2)-(3)可以得到利用模拟拉杆计算得到的材料动态应力应变关系。同时根据输入的材料参数（见表1和表2），利用方程(1)可以得到输入模型的应力应变曲线。 (a) v = 6m/s (b) v = 10m/s (c) v = 20m/s (d) result of 2D simulation 图5 模拟结果与输入模型的应力应变曲线 Fig.5 Stress-strain curves of simulation and input model 图8(a)-(c)是根据方程（2）-（4）利用模