# Realization of a Biped Robot Lower Limb Walking without Double Support Phase on Uneven Terrain

Authors: Hai-yan Wang; Yi-bin Li
Journal: Journal of Control Science and Engineering
 (2013)

---

## Abstract

Zero moment point (ZMP) is widely used in dynamical walking control of the biped robot, but it is hard to obtain the ZMP exactly. The paper describes a simple walking control method without using ZMP information directly. Firstly, the paper introduced a biped robot lower-limb prototype which is driven by linear hydraulic servocylinder. Then the paper simplifies the walking control in the lateral plane with a simple walking pattern generation method named “dynamic equilibrium method,” which is fit for active and underactuated biped robots. In the following section the paper provides the balance control methods without using ZMP information directly. Finally, simulation experiments with MD.DAMS and experiments in physical prototype are given. The experimental results confirm the effectiveness of the proposed control methods.


---

## Body

## 1. Introduction

To investigate the stability of human walking, Vukobratovic proposed a concept called zero moment point (ZMP) [1]. Although the ZMP algorithm and its variants have been used to realize stable walking for many bipedal robots, they are somewhat complex and computationally intensive since the full dynamics of the robot is considered. Kajita et al. introduced the single mass linear inverted pendulum model (LIPM), instead of using the complicated model where multibody mass is considered [2]. With three assumed condition, the ZMP equation for an n-DOF (degree of freedom) biped robot was simplified. Kajita also brought a successful result in the real experiment, and the maximum walking speed of HRP-2 was 2.5 km/h [3]. Kim introduced convolution sum [4], in which the proposed online algorithm has a similarity with the preview control [2] that needs a finite future reference, but it is not required to solve the Riccati equation to get gains. Kim et al. successfully applied this algorithm to the passenger-carrying biped robot HUBO FX-1 [5]. The control strategies provided by Kajita et al. and Kim et al. have no difference in the lateral plane and the sagittal plane. Walking pattern generation with LIPM only fits for the biped robot whose mass of legs can be neglected. But there are some robots that do not meet that assumed condition. The paper simplifies the walking control in the lateral plane with a simple walking pattern generation method named “dynamic equilibrium method.” It is fit for not only active biped but also underactuated robots, and there is no limit to the mass distribution between the body and the legs.Biped robots with weight loading capability are different from the general life-size humanoid robots, such as HRP [4, 5], ASIMO [6, 7], KHR [8], and BHR [9, 10]. Posture control of the upper body was helpful to realize the stable motion of the life-size humanoid robots. But the biped robots with the weight loading capability only have lower limb and a trunk which could be equipped with scaffolds for the loading or seats for the passenger as the users’ needs. The WL-16 has 6-DOF parallel mechanisms and was the first biped walking robot capable of carrying a human in the world [11, 12]. The weight of this original robot is 56 kgf, including 8 kgf battery weight, and the height is approximately 1.2 m. With the same load ability, the robots with electric motor driving system powered by batteries are usually large in volume and heavy in weight. Toyota released their “i-foot” robot in 2004 [13]. The height of the “i-foot” is 2.36 m, with a weight of 200 kgf. The maximum weight capacity of passenger was 60 kgf. Kim et al. proposed HUBO FX-1 [5] whose height is 1.987 m, with a weight of 130 kgf with the exception of the external covers and the cockpit. The maximum weight capacity of passenger is 100 kgf.Compared with electric motor driving system, hydraulic driving system has several important advantages [14], such as higher power-to-mass ratio, larger load capability, greater endurance, faster dynamic response, and better load stiffness. The paper introduced a hydraulic actuated biped robot lower limb with the theoretical maximum load capacity of 200 kgf, which is limited to the output force of the hydraulic cylinder. (In ADAMS simulation environment, there is no limit to the weight capability, and then we set a box on the robot trunk which is 100 kgf in weight.) The maximum height of the basic machine is about 882.4 mm, and the weight of the physical prototype robot is about 40 kg, which is mainly used in the no-load test.


## 2. Overview of the Hydraulic Actuated Biped Robot Lower Limb

