基于内聚力法则的高能硝酸酯增塑聚醚推进剂开裂过程细观模型

doi:10.3969/j.issn.1000-1093.2020.11.007

摘要/Abstract

摘要： 为从细观角度研究高能硝酸酯增塑聚醚（NEPE）推进剂的破坏机理，采用分子动力学算法生成细观颗粒填充模型，利用Python脚本语言在颗粒与基体界面及基体内部嵌入零厚度粘结单元。针对NEPE推进剂延展性失效特点，基于多项式内聚力法则建立一种多项式-梯形内聚力法则，并进行子程序VUMAT开发。对比考虑颗粒与基体界面脱湿及基体失效的数值模拟结果发现，NEPE推进剂颗粒与基体界面脱湿引起基体内部形成孔洞，孔洞周围形成的高应力区是导致推进剂开裂的关键。实验验证得知，多项式-梯形内聚力法则较双线性内聚力法则和多项式内聚力法则能更准确地描述推进剂的失效过程。

关键词: 高能硝酸酯增塑聚醚推进剂, 零厚度粘结单元, 多项式-梯形内聚力法则, 颗粒与基体界面脱湿, 基体失效, 数值模拟

Abstract: For the sake of acquiring the mesoscopic damage mechanism of nitrate ester plasticized polyether (NEPE) propellants, a random particle model based on the molecular dynamics theory is established. The zero-thickness cohesive elements are embedded in the matrix and interface in random particle model by using Python scripting language. Considering the ductile failure of NEPE propellant, a novel polynomial-trapezoidal cohesive zone model is presented based on polynomial cohesive zone model, and a vectorized user defined material subroutine VUMAT is developed. By comparing the numerical results of particle/matrix interface dewetting model and matrix damage model, it can be concluded that the particle/matrix interface dewetting of NEPE propellants causes formation of the holes in the matrix, and the high stress around the holes is the major factor that leads to the cracking of propellants. Experimental results validate that the polynomial-trapezoidal cohesive zone model can describe the cracking process of NEPE propellants more accurately compared to bilinear cohesive zone model and polynomial cohesive zone model.

Key words: nitrateesterplasticizedpolyetherpropellant, zero-thicknesscohesiveelement, polynomial-trapezoidalcohesivezonemodel, particle/matrixinterfacedewetting, matrixdamage, numericalsimulation

中图分类号:

V512⁺.3

侯宇菲，许进升，古勇军，周长省. 基于内聚力法则的高能硝酸酯增塑聚醚推进剂开裂过程细观模型[J]. 兵工学报, 2020, 41(11): 2206-2215.

HOU Yufei, XU Jinsheng, GU Yongjun, ZHOU Changsheng. Mesoscopic Model of Cracking Process of NEPE Propellant Based on Cohesive Zone Model[J]. Acta Armamentarii, 2020, 41(11): 2206-2215.

参考文献

