A New Kirchhoff-Love Beam Element and its Application to Polymer Mechanics

The novel finite element formulations fall into the category of geometrically exact Kirchhoff-Love beams. A prominent characteristic of this category is that the absence of shear deformation is strongly enforced by removing two degrees of freedom. Further, the corresponding beam theories exhibit not only translational but also rotational degrees of freedom and their configurations thus form a non-additive and non-commutative space. Sophisticated interpolation schemes are required that need to be tested not only for locking, spatial convergence behavior, and energy conservation, but also for observer invariance and path-independence. For the three novel beam element formulations all these properties are analytically and numerically studied and confirmed, if applicable. Two different rotation parameterization strategies are employed based on the well-known geodesic interpolation used in many Simo-Reissner beams and the lesser known split into the so-called extit{smallest rotation} and a torsional part. Application of the former parameterization results in a mixed finite element formulation intrinsically free of locking phenomena. Additionally, the first geometrically exact Kirchhoff-Love beam element is presented, which strongly enforces inextensibility by removing another degree of freedom. Furthermore, the numerical efficiency of the new beam formulations is compared to other beam elements that allow for or suppress shear deformation. When modeling very slender beams, the new elements offer distinct numerical advantages.