The reasons of nonlinear effects of cable-stayed bridges are investigated. A numerical simulation method for nonlinear behavior is developed in this work. By considering the girder-tower coupling effect, a new beam element is proposed. Based on the modified Lagrange method, a group of nonlinear finite element equations under earthquake excitation is built. The corresponding computer software is developed. The linear and nonlinear responses for Yueyang Dongtinghu cable-stayed Bridge have been computed. The results are discussed and compared. It can be concluded that the developed methodology ma...