A recursive rotational-coordinate-based formulation of a planar Euler–Bernoulli beam is developed, where large displacements, deformations, and rotations are considered. Different from the traditional rotational-coordinate-based formulations, relative rotational angles rather than absolute ones are used as generalized coordinates. The number of generalized coordinates is minimized, which is inherited from traditional rotational-coordinate-based formulations. A recursive scheme is used for adjacent elements of the beam that is considered as a chain-like structure. Integrals in the mass matrix and generalized force vector of each element of the beam can be analytically derived, and a numerical trick based on Taylor polynomial approximations is adopted to avoid numerical singularity. The current formulation entirely avoids the evaluation of integrals in each time step, which greatly improves the computational efficiency. Three widely-used examples are studied to illustrate the performance of the proposed method. Results indicate that the present formulation can achieve the same accuracy, as well as much higher efficiency, compared to some traditional formulations. In addition, calculation time of the current formulation almost linearly increases with the increasing number of elements of the beam, indicating that the computational complexity of the current formulation is $O\left ( N \right )$ .
Read full abstract