Walking Stability Compensation Strategy of a Small Humanoid Robot Based on the Error of Swing Foot Height and Impact Force

Abstract

Abstract In order to reduce the impact force of swing legs and improve walking stability when a small humanoid robot is walking, a set of impact dynamics equations based on the second kind Lagrange equation is produced, and an impact compensation control strategy with a BP network optimized by a particle swarm algorithm is designed. The core element of the compensation controller is replacing the error back propagation with a particle swarm algorithm. Due to the regulating joints of the knee, hip and ankle, the walking process is more stable than before. The experiment results show that when the left swing leg lands, the impact force drops by 2N and 1.5N respectively in the moments 4.5s and 10.5s. Therefore, the compensation strategy can reduce the impact force effectively and improve the walking stability.

Keywords

Humanoid Robot BP network Particle Swarm Optimization Compensating Controller Impact Dynamics

1. Introduction

With the promotion of robots' ability to simulate human behaviour, reducing the impact of the landing of the swing leg whilst walking is increasingly prominent for robot walking stability. Therefore, being able to establish impact dynamics equations and compensate for joint angles during the walking process in order to improve walking stability has become an urgent area of study around the world. At present, many institutes are doing research in this area and have achieved a considerable amount.

Based on the concept of a change in dynamics, Ken'ichiio N proposed a coordinated movement control [1]. By using this method, dynamics models can be generated synchronously to maintain balance for smooth walking. The effectiveness of the method has been validated by QRIO which was introduced by SONY. Honda created ASIMO which had a self-discipline movement technology [2], developed the I-WALK movement model, which can effectively reduce impact force and regulate walking processes by predicting the centre of gravity. Kenji H of Waseda University changed the foot structure in order to adapt to walking environments so as to reduce the impact force [3–4]. This special foot structure may absorb vibrations when the foot hits the floor. M2, which was made by MIT, can walk more naturally. A ball screw was designed to drive a spring joint for monitoring the torque at all times [5]. The compensation controller was installed on the THBIP-II made by Tsinghua University [6]. The height error and pace error of the swing leg are regarded as the input, whereas the three joint compensation angles are the output. The National University of Defense Technology designed four generations of humanoid robots with orthogonal axes [7]. The strategy of the BHR-2, a robot created by the Beijing Institute of Technology, is to utilize sensors to interact with the environment, thus reducing the impact with the use of slow foot regulation. As well as this, an approach to revising sensor system errors on feet was also proposed [8]. Moreover, other researchers make use of flexible foot structures to reduce impact force. For example, a cushion material was installed into the foot panel so as to offset impact energy [9–10]. A foot structure that can adapt to the walking environment can be directly designed without feedback components [11–13].

To sum up, there are mainly three means to compensate for the joint angle. (1) Utilizing passive flexibility to reduce impact force directly. (2) Revising moving parameters according to the environment. (3) Regulating joint angles based on dynamics models.

The above measures can modify walking characteristics to improve walking quality. Measure (1) seems simple but cannot be used in most fields. Measure (2) can increase the flexibility of a robot, but walking gait cannot solve every kind of impact problem. Measure (3) can revise walking parameters according to sensors; it can also be valid in complicated environments. Therefore, based on past achievements, this paper formulates a Particle Swarm Optimization and BP network (PSO-BP) to compensate for impact force in order to improve walking stability.

2. The Robot Prototype

The most notable difference between humanoid robots and other forms of robots is that humanoid robots' feet can perform comfortable walking. Its main structure is designed based on a human body structure for this version as is each part of the robot's size, quality distribution and DOFs design selection, etc.

Based on a human body structure, a 3D-model of the prototype was established by Solidworks. As shown in Figure 1, when the robot is standing upright, the height from head to feet is 375 mm, the width between the left and right shoulder is 289 mm, and the body weighs about 2.2 kg. To be precise, lower limb length (from feet to hip) is about 206 mm, ankle height is about 31 mm.

Figure 1.

The 3 dimensional model of a small humanoid robot

