In the seismic response analysis of large-scale reinforced concrete (RC) structures, such as nuclear power plant buildings, while considering the soil-structure interaction (SSI) using high-fidelity models, the application of high-performance computing (HPC) is indispensable owing to the number of the degrees of freedom (DOF). This study shows the development of a seismic response analysis method for large-scale RC structures with FEM utilizing HPC. Herein, we implemented a previously proposed constitutive relation of concrete to the HPC-FEM program, FrontISTR. In the HPC-FEM program, conjugate gradient method is usually used as the solver algorithm that requires a positive-definiteness of the stiffness matrix. However, the constitutive relation of concrete expresses the material damage by stress-strain curve softening; hence one cannot assure that the stiffness matrix is positive-definite. Several techniques have been proposed to overcome this limitation. Herein, we adopted the method to calculate the stiffness matrix using elastic tensor of the material by comparing the speed of computation. We applied the developed software to a seismic response analysis of a nuclear power plant with surrounding ground, which is approximately 150 million-DOF problem. This paper shows the results and discusses the analysis performance.