The hydraulic actuated biped robot lower limb consists of the original robot system and hydraulic driven system. Most humanoid robots has six freedoms in one leg, three in the hip, one in the knee, and two in the ankle [15]. According to the kinetic characteristics of the hydraulic actuated biped robot, the leg structure is simplified here. The biped robot lower limb only has four active freedoms in one leg (shown as Figure 1), two in the hip, one in knee, and one in ankle. The trunk is used to load weight or install scaffold, and there is no freedom on the trunk. Table 1 is the specifications of the hydraulic biped robot. For it is not easy to set the buffer mechanism in the ADAMS simulation model, so the rolling joints in the ankles are not locked. Each of them is installed a torsion spring with the stiffness coefficient of 106 N·mm/rad and a rotary damper with the damping coefficient of 106 N·mm·s/rad. The roll angle in the ankle is used to replace the deformation of the buffer mechanism fixed on the foot bottom.Table 1
Specifications of the hydraulic biped robot.
Weight (kg)
Trunk: 12.53A single leg: 13.61
Size (mm)
Height inZc: 882.36Width in Yc: 482Length in Xc: 440
DOF
Hip: 2Knee: 1Ankle: 1Figure 1
The physical prototype.Considering economy problems, all the joints in the biped lower limb are actuated by identical linear hydraulic servocylinder (shown as Figure2), which is composed of a hydraulic cylinder, a hydraulic servovalve, a linear displacement sensor and a force sensor. The maximum displacement of the linear hydraulic servocylinder is 60 mm.Figure 2
Linear hydraulic servocylinder.


## 3. Walking Pattern Generation

In the biped robotics research field, the gait trajectory, otherwise known as the walking pattern, generates the relative position trajectories of two feet with respect to the pelvis center [16]. The paper set the geometric center of the trunk upper surface as the origin point of the basic coordinate system; Xc points to the direction of forward motion, Zc is perpendicular to the trunk (shown as Figure 3), and then XcOcZc is the frontal plane and YcOcZc is the lateral plane.Figure 3
The diagram of the motion direction.


### 3.1.  “Mobile Equilibrium Method” in the Lateral Plane

The three-dimensional linear inverted pendulum mode (3D-LIPM) combined with zero moment point criterion are widely used in the biped robot control field. But the 3D-LIPM only fits for the biped robots with active ankle joints and whose legs are light compared with the upper part of the body. This paper proposes an easier walking pattern method which has no limits to the robot mass distribution. It also fits for the robot without active ankle joint in the lateral plane.We consider the robot’s center of mass (COM) as the research object. When feet move accelerated toward left relative to the trunk in the lateral plane, the inertia force in the robot’s COM points to right (shown as Figure4 left). If the acceleration is large enough and the ground can provide enough friction and support force, the force between the left foot and the ground becomes zero, and the trunk moves right. If the roll joint of the ankle in the lateral plane is active, in the ideal condition the trunk moves right without trunk roll motion (shown in Figure 4, middle); if the ankle does not have roll joints, the trunk’s roll angle is large when the trunk move right (shown as Figure 4 right); if the ankle has roll joints but they are passive, the trunk’s roll angle exists, but it is smaller.Figure 4
Force analysis in the lateral plane.The biped robot lower-limb introduced in this paper does not have roll joint in the ankle, that is to say the roll angle of the trunk must exist when the feet move accelerated to left/right relative to the trunk (suppose that the ground can provide enough friction force, this is the precondition). If the trunk’s roll angle changes periodically, it shows that the robot is stable in the lateral plane; otherwise, we should make some adjustments.When the feet move accelerated to left relative to the trunk with appropriate acceleration, the support force on the left foot decreases to zero, then the left foot lifts, and the trunk will not dump in the lateral plane. The theory is fit for the opposite direction. In order to keep the robot stable in the lateral plane, the following two conditions are needed.(1)
The initial acceleration inYc direction is adequate. If the acceleration is too small the swing foot cannot lift or the robot topples toward the swing leg; if the acceleration is too large, the robot topples toward the support leg.(2)
In the moving process, the resultant moment is zero which is created by the inertial force and the gravity relative to equilibrium pointO1 (the equilibrium is O1, when the robot COM moves to right; the equilibrium is O2, when the robot COM moves to left. Actually O1 and O2 are zero moment points).In the above process, the acceleration of the trunk keeps on changing with the displacement changes in the lateral plane, so the method is named as “mobile equilibrium method.”


### 3.2. Walking Pattern Generation inYc Direction

