1.一种含能材料冲击起爆的计算机模拟方法,其特征在于,包括以下步骤;
步骤1、用Materials Studio软件建立含能材料冲击起爆模拟体系的三维模型,并导出所有原子的三维坐标;
步骤2、确定模拟参数,包括初始温度、初始压力、迭代时间、时间步长和冲击速度,模拟计算并输出数据,具体包括:
101、根据玻尔兹曼分布随机抽样获得扰动初始速度;
102、步骤101获得原子的初始速度后,采用基于紧束缚近似的密度泛函SCC-DFTB方法计算迭代的第n步所有原子间的作用力来计算原子的速度;
103、完成步骤102计算所有原子间的作用力后,根据原子速度计算公式计算第n+1/2步所有原子的速度;
104、根据公式ri=ri+viΔt计算第n+1步的三维坐标位置,其中ri表示第i个粒子的位置,vi为上一步计算得到的速度,Δt为设定的时间步长;
105、再次根据根据原子速度计算公式计算第n+1/2步所有原子的速度;
106、根据分子动力学模拟体系的哈密顿量、能量关系式以及约束压力的方程计算并输出体系的压力、体积、能量和温度。
2.根据权利要求1所述的含能材料冲击起爆的计算机模拟方法,其特征在于,步骤101中扰动初始速度为:令初始位置在差分划分网格的格子上,初始速度则从玻尔兹曼分布随机抽样得到。