This work is focused on the development of numerical methods for the design and control of robots, with particular emphasis on joint elasticity. First, a general methodology is presented that is able to solve the problem of computing the inverse dynamics of a serial robot manipulator with an arbitrarily large number of elastic joints in a recursive numerical way. The solution algorithm is a generalized version of the standard Newton-Euler approach. The algorithm is presented with numerous extensions and variants, including the extension to variable-stiffness technologies and control applications. Then, an optimization framework is introduced for the design and analysis of biped walkers characterized by elastic joints, with comparative results demonstrating the scope of application of joint compliance in bipedal walking
