如题用abaqus中的sph方法算了好久了，一直求不出砰击力时程曲线。如果能求出楔形体表面某点的局部砰击压力就更好了。目前是通过修改inp的方法来生成sph粒子的。楔形体入水约束除下落方向的所有自由度，入水后虽然能够观察到位移和应力变化，但是表面和水体sph粒子的压力一直不知道输出方法。查了一下网上的，基本是要用固壁求解压力的方法来避免sph方法的负压问题，如果要这样求该怎么办呢，希望能够有哪位专家可以解惑。