复合固体推进剂颗粒与基体初始界面有无缺陷的细观模型对比

doi:10.3969/j.issn.1000-1093.2020.09.012

摘要/Abstract

摘要： 复合固体推进剂的制备工艺和环境因素会引起其颗粒与基体初始界面的浸润能力下降，进而导致材料内部存在随机分布的初始缺陷。为探究初始缺陷对推进剂力学性能的影响，基于分子动力学方法建立细观颗粒填充模型，利用通用有限元分析软件Abaqus子程序VUMAT开发指数内聚力模型，并建立颗粒与基体初始界面粘结强度相同的无缺陷模型，以及颗粒与基体初始界面粘结强度的密度函数呈Weibull函数分布的初始缺陷模型。通过对比两种模型的数值计算结果发现：初始缺陷模型的网格敏感性低于无缺陷模型；在初始缺陷模型中，形状参数越小，模型内初始缺陷含量越高，初始弹性模量越低。研究结果表明，初始缺陷模型较无缺陷模型能更准确描述推进剂在单轴拉伸条件下的损伤演化过程。

关键词: 复合固体推进剂, 细观颗粒填充模型, 指数内聚力法则, Weibull分布, 网格敏感性

Abstract: Impregnation capability at particle/matrix interface deteriorates during the preparation of composite solid propellant due to manufacturing and environmental factors，which leads to the initial defects in composite solid propellant. In order to analyze how initial defect affects the mechanical property of solid propellant，a microscopic particle packing model is created by applying the molecular dynamics method. An exponential cohesive model for particle/matrix interface is developed by vectorized user defined material subroutine(VUMAT). Moreover，an initial damage model，in which the distribution of density function of bond strength at propellant particle/matrix interface obeys Weibull function，and a damage-free model of particle/matrix interface with the same bond strength are established.Comparison of the numerical results of two models shows that the mesh sensitivity of the initial damage model is lower than that of the damage-free model.In the initial damage model，the smaller the shape parameter is，the higher the damage is，and the lower the initial modulus is. Through experiments and simulations，it could be seen that the initial damage model can be used to describe the damage characteristic more accurately in uniaxial tension compared to the damage-free model.

Key words: solidcompositepropellant, microscopicparticlepackingmodel, exponentialcohesivelaw, Weibulldistribution, meshsensitivity

中图分类号:

V512⁺.3

侯宇菲，许进升，周长省，陈雄，李宏文. 复合固体推进剂颗粒与基体初始界面有无缺陷的细观模型对比[J]. 兵工学报, 2020, 41(9): 1800-1808.

HOU Yufei， XU Jinsheng， ZHOU Changsheng， CHEN Xiong， LI Hongwen. Comparation of Solid Propellant Micromodels with and without Damage at Initial Particle/matrix Interface[J]. Acta Armamentarii, 2020, 41(9): 1800-1808.

参考文献 23

