A method of spherical indentation experiment based on the split Hopkinson pressure bar system
-
摘要: 提出一种基于霍普金森压杆装置的动态球形压痕实验方法,通过将硬质合金小球置于2个试件中间,实现加载过程中2个试件的同时压入,以得到准确的压痕力与位移的关系。利用有限元软件ABAQUS/Explicit对该实验方法进行了数值模拟,从实验结果的可行性、准确性等方面,对新实验方法与传统的动态压痕实验方法进行了比较;采用新方法对铝合金材料进行了实验,并得到了压入力-位移曲线等关系。研究结果表明:采用新实验方法能得到较准确的结果,能较真实地描述压痕过程。Abstract: A new technique for testing dynamic spherical indentation of materials was proposed.It was implemented by putting a WC alloy ball between two specimens, one of which was bonded with the input bar and the other was bonded with the output bar.Then the two speciemens were pressed in the process of indentation.With this simple arrangement we can accomplish a very accurate measurment of force and displacement during the indentation.A finite element simulation was carried out by using the ABAQUS/Explicit software to evaluate this new experimental technique.The results show that the new method gives a more accurate result than other dynamic indent experimental methods.In addition, by testing the alumimum alloy materials, the curves of the indentation force and depth were gotten based on the new technique.The investigation diaplays that the technique proposed can porvide more accuate results and refect the indentation process more truely.
-
Key words:
- solid mechanics /
- spherical indentation /
- SHPB /
- indentation experiment
-
利用压痕法测试材料的力学性能已有悠久的历史, 与其它的材料性能测试方法相比, 压痕法具有无损(微损)、便捷等优点[1-3]。传统的硬度检测就是压痕法的一种应用——通过分析压痕的最终几何量及荷载, 给出材料的硬度指标。硬度按照不同的测试方法又可分为布氏、维氏、洛氏、里氏、肖氏等硬度[4]。随着实验技术的进步, 已容易通过仪器化的压痕实验获取实验的压入荷载与位移(p-d)曲线, 如何从p-d曲线中获得更多的材料力学参数成为研究者们关注的问题。W.C.Oliver等[5]提出了由静态球形压痕力与位移的关系确定金属材料弹性模量的方法。此后, Y.P.Cao等[6]也进行了该方面的工作; Y.T.Cheng等[7-8]利用量纲分析方法对压痕实验进行了归纳, 并结合有限元数值模拟进行分析, 总结出压痕测试的基本原则。Y.P.Cao等[9]提出了一个分析框架, 从球形静压痕实验的p-d曲线进行反演, 获得幂强化材料的弹性模量、屈服应力及硬化指数等参数。
随着压痕法测试技术的发展, 研究者开始将其用于材料动态力学性能的研究。20世纪90年代霍普金森压杆(SHPB)被用于动态压痕测试。G.Subhash等[10]将SHPB入射杆作为加载部件, 杆与试件接触端套上锥形压头, 分别采用力、位移传感器测试压入力和位移, 获得相应的p-d曲线。M.Nilsson[11]采用一维应力波理论通过入射杆与透射杆上的信号来计算压头的压入位移。张新等[12]在上述装置上进行了铜、钛、铝合金动压痕实验研究。G.Subhash的方法[10]中采用悬臂梁测量位移时程曲线, 由于悬臂梁的振动模态比较复杂, 各个模态被激发的程度与加载速率有关, 应变测量值和端点的位移关系是不确定的, 所以无法的到准确的测量结果。M.Nilsson[11]采用了传统的实验方法[10], 利用入射、反射波信号来计算压头的位移, 但是锥形段的存在以及压头的安装导致其有效性存在疑问。
本文中提出一种采用双试件的球形压痕测试方法; 采用有限元软件ABAQUS/Explicit对该实验方法进行数值模拟, 以评估方法的有效性与准确性; 采用新测试方法对7075铝合金进行动态压痕实验, 研究7075铝合金的动态力学性能。
1. 动态压痕实验
1.1 实验装置
图 1为双试件动压痕实验装置示意图。将硬质合金小球放在2个试件中间(如图 1(b)所示), 以此取代传统霍普金森压杆实验装置中的试件。在入射应力脉冲作用下, 硬质合金小球同时压入2个试件, 小球对其两侧试件有相同的压入力, 因此能够方便地利用传统霍普金森压杆实验方法测出硬质合金小球压入试件过程中的压入位移和压入力。通过入射杆和透射杆的输出信号可得到材料的p-d曲线。
1.2 实验数据处理
将测量得到的入射波信号εi(t)、反射波信号εr(t)和透射波信号εt(t), 代入式(1), 即可获得小球对2个试件的压入速度、加速度和压力与时间的关系:
{˙u(t)=12c0[εi(t)−εr(t)−εt(t)]u(t)=12c0∫t0[εi(t)−εr(t)−εt(t)]dtp(t)=A0E0εt(t) (1) 式中:c0为应力波在杆件中的传播速度, E0为压杆的杨氏模量, A0为杆件的横截面面积, u(t)和
(t)为每一时刻硬质合金小球压入试件的深度与速度, p(t)为试件受到的压入力。
2. 动态压痕实验的数值模拟
2.1 计算模型
利用有限元软件ABAQUS/Explicit中的轴对称模型, 对M.Nilsson的方法[11]和本文中提出的新方法进行数值模拟, 并将模拟结果与实验结果进行了比对。计算中, 试件与杆件之间保持波阻抗匹配:ρscsAs=ρ0c0A0, 其中ρs、cs、和As分别为试件的密度、弹性波速和横截面积。子弹为尺寸为Ø14.5 mm×400 mm, 入射杆与透射杆尺寸均为Ø14.5 mm×1 000 mm。钢试件匹配的杆件材料为60Si2Mn, 铝合金试件匹配的杆件材料为铝合金。试件尺寸为Ø14.5 mm×6 mm, 试件材料模型选用Johnson-Cook(J-C)本构模型。硬质合金小球的材料为碳化钨, 直径为1.6 mm。材料具体参数见表 1~2, 其中:E为弹性模量, μ为泊松比, ρ为密度, A、B、C、n和m为试件的J-C本构常数。
表 1 杆件及试件的弹性参数Table 1. Elastic constants of the bars and the specimens材料 E/GPa μ ρ/(g·m-3) 60Si2Mn 206 0.29 7.80 45钢 210 0.29 7.80 铝合金 71 0.28 2.70 碳化钨 450 0.28 18.00 表 2 试件的J-C本构参数Table 2. J-C constants of specimens材料 A/MPa B/MPa C n m 45钢 1 150 739 0.014 0.26 1.03 铝合金 369 684 0.008 3 0.73 1.70 图 2为2种试件的加载方式示意图。在入射杆和透射杆中布置波形测量点, 在试件两端设置位移测量点。对于M.Nilsson的方法[11], 采取2个不同角度的压头:锥头(45°半锥角如图 2(a)所示)和平头(90°半锥角)。
2.2 实验的模拟分析
图 3为根据M.Nilsson的实验方法[11]模拟得到的压入深度时程曲线, 从图中可以看出, 压入深度的计算值与测量值存在较大差异。造成误差的原因主要是应力波在锥形段中的反射透射, 会影响最终计算结果。在实验中, 小球需要焊接在锥头上, 锥头与杆连接通常采用粘接的方式[12], 这使得局部的变形状况变得复杂, 也会导致计算结果与实验结果有差异。图 4为根据本文中实验方法模拟得到的压入深度时程曲线。由于排除了锥头与杆之间的波阻抗不匹配带来的影响, 提高了实验的准确性, 因此能够得到与位移测量点输出结果较一致的压头压入位移时程曲线。图 4(a)为铝合金计算结果, 与测量点测到的结果比较一致; 图 4(b)为45钢计算结果, 虽然与测量点的测试结果存在一定的误差, 但与真实压入过程吻合很好。
图 5为铝合金试件的Misses应力云图。由图 5可以看出, 在加载过程中, 入射杆端试件在应力脉冲作用下首先与硬质合金小球发生碰撞, 获得压痕, 然后压力脉冲再通过小球, 与透射杆端试件发生作用, 最终透射杆端试件得到压痕。图 5显示在撞击开始阶段, 前后试件的压痕存在较小的差异, 碰撞发生10 μs之后, Misses应力基本上达到对称, 两试件得到相同的压痕深度。
3. 铝合金动态力学性能的实验测量
采用双试件方法对7075铝合金试件进行动态压痕实验, 实验在Ø14.5 mm霍普金森压杆装置上进行, 选用规格为Ø14.5 mm×400 mm的子弹进行打击, 得到3.5、8和12 m/s等3种打击速度下的实验结果。图 6为实验测得的典型波形(电压信号), 由于透射波信号非常微弱, 透射杆上采用了半导体应变片。
图 7为波形处理得到的相关曲线, 其中应变率强化作用明显。由图 5也可观察到, 压痕范围内应变分布非常复杂, 在凹痕顶部应变由内向外依次变小(但未到零), 应变率难以给出。G.Subhash等[10]将压痕的压入速度与压痕几何尺寸的比值定义为平均应变率, 其中压痕几何尺寸为压痕深度或对角线的长度, 则上述实验的平均应变率分别约为3.5×104、4.0×104和4.2×104 s-1, 一般SHPB实验很难达到。
4. 结论
提出了双试件动态压痕实验方法, 通过数值模拟对其与现有实验方法进行比较, 结果表明采用双试件法能够得到较现有方法更准确的压头压入位移-时间曲线, 且双试件法的理想实验条件更易满足, 实验可靠性更高。在对铝合金材料进行实际实验后, 得到4×104 s-1左右的平均应变率估计值。该方法可在传统SHPB上轻易实现, 同时数据处理简单, 能够用于104~105 s-1应变率范围内材料动态性能的研究。
-
表 1 杆件及试件的弹性参数
Table 1. Elastic constants of the bars and the specimens
材料 E/GPa μ ρ/(g·m-3) 60Si2Mn 206 0.29 7.80 45钢 210 0.29 7.80 铝合金 71 0.28 2.70 碳化钨 450 0.28 18.00 表 2 试件的J-C本构参数
Table 2. J-C constants of specimens
材料 A/MPa B/MPa C n m 45钢 1 150 739 0.014 0.26 1.03 铝合金 369 684 0.008 3 0.73 1.70 -
[1] Ahn J H, Kwon D. Derivation of plastic stress strain relationship from ball indentions: Examination of strain definition and pileup effect[J]. Journal of Materials Research, 2001, 16(11): 3170-3171. [2] Hernot X, Bartier O, Bekouche Y, et al. Influence of penetration depth and mechanical properties on contact radius determination for spherical indentation[J]. International Journal of Solids and Structures, 2006, 43(14/15): 4236-4153. [3] Lee H, Lee J H, Pharr G M. A numerical approach to spherical indentation techniques for material property evaluation[J]. Journal of the Mechanics and Physics of Solids, 2005, 53(9): 2046-2063. [4] Taber D. A sample theory of static and dynamic hardness[J]. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Science, 1948, 192(1029): 247-274. [5] Oliver W C, Pharr G M. An improved technique for determining hardness and elastic modulus using load and displacement sensing indentation experiments[J]. Journal of Materials Research, 1992, 7(6): 1564-1583. [6] Cao Y P, Qian X Q, Lu J. On the determination of reduced Young's modulus and hardness of elasticplastic materials using a single sharp indenter[J]. Jounral of Materials Research, 2006, 21(1): 215-224. [7] Cheng Y T, Cheng C M. Relationships between hardness, elastic modulus, and the work of indentation[J]. Applied Physics Letters, 1998, 73(15): 614-616. [8] Cheng Y T, Cheng C M. Scaling, dimensional analysis and indentation measurements[J]. Materials Science and Engineering: R: Reports, 2004, 44(4): 91-149. [9] Cao Yan-ping, Lu Jian. A new method to extract the plastic properties of metal materials from a instrumented spherical indentation loading curve[J]. Act Materialia, 2004, 52(13): 4023-4032. [10] Subhash G, Koeppel B J, Chandra A. Dynamic indentation hardness and rate sensitivity in metals[J]. Journal of Engineering Materials and Technology, 1999, 121(3): 257-263. [11] Nilsson M. Dynamic hardness testing using a Split-Hopkinson pressure bar apparatus[R]. Sweden: Swedish Deference Research Agency, 2002. [12] 张新, 侯兵, 李玉龙.基于霍普金森压杆系统的动态压痕实验[J].爆炸与冲击, 2011, 31(3): 256-262.Zhang xin, Hou Bing, Li Yu-long. Dynamic indentation experiment based on the Split-Hopkinson pressure bar system[J]. Explosion and Shock Waves, 2011, 31(3): 256-262. 期刊类型引用(1)
1. 韩国庆,谈梦婷,邓佳杰,李逸,毛小建,张先锋. AlON透明陶瓷动态压痕加载响应研究. 力学学报. 2025(02): 473-487 . 百度学术
其他类型引用(8)
-