SHPB入射杆形状对动态断裂参量影响的数值方法研究.pdfVIP

结构强度研究 SHPB入射杆形状对动态断裂参量影响的 数值方法研究 陈华燕曾祥国钱晓军王清远许书生 (四川大学建筑与环境学院，成都) 擒要；利用分离式Hopkinson压杆(SmB)测量三点弯曲试样材料的动态应力强度因子，应力波载荷 的惯性效应和弥散效应对试验结果不可忽略的影响已被许多学者的实验结果证实．将不同SHPB入射杆 形状产生的载荷一时间历程曲线作为动态有限元模型的载荷信息，采用ANSYS／Ls．DmA有限元软件， 基于虚拟裂纹闭合技术(VCCT)计算得到三点弯曲试样的动态应力强度因子变化曲线，曲线上起裂 时间对应的值即为材料的动态断裂韧度。进一步考虑入射杆形状及惯性效应对冲击载荷曲线的影响， 研究了动态应力强度因子的变化规律，为Hopkinson压杆测试动态断裂韧度提供了合理，可行的方法． 关键词：Hopkinson压杆；三点弯曲试样；动态应力强度因子；虚拟裂纹闭合技术 1 引言 近年来材料的动态断裂参数得到了越来越多材料和力学工作者的重视【l-2]。动态断裂由于断裂 过程具有时间短，应变率高等特点例，不能不考虑应力波传播和惯性效应的影响，一方面使得动 态断裂的理论研究存在极大的困难，另一方面，动态断裂韧度的测试也非常困难，尚未形成统一 的测试规范，其测试难点主要体现在：1)动态载荷的施加和测试；2)动态应力强度因子的实时 测试；3)起裂时问的确定。目前，动态载荷的施加和测试技术主要有落锤、摆锤冲击加载和 根据一维应力波理论得到完整准确的加载波以及确定起裂时间，得到越来越多的应用14J．动态应 力强度因子的实时测试，可采用试验方法，例如应变片监测；也可采用近似计算方法[5-61；或采用 数值计算方法【7J。本文采用数值方法，利用有限元软件建立试样的动态模型嘲，并运用虚拟裂纹 闭合技术[9-14J首先计算出能量释放率响应曲线，再根据能量释放率与应力强度因子的关系得到应 力强度因子变化曲线，曲线上起裂时间对应的值即为材料的动态断裂韧度。随着计算机硬件和软 件快速且广泛的发展，该方法在进行动态断裂分析方面正成为一种经济可行的方法，可以对各种 新奇材料组成的、几何特征复杂和各种不同加载历程下的实际结构进行分析。 2基于虚拟裂纹闭合技术的有限元算法及验证 图1显示3-维虚拟裂纹有限元法计算能量释放率的两个步骤。图1(a)和图l@分别为裂纹扩 展前和裂纹扩展后裂纹附近的节点和单元示意图。裂尖节点采用双重节点技术。由图1(a)可得裂 尖和裂尖之前共有三对节点：节点e、f为裂尖节点；节点g、h和节点i、j分别重合，表示裂纹尖端 为e(0节点，裂纹未向前扩展。图l嘞假设裂纹向前扩展一个微量，则此时i0)节点成为裂纹尖端， 节点e、f和节点g、h的约柬释放。计算模型(a)。可以得到裂尖附近的节点力{e，嘭，E，圪}；计 ·528· 结构强度研究 MPa·√磊，与传统奇异单元法结果相对误差小于o．3％，与试验 vccT方法计算得到的群。_26．53 值的相对误差小于1．2％，说Ij,0VCCT方法具有很高的精度。 3动态断裂算例 文献【6】采用HopkiIls∞压杆试验装置测试动态断裂韧度，入射杆加载端加工成同冲击摆锤相 似的楔形头．由于几何形状的改变必然影响应力波的传播，从而得到了不同的载荷时间曲线。文 献【6】比较了三种楔形头形状对动态断裂韧度测试造成的影响，楔形头长分别为3、15．5、40mm(图 5)。 图5三种楔形头尺寸 图6三种楔形头典型荷载时间曲线 图7三种楔形头置．(，)曲线 动态断裂三点弯曲试样与静态断裂韧度试样具有相同尺寸(图2)，选取平面应变3节点三角 将图6的载荷曲线作为有限元计算的载荷输入，采用虚拟裂纹闭合技术计算得到如图7所示 墨(f)曲线，进行比较可以得到：(1)三楔形头的荷载曲线对应于时间轴有明显的振荡，其中 随楔形头尺寸越大，荷载曲线振荡越厉害。通过VCCT方法计算所得动态应力强度因子墨(f)响 应曲线按体呈上升趋势，这与文献【6】中采用弹簧质量模型计算所得结果一致。墨(f)曲线随载荷 的振荡有微小波动，文献【5】、【9】等所得结果有相似现象。且楔形头尺寸越大，载荷曲线振荡越为 明显。(2)从0至5工心期