In addition, the number and allocation of the DOFs' setting are just like the human distribution of bone and muscle tissue. DOF distribution is shown in Figure 2:

Figure 2.

The layout of DOF of humanoid robot

This research is focused on the lower leg mechanism, which emphasizes the distribution of the lower limbs DOF: Hip 2: Pitch direction + Roll direction, and knee a Pitch direction, ankle 2: Pitch direction + Roll direction.

After the sensor test is finished, sensors are installed on the lower limbs of the robot; at the same time, the potentiometer shaft end and steering gear shaft inside are linked together to measure the real-time turning angle and a pressure sensor in the robot is pasted on the foot motherboard. The strain type pressure sensor and potentiometer can be seen from the picture on the right in figure 3.

Figure 3.

The prototype of a small humanoid robot

3. Collision Dynamics Equations

The Lagrange equation of the second category with the potential or conservative system is [14]:

\frac{d}{d t} \frac{\partial L}{\partial {\overset{•}{q}}_{j}} - \frac{\partial L}{\partial q_{j}} = 0 (j = 1, 2 \dots k)

(1)

Where:

L = T - V

(2)

L is the Lagrange function, T is the total kinetic energy of system, V is the total potential energy, q_j are generalized velocities and q_j is the generalized coordinates.

The small humanoid robot legs institutions involved in the process of moving can be treated as a rigid body. The position and orientation of the leg's modelling do not change in the process of the walking impact. Its shape and inertia will not change either. Momentum F towards the ground to swing the foot in the collision process is a quantity with a constant size and direction. Therefore, the momentum in this process can be powerful. It can be shown in a function which describes their potential relative to a basis point. So Equation (3) can be obtained from Equation (1).

Δ (\frac{\partial T}{\partial {\overset{•}{q}}_{j}}) + \frac{\partial V_{1}}{\partial q_{j}} = 0 (j = 1, 2 \dots, k)

(3)

V₁ is the impulse potential.

After the implicit equation introduced, the explicit equation should be established to make clear the parameter of the equations. Therefore, the D-H coordinate system has to be established, as is depicted in Figure 4.

Figure 4.

The D-H coordinate system of the humanoid robot

Based on the robot prototype body and the D-H coordinate system described above, coordinate transformation matrixes A_i+1 are created to describe the relationship between the adjacent coordinate system [15]. According to the structure of the humanoid robot and the established D-H coordinate system, the coordinate transformation matrixes between A₁ to A_i can be expressed as:

A_{i} = A_{1} \cdot A_{2} \cdot \cdot \cdot \cdot \cdot A_{i} (i = 1, 2, …, 10)

(4)

Each rigid bar i has its own centroid. Centroid position r_i can be described as below:

r_{i} = (\begin{matrix} x_{i} \\ y_{i} \\ z_{i} \\ 1 \end{matrix}) = {(x_{i}, y_{i}, z_{i}, 1)}^{T}

(5)

Then, each bar's centroid position r_i relative to the origin of the coordinates:

^{0} r_{i} =^{0} A_{i} \cdot r_{i}

(6)

Therefore, the centroid speeds relative to the origin of the coordinates:

^{0} V_{i} = \frac{d}{d t} (^{0} A_{i} \cdot r_{i}) = [\sum_{j = 1}^{i} \frac{\partial^{0} A_{i}}{\partial q_{j}} {\overset{•}{q}}_{j}] r_{i}

(7)

Set dT_i is on the quality of the micro elements dm kinetic energy, the components i(i = 1, 2, …, k) of kinetic energy:

{dT}_{i} = Tr (V_{i} V_{i}^{T}) dm / 2

(8)

Tr is the trace of the matrix.

Put V_i generation on the type and finishing. By introducing

W_{i} = \int r_{i} r_{i} ​^{T} d m_{​}

(9)

To simplify the type so that:

T_{i} = \frac{1}{2} Tr \int {(V}_{i} V_{i}^{T}) dm

(10)

