Search

The paper extends the formulation of a 2D geometrically exact beam element proposed in our previous paper [1] to curved elastic beams. This formulation is based on equilibrium equations in their integrated form, combined with the kinematic relations and sectional equations that link the internal forces to sectional deformation variables. The resulting first-order differential equations are approximated by the finite difference scheme and the boundary value problem is converted to an initial value problem using the shooting method. The paper develops the theoretical framework based on the Navier-Bernoulli hypothesis, with a possible extension to shear-flexible beams. Numerical procedures for the evaluation of equivalent nodal forces and of the element tangent stiffness are presented in detail. Unlike standard finite element formulations, the present approach can increase accuracy by refining the integration scheme on the element level while the number of global degrees of freedom is kept constant. The efficiency and accuracy of the developed scheme are documented by seven examples that cover circular and parabolic arches, a spiral-shaped beam, and a spring-like beam with a zig-zag centerline. The proposed formulation does not exhibit any locking. No excessive stiffness is observed for coarse computational grids and the distribution of internal forces is not polluted by any oscillations. It is also shown that a cross effect in the relations between internal forces and deformation variables arises, i.e., the bending moment affects axial stretching and the normal force affects the curvature. This coupling is theoretically explained in the appendix.

The paper presents a two-dimensional geometrically nonlinear formulation of a beam element that can accommodate arbitrarily large rotations of cross sections. The formulation is based on the integrated form of equilibrium equations, which are combined with the kinematic equations and generalized material equations, leading to a set of three first-order differential equations. These equations are then discretized by finite differences and the boundary value problem is converted into an initial value problem using a technique inspired by the shooting method. Accuracy of the numerical approximation is conveniently increased by refining the integration scheme on the element level while the number of global degrees of freedom is kept constant, which leads to high computational efficiency. The element has been implemented into an open-source finite element code. Numerical examples show a favorable comparison with standard beam elements formulated in the finite-strain framework and with analytical solutions.