个简化的MD程序结构 ●初始化: 读入运行参数(初始参数,时间步数,密度,粒子数, 时间步长) 生成或读取初始配置(粒子的位置和速度) ●循环时间步增量,直到t=t最后 一计算粒子受力 对运动方程积分 一若1>1甲断系统达到平衡状态 ●结果输出一个简化的MD程序结构 t t = 最后 t t > 平衡 z初始化： －读入运行参数（初始参数，时间步数，密度，粒子数， 时间步长） －生成或读取初始配置（粒子的位置和速度） －计算粒子受力 －对运动方程积分 －若 z结果输出 ，系统达到平衡状态 z循环时间步增量，直到