Two step periods are considered as one movement circle. Based on the “mobile equilibrium method”, we can obtain the walking pattern inYc direction.In the first step period, the feet move accelerated toward left relative to the trunk(1)mg(-Lf1-y)+my¨zc=0.
Then,
(2)y=c1et·g/zc+c2e-t·g/zc+Lf1,
where
(3)t∈[0,Tn],c1=-e-Tng/zc·Lf1eTng/zc-e-Tng/zc,c2=eTng/zc·Lf1eTng/zc-e-Tng/zc,t∈[Tn,Tc],c1=e-Tng/zc·Lf1eTcg/zc·e-Tng/zc-e-Tcg/zc·eTng/zc,c2=-eTng/zc·Lf1eTcg/zc·e-Tng/zc-e-Tcg/zc·eTng/zc.In the second step period, the feet move accelerated toward right relative to the trunk:(4)mg(Lf2-y)+my¨zc=0.
Then,
(5)y=c3et·g/zc+c4e-t·g/zc-Lf2.
While t∈[Tc,Tc+Tn],
(6)c3=e-(Tc+Tn)g/zc·Lf2e(Tc+Tn)g/zc·e-Tcg/zc-e-(Tc+Tn)g/zc·eTcg/zc,c4=-e(Tc+Tn)g/zc·Lf2e(Tc+Tn)g/zc·e-Tcg/zc-e-(Tc+Tn)g/zc·eTcg/zct∈[Tc+Tn,2Tc],c3=e-(Tc+Tn)g/zc·Lf2e(Tc+Tn)g/zc·e-2Tcg/zc-e-(Tc+Tn)g/zc·e2Tcg/zc,c4=-e(Tc+Tn)g/zc·Lf2e(Tc+Tn)g/zc·e-2Tcg/zc-e-(Tc+Tn)g/zc·e2Tcg/zc.m is the mass of the whole biped lower limb; g is the gravitational acceleration; y is the displacement of the robot COM relative to the ground, which is opposite to the displacement of the ankle joint relative to the body fixed coordinate system in Yc direction; y¨ is the acceleration of the robot COM relative to the trunk fixed coordinate system in Yc direction; zc is the absolute distance between the robot COM and the ground; Lf1 is the absolute distance between O1 and Oc in Yc direction; Lf2 is the absolute distance between O2 and Oc in Yc direction. If the robot is symmetrical, then Lf1=Lf2=Lf, actually if there is no sliding between the feet and the ground. Lf is the maximum displacement of the trunk COM in the lateral plane, and we named it as maximum offset in the lateral plane.


### 3.3. Walking Pattern Generation inXc Direction

The direction changes of the trunk’s acceleration will cause the loading sway back and forth. The larger the acceleration, the bigger the sway range is. In order to decrease this impact, the trunk is set to a constant velocity inXc direction (the COM adjustments change the velocity of the trunk slightly). That is to say, the ankle joint moves in an approximately constant velocity relative to the trunk fixed system, and then the trunk can keep a constant velocity relative to the ground. When the swing foot changes to the support foot, the swing foot should keep the same velocity as the trunk. With this constraint one of the easiest gait generation methods for swing leg is cubic polynomial fitting. One periodic motion contains two step periods.In the first step period, the trajectory of the ankle joint of the swing leg relative to the trunk fixed coordinate system inXc direction is shown as follows:
(7)xs=-8DsTc3t3+12Ds-DdTc2t2-2Ds-DdTct-Ds.The trajectory of the ankle joint of the support leg relative to the trunk fixed coordinate system inXc direction is shown as follows:
(8)x=-2DsTct+Ds,
where xs is the ankle displacement of the swing leg relative to the trunk fixed coordinate system in Xc direction; x is the ankle displacement of the support leg relative to the trunk fixed coordinate system in Xc direction; Ds is the step length of the ankle joint relative to the trunk fixed coordinate system in Xc direction; Dd is the increment of the step length between the two contiguous periods: if the step length increases, Dd is a positive, and if the step length decreases, Dd is a negative; Tc is the step period.The walking generation method in the second step period is the same as in the first step period, but the swing leg changes to the support leg and the support leg changes to the swing leg.


### 3.4. Walking Pattern Generation inZc Direction

In order to decrease the impact, when the swing foot hit the ground and changes to the support foot, the following two conditions are needed:(1)
the ideal dropping velocity inZc direction is zero;(2)
the raise velocity decreases to zero when the swing foot reaches the planned raise height.With the constraints, cubic polynomial fitting can be used here. In the first step period, the trajectory of the ankle in the swing leg relative to the trunk inZc direction is shown as follows:
(9)t∈[0,Tc2],Hz=-16HnTc3t3+12HnTc2t2,t∈[Tc2,Tc],Hz=16HnTc3t3-36HnTc2t2+24HnTct-4Hn,
where, Hz is the raise height of the swing leg; Hn is the expected maximum raise height of the ankle joint in the swing leg relative to the trunk fixed coordinate system in Zc direction; Tc is the step period.The support foot keeps on touch with the terrain. In the second step period, the swing leg changes to the support foot.


