混凝土高压状态方程实验与数值模拟研究

孙玉祥; 王杰; 武海军; 周婕群; 李金柱; 皮爱国; 黄风雷

doi:10.11883/bzycj-2020-0002

混凝土高压状态方程实验与数值模拟研究

doi: 10.11883/bzycj-2020-0002

孙玉祥^1,,
王杰^{1, 2},
武海军^1, ,,
周婕群³,
李金柱¹,
皮爱国¹,
黄风雷¹

1.
北京理工大学爆炸科学与技术国家重点实验室，北京 100081
2.
中国工程物理研究院化工材料研究所，四川绵阳 621999
3.
中国工程物理研究院总体工程研究所，四川绵阳 621999

基金项目: 国家自然科学基金联合基金（U1730128）；国家自然科学基金（11390362，11572048）

详细信息

作者简介:
孙玉祥（1993－　），男，硕士，1762566576@qq.com

通讯作者:
武海军（1974－　），男，博士，教授，wuhj@bit.edu.cn

中图分类号: O383
计量
- 文章访问数: 1635
- HTML全文浏览量: 1283
- PDF下载量: 185
- 被引次数: 0
出版历程
- 收稿日期: 2020-01-02
- 修回日期: 2020-06-09
- 刊出日期: 2020-12-05

Experiment and simulation on high-pressure equation of state for concrete

1.
State Key Laboratory of Explosion Science and Technology, Beijing Institute of Technology, Beijing 100081, China
2.
Institute of Chemical Materials, China Academy of Engineering Physics, Mianyang 621999, Sichuan, China
3.
Institute of Systems Engineering, China Academy of Engineering Physics, Mianyang 621999, Sichuan, China

摘要

摘要: 为研究高压下混凝土冲击压缩特性以及确定HJC本构模型状态方程参数，采用$\varnothing $58 mm火炮加载技术和多普勒探针系统（Doppler probe system，DPS）测速技术，对抗压强度为26.5 MPa（C25）和42.1 MPa（C40）的混凝土进行反向平板撞击实验研究与数值模拟。DPS探针记录TU1无氧铜靶自由面粒子速度历史，根据一维应变弹塑性波理论，计算撞击压力，拟合得到2～11 GPa高压条件下混凝土波速与粒子速度（u_s-u_p）、压力与体积应变（p-μ）关系。实验结果表明：高压条件下，混凝土波速-粒子速度呈线性关系；两种初始密度、孔隙率相近，强度不同的混凝土波速-粒子速度、压力-体积应变关系存在明显差异，相同压力下，混凝土试件强度越高，体积应变越小。基于实验结果，确定了两种强度混凝土HJC本构模型状态方程参数，利用LS-DYNA动力有限元分析软件对平板撞击实验进行了数值模拟，靶板自由面粒子速度历史与实验曲线吻合较好，仿真结果表明混凝土中冲击波的追赶卸载现象仅存在于低速撞击条件下。
- 混凝土 /
- 平板撞击 /
- 状态方程 /
- HJC /
- 无氧铜
Abstract: To study the dynamic compression characteristics of concrete under high hydrostatic pressure and to determine the equation of state parameters of the HJC constitutive model, inverse flyer-impact tests and numerical simulation analysis were conducted with two kinds of concrete flyers of which the compressive strengths were 26.5 MPa and 42.1 MPa, respectively. The concrete flyers were launched by $\varnothing $58 mm gun against TU1 copper targets. The particle velocity histories of the TU1 copper target free surface were measured by DPS (Doppler probe system). Based on the one-dimensional strain shock wave theory, the impact pressure was calculated. The relationships of shock velocity vs. particle velocity and pressure vs. volume strain in the pressure range of 2−11 GPa were fitted. The results show that the relationship between shock velocity and particle velocity of concrete is linear. The relationships of shock velocity vs. particle velocity and pressure vs. volume strain for concretes with similar initial density and porosity but different compressive strengths are obviously different. Under the same pressure, the higher the compressive strength of concrete, the smaller the volume strain is. According to the test results, the equation of state parameters of the HJC constitutive model were determined and the plate-impact tests were simulated by LS-DYNA. The simulated particle velocity histories of the TU1 target free surface were in good agreement with the experimental results. The simulation results show that the phenomenon of chasing and unloading of shock waves in concrete only exists under low velocity impact conditions.
- concrete /
- flyer-impact /
- equation of state /
- HJC /
- oxygen free copper

HTML全文

图 1 飞片撞击实验装置示意图

Figure 1. Flyer-impact experimental device

下载: 全尺寸图片幻灯片

图 2 DPS测速装置与靶室布置

Figure 2. DPS and target room arrangement