［1］	王哲君，强洪夫，王广，等.低温高应变率条件下HTPB推进剂拉伸力学性能研究[J].推进技术，2015，36(9):1426-1432.WANG Z J，QIANG H F，WANG G，et al. Tensile mechanical properties of HTPB propellant at low temperature and high strain rate[J].Journal of Propulsion Technology，2015，36(9):1426-1432.(in Chinese)
[2]	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.
[3]	王阳，李高春，张璇，等.基于SEM与数字图像相关方法的HTPB推进剂裂尖扩展过程分析[J].火炸药学报，2019，42(1):73-78.WANG Y，LI G C，ZHANG X，et al. Analysis of crack tip propagation process of HTPB propellant based on SEM and digital image correlation method[J].Chinese Journal of Explosives & Propellants，2019，42(1):73-78. (in Chinese)
[4]	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.
[5]	常武军，鞠玉涛，王蓬勃.HTPB推进剂脱湿与力学性能的相关性研究[J].兵工学报，2012，33(3):261-266.CHANG W J，JU Y T，WANG P B. Research on correlation between dewetting and mechanical property of HTPB propellant[J].Acta Armamentarii，2012，33(3):261-266.(in Chinese)
[6]	INGLIS H M. Modeling the effect of debonding on the constitutive response of heterogeneous materials[D].Champaign，IL，US: University of Illinois at Urbana-Champaign，2014.
[7]	刘承武，阳建红，邓凯，等.基于Mori-Tanaka有限元法的粘弹复合推进剂非线性界面脱粘[J]. 推进技术，2011，32(2): 225-229. LIU C W，YANG J H，DENG K，et al.Nonlinear interface debonding of viscoelastic composite propellant using Mori-Tanaka finite element method[J].Journal of Propulsion Technology，2011，32(2): 225-229.(in Chinese)
[8]	韩龙.复合固体推进剂细观损伤机理及本构模型研究[D].南京:南京理工大学，2017.HAN L.Research on the mesoscopic damage mechanism and nonlinear viscoelastic constitutive model of composite propellant[D].Nanjing: Nanjing University of Scienceand Technology，2017. (in Chinese)
[9]	DAI K，LU B，CHEN P，et al. Modelling microstructural deformation and the failure process of plastic bonded explosives using the cohesive zone model[J].Materials，2019，12(22): 3661-3682.
[10]	BALZER J E，SIVIOUR C R，WALLEY S M，et al. Behaviour of ammonium perchlorate-based propellants and a polymer-bonded explosive under impact loading[J].Proceedings of the Royal Society of London. Series A:Mathematical，Physical and Engineering Sciences，2004，460(2043):781-806.
[11]	刘晋湘，梁蓓，朱立勋，等.Ⅰ类AP初始缺陷对丁羟推进剂力学性能的影响[J].火炸药学报，2018，41(3):273-277.LIU J X，LIANG B，ZHU L X，et al. Effect of initial defect of type Ⅰ AP on the mechanical properties of HTPB propellant[J].Chinese Journal of Explosives & Propellants，2018，41(3):273-277.(in Chinese)
[12]	黄海峰，巨能攀，蓝康文，等.岩石统计损伤软化模型及其参数反演[J].长江科学院院报，2018，35(6):102-106.HUANG H F，JU N P，LAN K W，et al. Statistical damage softening model for rock and back analysis of its parameters[J]. Journal of Yangtze River Scientific Research Institute，2018，35(6):102-106.(in Chinese)
[13]	文志杰，田雷，蒋宇静，等.基于应变能密度的非均质岩石损伤本构模型研究[J].岩石力学与工程学报，2019，38(7):1332-1343.WEN Z J，TIAN L，JIANG Y J, et al. Research on damage constitutive model of inhomogeneous rocks based on strain energy density[J].Chinese Journal of Rock Mechanics and Engineering，2019，38(7): 1332-1343.(in Chinese)
[14]	BERNARD F，KAMALI-BERNARD S, PRINCE W.3D multi-scale modelling of mechanical behaviour of sound and leached mortar[J].Cement and Concrete Research，2008，38(4):449-458.
[15]	宋丹平.固体推进剂细观力学与本构关系研究[D].武汉：武汉理工大学，2008.SONG D P.Research on the micro-mechanics and constitutive relation for solid propellant[D].Wuhan:Wuhan University of Technology, 2008. (in Chinese)
[16]	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.
[17]	DUGDALE D S.Yielding of steel sheets containing slits[J].Journal of the Mechanics and Physics of Solids，1960，8(2):100-104.
[18]	BARENBLATT G I.The mathematical theory of equilibrium cracks in brittle fracture[J]. Advances in Applied Mechanics，1962，7(1): 55-129.
[19]	丁云龙. 基于多尺度内聚力模型的PRMMCs动态力学行为研究[D]. 武汉：武汉理工大学，2015.DING Y L. The study of dynamic behavior of particle reinforced metal matrix composites based on multi-scale cohesive zone model[D].Wuhan: Wuhan University of Technology, 2015. (in Chinese)
[20]	徐强，胡勤伟，沙宝林，等.NEPE固体推进剂含损伤的黏弹性本构研究[J].应用力学学报，2019，36(3): 652-657.XU Q，HU Q W，SHA B L，et al. Research on viscoelasticity model including damage for NEPE solid propellant[J].Chinese Journal of Applied Mechanics，2019，36(3): 652-657. (in Chinese)
[21]	ARORA H，CHARAKAMBIDS M N，TARLETON E，et al. An image based approach to modelling plastic bonded explosives (PBX) on the micro scale[C]∥Proceedings of International Conference on Composite Materials. Montreal, Canada:Canadian Association for Composite Structures and Materials，2013.
[22]	周清春，鞠玉涛，周长省.基于Hooke-Jeeves算法的挠性粘接件的高效内聚反演分析[J].工程力学，2015，32(4): 1-7.ZHOU Q C，JU Y T，ZHOU C S.An effective inverse analysis of cohesive parameters of flexible adhesive joints based on the Hooke-Jeeves algorithm[J].Engineering Mechanics，2015，32(4): 1-7. (in Chinese)
[23]	韩波，鞠玉涛，周长省. HTPB推进剂粘聚区本构模型反演识别研究[J].兵工学报，2012，33(11):1335-1341.HAN B，JU Y T，ZHOU C S.Inversion identification of constitutive model of cohesive zone in HTPB propellant[J].Acta Armamentarii， 2012，33(11):1335-1341.(in Chinese)

[1]	叶振威，余永刚. 脉冲点火射流与高氯酸铵/端羟基聚丁二烯固体推进剂耦合燃烧的试验研究及数值模拟[J]. 兵工学报, 2018, 39(8): 1507-1514.
[2]	邱明，周大威，周占生. 基于加速寿命试验的自润滑关节轴承可靠性分析[J]. 兵工学报, 2018, 39(7): 1429-1435.
[3]	王浩伟，滕克难，奚文骏. 非恒定环境下基于载荷谱的导弹部件寿命预测[J]. 兵工学报, 2016, 37(8): 1524-1529.
[4]	张详坡，尚建忠，陈循，张春华，汪亚顺. 三参数Weibull分布竞争失效场合加速寿命试验统计分析[J]. 兵工学报, 2013, 34(12): 1603-1610.
[5]	1：安伟光:2：周源泉. Weibull分布变差系数上限及蓼对数插值方法[J]. 兵工学报, 1997, 18(2): 129-133.