土结相互作用分析(群桩),通过自编FORTRAN根据NEWMARK隐式积分格式实现