Numerical simulation using finite element (FE) software packages has been widely adopted to investigate the seismic performance of large-scale structures, recently.Among these, Open System for Earthquake Engineering Simulation (OpenSees), as an open-source FE software program, has increasingly become one of the most influential packages.However, the computational efficiency of the matrix solvers currently cannot satisfy the demand for numerical simulation of large-scale structures.Consequently, two new, parallel, iterative linear solvers are proposed and implemented in OpenSees based on two graphic processing unit-based (GPU-based) libraries.The time history analysis of two typical large-scale structures are performed using the proposed solvers.The speed-up ratio of the proposed solvers are up to 9 to 15 (achieved with high accuracy in results).The research outcome can provide an effective computing technology for further numerical analysis of the seismic behavior of large-scale structures.