固体推进剂和高分子共混物的微观、介观和宏观多尺度模拟研究

固体推进剂和高分子共混物的微观、介观和宏观多尺度模拟研究随着计算机运算能力的不断提高，计算机模拟在材料设计和开发过程中的作用越来越重要，其可以揭示一些无法或很难从实验获得的微观机理和本质。高分子材料的特性是具有多个时间和空间尺度，到目前为止没有一种单一的模拟方法可以跨越多个时间及空间尺度。

本文采用分子动力学（MD）、介观动力学（MesoDyn）、耗散粒子动力学（DPD）和有限元分析方法（FEM）对几种高分子材料的结构与性能进行了研究，主要包括：高分子粘结剂/增塑剂、二元高分子共混物和高分子/粘土纳米复合材料。固体推进剂和塑性炸药（RBX）的力学性能在很大程度上依赖于配方中高分子粘结剂与增塑剂的相容性。

为预测这两种材料的相容性，本文对高分子粘结剂/增塑剂共混物进行了MD 和MesoDyn模拟研究，即在COMPASS力场条件下，对端羟基聚丁二烯（HTPB）/癸二酸二辛酯（DOS）和HTPB/硝化甘油（NG）共混物的的密度（ρ）、内聚能密度（CED）、溶度参数（δ）、Flory–Huggins相互作用参数（χ）、玻璃化转变温度（Tg）和力学性能等进行了模拟计算。结果表明通过比较溶度参数差值（Δδ）的大小、分子间径向分布函数g(r)或模拟前后体系密度变化情况均可以推测HTPB/DOS属于相容体系，而HTPB/NG则不然。

通过分析温度-体积曲线得到HTPB、HTPB/DOS和HTPB/NG的T_g分别为197.54 K、176.30 K和200.03 K，而HTPB/DOS中只出现一个T_g的情况则进一步表明其为相容体系；另外，通过MD模拟也可以得到共混物的弹性模量（E），体积模量（K）和剪切模量（G），添加DOS增塑剂后使E、K和G显著下降，柔性增强，共混物的力学性能得到有效改善。将通过MD模拟得到的Δδ转化为

MesoDyn模拟的输入参数，然后采用MesoDyn模拟方法对共混体系的介观形貌与动力学演变过程进行了计算。

计算所得的等密度图、有序度参数和自由能密度等也可用于判断共混体系的相容性。MD和MesoDyn模拟结果均表明：HTPB/DOS属于相容体系，而HTPB/NG 属于不相容体系，其结论与现有实验结果一致。

（2）通过MD、MesoDyn和FEM等模拟计算方法研究了聚丙烯（PP）/尼龙1（1PA11）和聚乳酸（PLA）/聚对苯二甲酸乙二醇酯（PET）共混物的相容性、介观结构和力学性能。分别建立了五个质量比例为10/90、30/70、50/50、70/30和90/10的PP/PA11和PLA/PET共混物。

通过MD模拟得到纯物质的δ与文献的实验值比较一致，通过比较模拟得到的χ与临界值χc的大小以及C-C原子对分子间g(r)的大小可以确定在90/10比例下PP/PA11共混具有一定相容性，其它比例下的PP/PA11共混物不相容，而PLA在任意比例下均与PET完全互溶。PP/PA11共混物出现两个T_g，每个值分别对应于各自组分的T_g，而PLA/PET共混物只有一个T_g，通过共混物的玻璃化转变过程也可以证明共混体系的相容与否。

为了进一步探究PP/PA11和PLA/PET共混物的介观结构，将MD模拟的组分间相互作用的χ参数转化为MesoDyn模拟的相互作用参数，采用MesoDyn模拟方法在介观水平研究了共混物的相分离动力学过程，共混物的相容性可以从其介观形态结构得到进一步证实。将通过MesoDyn模拟得到的周期性模拟盒子中基于网格的两相浓度分布转化为宏观模拟中FEM分析的输入结构，通过FEM分析得到共混体系的力学性能和局部应力分布情况。

结果表明：所有的共混物均为各向同性材料，PP/PA11共混物的模量（E，K

和G）随PP含量的增加线性降低，PLA/PET共混物的E和G随PLA含量的增加而增加，K随PLA含量的增加而降低，FEM预测的值与实验测试值吻合较好。（3）通过MD、DPD和FEM等计算模拟方法对PA11/季铵盐（Quat）/蒙脱土（MMT）纳米复合材料的微观分子结构、介观结构和宏观力学性能进行了研究。

MD模拟得到的平衡构象表明：Quat分子平铺在MMT上，几乎覆盖了MMT的整个表面，高分子链塌缩在季铵盐分子上，没有直接与MMT表面相接触。Quat 与MMT之间存在强烈的相互作用，Quat中的极性基团通过强烈的静电吸引保持在邻近的MMT层上，Quat中的非极性基团通过范德华（vdW）力与高分子链发生作用，Quat与PA11之间的结合能也比较强。

通过将MD模拟得到的非键相互作用能转化为介观模拟中的相互作用参数，采用DPD模拟方法研究了共混物的介观形态结构，DPD模拟得到的介观形貌结构与MD模拟结果一致。将通过DPD模拟得到的介观结构作为宏观模拟中FEM分析的输入结构，采用FEM方法预测了共混物的力学性能。

FEM分析结果表明，共混物为各向异性材料，在z方向的模量值与实验值一致，且比垂直与z方向（x和y方向）的低很多。