Among them W_i only depends on the quality distribution of component i, and has nothing to do with the position and movement speeds [16].

Therefore, the total kinetic energy of the system is:

T = \sum_{i = 1}^{n} T_{i}

(11)

Thus far, the total kinetic energy of the system has been established. The next step is to consider introducing the characteristics of the potential of the momentum.

Through the analysis of the characteristics of a humanoid robot prototype entity in the process of collision in walking, with regards to component 0 (the support legs) as a reference body, there is a force on component 6 with a known direction and strength, namely the impact on the swinging leg caused by the ground while landing.

P₆ is set as the impulse potential on component 6, namely the impact on the swinging leg caused by the ground while landing. Then,

P_{6} = - I_{6}^{0} A_{6} r_{6}

(12)

r₆ is the coordinates of the impact point. Therefore, in the moment that the swinging leg lands, the lower limbs of the robot are affected by only one impulse; so the total impulse V₁ is equal to P₆.

V_{1} = P_{6} = - I_{6}^{0} A_{6} r_{6}

(13)

From the formula, the kinetic energy and impulse potential of the lower limb of a small humanoid robot, respectively, can be obtained. They are applied in analyzing the form of the collision equation. So, the display equation can be achieved.

Δ (\frac{\partial T}{\partial {\overset{•}{q}}_{j}}) = \sum_{j = 1}^{n} \sum_{k = 1}^{j} T r (M_{j k} \cdot W_{j} \cdot M_{j i} ​^{T}) \cdot {\overset{​}{v}}_{k}

(14)

Among them:

{\overset{​}{v}}_{k} = {\overset{•}{q}}_{k^{0}} - {\overset{•}{q}}_{k}

(15)

${\overset{•}{q}}_{k}, {\overset{•}{q}}_{k^{o}},$ stand for the feet speed before the walking impact and after walking impact. It is obvious that after the walking impact the foot speed is 0, namely ${\overset{•}{q}}_{k^{0}} = 0$ . And M_ij = ∂⁰ A_i/∂q_j. Therefore,

{\overset{​}{v}}_{k} = - {\overset{•}{q}}_{k} = - {\overset{•}{q}}_{6}

(16)

So generations into the Lagrange form of the collision equation are:

\sum_{j = 1}^{n} \sum_{k = 1}^{j} T r (M_{j k} \cdot W_{j} \cdot M_{j i} ​^{T}) \cdot {\overset{​}{v}}_{k} - I_{6} M_{6 i} r_{6} = 0 (i = 1, 2, …, 6)

(17)

To conclude, all the theory collision dynamics of a humanoid robot equation have been established. The next section will involve the theoretical data needed to conduct research into the compensation controller.

4. Compensation Controller Design

Before designing the compensation controller, the most important step is to calculate the swing foot height in the moving plane. The limb body structure is simplified by the model shown in Figure 5.

Figure 5.

Marching and positive plane figure

According to Figure 5, the distance between each foot plane l₆ to the ground is deemed as the swing foot height. It is obvious from figure 5, l₆ is the ankle height (31mm), l₁ is the lower leg height (85mm), and l₂ is the higher leg height (90mm). Using a simplified model in the Figure 5, formula(18) is established to calculate the swing foot height during walking [17].

\begin{array}{l} h = (l_{1} \cdot \cos q_{A} + l_{2} \cdot \cos (q_{B} - q_{A})) \cdot \cos θ_{A} + l_{0} - l_{2} \cdot \cos (q_{C} - (q_{C 1} - q_{B} + q_{A})) \\ \cdot \cos (θ_{C} - θ_{B} + θ_{A}) + l_{1} \cdot \cos (q_{D} - q_{C} + (q_{C 1} - q_{B} + q_{A})) \cdot \cos (θ_{C} - θ_{B} + θ_{A}) \\ - l_{0} \cdot \cos (q_{E} - q_{D} + q_{C} - (q_{C 1} - q_{B} + q_{A})) \cdot \cos (θ_{D} - θ_{C} + θ_{B} - θ_{A}) \end{array}

