Design optimization of dynamic flexible multibody systems using the discrete adjoint variable method

Abstract The design space of dynamic multibody systems (MBSs), particularly those with flexible components, is considerably large. Consequently, having a means to efficiently explore this space and find the optimum solution within a feasible time-frame is crucial. It is well-known that for problems with several design variables, sensitivity analysis using the adjoint variable method extensively reduces the computational costs. This paper presents the novel extension of the discrete adjoint variable method to the design optimization of dynamic flexible MBSs. The extension involves deriving the adjoint equations directly from the discrete, rather than the continuous, equations of motion. This results in a system of algebraic equations that is computationally less demanding to solve compared to the system of differential algebraic equations produced by the continuous adjoint variable method. To describe the proposed method, it is integrated with a numerical time-stepping algorithm based on geometric variational integrators. The developed technique is then applied to the optimization of MBSs composed of springs, dampers, beams and rigid bodies, considering both geometrical (e.g., positions of joints) and non-geometrical (e.g., mechanical properties of components) design variables. To validate the developed methods and show their applicability, three numerical examples are provided.

[1]  Peter Betsch,et al.  Constrained integration of rigid body dynamics , 2001 .

[2]  L. Petzold,et al.  Adjoint sensitivity analysis for differential-algebraic equations: algorithms and software☆ , 2002 .

[3]  R. Y. Yakoub,et al.  Three Dimensional Absolute Nodal Coordinate Formulation for Beam Elements: Theory , 2001 .

[4]  Andrei Serban Schaffer,et al.  On the adjoint formulation of design sensitivity analysis of multibody dynamics cs , 2005 .

[5]  Adrian Sandu,et al.  Dynamic Response Optimization of Complex Multibody Systems in a Penalty Formulation using Adjoint Sensitivity , 2014, ArXiv.

[6]  L. Petzold,et al.  Sensitivity analysis of differential-algebraic equations: A comparison of methods on a special problem ✩ , 2000 .

[7]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[8]  R. Y. Yakoub,et al.  Three Dimensional Absolute Nodal Coordinate Formulation for Beam Elements: Implementation and Applications , 2001 .

[9]  J. Mayo,et al.  Describing Rigid-Flexible Multibody Systems Using Absolute Coordinates , 2003 .

[10]  Thomas Lauß,et al.  The discrete adjoint method for parameter identification in multibody system dynamics , 2017, Multibody System Dynamics.

[11]  D. García-Vallejo,et al.  Three-dimensional formulation of rigid-flexible multibody systems with flexible beam elements , 2008 .

[12]  L. Petzold,et al.  Numerical methods and software for sensitivity analysis of differential-algebraic systems , 1986 .

[13]  Li-Qun Chen,et al.  Second order adjoint sensitivity analysis of multibody systems described by differential–algebraic equations , 2007 .

[14]  Haiyan Hu,et al.  Structural optimization of flexible components in a flexible multibody system modeled via ANCF , 2016 .

[15]  Alexander Verbart,et al.  The equivalent static loads method for structural optimization does not in general generate optimal designs , 2018 .

[16]  Liping Chen,et al.  First order sensitivity analysis of flexible multibody systems using absolute nodal coordinate formulation , 2012 .

[17]  J. Marsden,et al.  DISCRETE MECHANICS AND OPTIMAL CONTROL: AN ANALYSIS ∗ , 2008, 0810.1386.

[18]  J. Marsden,et al.  Variational integrators for constrained dynamical systems , 2008 .

[19]  Peter Betsch,et al.  A rotationless formulation of multibody dynamics: Modeling of screw joints and incorporation of control constraints , 2009 .

[20]  D. Tortorelli,et al.  Structural optimization of multibody system components described using level set techniques , 2015 .

[21]  Michaël Bruyneel,et al.  Discussion on the optimization problem formulation of flexible components in multibody systems , 2013 .

[22]  Johannes Gerstmayr,et al.  Analysis of Thin Beams and Cables Using the Absolute Nodal Co-ordinate Formulation , 2006 .

[23]  Edward J. Haug,et al.  Second‐order design sensitivity analysis of mechanical system dynamics , 1982 .

[24]  J. Marsden,et al.  Variational Integrators and the Newmark Algorithm for Conservative and Dissipative Mechanical Systems , 2000 .

[25]  J. G. Jalón,et al.  Dynamic Analysis of Three-Dimensional Mechanisms in “Natural” Coordinates , 1987 .

[26]  A. Mikkola,et al.  Review on the Absolute Nodal Coordinate Formulation for Large Deformation Analysis of Multibody Systems , 2013 .

[27]  Edward J. Haug,et al.  Implicit Numerical Integration for Design Sensitivity Analysis of Rigid Multibody Systems , 2005 .

[28]  J. G. Jalón.  Twenty-five years of natural coordinates , 2007 .

[29]  Mathieu Desbrun,et al.  Discrete geometric mechanics for variational time integrators , 2006, SIGGRAPH Courses.

[30]  R. Haftka,et al.  Computational aspects of sensitivity calculations in transient structural analysis , 1989 .

[31]  Rudranarayan M. Mukherjee,et al.  A divide-and-conquer direct differentiation approach for multibody system sensitivity analysis , 2008 .

[32]  J. Marsden,et al.  Discrete mechanics and variational integrators , 2001, Acta Numerica.

[33]  Javier García de Jalón,et al.  Kinematic and Dynamic Simulation of Multibody Systems: The Real Time Challenge , 1994 .

[34]  A. Shabana,et al.  DEVELOPMENT OF SIMPLE MODELS FOR THE ELASTIC FORCES IN THE ABSOLUTE NODAL CO-ORDINATE FORMULATION , 2000 .

[35]  J. P. Dias,et al.  Sensitivity Analysis of Rigid-Flexible Multibody Systems , 1997 .