
tact force actively since the force/wrench was implicitly
manipulated. Previously, explicitly tracking desired contact
forces in quadrupedal systems was implemented in [9] but
required planning on both desired position trajectories and
desired contact forces at the same time if imposing higher
priority on force control than motion control, because desired
accelerations affected realizable contact forces. The tracking
of planned force profiles may not be suitable for tasks
requiring fast reaction (e.g., tracking desired interaction force
command from users); thus, a reactive control scheme is
preferable.
To this end, to approach an exact force output, we propose
a novel interaction force control scheme for underactuated
quadrupedal systems in the sense of reaction control. It
allows the system (e.g., using a foot) to apply a force as
precisely as possible to the environment, without requiring
force planning. Our scheme uses two QP designs that resolve
the underactuation problem by splitting it into two orthogonal
spaces: motion and constraint space, and then optimizing the
cost in each space in a hierarchical order. Physical limitations
(i.e., unilateral contact, torque limits, and friction cones)
are also considered as inequality constraints in the design.
To decouple force and motion control and accomplish the
exact constraint force control as much as possible, we apply
two selection matrices, in which one of them distributes
the desired force control task to the designated joints, and
another one selects the rest joints for the underlying motion
task. For instance, as shown in Fig. 1, a quadrupedal robot
uses its front right leg for force control, while the other three
legs support its base and conduct the motion task of the base.
These two selection matrices select the front right leg and the
other three legs for the force and motion control, respectively.
In summary, our control scheme is a reactive control (e.g.,
can respond to user-defined force inputs fast) that does not
rely on any motion planning techniques and force sensors
(FS), and induces minimal base motion.
The contributions of this paper are:
1) Presenting a novel interaction force control scheme for
underactuated quadrupedal systems that does not require
motion planning and FS.
2) For resolving the underactuation problem, we propose
a hierarchical QP structure that minimizes the cost
function for the motion and force control in motion and
constraint space, respectively. To reduce the coupling
effect between motion and force control, two selection
matrices are deployed, allowing us to decouple force
and motion to the greatest extent and thus approach the
exact force control.
II. BACKGROUND
A. Projected Inverse Dynamics with External Disturbance
The dynamics equation of a quadruped robot can be
expressed as:
M¨
q+h=Bτ +JT
cλ+JT
xFx,(1)
where q=qT
b,qT
jTis the generalized coordinate vec-
tor including unactuated floating base coordinates qb∈
SE(3) and actuated joint configuration qj∈Rnj,M∈
R(6+nj)×(6+nj)denotes the inertia matrix, h∈R(6+nj)is
the non-linear effect consisting of Coriolis, centrifugal and
gravitational forces, B=0nj×6,Inj×njTis the selection
matrix1of actuated joints, τ∈Rnjis the actuated joint
torques, Jc∈Rk×(6+nj)represents the constraint Jacobian
with k= 3nc(ncdenotes the number of legs in contact
with the environment and the contact is assumed as a point
contact), λ∈Rkdenotes the generalized constraint force
vector used to control unactuated qb,Jx∈Rne×(6+nj)is
the Jacobian at x, and Fx∈Rneis the external disturbances
due to the interaction from human or environments.
During locomotion, the support feet should not slip, i.e.,
the constraint Jc˙
q=0must be satisfied. This constraint
indicates that any admissible ˙
qlies in the constraint null
space N(Jc). According to [10], the dynamics equation (1)
can be projected into two subspaces by using the orthogonal
projection matrix Pand I−Prespectively as:
P M ¨
q+P h =P Bτ
|{z }
τm
+P JT
xFx,(2)
(I−P)(M¨
q+h)=(I−P)Bτ
| {z }
τc
+JT
cλ+ (I−P)JT
xFx,
(3)
where P=I−J+
cJcimplies that the orthogonal null
space projection matrix is computed from the Moore-Penrose
pseudoinverse of Jcand projects vectors into the null space
of the constraint, such that P=P2=PT,P JT
c=0, and
P˙
q=˙
qfor all ˙
q∈ N (Jc). Then I−Prepresents the
complementary projection into N⊥(Jc).
To solve the equation of λ,¨
qmust be computed first. How-
ever, ¨
qcannot be obtained directly through pre-multiplying
the inverse of P M in (2), as P M cannot be inverted
attributed to the rank deficiency of P. With the help of
additional equations (I−P)˙
q=0and its derivative
(I−P)¨
q=˙
P˙
qthat are derived from P˙
q=˙
q,¨
qcan
be solved by substituting the latter equation into (2) as:
¨
q=M−1
c(τm−P h +˙
P˙
q+P JT
xFx),(4)
where Mc=P M +I−P. Equations (4) shows that
only the motion torques τmcontributes to the motion of
the system, therefore N(Jc)is called the motion space as
described in [6]. Similarly, because in (3) the ¨
qis fixed
and constraint forces λare free to choose depending on the
constraint torques τc,N⊥(Jc)is named the constraint space.
Eventually, the constraint forces are obtained by inserting ¨
q
into (3) and yields:
λ= (JT
c)+h(I−P)[ ¯
M(τm−P h +˙
P˙
q) + h]−τc
+ (I−P)( ¯
M P −I)JT
xFxi,
(5)
with ¯
M=M M −1
c, and is the constraint inertia matrix and
is always invertible [10]. More details can be found in [6].
1Unlike the selection matrix in previous papers, the selection matrix here
is not square.