## 3.1.  “Mobile Equilibrium Method” in the Lateral Plane

The three-dimensional linear inverted pendulum mode (3D-LIPM) combined with zero moment point criterion are widely used in the biped robot control field. But the 3D-LIPM only fits for the biped robots with active ankle joints and whose legs are light compared with the upper part of the body. This paper proposes an easier walking pattern method which has no limits to the robot mass distribution. It also fits for the robot without active ankle joint in the lateral plane.We consider the robot’s center of mass (COM) as the research object. When feet move accelerated toward left relative to the trunk in the lateral plane, the inertia force in the robot’s COM points to right (shown as Figure4 left). If the acceleration is large enough and the ground can provide enough friction and support force, the force between the left foot and the ground becomes zero, and the trunk moves right. If the roll joint of the ankle in the lateral plane is active, in the ideal condition the trunk moves right without trunk roll motion (shown in Figure 4, middle); if the ankle does not have roll joints, the trunk’s roll angle is large when the trunk move right (shown as Figure 4 right); if the ankle has roll joints but they are passive, the trunk’s roll angle exists, but it is smaller.Figure 4
Force analysis in the lateral plane.The biped robot lower-limb introduced in this paper does not have roll joint in the ankle, that is to say the roll angle of the trunk must exist when the feet move accelerated to left/right relative to the trunk (suppose that the ground can provide enough friction force, this is the precondition). If the trunk’s roll angle changes periodically, it shows that the robot is stable in the lateral plane; otherwise, we should make some adjustments.When the feet move accelerated to left relative to the trunk with appropriate acceleration, the support force on the left foot decreases to zero, then the left foot lifts, and the trunk will not dump in the lateral plane. The theory is fit for the opposite direction. In order to keep the robot stable in the lateral plane, the following two conditions are needed.(1)
The initial acceleration inYc direction is adequate. If the acceleration is too small the swing foot cannot lift or the robot topples toward the swing leg; if the acceleration is too large, the robot topples toward the support leg.(2)
In the moving process, the resultant moment is zero which is created by the inertial force and the gravity relative to equilibrium pointO1 (the equilibrium is O1, when the robot COM moves to right; the equilibrium is O2, when the robot COM moves to left. Actually O1 and O2 are zero moment points).In the above process, the acceleration of the trunk keeps on changing with the displacement changes in the lateral plane, so the method is named as “mobile equilibrium method.”


## 3.2. Walking Pattern Generation inYc Direction

Two step periods are considered as one movement circle. Based on the “mobile equilibrium method”, we can obtain the walking pattern inYc direction.In the first step period, the feet move accelerated toward left relative to the trunk(1)mg(-Lf1-y)+my¨zc=0.
Then,
(2)y=c1et·g/zc+c2e-t·g/zc+Lf1,
where
(3)t∈[0,Tn],c1=-e-Tng/zc·Lf1eTng/zc-e-Tng/zc,c2=eTng/zc·Lf1eTng/zc-e-Tng/zc,t∈[Tn,Tc],c1=e-Tng/zc·Lf1eTcg/zc·e-Tng/zc-e-Tcg/zc·eTng/zc,c2=-eTng/zc·Lf1eTcg/zc·e-Tng/zc-e-Tcg/zc·eTng/zc.In the second step period, the feet move accelerated toward right relative to the trunk:(4)mg(Lf2-y)+my¨zc=0.
Then,
(5)y=c3et·g/zc+c4e-t·g/zc-Lf2.
While t∈[Tc,Tc+Tn],
(6)c3=e-(Tc+Tn)g/zc·Lf2e(Tc+Tn)g/zc·e-Tcg/zc-e-(Tc+Tn)g/zc·eTcg/zc,c4=-e(Tc+Tn)g/zc·Lf2e(Tc+Tn)g/zc·e-Tcg/zc-e-(Tc+Tn)g/zc·eTcg/zct∈[Tc+Tn,2Tc],c3=e-(Tc+Tn)g/zc·Lf2e(Tc+Tn)g/zc·e-2Tcg/zc-e-(Tc+Tn)g/zc·e2Tcg/zc,c4=-e(Tc+Tn)g/zc·Lf2e(Tc+Tn)g/zc·e-2Tcg/zc-e-(Tc+Tn)g/zc·e2Tcg/zc.m is the mass of the whole biped lower limb; g is the gravitational acceleration; y is the displacement of the robot COM relative to the ground, which is opposite to the displacement of the ankle joint relative to the body fixed coordinate system in Yc direction; y¨ is the acceleration of the robot COM relative to the trunk fixed coordinate system in Yc direction; zc is the absolute distance between the robot COM and the ground; Lf1 is the absolute distance between O1 and Oc in Yc direction; Lf2 is the absolute distance between O2 and Oc in Yc direction. If the robot is symmetrical, then Lf1=Lf2=Lf, actually if there is no sliding between the feet and the ground. Lf is the maximum displacement of the trunk COM in the lateral plane, and we named it as maximum offset in the lateral plane.


