This paper presents a systematic formulation for the kinematic and dynamic analysis of flexible multibody systems. The system equations of motion are derived in terms of relative and elastic coordinates using velocity transformation technique. The position transformation equations that relate the relative and elastic coordinates to the Cartesian coordinates for the two contiguous flexible bodies are derived. The velocity transformation matrix is derived systematically corresponding to the type of kinematic joints connecting the bodies and system path matrix. This matrix is employed to represent the equations of motion in relative coordinate space. Two examples are taken to test the method developed here.