Dynamics Study and Sensitivity Analysis of Flexible Multibody Systems With Interval Parameters

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.

Author(s):  
Alfonso Callejo ◽  
Valentin Sonneville ◽  
Olivier A. Bauchau

The gradient-based design optimization of mechanical systems requires robust and efficient sensitivity analysis tools. The adjoint method is regarded as the most efficient semi-analytical method to evaluate sensitivity derivatives for problems involving numerous design parameters and relatively few objective functions. This paper presents a discrete version of the adjoint method based on the generalized-alpha time integration scheme, which is applied to the dynamic simulation of flexible multibody systems. Rather than using an ad hoc backward integration solver, the proposed approach leads to a straightforward algebraic procedure that provides design sensitivities evaluated to machine accuracy. The approach is based on an intrinsic representation of motion that does not require a global parameterization of rotation. Design parameters associated with rigid bodies, kinematic joints, and beam sectional properties are considered. Rigid and flexible mechanical systems are investigated to validate the proposed approach and demonstrate its accuracy, efficiency, and robustness.


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):  
Kishor D. Bhalerao ◽  
Mohammad Poursina ◽  
Kurt S. Anderson

This paper presents a recursive direct differentiation method for sensitivity analysis of flexible multibody systems. Large rotations and translations in the system are modeled as rigid body degrees of freedom while the deformation field within each body is approximated by superposition of modal shape functions. The equations of motion for the flexible members are differentiated at body level and the sensitivity information is generated via a recursive divide and conquer scheme. The number of differentiations required in this method is minimal. The method works concurrently with the forward dynamics simulation of the system and requires minimum data storage. The use of divide and conquer framework makes the method linear and logarithmic in complexity for serial and parallel implementation, respectively, and ideally suited for general topologies. The method is applied to a flexible two arm robotic manipulator to calculate sensitivity information and the results are compared with the finite difference approach.


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.


Sign in / Sign up

Export Citation Format

Share Document