(18)

According to the formula (18) and the derivation of the definition, the expected gait programming data joints angle of the freedom and the type are substituted, and the swing expectations of the height of the feet movement tracking can be formulated.

We analyze the swing foot height error and the impact force error, which are caused by many servos' mutual reactions and the moving error. In addition, because of the inevitable error between the robot body and the simplified model as well as the servos' moving track error, a swing foot height error exists. Therefore, this paper focuses on compensating for the relative angle to eliminate the swing foot height error and impact force error.

This paper presents hip, knee and ankle joint angles to compensate for the swing foot height error and impact force error. According to impact dynamics equation and the formula (18), there is a nonlinear relationship between the feet height, the impact force and the steering joint angle. There will be a lot of complex inverse kinematics calculations to calculate the compensated joint angle. Therefore, the swinging leg compensating controller is designed based on the PSO optimization of the BP neural network in this paper, as shown in Figure 6.

Figure 6.

The controller structure of contact swing legs

The process of humanoid robot compensation for impact control is as follows:

Step 1: Servo control board output angle e_swknd to the knee joint, while in the process of walking, the real knee joint angle is e_swkna, e_swkn is the error between e_swknd and e_swkna, which is defined as a parameter of the control rules. The structure of the control rules will be described in the following text.

Step 2: The error ΔF is between the theoretical impact force F_S and the real impact force F during walking. ΔF is one of the BP controller's input parameters.

Step 3: In the calculation of the real swing foot height H₁ by function (18), H_S is the theoretical swing foot height. The error ΔH₁ is another input parameter of the BP controller.

Step 4: Calculating the adjustment knee incremental θ_rkn using the BP controller, whose inputs are ΔF and ΔH₁, is defined as another control rules parameter.

Step 5: After the “control rules” inputs are set; the actual compensated knee joint angle θ_rkn can be calculated through its control law; “Control rule” will be introduced in table 1.

Table 1.

The control rules

Relation Types	Actual Compensation Angle θ_tkn
e_swkn ≤ θ_rkn ≤ 0 0 ≤ θ_rkn ≤ e_swkn	θ_tkn = θ_swknd
θ_rkn ≤ 0 ≤ e_swkn e_swkn ≤ 0 ≤ θ_rkn	θ_tkn = θ_swknd + θ_rkn
0 ≤ e_swkn ≤ θ_rkn θ_rkn ≤ e_swkn ≤ 0	θ_tkn = θ_swknd + θ_rkn – e_swkn

Step 6: Similarly, e_swhip is the error between the designed angle of the servo control board to the hip joint θ_swhipd and the real hip joint angle θ_swhip, which is one of the control rules parameters.

Step 7: The error ΔF is the difference between the theoretical impact force F_S and the real impact force F during walking. ΔF is one of BP controller input parameters.

Step 8: After compensation of the knee joint, the real swing foot height is changed toH₂, while the theoretical values are still H_S. Consequently, the error is ΔH₂ for the BP controller's input parameter.

Step 9: Calculating the adjusted hip incremental θ_rhip using a BP controller, with the inputs ΔF and ΔH₂, which are defined as another control rules parameter.

Step 10: After the “control rules” inputs are set; the actual compensated hip joint angle θ_thip can be calculated using its control law.

Step 11: Owing robot's upper limb being upright, the ankle angle will be regulated by the compensated knee and hip angle in accordance with the function θ_tak = –(θ_thip + θ_tkh), and the compensated ankle angle will be directly calculated.

Step 12: The compensated knee, hip and ankle angle are directly entered into the servo control board as inputs, and the compensation process is accomplished.

The paper will present the control rules as follows.

From the above steps regarding the controller process during walking, the servo control board output angle e_swknd, and the actual compensated joint angle θ_rkn have the following 6 forms: corresponding to the eight kinds of form of “control rules” are the actual compensation angles. Among them: the hip control scheme and the knee control scheme are similar, only taking the knee joint compensation as an example, the control rules are given.

