The absolute nodal coordinate formulation (ANCF) is a nonlinear finite element approach proposed for the large deformation dynamics analysis of beam- and plate/shell-type structures. In the ANCF approach, elastic forces can be defined using three-dimensional elasticity-based continuum mechanics. This approach is often straightforward, and it makes it possible to use advanced material models in the ANCF framework. However, it has been pointed out in several studies that continuum ANCF-based elements with a full three-dimensional elasticity description can suffer from locking phenomena. In this study, a comparison between various combinations of locking alleviation techniques and their applicability to different ANCF beam variants is studied using numerical examples. Furthermore, the enhanced deformation gradient (EDG) technique, which has been proposed recently in finite element literature, is demonstrated for high-order ANCF beam elements. Based on the numerical tests, none of the currently available techniques are suitable for all types of ANCF elements. The paper also shows that the efficiency and accuracy of the techniques are case-dependent. For the ANCF beam element involving higher-order terms with respect to trapezoidal mode, however, the EDG-based techniques are preferable to reduce locking phenomena.

Experimental results have revealed the sophisticated Achilles tendon (AT) structure, including its material properties and complex geometry. The latter incorporates a twisted design and composite construction consisting of three subtendons. Each of them has a nonstandard cross-section. All these factors make the AT deformation analysis computationally demanding. Generally, 3D finite solid elements are used to develop models for AT because they can discretize almost any shape, providing reliable results. However, they also require dense discretization in all three dimensions, leading to a high computational cost. One way to reduce degrees of freedom is the utilization of finite beam elements, requiring only line discretization over the length of subtendons. However, using the material models known from continuum mechanics is challenging because these elements do not usually have 3D elasticity in their descriptions. Furthermore, the contact is defined at the beam axis instead of using a more general surface-to-surface formulation. This work studies the continuum beam elements based on the absolute nodal coordinate formulation (ANCF) for AT modeling. ANCF beam elements require discretization only in one direction, making the model less computationally expensive. Recent work demonstrates that these elements can describe various cross-sections and materials models, thus allowing the approximation of AT complexity. In this study, the tendon model is reproduced by the ANCF continuum beam elements using the isotropic incompressible model to present material features.