## 3.3. Walking Pattern Generation inXc Direction

The direction changes of the trunk’s acceleration will cause the loading sway back and forth. The larger the acceleration, the bigger the sway range is. In order to decrease this impact, the trunk is set to a constant velocity inXc direction (the COM adjustments change the velocity of the trunk slightly). That is to say, the ankle joint moves in an approximately constant velocity relative to the trunk fixed system, and then the trunk can keep a constant velocity relative to the ground. When the swing foot changes to the support foot, the swing foot should keep the same velocity as the trunk. With this constraint one of the easiest gait generation methods for swing leg is cubic polynomial fitting. One periodic motion contains two step periods.In the first step period, the trajectory of the ankle joint of the swing leg relative to the trunk fixed coordinate system inXc direction is shown as follows:
(7)xs=-8DsTc3t3+12Ds-DdTc2t2-2Ds-DdTct-Ds.The trajectory of the ankle joint of the support leg relative to the trunk fixed coordinate system inXc direction is shown as follows:
(8)x=-2DsTct+Ds,
where xs is the ankle displacement of the swing leg relative to the trunk fixed coordinate system in Xc direction; x is the ankle displacement of the support leg relative to the trunk fixed coordinate system in Xc direction; Ds is the step length of the ankle joint relative to the trunk fixed coordinate system in Xc direction; Dd is the increment of the step length between the two contiguous periods: if the step length increases, Dd is a positive, and if the step length decreases, Dd is a negative; Tc is the step period.The walking generation method in the second step period is the same as in the first step period, but the swing leg changes to the support leg and the support leg changes to the swing leg.


## 3.4. Walking Pattern Generation inZc Direction

In order to decrease the impact, when the swing foot hit the ground and changes to the support foot, the following two conditions are needed:(1)
the ideal dropping velocity inZc direction is zero;(2)
the raise velocity decreases to zero when the swing foot reaches the planned raise height.With the constraints, cubic polynomial fitting can be used here. In the first step period, the trajectory of the ankle in the swing leg relative to the trunk inZc direction is shown as follows:
(9)t∈[0,Tc2],Hz=-16HnTc3t3+12HnTc2t2,t∈[Tc2,Tc],Hz=16HnTc3t3-36HnTc2t2+24HnTct-4Hn,
where, Hz is the raise height of the swing leg; Hn is the expected maximum raise height of the ankle joint in the swing leg relative to the trunk fixed coordinate system in Zc direction; Tc is the step period.The support foot keeps on touch with the terrain. In the second step period, the swing leg changes to the support foot.


## 4. Balance Control

### 4.1. Motion Adjustment inYc Direction

From (3), (6), We can see that c1,c2,c3,  andc4 are constants, which are relative to zc,Tc, and Lf.As mentioned above,Tc is the step period; zc is the absolute distance between the robot COM and the ground, and we can obtain zc from experiments, but here we use the absolute distance between the origin point of the trunk fixed coordinate system and the ground z to replace, which is easier to obtain. The difference is compensated by adjusting Lf, with a different Lf, the trunk’s acceleration, velocity, and displacement in Yc direction changed accordingly. From the marking time (walking but the step length is zero) experiments, we may obtain the relationship between the loading weight, step period, and Lf, which could be retained in the following movements.When adjustingLf (shown as Figure 5), we care for the actual raise height Hz of the swing leg rather than ZMP position in Yc direction (that is to say, we do not care whether the equilibrium point O1 or O2 is in or out the convex polygon composed by the support foot). For the robot without active rolling joint in the ankle, the roll motion of the trunk is inevitable. The actual raise height Hz of the swing leg is larger than the planned He. We should guarantee that Hz≈k·He (k is more than one and obtained from experiments. It is related to the robot system and has little with the loading weight or the motion characteristics of the robot).(i)
IfLf is too small the swing leg cannot leave the ground normally or the robot topples toward the swing leg then the actual raise height Hz of the swing leg is less than the expected raise height k·He;(ii)
IfLf is too large, the robot topples toward the support leg Hz and the actual raise height is larger than k·He.Figure 5
Diagram of the motion adjustment inYc direction.