The algorithm of BP controller is drawn in Figure 7. Based on the analysis of the traditional BP neural network and the advantage of PSO algorithm, PSO-BP network is designed, in which PSO algorithm is used to optimize the parameter of BP network.

Figure 7.

PSO optimization algorithm of BP network flow chart

After the PSO-BP structure has been established, the relative parameters will be designed as follows.

First of all, the neural network hidden layer number is chosen. Many researchers have proved that if the hidden layer number of the neurons of a neural network is enough, the network can approach any continuous function with arbitrary precision. However, the research showed that two input parameters of the BP network produce an output fitting; so the single hidden layer of a network is applied in this paper.

Secondly, the number of hidden nodes is chosen. Because in different fields the designs of the neural network are different; experimental methods are usually involved in determining the optimal number of hidden nodes. Therefore, the Hecht-Nielsen law is used in this paper at first (one of them is a hidden layer node, the number of samples is the input dimension) and the number of hidden nodes increases, comparing different schemes of the training and test results.

Thirdly, a transfer function is chosen [18]. Since the final compensation angle of the neurons requires nonlinear fitting. Tan-Sigmoid functions are used in the hidden layer, whereas, linear functions (Purelin) are chosen in the output layer.

Finally, after data positive spreading, the error back propagation stage begins, in which the initialization network weights and threshold are optimized. The PSO algorithm is used to replace the traditional BP algorithm to optimize the neural network in weights and threshold. So, the performance defects that are brought about by the traditional BP algorithm can be avoided (such as: easy in minimum point, fitting, etc.).

From the above analysis, the single hidden layer with 5 nodes is the preliminary which is used to establish the BP network in this paper, and through the experimental proof it can be seen that the BP neural network has a better generalization ability with such a configuration.

The particle swarm parameters are set as follows:

Firstly, the swarm population scale is chosen. Owing to the general characteristics of this study, which does not belong to a difficult or certain type, the swarm population size is defined as 10.

Secondly, the length of the particles is connected to the network structure. After determining the BP network for the single hidden layer with the five nodes of the structure, particle length N = p×a + a×o + a + o = 21, can be determined, meaning that the length of each particle is 21 dimensional.

Next, v_max stands for the precision of the impact of the resolution between the current position and the best position. If v_max is too fast, the particle may miss the minimum point; If v_max is too slow, the particle may get stuck in the local extreme value area. Therefore, v_max is two based on experience.

Then the inertial factor wand learning factor c₁ and c₂ are chosen. The inertial factor w stands for the ability to extend the search space. c₁ and c₂ represent the acceleration weights when particles fly to P_best and G_best. Set the w = 4.1, and c₁ = c₂ = 2.

Finally, the termination conditions are set. When the error is lower than 0.005 and the maximum recycling times are 5000, the algorithm stops running.

Thus far, the parameters of BP controller have been set in Matlab software and the program is built after testing the performance analysis.

Based on the experimental data, an effective program is established. After the training of the network, a performance curve is acquired as shown in Figure 8.

Figure 8.

BP network performance training figure

Figure 8 shows that the expected precision can be achieved when the training time is 4 seconds and in 307 steps, the precision has declined. Before that, a BP controller is efficient and feasible. In order to further train the quantitative characteristics, we take a training of ten times as a benchmark and five as a hidden layer number of neurons, the specific parameters are listed in Table 2.

Table 2.

Training conditions of single hidden layer with 5 nodes BP network (MSE stands for the mean square error)

NO.	Time	Epoch	Speed	MSE
1	11s	1003	0.01	0.0367
2	9s	795	0.01	0.0611
3	12s	1109	0.01	0.0547
4	8s	691	0.01	0.0375
5	5s	485	0.01	0.0642
6	7s	583	0.01	0.0758
7	7s	610	0.01	0.0772
8	9s	773	0.01	0.0519
9	8s	703	0.01	0.0573
10	9s	802	0.01	0.0484

