The　method　employed　in　this　study　combines　large-scale　shaking　table　tests　with　fully　coupled　3D　non-linear　dynamic　analysis　to　investigate　the　site　response　and　structural　damage　mechanisms　of　laterally　spreading　site　cluster　pile　foundation　systems　under　weak　and　strong　seismic　excitation.　The　study　also　aims　to　verify　the　applicability　of　the　high-efficiency　and　universal　parallel　calculation　system　method,　which　involves　SolidWorks,　HyperMesh,　OpenSeesMP,　and　GiD,　with　a　fully　visualized　interface　operation　in　geotechnical　engineering,　based　on　the　test　database.　The　results　show　that　under　weak　seismic　excitation,　the　acceleration　and　displacement　in　the　free　field　and　near-pile　region　are　linearly　distributed　along　the　burial　depth,　and　the　motion　of　shallow　pile　foundation　is　mainly　controlled　by　the　inertial　effect　of　the　superstructure;　while　under　strong　seismic　excitation,　the　displacement　in　the　near-field　region　is　much　smaller　than　the　free　field　showing　significant　pile　pinning　effect,　and　the　displacement　and　acceleration　in　the　near-pile　domain　are　cosine-shaped　distributed　along　the　burial　depth,　with　the　lateral　spreading　of　site　liquefaction　intensifying.　At　the　intersection　of　clay　and　sand　layers,　a　weak　interlayer　is　present,　which　affects　the　structural　response.　The　kinematic　effect　of　the　site　dominates　the　response,　and　the　location　of　vulnerability　gradually　shifts　from　the　pile　bottom　to　the　pile　head.　©　2023,　The　Author(s),　under　exclusive　licence　to　Springer　Nature　B.V.