Hamilton’s equations can be used to define the dynamics of a tree configured flexible multibody system. Their states are the generalized coordinates and momenta . Numerical solution of these equations requires the time derivatives of the states be defined. Hamilton’s equations have the benefit that the time derivative of the system momenta are easy to compute. However, the generalized velocities need be solved from the system momenta as defined by to support the computation of and the propagation of . Because of the size of , the determination of by linear equation solution schemes requires order arithmetic operations, where is the number of bodies and is the number of mode shape functions used to model the body deformations. It has been shown that can be solved recursively from the momentum equations for rigid multibody systems (Naudet, Lefeber, and Terze, 2003, “Forward Dynamics of Open-Loop Multibody Mechanisms Using An Efficient Recursive Algorithm Based On Canonical Momenta,” Multibody Syst. Dyn., 10, pp. 45–59). This paper extends that result to flexible multibody systems. The overall arithmetic operations to solve for in this case is proportional to if the effort to solve for the flexible coordinate rates for each body is weighted the same as that for the joint rate. However, each time the flexible coordinates rate of a body is solved an order operations is incurred. Thus, the total computational effort for flexible multibody systems includes an additional order operations.