Abstract: Three-dimensional dynamic analysis of an open excavation subway station under earthquake is carried out. The space of subway station is complex, and the transverse section is variable. Major earthquake analysis must be performed. Soil-structure dynamic interaction is also need to be considered. The subway station is modeled and analyzed by means of the software ANSYS and ABAQUS respectively. The ANSYS model is introduced into ABAQUS based on a programming interface coded by the authors. A three-dimensional viscoelastic dynamic constitutive model is used to simulate the soils. Three-dimensional viscoelastic artificial boundary is applied. Seismic wave of bedrock is input as the acceleration. The results demonstrate that edge concrete fiber at the bottom of interior column enters into the plastic stage, and local damage but no collapse damage is found. The proposed method may help engineers in similar projects.