内爆驱动式超高速发射技术的初步研究

王马法; HIGGINSAndrew J; 焦德志; 黄洁; 柳森

doi:10.11858/gywlxb.20190870

内爆驱动式超高速发射技术的初步研究

doi: 10.11858/gywlxb.20190870

1.
中国空气动力研究与发展中心，四川绵阳　621000
2.
麦吉尔大学工程力学系，魁北克蒙特利尔　H3A2K6，加拿大

基金项目: 国家自然科学基金（11802330）

详细信息

作者简介:
王马法（1986－），男，博士，助理研究员，主要从事超高速发射与毁伤技术研究.E-mail：fujianwmf@163.com

通讯作者:
柳　森（1967－），男，博士，研究员，主要从事超高速碰撞与动能毁伤研究. E-mail：hvi@cardc.cn

中图分类号: O521.3
计量
- 文章访问数: 8918
- HTML全文浏览量: 2782
- PDF下载量: 46
出版历程
- 收稿日期: 2019-12-23
- 修回日期: 2020-01-20

Preliminary Simulation and Experimental Study on Implosion-Driven Hypervelocity Launching Technology

1.
China Aerodynamics Research and Development Center, Mianyang 621000, Sichuan, China
2.
Department of Mechanical Engineering, McGill University, Montreal H3A2K6, Quebec, Canada

摘要

摘要: 为获得10 km/s左右的超高速发射能力，以内爆发射器为研究对象，利用AUTODYN 2D软件对口径为8 mm的内爆发射器进行有限元仿真分析，获得了典型状态下的弹丸发射速度。研制了口径为8 mm的内爆发射器，并在压缩管中填充5 MPa氦气进行实验，分别获得了0.55 g铝合金弹丸7.95 km/s和0.37 g镁合金弹丸10.28 km/s的发射速度，与有限元仿真计算结果的速度偏差分别为15.3%和3.7%。结果表明，设计的内爆发射器具备10 km/s发射能力，满足空间碎片撞击和防护研究的超高速发射需求。
- 发射技术 /
- 超高速 /
- 内爆驱动式
Abstract: For achieving the hypervelocity launching of about 10 km/s, an implosion-driven launcher with the caliber of 8 mm diameter was analyzed using the AUTODYN 2D software. The projectile launching velocities under typical operation condition were obtained. Based on numerical simulation results, several tests of the implosion-driven launcher with the caliber of 8 mm diameter were carried out. The driven gas of helium with the pressure of 5 MPa was filled in the compressed pipe. The experimental results show that the 0.55 g aluminum and 0.37 g magnesium projectiles could be launched to the velocity of 7.95 km/s and 10.28 km/s, respectively, and the relative deviations between the numerical and experimental results are 15.3% and 3.7%, respectively. Consequently, the designed implosion-driven launcher can realize the launching of the projectiles to 10 km/s or even higher which could provide a new ground-test method for investigating the impact features of orbital debris and corresponding shield technologies.
- launching technique /
- hypervelocity /
- implosion-driven

HTML全文

1. 引言

旋耕刀是旋耕机作业最主要的受力部件, 也是最易受损的部件, 其性能直接影响旋耕机的作业效率和耕作质量。随着社会生产力的急剧发展, 大马力旋耕机的使用成为农业机械化发展的必然趋势, 而大耕深旋耕刀的加工制造成为关键技术难题。目前, 大耕深旋耕刀的制造工艺为同比例放大法, 即按比例将原有旋耕刀回转半径扩大至350 mm, 这就带来刀面是否加厚的问题。事实证明^[1], 当刀面厚度不增加时, 机具耕作100 h内的刀具折断率高达30%以上; 当刀面厚度同比例增加时, 机具耕作100 h内的刀具折断率仍在20%左右, 且两种情况下后者的能耗增加了15%。为节省能耗、减少刀具材料使用, 不宜增加刀面厚度, 通过合适的表面处理提高刀具性能是一种有效方法^[2-5]。

