The　saturated　soil　site-pile　foundation-superstructure　system,　with　large　degrees　of　freedom　or　strong　nonlinear　problems,　often　involves　a　large　amount　of　calculation　and　low　computational　efficiency.　In　this　paper,　a　fully　explicit　finite　difference　method　is　proposed　for　a　saturated　soil　site-pile　foundation-superstructure　system.　Since　the　proposed　method　has　the　advantages　of　decoupling　in　both　time　and　space,　it　does　not　need　to　solve　the　equations　simultaneously,　which　grants　it　high　computational　efficiency.　At　the　same　time,　the　method　is　implemented　on　the　open-source　software　OpenSees　and　is　used　to　compare　and　analyze　the　dynamic　responses　of　the　shaking　table　of　a　liquefiable　soil　site-pile　foundation-superstructure　system.　After　the　calculation　and　analysis,　the　numerical　solutions　were　found　to　be　in　good　agreement　with　the　experimental　solutions,　which　verifies　the　proposed　method　and　illustrates　that　the　proposed　method　can　reasonably　simulate　the　seismic　responses　of　the　whole　system.　In　addition,　the　proposed　calculation　platform　in　OpenSees　can　be　used　for　the　analysis　of　the　liquefaction　process　and　the　possible　large　deformation　of　soil　after　liquefaction,　as　well　as　for　analyzing　the　failure　mode　of　the　complex　and　nonlinear　saturated　soil　sites　and　structures　under　the　effects　of　an　earthquake.