### 4.2. Motion Control in the Frontal Plane

With the motion adjustment in the lateral plane, the robot walking system can walk front and back, left and right with a small step length. In order to increase the impact resistance, cylinders in the frontal plane are controlled by the hydraulic actuated forcefe which is used to servo the elongation lene of the hydraulic cylinder. Shown in Figure 6, finitial is used to meet the offset, and kp and kd are obtained from experiments.Figure 6
Position servo.


### 4.3. Balance Control on Uneven Terrain

Unlike human size biped robot, when the hydraulic actuated biped robot with load walks dynamically, the ground reaction force distributed unevenly and varied in a large dynamic range (shown as Figure7). The ground reaction force to the physical prototype is smaller in a certain extent, because there is elastic mechanism on the feet which provides a buffer.Figure 7
Diagram of the ground reaction force.ZMP information for the biped robot walking on the uneven terrain is important, but it is hard to calculate it precisely. Obviously, we can obtain it by six-dimensional force sensor. But the larger the necessary measure range is, the more expensive the sensor is.From the kinetic characteristics of the biped robot lower limb we find that trunk attitude, especially the pitch angle and the roll angle, could be used as controlled variable. In the actual environment, the roll and pitch movements of the trunk cannot be eliminated at least for the robot with foot, but without active roll joint in the ankle. We use this information to realize the stable walking of the biped robot lower limb (shown in Figure8).(1)
When the trunk rolls alongXc in an active angle, the length of the right leg is shorter than expected, while the left one is longer than expected, we should increase the length of the right leg and decrease the length of the left leg. But the adjustment amount Hδ is the same, which is relative to the roll angle γ and the angular velocity γ˙ of the trunk; when the trunk rolls along Xc in a negative angle, and we should decrease the length of the right leg and increase the length of the left leg accordingly.(2)
When the trunk rolls alongYc in an active angle, the trunk tends to fall forward; in particular the robot walks down slope. We should move the COM backward relative to the ground (the absolute displacement of the support leg x relative to the trunk fixed coordinate system increases Pδ); when the trunk turns along Yc in a negative angle, the trunk tends to fall backward, and in particular the robot walks upslope. Then we should move the COM forward (the absolute displacement of the support leg x relative to the trunk fixed coordinate system decreases Pδ). Pδ is relative to the pitch angle β and angular velocity β˙ of the trunk (shown in Figure 8). It is important that kp1,kd1, kp2, and kd2, are determined by the biped robot lower limb itself, they have little to do with the loading weight and the step period. kp1,kd1,kp2, and kd2 are obtained from simulation experiments and have some changes in the physical prototype experiments.Figure 8
Structure of balance controller on uneven terrain.


## 4.1. Motion Adjustment inYc Direction

From (3), (6), We can see that c1,c2,c3,  andc4 are constants, which are relative to zc,Tc, and Lf.As mentioned above,Tc is the step period; zc is the absolute distance between the robot COM and the ground, and we can obtain zc from experiments, but here we use the absolute distance between the origin point of the trunk fixed coordinate system and the ground z to replace, which is easier to obtain. The difference is compensated by adjusting Lf, with a different Lf, the trunk’s acceleration, velocity, and displacement in Yc direction changed accordingly. From the marking time (walking but the step length is zero) experiments, we may obtain the relationship between the loading weight, step period, and Lf, which could be retained in the following movements.When adjustingLf (shown as Figure 5), we care for the actual raise height Hz of the swing leg rather than ZMP position in Yc direction (that is to say, we do not care whether the equilibrium point O1 or O2 is in or out the convex polygon composed by the support foot). For the robot without active rolling joint in the ankle, the roll motion of the trunk is inevitable. The actual raise height Hz of the swing leg is larger than the planned He. We should guarantee that Hz≈k·He (k is more than one and obtained from experiments. It is related to the robot system and has little with the loading weight or the motion characteristics of the robot).(i)
IfLf is too small the swing leg cannot leave the ground normally or the robot topples toward the swing leg then the actual raise height Hz of the swing leg is less than the expected raise height k·He;(ii)
IfLf is too large, the robot topples toward the support leg Hz and the actual raise height is larger than k·He.Figure 5
Diagram of the motion adjustment inYc direction.