激光冲击强化(Laser Shock Peening, LSP)技术是一种新型的表面技术, 目前已广泛应用于材料表面改性的研究。激光冲击强化采用几十纳秒的短脉冲高峰值功率密度(> 10⁹ W/cm²)激光辐射向金属表层, 使金属表面涂覆的保护层吸收激光能量并发生爆炸性气化蒸发, 产生大于1 GPa高压等离子体冲击波, 利用冲击波的力效应使表层材料微观组织发生变化引入残余压应力, 从而提高金属材料抗疲劳、耐磨损和防应力腐蚀等性能^[6-8]。

旋耕刀切土时受到的土壤阻力大多集中在刀具的正切削刃上, 会导致旋耕刀刀柄处产生较大的集中应力, 刀柄在循环交变应力的作用下极易形成疲劳裂纹并失效^[9-12]。为解决这一问题, 可在刀柄应力集中部位引入残余压应力以抵消旋耕刀作业时的外界应力, 抑制或减缓疲劳裂纹的产生, 提高刀具使用寿命。本研究借助ANSYS分析出旋耕刀应力集中区域, 并模拟激光冲击试验得出残余应力引入情况, 同时采用激光冲击强化技术处理应力集中部位, 利用X射线衍射法对残余应力进行测试, 讨论激光冲击强化处理后刀面的残余应力的影响。

2. 理论分析

2.1 旋耕刀受力分析

旋耕机作业时, 旋转刀对土壤进行切削、破碎及抛掷, 土壤对旋耕刀的反作用力构成了土壤阻力^[13]。刀具在切削土壤时, 其耕作深度及切削土壤面均先由小到大, 后由大到小, 其在土壤中的位置也不断变化, 所以刀具受到土壤阻力的大小、方向和作用点在其切削土壤时是不断变化的。在切土过程中, 刀具基本不受轴向力的作用, 可将阻力F沿x、z两个坐标轴分解成F_x、F_z两分力, 则有

