Movement characteristics of high-velocity fragments in water medium: Numerical simulation using LS-DYNA
-
摘要: 基于大型有限元分析软件ANSYS/LS-DYNA,建立三维长方体高速破片在水介质中运动的有限元动力分析模型,采用ALE方法对破片在水下运动过程进行流固耦合数值模拟,获得了破片的速度衰减曲线。研究了速度衰减规律、破片墩粗变形规律以及冲击波传播过程。得到高速破片的侵彻能力随速度的变化规律:当初速度大于910~1 115m/s时破片头部将产生显著变形,并大大影响其侵彻阻力;当破片速度较小时,水中侵彻距离随破片初速的增大而增大,当破片速度达到某临界值以后,侵彻距离将随初始速度的增大而逐渐减小。Abstract: Based on the large-scale finite element analysis software ANSYS/LS-DYNA, the dynamic analysis model was established for the movement of a three-dimensional rectangular high-velocity fragment in water medium.The multi-material arbitrary Lagrangian Eulerian method was used to simulate the underwater movement process of the high-velocity fragment.And the velocity attenuation curve of the fragment was given.The velocity attenuation, the mushrooming rule and the shock wave propagation process were studied.The penetration performance of the high-velocity fragment was obtained by considering its velocity change.When the initial velocity is 910-1 115m/s, the mushrooming is definite, which will increase the penetration resistance for the fragment sufficiently.When the initial velocity is low, the penetration distance of the underwater fragment will increase with the increasing of its initial velocity.When the initial velocity arrives at some critical value, the penetration distance will decrease as the initial velocity of the fragment increases.
-
自1949年Kolsky建立分离式霍普金森压杆(split Hopkinson pressure bar, SHPB)实验技术以来,SHPB实验技术在此后的近百年里获得了长足发展,为了满足工程和科学研究的需要,也不断演变出其他动态加载(拉伸、扭转、剪切、围压等)和诊断技术,但SHPB始终是成熟度最高、使用最广泛的实验技术[1-10]。SHPB实验技术建立在3个基本假设基础上:(1)一维应力假设;(2)应力均匀性假设;(3)压杆弹性假设[10-11]。实验过程中需对压杆材质,压杆与样品的几何尺寸、表面状态、摆放姿态(同轴度)、端面摩擦效应消除和修正等提出较高要求,才能够满足以上假设,进而保证实验数据的准确性和可靠性[10-12]。如果实验设计和实施过程中出现失误,无法满足以上假设,将会导致实验数据出现较大偏差,甚至得出错误的研究结论[2, 10-11]。在以上众多影响SHPB实验结果准确性和可靠性的因素中,样品与压杆端面摩擦效应对材料塑性流动应力的影响是普遍关注的一个问题[2, 11, 13]。
金属样品外表面机加状态良好,粗糙度较低时,二硫化钼(MoS2)等润滑剂能够有效降低样品端面摩擦力,但在某些特殊工况下,样品外表面无法达到较高精度[14-16]。例如:铍作为聚变堆中的重要结构材料,普遍认为机械加工会在其外表面引入亚毫米级的机加损伤区(张鹏程等[15-16]认为该机加损伤区厚度可达0.5 mm),且常规退火无法完全消除机加损伤区内的损伤微结构和残余应力,因此,实验前采用腐蚀液化学腐蚀掉机加铍样品外表面的损伤区是常用的消除机加损伤区影响的方法[14-19]。然而,化学腐蚀都是从活性较高的晶界等缺陷位置开始,样品外表面的腐蚀速率无法统一,这必然导致化学腐蚀后样品尺寸精度降低,尤其是外表面粗糙度增加等问题出现[20-22]。腐蚀后样品表面粗糙度增加会提高样品端面摩擦力,因此,能否有效降低表面粗糙度较高样品的端面摩擦力,直接决定着SHPB实验结果的准确性和置信度。
本文中,采用SHPB实验技术,研究3种表面粗糙度紫铜样品在MoS2充分润滑和完全不润滑条件下的本构性能,以期获得样品表面粗糙度以及润滑效果对端面摩擦力和最终应力-应变曲线的影响。
1. 实 验
为了尽可能避免样品尺寸、加载状态浮动、应变率效应等因素对样品表面粗糙度和润滑效果对端面摩擦力的影响,SHPB实验中采用直径10 mm的马氏体钢压杆,压杆端面和外表面粗糙度(Ra)均为0.8 μm,加载弹速控制在(17±0.5)m/s,样品与压杆的同轴度优于0.1 mm(安装工艺保证),样品选用应变率效应不敏感的紫铜,样品尺寸统一为
∅ 4 mm×4 mm,如图1所示,机加端面粗糙度分别为0.8和1.6 μm。将端面粗糙度为1.6 μm的紫铜样品放入稀硝酸中腐蚀8 h后,获得端面粗糙度为3.2 μm的腐蚀样品。采用扫描电子显微镜(scanning electron microscope, SEM)和原子力探针[23-25]观测表面粗糙度为0.8、1.6 和3.2 μm样品的实验端面。显微分析结果证实,3种状态样品的表面粗糙度存在显著差异,如图2所示,图中上方为SEM图像,下方为原子力探针获得的表面图像。SEM显微分析结果表明,机加粗糙度为0.8 μm的样品端面存在非常均匀的圆形车刀痕迹,进刀量约为5 μm,样品中心无凸台;机加粗糙度为1.6 μm的样品端面的圆形车刀痕迹非常不均匀,进刀量约为10 μm,且样品中心存在凸台;酸蚀样品端面存在大量腐蚀坑,已无法观测到车刀痕迹。原子力探针显微分析的观测区域均为100 μm×100 μm,观测区域在半径约为1 mm的圆周上。将原子力显微观测结果导入Nanoscope Analysis后处理软件,通过软件自动分析获得3种状态样品观测区域的表面粗糙度实测值分别为73、534(该样品观测区域的部分测量点起伏超出测量仪器的最大量程4 μm,实际值应略大)、840 nm(该样品观测区域表面起伏太大,大量数据点超出了4 μm的仪器量程,实际值应更大)。样品端面粗糙度实测值显著低于机加要求的表面粗糙度数值,推测可能是由于观测区域太小以及部分测量点超出仪器量程导致的。但实测粗糙度的相对数值与要求的0.8、1.6和3.2 μm较接近,因此,依旧用0.8、1.6和3.2 μm表征3种样品的端面粗糙度。
2. 实验结果
采用MoS2作为润滑剂,对3种表面粗糙度的紫铜样品分别在完全润滑和完全不润滑(干摩擦)两种状态条件下开展SHPB动态压缩实验。实验结果如图3所示,图中Cu4×4-0.8-RH-1~Cu4×4-0.8-RH-3表示紫铜样品机加尺寸为
∅ 4 mm×4 mm,端面粗糙度为0.8 μm,在MoS2完全润滑的条件下开展SHPB实验的第1~3发实验;Cu4×4-0.8-WRH-1~Cu4×4-0.8-WRH-3表示紫铜样品机加尺寸为∅ 4 mm×4 mm,端面机加粗糙度为0.8 μm,在完全无润滑的条件下开展SHPB实验的第1~3发实验,依此类推。实验样品几何尺寸、弹速和回收试样形状汇总如表1所示。表 1 样品尺寸和弹速实测值以及回收样品形状Table 1. Measured values of sample size and velocity, and the shapes of recovered samples样品 直径/mm 长度/mm 弹速/(m∙s−1) 回收样品形状 Cu4×4-0.8-RH-1 4.003 4.014 17.07 标准圆柱 Cu4×4-0.8-RH-2 4.000 4.012 17.11 Cu4×4-0.8-RH-3 4.003 4.007 17.04 Cu4×4-0.8-WRH-1 4.004 4.010 17.03 鼔形 Cu4×4-0.8-WRH-2 4.000 4.014 17.00 Cu4×4-0.8-WRH-3 4.003 4.016 17.01 Cu4×4-1.6-RH-1 4.002 4.016 17.03 略带鼓形 Cu4×4-1.6-RH-2 3.998 4.022 17.02 Cu4×4-1.6-RH-3 3.998 4.018 17.04 Cu4×4-1.6-WRH-1 4.004 4.013 17.13 鼔形 Cu4×4-1.6-WRH-2 3.987 4.006 17.18 Cu4×4-1.6-WRH-3 3.990 4.004 17.13 Cu4×4-3.2-RH-1 3.975 3.982 17.00 显著鼔形 Cu4×4-3.2-RH-2 3.990 4.005 16.99 Cu4×4-3.2-RH-3 3.974 3.990 16.99 Cu4×4-3.2-WRH-1 3.986 4.002 17.15 显著鼔形 Cu4×4-3.2-WRH-2 3.980 3.998 17.32 Cu4×4-3.2-WRH-3 3.980 4.004 17.04 当样品端面机加粗糙度为0.8 μm时,完全润滑样品的应力-应变曲线实验结果重复性良好,无润滑样品应力-应变曲线的一致性也较好,但屈服强度和整体塑性流动应力水平更高,如图3(a)所示。这主要是由于端面无润滑时引入的较大端面摩擦力提高了样品的塑性流动应力水平,与此同时,样品端面摩擦力提高了样品的应力状态复杂程度,这也导致无润滑样品应力-应变曲线分散性略有增加,并使得回收样品为鼓形,如表1所示。当样品端面机加粗糙度为1.6 μm时,完全润滑样品的应力-应变曲线实验结果重复性依旧较好(回收样品略带鼓形),且优于无润滑状态样品(回收样品为鼓形),如图3(b)和表1所示。但二者的应力-应变曲线在塑性流动段几乎重合,且与端面粗糙度为0.8 μm的无润滑样品的塑性流动段也几乎重合,如图3(d)所示(汇总所有实验数据,并根据GJB3756-99中的B类评估方法,置信因子取1.732)。这表明MoS2已经无法有效降低端面粗糙度为1.6 μm的样品的端面摩擦力。当样品端面粗糙度为3.2 μm时,润滑和无润滑样品的应力-应变曲线整体几乎重合,难以区分,且各自应力-应变曲线的一致性均较差(回收样品的鼓形也更显著),如图3(c)~(d)和表1所示。图3(c)中的数据分散性大可能存在以下原因:(1)酸蚀后由较大表面粗糙度引入的端面摩擦力进一步强化了样品的应力状态复杂程度,导致样品的应力-应变曲线分散性也进一步增大;(2)由于酸蚀是由点到面的腐蚀过程,导致腐蚀后样品轮廓尺寸精度较低(几乎所有腐蚀样品直径在两个相互垂直方向上的实测值均相差0.02 mm左右,表明腐蚀后的样品已不是标准圆柱体),导致测量得到的样品尺寸值(二次测量平均值)与实际值存在较大误差,进而增加了实验数据的分散性。图3(d)中酸蚀后样品应力-应变曲线的塑性流动应力水平也整体较低,这可能是由于酸蚀在样品外壁面引入的腐蚀缺陷充当了预置损伤的角色,降低了材料的整体强度[22, 26-28]。
3. 分析与讨论
由以上3种不同表面粗糙度紫铜样品在完全润滑和完全无润滑两种工况条件下的常温SHPB实验结果可以看出,样品表面粗糙度越大,引入的端面摩擦力越大,数据的分散性也越大,置信度越低,如图3所示。表面粗糙度为0.8和1.6 μm的样品在充分润滑条件下的实验结果虽然存在一定偏差,但整体差异低于5%,这在不过分追求实验数据置信度的工况条件下是可以接受的。与此同时,部分材料(尤其存在织构和非均匀性较大的材料)自身性能的分散性可能大于5%[11, 29-31],且常规金属样品外表面很容易达到1.6 μm的机加粗糙度,这导致关于样品外表面粗糙度微小差异对SHPB实验结果影响的量化分析研究一直较少[11]。但部分特殊工况条件下,样品外表面难以达到1.6 μm的粗糙度,例如传统机加会在铍等材料的外壁面引入较大的机加损伤区,从而影响实验数据的准确性,因此实验前通常采用腐蚀液腐蚀掉一定厚度的样品外壁面来避免机加损伤区对材料性能的影响,这就导致最终样品的表面粗糙度较大[14-16]。基于以上不同粗糙度紫铜的SHPB实验结果可知,腐蚀会显著增加样品的表面粗糙度,进而提高样品的端面摩擦力,且MoS2并不能有效消减端面粗糙度大于1.6 μm的样品的端面摩擦力。采用SHPB实验技术研究表面粗糙度较大样品的本构关系时,需采取更有效的润滑措施消减样品端面的摩擦力,或在数据处理过程中扣除端面摩擦力的影响,对原始数据做必要的修正,才能够获得置信度较高的实验数据[11, 13, 29]。
相同加工工艺条件引入的机加损伤深度基本固定[15]。因此,研究铍等存在较强机加损伤行为的材料的本构关系时,可通过提高样品尺寸的方式降低机加损伤区域的体积占比,从而减小机加损伤区域对材料整体力学性能的影响,机加损伤区的体积占比为
1−(1−2δ2/D)2(1−2δ1/H) ,样品的初始设计尺寸(D 和 H )越大,机加损伤区的体积占比就越小,如图4所示。假定机加损伤区厚度恒定为0.5 mm,则机加损伤区∅ 2 mm×2 mm、∅ 4 mm×4 mm和∅ 6 mm×6 mm样品的体积占比分别为87.5%、57.8%和42.1%。但如果想量化评估机加损伤区对材料本构关系的影响,还需进一步评估样品尺寸效应、加载状态(应变率效应)等因素的影响,相关工作将在后续研究中作更详细的介绍。4. 结 论
通过机械加工和腐蚀液腐蚀的方式制备了3种不同典型表面粗糙度的紫铜样品,并基于SHPB实验技术在MoS2充分润滑和完全无润滑两种工况下,研究了样品粗糙度和端面润滑对紫铜SHPB动态压缩本构实验结果的影响,得出以下结论。
(1)样品表面粗糙度越大,引入的端面摩擦力越大,样品的受力状态越偏离一维应力状态,导致动态压缩应力-应变曲线数据的分散性越大。
(2)MoS2能够有效消减表面粗糙度为0.8 μm的样品的端面摩擦力,但对表面粗糙度大于1.6 μm的样品端面摩擦力的消减效果极为有限。
(3)开展表面粗糙度较大(>1.6 μm)样品的SHPB实验研究时,需采用比MoS2更有效的润滑措施消减端面摩擦力,或在数据处理时扣除端面摩擦力的影响,对实验数据进行修正才能够获得准确的实验数据。
(4)采用腐蚀液腐蚀材料外表面的方法虽然能够有效消除铍等材料外壁面的机加损伤区,但会显著降低样品尺寸精度,尤其是会提高样品端面粗糙度,进而提高端面摩擦力,降低实验数据的准确性和置信度。通过提高铍样品的几何尺寸,降低机加损伤区的体积占比,从而降低机加损伤区的影响程度或许是获取铍材料高精度动态本构关系参数的一个可行方法。
-
-
[1] 沈晓乐, 朱锡, 侯海量, 等.高速破片侵彻防护液舱试验研究[J].中国舰船研究, 2011, 6(3): 12-15. http://d.wanfangdata.com.cn/Periodical/zgjcyj201103003Shen Xiao-le, Zhu Xi, Hou Hai-liang, et al. Experimental study on penetration properties of high velocity fragment into safety liquid cabin[J]. Chinese Journal of Ship Research, 2011, 6(3): 12-15. http://d.wanfangdata.com.cn/Periodical/zgjcyj201103003 [2] 白金泽. LS-DYNA3D理论基础与实例分析[M].北京: 科学出版社, 2005. [3] 吴建萍, 魏志芳, 李柱.基于LS-DYNA的球形破片侵彻明胶的数值模拟[J].机电技术, 2011, 1: 28-30. http://www.cnki.com.cn/Article/CJFDTotal-JDJS201101011.htmWu Jian-ping, Wei Zhi-fang, Li Zhu. Numerical simulation of spherical fragment penetrating into gelatin based on ANSYS/LS-DYNA[J]. Mechanical & Electrical Technology, 2011, 1: 28-30. http://www.cnki.com.cn/Article/CJFDTotal-JDJS201101011.htm [4] Trevino T. Applications of arbitrary lagrangian eulerian(ALE)analysis approach to under water and air explosion problems[R]. AD-A384983, 2000. [5] 王肖钧, 赵新.高速钢球在水中贯穿过程的数值计算[J].爆炸与冲击, 1992, 12(3): 213-218. http://www.cnki.com.cn/Article/CJFDTotal-BZCJ199203003.htmWang Xiao-jun, Zhao Xin. Numerical study of a steel ball penetrating in water with high speed[J]. Explosion and Shock Waves, 1992, 12(3): 213-218. http://www.cnki.com.cn/Article/CJFDTotal-BZCJ199203003.htm [6] Tate A. A theory for the deceleration of long rods after impact[J]. Journal of the Mechanics and Physics of Solids, 1967, 15(6): 387-399. http://www.sciencedirect.com/science/article/pii/0022509667900105 期刊类型引用(0)
其他类型引用(2)
-