Based on Large-displacement Non-linear Elastic Generalized Variational Principle, coupling effect of axial and flexural action, shearing strain energy, torsional strain energy of stiffening girder were considered, the large-displacement incomplete generalized potential energy functional of space coupling free vibration of a three-span self-anchored cable-stayed suspension bridge was presented. By constraint variation, fundamental differential equations of vertical flexural vibration, lateral flexural vibration, longitudinal vibration and torisional vibration were formulated, also presented the equations for the main tower with respect to longitudinal and lateral vibration. The linear free vibration differential equation was obtained when the nonlinear items were discarded. This approach provides theoretical basis for analysis of natural vibration character of self-anchored cable-stayed suspension bridges.