$\left\{\begin{array}{l} \boldsymbol{F}=\boldsymbol{F}_{x}+\boldsymbol{F}_{z} \\ F=\sqrt{F_{x}^{2}+F_{z}^{2}} \end{array}\right.$

(1)

旋耕刀的回转半径为350 mm, 耕作深度为21 cm, 刀面切土宽度为6 cm, 土壤阻力计算公式为

$R_{x}=0.7 a b K$

(2)

式中:R_x为土壤阻力, N;a为刀面切土宽度, cm;b为耕作深度, cm;K为土壤比阻, N·cm^-2。一般土壤的K为4~5 N·cm^-2, 粘土的K为6~8 N·cm^-2, 本研究K取7 N·cm^-2, 根据公式(2)可求出旋耕刀所受土壤阻力R_x约为620 N。

2.2 有限元应力仿真

在PRO/E 5.0中建立回转半径为350 mm、刀面厚度为7 mm的旋耕刀三维模型, 并将模型导入ANSYS中进行有限元分析。65Mn旋耕刀的弹性模量E为208 GPa, 泊松比n为0.288。采用Solid 168单元将模型自动划分网格, 共计86 021个节点、53 902个单元。在旋耕刀柄两个侧面及顶面与连接孔处施加约束, 并在受力处施加620 N的集中载荷, 最后运算求解。

图 1(a)为旋耕刀的位移分布等值线图, 工作时旋耕刀的最大位移可以达到0.508 mm, 发生在旋耕刀切削土壤时的正切削刃和侧切削刃处, 而且位移变形是向两边等距递减的, 由于旋耕刀柄固定, 切土时正切削刃和侧切削刃受力, 因而此处变形较大。图 1(b)为旋耕刀的应力分布等值线图, 旋耕刀柄外弯角处应力最为集中, 最大拉应力达到303 MPa, 旋耕刀工作时该处极易折断。

图 1 (a) 位移分布等值线图(b)应力分布等值线图

Figure 1. (a) Isoline of displacement distribution (b) Isoline of stress distribution

下载: 全尺寸图片幻灯片

3. 激光冲击强化有限元仿真

3.1 材料本构方程确定

旋耕刀采用65Mn为材料, 其主要化学成分为(质量分数, %):C 0.65, Si 0.26, Mn 1.10, S 0.002, P 0.019, Cr 0.02, Ni 0.01, Cu 0.01, Fe余量。材料在激光冲击过程中产生动态变形, 因此在选择材料应力-应变本构方程时, 静态方程已经不能反映材料的真实响应, 必须用动态本构方程替代。一般材料的动态本构方程, 目前较多采用Johnson-Cook方程^[14], 其本构方程为

$\sigma=\left(A+B \varepsilon^{n}\right)\left[1+C \ln \left(\frac{\dot{\varepsilon}}{\dot{\varepsilon}_{0}}\right)\right]\left[1-\left(\frac{\theta-\theta_{\mathrm{r}}}{\theta_{\mathrm{m}}-\theta_{\mathrm{r}}}\right)^{m}\right]$

(3)

式中:ε、θ分别为应变和温度为应变率分别为参考应变率和参考温度; θ_m为材料熔点; A、B、n、C、m为待定系数, A、B、n表征材料应变强化项系数, C表征材料应变率强化项系数, m表征材料热软化系数。

3.2 激光冲击波等效载荷计算

段志勇等^[15]在实验基础上, 提出了一个与实际比较相符的半经验型激光冲击波估计模型公式

$p_{\max }=0.8 \sqrt{\rho I_{0}}$

(4)

式中:p_max为激光峰值压力, GPa; ρ为折合密度, g/cm³; I₀为激光功率密度, GW/cm²。折合密度公式为

$\frac{2}{\rho}=\frac{1}{\rho_{1}}+\frac{1}{\rho_{2}}$

(5)

式中:ρ₁为约束层密度, g/cm³; ρ₂为靶材密度, g/cm³。激光功率密度公式为

$I_{0}=\frac{E}{\pi r^{2} \tau}$

(6)

式中:E为激光能量, J;τ为激光脉宽, ns;r为光斑半径, mm。

已知水密度ρ₁=1 g/cm³, 靶材密度ρ₂=7.85 g/cm³, 算得折合密度ρ=1.774 g/cm³; 激光能量E=10 J, 脉冲宽度τ=30 ns, 光斑半径r=3 mm, 算得激光功率密度I₀=1.179 GW/cm², 激光峰值压力p_max=1.16 GPa。激光冲击波作用时间为激光脉宽的3倍左右甚至更高^[16], 实验采用的激光脉宽为30 ns, 取激光冲击波作用时间为100 ns, 其冲击波作用压力与时间的关系如图 2所示, 在50 ns时的压力最大, 最大压力值为1.16 GPa。

图 2 激光冲击波加载曲线

Figure 2. Shock wave of laser loading curve

下载: 全尺寸图片幻灯片

3.3 有限元仿真计算

采用Solid 70单元划分网格并定义材料属性, 设定环境温度为22 ℃, 设定对流及辐射边界条件, 并在应力集中区域的旋耕刀另一侧面设定绝热边界条件。在应力集中部位施加高斯表面热源, 由于区域面积较小, 高斯热源固定施加在待冲击区域中心, 设定载荷步选项和分析选项运算求解, 得到温度场模型。

进行应力分析时, 通过前处理模块定义温度场, 按图 2中激光冲击波加载曲线定义载荷函数, 采用Solid 45单元重新划分网格并定义材料属性, 在应力集中区域定义加载单元组元, 将定义的冲击压力载荷施加在定义的单元组元。同时在单元组元对应的模型各个表面施加透射条件将应力波透射, 防止过大的应力波使模型发生弯扭变形影响仿真结果。最后施加约束并设置求解选项进行运算, 输出表面残余应力结果。

4. 试验材料与方法

激光冲击强化前, 将激光冲击强化区域进行抛光, 利用超声波清洗器将试样在乙醇中清洗15 min, 除去表面油污等杂质, 最后烘干试样^[17]。采用铝箔作为激光冲击的金属涂敷层, 铝箔涂层长为25 mm, 宽为10 mm, 厚为0.1 mm。将铝箔粘贴在旋耕刀应力最集中的刀柄外弯角处, 如图 3(a)所示。使用专用夹具将贴有铝箔的旋耕刀夹于激光冲击强化处理工作台。用水作约束层, 水流厚度为0.5 mm。试验采用江苏大学Gaia-R系列高能量脉冲灯抽运YAG激光器, 激光器调整冲击参数为:激光能量10 J, 波长1.06 μm, 脉冲宽度30 ns, 频率0.1 Hz, 光斑直径6 mm。

图 3 (a) 旋耕刀冲击强化前准备(b)残余应力测点位置

Figure 3. (a) The preparation before LSP (b) Residual stress measuring point location

下载: 全尺寸图片幻灯片

利用超声波清洗器将试样在丙酮溶液中清洗15 min, 并用酒精反复冲洗, 最后吹干^[18-19]。采用MSF-3M型X射线应力测定仪进行残余应力测量, 应力测试取点如图 3(b)所示。X射线分析测试条件为V靶材, 所用衍射晶面为α-Fe(211)衍射晶面, 交相关方法定峰^[20]。

5. 结果分析与讨论

5.1 激光冲击仿真与试验实测对比

从ANSYS后处理输出各测试点的表面残余应力仿真值, 对比试验实测数据, 如图 4所示。各测点的仿真值与实测值曲线呈上下交替分布, 曲线拟合较好, 误差在±20 MPa内。仿真值呈下凹状态, 主要是因为ANSYS中施加的高斯热源为固定热源, 且集中施加在冲击区域的中心, 造成中心温度向边缘扩散, 结果表现为整个区域中心的残余压应力明显高于周边。同理, 实测值呈现下凹状态也是因为激光冲击时, 光斑对准区域中心, 中心受热过大, 从而残余应力的引入明显多于周边。

图 4 激光冲击仿真与实测残余应力

Figure 4. Residual stress datum of simulation and experiment after LSP

下载: 全尺寸图片幻灯片

利用ANSYS模拟激光冲击强化, 其在材料表层引入的残余应力与试验实测的结果吻合得较好, 说明ANSYS可完成激光冲击强化试验, 且误差在允许范围内, 通过此类仿真分析研究, 可对激光冲击强化的工艺参数进行优化, 进一步提高材料性能。

5.2 激光冲击前后残余应力比

图 5(a)、图 5(b)分别是激光冲击前后旋耕刀表面材料残余应力沿x、y轴的分布情况, 未经强化的平均残余压应力仅为146.90 MPa, 强化后的平均残余压应力达到390.70 MPa, 提高了166%。经激光冲击强化后, 最大残余压应力达到412.25 MPa, 在冲击区域的中心位置; 最小残余压应力为375.82 MPa, 在冲击区域的最边缘位置。同时, 冲击后沿x轴及y轴方向的残余应力误差大小在±20 MPa内, 说明整个冲击强化区域的残余应力分布均匀, 无应力集中现象, 说明激光冲击强化能有效提高材料表层残余应力。

图 5 冲击前、后的残余应力

Figure 5. Residual stress before and after LSP

下载: 全尺寸图片幻灯片

由于旋耕刀工作环境恶劣, 且工作时间较长, 易存在应力集中, 引发疲劳裂纹, 激光冲击强化处理引入的残余压应力可消除应力集中的影响, 减慢或抑制裂纹的扩展, 能起到提高板料疲劳寿命的作用, 进而提高旋耕刀使用寿命。

6. 结论

(1) 利用ANSYS对旋耕刀进行应力分析, 得出旋耕刀柄外弯角处应力最为集中, 最大拉应力达到303 MPa, 与日常刀具易受损折断处的位置基本吻合。

(2) ANSYS可模拟激光冲击强化工艺, 其算得的残余应力引入值与试验实测的结果吻合很好, 误差在±20 MPa内, 通过仿真分析可优化相关工艺参数, 进一步提高材料性能。

(3) 经激光冲击强化处理的旋耕刀表层材料残余压应力明显增大, 最大残余压应力达412.25 MPa, 可消除旋耕刀在恶劣工况下应力集中的影响, 减慢或抑制裂纹的扩展, 提高材料疲劳寿命。

图 1 内爆式超高速发射器结构示意图^[11]

Figure 1. Structural diagram of implosion-driven hypervelocity launcher^[11]

下载: 全尺寸图片幻灯片

图 2 内爆式超高速发射器工作原理^[11]

Figure 2. Working diagram of implosion launcher after ignition^[11]

下载: 全尺寸图片幻灯片

图 3 8 mm口径内爆式发射装置工作过程

Figure 3. Launch process of 8 mm caliber implosion launcher

下载: 全尺寸图片幻灯片

图 4 不同充气压力下铝弹丸的速度-时间历史

Figure 4. Velocity-time history of aluminum projectiles with different filling pressure

下载: 全尺寸图片幻灯片

图 5 不同材料弹丸速度-时间历史

Figure 5. Velocity-time history of projectiles with different materials

下载: 全尺寸图片幻灯片

图 6 内爆式超高速发射器实物

Figure 6. Image of implosion-driven hypervelocity launcher

下载: 全尺寸图片幻灯片

图 7 实验装置布局示意图

Figure 7. Layout of experimental equipment

下载: 全尺寸图片幻灯片

图 8 ILT08实验发射器装置回收

Figure 8. Recycle launcher of test ILT08

下载: 全尺寸图片幻灯片

图 9 0.55 g铝合金弹丸速度测试结果

Figure 9. Experimental muzzle velocities of 0.55 g aluminum projectile

下载: 全尺寸图片幻灯片

图 10 0.37 g镁合金弹丸速度测试结果

Figure 10. Experimental muzzle velocities of 0.37 g magnesium projectile

下载: 全尺寸图片幻灯片

图 11 弹丸序列激光阴影成像结果

Figure 11. Projectile photos shot by sequence laser shadowgraph imager

下载: 全尺寸图片幻灯片

表 1 不同弹丸材料和充气压力下仿真计算参数

Table 1. Simulation parameters of launchers with different projectile materials and filling pressure

No.	Material of projectile	Filling pressure/MPa	Diameter of launch tube/mm	Diameter of pump tube/mm	Muzzle velocity/（km·s^–1）
1	Aluminum alloy	4	8	16	8.62
2	Aluminum alloy	5	8	16	9.17
3	Aluminum alloy	6	8	16	9.25
4	Magnesium alloy	5	8	16	10.66

下载: 导出CSV

表 2 内爆发射器实验相关参数

Table 2. Parameters of the implosion-driven launchers in tests

No.	Material of projectile	Mass of projectile/g	Filling pressure/MPa	Diameter of launch tube/mm	Diameter of pump tube/mm	Maximum velocity/（km·s^–1）
ILT04	Aluminum alloy	0.55	5	8	16	7.26
ILT07						6.96
ILT14						7.95
ILT08	Magnesium alloy	0.37	5	8	16	9.73
ILT09						10.28
ILT11						9.36
ILT12						9.77

下载: 导出CSV

参考文献(15)

[1]	HUNEAULT J, LOISEAU J, HIGGINS A J. Coupled lagrangian gasdynamic and structural hydrocode solvers for simulating an implosion-driven hypervelocity launcher [C]//51st AIAA Aerospace Science Meeting Including the New Horizons Forum and Aerospace Exposition. Grapevine, 2013: 1–21.
[2]	王翔, 王青松, 彭建祥, 等. 三级炮超高速发射技术在空间碎片防护研究中的初步应用 [C]//第八届全国空间碎片学术交流会. 北京, 2015: 350–358. WANG X, WANG Q S, PENG J X, et al. Preliminary application of three-stage gas gun hypervelocity launcher techniques in space debris protection research [C]//8th National Symposium on Space Debris. Beijing, 2015: 350–358.
[3]	张旭平, 谭福利, 王桂吉, 等. 基于CQ4的磁驱动10 km/s以上超高速飞片发射 [C]//第八届全国空间碎片学术交流会. 北京, 2015: 385–390. ZHANG X P, TAN F L, WANG G J, et al. Magnetically driven flyer plates to velocities above 10 km/s on CQ4 [C]//8th National Symposium on Space Debris. Beijing, 2015: 385–390.
[4]	文尚刚, 赵锋, 王建, 等. 气炮加载下炸药强爆轰驱动技术的初步实验研究 [J]. 高压物理学报, 2011, 25(1): 36–40. doi: 10.11858/gywlxb.2011.01.006 WEN S G, ZHAO F, WANG J, et al. Primary experimental study on driving technique of strong detonation using gas gun [J]. Chinese Journal of High Pressure Physics, 2011, 25(1): 36–40. doi: 10.11858/gywlxb.2011.01.006
[5]	赵士操, 宋振飞, 姬广富, 等. 一种基于二级轻气炮平台的超高速弹丸发射装置设计 [J]. 高压物理学报, 2011, 25(6): 557–564. doi: 10.11858/gywlxb.2011.06.012 ZHAO S C, SONG Z F, JI G F, et al. A novel design of a hypervelocity launcher based on two-stage gas gun facilities [J]. Chinese Journal of High Pressure Physics, 2011, 25(6): 557–564. doi: 10.11858/gywlxb.2011.06.012
[6]	邢柏阳, 刘荣忠, 郭锐, 等. 强爆轰驱动超高速碎片发射装置设计因素分析 [J]. 国防科技大学学报, 2018(4): 151–158. XING B Y, LIU R Z, GUO R, et al. Analysis on design factors of hypervelocity fragment launcher using strong detonation drive [J]. Journal of National University of Defense Technology, 2018(4): 151–158.
[7]	林俊德, 张向荣, 朱玉荣, 等. 超高速撞击实验的三级压缩气炮技术 [J]. 爆炸与冲击, 2012, 32(5): 483–489. doi: 10.3969/j.issn.1001-1455.2012.05.006 LIN J D, ZHANG X R, ZHU Y R, et al. The technique of three-stage compressed-gas gun for hypervelocity impact [J]. Explosion and Shock Waves, 2012, 32(5): 483–489. doi: 10.3969/j.issn.1001-1455.2012.05.006
[8]	王青松, 王翔, 郝龙, 等. 三级炮超高速发射技术研究进展 [J]. 高压物理学报, 2014, 28(3): 340–345. WANG Q S, WANG X, HAO L, et al. Progress on hypervelocity launcher techniques using a three-stage gun [J]. Chinese Journal of High Pressure Physics, 2014, 28(3): 340–345.
[9]	MOORE J E T. Explosive hypervelocity launchers: PIFR-051 [R]. Physics International Company, 1968.
[10]	WATSON J D. High-velocity explosively driven guns: CR-1533 [R]. Physics International Company, NASA, 1970.
[11]	LOISEAU J, HUNEAULT J, HIGGINS A J. Development of a linear implosion-driven hypervelocity launcher [J]. Procedia Engineering, 2013, 58: 77–87. doi: 10.1016/j.proeng.2013.05.011
[12]	HUNEAULT J, LOISEAU J, HILDEBRAND M, et al. Down-bore velocimetry of an explosively driven light-gas gun [J]. Procedia Engineering, 2015, 103: 230–236. doi: 10.1016/j.proeng.2015.04.031
[13]	HILDEBRAND M, HUNEAULT J, LOISEAU J, et al. Down-bore two-laser heterodyne velocimetry of an implosion-driven hypervelocity launcher [J]. AIP Conference Proceedings, 2017, 1793: 160009.
[14]	田杨萌, 王莹. 炸药爆轰驱动高速击波管发射技术 [J]. 弹箭与制导学报, 2003, 23(3): 221–224. TIAN Y M, WANG Y. A propulsion technology of the fast shock tube driven by high explosive [J]. Journal of Projectiles Rockets Missiles and Guidance, 2003, 23(3): 221–224.
[15]	北京工业学院八系. 爆炸及其作用(上册) [M]. 北京: 国防工业出版社, 1979: 214–219. No.8 Department of Beijing Industrial College. Explosion and its application [M]. Beijing: National Defense Industrial Press, 1979: 214–219.