A finite element/quaternion/asymptotic numerical method for the 3D simulation of flexible cables

A flexible cable is modeled by a geometrically exact beam model with 3D rotations described using quaternion parameters. The boundary value problem is then discretized by the finite element method. The use of an asymptotic numerical method to solve the problem, quadratic equations, is well suited to the quaternion parametrization. This combination of methods leads to a fast, robust and accurate algorithm very well-adapted for the simulation of the assembly process of cables. This is proved by running many examples involving complicated solutions