Deformation field in 316L stainless steel by single shot peening
-
摘要: 运用金属材料表面纳米化试验机对单个弹丸撞击316L不锈钢表面进行了撞击实验；采用激光共聚焦显微镜观察了弹坑的三维形貌，测量不同振动频率下弹坑的直径及离面位移；采用云纹干涉法对弹坑周围的面内应变场进行测量，并分析振动频率及撞击方式对弹坑尺寸、塑性应变大小以及塑性应变区范围的影响；采用有限元方法对单个弹丸垂直撞击试件表面的应变场进行数值模拟，与实验结果进行比较，分析了弹坑周围残余应力的分布。结果表明：随振动频率的增加，弹坑直径和离面位移都增加，频率在50~55Hz，弹坑直径有突变，离面位移和振动频率呈线性关系；振动频率越大，塑性应变越大，塑性应变分布范围均大于弹坑直径的2倍；同一振动频率下弹丸垂直撞击比倾斜撞击的塑性应变大，而塑性应变分布范围相差不大；面内残余应变场的数值模拟结果和实验结果吻合较好，最大误差小于10%。Abstract: The experiment of the single shot impacting the 316L stainless steel surface was carried out using a surface nano-crystallization testing machine. Three-dimensional morphology of the dimple was observed with a laser scanning confocal microscope, and the dimple's diameter and off-plane displacement in different vibration frequencies were also measured. The in-plane strain around the dimple was measured by moiré interferometry. The effect of the vibration frequency and the way of impacting on the dimple size, the plastic strain size and the plastic strain zone were also analyzed. In comparison with the experimental result, the strain field was simulated using the finite element method, and the distribution of the residual stress around the dimple was also analyzed. The result showed that the crater diameter and the off-plane displacement increases with the increase of the vibration frequency. When the frequency is from 50 to 55 Hz, the crater diameter experiences mutations. When the shot impacts the surface vertically, the greater the vibration frequency, the greater the plastic strain and plastic strain zone, and the plastic strain zone are two times larger than the crater diameter. The plastic strain by the vertical impact is slightly greater than the plastic strain by the oblique impact, but it has little effect on the plastic strain zone. The experimental U field strain is in fairly good agreement with the numerical simulation result, and the maximum error is less than 10%.
-
增大装药的能量密度是提高炸药爆炸能量的重要手段。为了提高炸药的能量，在炸药中添加适量的如镁、锂、铝、硼等物质。最开始通过在炸药中添加铝粉，利用铝粉与炸药的爆轰产物反应释放能量以提高炸药的爆炸威力[1]。对含铝炸药的爆轰性能已经有大量研究，冯晓军等[2]通过实验得出铝粉含量对炸药爆炸加速能力的贡献有一最佳值，RDX基含铝炸药最佳含量为15%。周俊祥等[3]研究发现铝粉与爆轰产物的反应明显减缓了水中冲击波后压力的衰减，提高了冲击波的能量及冲量。还有一些实验探究了硼添加剂对炸药性能的影响，如Makhov[4]研究得出硼的加入能显著提高混合炸药的爆热和加速能力。黄亚峰[5]把硼粉加入RDX中，结果显示其可以提高炸药爆热、水下能量等。虽然对在炸药中加入各种添加剂的研究成果丰富，但对含钛炸药的研究很少。而金属钛具有比强度高，无磁性等优点，在地壳中含量约为0.64%, 其含量在金属中排在第四位[6]，其单位质量的钛放热比黑索金等传统炸药高出许多，且钛的密度比较大(钛的密度达到4.54 g/cm3, 而金属炸药中常用的含铝炸药中铝的密度只有2.7 g/cm3)，使得在战斗部狭小的空间内能够携带更多的能量，在军事上具有极其重要的意义。因此，本文中以直径为8 μm的钛纤维作为炸药金属添加剂，来研究钛纤维炸药水下爆炸能量。
1. 实验
RDX基含钛纤维炸药的配方见表 1。钛纤维直径8.0 μm，混药时把钛纤维长度截成3.0 cm，使钛纤维与RDX均匀混合，压成圆柱状药柱，药柱直径1.3 cm，高2.7 cm，如图 1所示。炸药总质量为5.0 g，密度为1.4 g/cm3，钛纤维的质量分数依次增大。采用铝壳包覆炸药，铝壳内径1.3 cm，外径1.6 cm，高5.0 cm。每组配方炸药重复做3次实验。
表 1 含钛纤维炸药配方Table 1. The formulation of Ti fiber explosive序号 φ/% RDX Ti wax 1 95 0 5 2 95 5 0 3 90 10 0 4 85 15 0 水下爆炸实验在水下爆炸箱中进行，如图 2所示。H1=5.0 m，为水的深度；H2=2.5 m，为传感器距离水面的深度；R1=1.0 m，为左侧传感器与药柱的距离；R2=1.2 m，为右侧传感器与药柱的距离；D=5.0 m，为水下爆炸罐的直径。传感器与药柱保持在同一水平面。由于冲击波在水中传播，衰减较慢且信号比较稳定，因此每个距离处只用1个传感器即可。测试前用雷管对测量系统进行现场标定。
2. 结果与分析
2.1 水下爆炸压力峰值分析
图 3所示为含钛纤维炸药水下2.5 m起爆后距爆心1.0 m处的经过低通滤波处理后的冲击波压力时程曲线图。图中冲击波后面的几个压力峰值是由于水下爆炸箱壁面反射形成的，从图中可以看出含钛纤维炸药气泡脉动周期约为51 ms。
图 3 实验得到的钛纤维炸药水下冲击波压力时程曲线Figure 3. History of pressure for underwater shock wave of Ti-fiber explosive by experiment学者们根据实验测试以及理论计算对炸药爆炸压力进行经验拟合，得出各种拟合函数表达式。比较经典的有奥尔连科[7]拟合的函数表达式：
p(t)={pme−tθ0<t<θpm0.368θ/tθ<t<(5∼10)θ (1) 式中：p(t)为冲击波压力随时间变化的函数，pm为压力峰值，θ是衰减时间常数，为pm到pm/e的时间差，t为时间。对压力时程曲线进行拟合前要确定压力峰值与时间常数，在对爆炸压力曲线进行分析的过程中，会发现压力时程曲线的压力峰值若直接读取会有较大误差，且有时不能准确判断峰值所处位置处，因此采用Ethridge[8]作图法，对压力时程曲线的压力取对数，得到压力对数时程曲线，如图 4所示。在压力对数时程曲线中，有一段呈线性，将线性段延伸，延伸段会与曲线上升段交于一点，取交点处对应的压力为该曲线的压力峰值，拟合得到的公式为：
lnp=2.307−0.059t (2) 
图 4 钛纤维炸药水下爆炸压力对数时程曲线Figure 4. History of logarithmic pressure for underwater explosion of Ti-fiber explosive根据式(1)对压力时程曲线进行拟合，如图 5所示。图 5中拟合曲线Ⅰ为t＜θ部分的拟合，拟合曲线Ⅱ为θ＜t＜10θ部分的拟合：对压力时程曲线初始段拟合效果较好，对后面曲线拟合效果不是十分理想，在17.9 μs之后，实验所得曲线在拟合曲线之上，这部分形成原因是钛的加入使得压力曲线相比于RDX有所提升，且在本实验中，药柱有铝壳约束，使钛纤维反应程度增强，加剧了实验压力曲线大于拟合压力曲线这种现象。图 5中拟合曲线Ⅰ的拟合公式为：
p=8.66e−t15.671 (3) 
图 5 钛纤维炸药水下爆炸压力时程拟合曲线Figure 5. Fitted history of pressure for underwater shock wave of Ti-fiber explosive由以上可知：钛纤维炸药的压力时程曲线初始衰减部分可以用式(1)拟合，但后半部分拟合效果不理想。
对压力峰值进行量纲分析[9]：决定水中爆炸冲击波强度的因素有炸药质量Q、装药密度ρe、单位质量炸药释放的化学能Ee、单位质量的钛释放的热量Et、爆炸产物的膨胀指数γe、钛纤维含量的质量分数φ、水的初始密度ρw、状态方程压力参数B、绝热指数γw、药柱离传感器测点处的距离R。水的状态方程可以表达为：
p−pa=B[(ρρw)γw−1] (4) 式中：ρw为水在常温下的密度，B=299.10 MPa，γw=7.15，pa为爆炸前水的压力，p为爆炸时水的压力，ρ为爆炸时水的密度。于是，水下冲击波的峰值压力pm应当为上述因素的函数，即:
pm=f(Q,ρe,Ee,Et,φ,γe,ρw,B,γw,R) (5) 可取Q、ρe、Ee作为基本量，则上式可以化为下面的无量纲关系：
pmρeEe=f(EtEe,φ,γe,ρwρe,BρeEe,γw,R(Q/ρe)13) (6) 由于ρe、Ee、Et、γe、ρw、B、γw这些参数是固定的，无量纲的超压峰值便可简化为:
pmρeEe=f(R(Q/ρe)13,φ) (7) 由于(Q/ρe)13表征炸药的特征长度，则水下爆炸的超压峰值服从几何相似律。利用实验结果计算压力峰值与各物理量之间的关系时，由于公式中的ρe、Ee保持不变，使得其对pm的大小影响保持不变，而且还由于对表达式中的R、Q及pm，都预先规定了所使用的单位，因此式(7)可以简化为:
pm=f(R/Q13,φ) (8) 即在该实验条件下，压力峰值只与R及φ有关，由于实验中在R变量上只有1.0、1.2 m两个数值，无法在R变量上进行深入分析，因此本文中主要对炸药随着φ的变化而分析各个性能的变化。对比冲击波能、比气泡能，总比能量进行量纲分析可以得出相似的结果。
图 6所示为不同钛纤维含量的炸药含钛纤维的质量分数与压力峰值之间的关系, 图中误差棒为每组数据的标准偏差。从图中可看出，在1.0、1.2 m距离处，拟合曲线分别为：
pm=2.002e−φ3.375+6.757 (9) pm=6.553−0.023φ (10) 
图 6 钛纤维炸药压力峰值与钛纤维质量分数关系Figure 6. Relationship between the shock wave peak pressure and the mass fraction of Ti fiber在1.0 m处，压力峰值随钛纤维含量的提高呈指数形式衰减，在1.2 m处，压力峰值衰减较平缓。1.0和1.2 m处，质量分数15%的钛纤维相比RDX峰值降低了22.95%、5.18%。一般情况下，炸药在水中爆炸时，其压力峰值与距离的关系为：
pm=A(RR0)α (11) 式中：pm为冲击波压力峰值，A、α为待拟合系数，R为药柱离传感器的距离，R0为药柱的半径。α的取值一般约为-1，可见压力峰值随比距离的增大，衰减较快。表 2所示为钛纤维炸药在不同含量和不同距离处压力峰值的大小。表 3所示为当加入钛纤维时，不同的钛纤维含量炸药相对于低含量钛纤维炸药压力峰值降低的比例。如钛纤维质量分数5%对应压力峰值降低比例计算过程如下：先得到质量分数0%时压力峰值与质量分数5%时的压力峰值之差，再拿差值除以质量分数0%时的压力峰值即得到压力峰值衰减比例。从表 3中可以看出，当加入钛纤维后，随比距离的增大炸药峰值衰减较慢，原因为钛纤维参与反应，延缓了压力峰值的衰减。且根据表 2中的数值可以计算出RDX随距离的增大压力峰值降低了25%，质量分数15%的钛纤维炸药随距离的增大压力峰值降低了7.70%。
表 2 含钛纤维炸药在不同距离处各质量分数下的压力峰值Table 2. The peak pressure of Ti-fiber explosive at different distances and different mass fractionR/m pm/MPa φ(Ti)=0% φ(Ti)=5% φ(Ti)=10% φ(Ti)=15% 1.0 8.76 7.21 6.91 6.75 1.2 6.57 6.41 6.29 6.23 表 3 含钛纤维炸药在不同距离处各质量分数下的压力峰值降低百分比Table 3. Percentage reduction in peak pressure of Ti-fiber explosive at different distances and different mass fractionR/m (pm(0)-pm(5%))/pm(0) (pm(5%)-pm(10%))/pm(5%) (pm(10%)-pm(15%))/pm(10%) φ(Ti)=5% φ(Ti)=10% φ(Ti)=15% 1.0 17.69% 4.16% 2.32% 1.2 2.44% 1.87% 0.95% 2.2 比冲击波能、比气泡能、总比能量、能量密度分析
比冲击波能Es和比气泡能Eb可以利用测得的冲击波压力时程曲线再经下面的公式计算得到[10]：
Es=4πR2Qρwcw∫6.7θ0Δp2(t)dt (12) Eb=(√1+4CTb−1)38C3k3Q (13) 式中：cw为水的音速，m/s；Tb为炸药气泡脉动时间；C为边界效应的校正系数，C在实验条件保持不变的情况下是常数，在与该实验相同的条件下测得C值为-0.919 5 s-1；k为关于ρw和测点处流体总静水压力的常数，经实验测试k取值为0.002 019。
图 7~8所示为不同钛纤维含量的炸药含钛纤维的质量分数与比冲击波能Es、比气泡能Eb之间的关系。在图 7中，随钛纤维质量分数的增大，比冲击波能呈递减趋势，质量分数为15%的钛纤维炸药在1.0、1.2 m处相比RDX比冲击波能分别降低6.25%、5.04%。在图 8中，随钛纤维质量分数的增大，比气泡能呈增大趋势，质量分数为15%的钛纤维炸药在1.0、1.2 m处相比RDX比气泡能分别增大了24.33%、24.67%。比冲击波能在1.0、1.2 m距离处的拟合曲线分别为：
Es=1.436−0.006φ (14) Es=1.186−0.004φ (15) 
图 7 比冲击波能与钛纤维质量分数关系Figure 7. Relationship between the shock wave energy per unit mass and the mass fraction of Ti fiber
图 8 比气泡能与钛纤维质量分数关系Figure 8. Relationship between the bubble energy per unit mass and the mass fraction of Ti fiber比气泡能在1.0、1.2 m距离处的拟合曲线分别为：
Eb=3.021+0.049φ (16) Eb=3.028+0.049φ (17) 钛纤维作为一种金属添加剂，只有极少部分参与了炸药的爆轰反应，大部分钛纤维是与炸药的爆轰产物反应，释放热量。钛纤维炸药的比气泡能大于RDX的，而比冲击波能小于RDX的，这是由于大部分钛纤维在炸药爆轰过后才参与反应的。
钛纤维炸药在反应的开始阶段，在爆轰波压缩下炸药发生化学反应放热，纯钛的熔点在1 670 ℃，RDX炸药的爆温在3 400 ℃[11]，二氧化钛的熔点约为1 800 ℃[12]。因此，钛纤维表面的氧化膜会因为高温破裂，纯钛外露，周围的炸药与钛发生反应。另一方面，爆轰波经过炸药，相当于对炸药进行绝热压缩，前一部分炸药爆炸产生的爆轰波会破坏钛纤维表面的氧化层，使部分未被氧化的钛释放出来，进而与炸药产生化学反应。由于炸药反应速度很快，这就使得只有很少的钛纤维直接参与炸药的爆炸反应，大部分钛是与炸药的爆轰产物反应。
RDX炸药爆炸产生的气泡脉动时间约在46 ms，在这段时间内，爆轰产物在高温状态下迅速向外膨胀，钛纤维与爆轰产物充分混合，其反应速率大大提高，钛纤维燃烧放出的热量又反过来加速产物的热膨胀，使整个爆炸过程持续下去，表现为气泡脉动时间随着钛纤维质量分数的增大而增加，气泡脉动时间的增加又表现在比气泡能的增大。钛纤维的加入，使得对前期压力做主要贡献的RDX的质量分数减少，大部分钛纤维还没有燃烧放热，而冲击波冲量是对压力的平方积分，从而含钛纤维炸药的比冲击波能随钛纤维质量分数的增大而降低。
炸药水下爆炸产生的总比能量计算公式：
E=Kf(μEs+Eb) (18) 式中：μ=1+1.332 8×10-1pCJ-6.577 5×10-3pCJ2+1.259 4×10-4pCJ3，pCJ=14ρeD2。其中：E为总比能量，MJ/kg；μ为冲击波损失因数；Kf为炸药的形状参数，对于球形，取1.00；非球形，取1.02~1.10；本实验中Kf取1.05。ρe为药柱的密度，kg/m3；D为炸药的爆速，m/s；pCJ为C-J压力，GPa。
图 9所示为钛纤维炸药爆炸总比能量与钛纤维质量分数关系。图中可以看出含钛纤维炸药的总比能量随钛纤维的加入呈增大趋势，加入质量分数为10%的钛纤维后，总比能量在1.0、1.2 m处分别增大了3.72%、5.65%；加入质量分数为15%钛纤维后，在1.0、1.2 m处分别增大了6.79%、9.36%。钛纤维炸药在较远距离处的能量增大比较高，是因为在较远测点处，钛纤维得到充分燃烧，所释放能量更多。钛纤维加入使得炸药比冲击波能降低，比气泡能增大，相当于把炸药的总比能量在水中以气泡能形式释放出来，增强了炸药在水中的破坏效应。
图 9 总比能量与钛纤维质量分数关系Figure 9. Relationship between the total energy per unit mass and the mass fraction of Ti fiber图 9中在1.0、1.2 m距离处，拟合曲线分别为：
E=6.201+0.027φ (19) E=5.675+0.035φ (20) 对RDX而言，因距离的增大，能量降低了8.56%。对质量分数15%的钛纤维炸药而言，因距离的增大，能量降低了6.35%，即含钛纤维炸药在传播的过程中能量衰减较慢。
在本实验中，之所以加入铝壳是因为在前期实验中发现，若不加铝壳则在炸药爆炸时，炸药与水的接触面处反射的稀疏波会使钛纤维炸药中的钛纤维发生破坏，影响整个炸药的反应过程。而加入铝壳后，由于铝壳的约束作用使得钛纤维炸药能在炸药爆炸的瞬间反应，提高钛纤维炸药的能量利用率。铝壳的存在会使炸药的能量一部分用于破坏铝壳，且使铝壳产生一定的速度，从而消耗炸药的能量，使所测得的能量降低。由文献[13]可知：在有外壳约束的情况下，炸药爆炸产生破片的速度为：
Vf=√2Eg(MQ+12)−12 (21) 式中：Eg为壳体破碎形成破片时Gurney能量，M为壳体质量。在本文中铝壳的质量M是相同的，炸药的质量Q也是相同的，破片速度的差别取决于Eg的数值。由文献[14]可知：
√2Eg=0.739+0.622√ρeΦ (22) Φ=0.4889N√q(1−ξ) (23) 式中：Ф为混合炸药的能量示性数，N为每克炸药所产生的气体爆炸产物的物质的量，mol/g；q为每摩尔炸药的爆炸反应热，J/mol；ξ为添加物占炸药的质量比，N、q的值可在文献中查找到。Ф值根据上面公式计算钛纤维炸药的结果为6.47，RDX对应的数值可在文献中直接查找，为6.78，则根据Eg计算公式可以分别计算出钛纤维炸药、RDX炸药的√2Eg值分别为2.63、2.65，其差别仅为0.75%。由于铝壳各种性能是一样的，因此其只破坏铝壳所消耗的能量是相同的，而破碎的铝壳产生速度所消耗的能量差值只有0.75%，且计算过程繁琐，因此本文中做了简化处理，简化过程对最终结果影响不大。
在比气泡能的计算过程中，由文献[15]可知：Eb在实验条件相同的情况下只与炸药的气泡脉动时间Tb有关，而Tb=aW13+bW23，其中a、b在相同的实验条件下为常数，本文实验中a=0.284 89，b=-0.074 63，W为等效的炸药质量，即把含钛炸药按爆热换算成TNT的质量。根据爆热转换，5 g RDX转换成6.91 g TNT，利用Tb的计算公式可求出Tb=0.051 56 s，由前面可知，因铝壳存在使得两种炸药消耗的能量差别为0.75%，则含钛炸药换算成TNT的等效质量应为6.858 g，计算出的Tb=0.051 43 s, 其差别只有0.25%，其数值比较小，忽略其计算对最终结果影响不大。
对于比冲击波能，比冲击波能的计算公式由式(12)可知，Es的变化只与距离和压力峰值有关，当距离不变时只与压力峰值有关。本实验中铝壳对比冲击波能的影响在于铝壳使炸药压力峰值产生变化，而在文献[16]中：有外部约束时(约束外壳厚度较低时，如本文中铝壳厚度为0.15 cm)，水下爆炸测得的压力峰值会稍有增大，即铝壳会影响炸药的比冲击波能，而具体影响范围的大小可以对比前面铝壳对总比能量和比气泡能的影响这两部分进行判断。根据前文中式(18)，可知：Es和Eb的综合变化导致了总比能量E的变化，而E的变化范围只有0.75%，Eb的变化范围只有0.25%，且Es和Eb是同方向变化的，则Es的变化范围小于0.75%。由此可见，铝壳对比冲击波能的影响在实验中也可以忽略。
根据前文中计算得到的炸药总比能量、RDX炸药的密度、钛的密度、钛的质量分数，可以得到炸药的能量密度为：
e=Eaρe+bρTi (24) 式中：a为单位质量钛纤维炸药中RDX炸药的质量，ρe为RDX的密度，g/cm3；b为单位质量钛纤维炸药中钛的质量；ρTi为钛的密度，g/cm3。
图 10所示为钛纤维炸药能量密度与钛纤维质量分数关系。由图 10可知：含钛纤维炸药的能量密度随钛纤维质量分数的增大而增大，在1.0、1.2 m处质量分数15%的钛纤维炸药的能量密度分别增大了19.15%和22.02%。含钛纤维炸药能量密度比纯RDX炸药增大的较多，这是因为钛的密度比RDX大，单位体积所携带的能量比RDX高。
图 10 能量密度与钛纤维质量分数关系Figure 10. Relationship between the energy per unit volume and the mass fraction of Ti fiber图 10中的拟合曲线分别为：
e=8.657+0.108φ (25) e=7.923+0.115φ (26) 对RDX而言，因距离的增大，所测得的能量降低了8.56%。对质量分数15%的钛纤维炸药而言，因距离的增大，所测得的能量降低了6.35%。
从上文中分析可知：1.2 m处钛纤维炸药能量密度相对于RDX增大了22.02%。文献中，在实验条件相同的情况下(1.2 m距离处)10%质量分数的铝纤维炸药总比能量为5.290 MJ/kg[15]，把该比能量按照与钛纤维类似的方法转换成能量密度时，e=5.290.12.7+0.91.6 kJ/cm3=8.82 kJ/cm3, 本文中计算的钛纤维能量密度为9.66 kJ/cm3，比铝纤维能量密度提高了9.52%，这在军用炸药的使用中具有极其重要的意义，在战斗部空间普遍狭窄的情况下，提高单位体积的能量能使战斗部威力得到大幅度的提高。
2.3 对炸药爆炸反应方程式的理论分析
以下根据本实验中炸药所测得的能量，用理论方法确定炸药的爆炸反应方程式[17]。
对于CaHbOcNd类炸药，产物组分多，情况复杂，在这里只考虑常见的10种爆炸产物，因此这里纯药柱爆炸变化方程的一般形式为:
CaHbOcNd=xCO2+yCO+zC+uH2O+wN2+hH2+iO2+jNO+NH3+HCN (27) RDX是一种负氧平衡炸药，这类炸药的含氧量虽不能完全氧化炸药组分中的可燃元素，但足以使其生成完全汽化的产物，不含游离的固体碳，此时产物的主要成分为CO2、CO、H2O、N2、H2，加入钛纤维后，会生成钛的氧化物和很少量的氮化物，这里为简化处理，从能量角度把钛的各种氧化物统一处理成TiO2，忽略钛的氮化物，爆炸变化方程为：
mTi+CaHbOcNd=mTiO2+xCO2+yCO+uH2O+wN2+hH2 (28) 根据RDX炸药在1.0 m处总比能量为6.19 kJ/g，质量分数15%的钛纤维炸药在1.0 m处总比能量为6.61 kJ/g，其差别是由于总质量的15%被换成了钛纤维，而上面忽略了钛发生的复杂反应，按钛纤维放热角度把钛的氧化物统一处理成TiO2，这样得出由于钛纤维的加入，5 g炸药总能量提高了2.1 kJ，而钛纤维反应放出的热量为19.57 kJ/g，RDX释放的能量为5.54 kJ/g，则反应的钛的质量计算过程如下：先计算出被钛纤维替代的总质量为5 g的RDX的15%部分释放的能量为4.155 kJ, 则5 g炸药中参与反应的钛释放的能量为6.255 kJ，则参与反应的钛的质量为0.319 6 g，取值为0.32 g。根据钛的相对分子质量为48，RDX的相对分子质量为222，反应的RDX为4.25 g，则钛与RDX反应的摩尔比为1：2.87，爆炸变化方程变成：
Ti+2.87C3H6O6N6=TiO2+xCO2+yCO+uH2O+wN2+hH2 (29) 产物间发生的平衡反应为：
CO+H2O⇌CO2+H2 (30) 根据质量守恒原则有：
{x+y=8.61u+h=8.612x+y+u=15.22w=8.61 由于式(29)在爆炸化学反应中可以认为能瞬间完成，建立化学反应平衡，先假设该爆炸反应的爆温为4 500 K，查表可得分压平衡常数为100.93，利用迭代法求解这些未知数为：
{x=1.34y=7.27u=5.27w=8.61h=3.34 把RDX的系数化为1，得到如下方程式：
0.35Ti+C3H6O6N6=0.35TiO2+0.47CO2+2.53CO+1.84H2O+3N2+1.16H2 (31) 利用爆炸产物的平均热容计算爆温，先计算爆炸产物的热容，爆炸产物的平均摩尔定容热容与温度的关系一般为：
ˉCi=ai+biT+ciT2+diT3 (32) 式中：Ci为第i种爆炸产物的平均摩尔定容热容，J/(mol·K)；ai、bi、ci、di为与产物组分有关的常数。
计算时一般取前两项，即认为平均摩尔定容热容与温度之间呈直线关系：
ˉCi=ai+biT (33) 对于双原子气体(CO、N2、H2)，ai、bi分别为20.08、0.001 883；对于CO2，ai、bi分别取37.66、0.002 427；对于H2O，ai、bi分别取16.74、0.008 996；对于TiO2，ai、bi分别取56.16、0。
则由式(31)、(33)可得：
ˉC=∑niˉCi=202.53+0.03022T (34) 式中：ni为第i种爆炸产物的物质的量，C为炸药全部爆炸产物在温度间隔T内的平均热容。
假设在爆炸过程中：(1)爆炸过程反应热全部用来加热爆炸产物；(2)爆炸产物处于化学平衡和热力学平衡状态，热容只是温度的函数，则有：
q=ˉC(Te−T0)=ˉCΔT (35) 式中：q为每摩尔炸药的爆炸反应热，J/mol；Te为炸药的爆温，K；T0为炸药的初温，取298 K；ΔT为爆炸产物净增温度。
{ˉCi=ai+biΔTˉC=A+BΔTA=∑niai,B=∑nibi (36) 把式(35)与式(36)联立，可得：
q=AΔT+BΔT2 (37) Te=−A+√A2+4Bq2B+298 K (38) q可通过盖斯定律由各物质的定压生成热计算得出，质量分数为15%的钛纤维炸药爆热计算结果为1 334.18 kJ/mol，即在该计算中取1 mol的RDX反应释放的能量，为1 334.18 kJ。由式(34)可知A、B分别为202.53、0.030 22，将其代入式(38)可求出Te为4 389 K，与假设的4 500 K相差111 K，在计算爆温时能满足要求，式(31)即为质量分数为15%的钛纤维的爆炸反应方程式，类似计算可以求出其他不同钛纤维含量的爆炸反应方程式。
3. 结论
本文中通过对加入不同质量分数钛纤维的RDX炸药做水下爆炸实验，分析压力峰值、比冲击波能、比气泡能、总比能量及能量密度变化情况，得出如下结论：
(1) 钛纤维炸药为非理想炸药，在铝壳约束的情况下，压力时程曲线衰减的初始部分可用指数函数拟合，但压力衰减后半部分拟合效果不理想；
(2) 含钛纤维炸药的压力峰值随钛纤维含量的提高而降低，随距离的增大，压力峰值相对RDX衰减明显变慢，RDX随距离增大压力峰值降低25.00%，质量分数为15%的钛纤维炸药随距离增大压力峰值降低7.70%；
(3) 含钛纤维炸药比冲击波能随钛纤维含量的提高而降低，含钛纤维炸药的比气泡能随钛纤维含量的提高而增大，含钛纤维炸药比冲击波能、比气泡能在不同距离处随钛纤维含量变化趋势基本一致；
(4) 含钛纤维炸药总比能量随钛纤维质量分数的增大而增大，质量分数为10%的钛纤维炸药的能量密度相比质量分数为10%的铝纤维炸药提升了9.52%，钛纤维炸药的能量密度相对于铝纤维炸药的能量密度有较大提升；
(5) 通过实验测得的炸药总比能量，用理论分析方法从能量角度确定了炸药的爆炸反应方程式。
-
图 4 弹丸撞击试件的变形示意图
Figure 4. Schematic diagram of deformation after shot impact substrate
图 5 频率为50 Hz时弹坑的三维形貌图
Figure 5. Three-dimensional topography of crater at frequency of 50 Hz
图 10 不同频率时U场应变分布
Figure 10. Strain distribution in U-displacement at different frequencies
图 11 不同频率下U场和V场的应变分布规律
Figure 11. Strain distribution regularities in U-displacement and V-displacement at different frequencies
图 12 弹丸垂直撞击和倾斜撞击时相同频率下U场云纹图
Figure 12. Contour maps in U-displacement at same frequency under vertical or oblique impact
图 13 相同撞击频率不同撞击角度时应变变化规律
Figure 13. Strain regularities at the same impact frequency and the different impact angles
图 14 x方向应变沿y方向变化规律
Figure 14. Change of strain regularity in x direction along the y direction
-
[1] Marteau J, Bigerelle M, Mazeran P E, et al. Relation between roughness and processing conditions of AISI 316L stainless steel treated by ultrasonic shot peening[J]. Tribology International, 2015, 82:319-329. doi: 10.1016/j.triboint.2014.07.013 [2] Ganesh B K C, Sha W, Ramanaiah N, et al. Effect of shotpeening on sliding wear and tensile behavior of titanium implant alloys[J]. Materials and Design, 2014, 56(4):480-486. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=d2428b4f6fa62767ff66efd3f8f0b4cd [3] Benedetti M, Fontanari V, Santus C, et al. Notch fatigue behavior of shot peened high-strength aluminium alloys: Experiments and predictions using a critical distance method[J]. International Journal of Fatigue, 2010, 32(10):1600-1611. doi: 10.1016/j.ijfatigue.2010.02.012 [4] 栾伟玲, 涂善东.喷丸表面改性技术的研究进展[J].中国机械工程, 2005, 16(15):1405-1049. doi: 10.3321/j.issn:1004-132X.2005.15.023Luan Weiling, Tu Shandong. Recent trends on surface modification technology of shot peening[J]. China Mechanical Engineering, 2005, 16(15):1405-1049. doi: 10.3321/j.issn:1004-132X.2005.15.023 [5] Al-Obaid Y F. Shot peening mechanics: experimental and theoretical analysis[J]. Mechanics of Materials, 1995, 19(2/3):251-260. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ0211961014/ [6] Menig R, Pintschovius L, Schulze V, et al. Depth profiles of macro residual stresses in thin shot peened steel plates determined by X-ray and neutron diffraction[J]. Scripta Materialia, 2001, 45(8):977-983. doi: 10.1016/S1359-6462(01)01063-6 [7] Xing Y M, Lu J. An experimental study of residual stress induced by ultrasonic shot peening[J]. Journal of Materials Processing Technology, 2004, 152(1):56-61. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=f4518f2da0effa1755be7dca1bc56be4 [8] 张洪伟, 张以都, 吴琼.喷丸强化过程及冲击效应的数值模拟[J].金属学报, 2010, 46(1):111-117. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=CAS201303040000302556Zhang Hongwei, Zhang Yidu, Wu Qiong. Numerical simulations of shot-peening process and impact effect[J]. Acta Metallurgica Sinica, 2010, 46(1):111-117. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=CAS201303040000302556 [9] Taehyung K, Hyungyil L, Hong C H, et al. Effects of Rayleigh damping, friction and rate-dependency on 3D residual stress simulation of angled shot peening[J]. Materials and Design, 2013, 46(4):26-37. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=1ab77b975ca1cdf02a33903920abb275 [10] Taehyung K, Hyungyil L, Minsoo K, et al. A 3D FE model for evaluation of peening residual stress under angled multi-shot impacts[J]. Surface & Coatings Technology, 2012, 206(19/20):3981-3988. http://www.sciencedirect.com/science/article/pii/S0257897212002599 [11] Sheng X F, Xia Q X, Cheng X Q, et al. Residual stress field induced by shot peening based on random-shots for 7075 aluminum alloy[J]. Transactions of Nonferrous Metals Society of China, 2012, 22:261-267. doi: 10.1016/S1003-6326(12)61717-8 [12] Mylonas G I, Labeas G. Numerical modeling of shot peening process and corresponding produces: Residual stress, surface roughness and cold work prediction[J]. Surface & Coatings Technology, 2011, 205(19):4480-4494. http://www.sciencedirect.com/science/article/pii/S0257897211002696 [13] Watanabe M, Kishimoto S, Xing Y M, et al. Evaluation of strain field around impacted particles by applying electron Moiré method[J]. Journal of Thermal Spray Technology, 2007, 16(5):940-946. doi: 10.1007/s11666-007-9129-1 [14] Schiffner K, Helling C D. Simulation of residual stresses by shot peening[J]. Computers & Structures, 1999, 72(1/2/3):329-340. http://www.sciencedirect.com/science/article/pii/S0045794999000127 [15] Meguid S A, Shagal G, Stranart J C, et al. Three-dimensional dynamic finite element analysis of shot-peening induced residual stresses[J]. Finite Elements in Analysis and Design, 1999, 31(3):179-191. doi: 10.1016/S0168-874X(98)00057-2 [16] Meo M, Vignjevic R. Finite element analysis of residual stress induced by shot peening process[J]. Advances in Engineering Software, 2003, 34(03):569-575. doi: 10.1016-j.clon.2010.02.005/ [17] Umbrello D, Saoubi R M, Outeiro J C. The influence of Johnson-Cook material constants on finite element simulation of machining of AISI 316L steel[J]. International Journal of Machine Tools & Manufacture, 2007, 47(3/4):462-470. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=da2b63f7ce1c3f04be8bf54f29588b1e [18] Wang J M, Liu F H, Yu F, et al. Shot peening simulation based on SPH method[J]. The International Journal of Advanced Manufacturing Technology, 2011, 56(5/6/7/8):571-578. http://d.old.wanfangdata.com.cn/Periodical/sdgydxxb201006013 期刊类型引用(2)
1. 阚明先，王刚华，刘利新，南小龙，计策，何勇，段书超. 带窗口磁驱动准等熵压缩实验模拟. 强激光与粒子束. 2021(05): 121-125 . 
百度学术2. 阚明先，段书超，王刚华，肖波，赵海龙. 磁驱动飞片发射实验结构系数初步研究. 强激光与粒子束. 2020(08): 117-122 . 百度学术其他类型引用(0)
-
下载:
下载:
百度学术
计量
- 文章访问数: 4464
- HTML全文浏览量: 1716
- PDF下载量: 475
- 被引次数: 2