Bursting Performance Optimization of Reverse-Arched Bursting Discs Based on Variable Fidelity Surrogate Models
-
摘要: 根据反拱爆破片(reverse-arched bursting discs,RABDs)爆破性能的高、低精度有限元计算结果,通过构建分层克里金(H-Kriging)代理模型,实现了反拱爆破片爆破压力的快速预报,建立了反拱爆破片爆破性能的数学模型,基于此进行反拱爆破片的结构优化设计。研究结果表明,基于高、低精度有限元模型所构建的爆破压力-结构参数的分层克里金代理模型,能在显著节约计算成本的前提下,准确预报反拱爆破片的爆破压力。针对反拱爆破片结构的初始设计方案,采用遗传算法进行优化设计,优化方案能够在考虑爆破片厚度加工公差的情形下,使爆破压力的波动幅度降低58.8%,从而大大降低反拱爆破片爆破压力对厚度加工误差的敏感程度,具有较高的工程参考价值。Abstract: To address the optimization design problem of the bursting performance of reverse-arched bursting discs (RABDs), a hierarchical Kriging (H-Kriging) surrogate model was constructed based on both high- and low-fidelity finite element analysis results. This model enables the rapid prediction of the burst pressure of RABDs, facilitating the development of a mathematical model for performance optimization and structural improvement. The results show that the H-Kriging surrogate model relating burst pressure to structural parameters based on high- and low-fidelity finite element models can significantly reduce computational cost while accurately predicting the burst pressure of RABDs. For the initial structural design scheme of RABDs, optimization was carried out using a genetic algorithm, with the optimized design accounting for manufacturing tolerance in disc thickness. This resulted in a 58.8% reduction in burst pressure fluctuation, significantly reducing the sensitivity of burst pressure to thickness manufacturing errors and providing valuable engineering reference.
-
表 1 反拱爆破片数值模拟计算参数
Table 1. Parameters of RABDs in numerical simulations
A/MPa B/MPa n c m 340.68 116 0.61 0.01 0.517 表 3 助爆环锥爆破片结构的无量纲几何尺寸
Table 3. Dimensionless geometric dimensions of bursting discs with aid-bursting ring cones
Sample H/rD L1/rD d1/rD d2/rD b/rD L/rD 1 0.0125 0.31 1.5 1.25 0.01 0.0375 表 4 设计变量的样本点取值
Table 4. Sample values of design variables
L1/mm H/mm d1/mm R/mm 1.8 0.4 6.0 0.5 1.9 0.5 8.0 0.8 2.0 0.6 9.3 1.0 2.1 0.7 10.0 1.3 2.2 0.8 11.5 1.6 12.6 2.0 Note:considering that the variations in arch height and thickness are relatively small, fewer sample points were selected. 表 5 交互项影响分析
Table 5. Analysis of interaction effects
Interaction term P Result L1:H 0.203 8 Not statistically significant L1:d1 0.673 7 Not statistically significant H:d1 0.000 3 Statistically significant effect L1:R 0.069 3 Marginally significant effect H:R 0.016 2 Statistically significant effect d1:R 0.010 7 Statistically significant effect 表 6 优化目标和约束条件
Table 6. Optimization objectives and constraints
Proxy model $ p=f\left(H,\ L_1,\ d_1,\ R\right) $ Optimization objective $ \left[f\left(H+\varDelta,\ L_1,\ d_1,\ R\right)-f\left(H-\varDelta,\ L_1,\ d_1,\ R\right)\right] $ Constraint conditions $ p_0-0.5 \lt p \lt p_0+0.5 $ $ 0.04< H< 0.08 $ $ 1.8 \lt L_1 \lt 2.2 $ $ 6 \lt d_1 \lt 12.6 $ $ 0.5< R< 2 $ Note: $ {p}_{0} $=1.5 MPa,ideal blasting pressure;$ \varDelta=0.02\mathrm{\ mm} $,processing tolerance; the allowable fluctuation range of the blasting pressure is $ {p}_{0}-0.5< p< {p}_{0}+0.5 $; the units of H, L1, d1, R are mm. 表 7 不同设计方案的有限元分析结果
Table 7. Finite element results of different design schemes
Scheme Structure parameters/mm Blasting pressure/MPa Error/% L1 H d1 R High-precision finite element Proxy model P0 2.0 $ 0.06-\varDelta $ 9.3 1 0.64 0.61 −4.7 $ 0.06+\varDelta $ 2.09 2.09 0 P1 2.2 $ 0.06-\varDelta $ 9.3 2 0.94 1.00 6.3 $ 0.06+\varDelta $ 1.71 1.61 −5.8 表 8 优化结果
Table 8. Optimization results
Scheme Structure parameters/mm Target variable Optimization effect/% L1 H d1 R Fluctuation of blasting pressure/MPa P0 2.0 0.06 9.3 1 1.48 P1 2.2 0.06 9.3 2 0.61 58.8 -
[1] API. API Standard 521 Pressure-relieving and depressuring systems [S]. Washington: American Petroleum Institute, 2014. [2] ISO 4126 Safety devices for protection against excessive pressure [S]. 2003. [3] 李志义, 喻健良. 爆破片技术及应用 [M]. 北京: 化学工业出版社, 2006.LI Z Y, YU J L. Bursting disc technology and application [M]. Beijing: Chemical Industry Press, 2006. [4] 刘中华, 苏晨, 汪军安, 等. 气路膜片设计研究 [J]. 火箭推进, 2015, 41(5): 95–100. doi: 10.3969/j.issn.1672-9374.2015.05.018LIU Z H, SU C, WANG J A, et al. Design and study of pneumatic diaphragm in gas circuit [J]. Journal of Rocket Propulsion, 2015, 41(5): 95–100. doi: 10.3969/j.issn.1672-9374.2015.05.018 [5] 袁晨星. 反拱型爆破片结构参数对爆破压力及泄放特性影响研究 [D]. 徐州: 中国矿业大学, 2022.YUAN C X. Influence of structural parameters of reverse arched rupture disc on bursting pressure and relief characteristics [D]. Xuzhou: China University of Mining and Technology, 2022. [6] 李岳, 闫兴清, 喻健良, 等. 升压速率对反拱带槽型爆破片爆破压力影响的试验研究 [J]. 压力容器, 2020, 37(1): 8–13, 19. doi: 10.3969/j.issn.1001-4837.2020.01.002LI Y, YAN X Q, YU J L, et al. An experimental study for the effects of pressure rise rate on bursting pressure of reverse domed scored bursting disc [J]. Pressure Vessel Technology, 2020, 37(1): 8–13, 19. doi: 10.3969/j.issn.1001-4837.2020.01.002 [7] 李岳, 喻健良, 姜一昌, 等. 夹持条件对反拱形爆破片失稳压力影响研究 [J]. 压力容器, 2016, 33(7): 53–57. doi: 10.3969/j.issn.1001-4837.2016.07.009LI Y, YU J L, JIANG Y C, et al. Study of effect of clamping condition of reverse domed bursting disc on its buckling pressure [J]. Pressure Vessel Technology, 2016, 33(7): 53–57. doi: 10.3969/j.issn.1001-4837.2016.07.009 [8] 闫照锋. 金属爆破片爆破性能及影响因素数值研究 [D]. 大连: 大连理工大学, 2012.YAN Z F. Numerical study on explosion performance and influential factors of the metallic rupture discs [D]. Dalian: Dalian University of Technology, 2012. [9] YANG H T Y, HEINSTEIN M, SHIH J M. Adaptive 2D finite element simulation of metal forming processes [J]. International Journal for Numerical Methods in Engineering, 1989, 28(6): 1409–1428. doi: 10.1002/nme.1620280612 [10] AHMED M, HASHMI M S J. Comparison of free and restrained bulge forming by finite element method simulation [J]. Journal of Materials Processing Technology, 1997, 63(1/2/3): 651–654. doi: 10.1016/S0924-0136(96)02701-X [11] 张龙, 尉成果, 刘秉斌. 基于遗传算法的大型矩形压力容器刚强度优化设计 [J]. 压力容器, 2024, 41(2): 58–65. doi: 10.3969/j.issn.1001-4837.2024.02.008ZHANG L, WEI C G, LIU B B. Stiffness and strength optimization design of large rectangular pressure vessels based on genetic algorithm [J]. Pressure Vessel Technology, 2024, 41(2): 58–65. doi: 10.3969/j.issn.1001-4837.2024.02.008 [12] MURTY D V R. Finite element analysis of rupture disc [J]. Scientist EI, Design & Engineering Division, Indian Institute of Chemical Technology. [13] ISEKI H, JIMMA T, MUROTA T. Finite element method of analysis of the hydrostatic bulging of a sheet metal (part 1) [J]. Bulletin of JSME, 1974, 17(112): 1240–1246. doi: 10.1299/jsme1958.17.1240 [14] 詹晶. 普通反拱形爆破片爆破压力影响因素研究 [D]. 上海: 华东理工大学, 2017.ZHAN J. Burst pressure research of influence factors on reverse arch rupture disc [D]. Shanghai: East China University of Science and Technology, 2017. [15] 王丛飞, 许光, 满满, 等. 反拱带槽型爆破片结构参数研究 [J]. 真空与低温, 2021, 27(5): 485–489. doi: 10.3969/j.issn.1006-7086.2021.05.012WANG C F, XU G, MAN M, et al. Research on parameter characteristics of reverse domed scored bursting disc [J]. Vacuum and Cryogenics, 2021, 27(5): 485–489. doi: 10.3969/j.issn.1006-7086.2021.05.012 [16] 钱志强, 章碧成, 施威, 等. 反拱形爆破片翻转行为数值模拟研究 [J]. 压力容器, 2022, 39(12): 30–37. doi: 10.3969/j.issn.1001-4837.2022.12.005QIANG Z Q, ZHANG B C, SHI W, et al. Numerical simulation of overturning behavior of reverse domed bursting disc [J]. Pressure Vessel Technology, 2022, 39(12): 30–37. doi: 10.3969/j.issn.1001-4837.2022.12.005 [17] ZHANG P, CHENG Y S, LIU J, et al. Experimental and numerical investigations on laser-welded corrugated-core sandwich panels subjected to air blast loading [J]. Marine Structures, 2015, 40: 225–246. doi: 10.1016/j.marstruc.2014.11.007 [18] 鲁寨军, 冯书颖, 刘杰夫, 等. 十字开槽正拱型爆破片爆破性能影响因素研究 [J]. 工程爆破, 2023, 29(4): 114–121. doi: 10.19931/j.EB.20220279LU Z J, FENG S Y, LIU J F, et al. Study on the factors influencing the bursting performance of domed cross-scored rupture discs [J]. Engineering Blasting, 2023, 29(4): 114–121. doi: 10.19931/j.EB.20220279 -

下载: