Parameter inversion of rock RHT constitutive model using PAWN global sensitivity analysis and intelligent optimization algorithm

TIAN Haofan; SHAO Zekai; YU Ji; YOU Shuai; WANG Zhengzheng

doi:10.11883/bzycj-2025-0254

Article Contents

Article Navigation > Explosion And Shock Waves > 2026 > Accepted Manuscript

TIAN Haofan, SHAO Zekai, YU Ji, YOU Shuai, WANG Zhengzheng. Parameter inversion of rock RHT constitutive model using PAWN global sensitivity analysis and intelligent optimization algorithm[J]. Explosion And Shock Waves. doi: 10.11883/bzycj-2025-0254

Citation:

TIAN Haofan, SHAO Zekai, YU Ji, YOU Shuai, WANG Zhengzheng. Parameter inversion of rock RHT constitutive model using PAWN global sensitivity analysis and intelligent optimization algorithm[J]. Explosion And Shock Waves. doi: 10.11883/bzycj-2025-0254

Citation:

TIAN Haofan, SHAO Zekai, YU Ji, YOU Shuai, WANG Zhengzheng. Parameter inversion of rock RHT constitutive model using PAWN global sensitivity analysis and intelligent optimization algorithm[J]. Explosion And Shock Waves. doi: 10.11883/bzycj-2025-0254

PDF( 8731 KB)

Parameter inversion of rock RHT constitutive model using PAWN global sensitivity analysis and intelligent optimization algorithm

doi: 10.11883/bzycj-2025-0254

1.
School of Infrastructure Engineering, Dalian University of Technology, Dalian 116024, Liaoning, China
2.
School of Civil and Resource Engineering, University of Science and Technology Beijing, Beijing 100083, China

Received Date: 2025-08-05
Rev Recd Date: 2025-11-09

Available Online: 2025-11-13

Abstract

Abstract

The Riedel-Hiermaier-Thoma (RHT) constitutive model has been widely applied in tunnel blasting, impact-resistant structural design, and underground protective engineering due to its strong capability to describe the mechanical behavior of brittle materials such as rock and concrete under high-strain-rate and high-pressure conditions. However, the model involves a large number of nonlinear parameters, some of which are difficult to determine experimentally because of the high cost of testing. These key parameters are often adjusted through trial-and-error methods, which limit both modeling efficiency and simulation accuracy. To overcome these limitations, a comprehensive parameter inversion framework was developed for 16 difficult-to-calibrate parameters of the RHT model. The framework integrated the PAWN (Pianosi-Wagener) global sensitivity analysis method with intelligent optimization algorithms and coupled MATLAB with the ANSYS/LS-DYNA simulation platform. The area difference of the stress-strain curve was introduced as the core evaluation metric, and a batch result-extraction and automated three-wave alignment technique was developed. Based on these developments, an efficient and reliable RHT parameter inversion system was established, achieving, for the first time, a global sensitivity analysis (GSA) and automated inversion of key parameters in the RHT model. The results show that, among the 16 parameters analyzed, only eight exert a significant influence on the model response. The intelligent optimization–based inversion achieved relative errors ranging from 0.23% to 9.28%, and the reliability of the calibrated parameters was verified through Semicircular Bend Split Hopkinson Pressure Bar (SCB-SHPB) tests and scaled blasting experiments. The proposed method significantly enhances both the efficiency and accuracy of RHT parameter calibration without the need to construct large sample datasets, and it is applicable to a wide range of loading conditions. Compared with traditional calibration approaches, the required inversion accuracy was achieved in fewer than 15 iterations, meeting the dual demands of computational efficiency and precision. Overall, the proposed framework provides a new and effective approach for sensitivity analysis and parameter calibration of dynamic constitutive models, demonstrating strong engineering applicability and practical value.
- calibration of RHT constitutive parameters,
- PAWN method,
- intelligent optimization algorithm,
- global sensitivity,
- parameter inversion

FullText(HTML)

References(53)

References

