This paper investigates the elastic-plastic seismic behavior of long span cable-stayed steel bridges through the plane finite-element model. Both geometric and material nonlinearities are involved in the analysis. The geometric nonlinearities come from the stay cable sag effect, axial force-bending moment interaction, and large displacements. Material nonlinearity arises when the stiffening steel girder yields. The example bridge is a cable-stayed bridge with a central span length of 605 m. The seismic response analyses have been conducted from the deformed equilibrium configuration due to dead loads. Three strong earthquake records of the Great Hanshin earthquake of 1995 in Japan are used in the analysis. These earthquake records are input in the bridge longitudinal direction, vertical direction, and combined longitudinal and vertical directions. To evaluate the residual elastic-plastic seismic response, a new kind of seismic damage index called the maximum equivalent plastic strain ratio is proposed. The results show that the elastic-plastic effect tends to reduce the seismic response of long span cable-stayed steel bridges. The elastic and elastic-plastic seismic response behavior depends highly on the characteristics of input earthquake records. The earthquake record with the largest peak ground acceleration value does not necessarily induce the greatest elastic-plastic seismic damage.