Reduced Isogeometric Analysis Models for Impact Simulations

2021 ◽  
Author(s):  
Tobias Rückwald ◽  
Alexander Held ◽  
Robert Seifried

Abstract Detailed impact simulations in flexible multibody systems are usually based on isoparametric finite element models. For modeling the dynamics of an impact, a precise representation of the geometry is essential. However, isoparametric finite elements involve the discretization of the geometry. This work tests the isogeometric analysis (IGA) as an alternative approach in flexible multibody systems. The IGA enables the exact representation of the geometry by using non-uniform rational basis splines (NURBS) as element shape functions. In the context of an efficient impact simulation a model reduction and a possible inclusion of the floating frame of reference formulation is beneficial. The degrees of freedom of the flexible bodies are reduced using component mode synthesis to save computation time in the multibody simulation. For the precise description of deformations and stresses in the contact area as well as elastodynamic effects, a large number of global shape functions is required. As testing examples, the impact of two elastic spheres and a multibody multicontact problem including wave propagation in a long elastic rod are simulated and compared to reference solutions.

Author(s):  
Tobias Rückwald ◽  
Alexander Held ◽  
Robert Seifried

AbstractUsually detailed impact simulations are based on isoparametric finite element models. For the inclusion in multibody dynamics simulation, e.g., in the framework of the floating frame of reference, a previous model reduction is necessary. A precise representation of the geometry is essential for modeling the dynamics of the impact. However, isoparametric finite elements involve the discretization of the geometry. This work tests isogeometric analysis (IGA) models as an alternative approach in the context of impact simulations in flexible multibody systems. Therefore, the adaption of the flexible multibody system procedure to include IGA models is detailed. The use of nonuniform rational basis splines (NURBS) allows the exact representation of the geometry. The degrees of freedom of the flexible body are afterwards reduced to save computation time in the multibody simulation. To capture precise deformations and stresses in the area of contact as well as elastodynamic effects, a large number of global shape functions is required. As test examples, the impact of an elastic sphere on a rigid surface and the impact of a long elastic rod are simulated and compared to reference solutions.


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.


2015 ◽  
Vol 10 (2) ◽  
Author(s):  
Wenhao Guo ◽  
Tianshu Wang

By means of a recursive formulation method, a generalized impulse–momentum-balance method, and a constraint violation elimination (CVE) method, we propose a new global simulation method for flexible multibody systems with kinematic structure changes. The constraint equations of a pair of adjacent bodies, considering body flexibility in Cartesian space, are derived for a recursive formulation. Constraint equations in configuration space, which are obtained from the constraints presented in this paper via recursive formulation, are very useful for modeling different kinematic structures and impacting governing equations. The novelty is that the impact governing equations, which calculate the jumps of generalized velocities, are modified by taking velocity-level CVE into consideration. Numerical examples are given to validate the presented method. Simulation results show that the new method can effectively suppress constraint drifts at the velocity level and stabilize constraint violations at the position level.


Author(s):  
Lars Ku¨bler ◽  
Peter Eberhard ◽  
Johannes Geisler

In this paper a formulation for flexible Multibody Systems (MBS) is proposed where flexible bodies are described using absolute coordinates for isoparametric brick elements. The use of absolute coordinates allows for large deformations and provides an accurate description of rigid body motion and inertia in the case of large rotations without additional considerations. Further, constant mass matrices are obtained for isoparametric elements. Brick elements are important, e. g. if general solid bodies with low stiffness, i. e. not negligible large deformations, are part of the MBS and cannot be modeled using beam, plate, or shell elements. Since only nodal translational degrees of freedom are used for brick elements additional questions arise. For example, imposing joint constraints for relative rotations between two bodies requires a nodal reference frame at connection points. An approach is proposed to define such a reference system utilizing displacement information of three finite element nodes.


Author(s):  
Thomas Gorius ◽  
Robert Seifried ◽  
Peter Eberhard

In many cases, the design of a tracking controller can be significantly simplified by the use of a 2-degrees of freedom (DOF) control structure, including a feedforward control (i.e., the inversion of the nominal system dynamics). Unfortunately, the computation of this feedforward control is not easy if the system is nonminimum-phase. Important examples of such systems are flexible multibody systems, such as lightweight manipulators. There are several approaches to the numerical computation of the exact inversion of a flexible multibody system. In this paper, the singularly perturbed form of such mechanical systems is used to give a semianalytic solution to the tracking control design. The control makes the end-effector to even though not exactly, but approximately track a certain trajectory. Thereby, the control signal is computed as a series expansion in terms of an overall flexibility of the bodies of the multibody system. Due to the use of symbolic computations, the main calculations are independent of given parameters (e.g., the desired trajectories), such that the feedforward control can be calculated online. The effectiveness of this approach is shown by the simulation of a two-link flexible manipulator.


2013 ◽  
Vol 30 (1) ◽  
pp. 13-35 ◽  
Author(s):  
Maria Augusta Neto ◽  
Jorge A. C. Ambrósio ◽  
Luis M. Roseiro ◽  
A. Amaro ◽  
C. M. A. Vasques

1999 ◽  
Vol 122 (4) ◽  
pp. 498-507 ◽  
Author(s):  
Marcello Campanelli ◽  
Marcello Berzeri ◽  
Ahmed A. Shabana

Many flexible multibody applications are characterized by high inertia forces and motion discontinuities. Because of these characteristics, problems can be encountered when large displacement finite element formulations are used in the simulation of flexible multibody systems. In this investigation, the performance of two different large displacement finite element formulations in the analysis of flexible multibody systems is investigated. These are the incremental corotational procedure proposed in an earlier article (Rankin, C. C., and Brogan, F. A., 1986, ASME J. Pressure Vessel Technol., 108, pp. 165–174) and the non-incremental absolute nodal coordinate formulation recently proposed (Shabana, A. A., 1998, Dynamics of Multibody Systems, 2nd ed., Cambridge University Press, Cambridge). It is demonstrated in this investigation that the limitation resulting from the use of the infinitesmal nodal rotations in the incremental corotational procedure can lead to simulation problems even when simple flexible multibody applications are considered. The absolute nodal coordinate formulation, on the other hand, does not employ infinitesimal or finite rotation coordinates and leads to a constant mass matrix. Despite the fact that the absolute nodal coordinate formulation leads to a non-linear expression for the elastic forces, the results presented in this study, surprisingly, demonstrate that such a formulation is efficient in static problems as compared to the incremental corotational procedure. The excellent performance of the absolute nodal coordinate formulation in static and dynamic problems can be attributed to the fact that such a formulation does not employ rotations and leads to exact representation of the rigid body motion of the finite element. [S1050-0472(00)00604-8]


Sign in / Sign up

Export Citation Format

Share Document