一种基于材料细观特征量的PBX拉伸强度理论模型_李克武.pdf

DOI：10.11883/bzycj-2021-0514一种基于材料细观特征量的PBX 拉伸强度理论模型*李克武1，胡秋实1，郑贤旭1，李涛1，傅华1，唐维2（1.中国工程物理研究院流体物理研究所，四川绵阳621999；2.中国工程物理研究院化工材料研究所，四川绵阳621999）摘要：为给塑性黏结炸药（PBX）的力学强度设计提供支撑、探索材料细观特征量与材料强度之间的定量规律，应用微裂纹扩展区理论，将 PBX 炸药的单轴拉伸过程中力学响应特征的变化归结为扩展裂纹取向角度的增加，将扩展裂纹最大取向角与拉伸强度相关联，构建了基于材料细观特征量的拉伸强度理论模型，并采用不同温度的单轴拉伸实验验证了该理论模型的有效性。研究表明：该拉伸强度理论模型可以实现对 PBX 炸药拉伸强度与炸药微裂纹密度、颗粒/黏结剂界面性能以及颗粒/黏结剂体系的表观杨氏模量、泊松比等细观特征量之间关系的定量描述。关键词：塑性黏结炸药；拉伸强度；理论模型；细观损伤；微裂纹扩展区理论中图分类号：O383国标学科代码：13035文献标志码：AA theoretical model of PBXs tensile strengthbased on meso-structure parametersLIKewu1,HUQiushi1,ZHENGXianxu1,LITao1,FUHua1,TANGWei2(1.Insititute of Fluid Physics,China Academy of Engineering Physics,Mianyang 621999,Sichuan,China;2.Institute of Chemical Materials,China Academy of Engineering Physics,Mianyang 621999,Sichuan,China)Abstract:In practical applications,plastic bonded explosive(PBX)explosives are often used as load bearing structuralcomponents.Hence,themechanicalstrengthisanimportantparameterinPBXsdesign.ItisofgreatinteresthowtoadjustthemicroscopiccharacteristicsofthematerialinthemanufacturingprocesstoobtainanPBXwiththerequiredstrength.PBXconsists of a large portion of energetic particles with a small portion of binder.Therefore,a large number of randomlydistributedmicrocracksexistinside.InspiredbythesuccessfulapplicationofmicrocrackextensionmechanisminthePBXsmechanical response simulation in recent years,the theory related to microcrack extension was applied to the strengthmodelingstudy.Thedomainofmicrocrackgrowth(DMG)theorywasappliedtoanalyzePBXsuniaxialtension.Theresultsshowsthattheincreaseintensileloadresultedinmicrocrackextensionatdifferentorientationangles,aswellasadecreaseinthecrackinterval.Meanwhile,recentstudiesonthecascadingbehaviorofrandomlydistributedmicrocracksshowthatthecrackcascadingbehaviorwitheachotheriscontrolledbytheorientationangleandcrackintervalofadjacentcracks.ThePBXtensilefracturecomesfromtheinstableextensionofmacroscopiccracks,andthemacroscopiccrackformationandextensionbothcomefromthecascadingofmicrocracks.Bytreatingthestrengthastheminimumloadrequiredforcatastrophiccrackcascadeextension,thereisaone-to-onecorrespondencebetweenthestrengthandthemaximumorientationangleofextendedcrack.BasedontheDMGtheory,thisone-to-onecorrespondenceisexpressedasthetheoreticalmodeloftensilestrength.Thistheoreticalmodelestablishestheconnectionbetweenthetensilestrengthandthemeso-structureparameterssuchasstress*收稿日期：2021-12-15；修回日期：2022-09-25 基金项目：国家自然科学基金（11802288，11802283，11902306）第一作者：李克武（1983），男，博士，助理研究员，第43卷第1期爆炸与冲击Vol.43,No.12023年1月EXPLOSIONANDSHOCKWAVESJan.,2023013106-1intensityfactor,microcrackdiameterandmicrocrackdensity.Comparisonoftensilestrengthmodelpredictionswithwidetemperaturerange(from40to45)experimentaldata,indicatedthatthismodeliscapabletodescribethePBXstensilestrength.SuggestionsforthedesignofstrengthenhancementofPBXscanbeprovidedbythistheoreticalmodel.Keywords:PBX;tensilestrength;theoreticalmodel;mesodamage;DMGtheory塑性黏结炸药（plasticbondedexplosive,PBX）在实际应用中除作为毁伤能量的主要提供部件外，通常还作为结构件承担载荷，因此炸药的力学强度是一个重要指标。一方面，强度是材料的一种力学特征量，通过力学实验确定 PBX 在不同应变率与不同温度下的强度，可以确定其安全使用范围。Stevens1基于大量的力学实验数据，分析了环境温度和加载应变率对PBX-9501 强度的影响，利用Williams-Landel-Ferry 时温转换函数构建了能够定量描述温度和应变率影响的强度模型，确定了 PBX-9501的安全载荷边界。另一方面，强度也是多个材料参量共同作用下的综合指标，通过各种方法调控相关参数，可以实现材料强度的提高。例如，肖磊等2开展了 PBX 炸药颗粒级配研究，发现通过添加纳米级炸药晶体颗粒可以有效提高材料拉伸强度，增幅最高可达 16.7%；黄辉3针对浇注固化炸药，将颗粒级配与悬浮体系的流变性能相关联，建立了体现颗粒级配影响的流变模型，提出了一种降低表观黏度提升力学性能的最佳颗粒级配计算方法；Lv 等4则对 PBX 制备过程中的超声波辅助工艺进行了研究，认为在此工艺过程中较大粒径颗粒破碎引起粒径级配变化，可以提高材料强度。此外，通过强化颗粒/黏结剂界面也可以提升材料强度，如黄辉等5分析了多种偶联剂对浇注固化炸药强度提高的作用，通过实验确定了酰胺类化合物为最佳偶联剂；Li 等6研究了 4 种含硼偶联剂对 PBX 炸药界面的强化作用，实验结果证实偶联剂均可有效提高界面黏附功，从而提高材料强度；除偶联剂外，Yang 等7还建立了在 PBX 炸药晶体外添加分子涂层以提高颗粒/黏结剂界面性能的方法，力学实验证实该种方法可以将拉伸强度提高 20%以上。然而，开展各种工艺探索以有效提高材料强度的同时，还迫切需要深入了解各材料参量对材料强度的影响机制，确定各材料参量对材料强度的定量影响规律，以方便 PBX 炸药性能的设计规划。本文中，在细观损伤理论的框架内，应用微裂纹扩展区（domainofmicrocrackgrowth，DMG）理论对PBX 炸药的单轴拉伸过程进行分析，将 PBX 炸药拉伸响应特征通过扩展裂纹取向角度进行描述；同时引入随机分布裂纹相互连接行为的研究成果，将扩展裂纹取向角的最大值与裂纹灾难性延伸相关联；考虑到拉伸强度是裂纹灾难性延伸的重要因素，将拉伸强度与扩展裂纹最大取向角对应，构建一种基于细观结构参数的拉伸强度理论模型，以实现细观结构参量对 PBX 炸药拉伸强度影响的定量描述。1 微裂纹扩展区理论简介微裂纹扩展区理论是冯西桥等8在细观损伤理论框架下针对准脆性材料准静态力学响应特征而发展起来的一套理论。在细观损伤理论框架下，由于非均匀材料的宏观力学性能仍然是统计均匀的，因此常引入代表性体积单元（representativevolumeelement,RVE）的概念，用以分析非均匀材料的有效本构关系。在细观尺度上，该 RVE 尺寸足够大，包含有足够多的细观结构，从而能够反映材料的统计平均性质；同时，在宏观尺度上，RVE 尺寸又足够小，作用在 RVE 上的宏观应力与应变可认为是均匀的。Ox1x2x3微裂纹扩展区理论假定了在 RVE 内分布大量币形微裂纹，所有微裂纹都具有相同的统计平均半径a0，各条微裂纹之间无相互影响，并且微裂纹取向随机。对任意一条微裂纹均可建立整体坐标系Ox1x2x3与局部坐标系，如图 1 所示。基于坐标系变换准则，微裂纹的取向使用角度参数 和实现唯一定义，其中 的取值范围为 0,/2)，的取值范围为 0,2)。在拉伸载荷 22作用下，当裂纹的应力强度因子满足下式时，裂纹发生自相似扩展，并在遇到能障后停止8。扩展后的微裂纹半径为 au。第43卷李克武，等：一种基于材料细观特征量的PBX拉伸强度理论模型第1期013106-2KKC2+KKC2=1K=222a0K=42a0(221+223)(1)KKij式中：上标“”表示该物理量基于局部坐标系，和分别为型和型临界应力强度因子，KC和 KC为裂纹临界应力强度因子，为局部坐标系下应力张量的各个分量。对于单轴拉伸状态，由于裂纹取向角相对于拉伸载荷的对称性，该参量不影响裂纹扩展行为，式(1)可以化简为：KCKC2222cos4sup+22v2222cos2supsin2sup4a0K2C=0(2)式中：v 为材料的泊松比，sup为扩展裂纹取向角的上限。KC/(4a0)由式(2)可知，当拉伸载荷增至时，法线平行于载荷方向的裂纹（=0）首先发生扩展。并且，随着拉伸载荷的持续增加，更大角度范围的裂纹发生扩展，sup持续增大。当载荷趋向无穷大时，sup随载荷趋向/2，所有取向的裂纹均发生扩展。由于式(2)中的 KC和 KC均为材料常数，因此拉伸载荷和 sup一一对应。随着载荷的增加，越来越多的能量耗散在裂纹扩展行为中，宏观表现为材料刚度下降损伤增加。该理论通过扩展裂纹取向角上限 sup来描述材料损伤状态的变化。2 PBX 的拉伸强度理论模型PBX 由大量含能颗粒与少量黏结剂组成，存在着大量的颗粒/黏结剂界面。受载过程中的大部分损伤破坏为沿这些界面的裂纹扩展9，在单轴拉伸载荷下表现为典型的准脆性特征10-11。针对 PBX 炸药的准脆性响应，许多研究者们基于 RVE 概念，将真实 PBX 细观结构缺陷简化为均质体内的微裂纹体系，建立了诸如成核与增长（nucleationandgrowth,NAG）模型12、统计微裂纹（statisticalcrackmechanics,SCRAM）模型13以及主裂纹（dominantcrackalgorithm,DCA）模型14，实现了 PBX 力学响应特征的准确描述。依据力学实验与模型计算的对比分析结果，细观损伤理论框架下的微裂纹扩展理论在描述 PBX 力学响应方面表现出足够的准确度15。因此，本文中也将 PBX 视为微裂纹体，应用前述的微裂纹扩展区理论分析 PBX 的拉伸断裂行为。真实的拉伸断裂来自于宏观失稳裂纹的扩展传播16，而 PBX 宏观主裂纹形成与扩展均来自于微裂纹的连接17，因此假定 PBX 的拉伸断裂过程如下：在 PBX 发生拉伸断裂前，最大扩展裂纹取向角 sup随拉伸载荷持续增大，PBX 内部大量裂纹发生扩展，与相邻裂纹的裂尖间距减小，而当载荷增至某一阈值时，裂尖间距减小的相邻微裂纹间发生相互连接，成为失稳大裂纹并在载荷作用下持续与相邻微裂纹串接，形成灾难性裂纹延伸。对于代表性体积单元，其内部各种微裂纹状态均存在。换言之，最易于发生微裂纹相互连接形成失稳扩展的微裂纹分布状态也是存在的。因此，拉伸强度即为微裂纹发生失稳连接的最小载荷阈值。相邻裂纹间相互连接行为如图 2 所示。Zhang 等18研究发现，能量释放率达到阈值时，裂纹即与相邻裂纹连接，具体表示为：G(a/d,)K2E=GC(3)x2x2aOx1x3x1x3图1微裂纹的整体坐标系与局部坐标系8Fig.1Globecoordinateandlocalcoordinateofmicro-crack8第43卷李克武，等：一种基于材料细观特征量的PBX拉伸强度理论模型第1期013106-3K式中：a 为裂纹半径；d 为裂尖至相邻裂纹裂心的距离；为裂尖至相邻裂纹裂心连线与裂纹方向的夹角；为扩展裂纹取向角；为远场应力强度因子；对平面应变状态，E=E/(1v2)，对平面应力状态 E=E，其中 E 为材料的杨氏模量；GC为临界能量释放率。G(a/d,)Gmax(a/d,)max对于式(3)内的各个变量，据微裂纹扩展区理论，扩展后的微裂纹半径为 au，因此 a=au。au和 d 为材料常数，au与材料的颗粒级配相关8，d 为单位体积裂纹数 n 的函数。由于代表性体积单元内部的各种微裂纹状态均存在，因此总是可以取到令为最大值的。此外，杨氏模量与临界能量释放率均为材料常数。因此是否会形成微裂纹相互连接，受拉伸载荷与扩展裂纹取向角控制。f据微裂纹扩展区理论，拉伸载荷与扩展裂纹取向角存在一一对应关系。在材料拉伸破坏前，扩展裂纹取向角上限 sup随拉伸载荷增加而持续增加。sup增至 max使得式(3)成立时，材料拉伸破坏，拉伸强度即为增至 max时的拉伸载荷。因此，由式(2)有：f=KC2a0/cos4max+4cos2maxsin2maxK2C(2v)2K2C(4)式(4)即为基于材料细观特征量的 PBX 拉伸强度理论模型，其中，KC和 KC反映 PBX 炸药晶体颗粒/黏结剂界面性能，因此本模型描述强化界面对材料强度提高的规律。而 n、a0和 au与炸药晶体颗粒级配相关8，因此本模型还可以用于分析颗粒级配对材料强度的影响。3 算例应用拉伸强度理论模型，对4545 的 PBX-3 炸药准静态拉伸试验进行分析。该 PBX-3以奥克托金（HMX）为基，HMX、TATB（三氨基三硝基苯）、黏结剂和钝感剂的质量分数分别为 87.0%、7.0%、4.2%和 1.8%，在模具中热压成型。其拉伸响应表现出明显的准脆性特征，应力-应变曲线如图 3 所示。本文中使用该应力和应变数据确定 PBX-3 在不同温度下的细观特征量，并将实验测得的拉伸强度与模型预测强度进行对比。据微裂纹扩展区理论8，材料单轴拉伸响应过程可以分为弹性响应阶段与非弹性响应阶段。2 个阶段的应力-应变关系为=E1+1645AKC2a0E1+1645A+1645B(103v10cos3max+3vcos5max)KC2a0(5)式中：为应力，为应变，A=na03(1v2)(103v)/(2v)，B=n(au3a03)(1v2)(103v)/(2v)。应用式(5)对该实验获得的 7 条不同温度应力-应变曲线进行全局拟合，其中 n、a0和 au不受温度影响，拟合时共用同一套参数。而 KC和 KC随温度变化，拟合时每个温度设置一套参数。杨氏模量 E 由xxyy2audKxConnectingxyymax2auK图2相邻裂纹的连接Fig.2Connectingofneighboringcracks第43卷李克武，等：一种基于材料细观特征量的PBX拉伸强度理论模型第1期013106-4应力-应变曲线直接测量获得。拟合结果如图 3 所示。结果表明，当环境温度在4045 范围时，微裂纹扩展区模型可以较好地描述 PBX 的拉伸响应行为，拟合得到的细观特征量较为准确。所确定的 PBX-3 各细观特征量如表 12 所示。由于扩展裂纹最大取向角 max与温度不直接相关，因此将根据25 时拉伸强度所确定的 max=1.317rad 代入理论模型，预测其余温度点的拉伸强度。与实验测得的拉伸强度对比如图 4 所示，可以看出最大相对误差不超过 10%，表明该拉伸强度理论模型可以用于 PBX 炸药拉伸强度分析。0.020.040.060.08012345678Stress/MPa012345678Stress/MPa012345678Stress/MPa012345678Stress/MPa012345678Stress/MPa012345678Stress/MPa Experimental value Fitting data Experimental value Fitting data Experimental value Fitting data Experimental value Fitting data Experimental value Fitting data Experimental value Fitting data(a)40 0.020.040.060.020.040.060.020.040.060.080.020.040.060.080.020.040.06Strain/%(b)30 Strain/%(c)20 Strain/%(d)10 Strain/%(e)25 Strain/%(f)45 Strain/%图3拟合结果Fig.3Fittingresults表 1 与温度无关的 PBX-3 细观特征量Table 1 Temperature independent mesoscopic characteristics of PBX-3n/cm3a0/cmau/cmv13720.0063520.037030.3表 2 与温度相关的 PBX-3 细观特征量Table 2 Temperature dependent mesoscopic characteristics of PBX-3T/KC/（MPacm1/2）KC/（MPacm1/2）400.16280.1911300.15020.1763250.15560.1827200.15200.1784100.14570.1711250.14570.1711450.13220.1552第43卷李克武，等：一种基于材料细观特征量的PBX拉伸强度理论模型第1期013106-5由拉伸强度理论模型可知，PBX 炸药的拉伸强度受颗粒级配与晶体颗粒/黏结剂界面性能控制。由于颗粒级配不受温度影响，因此 max基本恒定不变。但是随着温度的上升，黏结剂受热软化，颗粒更易于错动，模型中表现为临界应力强度因子随温度持续下降。据强度理论模型，随着温度的上升，模型预测的拉伸强度持续下降，如图 4 所示，与实验测量结果变化趋势一致。4 讨论本文中，以 PBX 炸药为研究目标，围绕其拉伸强度及细观特征量影响机理展开研究，提出了一个基于材料细观特征量的拉伸强度理论模型。相较于已有 PBX 拉伸强度的研究成果2-7，本文中所获得的结果着重从理论上解释各个细观特征量对宏观力学强度的影响规律，并给出了定量的理论模型。基于本文中所获得的理论模型，可以为 PBX 炸药性能的设计规划提供支持。例如，强化颗粒/黏结剂界面以提高材料力学强度的工艺，是通过调整材料细观特征量中的型/型临界应力强度因子以实现 PBX 强度提高。由于一部分黏结剂在提高型临界应力强度因子更有效，一部分黏结剂提高型临界应力强度因子更有效，在面对黏结剂种类选择时，可应用强度理论模型辅助确定。以 PBX-3 为例，基于所获得的细观特征量 n、a0、au以及泊松比 v，应用拉伸强度理论模型（式(4)）可以给出强度变化曲线，如图 5 所示。据图 5(a)，单纯提高型临界应力强度因子存在边际效应，当 KC0.02MPacm1/2时，提升材料强度的效果明显减弱。据图 5(b)，KC对材料强度的提升几乎没有边际效应，因此选择能够能有效提高 KC的黏结剂是一个较为合理的选择。50 40 30 20 10010203040500123456789Tension strength/MPaT/Experimental value Model prediction图4理论模型预测与实验结果对比Fig.4Comparisonofmodelpredictionandexperimentalvalue0.040.080.120.160.200123456789Tension strength/MPaTension strength/MPaKC/(MPacm1/2)KC/(MPacm1/2)KC=0.14 MPacm1/2KC=0.16 MPacm1/2KC=0.18 MPacm1/2KC=0.20 MPacm1/2KC=0.22 MPacm1/2KC=0.12 MPacm1/2KC=0.14 MPacm1/2KC=0.16 MPacm1/2KC=0.18 MPacm1/2KC=0.20 MPacm1/20.040.080.120.160.20012345678图5两种应力强度因子（KC和 KC）对材料强度的影响Fig.5Curvesoftensilestrengthwithtwostressintensityfactors(KCandKC)第43卷李克武，等：一种基于材料细观特征量的PBX拉伸强度理论模型第1期013106-6此外，依据强度理论模型，还能够直接推算强度高于指定值所需要的细观特征量范围。仍以 PBX-3为例，若计划通过调整颗粒/黏结剂界面性能的技术途径来增强该 PBX 的拉伸强度至 10MPa。那么将拉伸强度 10MPa 代入拉伸强度理论模型式(4)中，基于所获得的细观特征量 n、a0、au以及泊松比 v，可以给出材料强度不低于 10MPa 时，KC和 KC需满足：1K2C+20.6K2C1(0.0566 MPacm1/2)2(6)进而，可以依据式(6)所确定的材料参数范围，选择合理的具体技术手段以实现力学强度设计目标。本研究也存在一定欠缺，主要表现在该强度理论模型只适用于脆性或准脆性材料。这是因为本模型基于微裂纹扩展区理论，其非弹性响应被认为源于微裂纹扩展导致的能量耗散。当 PBX 炸药非弹性响应主要因为塑性滑移耗散时，本理论模型不再适用。5 结论基于微裂纹扩展区理论与裂纹连接准则，分析了 PBX 炸药单轴拉伸断裂行为，提出了一种基于细观特征量的拉伸强度理论模型，展示了各细观特征量对材料强度的影响机制。并通过 PBX-3 的拉伸实验数据对模型进行了验证，证实该模型可应用于 PBX 炸药拉伸强度研究。主要结论有：f=KC2a0 cos4max+4cos2maxsin2maxK2C(2v)2K2C(1)基于细观特征量的PBX炸药拉伸强度理论模型为；(2)PBX 炸药的拉伸强度与材料的炸药颗粒粒径分布、颗粒/黏结剂界面性能以及颗粒/黏结剂体系的表观杨氏模量、泊松比相关；(3)基于实验对比结果，本拉伸强度理论模型可以较好地描述 PBX-3 材料的强度变化特征。参考文献：STEVENS R.A strength model and service envelope for PBX 9501:LA-UR-14-20696 R.Los Alamos:Los AlamosNationalLaboratory,2014.1肖磊,刘杰,郝嘎子,等.微纳米 RDX 颗粒级配对压装 PBX 性能影响J.含能材料,2016,24(12):11931197.DOI:10.11943/j.issn.1006-9941.2016.12.011.XIAOL,LIUJ,HAOGZ,etal.Effectsofnano-/micrometerRDXparticlegradationonthepropertyofPBXJ.ChineseJournalofEnergeticMaterials,2016,24(12):11931197.DOI:10.11943/j.issn.1006-9941.2016.12.011.2黄辉.颗粒级配技术及其在含能材料中的应用J.含能材料,2001,9(4):161164.DOI:10.3969/j.issn.1006-9941.2001.04.005.HUANGH.ParticlegradetechniqueandapplicationonenergeticmaterialsJ.EnergeticMaterials,2001,9(4):161164.DOI:10.3969/j.issn.1006-9941.2001.04.005.3LVKZ,YANGK,ZHOUB,etal.Thedensificationandmechanicalbehaviorsoflarge-diameterpolymer-bondedexplosivesprocessedbyultrasonic-assistedpowdercompactionJ.Materials&Design,2021,207:109872.DOI:10.1016/j.matdes.2021.109872.4黄辉,王晓川.偶联剂在 HMX 基浇注固化炸药中的作用J.含能材料,2000,8(1):1317.DOI:10.3969/j.issn.1006-9941.2000.01.004.HUANGH,WANGXC.BehaviorofcouplingagentinHMX-basedextrudablecastexplosiveJ.EnergeticMaterials,2000,8(1):1317.DOI:10.3969/j.issn.1006-9941.2000.01.004.5LIF,YEL,NIEFD,etal.Synthesisofboron-containingcouplingagentsanditseffectontheinterfacialbondingoffluoropolymer/TATBcompositeJ.JournalofAppliedPolymerScience,2007,105(2):777782.DOI:10.1002/app.26117.6YANGXL,GONGFY,ZHANGK,etal.EnhancedcreepresistanceandmechanicalpropertiesforCL-20andFOX-7basedPBXsbycrystalsurfacemodificationJ.Propellants,Explosives,Pyrotechnics,2021,46(4):572578.DOI:10.1002/prep.202000277.7第43卷李克武，等：一种基于材料细观特征量的PBX拉伸强度理论模型第1期013106-7冯西桥,余寿文.准脆性材料细观损伤力学M.北京:高等教育出版社,2002:2041.FENGXQ,YUSW.Damagemicromechanicsofquasi-brittlematerialsM.Beijing:HigherEducationPress,2002:2041.8YANGYZ,LIJL,CHENJK.ExperimentalandnumericalanalysisofshearprocessofahighparticlecontentbondingmaterialJ.JournalofAppliedMechanics,2021,88(6):061003.DOI:10.1115/1.4050277.9LIUC,THOMPSONDG.CrackinitiationandgrowthinPBX9502highexplosivesubjecttocompressionJ.JournalofAppliedMechanics,2014,81(10):101004.DOI:10.1115/1.4028087.10LIUM,HUANGXC,WUYQ,etal.Numericalsimulationsofthedamageevolutionforplastic-bondedexplosivessubjectedtocomplexstressstatesJ.MechanicsofMaterials,2019,139:103179.DOI:10.1016/j.mechmat.2019.103179.11CURRANDR,SEAMANL,SHOCKEYDA.DynamicfailureofsolidsJ.PhysicsReports,1987,147(5/6):253388.DOI:10.1016/0370-1573(87)90049-4.12BENNETTJG,HABERMANKS,JOHNSONJN,etal.Aconstitutivemodelforthenon-shockignitionandmechanicalresponseofhighexplosivesJ.JournaloftheMechanicsandPhysicsofSolids,1998,46(12):23032322.DOI:10.1016/s0022-5096(98)00011-8.13ZUOQH,ADDESSIOFL,DIENESJK,etal.Arate-dependentdamagemodelforbrittlematerialsbasedonthedominantcrackJ.InternationalJournalofSolidsandStructures,2006,43(11/12):33503380.DOI:10.1016/j.ijsolstr.2005.06.083.14RAGASWAMYP,LEWISMW,LIUC,etal.ModelingthemechanicalresponseofPBX9501:LA-UR-10-01416R.LosAlamos:LosAlamosNationalLaboratory,2010.15PIJAUDIER-CABOTG,BAANTZP.NonlocaldamagetheoryJ.JournalofEngineeringMechanics,1987,113(10):15121533.DOI:10.1061/(asce)0733-9399(1987)113:10(1512.16LIU C.Damage and fracture in high explosive materials:LA-UR-11-05618 R.Los Alamos:Los Alamos NationalLaboratory,2011.17ZHANGSL,YANGW.MacrocrackextensionbyconnectingstatisticallydistributedmicrocracksJ.InternationalJournalofFracture,1998,90(4):341353.DOI:10.1023/A:1007440012284.18(责任编辑王小飞)第43卷李克武，等：一种基于材料细观特征量的PBX拉伸强度理论模型第1期013106-8