A Novel Approach to Real-Time Flexible Multibody Simulation: Sub-System Global Modal Parameterization

Author(s):  
Frank Naets ◽  
Gert H. K. Heirman ◽  
Wim Desmet

This paper introduces a novel model reduction technique, namely Sub-System Global Modal Parameterization (SS-GMP), for real-time simulation of flexible multibody systems. In the past, other system-level model reduction techniques have been proposed for this purpose, but these were limited in applicability due to the large storage requirements for systems with many rigid degrees-of-freedom (DOFs). However, in the SS-GMP approach, the motion of a mechanism is split up into a global motion and a relative motion of the (sub-)system. The relative motion is then reduced according to the Global Modal Parameterization, which is a model reduction procedure suitable for closed chain flexible multibody systems. In combination with suitable explicit solvers, the SS-GMP approach enables (hard) real-time simulations due to the strong reduction in the number of DOFs and the conversion of a system of differential-algebraic equations into a system of ordinary differential equations. The proposed approach is validated numerically with a quarter-car model. This fully flexible mechanism is simulated faster than real-time on a regular PC with the SS-GMP approach while providing accurate results.

Author(s):  
Gert H. K. Heirman ◽  
Wim Desmet

The presence of both differential and algebraic equations in the model equations, as well as the number of degrees of freedom needed to accurately represent flexibility, prohibit fast simulation of flexible multibody systems (e.g. real-time). In this research, Global Modal Parametrization, a model reduction technique for flexible multibody systems is further developed to speed up simulation of flexible multibody systems. The reduction of the model is achieved by projection on a curvilinear subspace instead of a fixed vector space, requiring significantly less degrees of freedom to represent the system dynamics with the same level of accuracy. The complexity of simulation of the reduced model equations is estimated. In a numerical experiment, simulation results for the original model equations are compared with simulation results for the model equations obtained after model reduction, showing a good match. The dominant sources of error of the proposed methodology are illustrated and explained.


Author(s):  
Robert Seifried ◽  
Markus Burkhardt

This paper presents inversion based feedforward control design for flexible multibody systems with kinematic loops and end-effector contact. The inverse model provides for a given desired output trajectories, e.g. end-effector point and contact force, the required control inputs for exact output reproduction. A very appealing and efficient model inversion approach for such multibody systems is the use of so-called servo-constraints. These can be seen as an extension of classical mechanical constraints and yield a set of differential-algebraic equations. This allows an efficient numerical solution without burdensome symbolic manipulations. In addition, the use of servo-constraints allows the straight-forward treatment of flexible multibody systems with various topologies. The arising set of differential-algebraic equations describes the inverse model. The inverse model might be purely algebraic or include a dynamical part, which is called internal dynamics in nonlinear control theory. For its numerical solution it is advisable to transform the set of differential-algebraic equations to its underlying set of ordinary differential equations. The solution method for this internal dynamics depends then on its stability. For systems with unstable internal dynamics, as considered in this paper, a solution can be computed from a boundary-value problem. The efficiency of this approach is demonstrated for a flexible multibody system with a kinematic loop and a closed end-effector contact.


Author(s):  
Wang Zhe ◽  
Qiang Tian ◽  
Hiayan Hu

The dynamics of flexible multibody systems with interval parameters is studied based on a non-intrusive computation methodology. The Absolute Nodal Coordinate Formulation (ANCF) is used to model the rigid-flexible multibody system, including the finite elements of the ANCF and the ANCF Reference Nodes (ANCF-RNs). The Chebyshev sampling methods including Chebyshev tensor product (CTP) sampling method and Chebyshev collocation method (CCM), are utilized to generate the Chebyshev surrogate model for Interval Differential Algebraic Equations (IDAEs). For purpose of preventing the interval explosion problem and maintaining computation efficiency, the interval bounds of the IDAEs are determined by scanning the deduced Chebyshev surrogate model. To further improve the computation efficiency, OpenMP directives are also used to parallelize the solving process of the Differential Algebraic Equations (DAEs) by fixing the uncertain interval parameter at the given sampling points. The sensitivity analysis of flexible multibody systems with interval parameters is initially performed by using the direct differentiation method. The direct differentiation method differentiates the dynamic equations with respect to the design variable, which yields the system sensitivity equations governed by DAEs. The generalized alpha method is introduced to integrate the sensitivity DAEs. The sensitivity equations of flexible multibody systems with interval parameters are also described by the IDAEs. Based on the continuum mechanics, the computational efficient analytical formulations for the derivative items of the system sensitivity equations are deduced. Three examples are studied to validate the proposed methodology, including the complicated spatial rigid-flexible multibody systems with a large number of uncertain interval parameters, the flexible system with uncertain interval clearance size joint, and the first order sensitivity analysis of flexible multibody systems with interval parameters. Firstly, the dynamics analysis of a six-arm space robot with six interval parameters is performed. For this case study, the interval dynamics cannot be obtained by directly scanning the IDAEs because extremely huge sets of DAEs with deterministic samples have to be solved. The estimated total computational time for solving the scanned IDAEs will be 1850 days! However, the computational time for solving the scanned Chebyshev surrogate model is 9796.97 seconds. It shows the effectiveness of the proposed computation methodology. Then, the nonlinear dynamics of a planar slider-crank mechanism with uncertain interval clearance size joint is studied in this work. The kinetics model of the revolute clearance joints is formulated under the ANCF-RN framework. Moreover, the influence of the LuGre and the modified Coulomb’s friction force models on the system’s dynamic response is investigated. By analyzing the bounds of dynamic response, the bifurcation diagrams are observed. It must be highlighted that with increasing the size of clearance, it does not automatically lead to unstable behaviors. Finally, the first order sensitivity analysis of flexible multibody systems with interval parameters is also studied in this work. The third one of a flexible mechanism with interval parameters is used to perform the sensitivity analysis.


Sign in / Sign up

Export Citation Format

Share Document