You need something much more sophisticated than a general purpose piping program. Basic 3D beam theory, with discrete restriants (simulating a continuous beam on an elastic foundation) has too many course approximations to yield reliable results in this situation. FEA is the way to go, assuming you can model (mesh) the soil also.