A double-striker electromagnetic driving SHPB system for soft materials
-
摘要: 软材料的动态力学性能研究一直备受关注,目前分离式Hopkinson压杆(split Hopkinson pressure bar, SHPB)技术是其最重要的测试手段,然而在测试超软材料时实验装置设计方面仍存在许多有待改进之处。本文中研制了一套双子弹电磁驱动SHPB系统,使用聚碳酸酯作为杆件材料以克服软材料试件带来的诸多困难,引入了双子弹设计方案解决了电磁驱动方式难以应用于非铁磁材料的问题,并有效保证了子弹速度的准确控制。使用双子弹电磁驱动SHPB系统和传统金属SHPB装置同时对硅胶材料的动态力学性能进行了测试,实验结果的吻合性验证了本套系统的可靠性。应用双子弹电磁驱动SHPB系统开展了聚乙烯醇(polyvinyl alcohols, PVA)水凝胶这种超软材料在高应变率下的实验,成功表征出其动态力学性能。
-
关键词:
- 分离式Hopkinson压杆(SHPB) /
- 电磁驱动 /
- 双子弹 /
- 软材料
Abstract: The dynamic mechanical properties of soft materials have attracted great attention, and the separated Hopkinson pressure bar (SHPB) technique is the most important test method at present. However, the SHPB technique still needs some improvements when it is applied to extra-soft materials. So, a double-striker electromagnetic driving SHPB system was developed. In the developed system, the bars were made of polycarbonate, and the double-striker scheme was designed to precisely control the velocities of the strikers. Both the aluminum SHPB facility and the double-striker electromagnetic driving SHPB system were employed to carry out dynamic experiments on the silicone, and the agreement between the test results by these two systems indicated the reliability of the innovative system. And dynamic mechanical properties of PVA hydrogels were successfully tested by the double-striker electromagnetic driving SHPB system. -
可燃气体着火过程的影响因素、燃烧特性、传播特性等一直是安全领域的研究热点,已取得了许多成果。其中对单质可燃气体的研究成果较多,对混合可燃气体的研究较少,如油气。油气是一种主要由轻质烷烃组成的混合物。混合物的研究较困难,如对汽油燃烧特性的研究,往往采用其替代物进行研究[1-3]。Du等[4]已通过实验研究了油气着火的影响因素,欧益宏等[5]、杜扬等[6]、Ou等[7]已通过实验研究了受限空间爆炸的传播特性。但受实验条件限制,某些研究还有待深入,如温度、压力、火焰速度、湍流变化速度、组分变化速度等参数的演变特征。数值模拟是研究可燃气体着火、爆炸过程的一种有效方法,如Bi等[8]、Wang等[9]、Skjold等[10]研究了具有某种特殊几何结构的封闭管道内的气体爆炸过程,Sarli等[11]研究了有障碍物直管道内的气体爆炸过程。通过数值模拟能获得采用实验方法不能测试和观察到的一些信息,而且还具有安全、经济的特点。但是数值模拟的关键是建立恰当的模型、反应机理和初始条件。吴松林等[12]对油气的化学反应机理进行了简化,为油气着火、爆炸的数值模拟奠定了基础。本文中,基于油气简化机理[12]和实验结果[13],将油气化学动力学、辐射模型和传质相互耦合,通过建立统一控制方程组模型,选择适当的数值解法和初始条件,对持续热壁条件下油气热着火的发生进行数值模拟,以期获得的相关数值结果可作为实验结果的有益补充和启发进一步探讨。
1. 统一控制方程组模型和数值解法
基于实验背景,对受限空间中油气热着火的发生作如下假设:(1)混合气体处于绝热环境中,受限空间与外界的热量传递可忽略;(2)混合气体处于可压缩、低马赫数的二维流动状态;(3)黏性耗散、压缩功的影响可忽略。
油气着火过程是流动过程和组分化学反应过程的耦合。在受限空间持续热壁加热情况下,还要同时考虑热壁与气体之间的对流换热、气体间能量的辐射。因此,对油气热着火需要同时建立化学反应模型、对流换热模型、辐射模型以及综合上述模型的流场控制方程组模型。对化学反应采用吴松林等[12]提出的47个组分、100个基元反应的油气简化机理;对辐射模型采用P-1模型;对加热热壁与可燃气体的对流采用牛顿冷却公式;其他的控制方程有连续方程、动量方程、能量方程、组分方程、k方程、ε方程。化学反应和流场方程通过组分和能量关系直接耦合,辐射模型、对流换热方程与能量方程直接耦合。
统一模型中基本方程组可表示为瞬态项、对流项、扩散项和源项组成通用输运方程:
∂∂t(ρφ)+div(ρuφ)=div(Γ⋅∂φ∂x)+S 式中:φ代表温度、质量分数、单位体积的能量和单位体积的动量等控制变量,Γ为对应的交换系数,S为对应的源项,u为气体的速度。
由于这类方程相互耦合,有强烈的非线性,特别是源项, 因此,采用迭代法进行数值求解,即将方程表示为节点的迭代关系,把连续的基本方程离散化,获得迭代方程。在计算方法上采用有限体积法,在体积域采用混合格式,时间域上采用向后差分格式对控制方程进行离散。控制方程的分离式解法如图 1所示。
2. 初始条件及结果验证
数值模拟区域采用和实验[13]一致的1 700 mm×400 mm绝热受限空间,如图 2所示。模拟区域共分1 400个单元、2 910个内部面、1 511个点。初始压强为一个标准大气压,环境温度为300 K,热壁温度为2 000 K。初始组分如表 1所示。以下通过油气热着火发生过程的流场特征和某些点参数的变化曲线来分析着火特征,选取的7个点位置见图 2。图 3中给出了通过数值模拟获得的热壁表面3个点的温度变化曲线。当热爆炸发生时,热壁表面的临界温度为880 K,这与实验结果877 K[13]基本一致。
表 1 初始组分Table 1. Initial components初始组分 质量分数 摩尔分数 摩尔浓度/(mol·cm-3) CH4 0.01 0.018 2 7.45×10-7 C2H6 0.01 0.009 7 3.98×10-7 C3H8 0.01 0.006 6 2.71×10-7 C4H10 0.01 0.005 0 2.05×10-7 n-C7H16 0.01 0.002 9 1.19×10-7 i-C8H18 0.01 0.002 6 1.04×10-7 O2 0.21 0.192 1 7.85×10-6 N2 0.73 0.763 0 3.12×10-5 3. 数值结果及分析
3.1 温度
3.1.1 温度场变化过程
选取有代表性的6个时间点来刻画温度场和参数的变化过程,如表 2所示,表中t为时间,n为计算步数,Tmax为最高温度,pmax为最高压力,vx, max为最高横向速度,vy, max为最高纵向速度, I为最大扰动强度,其温度场模拟结果如图 4所示。
表 2 热着火过程的6个时间点的基本参数Table 2. Basic parameters of the thermal ignition process at 6 time points时间点 t/ms n Tmax/K pmax/Pa vx, max/(m·s-1) vy, max/(m·s-1) I/% No.1 100.09 10 000 452.12 10 484.3 4.93 4.76 280.551 No.2 413.06 63 000 581.78 92 218.8 66.89 52.59 2 435.931 No.3 413.53 64 800 844.91 127 912.0 109.51 55.40 2 812.672 No.4 413.53 80 000 872.86 144 851.9 100.20 42.66 2 266.638 No.5 413.53 100 000 914.99 179 834.6 99.01 40.38 2 195.186 No.6 414.06 587 200 1271.99 381 857.6 306.55 142.44 14 418.030 图 4(a)显示热对流换热和热辐射后,加热壁上空温度场呈笔尖状,气体流动和升温同时发生;图 4(b)温度场显然受到了流动的影响,其形状有所变化,向一侧倾斜,这同实验的可视化结果[13]完全一致;图 4(c)~(e)显示的现象几乎是瞬间完成的,热着火已经发生,形成燃烧核,同时显示了温度场变化过程也是化学反应和流动的耦合过程,使得燃烧向上和向两侧扩张;图 4(f)显示已经形成了次级燃烧火焰,向两侧扩张,弱冲击波也已经形成。
3.1.2 不同位置的温度突变
(1) 温度存在突变现象,整个热着火发生过程可分成两段,分别为温度受加热壁和缓慢氧化放热影响的温度缓慢上升过程,以及热着火发生后温度跃升过程;
(2) 在温度缓慢上升过程阶段,7个空间点的温度基本呈现与空间相关上升,空间对称两侧的温度并不一定对称;
(3) 在温度缓慢上升过程阶段,受流场的影响,7个点的温度稳定上升;
(4) 从图 5(b)可以发现,尽管起燃时间很短,但还存在2个小的温度突变阶段,说明起燃中存在非线性特性,其原因需要深入研究;
(5) 受流场影响,燃烧过程呈现不对称性。
3.2 压力
3.2.1 压力变化过程
选取与温度场同样的6个时间点,作受限空间的压力等值线, 如图 6所示。图 6(a)中,随着气体受热壁影响,压力开始在热壁上方上升,压缩受限空间的气体,使得受限空间的气体开始向四周流动。从图 6(b)、(c)开始,热着火发生了,压力等值线开始变得稠密。模拟数据结果显示最大压力有1个多大气压,这与实验结果基本一致[13]。图 6(d)中,压力等值线较弯曲,压缩两侧的气体,形状不规则。图 6(e)和(f)的等值线在两侧非常稠密,说明随着气体燃烧强度的增大,两侧的弱冲击波已经形成。但是,燃烧火焰还主要在热壁的中上部,最高压力还处于中间,燃烧化学反应和气体流动的强烈耦合还需要更大的空间,所以总的压力并不高,模拟数据显示最高压力有3个多大气压。
3.2.2 不同位置压力突变
图 7是7个点的压力变化曲线,从图 7可知:(1)热着火发生前,7个点的压力变化不大;(2)热着火发生后,压力也存在阶段性突变特征。在同一时刻,7个点的压力不同,存在很大的差别;(3)从模拟的最后时刻的压力来看,第3~6点的压力很接近,而第7点压力很低,这也说明弱冲击波正在形成。
3.3 主要组分质量变化
3.3.1 主要反应物变化
图 8是1号位置的主要反应物质量变化曲线。从图 8可知,t∈[0, 0.36] s时,主要反应物相差不大,它们的质量分数从0.01还降到了0.009 6,这主要是受热壁加热影响,密度下降的原因。t∈[0.36,0.415] s时,主要反应物呈现了明显的阶段性特征。第1阶段,从0.36 s到0.39 s,主要反应物发生了缓慢氧化反应,浓度呈现下降趋势,特别正庚烷发生了热裂解,质量分数几乎变为了零。第2阶段,从0.39 s到0.40 s,主要反应物质量分数有所增大。主要原因是密度回升,流场影响相对减弱,大分子热裂解成低碳烷烃的因素增加,造成了C1~C4的质量分数增大。第3阶段,从0.40 s到0.41 s,主要反应物发生反应被消耗,浓度降低。0.41 s后,热着火发生,反应速度增大,流场变化加剧,主要反应物质量分数呈现波动。
3.3.2 不同位置的组分质量分数
图 9~13为不同位置主要组分的质量分数变化曲线。从以上反应物在不同阶段的质量分数变化曲线可知:(1)t∈[0, 0.15] s阶段,主要生成物CO、H2O基本没有变化,但随后到0.39 s,这些生成物有较小量级增加,说明随着温度的增加,主要反应物发生了非常缓慢的氧化反应,特别是1号位置表现明显。(2)t∈[0.39, 0.41] s阶段,发生缓慢氧化阶段,生成物不断增加,组分的质量分数与反应强度,流场变化特征、空间位置关系紧密,CO、H2O等组分的质量分数总体增加,O2、轻质烷烃等组分质量分数总体减小,但波动较明显。(3)热着火发生后,主要受燃烧反应强度的影响,反应物减少,生成物增加。从OH基团的变化可以发现不同位置反应强度不同。
综上,受限空间油气热着火发生过程中,无论是缓慢氧化阶段,还是热着火发生过程中,整个受限空间内气体都受到了化学反应和流动的耦合作用,气体间化学反应不断加剧,流动不断增强。但从所起的主导作用来看,加热初期,辐射换热为主导作用;在加热中期,辐射换热和缓慢氧化反应为主导作用;在热着火发生阶段,化学反应和流动强烈耦合,同时占主导作用。
3.4 火焰速度
3.4.1 层流火焰速度
图 14~17为不同位置流场速度变化曲线。由横轴方向层流速度来看:热着火发生前,速度很小,仅有每秒几个厘米,左右流动也基本对称;受受限空间影响,气体流动发生了回流,两侧气体呈现向中间流动的现象。在0.405 s后,热着火发生,形成弱冲击波,气体向两侧流动。速度超过了100 cm/s。从第1和第2点的流动方向来看,热着火发生前,火焰基本呈振动状态,方向不断摆动,热着火发生后火焰偏向右侧,说明有湍流影响流动方向。由纵轴方向的层流速度来看:与横向速度比较,开始加热阶段,最大速度略大,而后受上壁面影响,速度相对减小。起燃阶段,速度方向总体向上,呈现震荡特征。
3.4.2 湍流火焰速度
图 18为不同位置湍流变化曲线。湍流也出现了2个阶段的特征:在加热开始阶段,湍流突然增大,而后略有降低;当热着火发生时,湍流急剧增大。
4. 结论
油气在热壁下的热着火发生过程主要受辐射放热、化学反应和流动的耦合作用,呈现出阶段性。在高温热壁加热初始阶段,加热使得气体膨胀,同时气体温度升高。流动改变了流场流动特征,造成受限空间中组分浓度、流速等参数的变化,辐射放热和流动成为主导作用;在热着火发生前,随着辐射放热进行,热壁上方温度升高,油气缓慢氧化反应放热逐渐成为影响流场特征的主导作用。在热着火发生后,剧烈化学反应和高强度湍流相互耦合,弱冲击波瞬间形成,化学反应和流动同时占主导作用,气体的温度、压力、质量分数、湍流强度等参数均呈现出一定的震荡和非线性特性。
-
表 1 聚碳酸酯杆件的长度
Table 1. Lengths of polycarbonate bars
组号 子弹长度/mm 入射杆长度/mm 透射杆长度/mm 1 200 600 300 2 400 1 000 600 表 2 子弹速度与充电电压的关系
Table 2. Relation between striker velocity and charging voltage
充电电压/V 长200 mm子弹的
速度/(m·s−1)长400 mm子弹的
速度/(m·s−1)充电电压/V 长200 mm子弹的
速度/(m·s−1)长400 mm子弹的
速度/(m·s−1)100 1.52±0.10 1.61±0.08 260 2.86±0.09 4.37±0.08 120 1.82±0.09 1.91±0.07 280 2.90±0.10 4.44±0.09 140 1.98±0.09 2.26±0.07 300 3.13±0.11 4.64±0.08 160 2.14±0.09 2.60±0.08 320 3.19±0.11 4.85±0.10 180 2.30±0.10 2.97±0.08 340 3.28±0.10 5.07±0.09 200 2.46±0.10 3.35±0.07 360 3.32±0.11 5.28±0.09 220 2.55±0.09 3.98±0.07 380 3.52±0.12 5.57±0.10 240 2.67±0.09 4.17±0.08 -
[1] TAYLOR G I. The testing of materials at high rates of loading [J]. Journal of the ICE, 1946, 26(8): 486–519. [2] 胡时胜, 王礼立, 宋力, 等. Hopkinson压杆技术在中国的发展回顾 [J]. 爆炸与冲击, 2014, 34(6): 641–657. DOI: 10.11883/1001-1455(2014)06-0641-17.HU Shisheng, WANG Lili, SONG Li, et al. Review of the development of Hopkinson pressure bar technique in China [J]. Explosion and Shock Waves, 2014, 34(6): 641–657. DOI: 10.11883/1001-1455(2014)06-0641-17. [3] CHEN W, LU F, FREW D J, et al. Dynamic compression testing of soft materials [J]. Journal of Applied Mechanics, 2002, 69(3): 214–223. [4] 卢芳云, CHEN W, FREW D J. 软材料的SHPB实验设计 [J]. 爆炸与冲击, 2002, 22(1): 15–19.LU Fangyun, CHEN W, FREW D J. A design of SHPB experiments for soft materials [J]. Explosion and Shock Waves, 2002, 22(1): 15–19. [5] SONG B, GE Y, CHEN W W, et al. Radial inertia effects in Kolsky bar testing of extra-soft specimens [J]. Experimental Mechanics, 2007, 47(5): 659–670. DOI: 10.1007/s11340-006-9017-5. [6] ZHAO H, GARY G, KLEPACZKO J R. On the use of a viscoelastic split Hopkinson pressure bar [J]. International Journal of Impact Engineering, 1997, 19(4): 319–330. DOI: 10.1016/s0734-743x(96)00038-3. [7] CHEN W, ZHANG B, FORRESTAL M J. A split Hopkinson bar technique for low-impedance materials [J]. Experimental Mechanics, 1999, 39(2): 81–85. DOI: 10.1007/BF02331109. [8] 林玉亮, 卢芳云, 卢力. 石英压电晶体在霍普金森压杆实验中的应用 [J]. 高压物理学报, 2006, 19(4): 299–304.LIN Yuliang, LU Fangyun, LU Li. The application of quartz transducer technique in SHPB [J]. Chinese Journal of High Pressure Physics, 2006, 19(4): 299–304. [9] 宋博, 姜锡权, 陈为农. 霍普金森压杆实验中的脉冲整形技术 [C] // 第三届全国爆炸力学实验技术交流会. 合肥: 2004. [10] SONG B, CHEN W. One-dimensional dynamic compressive behavior of EPDM rubber [J]. Journal of Engineering Materials and Technology, 2003, 125(3): 294–301. [11] 王宝珍, 胡时胜. 软质聚氨酯泡沫的冲击力学性能 [J]. 高分子材料科学与工程, 2009, 25(12): 48–51.WANG Baozhen, HU Shisheng. Mechanical properties of flexible polyurethane foams under impact loading [J]. Polymer Materials Science and Engineering, 2009, 25(12): 48–51. [12] SHARMA A, SHUKLA A, PROSSER R A. Mechanical characterization of soft materials using high speed photography and split Hopkinson pressure bar technique [J]. Journal of Materials Science, 2002, 37(5): 1005–1017. [13] 郭伟国, 赵融, 魏腾飞, 等. 用于Hopkinson压杆装置的电磁驱动技术 [J]. 实验力学, 2010(6): 682–689.GUO Weiguo, ZHAO Rong, WEI Tengfei, et al. Electromagnetic driving technique applied to split-Hopkinson pressure bar device [J]. Journal of Experimental Mechanics, 2010(6): 682–689. [14] LIU Zhanwei, CHEN Ximin, LV Xintao, et al. A mini desktop impact test system using multistage electromagnetic launch [J]. Measurement, 2014, 49: 68–76. [15] 刘战伟, 吕新涛, 陈喜民, 等. 基于多级电磁发射的mini-SHPB装置 [J]. 实验力学, 2013, 28(5): 557–562.LIU Zhanwei, LV Xintao, CHEN Ximin, et al. On the mini-SHPB device based on multi-level electromagnetic emissions [J]. Journal of Experimental Mechanics, 2013, 28(5): 557–562. [16] BACON C. An experimental method for considering dispersion and attenuation in a viscoelastic Hopkinson bar [J]. Experimental Mechanics, 1998, 38(4): 242–249. DOI: 10.1007/bf02410385. [17] 刘孝敏, 胡时胜, 陈智. 粘弹性Hopkinson压杆中波的衰减和弥散 [J]. 固体力学学报, 2002, 23(1): 81–86.LIU Xiaomin, HU Shisheng, CHEN Zhi. The wave propagation attenuation and dispersion in a viscoelastic Hopkinson pressure bar [J]. Acta Mechanica Solida Sinica, 2002, 23(1): 81–86. [18] YANG L M, SHIM V P W. An analysis of stress uniformity in split Hopkinson bar test specimens [J]. International Journal of Impact Engineering, 2005, 31(2): 129–150. DOI: 10.1016/j.ijimpeng.2003.09.002. [19] 王宝珍, 胡时胜. 猪后腿肌肉的冲击压缩特性实验 [J]. 爆炸与冲击, 2010, 30(1): 33–38. DOI: 10.11883/1001-1455(2010)01-0033-06.WANG Baozhen, HU Shisheng. Dynamic compression experiments of porcine ham muscle [J]. Explosion and Shock Waves, 2010, 30(1): 33–38. DOI: 10.11883/1001-1455(2010)01-0033-06. [20] GORHAM D A. Specimen inertia in high strain-rate compression [J]. Journal of Physics D: Applied Physics, 1989, 22(12): 1888–1893. DOI: 10.1088/0022-3727/22/12/014. [21] GORHAM D A. The effect of specimen dimensions on high strain rate compression measurements of copper [J]. Journal of Physics D: Applied Physics, 1991, 24(8): 1489–1492. DOI: 10.1088/0022-3727/24/8/041. -