激波与固体颗粒群的相互作用问题是超音速气固两相流研究领域的一个重要课题，其研究成果在航空航天、医疗卫生、安全防控等领域都有着重要的应用。所以对于激波与颗粒群相互作用机理的进一步研究有助于超音速气固两相流理论的发展和完善以及工业生产相关技术的改进。　　本文基于ICEM、Fluent等软件组成的CFD计算平台，以及Origin、Tecplot等后处理软件，数值模拟了激波加载下三种周期性边界单元体模型球阵的周围流场和非稳态阻力。通过对压力云图和非稳态阻力系数变化曲线的同步对比，揭示了单元体模型球阵非稳态阻力的产生机理和变化规律。　　本文研究的主要内容和结论如下:　　1、对单元体内不同球的阻力系数曲线进行对比分析，发现单元体内阻力系数变化规律具有稳定性、相似性，证明了将阻力系数模型扩展到颗粒群中的可行性。　　2、对于目标单元体中的单个球，其阻力系数变化规律受到入射激波、前排球衍射激波、同排球和后排球反射激波以及激波-激波相互作用后形成的各种波形结构的影响。所以要比单球、双球、四球等模型受力复杂得多。阻力系数变化规律为先急剧增到到峰值，在单调下降到一个谷值，随后继续增到到第二个峰值，再波动下降出现一个最小值波谷，最后在0附近震荡并逐渐趋于稳定。　　3、在入射激波与单元体相互作用阶段（大概对应无量纲时间t=1～3），入射激波阵面与单元体内球的相对位置决定了阻力系数的变化趋势，即当激波阵面处于球前驻点到赤道位置时阻力系数上升，当激波阵面处于赤道到后主点位置时阻力系数下降。　　4、对于单元体A，阻力系数峰值随着空隙比（无量纲间距）的增大而减小，并且空隙比越大，在入射激波与单元体相互作用结束之后阻力系数越容易趋于稳定。　　5、对于A、B、C三种单元体来说，阻力系数都随着马赫数的增大而减小，且曳马赫数越大，入射激波与单元体相互作用结束后阻力系数越容易趋于稳定。　　6、在模型球阵中的前后排球（垂直与激波来流方向）对周围流场和激波结构的影响要比同排相邻球大。　　7、在所有工况条件下，单元体阻力系数变化曲线的第二波峰都高于第一波峰，造成这现象的主要原因是因为激波的二次反射。