AbstractA general substructure method for analysis of response of structures to earthquake ground motion, including the effects of structure‐soil interaction, is presented. The method is applicable to complex structures idealized as finite element systems and the soil region treated as either a continuum, for example as a viscoelastic halfspace, or idealized as a finite element system. The halfspace idealization permits reliable analysis for sites where essentially similar soils extend to large depths and there is no rigid boundary such as soil‐rock interface. For sites where layers of soft soil are underlain by rock at shallow depth, finite element idealization of the soil region is appropriate; in this case, the direct and substructure methods would lead to equivalent results but the latter provides the better alternative. Treating the free field motion directly as the earthquake input in the substructure method eliminates the deconvolution calculations and the related assumption—regarding type and direction of earthquake waves—required in the direct method.Spatial variations in the input motion along the structure‐soil interface of embedded structures or along the base of long surface supported structures are included in the formulation.The substructure method is computationally efficient because the two substructures—the structure and the soil region—are analysed separately; and, more important, it permits taking advantage of the important feature that response to earthquake ground motion is essentially contained in the lower few natural modes of vibration of the structure