A nonlinear finite element framework for flexible multibody dynamics

We present a uniform treatment of rigid body dynamics and nonlinear structural dynamics. The advocated approach is based on a rotationless formulation of rigid bodies, nonlinear beams and shells. In this connection, the specific kinematic assumptions are taken into account by the explicit incorporation of holonomic constraints. This approach facilitates the straightforward extension to flexible multibody dynamics by including additional constraints due to the interconnection of rigid and flexible bodies. We further address the design of energy-momentum schemes for the stable numerical integration of the underlying finite-dimensional Hamiltonian systems. To demonstrate the superior numerical performance of the proposed methodology, the numerical examples deals with a multibody system containing both rigid and flexible bodies undergoing large deformations. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Natural co-ordinates for control applications

Natural coordinates have emerged to be well-suited for both rigid and flexible multibody dynamics. Especially the combination of structural elements and energy-momentum consistent time stepping schemes leads to superior numerical stability as well as an automatable assembly, resulting in both excellent run-time behaviour as well as moderate modelling effort (see [1]). Incorporation of modern methods for finite-element simulations, such as mortar methods for contact or domain decomposition both for structural elements as well as continuum elements is straightforward ([2]). Augmentation techniques allow a systematic integration of both mechanical and non-mechanical quantities for simulation (see [3] and [4]), which makes this approach suitable especially for emulation and simulation of mechatronic systems. We will present an approach for evaluating forward control strategies with flexible multibody systems. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A flexible approach to the simulation of hybrid multibody systems

The present works deals with the incorporation of both flexible beam and shell structures into the realm of flexible multibody dynamics. Geometrically exact beam formulations based on classical Simo-Reissner kinematics are suitable for modelling beam-type flexible components in the context of finite-deformation multibody dynamics. So geometrically exact shell formulations are based on Reissner-Mindlin kinematics. In [2], a flexible framework for dealing with flexible structural elements in a multibody context is described. A specific isoparametric finite element discretization of a shell formulation leads to semi-discrete equations of motions assuming the form of differential-algebraic equations (DAEs). A compatible isoparametric formulation of beams has already been developed in [1]. The uniform DAE framework makes possible the incorporation of alternative finite element formulations. In addition to that, various time-stepping schemes such as energy-momentum methods or variational integrators can be applied. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Real-Time Simulation of Elastic Multibody Systems with Application in Vehicle Dynamics

Real-time simulation models are widely used for vehicle development, usually built up as rigid multibody systems. However, since lightweight structures are commonly used, body deformation is no longer negligible and rigid multibody simulations may be inaccurate. This work presents a real-time capable full vehicle model with a flexible car body, derived from a finite element model, whose performance has been improved by model order reduction. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Maxwell-Contact

In multibody systems, two alternative approaches exist to compute contact forces between bodies. The unilateral constraint contact and the unilateral regularized contact do not take into account the influence of the deformation due to other contacts on a body using rigid body dynamics. In this paper, a third alternative, i.e. the Maxwell-Contact, is derived coupling the deformations of different contacts on one body quasi-statically. An academic example validates the fundamental properties of the contact model and the application in a simulation of a pushbelt CVT shows improved results. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Multilayer beams: A geometrically exact formulation

Summary We review and extend our recent work on a new theory of multilayer structures, with particular emphasis on sandwich beams/1-D plates. Both the formulation of the equations of motion in the general dynamic case and the computational formulation of the resulting nonlinear equations of equilibrium in the static case based on a Galerkin projection are presented. Finite rotations of the layer cross sections are allowed, with shear deformation accounted for in each layer. There is no restriction on the layer thickness; the number of layers can vary between one and three. The deformed profile of a beam cross section is continuous, piecewise linear, with a motion in 2-D space identical to that of a planar multibody system that consists of three rigid links connected by hinges. With the dynamics of this multi (rigid/flexible) body being referred directly to an inertial frame, the equations of motion are derived via the balance of (1) the rate of kinetic energy and the power of resultant contact (internal) forces/couples, and (2) the power of assigned (external) forces/couples. The present formulation offers a general method for analyzing the dynamic response of flexible multilayer structures undergoing large deformation and large overall motion. With the layersnot required to have equal length, the formulation permits the analysis of an important class of multilayer structures with ply drop-off. For sandwich structures, an approximated theory with infinitesimal relative outer-layer rotations superimposed onto finite core-layer rotation is deduced from the general nonlinear equations in a consistent manner. The classical linear theory of sandwich beams/1-D plates is recovered upon a consistent linearization. Using finite element basis functions in the Galerkin projection, we provide extensive numerical examples to verify the theoretical formulation and to illustrate its versatility. Dedicated to the memory of Professor Juan Carlos Simo, whose early demise is a great loss for the applied and computational mechanics community This paper was solicited by the editors to be part of a volume dedicated to the memory of Juan Simo.     

Using Multibody Systems for the Investigation of Dynamic Aberrations in High Precision Optics

The coupling of multibody system dynamics and optical simulations by using ray tracing through moving lens systems will be summarized for both rigid and flexible lens systems. Furthermore, a method will be introduced for efficient simulations in the time domain. This method provides a direct integration of the optical simulation into the equations of motion of a multibody system. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Non-linearity of multibody dynamic equations with respect to Lagrange multipliers: application to railway dynamics

