一种含能材料细观热点物理模型的建模方法
【技术领域】
[0001] 本发明属于含能材料热安全技术领域，涉及一种含能材料细观热物理模型的建模 方法。
【背景技术】
[0002] 对于含能材料热分解，通常伴随有热量的释放。如果热分解所释放的热量不能扩 散到周围环境，将有可能引发自燃甚至爆炸事故。
[0003] 目前对于热安全物理模型主要是基于宏观均质假设的化学反应-传热模型，这一 模型在仿真宏观烤燃过程中得到较好的应用，但采用宏观均质化的仿真模型不能揭示含能 材料热安全的深层次的机制，对于含能材料热安全配方设计和危险性起不到预见作用。
[0004] 含能材料的细观尺度决定其热安全和热爆炸机理，在细观尺度上含能材料是非均 质性的，可以看作是由氧化剂颗粒、粘合剂基体以及两者之间的粘结界面带所组成的三相 非均质复合材料。各组分不一致的传热性能和分解放热性能、组份间界面以及各组份在热 分解过程中所遵循的不一致的化学反应机理，以及在热分解过程中各组份内部及组份间界 面所发生的复杂的结构演化，这些都是决定了含能材料热安全的复杂性。
[0005] 为了揭示含能材料热安全深层次的机制，本发明综合多种宏细观试验方法，从含 能材料及其组分的热分解动力学参数和机理函数、含能材料及其组分的微结构表面形貌演 化规律、含能材料及其组分的微结构内部形貌演化规律、对含能材料组分热物理常数的测 定四个方面展开，以建立正确的含能材料细观热点物理模型。
[0006] 在本发明以前的现有技术中，所检索到的文献有：[1]刘礼斌等温度对发射药安全 和能量影响的数值计算[D].弹道学报.2004:2; [2]陶文锉.数值传热学[M].西安交通大学 出版社，2001;[ 3]杜志明.有限空间内化学放热系统的热点火[D].北京：北京理工大学， 1993; [4]荆松吉.凝聚炸药烤燃机理研究及二维数值模拟[D].长沙：国防科技大学，2004; 上述文献是与本发明技术主题相关的同类别技术，但是，从所公开的内容看，尚未发现有关 含能材料细观热物理模型构建方法的报道。
【发明内容】
[0007] 根据上述【背景技术】，本发明的目的在于，为研究含能材料的热安全机理和热安全 临界温度及爆炸延迟时间的预测提供依据，为改善含能材料热稳定性配方提供一种能够构 建高质量高效率的含能材料细观热点物理模型的建模方法。
[0008] 现将本发明构思及技术解决方案叙述如下：
[0009] 本发明的基本构思是，在考虑细观结构形成及演变和温度对组分材料的热物理性 能的影响，考虑的影响因素全面合理，且形式简单、使用方便的前提下，对其细观热点物理 模型的建立方法进行描述，该方法包括:热分解机理函数和参数的确定;对含能材料微结构 表面形貌演化规律的研究;对含能材料微结构内部形貌演化规律的研究;对含能材料组分 热物理常数的测定四个部分，具体包括以下步骤：
[0010] 步骤1:通过DSC法和激光导热法确定组分材料的热容和导热率；
[0011] 步骤2:基于DSC-TG试验确定热分解机理函数和参数；
[0012] 步骤3:通过扫描电镜试验获取含能材料表面形貌随温度的演化规律，进行定量分 析，得出表面形貌随温度演化规律；
[0013] 步骤4:通过yCT试验获取含能材料内部形貌随温度的演化规律，并进行定量分析， 与表面形貌演化规律相比较，为数值模型的建立提供基础；
[0014] 步骤5:总结以上四步骤的结果，基于传热学方程(3)和结构演化规律得出含能材 料的热物理模型：
[0015] pC- = ?N2T + S ? (3)
[0016] 其中ρ为密度、C为比热容、T为温度、t为时间、λ为导热系数、S为化学反应放热项； // (f
[0017] 对于化学反应放热项S有以下普适表达式^ = 丁 ell .
[0018] 其中Q为单位质量的化学反应放热、α为含能材料已反应掉的质量百分数。
[0019] 本发明进一步提供一种含能材料细观热物理模型的建模方法，其特征在于:步骤1 中所述的"通过DSC法和激光导热法确定组分材料的热容和导热率"的具体方法为：
[0020] 步骤1.1采用差示扫描量热DSC法获取含能颗粒和粘合剂基体的热容
[0021]高纯氮气气氛，流速为30-50ml/min ;粘合剂基体试样预先加工成厚度为0.5-1.5mm，直径为5-7_的平整圆片型试件;含能颗粒粉末状样品直接装入样品盘中，装样品时 要轻震样品盘，以便样品间、样品与盘之间接触良好。从室温25摄氏度，按照一定的升温速 率(一般1~10摄氏度/min)加热试样，直至试样完全挥发，得到试样的DSC曲线。根据曲线计 算得到样品的热容；
[0022] 步骤1.2:采用激光导热法获取含能颗粒和粘合剂基体的导热率。
[0023] 将含能颗粒和粘合剂基体预先加工成直径为10-14mm，误差小于1.5mm;厚度0.5_ 1.5 mm，误差小于0.5 mm的平整圆片型试件。将制作好的试件放入导热仪中，按照一定的升温 速率(一般1~10摄氏度/min)加热试件，至400摄氏度，每升温50摄氏度，通过探针获取其导 热系数，每试件获取5个值，通过取平均值，得其导热率。
[0024]本发明进一步提供一种含能材料细观热物理模型的建模方法，其特征在于:步骤2 中所述的"通过DSC法和激光导热法确定组分材料的热容和导热率"的具体方法为：
[0025]步骤2.1获取含能材料的组分材料(含能颗粒和粘合剂)在四种升温速率下的DSC 曲线，四种温度分别为l°C/min、5°C/min、10°C/min、15°C/min，用于计算组分的阿伦尼乌斯 参数，包括反应活化能、指前因子、机理函数，计算方法如下：
[0026] 采用非等温动力学Kissinger方程（1)和Coats-Redfern方程（2)共同确定含能材 料组分的化学反应动力学参数和机理函数:首先从所测试样的四种升温速率的TG曲线获得 DTG峰温和反应深度值，用Kissinger方程得到反应活化能和指前因子，再通过Coats-Redfern 方程确定g (α) 函数。 ' β ' 、AR Εα
[0027] 1，1(^ 、 ρ α ρ (1)
[0028] ln(-7-，= ln^[1 / u (2)
[0029] 其中为β升温速率，^为峰温，T为在热分解过程中任一温度，Ea*表观活化能，A为 指前因子，R为气体常数，8.314J/K. mo 1，g (α)为机理函数积分形式，α为反应深度；
[0030] 步骤2.2获取含能材料组分材料(含能颗粒和粘合剂)在升温速率为5°C/min下的 TG曲线，用于对比计算材料随温度几何形貌的演化。
[0031] 本发明进一步提供一种含能材料细观热物理模型的建模方法，其特征在于:步骤3 中所述的"通过扫描电镜试验获取含能材料表面形貌随温度的演化规律，进行定量分析，得 出表面形貌随温度演化规律"的具体方法为：
[0032] 步骤3.1获取含能颗粒在不同温度下的细观表面形貌图片，为计算含能颗粒形貌 随温度演化曲线提供基础；
[0033] 步骤3.2获取粘合剂胶片在不同温度下的细观表面形貌图片，为研究粘合剂胶片 形貌随温度演化规律提供依据。
[0034]本发明进一步提供一种含能材料细观热物理模型的建模方法，其特征在于:步骤3 中所述的"通过yCT试验获取含能材料内部形貌随温度的演化规律，并进行定量分析，与表 面形貌演化规律相比较"的具体方法为：
[0035] 步骤4.1获取含能颗粒在不同温度下的细观内部形貌图片，为计算颗粒形貌随温 度演化曲线提供基础；
[0036] 步骤4.2获取含能材料在不同温度下的细观内部形貌图片，为研究含能材料形貌 随温度演化提供依据。
[0037] 本发明同现有技术相比的优越性在于:该方法充分考虑了含能材料自身损伤所特 有的特性，为含能材料细观损伤数值仿真采用的界面模型参数的选择提供了方便性和通用 性，以全面模拟预测实际的特征和过程，建立了物理模拟准则体系和满足准则的细观热物 理模型，功能全面、易于实现。
【具体实施方式】
[0038]以下将以丁羟推进剂(含能颗粒为AP颗粒，粘合剂基体为丁羟HTPB)为例，通过对 其细观热点物理模型的构建说明细观热物理模型建模方法。
[0039]步骤1:通过DSC法和激光导热法确定组分材料的热容和导热率。
[0040] 步骤1.1采用差示扫描量热DSC法获取含能颗粒和粘合剂基体的热容
[0041] 高纯氮气气氛，流速为40ml/min;粘合剂基体试样预先加工成厚度为1mm，直径为 6mm的平整圆片型试件;含能颗粒粉末状样品直接装入样品盘中，装样品时要轻震样品盘， 以便样品间、样品与盘之间接触良好。从室温25摄氏度，按照一定的升温速率(一般1~10摄 氏度/min)加热试样，直至试