Abstract: Based on Harichandran-Vanmarcke coherence function model, modified Kanai-Tajimi filtered white noise power spectrum model and Amin-Ang non-stationary function model, the Fast Fourier Transform (FFT) algorithm are adopted to compile the non-stationary Artificial acceleration time history with the spatial variation of seismic ground motions which is based on Geodatis-innovated spectral representation method by using program of MATLAB. A linear acceleration integral method and the FFT frequency-domain integral method are applied to transform acceleration time history to displacement time history respectively, which will be the ground motion excitations on the bridge structure. Importing the bilinear restoring force models , and using the Combin40 element in program ANSYS to simulate the nonlinear property of materials on plastic hinge sections in the pier top and pier bottom, simultaneously taking the geometrical nonlinearity due to large displacement of piers into account, the double nonlinear Finite Element Method model of the bridge has been set up, and the nonlinearity and ductility of the plastic hinge sections of the high-rise rigid frame bridges under rear earthquake has been obtained. The approach and achievement represented in this paper have some certain significance as a reference for the rear earthquake-resistant design and research for other similar bridges.