In this paper, the dynamics of multibody systems with closed kinematical chains of bodies is considered. The main focus is set on non-linearity of the multibody equations with respect to the Lagrange multipliers. When closed chains are considered, loop cutting procedure is a solution to express the constraint equations associated with the loops. Dynamic equations of the multibody tree-like structure are thus completed with the constraint forces via the Lagrange multipliers. In the considered case of railway vehicles, constraints arise from the contact between the rigid wheels and the rails. Corresponding contact forces applied to the wheels appears via the Lagrange multipliers λ and the tangent creep forces as well. Resulting differential-algebraic equations can be transformed into an ODE system and then time-integrated using the coordinate partitioning method [3], when the system is linear with respect to λ. This paper presents an algorithm allowing us to solve this system in case of nonlinearities with respect to λ, which is typical of wheel/rail contact force models. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Optimization of Nonlinear Mechanical Systems under Constraints with the Particle Swarm Method

In recent years, stochastic optimization methods have gained increasing attention in parameter optimization of mechanical systems. Most popular techniques are Evolutionary Computation and the Simulating Annealing algorithm, which are applied more frequently to mechanical problems due to the increasing computing resources available now. Since theses methods do not require any gradient information, they are well suited for non‐smooth or discontinuous optimization tasks occurring in nonlinear multibody systems. In addition to these techniques, Kennedy and Eberhart [5] introduced the Particle Swarm Optimization method (PSO) based on the simulation of bird flocking. In this work, the efficiency of an extended PSO algorithm has been compared with an Evolutionary Strategy (ES) [6] and an Adapted Simulated Annealing method (ASA) [4]. In order to solve optimization tasks with both equality and inequality constraints the PSO algorithm has been extended by the Augmented Lagrangian Multiplier Method [2]. The proposed method shows often superior results and is quite simple to implement. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Inverse dynamics of underactuated multibody systems

The present work deals with controlled mechanical systems subject to holonomic constraints. In particular, we focus on underactuated systems, defined as systems in which the number of degrees of freedom exceeds the number of inputs. The governing equations of motion can be written in the form of differential-algebraic equations (DAEs) with a mixed set of holonomic and control constraints. The rotationless formulation of multibody dynamics will be considered [1]. To this end, we apply a specific projection method to the DAEs in terms of redundant coordinates. A similar projection approach has been previously developed in the framework of generalized coordinates by Blajer & Kołodziejczyk [2]. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A numerical method for the servo constraint problem of underactuated mechanical systems

Servo constraints are used in inverse dynamics simulations of discrete mechanical systems, especially for trajectory tracking control problems [1], whose desired outputs are represented by state variables and treated as servo constraints [2]. Servo constraint problems can be classified into fully actuated and underactuated multibody systems, and the equations of motion take the form of differential algebraic equations (DAEs) including holonomic and servo constraints. For fully actuated systems, control inputs can be solved from the equations by model inversion, as the input distribution matrix is nonsingular and invertible. However, underactuated systems have more degrees of freedom than control inputs. The input distribution matrix is not invertible, and in contrast to passive constraints, the realization of servo constraints with the use of control forces can range from orthogonal to tangential [3]. Therefore, it is challenging for the determination of control inputs which force the underactuated system to realize the partly specified motion. For differentially flat underactuated systems, the differentiation index of DAEs may exceed three. Hence we need to apply specific index reduction techniques, such as the projection approach applied in [3], [4], and [6]. The present work applies index reduction by minimal extension [5] to differentially flat underactuated crane systems and shows that the index can be reduced from five to three and even to one. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Influence of model reduction techniques on the impact force calculation of two flexible bodies

One important issue for the simulation of flexible multibody systems is the reduction of the flexible body's degrees of freedom. For the reduction process finite element data and user inputs are necessary. The model reduction program for elastic multibody systems MOREMBS, which is developed at the ITM, has an easy-to-use interface and the data can be gained from the programs ABAQUS or ANSYS. In this work, the simulation of a fuel injection process is investigated with MOREMBS. We focus on the interaction between valve and armature. These two bodies impact in every injection circle. The impacting bodies are modeled as flexible and the contact force is calculated by a penalty approach. One essential part of this work is the investigation of the influence of different model reduction techniques on the impact force calculation of the flexible multibody system. The main reduction techniques modal reduction, Krylov-subspace based and Gramian matrix based techniques are compared. The results achieved with modal reduction, the state of the art reduction method, are not acceptable here. Krylov-subspace based techniques are especially well-suited for large sparse systems but are not error controlled. However, by choosing appropriate moment-matching properties the impact force calculation is nearly as good as with a full finite element model. The Gramian matrix based reduction techniques can be fully automated and are error controlled but require high computational effort. Hence, appropriate approximation schemes have to be used for them. With Gramian matrix based methods we can even further reduce the size of the subsystems compared to Krylov-subspace based methods and still have an impact force calculation nearly as good as with finite element results, but we gain a simulation speedup by the factor 4000. In addition, a parameter study of the parameters involved in the model reduction process is presented. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical integration of differential-algebraic equations with mixed holonomic and control constraints

The present work aims at the incorporation of control (or servo) constraints into finite–dimensional mechanical systems subject to holonomic constraints. In particular, we focus on underactuated systems, defined as systems in which the number of degrees of freedom exceeds the number of inputs. The corresponding equations of motion can be written in the form of differential–algebraic equations (DAEs) with a mixed set of holonomic and control constraints. Apart from closed–loop multibody systems, the present formulation accommodates the so–called rotationless formulation of multibody dynamics. To this end, we apply a specific projection method to the DAEs in terms of redundant coordinates. A similar projection approach has been previously developed in the framework of generalized coordinates by Blajer & Kołodziejczyk [1]. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)