[1]	CHEN C, ZHANG X H, HAO H, et al. Discussion on the suitability of dynamic constitutive models for prediction of geopolymer concrete structural responses under blast and impact loading [J]. International Journal of Impact Engineering, 2022, 160: 104064. DOI: 10.1016/j.ijimpeng.2021.104064.
[2]	YAN L, CHEN L, CHEN B Y, et al. Analysis and evaluation of suitability of high-pressure dynamic constitutive model for concrete under blast and impact loading [J]. International Journal of Impact Engineering, 2025, 195: 105145. DOI: 10.1016/j.ijimpeng.2024.105145.
[3]	RIEDEL W, THOMA K, HIERMAIER S, et al. Penetration of reinforced concrete by BETA-B-500 numerical analysis using a new macroscopic concrete model for hydrocodes [C]//Proceedings of the 9th International Symposium on the Effects of Munitions with Structures: Vol. 315. Berlin-Strausberg Germany, 1999: 315–322.
[4]	BANADAKI M M D. Stress-wave induced fracture in rock due to explosive action [D]. Toronto: University of Toronto, 2010.
[5]	XIE L X, LU W B, ZHANG Q B, et al. Damage evolution mechanisms of rock in deep tunnels induced by cut blasting [J]. Tunnelling and Underground Space Technology, 2016, 58: 257–270. DOI: 10.1016/j.tust.2016.06.004.
[6]	PAN C, XIE L X, LI X, et al. Numerical investigation of effect of eccentric decoupled charge structure on blasting-induced rock damage [J]. Journal of Central South University, 2022, 29(2): 663–679. DOI: 10.1007/s11771-022-4947-3.
[7]	TIAN H F, BAO T, LI Z, et al. Determination of constitutive parameters of crystalline limestone based on improved RHT model [J]. Advances in Materials Science and Engineering, 2022, 2022: 3794898. DOI: 10.1155/2022/3794898.
[8]	LI S L, LING T L, LIU D S, et al. Determination of rock mass parameters for the RHT model based on the Hoek–Brown criterion [J]. Rock Mechanics and Rock Engineering, 2023, 56(4): 2861–2877. DOI: 10.1007/s00603-022-03189-9.
[9]	WANG H C, WANG Z L, WANG J G, et al. Effect of confining pressure on damage accumulation of rock under repeated blast loading [J]. International Journal of Impact Engineering, 2021, 156: 103961. DOI: 10.1016/j.ijimpeng.2021.103961.
[10]	NI Y, WANG Z L, LI S Y, et al. Numerical study on the dynamic fragmentation of rock under cyclic blasting and different in-situ stresses [J]. Computers and Geotechnics, 2024, 172: 106404. DOI: 10.1016/j.compgeo.2024.106404.
[11]	李洪超, 刘殿书, 赵磊, 等. 大理岩RHT模型参数确定研究 [J]. 北京理工大学学报, 2017, 37(8): 801–806. DOI: 10.15918/j.tbit1001-0645.2017.08.006. LI H C, LIU D S, ZHAO L, et al. Study on parameters determination of marble RHT model [J]. Transactions of Beijing Institute of Technology, 2017, 37(8): 801–806. DOI: 10.15918/j.tbit1001-0645.2017.08.006.
[12]	LEVASSEUR S, MALÉCOT Y, BOULON M, et al. Soil parameter identification using a genetic algorithm [J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2008, 32(2): 189–213. DOI: 10.1002/nag.614.
[13]	ROKONUZZAMAN M, SAKAI T. Calibration of the parameters for a hardening–softening constitutive model using genetic algorithms [J]. Computers and Geotechnics, 2010, 37(4): 573–579. DOI: 10.1016/j.compgeo.2010.02.007.
[14]	PARK D, PARK E S. Inverse parameter fitting of tunnels using a response surface approach [J]. International Journal of Rock Mechanics and Mining Sciences, 2015, 77: 11–18. DOI: 10.1016/j.ijrmms.2015.03.026.
[15]	SUN J L, WU S C, WANG H, et al. Inversion of surrounding rock mechanical parameters in a soft rock tunnel based on a hybrid model EO-LightGBM [J]. Rock Mechanics and Rock Engineering, 2023, 56(9): 6691–6707. DOI: 10.1007/s00603-023-03387-z.
[16]	CUI J Q, WU S C, CHENG H Y, et al. Composite interpretability optimization ensemble learning inversion surrounding rock mechanical parameters and support optimization in soft rock tunnels [J]. Computers and Geotechnics, 2024, 165: 105877. DOI: 10.1016/j.compgeo.2023.105877.
[17]	ZHOU C J, GAO B, YAN B, et al. A combined machine learning/search algorithm-based method for the identification of constitutive parameters from laboratory tests and in-situ tests [J]. Computers and Geotechnics, 2024, 170: 106268. DOI: 10.1016/j.compgeo.2024.106268.
[18]	ÖZEL T, KARPAT Y. Identification of constitutive material model parameters for high-strain rate metal cutting conditions using evolutionary computational algorithms [J]. Materials and Manufacturing Processes, 2007, 22(5): 659–667. DOI: 10.1080/10426910701323631.
[19]	ANDRADE-CAMPOS A, THUILLIER S, PILVIN P, et al. On the determination of material parameters for internal variable thermoelastic–viscoplastic constitutive models [J]. International Journal of Plasticity, 2007, 23(8): 1349–1379. DOI: 10.1016/j.ijplas.2006.09.002.
[20]	ZHAI Y, ZHAO R F, LI Y B, et al. Stochastic inversion method for dynamic constitutive model of rock materials based on improved DREAM [J]. International Journal of Impact Engineering, 2021, 147: 103739. DOI: 10.1016/j.ijimpeng.2020.103739.
[21]	MENG K, YE G G. An inversion study of constitutive parameters for powder liner and hard rock based on finite element simulation [J]. Applied Sciences, 2025, 15(6): 3065. DOI: 10.3390/app15063065.
[22]	宋力, 胡时胜. SHPB数据处理中的二波法与三波法 [J]. 爆炸与冲击, 2005, 25(4): 368–373. DOI: 10.11883/1001-1455(2005)04-0368-06. SONG L, HU S S. Two-wave and three-wave method in SHPB data processing [J]. Explosion and Shock Waves, 2005, 25(4): 368–373. DOI: 10.11883/1001-1455(2005)04-0368-06.
[23]	WANG Z L, HUANG Y P, LI S Y, et al. SPH-FEM coupling simulation of rock blast damage based on the determination and optimization of the RHT model parameters [J]. IOP Conference Series: Earth and Environmental Science, 2020, 570(4): 042035. DOI: 10.1088/1755-1315/570/4/042035.
[24]	刘殿柱, 刘娜, 高天赐, 等. 应用正交试验法的RHT模型参数敏感性研究 [J]. 北京理工大学学报, 2019, 39(6): 558–564. DOI: 10.15918/j.tbit1001-0645.2019.06.002. LIU D Z, LIU N, GAO T C, et al. Study on the parameter sensitivity of RHT concrete model by orthogonal test technique [J]. Transactions of Beijing Institute of Technology, 2019, 39(6): 558–564. DOI: 10.15918/j.tbit1001-0645.2019.06.002.
[25]	ROUHANI H, FARROKH E. Failure analysis of Nehbandan granite under various stress states and strain rates using a calibrated Riedel–Hiermaier–Thoma constitutive model [J]. Geomechanics and Geophysics for Geo-Energy and Geo-Resources, 2024, 10(1): 157. DOI: 10.1007/s40948-024-00876-5.
[26]	马凯, 任高峰, 葛永翔, 等. 硬石膏矿RHT动力学模型参数研究 [J]. 爆破, 2024, 41(4): 35–44. DOI: 10.3963/j.issn.1001-487X.2024.04.005. MA K, REN G F, GE Y X, et al. Study on parameters of anhydrite RHT dynamic model [J]. Blasting, 2024, 41(4): 35–44. DOI: 10.3963/j.issn.1001-487X.2024.04.005.
[27]	李洪超, 陈勇, 刘殿书, 等. 岩石RHT模型主要参数敏感性及确定方法研究 [J]. 北京理工大学学报, 2018, 38(8): 779–785. DOI: 10.15918/j.tbit1001-0645.2018.08.002. LI H C, CHEN Y, LIU D S, et al. Sensitivity analysis determination and optimization of rock RHT parameters [J]. Transactions of Beijing Institute of Technology, 2018, 38(8): 779–785. DOI: 10.15918/j.tbit1001-0645.2018.08.002.
[28]	LIU W L, DING L Y. Global sensitivity analysis of influential parameters for excavation stability of metro tunnel [J]. Automation in Construction, 2020, 113: 103080. DOI: 10.1016/j.autcon.2020.103080.
[29]	WANG Y K, HUANG J S, TANG H M. Global sensitivity analysis of the hydraulic parameters of the reservoir colluvial landslides in the Three Gorges Reservoir area, China [J]. Landslides, 2020, 17(2): 483–494. DOI: 10.1007/s10346-019-01290-9.
[30]	ZHOU S T, ZHANG Z X, LUO X D, et al. Predicting dynamic compressive strength of frozen-thawed rocks by characteristic impedance and data-driven methods [J]. Journal of Rock Mechanics and Geotechnical Engineering, 2024, 16(7): 2591–2606. DOI: 10.1016/j.jrmge.2023.09.017.
[31]	SALTELLI A, ANNONI P, AZZINI I, et al. Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index [J]. Computer Physics Communications, 2010, 181(2): 259–270. DOI: 10.1016/j.cpc.2009.09.018.
[32]	PIANOSI F, WAGENER T. A simple and efficient method for global sensitivity analysis based on cumulative distribution functions [J]. Environmental Modelling and Software, 2015, 67: 1–11. DOI: 10.1016/j.envsoft.2015.01.004.
[33]	HERRMANN W. Constitutive equation for the dynamic compaction of ductile porous materials [J]. Journal of Applied Physics, 1969, 40(6): 2490–2499. DOI: 10.1063/1.1658021.
[34]	MEYERS M A. Dynamic behavior of materials [M]. New York: John Wiley & Sons, Inc. , 1994. DOI: 10.1002/9780470172278.
[35]	辛春亮, 薛再清, 涂建, 等. 有限元分析常用材料参数手册 [M]. 北京: 机械工业出版社, 2020.
[36]	汪大为, 王志亮, 汪书敏. 基于RHT本构的隧道爆破参数优化与振动损伤数值研究 [J]. 水文地质工程地质, 2023, 50(6): 129–136. DOI: 10.16030/j.cnki.issn.1000-3665.202210009. WANG D W, WANG Z L, WANG S M. Numerical study of tunnel blasting parameter optimization and vibration damage based on the RHT constitutive model [J]. Hydrogeology & Engineering Geology, 2023, 50(6): 129–136. DOI: 10.16030/j.cnki.issn.1000-3665.202210009.
[37]	PATTAJOSHI S, RAY S. Dynamic fracture analysis of multi-layer composites using an improved RHT model [J]. Engineering Fracture Mechanics, 2025, 325: 111220. DOI: 10.1016/j.engfracmech.2025.111220.
[38]	SONG S, XU X Y, REN W J, et al. Determination and application of the RHT constitutive model parameters for ultra-high-performance concrete [J]. Structures, 2024, 69: 107488. DOI: 10.1016/j.istruc.2024.107488.
[39]	WANG Z L, NI Y, WANG J G, et al. Improvement and performance analysis of constitutive model for rock blasting damage simulation [J]. Simulation Modelling Practice and Theory, 2025, 138: 103043. DOI: 10.1016/j.simpat.2024.103043.
[40]	WANG Z X, HUANG J H, CHEN Y L, et al. Dynamic mechanical properties of different types of rocks under impact loading [J]. Scientific Reports, 2023, 13(1): 19147. DOI: 10.1038/s41598-023-46444-x.
[41]	WANG Z L, WANG H C, WANG J G, et al. Finite element analyses of constitutive models performance in the simulation of blast-induced rock cracks [J]. Computers and Geotechnics, 2021, 135: 104172. DOI: 10.1016/j.compgeo.2021.104172.
[42]	邱薛, 刘晓辉, 胡安奎, 等. 煤岩动态RHT本构模型数值模拟研究 [J]. 煤炭学报, 2024, 49(S1): 261–273. DOI: 10.13225/j.cnki.jccs.2023.0540. QIU X, LIU X H, HU A K, et al. Research on numerical simulation of coal dynamic RHT constitutive model [J]. Journal of China Coal Society, 2024, 49(S1): 261–273. DOI: 10.13225/j.cnki.jccs.2023.0540.
[43]	ZHANG J W, CHEN Z W, SHAO K. Numerical investigation on optimal blasting parameters of tunnel face in granite rock [J]. Simulation Modelling Practice and Theory, 2024, 130: 102854. DOI: 10.1016/j.simpat.2023.102854.
[44]	ALTMAN D G, BLAND J M. Statistics notes: quartiles, quintiles, centiles, and other quantiles [J]. BMJ, 1994, 309(6960): 996–996. DOI: 10.1136/bmj.309.6960.996.
[45]	LUO D H, WAN X, LIU J M, et al. Optimally estimating the sample mean from the sample size, median, mid-range, and/or mid-quartile range [J]. Statistical Methods in Medical Research, 2018, 27(6): 1785–1805. DOI: 10.1177/0962280216669183.
[46]	聂铮玥, 彭永, 陈荣, 等. 侵彻条件下岩石类材料RHT模型参数敏感性分析 [J]. 振动与冲击, 2021, 40(14): 108–116. DOI: 10.13465/j.cnki.jvs.2021.14.015. NIE Z Y, PENG Y, CHEN R, et al. Sensitivity analysis of RHT model parameters for rock materials under penetrating condition [J]. Journal of Vibration and Shock, 2021, 40(14): 108–116. DOI: 10.13465/j.cnki.jvs.2021.14.015.
[47]	XUE J K, SHEN B. A novel swarm intelligence optimization approach: sparrow search algorithm [J]. Systems Science & Control Engineering, 2020, 8(1): 22–34. DOI: 10.1080/21642583.2019.1708830.
[48]	MIRJALILI S, MIRJALILI S M, LEWIS A. Grey wolf optimizer [J]. Advances in Engineering Software, 2014, 69: 46–61. DOI: 10.1016/j.advengsoft.2013.12.007.
[49]	MIRJALILI S, LEWIS A. The whale optimization algorithm [J]. Advances in Engineering Software, 2016, 95: 51–67. DOI: 10.1016/j.advengsoft.2016.01.008.
[50]	ARORA S, SINGH S. Butterfly optimization algorithm: a novel approach for global optimization [J]. Soft Computing, 2019, 23(3): 715–734. DOI: 10.1007/s00500-018-3102-4.
[51]	HEIDARI A A, MIRJALILI S, FARIS H, et al. Harris hawks optimization: algorithm and applications [J]. Future Generation Computer Systems, 2019, 97: 849–872. DOI: 10.1016/j.future.2019.02.028.
[52]	LIAN J B, HUI G H, MA L, et al. Parrot optimizer: algorithm and applications to medical problems [J]. Computers in Biology and Medicine, 2024, 172: 108064. DOI: 10.1016/j.compbiomed.2024.108064.
[53]	DAI F, XIA K W, TANG L Z. Rate dependence of the flexural tensile strength of Laurentian granite [J]. International Journal of Rock Mechanics and Mining Sciences, 2010, 47(3): 469–475. DOI: 10.1016/j.ijrmms.2009.05.001.