We present an efficient numerical scheme based on asymptotic numerical method for continuation of spatial equilibria of special Cosserat rods. Using quaternions to represent rotation, the equations of static equilibria of special Cosserat rods are posed as a system of thirteen first order ordinary differential equations having cubic nonlinearity. The derivatives in these equations are further discretized to yield a system of cubic polynomial equations.
Recent comments