[论文解读] Micromechanics-based prediction of thermoelastic properties of high energy materials

B. Banerjee, Daniel O. Adams|arXiv (Cornell University)|Jan 11, 2012

Composite Material Mechanics参考文献 11被引用 13

一句话总结

本文提出一种基于细观力学的方法，结合递归单元法（RCM）与广义细胞法（GMC），用于预测高能材料（如PBX 9501）的热弹性性能。PBX 9501是一种高颗粒体积分数（92%）且剪切模量对比极端（高达20,000）的聚合物粘结炸药。RCM能够以计算高效且高度准确的方式预测结果，接近实验值；而GMC由于未能有效模拟应力传递及法向-剪切响应的耦合效应，显著低估了刚度。

ABSTRACT

High energy materials such as polymer bonded explosives are commonly used as propellants. These particulate composites contain explosive crystals suspended in a rubbery binder. However, the explosive nature of these materials limits the determination of their mechanical properties by experimental means. Therefore micromechanics-based methods for the determination of the effective thermoelastic properties of polymer bonded explosives are investigated in this research. Polymer bonded explosives are two-component particulate composites with high volume fractions of particles (volume fraction $>$ 90%) and high modulus contrast (ratio of Young's modulus of particles to binder of 5,000-10,000). Experimentally determined elastic moduli of one such material, PBX 9501, are used to validate the micromechanics methods examined in this research. The literature on micromechanics is reviewed; rigorous bounds on effective elastic properties and analytical methods for determining effective properties are investigated in the context of PBX 9501. Since detailed numerical simulations of PBXs are computationally expensive, simple numerical homogenization techniques have been sought. Two such techniques explored in this research are the Generalized Method of Cells and the Recursive Cell Method. Effective properties calculated using these methods have been compared with finite element analyses and experimental data.

研究动机与目标

• 开发计算高效的细观力学方法，用于预测高能材料（如PBX 9501）的有效热弹性性能。
• 解决传统细观力学方法在处理高颗粒体积分数（>90%）与极端模量对比（5,000–10,000）时的局限性。
• 将数值预测结果与PBX 9501的实验数据进行验证，尤其关注其弹性模量与热膨胀系数。
• 评估RCM与GMC在有限元法（FEM）模拟与实验基准下的性能表现。
• 识别并修正现有均质化技术中的缺陷，特别是GMC对刚度分量的严重低估问题。

提出的方法

• 采用递归单元法（RCM），一种多级均质化技术，通过递归地将代表性体积元（RVE）划分为更小的块体，以高效模拟高对比度复合材料。
• 应用广义细胞法（GMC），一种数值均质化方案，用于计算有效性能，但识别出其在捕捉应力传递与耦合效应方面的局限性。
• 采用有限元法（FEM）并使用256×256结构化网格模拟PBX 9501的RVE，其中颗粒为圆形，结果与RCM和GMC进行对比。
• 基于干混颗粒分布生成五个具有统计代表性的微观结构，以评估变异性与对称性。
• 计算有效刚度矩阵（C*），包括C₁₁、C₂₂、C₁₂与C₆₆，并与实验值及理论边界进行比较。
• 通过计算成本、精度与不同RVE尺寸及离散化水平下的收敛性，评估方法性能。

实验结果

研究问题

• RQ1 与实验数据相比，递归单元法（RCM）在预测PBX 9501的有效弹性模量与热膨胀系数方面有多准确？
• RQ2 为何广义细胞法（GMC）在高模量对比与高颗粒体积分数条件下无法准确预测PBX 9501的刚度？
• RQ3 对具有圆形颗粒的RVE进行有限元模拟，在多大程度上能反映真实PBX 9501的各向同性行为与有效性能？
• RQ4 不同微观结构实现下的预测有效性能如何变化？RVE尺寸与网格分辨率的影响是什么？
• RQ5 RCM是否能在显著降低计算成本的前提下，实现与全FEM模拟相当的高精度，适用于高对比度颗粒复合材料？

主要发现

• FEM模拟结果较PBX 9501的实验值（1.60 GPa）高估了法向刚度分量（C₁₁与C₂₂）约15–25%。
• RCM分别预测C₁₁与C₂₂为3.20 GPa与2.99 GPa，较FEM结果高出约60%，且为实验值的两倍，但仍在高对比度复合材料可接受范围内。
• GMC严重低估了法向刚度分量（C₁₁与C₂₂），仅预测为0.15 GPa，不足实验值的10%，原因在于对应力传递效应建模不佳。
• RCM预测C₁₂为0.54 GPa，与FEM结果一致；而GMC仅预测0.12 GPa，约为FEM值的1/7，表明其在捕捉面内耦合效应方面完全失败。
• FEM对剪切刚度（C₆₆）的预测值为0.51 GPa，与实验值相差在10%以内；RCM低估60%（0.27 GPa）；GMC预测值仅为0.01 GPa，接近Hashin-Shtrikman下界。
• RCM在计算成本仅为FEM约三分之一的情况下，实现了可接受误差范围内的结果，展现出在高对比度复合材料中卓越的效率与精度。