如何用SPH程序来计算？计算的思路是什么？

1用粒子描述问题域，再对边界条件进行适当处理则可以得到初始瞬时t时刻离散粒子上的计算信息；2用适当的粒子近似法进行离散控制方程的导数或者积分；3从得到的速度或者位置坐标计算应变和应变率，然后计算每个粒子在瞬时t的应力；4用计算出来的应力计算每个粒子的加速度；5用瞬时t的加速度计算t+Δt时刻的速度和位置新坐标；6用新速度或者新位置坐标计算t+Δt瞬时新的应变率和新应变，然后计算t+Δt瞬时的新应力，重复第四、第五和第六步，直到指定的瞬间。