A finite element aerodynamic model that can be used to analyse flutter instability of long span bridges in the time domain is presented. This approach adopts a simplified quasi-steady formulation of the wind forces neglecting the vortex shedding effects. The governing equations used are effective only for reduced velocities $V^*$ sufficiently great: this is generally acceptable for long-span suspension bridges and, then, the dependence of the wind forces expressions of the flutter derivatives can be neglected. The procedure describes the mechanical response in an accurate way, taking into account the non-linear geometry effects (large displacements and large strains) and considering also the compressed locked coil strands instability. The time-dependence of the inertia force due to fluid structure interaction is not considered. The numerical examples are performed on the three-dimensional finite element model of the Great Belt East Bridge (DK). A mode frequency analysis is carried out to validate the model and the results show good agreement with the experimental measurements of the full bridge aeroelastic model in the wind tunnel tests. Significant parameters affecting bridge response are introduced and accurately investigated.