## 4.2. Motion Control in the Frontal Plane

With the motion adjustment in the lateral plane, the robot walking system can walk front and back, left and right with a small step length. In order to increase the impact resistance, cylinders in the frontal plane are controlled by the hydraulic actuated forcefe which is used to servo the elongation lene of the hydraulic cylinder. Shown in Figure 6, finitial is used to meet the offset, and kp and kd are obtained from experiments.Figure 6
Position servo.


## 4.3. Balance Control on Uneven Terrain

Unlike human size biped robot, when the hydraulic actuated biped robot with load walks dynamically, the ground reaction force distributed unevenly and varied in a large dynamic range (shown as Figure7). The ground reaction force to the physical prototype is smaller in a certain extent, because there is elastic mechanism on the feet which provides a buffer.Figure 7
Diagram of the ground reaction force.ZMP information for the biped robot walking on the uneven terrain is important, but it is hard to calculate it precisely. Obviously, we can obtain it by six-dimensional force sensor. But the larger the necessary measure range is, the more expensive the sensor is.From the kinetic characteristics of the biped robot lower limb we find that trunk attitude, especially the pitch angle and the roll angle, could be used as controlled variable. In the actual environment, the roll and pitch movements of the trunk cannot be eliminated at least for the robot with foot, but without active roll joint in the ankle. We use this information to realize the stable walking of the biped robot lower limb (shown in Figure8).(1)
When the trunk rolls alongXc in an active angle, the length of the right leg is shorter than expected, while the left one is longer than expected, we should increase the length of the right leg and decrease the length of the left leg. But the adjustment amount Hδ is the same, which is relative to the roll angle γ and the angular velocity γ˙ of the trunk; when the trunk rolls along Xc in a negative angle, and we should decrease the length of the right leg and increase the length of the left leg accordingly.(2)
When the trunk rolls alongYc in an active angle, the trunk tends to fall forward; in particular the robot walks down slope. We should move the COM backward relative to the ground (the absolute displacement of the support leg x relative to the trunk fixed coordinate system increases Pδ); when the trunk turns along Yc in a negative angle, the trunk tends to fall backward, and in particular the robot walks upslope. Then we should move the COM forward (the absolute displacement of the support leg x relative to the trunk fixed coordinate system decreases Pδ). Pδ is relative to the pitch angle β and angular velocity β˙ of the trunk (shown in Figure 8). It is important that kp1,kd1, kp2, and kd2, are determined by the biped robot lower limb itself, they have little to do with the loading weight and the step period. kp1,kd1,kp2, and kd2 are obtained from simulation experiments and have some changes in the physical prototype experiments.Figure 8
Structure of balance controller on uneven terrain.


## 5. Experimental Verification

### 5.1. Simulation Experiments

It is assumed that the terrain information has been contained in the walking generation. If the step height is less than 20 mm or the gradient of the slope is less than±3∘, the road could be seen as flat road. But we should know when the robot walks up/down steps, because the robot does not have toe, and then it is forbidden that less than half of the support foot is on the step.The robot loads a block whose weight is 100 kgf. Walking across a step which is 20 mm in height and two slopes (the gradients are3∘and -3∘, resp.). At the same time there are two disturbance forces on the center of the loading pointing to Xc and Yc, respectively. The disturbances are added manually by two cosine waves whose amplitude is 40 N. With the proposed adjusting method, the biped robot lower limb walks on the uneven terrain under disturbance successfully (shown in Figure 9). When walking without the balance control method proposed in this paper the robot cannot even pass the step which is 20 mm in height (shown in Figure 10). For the robot having big foot which is 270 mm in length, when the robot walks on steps, it is advised to use bigger step length, but when the robot walk on slopes, it is advised to use smaller step length. For example, when the robot walking on slopes the step length is 240 mm and the step period is 0.4 s (shown in Figures 9 and 10); when walking up steps which is 60 mm, the step length is 300 mm and the step period is 0.4 s (shown in Figure 11).Figure 9
Walking with the balance control method on uneven terrain.Figure 10
Walking without the balance control method on uneven terrain.Figure 11
Walking up steps.


### 5.2. Physical Prototype Experiments

