Advanced beam theory for multibody dynamics

In flexible multibody systems, many components are often approximated as beams or shells. More often that not, classical beam theories, such as Euler-Bernoulli beam theory, form the basis of the analytical development for beam dynamics. The advantage of this approach is that it leads to a very simple kinematic representation of the problem: the beam's section is assumed to remain plane and its displacement field is fully defined by three displacement and three rotation components. While such approach is capable of capturing the kinetic energy of the system accurately, it cannot represent the strain energy adequately. For instance, it is well known from Saint-Venant's theory for torsion that the cross-section will warp under torque, leading to a threedimensional deformation state that generates a complex stress state. To overcome this problem, sectional stiffnesses are computed based on sophisticated mechanics of material theories that evaluate the complete state of deformation. These sectional stiffnesses are then used within the framework of an Euler-Bernoulli beam theory based on far simpler kinematic assumptions. While this approach works well for simple cross-sections made of homogeneous material, very inaccurate predictions result for realistic sections, specially for thin-walled beams, or beams made of anisotropic materials. This paper presents a different approach to the problem. Based on a finite element discretization of the cross-section, an exact solution of the theory of threedimensional elasticity is developed. The only approximation is that inherent to the finite element discretization. The proposed approach is based on the Hamiltonian formalism and leads to an expansion of the solution in terms of extremity and central solutions, as expected from Saint-Venant's principle.