It can be seen from Table 2, when the learning rate is 0.01, the training times vary from 5 s to 12 s, the number of times of the training also differ from 485 times to 1109 times, and the fluctuation is obvious. Through calculating the square error of the testing results and the expected output, it can be seen that the mean square errors are all smaller than 0.1. The minimum mean square error is only 0.0367, the maximum mean square error is 0.0772. This suggests that the network structure design is ideal and the expected results can be obtained.

5. Walking Testing

Through the experimental testing, the control strategy is embedded in a humanoid robot during walking. At the same time, we can obtain the curve of the three joints changing with time, compensated before and after. Figures 9 to 11 are views of the curve of joint points of the knee, hip and ankle of the swinging left leg changing with time compensated before and after. Figures 12 to 14 are views of the curve of the points of three joints, changing with time compensated before, when the right leg swings during walking.

Figure 9.

Left knee angle curve of before & after compensation

Figure 10.

Left hip angle curve of before & after compensation

Figure 11.

Left ankle angle curve of before & after compensation

Figure 12.

Right knee angle curve of before & after compensation

Figure 13.

Right hip angle curve of before & after compensation

Figure 14.

Right ankle angle curve of before & after compensation

Figures 9 to 11 indicate that the three joints have a different degree of volatility because of collision impact, when the swinging leg lands (4.5s&10.5s), both before & after compensation. Joint fluctuations are obvious before the compensation; for example, the change in the angle of the knee-joint is 2 ˜ 3 °or so, in the second step in a crash landing; after compensation, we can see that the knee joint angle evidently decreased about 1° due to the influence of the impact. Similarly, the hip and ankle joint angle values have a small variation when the legs come into contact with the ground after compensation. This indicates that the influence of the impact is less when the left leg comes into contact with the ground.

In the same gait, the impact force influences the stability of the walk only when the right leg comes into contact with the earth at 7.5 s. So we should compensate for the three joint angles at 7.5 s particularly. In the same case as a left swinging leg, in figures 12 to 14 we can draw the conclusion that swinging the right leg at the time of 7.5 s will have obvious signs of compensation; there are different changes in the joint angles as well after a crash landing.

By compensating for the lower limbs of the robot's six joints, fluctuations of robot joint angles can be effectively reduced, making steering gears rotate more smoothly and not change largely when legs have a landing collision.

In the process of the gait of a walk after compensation, acquisition of the impact force of the feet takes place as shown in figures 15 and 16:

Figure 15.

Left foot contact force curve of before & after compensation

Figure 16.

Right foot contact force curve of before & after compensation

Through the collision impact collection of the swinging feet, the change of collision power before and after compensation can be seen; in the moments at 4.5 s and 10.5 s when the left swinging leg lands, the collision impact rises from zero to peak abruptly both before and after compensation; for example, the impact force decreases from about 17N before compensation to about 15.5N in the first step when the left swinging leg lands. After a short adjustment, it comes to 15N; this indicates that the walking robot has a smaller collision compact when landing after compensation and institutions ontology has a smaller effect than before. In the same way, the collision impact of the second step is even larger before compensation: about 17.5N; however, after compensation the peak is confined to about 16N.

6. Conclusions

The paper introduces research into a compensation control for the impact force of a small humanoid robot BJHR-II, which is aimed at solving the problem that the robot cannot walk stably after the swing leg collides with the ground. In order to explicitly reflect the information of the position and pose of the joints of the robot's lower limbs, the implicit and explicit equations for walking impact are established based on the second kind of Lagrange equations. Furthermore, a compensation strategy for impact of the swing legs during landing is derived from the error analysis of humanoid robots. In this strategy, the error of the swing leg is processed with BP networks and PSO, in which joint angles of the knee and hip can be compensated for based on the control rules of each joint. Based on the experiments, it is proved that the compensation strategy has a good performance on increasing the walking stability of humanoid robots.

7. Acknowledgments

This research is supported by the Fundamental Research Funds for the Central Universities (No.2013JBM080).

References

Nagasaka

