基于Python二次开发的三相混凝土侵彻模拟与参数分析

Three-phase Concrete Penetration Simulation and Parametric Analysis Based on Python Secondary Development

摘要

混凝土的细观结构对其抗侵彻性能影响显著，但界面过渡层厚度极小，难以在数值模型中精确表征。为此，基于Python语言对Abaqus进行二次开发，提出一种三相混凝土细观力学建模方法。该方法通过在骨料与砂浆基体之间插入零厚度Cohesive单元模拟界面过渡层，以表征其力学行为；同时编写考虑孔隙状态方程、强度面、拉压损伤演化及应变率效应的材料子程序。将细观模型与材料子程序结合，模拟卵头钢弹侵彻混凝土靶的过程，计算结果与实验数据吻合良好。基于该模型，进一步分析粗骨料强度、含量及界面过渡层对侵彻深度的影响。结果表明：当骨料体积分数由0%增至40%时，在650 m/s高速撞击下侵彻深度降低18.0%；骨料强度由100 MPa提升至200 MPa时，侵彻深度降低28.6%；当界面相对强度比值由20%提升至100%时，侵彻深度降低17.1%。研究揭示了界面过渡层、骨料含量与强度等细观参数对抗侵彻性能的调控机制，可为抗侵彻混凝土的细观设计提供理论依据。

Abstract

The meso-structure of concrete was found to have a significant influence on its anti-penetration performance. However, the interfacial transition zone (ITZ) is extremely thin, making it difficult to be accurately characterized in numerical models. To address this issue, a three-phase meso-mechanical modeling method for concrete was proposed based on the secondary development of Abaqus using the Python language. Zero-thickness cohesive elements were inserted between the aggregates and the mortar matrix to simulate the ITZ, through which its mechanical behavior was characterized. Meanwhile, a user-defined material subroutine was developed, in which a pore equation of state, a strength surface, tensile/compressive damage evolution laws, and a strain rate effect were taken into consideration. The proposed meso-scale model was then combined with the material subroutine to simulate the penetration process of an ogive-nose steel projectile into a concrete target. The results indicate that the depth of penetration (DOP) is reduced by 18.0% when the aggregate volume fraction is increased from 0% to 40% under an impact velocity of 650 m/s. A reduction of 28.6% in DOP is observed when the aggregate strength is increased from 100 MPa to 200 MPa. Furthermore, a reduction of 17.1% in DOP is obtained when the relative strength ratio of the ITZ is increased from 20% to 100%. The regulatory effects of meso-parameters—namely the ITZ, aggregate content, and aggregate strength—on the anti-penetration performance are thus revealed. This study provides a theoretical basis for the meso-structural design of anti-penetration concrete.

关键词

细观混凝土 / Python二次开发 / Cohesive单元 / 界面过渡层 / 材料子程序 / 卵头钢弹 / 混凝土侵彻

Key words

meso-concrete / Python-based secondary development / Cohesive element / interface transition layer / material subroutine / ogive-nosed steel / concrete penetration

参考文献

[1] WANG Y, ZHANG W, WANG J, et al. Effects of coarse aggregate size on thickness and micro—properties of ITZ and the mechanical properties of concrete[J]. Cement and Concrete Composites, 2024, 154: 105777.

[2] 韩晶, 王华, 陈智刚. 分体式动能子弹对混凝土靶的侵彻与毁伤特性分析[J]. 高压物理学报, 2016, 30(5): 399-405.

[3] HAN J, WANG H, CHEN Z G. Penetrating and damaging effects of a split kinetic energy projectile on concrete targets[J]. Chinese Journal of High Pressure Physics, 2016, 30(5): 399-405.

[4] WANG W, SONG X, YANG J, et al. Experimental and numerical research on the effect of ogive—nose projectile penetrating UR50 ultra—early—strength concrete[J]. Cement and Concrete Composites, 2023, 136: 104902.

[5] YAO Z, QIN F, LI J, et al. Experimental and 3D mesoscopic numerical simulation study of kinetic projectile penetrating into concrete[J]. International Journal of Impact Engineering, 2025, 195: 105140.

[6] LIU C, ZHANG X, CHEN H, et al. Experimental and theoretical study on steel Long—rod projectile penetration into concrete targets with elevated impact velocities[J]. International Journal of Impact Engineering, 2020, 138: 103482.

[7] DONG H, LIU Z H, WU H J, et al. Study on penetration characteristics of high—speed elliptical cross—sectional projectiles into concrete[J]. International Journal of Impact Engineering, 2019, 132: 103311.

[8] ZHANG J, WANG Z, WANG Z, et al. Impact mechanisms of aggregate on rigid projectile normal penetration into concrete target[J]. Scientia Sinica Technologica, 2021, 51(3): 272-280.

[9] LAI J, YANG H, WANG H, et al. Properties and modeling of ultra—high—performance concrete subjected to multiple bullet impacts[J]. Journal of Materials in Civil Engineering, 2018, 30(10): 04018256.

[10] WERNER S, THIENEL K C, KUSTERMANN A. Study of fractured surfaces of concrete caused by projectile impact[J]. International Journal of Impact Engineering, 2013, 52: 23-27.

