Abstract: The evaluation method of ship's explosion shock resistance is challenged by some key mechanical problems, such as strong nonlinear fluid-structure coupling, large-deformation and failure evolution of solid structure. By coupling the respective advantages of peridynamics (PD) and smoothed particle hydrodynamics (SPH), an efficient PD-SPH numerical model suitable for underwater explosion shock simulations is developed. The SPH method is employed to simulate underwater shock wave propagation and fluid-structure interaction, while the PD method accurately characterizes the complete mechanical behavior of solid structures from elastic deformation to progressive damage failure. A PD-SPH numerical model is established for non-medicinal underwater shock wave loading devices. To improve the computational efficiency in large-scale simulations, a multi-GPU parallel computing framework based on domain decomposition and data-communication mechanisms is established. Model validation and parallel efficiency tests demonstrate that the proposed method can accurately predict shock wave wall pressure and target dynamic deformation, successfully reproduce typical crack propagation patterns in thin-plate structures, and simulate the entire damage process of complex grid sandwich structure. In complex fluid-structure coupling scenarios with more than 5 million particles, the actual calculation time can be compressed to nearly 1 hours. The research outcomes provide a high-precision and efficient numerical analysis tool for the design of explosion-resistant naval structures, offering significant reference value for engineering applications of fluid-structure interaction in underwater explosion problems.