“…The equivalent mass M rj in the j-th Cartesian direction (j = x, y, z), depending on the robot configuration q, is given by [39] M rj = 1 e T j [J(q)B −1 (q)J T (q)] e j where J(q) and B(q) are the geometric Jacobian and the generalized mass matrix of the robot, respectively, q is the joint position vector, and e j is a unitary vector defined as follows Similarly, the equivalent stiffness K rj in the j-th Cartesian direction (j = x, y, z), depending on the robot configuration q, is given by [40] K rj = 1 e T j C ee e j where the end-effector compliance C ee is composed by two different contributions, one related to lumped joint compliance and the other one due to distributed arm stiffness, i.e.,…”