This paper is about an application of an energy approach to computational aeroelasticity. A frequency-domain calculation of aerodynamic work is presented in a form that has not been previously discussed. In large aeroelastic systems such as aircraft flutter models, the obtained expression allows us to quantify the roles played in flutter by the generalized coordinates, the phases between them, and the generalized aerodynamic forces. This is exemplified in a body-freedom flutter analysis of a flying-wing aircraft: the X-56A. Another interesting feature that is proposed is a diagram of the aerodynamic work as a function of airspeed. Such functions allow an observance of an evolution with the airspeed of the terms that are most reflective of the aeroelastic stability changes: for example, a phase between the two dominant modes. These diagrams can complement typical flutter trends in analysis documentation and aid in flutter suppression. The approach also permits a perspective on flutter as an interaction of aircraft surfaces rather than vibrational modes. Once a sensitivity of an aeroelastic eigenvalue to a surface area is measured with the presented approach, it can be used in the aircraft design to mitigate flutter and other undesirable aeroelastic responses associated with lightly damped eigenvalues. An illustration of this idea is provided. Finally, energy-based computations allow posing energy-efficient active flutter suppression problems. This has been presented before in the literature. Examples of this aspect are made here with aircraft models (for the first time, as far as the author knows).