Without power producer, the weight of the original robot is 40 kg, but the mass of one single leg is about 14 kg, and the leg’s mass cannot be ignored. Obviously, it is not fit for the assumed condition of LIPM. With the “dynamic equilibrium method,” the biped lower-limb can walk along, walk back, and walk in the lateral plane successfully. With the proposed control method in the paper, the robot pass through a carpet which is about 10 mm in thickness, a rubber blanket which is about 8 mm in thickness, and a plank which is about 13 mm in thickness successfully (shown in Figure12). For the robot walking without load, the robot mass is concentrated on the legs. Then the position changes of the leg may cause body COM position changing obviously, so the step length of the physical prototype in frontal plane is 60 mm and the step period is 0.45 s.Figure 12
Physical prototype experiments.


## 5.1. Simulation Experiments

It is assumed that the terrain information has been contained in the walking generation. If the step height is less than 20 mm or the gradient of the slope is less than±3∘, the road could be seen as flat road. But we should know when the robot walks up/down steps, because the robot does not have toe, and then it is forbidden that less than half of the support foot is on the step.The robot loads a block whose weight is 100 kgf. Walking across a step which is 20 mm in height and two slopes (the gradients are3∘and -3∘, resp.). At the same time there are two disturbance forces on the center of the loading pointing to Xc and Yc, respectively. The disturbances are added manually by two cosine waves whose amplitude is 40 N. With the proposed adjusting method, the biped robot lower limb walks on the uneven terrain under disturbance successfully (shown in Figure 9). When walking without the balance control method proposed in this paper the robot cannot even pass the step which is 20 mm in height (shown in Figure 10). For the robot having big foot which is 270 mm in length, when the robot walks on steps, it is advised to use bigger step length, but when the robot walk on slopes, it is advised to use smaller step length. For example, when the robot walking on slopes the step length is 240 mm and the step period is 0.4 s (shown in Figures 9 and 10); when walking up steps which is 60 mm, the step length is 300 mm and the step period is 0.4 s (shown in Figure 11).Figure 9
Walking with the balance control method on uneven terrain.Figure 10
Walking without the balance control method on uneven terrain.Figure 11
Walking up steps.


## 5.2. Physical Prototype Experiments

Without power producer, the weight of the original robot is 40 kg, but the mass of one single leg is about 14 kg, and the leg’s mass cannot be ignored. Obviously, it is not fit for the assumed condition of LIPM. With the “dynamic equilibrium method,” the biped lower-limb can walk along, walk back, and walk in the lateral plane successfully. With the proposed control method in the paper, the robot pass through a carpet which is about 10 mm in thickness, a rubber blanket which is about 8 mm in thickness, and a plank which is about 13 mm in thickness successfully (shown in Figure12). For the robot walking without load, the robot mass is concentrated on the legs. Then the position changes of the leg may cause body COM position changing obviously, so the step length of the physical prototype in frontal plane is 60 mm and the step period is 0.45 s.Figure 12
Physical prototype experiments.


## 6. Conclusion

Human size biped robot actuated by electrical motor developed rapidly and became mature gradually. BigDog [17] and Petman [18] made by Boston Dynamics bring hydraulic driving technology to robot study again. The paper gives an outline of a hydraulic actuated biped robot lower limb with load ability. The biped robot lower limb only has 8 active freedoms, for each leg there are four, two in the hip, one in the knee and one in the ankle.3D-LIPM with ZMP criterion is widely used in the biped robot control. Considering the boundedness of the LIPM model and the complexity of ZMP, the paper proposes an easy walking generation method and balance control method. “Dynamic equilibrium method” and the cubic polynomial fitting method are used in the walking generation. This method is fit for both the active and underactive biped robots with foot. Not as LIMP, “dynamic equilibrium method” proposed in this paper has no limit to the mass distribution between the trunk and the legs. The paper also proposes a simple control method without using the ZMP information directly. The biped robot lower-limb adjusts the maximum offset in the lateral plane according to the actual raise height of the swing leg in the marking time; the biped robot adjusts the COM position inXc direction, according to the pitch angle and angular velocity without using ZMP information directly.With the proposed method the biped robot lower limb with a box weight 100 kgf can walk stably on uneven terrain under the disturbance whose amplitude is 40 N. The control method proposed in this paper is also effective to the hydraulic actuated biped robot lower limb which does not fit for the assumed conditions of LIPM. When weighting empty, with the proposed method in the paper, the hydraulic actuated biped robot lower limb can walk on the uneven terrain successfully.



---

*Source: 429458-2013-09-01.xml*