下载: 全尺寸图片幻灯片

图 3 混凝土飞片试件

Figure 3. Concrete flyer specimen

下载: 全尺寸图片幻灯片

图 4 飞片与弹托

Figure 4. Flyer and sabot

下载: 全尺寸图片幻灯片

图 5 TU1无氧铜靶片试件

Figure 5. TU1 copper specimen

下载: 全尺寸图片幻灯片

图 6 冲击波传播X-t图

Figure 6. Shock wave propagation

下载: 全尺寸图片幻灯片

图 7 TU1靶板自由面粒子速度历史示意图

Figure 7. Particle velocity history curve on TU1 target

下载: 全尺寸图片幻灯片

图 8 6路DPS探针原始影像图

Figure 8. Images by six DPS probes

下载: 全尺寸图片幻灯片

图 9 速度为989 m/s的C40混凝土飞片撞击TU1靶板的自由面粒子速度历史

Figure 9. Particle velocity of TU1 free surface in C40 concrete flyer impact test at 989 m/s

下载: 全尺寸图片幻灯片

图 10 C25混凝土飞片撞击TU1靶板自由面粒子速度历史

Figure 10. Particle velocity of TU1 free surface in C25 concrete flyer impact tests

下载: 全尺寸图片幻灯片

图 11 C40混凝土飞片撞击TU1靶板自由面粒子速度历史

Figure 11. Particle velocity of TU1 free surface in C40 concrete flyer impact tests

下载: 全尺寸图片幻灯片

图 12 混凝土波速-粒子速度（u_s-u_p）关系拟合

Figure 12. Shock velocity and particle velocity linear fitting

下载: 全尺寸图片幻灯片

图 13 混凝土压力-密度

Figure 13. Pressure-density data for concrete

下载: 全尺寸图片幻灯片

图 14 混凝土压力-体积应变

Figure 14. Pressure-volumetric strain data for concrete

下载: 全尺寸图片幻灯片

图 15 HJC模型状态方程^[2]

Figure 15. EOS of the HJC model

下载: 全尺寸图片幻灯片

图 16 数值模拟模型

Figure 16. Numerical simulation model

下载: 全尺寸图片幻灯片

图 17 C25混凝土粒子速度数值模拟与实验结果对比

Figure 17. Comparison of simulation and test results for C25 concrete

下载: 全尺寸图片幻灯片

图 18 C40混凝土粒子速度数值模拟与实验结果对比

Figure 18. Comparison of simulation and test results for C40 concrete

下载: 全尺寸图片幻灯片

图 19 速度为479 m/s C25混凝土飞片压力时程曲线

Figure 19. Pressure-time curves of C25 flyer at 479 m/s

下载: 全尺寸图片幻灯片

图 20 速度为479 m/s C25混凝土飞片压力云图

Figure 20. Pressure cloud in C25 concrete flyer at 479 m/s

下载: 全尺寸图片幻灯片

图 21 速度为1017 m/s C25混凝土飞片压力时程曲线

Figure 21. Pressure-time curve of C25 flyer at 1017 m/s

下载: 全尺寸图片幻灯片

表 1 混凝土平板撞击实验数据

Table 1. Flyer-impact test data of concrete

编号	类别	v₀/(m·s⁻¹)	u_p1/(m·s⁻¹)	u_p2/(m·s⁻¹)	u_s2/(m·s⁻¹)	σ_H/MPa
C25v1017	C25	1017	343.9	845.0	3469.7	6450.4
C25v1259	C25	1259	441.0	1038.5	3682.8	8414.2
C25v1523	C25	1523	575.6	1235.2	4136.2	11240.0
C25v479	C25	479	112.6	422.7	2177.6	2025.1
C40v834	C40	834	276.3	695.8	3383.8	5121.2
C40v989	C40	989	355.5	811.3	3786.7	6681.6
C40v1279	C40	1279	476.1	1040.9	4037.4	9140.6
C40v1484	C40	1484	580.6	1193.7	4370.9	11348.0

下载: 导出CSV

表 2 混凝土HJC模型状态方程参数

Table 2. EOS parameters of the HJC model for concrete

混凝土	K₁/GPa	K₂/GPa	K₃/GPa	p_L/MPa	μ_L
C25	56.34	−363.96	1689.50	0.9	0.171
C40	137.37	−971.90	3483.43	1.2	0.187

下载: 导出CSV

表 3 混凝土HJC模型主要参数

Table 3. Parameters of HJC model for concrete

