A new numerical method for large-sized buildings seismic resistance modeling is proposed. It assumes various factors: contact interaction with soil, gravity field, ground nonhomogeneity and earthquake hypocenter location. Deep and shallow earthquakes modeling methods are proposed additionally. It is shown that the location of the seismic source has different effects on the behavior of different types of structures during earthquakes.