[1] KIM S H, IM Y T. Experimental study of material behavior of AP-HTPB base composite solid propellant[J]. Journal of Mechanical Science and Technology, 2019, 33(7): 3355-3361.
[2] 韩龙. 复合固体推进剂细观损伤机理及本构模型研究[D].南京: 南京理工大学, 2017.
HAN L. Research on the mesoscopic damage mechanism and nonlinear viscoelastic constitutive model of composite propellant[D]. Nanjing: Nanjing University of Science and Technology, 2017. (in Chinese)
[3] VAN RAMSHORST M C J, DI BENEDETTO G L, DUVALOIS W, et al. Investigation of the failure mechanism of HTPB/AP/Al propellant by in-situ uniaxial tensile experimentation in SEM[J]. Propellants, Explosives, Pyrotechnics, 2016, 41(4): 700-708.
[4] 刘新国, 刘佩进, 强洪夫, 等. 基于微CT技术的丁羟推进剂脱湿定量表征方法研究[J]. 推进技术, 2019, 40(5): 1162-1168.
LIU X G, LIU P J, QIANG H F, et al. Quantificational method of dewetting damage of HTPB propellant based on micro-CT detection[J]. Journal of Propulsion Technology, 2019, 40(5):1162-1168. (in Chinese)
[5] 高红俐, 刘欢, 齐子诚, 等. 基于高速数字图像相关法的疲劳裂纹尖端位移应变场变化规律研究[J]. 兵工学报, 2015, 36(9): 1772-1781.
GAO H L, LIU H, QI Z C, et al. Research on variation law of fatigue crack-tip displacement and strain fields based on high-speed digital image correlation method[J]. Acta Arrnamentarii, 2015, 36(9): 1772-1781. (in Chinese)
[6] 陈煜, 刘云飞, 夏吉东, 等. NEPE推进剂应力分布的数值模拟及损伤破坏趋势分析[J]. 含能材料, 2009, 17(1): 87-90.
CHEN Y, LIU Y F, XIA J D, et al. Numerical simulation of stress distribution and analysis of damage trend of NEPE propellants[J]. Chinese Journal of Energetic Materials, 2009, 17(1): 87-90. (in Chinese)
[7] MATOU K, INGLIS H M, GU X, et al. Multiscale modeling of solid propellants: from particle packing to failure[J]. Composites Science & Technology, 2007, 67(7/8): 1694-1708.
[8] INGLIS H M. Modeling the effect of debonding on the constitutive response of heterogeneous materials[D]. Urbana-Champaign, IL, US: University of Illinois at Urbana-Champaign, 2014.
[9] CUI H R, LI H Y, SHEN Z B. Cohesive zone model for mode-I fracture with viscoelastic-sensitivity[J]. Engineering Fracture Mechanics, 2019, 221: 106578.
[10] ZHANG L, ZHI S, SHEN Z. Research on tensile mechanical properties and damage mechanism of composite solid propellants[J]. Propellants, Explosives, Pyrotechnics, 2018, 43(3): 234-240.
[11] 周水平, 唐根, 庞爱民，等. GAP复合固体推进剂细观结构演变特性[J]. 固体火箭技术, 2018, 41(5): 580-585,592.
ZHOU S P, TANG G, PANG A M, et al. Evolution of meso-structure of GAP composite solid propellants[J]. Journal of Solid Rocket Technology, 2018, 41(5): 580-585,592. (in Chinese)
[12] DE FRANCQUEVILLE F, GILORMINI P, DIANI J. Represen- tative volume elements for the simulation of isotropic composites highly filled with monosized spheres[J]. International Journal of Solids and Structures, 2019, 158: 277-286.
[13] ZHANG H, AVIJA B, SCHLANGEN E. Towards understanding stochastic fracture performance of cement paste at micro length scale based on numerical simulation[J]. Construction and Building Materials, 2018, 183: 189-201.
[14] VAN RAMSHORST M C J, DI BENEDETTO G L, DUVALOIS W, et al. Investigation of the failure mechanism of HTPB/AP/Al propellant by in-situ uniaxial tensile experimentation in SEM[J]. Propellants, Explosives, Pyrotechnics, 2016, 41(4): 700-708.
[15] MERCIER S, KOWALCZYK- GAJEWSKA K, CZARNOTA C. Effective behavior of composites with combined kinematic and isotropic hardening based on additive tangent Mori-Tanaka scheme[J]. Composites Part B: Engineering, 2019, 174: 107052.
[16] 马昌兵. 复合固体推进剂细观结构建模及其力学行为数值模拟[D]. 西安: 第二炮兵工程学院, 2011.
MA C B. Research on mesostructure modeling and mechanical behaviors numerical simulation of composite solid propellant[D]. Xi'an:The Second Artillery Engineering University, 2011. (in Chinese)
[17] BARENBLATT G I. The mathematical theory of equilibrium cracks in brittle fracture[J]. Advances in Applied Mechanics, 1962, 7: 55-129.
[18] DUGDALE D S. Yielding of steel sheets containing slits[J]. Journal of the Mechanics and Physics of Solids, 1960, 8(2): 100-104.
[19] 张军, 贾宏, 田阳. 粘接界面弹塑性内聚力模型子程序开发[J]. 郑州大学学报 (工学版), 2014, 35 (1): 77-80.
ZHANG J, JIA H, TIAN Y. Development of subroutine for elastic-plastic cohesive zone model of bonded interface[J]. Journal of Zhengzhou University (Engineering Science), 2014, 35(1): 77-80. (in Chinese)
[20] 闫亚宾, 尚福林. PZT 薄膜界面分层破坏的内聚力模拟[J]. 中国科学(G辑:物理学力学天文学), 2009, 39(7): 1007- 1017.
YAN Y B, SHANG F L. Cohesive zone modelling of interface delamination in PZT thin films[J].SCIENTIA SINICA Physica,Mechanica & Astronomica, 2009, 39(7): 1007-1017. (in Chinese)
[21] TVERGAARD V, HUTCHINSON J W. The relation between crack growth resistance and fracture process parameters in elastic-plastic solids[J]. Journal of the Mechanics and Physics of Solids, 1992, 40(6): 1377-1397.
[22] NEEDLEMAN A.Micromechanical modelling of interfacial ecohesion［J］.Ultramicroscopy,1992,40(3):203-214.
[23] 周清春, 鞠玉涛, 周长省. 基于Hooke-Jeeves算法的挠性粘接件的高效内聚反演分析[J].工程力学, 2015, 32(4): 1-7.
ZHOU Q C, J Y T, ZHOU C S. An effective inverse analysis of cohesive parameters of flexible adhesive joints based on the Hook-Jeeves algorithm[J]. Engineering Mechanics, 2015, 32(4): 1-7. (in Chinese)