Characteristics of wall pressure and cavitation on the plate subjected to underwater explosion shockwaves at any angle of incidence
-
摘要: 以平面波理论为基础,推导了无限平板在全入射角度下的冲击波壁压载荷计算公式,利用实验数据对壁压公式进行修正,提出了一种适用于计算有限尺度平板壁压的经验公式;分析了不同入射角度下壁压载荷的变化特性,初步研究了壁压载荷负压特性对平板局部空化的影响。结果表明:修正后的壁压曲线与实际壁压曲线吻合较好;入射角度的增大会加快壁压衰减过程,并使最低壁压的绝对值减小;随着药量或平板厚度的增加,壁压最低负压的绝对值增大,形成局部空化的能力增强;局部空化仅在一定条件范围内才会形成,空化范围受局部空化形成压力及冲击强度等因素的影响较大。Abstract: Wall pressure as a load parameter should be confirmed exactly, and then offered for investigating dynamic response of structure such as a plate subjected to underwater explosion.On the basis of plane shock wave hypothesis, a theoretical formula was deduced and proposed to calculate the wall pressure values of an infinite plate subjected to underwater explosion shockwaves at any angle of incidence.By using the wall pressure data derived from underwater explosion experiments, the theoretical formula of wall pressure was revised to be suitable for plates with finite dimensions.With the revised formula, the characteristics of the wall pressure on the plate were analyzed by considering various incidence angles, and the effects of the negative wall pressure values on local cavitation on the plate were discussed primarily.The results show that the wall pressure values derived from the revised formula agree well with those from the experiments.As the incidence angle increases, the wall pressure decays more rapidly, and the absolute value of the minimum wall pressure decreases.With the more explosive or the thicker plate, the absolute value of the minimum wall pressure is larger, and local cavitation takes place more likely.Local cavitation can only come into being in certain conditions, and the cavitation range will be affected by the factors such as cavitation pressure criterion, shock intensity, and so on.
-
Key words:
- mechanics of explosion /
- plate wall-pressure /
- underwater explosion /
- shockwaves /
- local cavitation /
- oblique incidence
-
漂珠空心颗粒(cenospheres, CPs)材料不仅能够作为填充相使材料轻质化,又能在单体破碎后释放塌缩空间,作为缓冲材料时具有良好的可压性。Wang等[1]利用石英砂堆积体取代传统的缓冲材料对弹体在高g值冲击时进行缓冲可发现石英砂堆积体对加速度的降低具有明显的效果。为研究脆性空心颗粒材料堆积体在缓冲防护中的应用,所以以典型的脆性空心颗粒材料粉煤灰漂珠为研究对象考察其在冲击载荷下的力学行为及缓冲特性。传统的复合泡沫材料[2-3]是将脆性空心颗粒材料填充到高聚物中制成,复合泡沫材料具有明显的应变率效应,为考察填充材料对复合泡沫应变率效应的影响机制。需要对脆性空心颗粒的应变率效应进行分析。所以对脆性空心颗粒堆积体进行静动态压缩实验研究。然而,目前对松散颗粒的研究大多局限于实心颗粒材料,空心颗粒材料的力学性能及应变率效应的研究却鲜有涉及。
松散脆性颗粒堆积材料的应变率效应问题虽然经历了较长时间的研究,但目前有关机理的认识尚存在分歧。且目前研究多集中于干燥砂等一类实心颗粒材料,Bragov等[4]利用Kolsk实验对干燥砂进行研究表明,在较大的应变率范围内材料力学行为并无明显变化。Song等[5]通过SHPB实验考察干燥砂在450~1 400 s−1应变率范围内的动态力学性能,发现其对应变率并不敏感。但是Farr[6]和Yamamuro等[7]观察到10−2~103 s−1应变率范围内干燥砂的应力应变曲线有显著区别,具有典型的应变率效应。然而,对于空心颗粒材料而言,颗粒破碎后材料的体积可压缩性使得问题更为复杂,目前尚无文献对脆性空心颗粒材料在不同应变率下的力学性能进行的研究。
对于实心颗粒的应变率效应,Huang等[8-9]从颗粒破碎这一微细观角度研究了干燥砂颗粒材料的动态力学行为和应变率效应,发现颗粒材料宏观应变率效应机理的本质是颗粒破碎特征对冲击速率的依赖性,实验结果表明同应力加载水平下干燥沙准静态加载时的颗粒破碎量和破碎程度明显高于动态加载,且颗粒破碎所需能量减少,即产生相同的颗粒破碎量准静态加载时所需的能量更少,因此动静态加载下颗粒破碎能量利用效率不同导致了宏观上的应变率效应。另外,Huang等[10]研究指出对于均匀级配的颗粒材料,颗粒尺寸小的材料的屈服应力比颗粒尺寸大的材料髙,相同应力水平下颗粒破碎量更低,而级配良好的颗粒材料颗粒破碎量比级配均匀的颗粒材料的低。
本文以具有不同粒径分布的粉煤灰漂珠为研究对象,利用筒压实验考察空心颗粒在低速冲击下的力学性能。考虑到空心颗粒的体积可压缩性和存在应力平台阶段,本文采用等位移加载的方式,通过加载套筒限制颗粒材料的单次加载压缩应变为50%;然后对不同加载速率下所得的破碎颗粒进行粒径分析,基于颗粒破碎率和破碎程度分析粉煤灰漂珠的材料宏观应变率效应;最后通过Hardin破碎势理论分析粉煤灰漂珠的能量吸收与颗粒破碎之间的关系。
1. 实验材料
图1为实验所用的大粒径粉煤灰漂珠(large cenospheres, LGs)和小粒径粉煤灰漂珠(little cenospheres, LTs)的宏观图像,图2为LGs和LTs的扫描电镜(scanning electron microscope, SEM)图片,并通过碎片(如图3所示)可看出两者是具有一定厚度的空心颗粒。本文使用的颗粒材料为同一种粉煤灰漂珠,LGs的粒径为800~1 000 μm壁厚为(50±10) μm,LTs的粒径为500~600 μm壁厚为(30±10) μm。总体来说,空心微球的壁厚大约为球体外部半径的0.1倍,内部球体空腔大约占球体体积的72%,由SEM图可知漂珠的球壁也是多孔的,球壁材料孔隙率约为10%。文献[11]中对粉煤灰漂珠进行研究发现与文中所用的漂珠具有相同的性质。LGs的颗粒散体的堆积密度为0.33 g/cm3,LTs的颗粒散体的堆积密度为0.35 g/cm3。
在实验中选择在限制应变为50%下的颗粒材料堆积体的强度进行分析。进行准静态实验时采用最大加载值为100 kN的Suns万能试验机,采用刚性套筒对颗粒材料的堆积体进行约束,此时应变率为0.001 s−1。对漂珠材料进行中应变率实验时采用Suns全自动落锤试验机,此时应变率为150和300 s−1。在准静态和落锤冲击时采用尼龙加载装置,试样直径为 48.5 mm,高度为20 mm。预留加载的位移为10 mm,如图4(a)所示。并通过 SHPB (split Hopkinson pressure bar)装置对粉煤灰漂珠在高应变率下的力学性能进行测试,此时的应变率约为800 s−1,采用DHPB (direct-impact Hopkinson pressure bar)实验方法即质量块直接撞击试样,通过背部透射杆的信号获得试样的压缩载荷时间曲线,通过高速摄影对试样的位移进行统计,试样直径为35 mm,高度为10 mm,预留高度为5 mm,如图4(b)所示。在实验中将每次的实验进行3次重复试验求出平均值对实验的结果进行分析。
为了利用颗粒破碎率和破碎程度对颗粒堆积体的应变率效应进行解释,所以在实验中对颗粒材料堆积体在应变为50%时的破碎率和粒径进行统计分析。通过利用振荡筛分法对粒径进行粗略的筛分对颗粒的破碎率进行统计,然后通过激光粒度仪进行分析得出颗粒材料破碎后的粒径级配曲线。
2. 实验结果
图5为两种漂珠颗粒堆积体在不同应变率加载下的工程应力应变曲线,由图可知,该颗粒材料在准静态和动态冲击下力学性能有显著区别,但是在当前动态加载应变率范围内,材料动态力学行为变化不大。在低速冲击下,两种漂珠颗粒堆积体均表现出应变率效应,相比于准静态加载,材料的动态强度和平台段的应力均有显著提高。
考虑到粉煤灰漂珠的主要吸能途径是颗粒破碎。所以需要通过测量破碎颗粒的质量分数来表征空心颗粒的破碎率,然后检测碎片的大小分布。将压缩实验后的漂珠颗粒进行回收,通过振筛机将破碎后的试样按照粒径分布分为若干份,然后统计不同粒径分布内颗粒碎片的质量分数,结果如图6所示。以图6(a)中LTs试样颗粒破碎规律为例,试样初始粒径分布范围为500~600 μm,因此认为粒径500 μm以上部分为未破碎颗粒,颗粒在200~500 μm范围内分布较少,表明单个漂珠颗粒对径劈裂为两部分的破碎模式发生概率较小。颗粒破碎后碎片粒径主要集中于200 μm以内,表明大部分颗粒破碎为4部分或更多碎片。另外,随着应变率增大,压缩后完整颗粒的质量分数减小,但破碎程度较高的碎片质量分数越大,如图6(a)所示。结果表明,随着冲击速度的增大,漂珠颗粒的总破碎量以及单个颗粒的破碎程度均有所提高,即漂珠颗粒的破碎机理出现明显的冲击速率敏感性。
图6(b)为LGs试样的粒径分布统计图,其基本规律与LTs的实验结果一致。对比可知,颗粒材料在动态加载下每种粒径分布范围内质量分数分布情况并无明显差异,300和800 s−1应变率下颗粒材料的碎片分布规律基本一致,同样材料的强度等宏观力学性能也比较接近,说明漂珠颗粒材料的力学行为与颗粒的破碎机理密切相关。为精确统计漂珠颗粒碎片的分布规律,将各组分的颗粒碎片通过激光粒度仪进行分析,结合分样筛所得结果可获得颗粒材料在破碎后的各个粒径范围内的碎片的分布曲线,如图7所示。
由颗粒破碎后不同粒径级配曲线可以看出在0.001和300 s−1时颗粒的破碎率越来越高,颗粒的破碎程度越来越大,得出的规律与用分样筛得出的规律相同。在300和800 s−1应变率下的颗粒碎片的粒径分布曲线基本相同,表明颗粒的破碎行为在低速冲击范围内表现出显著的冲击速率敏感性。
3. 讨论与分析
3.1 力学性能
Huang等[10]研究指出对于级配均匀的颗粒材料,平均颗粒尺寸较小的材料屈服应力比平均颗粒尺寸大的材料髙。将LTs和LGs在相同应变率(0.001 s−1)下的力学性能进行对比,如图8所示,可发现在相同的应变率下LTs的力学性能较LGs的力学性能好。在相同应变率下小颗粒较大颗粒的强度提高35%~40%。由于小颗粒的堆积体密度较大颗粒堆积体的密度较高,小颗粒堆积体孔隙率较大颗粒的孔隙率低,小颗粒堆积体更加密实这是造成小颗粒堆积体强度增高的重要原因。
参考泡沫材料力学性能的分析方法,取应变为5%时对应的应力作为颗粒堆积体的抗压强度,将LGs和LTs在不同应变率下的抗压强度统计如图9所示,可以发现在150和300 s−1范围内大小粒径漂珠的强度分别增强39%和51.5%,在300和800 s−1内大小粒径漂珠的强度均未发生明显变化,在0.001和150 s−1范围内大小粒径的颗粒材料的强度分别增强200%和195%,表明漂珠颗粒堆积体具有显著的应变率效应。
3.2 颗粒破碎率分析
本实验控制材料的压缩应变均为50%,通过对等压缩量时、不同应变率下颗粒破碎规律进行统计,分析材料宏观的应变率效应。首先,对颗粒的破碎量进行统计分析,如图10所示,可以看出随着应变率的提高颗粒破碎量会提高。文献[12]提出颗粒尺寸越大级配越差,颗粒越容易产生破碎,因为颗粒尺寸越大初始缺陷越多,强度也越低;颗粒级配越差,则单个颗粒拥有的接触点数越少,颗粒越容易形成对径受集中力压缩的状态而产生破碎[13]。通过对相同应变率下不同粒径的颗粒破碎量进行对比可见,在相同应变率下大颗粒的破碎率较小颗粒的破碎率高。
倪素环等[14]指出能量吸收是决定破碎比的主要因素。通过式(1)得出不同应变下的能量吸收,根据实验得出随着应变率的升高颗粒材料所吸收的能量增高,在相同应变率下小颗粒的吸能效果较大颗粒的的吸能效果提高35%~48%。如图11所示。
W=∫εd0σdε (1) 式中:σ为应力,ε为应变,εd为最终应变。
分析可知,颗粒材料产生相同压缩量时,LTs比LGs的能量吸收高,动态下颗粒破碎量多于静态,因此其动态耗能高于准静态。然而,由于颗粒材料的压缩体积是相同的,因此同等能量耗散时,动态冲击下的应力必然高于准静态,即材料的强度和流动应力均表现出应变率敏感性,而动态冲击下颗粒破碎量的增多是粉煤灰漂珠应变率效应产生的主要原因。
通过对应变为50%时颗粒的碎片进行统计分析,得出高应变率下能量吸收显著增强的原因有两个:(1) 单个空心颗粒破碎程度的提高;(2) 空心颗粒损伤百分比的增加。
3.3 Hardin破碎势分析
为了量化颗粒破碎与能量耗散之间的关系,提出了破碎率、粒度和破碎势等指标。Hardin[15]提出利用颗粒破碎前后的粒径分布曲线和D=0.074 mm的界限围成的面积来对颗粒破碎进行衡量,围成的面积称为破碎势。Hardin[15]认为小于0.074 mm的颗粒很难再产生破碎,形成更小粒径的颗粒,因此将D=0.074 mm作为颗粒破碎的界限,用粒径分布曲线和该界限围成的面积作为破碎势Bp。颗粒破碎前后的粒径分布曲线围成的面积定义为整体破碎势Bt。实际上Bt就是破碎前的破碎势与破碎之后破碎势之差。Hardin提出了相对破碎势的概念
Br=Bt/Bp0 ,其中Bp0指的是初始破碎势。图12为破碎势示意图。池昌江[16]对准脆性颗粒在静态压缩下的颗粒破碎通过Hardin破碎势对颗粒的破碎进行分析,得出相对破碎量与能量吸收之间的关系。黄俊宇[17]也通过相对破碎势对石英砂在不同应变率下的颗粒破碎进行分析。利用Hardin提出的破碎势和相对破碎势的概念进行分析。图13为LTs和LGs在实验前后的粒径分布曲线。
利用Hardin相对破碎势概念可以计算得到LTs和LGs在不同应变率实验前后的相对破碎势,并且可以利用位移-载荷曲线求出每次实验所输入的能量,得出单位能量下的相对破碎势。如表1所示。
表 1 粉煤灰漂珠破碎势Table 1. Cenospheres breaking potential颗粒大小 应变率/s−1 Bp Bt Br 输入总能量/J 单位能量相对破碎势/J−1 LGs 0.001 0.826 0.430 0.520 4.610 0.113 150 0.490 0.593 10.600 0.056 300 0.530 0.642 11.670 0.055 800 0.530 0.641 11.630 0.055 LTs 0.001 0.476 0.237 0.497 6.085 0.082 150 0.268 0.562 14.770 0.038 300 0.294 0.617 17.062 0.036 800 0.290 0.609 17.170 0.036 综上所述,在低速冲击情况下随着压缩速度的增加,相对破碎势增大,这意味着在动态压缩过程中发生了更多的颗粒破碎。但是,考虑到能量耗散也随着压缩速度的增加而增加。因此,研究了单位能量耗散所产生的相对破坏势。结果表明,LTs和LGs均有明显的下降趋势。结果表明,随着压缩速度的增加,单位能量输入下的颗粒破碎率降低。相反,为了产生相同的破坏程度,在动态条件下需要更高的能量输入。简而言之,与静态压缩相比,动态颗粒破碎时的能量利用率有所降低,因此,当单位压缩程度或破碎产生时,能量利用率的降低导致颗粒系统的能量输入增加,从而导致动态压缩下的比能量吸收增强。显然,能量吸收越高,动力响应上的应力水平越高。得出颗粒破碎特性的速度依赖性和能量利用率是产生名义应变率效应的根本原因。
4. 结 论
通过对不同粒径的漂珠堆积体在准静态和动态下的力学性能、粒径分布和Hardin破碎势进行分析,可以得到如下结论:
(1)不同粒径的漂珠在不同冲击速度下表现出明显的应变率效应,动态加载下颗粒材料的强度明显提升,LTs和LGs在0.001~800 s−1范围内材料的强度分别提高294%和320%;
(2)在相同加载速度下,粒径较小的颗粒强度较大粒径的强度和吸能效果都有明显提升,提升量分别为35%~40%和35%~48%,在限制应变为50%情况下大粒径的破碎率较小颗粒的破碎率高;
(3)随着应变率的提高,颗粒的破碎率和吸能效果都会明显提高,可以得出能量吸收是决定破碎比的主要影响因素,颗粒破碎特性的速率依赖性导致材料宏观应变率效应;
(4)通过Hardin破碎势对颗粒材料在准静态和动态压缩后粒径分布,可以发现相对破碎势与能量吸收有直接关系,动态冲击下颗粒破碎能量利用率降低,从而导致了堆积体宏观能量耗散和应力水平高于准静态压缩。
-
表 1 不同爆炸工况下冲击波壁压峰值实验值及理论值对比
Table 1. Comparison of theoretical and experimental results of wall pressure in test cases
No. me/kg R/m (√me/R)/(kg0.5⋅m−1) mp/(kg·m-2) α /(°) pwt/MPa pwe/MPa Cp 1 0.005 0.10 0.707 7.8 3 191.680 124.960 0.652 2 0.005 0.15 0.471 7.8 5 121.220 80.010 0.660 3 0.050 1.00 0.224 7.8 6 33.820 22.560 0.667 4 6.000 19.55 0.125 93.6 6 7.420 5.580 0.752 5 6.000 10.21 0.240 93.6 20 15.520 12.240 0.789 6 6.000 15.23 0.100 23.4 67 0.154 0.068 0.442 7 6.000 15.23 0.100 23.4 4 0.170 0.127 0.747 8 1.000 12.59 0.079 23.4 83 5.910 0.990 0.168 9 1.000 12.99 0.077 23.4 64 5.770 2.250 0.390 10 1.000 8.83 0.113 23.4 69 8.920 4.890 0.548 11 1.000 5.33 0.188 23.4 65 15.780 9.050 0.574 12 1.000 4.69 0.213 23.4 27 18.240 14.140 0.775 13 6.000 16.56 0.148 23.4 77 8.600 2.240 0.260 14 6.000 16.50 0.148 23.4 69 8.640 3.600 0.417 15 6.000 15.07 0.163 23.4 50 9.580 4.200 0.438 注:工况1~3、5~8、13~15为函数f(α)拟合样本点工况;工况4、9~12为考查点工况。 表 2 不同爆炸工况下冲击波衰减常数实验值及理论值对比
Table 2. Comparison of theoretical and experimental results of θ in test cases
No. me/kg R/m (√me/R)/(kg0.5⋅m−1) mp/(kg·m-2) α /(°) θpt/μs θpe/μs ε/% 1 0.005 0.10 0.707 7.8 3 3.5 3.2 8.6 2 0.005 0.15 0.471 7.8 5 3.6 3.6 0.0 3 0.050 1.00 0.224 7.8 6 4.5 4.5 0.0 4 6.000 10.21 0.240 93.6 20 13.6 13.3 2.2 5 2.000 15.21 0.093 23.4 81 2.5 2.3 8.0 -
[1] Jiang J, Olson M D. Rigid-plastic analysis of underwater blast loaded stiffened plates[J]. International Journal of Mechanic and Science, 1995, 37(8): 843-859. doi: 10.1016/0020-7403(94)00100-X [2] Makinen K. Cavitation models for structures excited by a plane shock wave[J]. Journal of Fluids and Structures, 1998, 12: 85-101. doi: 10.1006/jfls.1997.0120 [3] Santiago L D. Fluid-interaction and cavitation effects on a surface ship model due to an underwater explosion[R]. California: Naval Postgraduate School, 1997. [4] 金辉, 周学滨, 周华, 等.水下爆炸中自由场压力和船体壁压的测量与分析[J].海军工程大学学报, 2009, 21(5): 82-87.Jin Hui, Zhou Xue-bin, Zhou Hua, et al. Measurement and analysis of free-field pressure and ship hull pressure under water explosion[J]. Journal of Naval University of Engineering, 2009, 21(5): 82-87. [5] 程素秋, 樊宝顺, 薛飞, 等.水下非接触爆炸作用下舱段模型的动态响应[J].爆炸与冲击, 2008, 28(4): 360-366. doi: 10.3321/j.issn:1001-1455.2008.04.013Cheng Su-qiu, Fan Bao-shun, Xue Fei, et al. Dynamical response measurement of a cabin model subjected to noncontact underwater explosion[J]. Explosion and Shock Waves, 2008, 28(4): 360-366. doi: 10.3321/j.issn:1001-1455.2008.04.013 [6] 樊宝顺, 程素秋, 韩峰.舰船舱段模型在水下爆炸作用下的壁压分析[J].中国舰船研究, 2009, 4(5): 20-22.Fan Bao-shun, Cheng Su-qiu, Han Feng. Wall pressure analysis of cabin model subjected to underwater explosion[J]. Chinese Journal of Ship Research, 2009, 4(5): 20-22. [7] 李海涛, 朱锡, 张振华.水下爆炸球面冲击波作用下船体梁的刚塑性动响应特性[J].工程力学, 2010, 27(10): 202-207.Li Hai-tao, Zhu Xi, Zhang Zhen-hua. Dynamic rigid-plasitc response of ship-like beam subjected to underwater spherical shockwaves[J]. Engineering Mechanics, 2010, 27(10): 202-207. [8] 顾文彬, 苏青笠, 刘建青, 等.水下平面爆炸冲击波作用下空化区域形成及其特征[J].爆破, 2004, 21(4): 8-11.Gu Wen-bin, Su Qing-li, Liu Jian-qiang, et al. Formation process of cavitation region affected by underwater plane blast shock wave and its characteristics[J]. Blasting, 2004, 21(4): 8-11. -