[11] CHEN X G, LU F Y, ZHANG D. Penetration trajectory of concrete targets by ogived steel projectiles: experiments and simulations[J]. International Journal of Impact Engineering, 2018, 120: 202-213.

[12] TIAN F, CHEN B. An FDEM meso—scale fracture mechanism analysis of concrete under compression based on cohesive zone model[J]. Structures, 2025, 82: 110822.

[13] SU H, LI H, HU B, et al. A research on the macroscopic and mesoscopic parameters of concrete based on an experimental design method[J]. Materials, 2021, 14(7): 1627.

[14] LIU J, WU C, LI J, et al. Experimental and numerical study of reactive powder concrete reinforced with steel wire mesh against projectile penetration[J]. International Journal of Impact Engineering, 2017, 109: 131-149.

[15] PENG Y, WU H, FANG Q, et al. Geometrical scaling effect for penetration depth of hard projectiles into concrete targets[J]. International Journal of Impact Engineering, 2018, 120: 46-59.

[16] ZHANG J, CHEN W S, HAO H, et al. Performance of concrete targets mixed with coarse aggregates against rigid projectile impact[J]. International Journal of Impact Engineering, 2020, 141: 103565.

[17] 张义彬, 文玉雷, 朱敏, 等. 混杂纤维增强混凝土抗侵彻性能研究进展[J]. 兵器装备工程学报, 2025, 46(11): 307-318.

[18] ZHANG Y B, WEN Y L, ZHU M, et al. Research progress on penetration resistances of hybrid fiber reinforced concrete[J]. Journal of Ordnance Equipment Engineering, 2025, 46(11): 307-318.

[19] 周冠韬, 徐志洪, 刘聪. 骨料不规则程度对混凝土单轴受压性能影响的数值研究[J]. 应用力学学报, 2025, 42(6): 1292-1300.

[20] ZHOU G T, XU Z H, LIU C. Mesoscale study on the influence of aggregate irregularityon uniaxial compressive properties of concrete[J]. Chinese Journal of Applied Mechanics, 2025, 42(6): 1292-1300.

[21] XU H, WEN H M. A computational constitutive model for concrete subjected to dynamic loadings[J]. International Journal of Impact Engineering, 2016, 91: 116-125.

[22] XU P B, XU H, WEN H M. 3D meso—mechanical modeling of concrete spall tests[J]. International Journal of Impact Engineering, 2016, 97: 46-56.

[23] HERRMANN W. Consitutive equation for the dynamic compaction of ductile porous materials[J]. Journal of Applied Physics., 1969, 40(6): 2490-2499.

[24] 黄晨瑞, 穆朝民, 刘安坤, 等. 高速冲击下混凝土动力学性质和动态温度研究[J]. 爆炸与冲击, 2025, 45(5): 42-59.

[25] HUANG C R, MU C M, LIU A K, et al. Study on dynamic properties and dynamic temperature of concrete under high—speed impact[J]. Explosion and Shock Waves, 2025, 45(5): 42-59.

[26] 马芹永, 简伟强, 张波波. 石墨烯增强模拟月壤地聚合物的动态压缩力学性能[J]. 安徽工业大学学报(自然科学版), 2025, 42(5): 534-545.

[27] MA Q Y, JIAN W Q, ZHANG B B. Dynamic compressive mechanical properties of graphene—reinforced simulated lunar soil geopolymer[J]. Journal of Anhui University of Technology (Natural Science), 2025, 42(5): 534-545.

[28] HARTMANN T, PIETZSCH A, GEBBEKEN N. A hydrocode material model for concrete[J]. International Journal of Protective Structures, 2010, 1(4): 443-468.

[29] MALVAR L J, CRAWFORD J E, WESEVICH J W, et al. A plasticity concrete material model for DYNA3D[J]. International Journal of Impact Engineering, 1997, 19(9/10): 847-873.

[30] GEBBEKEN N, GREULICH S. A new material model for SFRC under high dynamic loadings[C]// Proceedings of the 11th international Symposium on Interaction of the Effects of Munitions with Structures (ISIEMS) Mannheim, Germany. 2003.

[31] 徐浩. 混凝土动态计算本构新模型[D]. 合肥: 中国科学技术大学, 2013.

[32] XU H. A New Computational Constitutive Model for Concrete Subjected to Dynamic Loadings[D]. Hefei: University of Science and Technology of China, 2013.

[33] FORRESTAL M J, FREW D J, HANCHAK S J, et al. Penetration of grout and concrete targets with ogive—nose steel projectiles[J]. International Journal of Impact Engineering, 1996, 18(5): 465-476.

[34] SUN W, ZHANG W, YUAN J, et al. Multi—scale study on penetration performance of steel fiber reinforced ultra—high performance concrete[J]. Construction and Building Materials, 2024, 422: 135846.

[35] 周荣欣, 刘页. 冲击荷载下橡胶超材料混凝土的数值模拟[J]. 高压物理学报, 2025, 39(7): 38-54.

[36] ZHOU R X, LIU Y. Numerical simulation of rubber metamaterial concrete under impact load[J]. Chinese Journal of High Pressure Physics, 2025, 39(7): 38-54.

基金资助

国家自然科学基金项目(12102001)

产学研合作项目(HYB20240205)