In the aeronautical industry, the robustness properties of an aircraft are usually assessed using intensive and time-consuming simulations. Fortunately, several optimization techniques can be implemented and applied for clearance of flight control laws, in order to improve the efficiency and reliability of the certification process. Some of them, such as Lyapunov-based analysis and μ-analysis, require the considered models to be written as Linear Fractional Representations (LFR), which is unfortunately not the case of the aeroelastic models available in an industrial context. A whole methodology is thus proposed in this paper to convert a set of numerical flexible aircraft models into a suitable LFR which depends on the aircraft configuration and the flight conditions. This is a challenging issue, since the size of the initial models is very large and the state vector does not have the same physical meaning for the whole model set. Nevertheless, an LFR is obtained, which is representative of the aircraft behavior in the sense that its eigenvalues and frequency responses almost exactly match those of the initial flexible models. Moreover, its complexity proves compatible with the use of robustness analysis tools. Numerical results indeed show that the considered clearance problem is computationally tractable., clearance of flight control laws, LFT modeling and reduction techniques, aeroelastic models, robustness analysis.