Ken'ichiio

Kuroki

Yoshihiiro

Suzuki

Shin'ya

(2004) Integrated Motion Control for Walking, Jumping and Running on a Small Bipedal Entertainment Robot. IEEE Int. Con. on Robotics and Automation, USA, pp. 3189–3194.

Sakagami

Yoshiaki

Chiaki

Aoyama Ryujin Watanabe

(2002) The intelligent ASIMO: System overview and integration. IEEE/RSJ Int. Con. on Intelligent Robots and Systems, Switzerland, pp. 2478–2483.

Yoshikawa

Taizo

Khatib

Oussama

(2009) Compliant Humanoid Robot Control by the Torque Transformer. IEEE/RSJ Int. Con. on Intelligent Robots and Systems, USA, pp. 3011–3018.

Kaneko

Kanehiro

Kajita

S.H.

(2004) Humanoid Robot HRP-2. IEEE/RSJ Int. Con. on Robotics and Automation, USA, pp. 1083–1090.

Lei

Xusheng

Wang

Tianmaio

Chen

Dianshen

(2007) The Design and Control of a Low Cost Humanoid Robot. IEEE Int. Con. on Robotics and Biomimetics, China, pp. 763–768.

Zhangguo

Huang

Qiang

Jianxi

Chen

Xuechao

Kejie

(2008) Computer Control System and Walking Pattern Control for a Humanoid Robot. IEEE/ASME Int. Con. on Advanced Intelligent Mechatronics, China, pp. 1018–1023.

Park

Ill-Woo

Kim

Jung-Yup

Jun-Ho

(2006) Online Biped Walking Pattern Generation for Humanoid Robot KHR-3(KAIST Humanoid Robot - 3: HUBO). IEEE Int. Con. on Robotics and Automation, Italy, pp. 398–403.

Zhang

Hong-Yuan

Xing-Guo

(2013) Crank Round Slider Engine Multi-Flexible-Body Dynamics Simulation. Journal of Theoretical and Applied Information Technology, 47(2), pp. 575–579.

Jiandong

Zhao

Zishang

Wang

Chenglong

Ken

Chen

(2009) Compensating Control Based on Swing Foot Locomotion Error of Humanoid Robot for Diminishing Impact. Journal of Beijing Jiaotong University, 33(1), pp. 45–49.

10.

Ayaz

Konno

Munawar

(2011) A Human-Like Approach Towards Humanoid Robot Footstep Planning. International Journal of Advanced Robotic Systems, 8 (4). pp.98–109.

11.

Hongxu

Peng

Zhang

Lianqi

Zhang

(1997) Posture stability and Posture Control of Biped Dynamic Walking. Journal of Robot, 19(3), pp. 180–186.

12.

Leibin

Qinxin

Cao

Yijun

Sun

(2005) Discrete Trajectory Controlling of Humanoid Robot Based on Center of Pressure. Journal of Shanghai Jiaotong University, 41(8), pp. 1263–1266.

13.

Zongying

Shi

(2005) Research on gait planning and decentralized robust tracking control of biped robots. Beijing, Tsinghua University.

14.

Fengxiang

Mei

Jiping

Zhou

Xiaoping

Shui

(2005) Engineering Mechanics. High Education Press.

15.

Qiang

Ken

Chen

Liu

Chenglong

(2009) Research and Development of Kid-size Humanoid Robot THBIP-II. Journal of Robot, 31(6), pp. 586–593.

16.

Dingguo

Zhang

(1997) The Equations of External Impacted Dynamics between Multi-Rigibody Systems. Applied Mathematics and Mechanics, 18(6), pp. 551–555.

17.

Jiandong

Zhao

(2004) Research on Self-tuning Fuzzy Control of the Humanoid Robot Based on Walking Error. Beijing, Tsinghua University.

18.

Cai

Yao

Zhan

Qiang

(2011) Neural network control for the linear motion of a spherical mobile robot. International Journal of Advanced Robotic Systems, 8(4), pp. 79–87.