In Korea, the underground silo structure for low- and intermediate-level radioactive waste disposal facilities has been constructed and operated since 2014. Large-scale earthquakes occurred in 2016 and 2017, respectively, in Gyeongju and Pohang areas near the underground silo structures, and interest in the stability of the underground silo increased significantly. In this paper, one-dimensional free-field analyses have been carried out before the three-dimensional silo dynamic analyses subjected to earthquake loadings. As an additional study, a new form of the finite element equilibrium equation is derived in terms of relative motions, which is essentially the same equation expressed in terms of total motions where the base shear force is applied to the earthquake load. The accuracy of conventional finite element solutions is evaluated by directly comparing them with closed-form solutions by frequency domain analysis such as SHAKE91.