An integrated numerical method is proposed for analyzing the nonlinear seismic response of near-fault building clusters, comprising three algorithms: (1) a structural investigated lump algorithm for elastoplastic dynamic response of structure; (2) a connecting investigated lump algorithm for bidirectional wave propagation between the site and elastoplastic building clusters; (3) a geomedia investigated lump algorithm for seismic wave propagation with an improved viscoelastic constitutive model, which allows independent definition of P/S-wave quality factors to characterize geomedia attenuation. Validated for its capability in simulating site-city dynamic interaction problems via a shaking table test, the method is applied to study the seismic response of near-fault building clusters in Xichang City under a hypothetical Mw6.8 earthquake. It is shown that irrespective of whether shallow geological structures are considered, clusters (c2–c4) situated in rupture-forward surface area within ~1.5 km of the fault trace entered the elastoplastic stage, while others (c1, c5) remained elastic. Shallow geological structures may reverse locally hanging-wall/footwall effects of both near-fault structural seismic response and ground motion. A notable seismic-response characteristic of near-fault structures undergoing the elastoplastic stage is that the permanent structural motion displacement (PSMD) at the slab of a specific floor incorporates not only the non-zero permanent ground motion displacement (PGMD) but also the non-zero final structural residual displacement (FSRD) relative to the supporting ground. The developed method could provide support for seismic damage assessment, site selection, and structural optimization design of near-fault building clusters.
Zhong et al. (2026) studied this question.