Development of flutter constraints for high-fidelity aerostructural optimization

High fidelity computational modeling and optimization of aircraft has the potential to allow engineers to produce more efficient designs requiring fewer unforeseen design modifications late in the design process. In order for the optimization algorithm to generate a useful design, all the relevant physics must be considered, including flutter. This is especially important for the high-fidelity aerostructural optimization of commercial aircraft, which is likely to result in wing designs that are prone to flutter. To address this issue, we developed a flutter constraint formulation suitable for gradient-based optimization. This paper investigates the feasibility of using a Doublet-Lattice Method (DLM) based flutter constraint for high-fidelity aerostructural optimization. The p-k flutter equation is solved using a determinant iterative method to obtain the eigenvalues. The Kreisselmeier–Steinhauser (KS) function is used to aggregate the damping values for individual modes into a single value that is used as the flutter constraint. To study the behavior of the flutter constraints using this method, we optimize a simple flat plate problem and perform a flutter analysis for a full transport aircraft using the uCRM configuration. We compute accurate and efficient derivatives for the DLM and the coupled derivatives with respect to structural sizing variables, as well as wing shape variables. These derivatives are computed using an automatic differentiation method and validated using the complex-step method. However, it is found that the current formulation using determinant iteration root finding method, is not adequate or robust, even for the simplest problems and reformulation is required.

[1]  Valéerie Pascual,et al.  Extension of TAPENADE toward Fortran 95 , 2006 .

[2]  David M. Schuster,et al.  Computational Aeroelasticity: Success, Progress, Challenge , 2003 .

[3]  C. Mader,et al.  Derivatives for Time-Spectral Computational Fluid Dynamics Using an Automatic Differentiation Adjoint , 2012 .

[4]  Carol D. Wieseman,et al.  Aeroelastic Tailoring of Transport Wings Including Transonic Flutter Constraints , 2015 .

[5]  Joaquim R. R. A. Martins,et al.  A parallel aerostructural optimization framework for aircraft design studies , 2014 .

[6]  M. Gibbons,et al.  Aeroelastic Calculations Using CFD for a Typical Business Jet Model , 1996 .

[7]  Joaquim R. R. A. Martins,et al.  A CAD-Free Approach to High-Fidelity Aerostructural Optimization , 2010 .

[8]  Joaquim R. R. A. Martins,et al.  pyOpt: a Python-based object-oriented framework for nonlinear constrained optimization , 2011, Structural and Multidisciplinary Optimization.

[9]  Thomas W. Sederberg,et al.  Free-form deformation of solid geometric models , 1986, SIGGRAPH.

[10]  Joaquim R. R. A. Martins,et al.  Extensions to the design structure matrix for the description of multidisciplinary design, analysis, and optimization processes , 2012, Structural and Multidisciplinary Optimization.

[11]  Joaquim R. R. A. Martins,et al.  Aerostructural optimization of the common research model configuration , 2014 .

[12]  Michael A. Saunders,et al.  SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization , 2002, SIAM J. Optim..

[13]  K. Bathe Finite Element Procedures , 1995 .

[14]  Her Mann Tsai,et al.  Calculation of Wing Flutter by a Coupled Fluid-Structure Method , 2001 .

[15]  John C. Vassberg,et al.  Influence of Shape Parameterization on Aerodynamic Shape Optimization , 2014 .

[16]  Laurent Hascoët,et al.  TAPENADE 2.1 user's guide , 2004 .

[17]  Joaquim R. R. A. Martins,et al.  RANS-based aerodynamic shape optimization investigations of the common research modelwing , 2014 .

[18]  William P. Rodden,et al.  Further Refinement of the Subsonic Doublet-Lattice Method , 1998 .

[19]  Rudolph N. Yurkovich,et al.  Status of Unsteady Aerodynamic Prediction for Flutter of High-Performance Aircraft , 2003 .

[20]  John C. Vassberg,et al.  Development of a Common Research Model for Applied CFD Validation Studies , 2008 .

[21]  Joaquim R. R. A. Martins,et al.  High-fidelity aerostructural optimization considering buffet onset , 2015 .

[22]  Peter D. Dunning,et al.  Optimal Topology of Aircraft Rib and Spar Structures Under Aeroelastic Loads , 2014 .

[23]  Joaquim R. R. A. Martins,et al.  Multipoint High-Fidelity Aerostructural Optimization of a Transport Aircraft Configuration , 2014 .

[24]  G. Kreisselmeier,et al.  SYSTEMATIC CONTROL DESIGN BY OPTIMIZING A VECTOR PERFORMANCE INDEX , 1979 .

[25]  Max Blair,et al.  A Compilation of the Mathematics Leading to the Doublet Lattice Method , 1992 .

[26]  W. Rodden,et al.  A doublet-lattice method for calculating lift distributions on oscillating surfaces in subsonic flows. , 1969 .

[27]  Nancy Wilkins-Diehr,et al.  XSEDE: Accelerating Scientific Discovery , 2014, Computing in Science & Engineering.

[28]  Laurent Hascoët,et al.  The Tapenade automatic differentiation tool: Principles, model, and specification , 2013, TOMS.

[29]  Andreas Griewank,et al.  Evaluating derivatives - principles and techniques of algorithmic differentiation, Second Edition , 2000, Frontiers in applied mathematics.

[30]  John W. Edwards,et al.  AIAA 98-2421 An Overview of Recent Developments in Computational Aeroelasticity , 1998 .

[31]  Joaquim R. R. A. Martins,et al.  An adaptive approach to constraint aggregation using adjoint sensitivity analysis , 2007 .

[32]  Joaquim R. R. A. Martins,et al.  Parallel Solution Methods for Aerostructural Analysis and Design Optimization , 2010 .

[33]  Hermann J. Hassig,et al.  An Approximate True Damping Solution of the Flutter Equation by Determinant Iteration , 1971 .

[34]  Gregory A. Wrenn,et al.  An indirect method for numerical optimization using the Kreisselmeir-Steinhauser function , 1989 .

[35]  Joaquim R. R. A. Martins,et al.  The complex-step derivative approximation , 2003, TOMS.

[36]  Joaquim R. R. A. Martins,et al.  A Scalable Parallel Approach for High-Fidelity Aerostructural Analysis and Optimization , 2012 .