Lagrangian forward analysis in data processing of ramp wave compression experiments
-
摘要: 提出了一种联合使用Lagrange方法和转换函数方法来处理斜波压缩实验数据的新途径,和传统的数据处理方法相比,该方法更适合处理材料的复杂力学响应并具有较好的精度,同时该方法使用过程中对样品材料参数的初猜值精度要求更低。分析了转换函数方法的可靠性和健壮性,讨论了转换函数和Lagrange方法联合在斜波压缩强度实验数据处理中的应用,获得了可靠的结果。
-
关键词:
- Lagrange方法 /
- 转换函数 /
- 数据处理 /
- 斜波压缩
Abstract: In the present work, we developed a novel method combining the forward Lagrange method and the transfer function method to process the data of ramp wave compression experiments. Compared with traditional methods, this method is more suitable to process the data containing materials' complex behaviors and produces more accurate results. Meanwhile, this method has lower precision requirement regarding the initial gauss of testing material parameters. The feasibility and robustness of the transfer function method was analyzed, and the application of the forward Lagrange method and the transfer function method were examined in forward data processing of strength measurement experiments under ramp wave compression.-
Key words:
- Lagrangian method /
- transfer function method /
- data processing /
- ramp wave compression
-
空腔爆炸中,爆后空腔内气体的温度及压力变化历程对于分析评估爆炸实施效果及气体封闭和扩散规律具有重要意义。准确测量空腔内温度、压力变化历程存在以下3个难点:(1)爆炸瞬时产生的高压冲击波及高温灼烧,对传感器生存和正常工作构成了极大威胁,传感器的长期生存率较低;(2)爆炸产生的强电磁脉冲,会对采集系统造成严重损伤,测量系统的抗电磁防护要求较高;(3)要完整取得爆后空腔内气体的温度、压力变化历程,传感器的量程选择存在一定难度,既要保证爆炸瞬时能够生存,又要求量程不能过高,以免温度压力变化历程测量后期的数据不确定度过大。在相似领域内,李媛媛等[1]通过压力传感器和热电偶进行了半密闭条件下爆炸场的温度和压力测量,获得了不同爆炸条件下爆炸场温度和压力的响应特征,该系统适用于半密闭条件下爆炸温度、压力测量。王长利等[2]等使用锥体的形式安装侵蚀热电偶,获得小比距爆炸瞬时温度。李芝绒等[3-4]、马红等[5]、王代华等[6]、王等旺等[7]、黄亚峰等[8]、王铁良等[9]分别对密闭内爆炸条件下温度及准静态压力的特征进行了深入研究,对空腔爆炸温度、压力变化历程测量具有较大的指导意义。本文中在参考已有研究成果的基础上,基于铠装K型热电偶和压阻式绝压传感器,设计密封隔热防护装置,建立空腔爆炸静态温度、压力测量系统,以期实现小比距离密闭空腔爆炸爆后气体温度、压力变化历程的完整测量。
1. 测量系统设计
1.1 测量系统基本组成
如图 1所示,系统前端为温度、压力传感器,通过安全隔离后,信号被数据采集单元采集并进行转换;中央控制单元接收转换后的数字信号,同时根据监测实施程序,远程控制数据采集单元和前端信号的通断,以保护后端采集系统免受爆炸瞬时电磁脉冲冲击。同时中央控制单元还完成远程数据传输显示功能,并接收远程终端的指令控制。
1.2 安全隔离单元设计
在传感器单元与系统采集单元之间,应采取隔离防护措施[10],当前端在测量过程中出现瞬间的电流或电压冲击时,能够阻断其传导至后端采集和控制系统,保证系统核心单元的生存安全。该防护采用的技术方法是:安全隔离单元能够接收来自前端传感器的正常测量信号,同时在经过信号隔离转换后,能够无失真地保持原输入信号的大小,并传输至信号采集系统。同时,当信号线上出现冲击电流或电压时,安全隔离单元能够有效地阻止该冲击,从而保证后端采集和控制系统的正常运行。在本系统中,选用上海晨竹生产的FBE054-AAAA型信号隔离安全栅对系统采集仪器进行隔离保护,每块安全栅共有2个通道,可以对2个传感器输出的标准电流信号进行无失真隔离传输,从而有效地保护了后端的采集控制仪器,如图 2所示。
2. 密封隔热防护装置设计与分析
2.1 安装防护结构设计
爆炸瞬时会产生高温高压冲击波,尤其在小比距、大当量爆炸情况下,爆炸瞬时产生的高温高压冲击波可能会直接摧毁传感器,为确保测量取得完整的气体温度、压力信号曲线,需对传感器进行抗爆防护,设计如图 3所示传感器安装防护结构。该结构安装在爆炸迎爆面,主要由锥形引气装置、保护挡板、引气管、安装隔板组成,整个防护结构浇注在混凝土中。爆炸产生的高温高压气体通过锥形引气装置和引气管进入密闭气室,密闭气室和传感器安装室之间有一道密封隔板,将两个气室隔开。通过密封螺母将温度、压力传感器固定在安装隔板上,安装完成后,传感器的敏感面与爆室连通,信号调理模块与爆室隔开,可以有效避免爆炸产生的高温高压气体损坏传感器的信号调理模块。
2.2 锥形导气装置抗冲击性能分析
图 4(a)所示的锥形导气装置为防护结构最关键部件,该装置的抗冲击性能直接关系传感器的生存与否。其设计参数如下:最大外径10 cm,内孔直径2 cm,总高15 cm,侧壁通孔直径1 cm,锥角120°。装置采用Q345钢整体加工,通过焊接的方式固定在厚度不低于30 cm的保护挡板板上。材料参数如表 1所示,其中:E为杨氏模量;EP为硬化模量;ν为泊松比;σy为屈服强度;ρ为密度。
表 1 材料参数Table 1. Mechanical parameters of relevant materials材料 E/GPa EP/GPa ν σy/MPa ρ/(g·cm-3) Q345钢 210 10 0.28 345 7.85 为分析该装置在小比距、大当量爆炸条件下的抗冲击性能,按照其实际尺寸建立三维计算模型,网格平均尺寸为2.5 mm,网格划分效果如图 4(b)所示。由于装置所处环境比较复杂,分析时没有考虑爆炸冲击波在爆室壁面的反射会聚作用,仅考虑小比距下的正反射冲击波。根据Kingery-Bulmash公式,假设装置的比距离为0.86 m/kg1/3。引爆后,冲击波约到达锥顶的反射压力峰值为6 MPa,正相持续时间约48 ms。计算时,在支座的锥形顶部施加爆炸载荷,由Load-Blast关键字施加,底面施加位移约束,计算时长8 ms, 爆炸冲击作用下不同时刻装置的von Mises应力分布如图 5所示。
由图 5可以看出,爆炸波到达时刻,von Mises应力最大值出现在锥顶,峰值约8 MPa,0.17 ms后爆炸波传至支座底部,在圆孔周围von Mises应力最大,约为21 MPa。随着时间的增加,应力值始终未超过21 MPa,说明装置在爆炸波作用下没有出现塑性变形,结构强度满足要求。
考虑到载荷作用时间较长,类似准静态压力的作用,根据计算,装置的等效静载系数约为1.1,按安全系数为1.2考虑,取静态压力8 MPa作用于支座顶部,底部施加位移约束,分析防护装置在静力作用下的变形。
由图 6可知,模型最大位移位于支座锥顶靠近开孔位置,以及圆孔所对的圆柱侧壁,位移峰值约0.05 mm;模型最大应变位于支座开孔位置,应变峰值约为1.88×10-4,其他部位应变小于1.00×10-4,没有出现塑性变形。由计算结果可知,结构刚度也能够满足冲击防护要求。
3. 实验测试及结果分析
3.1 传感器选型及量程选择
某密闭空腔爆炸实验,比距离为0.86 m/kg1/3,假设爆炸后气体初态温度为500 ℃,爆后爆室内的静态气压峰值计算值约为1.5 MPa。综合考虑性能特征及安装便捷性,选用沈阳东大传感技术有限公司生产的铠装K型热电偶作为测温器件,其基本性能参数为:温度测量范围0~1 250 ℃,响应时间0.73 s,精确度±2.2 ℃;选用宝鸡麦克传感技术有限公司生产压阻式绝压传感器作为测压器件,其基本性能参数为:压力测量范围0~2 MPa,精度为全量程±0.25%,过载压力为全量程的1.5倍。
3.2 测量结果及分析
图 7是爆后空腔内气体温度变化历程监测结果,图 8是气体压力变化历程监测结果。
经过数据处理,爆后17 s, 热电偶获得爆后气体温度测量峰值708.6 ℃,爆后25 s, 压力传感器测量获得爆室内气体压力测量峰值1518.5 kPa。随后,爆室内气体温度、压力呈快速下降趋势,最终传感器输出可恢复至初始值,说明在测量过程中,传感器并未发生损坏,测量数据能够真实反应爆室内气体温度、压力变化趋势。
4. 结论
(1) 设计的温度、压力测量系统能够完整的测量取得小比距、密闭条件下爆后气体温度、压力变化历程。通过锥形导气结构,有效的减小了爆炸瞬时高温高压对传感器的损伤,通过密封隔热设计,可确保传感器在复杂环境中生存,该方法可以应用到其他类似实验场合中。
(2) 密闭空腔爆后压力下降是气体温度下降和气体总量减小2个因素共同作用的结果, 在爆后空腔内气体总量分析时,需引入温度测量结果进行修正。
(3) 通过合理有效的防护设计,可以通过降低压力传感器的量程,提高密闭空腔爆后气体压力测量精度。
-
-
[1] 孙承纬, 赵剑衡, 王桂吉, 等.磁驱动准等熵平面压缩和超高速飞片发射实验技术原理、装置及应用[J].力学进展, 2012, 42(3):206-219. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201200497977Sun Chengwei, Zhao Jianheng, Wang Guiji, et al. Progress in magnetic loading techniques for isentropic compression experiments and ultra-high velocity flyer launching[J]. Advances in Mechanics, 2012, 42(3):206-219. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201200497977 [2] 李牧, 孙承纬, 赵剑衡.固体材料高功率激光斜波压缩研究进展[J].爆炸与冲击, 2015, 35(2):145-156. http://www.bzycj.cn/CN/abstract/abstract9440.shtmlLi Mu, Sun Chengwei, Zhao Jianheng. Progress in high-power laser ramp compression of solid[J]. Explosion and Shock Waves, 2015, 35(2):145-156. http://www.bzycj.cn/CN/abstract/abstract9440.shtml [3] Wang G, Zhao J, Zhang H, et al. Advances in quasi-isentropic compression experiments at institute of fluid physics of CAEP[J]. The European Physical Journal Special Topics, 2012, 206(1):163-172. doi: 10.1140/epjst/e2012-01597-y [4] 王桂吉, 赵剑衡, 孙承纬, 等.磁驱动准等熵加载装置CQ-4的加载能力及主要应用[J].实验力学, 2015, 30(2):252-263. http://d.old.wanfangdata.com.cn/Periodical/sylx201502016Wang Guiji, Zhao Jianheng, Sun Chengwei, et al. On the loading capability and main application of magnetically driven quasi-isentropic compression device CQ-4[J]. Journal of Experimental Mechanics, 2015, 30(2):252-263. http://d.old.wanfangdata.com.cn/Periodical/sylx201502016 [5] Wang G, Luo B, Zhang X, et al. A 4 MA, 500 ns pulsed power generator CQ-4 for characterization of material behaviors under ramp wave loading[J]. Review of Scientific Instruments, 2013, 84(1):015117. doi: 10.1063/1.4788935 [6] Davis J P, Deeney C, Knudson M D, et al. Magnetically driven isentropic compression to multimegabar pressures using shaped current pulses on the Z accelerator[J]. Physics of Plasmas, 2005, 12(5):056310. doi: 10.1063/1.1871954 [7] Davis J P. Experimental measurement of the principal isentrope for aluminum 6061-T6 to 240 GPa[J]. Journal of Applied Physics, 2006, 99(10):103512-103512-6. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=77753de7c645e2e3af9b3c6e90d85bd8 [8] Davis J P, Brown J L, Knudson M D, et al. Analysis of shockless dynamic compression data on solids to multi-megabar pressures: Application to tantalum[J]. Journal of Applied Physics, 2014, 116(20):204903. doi: 10.1063/1.4902863 [9] Smith R F, Eggert J H, Jeanloz R, et al. Ramp compression of diamond to five terapascals[J]. Nature, 2014, 511(7509):330-333. doi: 10.1038/nature13526 [10] Luo Binqiang, Wang Guiji, Mo Jianjun, et al. Verification of conventional equations of state for tantalum under quasi-isentropic compression[J]. Jounal of Applied Physics, 2014, 116(19):193506. doi: 10.1063/1.4902064 [11] Smith R F, Eggert J H, Jankowski A, et al. Stiff Response of aluminum under ultrafast shockless compression to 110 GPa[J].Physical Review Letters, 2007, 98(6):065701. doi: 10.1103/PhysRevLett.98.065701 [12] Asay J R, Ao T, Davis J, et al. Effect of initial properties on the flow strength of aluminum during quasi-isentropic compression[J]. Journal of Applied Physics, 2008, 103(8):083514. doi: 10.1063/1.2902855 [13] Ao T, Knudson M D, Asay J R, et al. Strength of lithium fluoride under shockless compression to 114 GPa[J]. Journal of Applied Physics, 2009, 106(10):103507-103507-12. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=1e23ba0afada6fde67baacf6ad99e237 [14] Brown J L, Alexander C S, Asay J R, et al., Flow strength of tantalum under ramp compression to 250 GPa[J].Journal of Applied Physics, 2014.115(4):043530. doi: 10.1063/1.4863463 [15] 罗斌强, 王桂吉, 谭福利, 等.磁驱动准等熵压缩下LY12铝的强度测量[J].力学学报, 2014.46(2):241-247. http://www.cnki.com.cn/Article/CJFDTOTAL-LXXB201402009.htmLuo Binqiang, Wang Guiji, Tan Fuli, et al. Measurement of dynamic strength of LY12 aluminum under magnetically driven quasi-isentropic compression[J]. Chinese Journal of Theoretical and Applied Mechanics, 2014, 46(2):241-247. http://www.cnki.com.cn/Article/CJFDTOTAL-LXXB201402009.htm [16] Smith R F, Eggert J H, Swift D C, et al. Time-dependence of the alpha to epsilon phase transformation in iron[J]. Journal of Applied Physics, 2013, 114(22):223507. doi: 10.1063/1.4839655 [17] Rigg P A, Greeff C W, Knudson M D. Influence of impurities on the a to w phase transition in zirconium under dynamic loading conditions[J]. Journal of Applied Physics, 2009, 106(12):123532-1-123532-9. doi: 10.1063/1.3267325 [18] 种涛, 王桂吉, 谭福利, 等.磁驱动准等熵压缩下铁的相变[J].中国科学:物理学力学天文学, 2014, 44(6):630-636. http://www.cnki.com.cn/Article/CJFDTotal-JGXK201406011.htmZhong Tao, Wang Guiji, Tan Fuli, et al. Phase transition of iron under magnetically driven quasi-isentropic compression[J]. Scientia Sinica: Physica, Mechanica & Astronomica, 2014, 44(6):630-636. http://www.cnki.com.cn/Article/CJFDTotal-JGXK201406011.htm [19] 孙承纬. 磁驱动等熵压缩和高速飞片技术, 动高压原理与技术[M]. 经福谦, 陈俊祥, 主编. 北京: 国防工业出版社, 2006: 221 [20] Hayes D. Backward integration of the equation of motion to correct for free surface perturbations: SAND2001-1440[R]. Albuquerque, New Mexico, USA: Sandia National Laboratorys, 2001. https://www.osti.gov/biblio/783087/ [21] 张红平, 孙承纬, 李牧, 等.准等熵实验数据处理的反积分方法研究[J].力学学报, 2011, 43(1):105-111. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201100042122Zhang Hongping, Sun Cheng wei, Li Mu, et al. Backward integration method in data processing of quasi-isentropic compression experiments[J]. Chinese Journal of Theoretical and Applied Mechanics, 2011, 43(1):105-111. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201100042122 [22] Xue Quanxi, Wang Zhebin, Jiang Shaoen, et al. Characteristic method for isentropic compression simulation[J].AIP Advances, 2014, 4(5):057127. doi: 10.1063/1.4880039 [23] Davis J P. CHARICE 1. 0: An IDL application for characteristics-based inverse analysis of isentropic compression experiments: SAND2007-4984[R]. Albuquerque, New Mexico, USA: Sandia National Laboratorys, 2007. [24] Vogler T J, Ao T, Asay J R. High-pressure strength of aluminum under quasi-isentropic loading[J]. International Journal of Plasticity, 2009, 25(25):679-694. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=c0f72dcaf714e7ca22c8da7485c488a3 [25] Brown J L, Alexander C S, Asay J R, et al. Extracting strength from high pressure ramp-release experiments[J]. Journal of Applied Physics, 2013, 114(22):223518-223518-16. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ0231766782/ [26] Brown J L, Knudson M D, Alexander C S, et al. Shockless compression and release behavior of beryllium to 110 GPa[J]. Journal of Applied Physics, 2014, 116(3):033502. doi: 10.1063/1.4890232 期刊类型引用(1)
1. 徐鸿飞,王放,武郁文,翁春生. 基于小尺寸模型分析氢氧旋转爆轰波传播的不稳定性机制. 爆炸与冲击. 2025(01): 30-42 . 本站查看
其他类型引用(0)
-