A Lie group variational integration approach to the full discretization of a constrained geometrically exact Cosserat beam model

In this paper, we will consider a geometrically exact Cosserat beam model taking into account the industrial challenges. The beam is represented by a framed curve, which we parametrize in the configuration space $\mathbb{S}^{3}\ltimes \mathbb{R}^{3}$ S 3 ⋉ R 3 with semi-direct product Lie group structure, where $\mathbb{S}^{3}$ S 3 is the set of unit quaternions. Velocities and angular velocities with respect to the body-fixed frame are given as the velocity vector of the configuration. We introduce internal constraints, where the rigid cross sections have to remain perpendicular to the center line to reduce the full Cosserat beam model to a Kirchhoff beam model. We derive the equations of motion by Hamilton’s principle with an augmented Lagrangian. In order to fully discretize the beam model in space and time, we only consider piecewise interpolated configurations in the variational principle. This leads, after approximating the action integral with second order, to the discrete equations of motion. Here, it is notable that we allow the Lagrange multipliers to be discontinuous in time in order to respect the derivatives of the constraint equations, also known as hidden constraints. In the last part, we will test our numerical scheme on two benchmark problems that show that there is no shear locking observable in the discretized beam model and that the errors are observed to decrease with second order with the spatial step size and the time step size.

Flexible multibody impact simulations based on the isogeometric analysis approach

Usually 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.

A review of flexible multibody dynamics for gradient-based design optimization

Design optimization of flexible multibody dynamics is critical to reducing weight and therefore increasing efficiency and lowering costs of mechanical systems. Simulation of flexible multibody systems, though, typically requires high computational effort which limits the usage of design optimization, especially when gradient-free methods are used and thousands of system evaluations are required. Efficient design optimization of flexible multibody dynamics is enabled by gradient-based optimization methods in concert with analytical sensitivity analysis. The present study summarizes different formulations of the equations of motion of flexible multibody dynamics. Design optimization techniques are introduced, and applications to flexible multibody dynamics are categorized. Efficient sensitivity analysis is the centerpiece of gradient-based design optimization, and sensitivity methods are introduced. The increased implementation effort of analytical sensitivity analysis is rewarded with high computational efficiency. An exemplary solution strategy for system and sensitivity evaluations is shown with the analytical direct differentiation method. Extensive literature sources are shown related to recent research activities.

Dynamics of a flexible body: a two-field formulation

A two-field formulation of the nonlinear dynamics of an elastic body is presented in which positions/orientations and the resulting velocity field are treated as independent. Combining a nonclassical description of elastic velocity that includes the convection velocity due to elastic deformation with floating reference axes minimizing the relative kinetic energy due to elastic deformation provides a fully uncoupled expression of kinetic energy. A transformation inspired by the classical Legendre transformation concept is introduced to develop the motion equations in canonical form. Finite element discretization is achieved using the same shape function sets for elastic displacements and velocities. Specific attention is brought to the discretization of the gyroscopic forces induced by elastic deformation. A model reduction strategy to construct superelement models suitable for flexible multibody dynamics applications is proposed, which fulfills the essential condition of orthogonality between a rigid body and elastic motions. The problem of expressing kinematic connections at superelement boundaries is briefly addressed. Two academic examples have been developed to illustrate some of the concepts presented.

On design sensitivities in the structural analysis and optimization of flexible multibody systems

The structural analysis and optimization of flexible multibody systems become more and more popular due to the ability to efficiently compute gradients using sophisticated approaches such as the adjoint variable method and the adoption of powerful methods from static structural optimization. To drive the improvement of the optimization process, this work addresses the computation of design sensitivities for multibody systems with arbitrarily parameterized rigid and flexible bodies that are modeled using the floating frame of reference formulation. It is shown that it is useful to augment the body describing standard input data files by their design derivatives. In this way, a clear separation can be achieved between the body modeling and parameterization and the system simulation and analysis.

A mortar formulation for frictionless line-to-line beam contact

This paper describes the quasi-static formulation of frictionless line contact between flexible beams by employing the mortar finite element approach. Contact constraints are enforced in a weak sense along the contact region using Lagrange multipliers. A simple projection appropriate for thin beams with circular cross-sections is proposed for the computation of contact regions. It is combined with the geometrically exact beam formalism on the Lie group $SE(3)$ S E ( 3 ) . Interestingly, this framework leads to a constraint gradient and a tangent stiffness invariant under rigid body transformations. The formulation is tested in some numerical examples.