混凝土	ρ/(g·cm⁻³)	G/GPa	f_c/MPa	T/MPa	p_C/MPa	μ_C/10⁻³
C25	2.202	9.517	26.5	3.19	8.83	0.696
C40	2.177	11.79	42.1	4.02	14.03	0.893
注：G为剪切模量，f_c为单轴抗压强度，T为最大静水拉伸强度，p_C为压碎压力，μ_C为压碎体积应变。

下载: 导出CSV

参考文献(15)

[1]	高飞, 王明洋, 张先锋, 等. 水泥砂浆的平板撞击实验与高压状态方程研究 [J]. 振动与冲击, 2018, 37(12): 41–47. DOI: 10.13465/j.cnki.jvs.2018.12.007. GAO F, WANG M Y, ZHANG X F, et al. A study on planar impact and equation of state for cement mortar [J]. Journal of Vibration and Shock, 2018, 37(12): 41–47. DOI: 10.13465/j.cnki.jvs.2018.12.007.
[2]	HOLMQUIST T J, JOHNSON G R. A computational constitutive model for concrete subjected to larger strains, high strain rates and high pressure [C] // The 14th International Symposium on Ballistics. Quebec, Canada: American Defense Preparedness Association, 1995: 591−600.
[3]	GRADY D E. Impact compression properties of concrete [C] // The 6th International Symposium on Interaction of Nonnuclearmunitions with Structures. Panama City, Florida, 1993:172−175.
[4]	GEBBEKEN N, GREULICH S, PIETZSCH A. Hugoniot properties for concrete determined by full-scale detonation experiments and flyer-plate-impact tests [J]. International Journal of Impact Engineering, 2006, 32(12): 2017–2031. DOI: 10.1016/j.ijimpeng.2005.08.003.
[5]	GEBBEKEN N, GREULICH S, PIETZSCH A. Equation of state data for concrete determined by full-scale experiments and flyer-plate-impact tests [C] // European Conference on Computational Mechanics. Cracow Poland, 2001.
[6]	RIEDEL W, WICKLEIN M, THOMA K. Shock properties of conventional and high strength concrete: experimental and mesomechanical analysis [J]. International Journal of Impact Engineering, 2008, 35(3): 155–171. DOI: 10.1016/j.ijimpeng.2007.02.001.
[7]	RIEDEL W, THOMA K, HIERMAIER S. 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. Berlin, 1999: 315-322
[8]	王永刚, 张远平, 王礼立. C25混凝土冲击绝热关系和Grüneisen型状态方程的实验研究 [J]. 物理学报, 2008, 57(12): 7789–7793. DOI: 10.7498/aps.57.7789. WANG Y G, ZHANG Y P, WANG L L. Experimental study on the shock Hugoniot relationship and the Grüneisen-type equation of state for C25 concrete [J]. Acta Physica Sinica, 2008, 57(12): 7789–7793. DOI: 10.7498/aps.57.7789.
[9]	GRADY D E. Shock equation of state properties of concrete[R]. Office of Scientific & Technical Information Technical Reports, 1996. DOI: 10.2495/SUSI960371.
[10]	KIPP M E, CHHABILDAS L C, REINHART W D. Elastic shock response and spall strength of concrete[C]// AIP Conference. The tenth American Physical Society Topical Conference on Shock Compression of Condensed Matter. Amherst, Massachusetts, 1998: 557−560. DOI: 10.1063/1.55664.
[11]	HALL C A, CHHABILDAS L C, REINHART W D. Shock Hugoniot and release in concrete with different aggregate sizes from 3 to 23 GPa [J]. International Journal of Impact Engineering, 1999, 23(1): 341–351. DOI: 10.1016/S0734-743X(99)00085-8.
[12]	TSEMBELIS K, MILLETT J C F, PROUD W G, et al. The shock Hugoniot properties of cement paste up to 5GPa [J]. AIP Conference Proceedings, 2000, 505(1): 1267–1270. DOI: 10.1063/1.1303692.
[13]	TSEMBELIS K, PROUD W G, WILLMOTT G R, et al. The shock Hugoniot properties of cement paste & mortar up to 18 GPa [J]. AIP Conference Proceedings, 2004, 706(1): 1488–1491. DOI: 10.1063/1.1780520.
[14]	蒋国平, 浣石, 郝洪, 等. 钢纤维高强混凝土材料的气体炮试验研究 [J]. 物理学报, 2013(1): 351–356. DOI: CNKI:SUN:WLXB.0.2013-01-054. JIANG G P, HUAN S, HAO H, et al. Performance of steel reinforced high strength concrete investigated in the gas gun experiment [J]. Acta Physica Sinica, 2013(1): 351–356. DOI: CNKI:SUN:WLXB.0.2013-01-054.
[15]	马晓青. 冲击动力学[M]. 北京: 北京理工大学出版社, 1992: 168−169.