Locomotion approach of bipedal robot utilizing passive wheel without swing leg based on stability margin maximization and fall prevention functions

This paper proposes a locomotion approach of leg-wheel robot utilizing passive wheel attached to the foot of bipedal robot. The key feature of this approach is bipedal mobility without swing leg. This mobility contributes the stability based on expansion of support polygon during locomotion, the robustness for obstacles and stopping to prevent fall, and the adaptability by prevention of body swing sideways. To achieve these, we propose the stability margin maximization to optimize center of gravity projection for support polygon and the fall prevention functions for real environment that is a difficult situation to prevent unexpected fall by the only planning. Finally, we apply the proposed methods to leg-wheel phases through locomotion and verify the contribution by experiments using real bipedal robot.


Introduction
There are some studies of leg-wheel robots that are composed of the integration of legged robots and wheeled robots [1][2][3]. Legged robots have the advantage of versatile gait locomotion for various environments and wheeled robots have the advantage of stable and efficient locomotion. As the previous researches related to leg-wheel robots, the quadruped type [1] and the biped type [2,3] are proposed. The biped type leg-wheel robots enable the transformation between bipedal walking and wheel locomotion by utilizing the stewart platform mechanisms of legs [2] and having the additional actuators in knees as active driven wheels [3]. In this paper, we aim the establishment of the locomotion approach of legwheel robot that slides the foot by rolling passive wheel attached to the heel. Without the complex linkage mechanisms and the additional actuators, this locomotion approach is simply realized by only attaching the passive wheel to the heel of bipedal robot as shown in Fig. 1. There are several merits to utilize passive wheel for legwheel robot as well as the simplicity of attachment. Even if the general leg-wheel robot is unable to drive its wheels actively due to the accident of wheel drive system and the road surface conditions, the leg-wheel robot equipped with passive wheel has the possibility of locomotion by driving other joints of legs. Furthermore, integrating passive wheel into robot system enables the compliance locomotion for external environments because passive wheel joints relieve the internal force applied inside the closed link system consisting of legs and ground.
Having no swing leg is the key feature of the proposed locomotion approach of bipedal robot utilizing passive wheel. This feature contributes (1) expansion of support polygon during locomotion, (2) robustness and stopping for obstacles, and (3) prevention of body swing sideways.
By virtue of utilizing passive wheel and having no swing leg, a bipedal robot benefits from the features of locomotion as shown in Fig. 2.

Expansion of support polygon during locomotion
The locomotion rolling passive wheel on the ground expands support polygon of robot because of no swing leg as shown in Fig. 2-I. Figure 2-I is the bottom view of bipedal robot and the support polygon is generated by applying convex hull to the contact region of robot and ground. To maintain support polygon extensively implies  [4], and Capture Point [5] can exist. Capture Point is utilized as an indicator not to fall. To take these points into consideration is necessary to keep balance of robot and contributes the stable locomotion.

Robustness and stopping for obstacles
As shown in Fig. 2-II, we compare the case of having swing leg with the case of having no swing leg in the situation where robot collides with obstacle during locomotion. In the case that swing leg of robot collides with obstacle in the left of Fig. 2-II, the external force by impact is applied above the level of support region and robot tends to fall due to the large overturning moment around support region. By contrast, in the case of locomotion without swing leg by sliding passive wheel in the right of Fig. 2-II, the point of external force applied by impact can be always kept in the lower position. Therefore, even if the external force on the same level with the former case is applied, the overturning moment around support region is reduced smaller and robot acquires the robustness.
In the case of dynamic walking in the left of Fig. 2-II, robot is unable to stop immediately after the collision with obstacle is detected, and the complex processing for fall avoidance such as foot step modification [6] is required. On the contrary, the static locomotion approach in the right of Fig. 2-II has the feature that robot is able to stop with holding its pose immediately after the collision with obstacle occurs. As other features, to put swing leg always on the ground contributes the avoidance from the accident that robot steps on obstacle and falls over.

Prevention of body swing sideways
The walking of general bipedal robot has the problem that the body of robot swings sideways. This is caused by the sideways trajectory planning and generation of center of gravity projection and ZMP so that these points are included in support polygon as shown in the left of Fig. 2-I while the support polygon for right foot and left foot switches alternately. Fig. 2-III shows the situation that bipedal robot passes through narrow space between walls. As shown in the left of Fig. 2-III, the normal bipedal walking causes the body swing sideways and consequently, the body links of robot interfere with walls. In contrast, the center of gravity projection onto the expanded support polygon in the right of Fig. 2-I enables robot to prevent the lateral swing of body and achieve the passing through narrow space as shown in the right of Fig. 2-III.
Besides the above effectiveness, the prevention of body swing sideways contributes the fall prevention of the load mounted on robot and the prevention of the horizontal blur of robot vision.
Utilizing the slipperiness of passive wheel reduces the high friction and the high load such as the shuffle motion [7] of bipedal robot and enables the smaller shaking locomotion than the shuffle motion requiring the complex distribution of sole contact force and frictional force.

Stability margin maximization
In this section, we describe the optimization planning how to project the center of gravity of bipedal robot onto the support polygon expanded by utilizing passive wheel.

Definition
In general, when a center of gravity projection point of robot exists in a support polygon, robot is statically stable. Here, the center of gravity projection point is a foot of a perpendicular from the center of gravity of robot to the ground. Accordingly, when the center of gravity of robot is projected onto a certain point located within a support polygon, the planning that stability is guaranteed can be implemented. We consider the optimal position of a certain point located within a support polygon. This optimal position is the target of center of gravity projection. The indicator to optimize in this paper is the stability and the idea of stability margin [8] is applied as stability criterion and quantitative evaluation. Therefore, we define stability margin maximization to maximize the stability margin as objective function.
Stability margin is represented as the distance between center of gravity projection point and boundary line of support polygon. Namely, stability margin S(x, y) equals the minimum value among the distance d i (x, y) between center of gravity projection point (x, y) and each boundary line i of support polygon in the following formula.
For the static state, the center of gravity projection point (x, y) is in the set P composed of all points within support polygon. Let stability margin maximum point p max be the optimum point where stability margin S(x, y) is maximum, and p max is defined as follows.
Stability margin maximization is to provide the stability margin maximized by solving p max in Eq. (2).

Solution
To apply stability margin maximization to real robot, we need to solve the problem Eq. (2). When a support polygon has an inscribed circle, p max in Eq. (2) geometrically corresponds to the center of this inscribed circle. However, the case that the support polygon of real robot does not have an inscribed circle exists. Therefore, we formalize Eq. (2) into the generalized algorithm that is possible to be solved by computer. We generalize the support polygon of robot as shown in Fig. 2-I and show n-sided support polygon in Fig. 3. Let i be i = 0, 1, . . . , n − 1.
Equation (2) can be interpreted as the general optimization problem composed of objective function and constraint conditions. First, we assume that stability margin maximization is achieved with respect to a certain boundary line i. Let p i max be the stability margin maximum point for i. p i max is obtained as 2 variables x, y to maximize d i (x, y) (Objective Function). Meanwhile, from the definition Eq. (1) of stability margin, d i (x, y) should be smaller than any other distance d j (x, y), ∀j / ∈ i (Constraint Condition 1). Moreover, the inclusive condition (x, y) ∈ P for support polygon should be satisfied, which is provided as the precondition in Eq. (2). Based on Fig. 3, this condition is equivalent to the condition that u k × v k (x, y) ≥ 0 is satisfied for all boundary lines ∀k = 0, 1, . . . , i, . . . , n − 1 (Constraint Condition 2). Here, u k represents the vector from the start point to the end point of the boundary line k when turning clockwise round the boundary lines of support polygon, and v k (x, y) represents the vector from the start point of the boundary line k to an arbitrary point (x, y). Each component of the vectors follows the coordinate system in Fig. 3.
From the above, the optimization problem with respect to p i max is as follows.
p i max represents the center of the largest circle that contacts with the boundary line i and is included in the support polygon. d i (x, y) is denoted by u i and v i (x, y) from the geometric relationship in Fig. 3 and its absolute value can be removed from the Constraint Condition 2 in Eq.
(3). Besides, note that u i × v i (x, y) is the linear expression for x, y. d i (x, y) is as follows.
a i , b i , and c i are the constants particular to support polygon. By substituting Eq. (4) for Eq. (3), the optimization problem with respect to p i max is eventually summarized as follows.
(3) In the optimization problem Eq. (5), both objective function and constraint conditions are represented as the linear expression for 2 variables x, y, and therefore this problem can be solved by Linear Programming (LP).
In this paper, we utilize qpOASES [9] as the LP solver library. Finally, we obtain p max from p i max as follows.
represents the distance between p i max and the boundary line i and geometrically implies the radius of the largest circle that contacts with the boundary line i and is included in the support polygon.

Application to leg-wheel phases
We practically apply stability margin maximization to bipedal robot with passive wheel. We clarify the definition of support polygon for this leg-wheel robot.
As shown in Fig. 4, two types of support polygons are assumed for this robot. The support polygon P is generated by applying convex hull to the contact region including wheels and the support polygon Q is composed of the only feet soles. We compare P with Q in terms of the stability when the center of gravity (x, y) of robot moves posterior to the stance foot. In Fig. 4, four poses 1 � ∼ 4 � of robot are generated by inverse kinematics for the center of gravity (x, y). The relationship between (x, y) and P, Q is ∈ Q , and 4 : (x, y) ∈ bd(P), (x, y) / ∈ Q , respectively. Here, bd( * ) represents the boundary of set * . For the poses 1 � ∼ 3 � , the robot can keep the posture on the slippery tile floor, especially the robot has not fallen in 3 even though the center of gravity of robot locates posterior to the edge of the feet soles ( (x, y) / ∈ Q ). For the pose 4 , the robot falls backward because the center of gravity locates in the boundary of P ( (x, y) ∈ bd(P) ). Therefore, the support polygon that the center of gravity of robot can exist without falling is defined as P in this paper.
We define the phase that one leg of bipedal robot lands a sole and the other leg lands a passive wheel as leg-wheel phase. As opposed to the phase that both legs land soles, Fig. 4 Comparison of support polygon P including contact of wheels and support polygon Q composed of only feet soles when center of gravity (x, y) of robot moves posterior to stance foot the effectiveness of stability margin maximization for legwheel phase is remarkable because of the difference of right and left contact area. We show the results that stability margin maximization is applied to each leg-wheel phase in Fig. 5. The leg that lands a passive wheel in legwheel phase lifts up the foot around pitch by 20° keeping the height of the passive wheel axle as shown in the bottom right of Fig. 1. Keeping the height of axle constant guarantees the contact between wheel and ground during swing movement not to float and press the wheel against the ground. By applying stability margin maximization to the support polygons formed from sole and passive wheel on the ground, the stability margin maximum points are obtained as shown in Fig. 5a1-a4. Here, in Fig. 5a1-a4, the black coordinate systems are the world coordinate systems referenced in Fig. 3, the yellow regions are the support polygons, the red points are the stability margin maximum points, and the green circles are the stability margin circles whose radius equals stability margin S(x, y). The stability margin circles in Fig. 5a1-a4 are the circles that have the maximized stability margin. It is observed that the stability margin circles maximized by stability margin maximization represent the largest circles that are included in the support polygons and these circles contact with the several boundary lines of the support polygons. When these circles contact with all boundary lines of the support polygons, these circles correspond to inscribed circles.
By solving inverse kinematics for the center of gravity of robot to be projected onto the obtained stability margin maximum points, the reference joint angles of robot are determined. We show the results of dynamics simulation when these reference joint angles are commanded to bipedal robot on the simulator Choreonoid [10] in Fig. 5b1-b4. It is observed that robots keep balance without falling in all phases. Furthermore, following the sequential commands from Phase1 to Phase2 and from Phase3 to Phase4, robots are able to slide their passive wheels keeping balance. The swing leg trajectory that rolls passive wheel is determined as simple straight line trajectory because the height of passive wheel axle is kept constant during the sliding of foot. To guarantee the smoothness of the trajectory for the center of gravity of robot through the time, minimumjerk model [11] is adopted as the trajectory algorithm. The minimum-jerk model enables the motion trajectory of robot to be interpolated smoothly without loss of balance for sudden change of stability margin maximum point p max due to the change of support polygon. Figure 6 represents the single leg stance on Phase3 where the center of gravity of robot is projected onto the center of sole without stability margin maximization. Fig. 6a shows the center of gravity projection point (the red point) and the stability margin circle (the green circle). Considering the radius of stability margin circle (the stability margin) S, the lateral amplitude |y|, and the normal force distribution applied to right foot and left foot F R z : F L z , we make a comparison with the result of Fig. 5a3. We show the comparison result in Table 1. The larger the number of S is, the higher the stability is. This effectiveness is observed in "Stability margin maximization (Fig. 5a3)" in Table 1. The smaller the number of |y| is, the more effective the prevention of body swing sideways is. Similarly, this prevention is observed in "Stability margin maximization (Fig. 5a3)" in Table 1. This result regarding the lateral amplitude is also observed from the difference of simulation results between Fig. 5b3 and Fig. 6b. Finally, from the result of F R z : F L z , it is observed that stability margin maximization distributes the robot's own weight to the foot whose wheel is grounded.

Stability margin maximization with lateral constraint
Although stability margin maximization enables the prevention of body swing sideways, the further prevention is required so that robot can pass through narrower space. In this paper, we propose stability margin maximization with lateral constraint by adding the constraint with respect to y to the constraint conditions of stability margin maximization. This strategy is feasible based on the addition of the following conditional expression to the constraint conditions in Eq. (5).
ε represents the parameter regarding lateral tolerance limit.
For example Phase3 of leg-wheel phases, we show the comparison result between (A) Without lateral  constraint and (B) With lateral constraint Eq. (7) in Fig. 7. In Fig. 7, the walls of narrow space are highlighted in red and we set the distance between walls as 750 mm. When bipedal robot passes through this space, normal stability margin maximization causes the collision between robot links and wall as shown in Fig. 7a. By contrast, stability margin maximization with lateral constraint of ε = 35 mm enables robot to pass through this space without collision based on the smaller swing sideways than (A). We show the result of S, |y|, and F R z : F L z for Fig. 7b in Table 1. In comparison with the result for normal stability margin maximization (Fig. 5a3), the stability margin S decreases, whereas the lateral amplitude |y| gets smaller and is given as the value of ε . From this observation, it is found out that stability margin maximization with lateral constraint is the trade-off strategy that prevents the body swing sideways more effectively by reducing the stability margin. From the result of F R z : F L z , it is also observed that the ratio of weight load distribution to the foot whose wheel is grounded increases besides the prevention of body swing sideways.

Fall prevention functions
In this section, we describe the online-based prevention functions against the fall factors that are difficult to be prevented by only offline planned stability margin maximization.

Emergency Stopper
Because bipedal robots such as humanoid robots tend to lose balance and fall, the methods to prevent the fall accident are proposed [12][13][14]. Emergency Stop [12] is one of those methods. Emergency Stop is utilized for fall detection and motion reverting not to fall based on the measurement of Capture Point (CP) [5] ξ CP x,y defined as the following formula.
ξ A x,y represents the vector of (x, y) elements of point A and ξ A z represents the value of z element of point A. Regarding the point A, CP represents Capture Point, COG represents center of gravity of robot, and ZMP represents Fig. 7 Comparison at narrow space on Phase3

Fig. 8 Support polygon and safe region
Zero Moment Point, respectively. g represents gravitational acceleration.
In this paper, we define the region with margin from support polygon as Safe Region as shown in Fig. 8 and fall detection works in the moment CP gets out of the Safe Region. And we define the function to stop active motion immediately after fall detection is detected as Emergency Stopper. Safe Region in Fig.8 is obtained by applying convex hull to the region reduced scale α ( 0 < α ≤ 1 ) for the grounded region of foot and wheel.

Fall prevention function for stumbling
We describe the problem when bipedal robot stumbles over obstacle during locomotion and propose the strategy to solve this problem. Fig. 9a shows the snapshot when the right foot of real bipedal robot stumbles over the concrete block on the road during the locomotion from Phase3 to Phase4 in Fig. 5. The real robot falls immediately after this stumbling. Figure 9b shows the state of CP and Safe Region at this time. The real state of robot measured from joints encoder and attitude sensor is reflected in Fig. 9b and we set α = 0.9 . While the real robot falls, it is observed that CP can not break out of the Safe Region and Emergency Stopper does not work. According to Eq. (8), it is inferred that the impact of stumbling is not transmitted to the change of center of gravity.

Force translation of Safe Region
As described in Problem, because Emergency Stopper is valid only when CP breaks out of Safe Region, the novel strategy is required for the prevention of fall due to stumbling. Therefore, besed on the inference in Problem, we assume that the virtual force f is applied to center of gravity, which f is the stumbling force applied to foot of robot and can be measured by attached force sensor. And we propose force translation of Safe Region based on the virtual change of center of gravity as shown in Fig. 10. Figure 10 represents the model that robot is approximated by inverted pendulum. Force translation of Safe Region is the technique to translate Safe Region by the change of center of gravity x SR that can keep the balance when the horizontal virtual force f is applied to the center of mass

Fall prevention function for leg-wheel transition
To achieve the locomotion without swing leg, the legwheel transition from two legs to one leg and one wheel is essential. We describe the fall problem in leg-wheel transition and propose two strategies to solve this problem. Figure 11 shows the block diagram of the legwheel transition and the approach to solve the problem that real robot falls in spite of the unchanged target COG (center of gravity). The leg-wheel transition is executed between the stance phase and the phase to which stability margin maximization with lateral constraint ε = 0 is applied. The foot sole of robot is rotated around pitch by 20° with passive wheel on the ground as shown in the upper left of Fig. 11. Note that the robot keeps the height of the passive wheel axle during the leg-wheel transition not to float and press the wheel against the floor. In the stance phase, the target COG is determined in advance from the result of stability margin maximum point p max calculated in the phase with ε = 0 . Therefore, the target COG is not changed in the leg-wheel transition between two phases. Even though the target COG in the robot model is not moved, the actual COG of the real robot can be moved and the real robot falls in the real environment. There are several causes of fall such as modeling error as described below. This fall problem of real robot can be solved by the detection of CP that gets out of Safe Region because the velocity of falling affects CP as shown in Eq. (8).
We propose two strategies 1 Online search of lateral constraint and 2 Longitudinal shift of stability margin maximum point to correct the stability margin maximum point based on the detection of CP so that real robot does not fall.  Fig. 12b shows the dynamics simulation result, Fig. 12c shows the real experimental result, and Fig. 12d shows the state of CP and Safe Region. While the robot succeeds in transition without falling in (B) Simulation, the robot falls in (C) Experiment. From the data (D) in this moment, fall detection is confirmed because CP gets out of the Safe Region. After Emergency Stopper halts the transition motion based on fall detection, the execution of motion reverting [12] enables the fall prevention. However, because leg-wheel transition is not accomplished by only motion reverting, additional strategies are required. In particular, online-based behavior strategies are desired against the reduction of stability margin S due to the too strong constraint such as ε = 0 , the deflection of joints of real robot, the error of end effector, and the modeling error regarding contact area of ground and wheel.

Online search of lateral constraint
This is the online-based technique to search the minimum limit value ε in lateral constraint Eq. (7) required for achievement of leg-wheel transition. Let ε 0 = 0 be initial value of ε , and ε k , which is k-th search result for ε , is obtained sequentially based on the update rule Eq. (10). Note that k is updated to k + 1 every time CP gets out of Safe Region.
Let η > 0 , T Start k , and T Fall k be the learning rate, the start time of leg-wheel transition for ε = ε k , and the time CP gets out of Safe Region, respectively. T Fall k − T Start k is equivalent to the elapsed time (s) from the start of transition to the moment robot is likely to fall. The larger ε is, the weaker lateral constraint is. Eventually, ε converges as T Fall k − T Start k → ∞ , which denotes the state robot does not fall. As a device to hasten the convergence, the scale factor α k of Safe Region is changed adaptively for the value of ε k by the following rule.
Let α 0 and α max be lower limit and upper limit of α , respectively. In this paper, we set α 0 = 0.6, α max = 1.0 . ε max represents the value of |y| of stability margin maximum point without lateral constraint. According to Eq. (11), Safe Region is expanded as search for ε k proceeds, and T Fall k − T Start k also grows larger.

Longitudinal shift of stability margin maximum point
As opposed to the preceding paragraph, this is the onlinebased technique to achieve leg-wheel transition without changing the constraint ε = 0 and without lateral movement. In this technique, without shifting y element of stability margin maximum point (x, y), x element regarding longitudinal direction is shifted by the following shift amount x shift .

Locomotion experiment based on stability margin maximization
We show the locomotion experiment by applying stability margin maximization with lateral constraint to a real bipedal robot in Fig. 13. The weight of bipedal robot is 63.9 kg. The feet of robot are equipped with passive wheels without actuators as shown in Fig. 1. The radius of wheel is 44 mm, the width of wheel is 90 mm, and the surface of wheel is made of rubber. The radius of wheel is chosen as small as possible so that the wheel can be stored inside the foot and the robot is safer against slipping, whereas it needs to be chosen not too small so that the contact area between wheel and ground is larger enough to secure support polygon. Because we aim the minimum configuration that wheels can be attached simply to the feet of bipedal robot, the number of passive wheels is defined as one wheel for one foot. And we determine the width of wheel to be the same as the width of foot so that the bipedal robot does not become unstable around roll. The upper part of Fig. 13 shows each phase for locomotion planning of center of gravity projection marked in red point on the bottom view of robot. The lower part of Fig. 13 shows the snapshots of real experiment on the flat carpet floor and the slippery tile floor for the corresponding phases. In the lower part of Fig. 13, timestamps are shown when we set the start time on Phase1 as 0.0 s. In this experiment, we set lateral constraint as ε = 55 . One cycle of this locomotion is composed of 18 phases and stability margin maximization with lateral constraint is applied on Phase3, 4,7,8,11,12,15,16. On the phases that both feet of robot are grounded before and after these leg-wheel phases, the transition is executed without changing center of gravity projection point such as 2 → 3 and 4 → 5 . The achievement of locomotion based on the transition of leg-wheel phases is observed in the snapshots of real experiment. The passive wheels of real bipedal robot slide 0.36 m in 1.8 s for 3 → 4 and 15 → 16 , and slide 0.72 m in 3.6 s for 7 → 8 and 11 → 12.
To consider the efficiency of locomotion, we compare the electric power and energy consumed by the locomotion based on stability margin maximization with utilizing passive wheel with by the locomotion based on single leg stance as shown in Fig. 6 without utilizing passive wheel. Figure 14 shows the change of electric power consumed by the power supply for robot servo during one cycle locomotion from Phase1 to Phase18. Figure 14a shows the power consumption P Single for the locomotion by switching the right and left single stance phases without utilizing passive wheel. Figure 14b shows the power consumption P SMM for the locomotion by applying stability margin maximization with lateral constraint of Fig. 13 Locomotion based on stability margin maximization with lateral constraint ε = 55 ε = 55 with utilizing passive wheel. These are compared under the same conditions of the locomotion distance of robot and the transition time among each phase. On Phase5-Phase6, Phase9-Phase10, Phase13-Phase14, and Phase17-Phase18, the unnecessary power consumption is observed in (A) P Single , which requires the large movement of center of gravity to switch single stance phases. By contrast, in (B) P SMM based on stability margin maximization, the cooperative use of support regions of foot and wheel contributes to the achievement of efficient locomotion without the unnecessary movement of center of gravity. On Phase3-Phase4, Phase7-Phase8, Phase11-Phase12, and Phase15-Phase16 that are the phases to slide wheels, it is observed that P SMM is smaller than P Single . It is inferred that this effect reducing power consumption is brought by the weight load distribution to grounded wheel based on stability margin maximization as shown in the result of F R z : F L z in Table 1.
To evaluate these locomotion, we compare the energy consumption E {Single,SMM} during locomotion under the conditions of the same distance and the same time as calculated by the following formula.
T cycle represents the locomotion time for one cycle from Phase1 to Phase18. According to the comparison conditions, we set T cycle = 43 s for both P Single and P SMM . By calculating the value of E {Single,SMM} in Eq. (13) from the results of power consumption P {Single,SMM} in Fig. 14, we obtain the results of E Single = 5499.49 J and E SMM = 4093.29 J. These results denote E Single > E SMM and stability margin maximization succeeds in reducing energy consumption during locomotion by approximately 1400 J.
Leg-wheel robots have the road ability of legged robots to travel over various environments and the stable and efficient mobility of wheeled robots. The leg-wheel robot to utilize passive wheels based on the proposed method can achieve these ability and mobility by locomotion with pushing obstacles in the environment where disturbance is mixed. To clarify whether this leg-wheel robots rather than general bipedal robots are useful for specific task, we compare (A) General bipedal walking with (B) Proposed locomotion approach for the locomotion tasks with pushing obstacles out in the debris environment as shown in Fig. 15. Figure 15 shows the dynamics simulation results that the bipedal robot and the leg-wheel robot try to pass through the debris environment where 5 kg of cuboid obstacle, 5 kg of cylinder obstacle, and 3 kg of cuboid obstacle are mixed. In Fig.15a, the robot starts bipedal walking with swing leg and tries to push the obstacles by walking. However, the robot loses balance while pushing three heavy obstacles by swing leg in Fig. 15a Fig. 15a-6. By contrast, in Fig. 15b, the robot starts proposed locomotion utilizing passive wheel without swing leg and succeeds in pushing the debris obstacles out. The leg-wheel robot is able to turn by bipedal walking in Fig. 15b-5 and restart this locomotion toward the different direction in Fig. 15b-6. These comparison results demonstrate the robust locomotion skill in a debris environment for the proposed locomotion approach utilizing passive wheel. The locomotion ability of this leg-wheel robot to push obstacles away is demonstrated in real environment. In Fig. 16, 3.4 kg, 3.8 kg, 5.5 kg, and 8.3 kg of the mixed obstacles are placed randomly on the road. It is observed that the real leg-wheel robot moves while pushing obstacles away. Stability margin maximization contributes to the response to disturbance during locomotion and the improvement of robustness.

Fall prevention experiment for stumbling
Although the robustness is secured by stability margin maximization, robot stumbles in the case that the larger disturbance exists. We show the fall prevention experiment for this case in Fig. 17. In Fig. 17-1, 13.9 kg * 2 pieces of the concrete blocks are placed side by side on the road and the robot can not push these obstacles away due to the heavy weight. For this situation, Fig. 17-2 shows the state that force translation of Safe Region works. In Fig. 17-2, the Safe Region is translated forward by x SR = 243.2 mm and Emergency Stopper works since CP gets out of the Safe Region. After Emergency Stopper works, the robot slides the passive wheel of stumbling foot backward and the emergency stopping to prevent fall is achieved so that CP gets inside the Safe Region as shown in Fig. 17-3.

Online search experiment of lateral constraint
We set the initial value of lateral constraint as ε 0 = 0 and show the result of online search of the value ε so that the real robot can achieve leg-wheel transition on the carpet floor without falling in Fig. 18. Figure 18 shows the results of stability margin maximization for ε = ε k and the moments CP gets out of the Safe Region, where CP is inside the Safe Region when the online search is converged. The elapsed time T k = T Fall k − T Start k from the start of leg-wheel transition to these moments is also shown in Fig. 18. Note that the leg-wheel transition is expected to be executed as Phase2 → 3 in Fig. 13 for 5 s. Regarding the way to determine η in the update rule Eq. (10), we can set η = 138.542 from the relationship of Eq. (10) so that ε k+1 − ε k , which is the change of ε , equals 1 % of the maximum range ε max − ε 0 = 93.3491 when the sufficient time T k = 5 s for transition elapses. From the search experiment under the above conditions, we obtain the result that ε eventually converges with ε 5 = 54.0728 after 6 trials. This convergence value is close to the setting value ε = 55 applied in Fig. 13 and the real bipedal robot is able to achieve locomotion as a result.
We also show the experiment of fall prevention for leg-wheel transition and then moving on the slippery tile floor in Fig. 19. In Fig. 19(1)-(4), the online search of lateral constraint ε is executed on the tile floor and ε converges with ε 3 = 23.095 after only 4 trials. In Fig.19(5)- (8), stability margin maximization with the obtained lateral constraint ε = 23.095 enables the robot to succeed in moving so that CP does not get out of the Safe Region. It is inferred that the smaller convergence value of ε and fewer trials than the previous experimental result shown in Fig. 18 are caused by the floor conditions (the difference between elastic carpet floor and rigid tile floor). The leg-wheel transition is achieved without

Longitudinal shift experiment of stability margin maximum point
In order to achieve leg-wheel transition without falling, besides the approach to change lateral constraint, it is effective to shift the stability margin maximum point which center of gravity of robot is projected onto in the longitudinal direction as shown in Fig. 20. Figure 20 shows the plot regarding the change of longitudinal position ξ COG x of center of gravity. At time 0 s, the legwheel transition with ε = 0 starts. At time T 1 = 2.24 s, Emergency Stop works since CP gets out of the Safe Region. From time T 1 = 2.24 s to time T 2 = 6.08 s, the motion reverting to previous phase is executed. At time T 3 = 12.79 s, the longitudinal shift of stability margin maximum point is achieved so that CP gets inside the Safe Region. The longitudinal shift amount x shift = −42.47 mm is obtained during the period from time 0 s to time T 1 = 2.24 s based on Eq. (12). After the stability margin maximum point is shifted by x shift at time T 2 = 6.08 s when motion is reverted to previous phase, the leg-wheel transition with ε = 0 is executed without falling during the period from time T 2 = 6.08 s to time T 3 = 12.79 s.
Locomotion with ε = 0 is also achieved by utilizing the obtained longitudinal shift amount x shift . Fig. 21 shows the plots regarding the change of lateral position of center of gravity during the locomotion applied by stability margin maximization without lateral constraint ("no constraint"), stability margin maximization with lateral constraint of ε = 55 ("with ε = 55"), and longitudinal shift of stability margin maximum point with ε = 0 ("with ε = 0"), respectively. It is observed that the effectiveness for prevention of body swing sideways is heightened in this order. Utilizing the locomotion with ε = 0 , the real bipedal robot is able to achieve the passing through narrow space as shown in Fig. 22. In Fig. 22a applied by stability margin maximization with lateral constraint of ε = 55 , the interference between the real robot and the wall of narrow space is observed. By contrast, in Fig. 22b during the locomotion with ε = 0 , minimizing the lateral amplitude enables the real robot to pass through narrow space without the collision with walls.

Conclusion
This paper established the locomotion approach without swing leg by utilizing passive wheels attached to the feet of bipedal robot without wheel actuators. The proposed locomotion approach provided the effectiveness of the stability based on expansion of support polygon during locomotion, the robustness for external obstacles and emergency stopping not to fall, and the adaptability to prevent body swing sideways. To achieve these features, this paper proposed stability margin maximization and fall prevention functions. Stability margin maximization is the optimization approach how to project the center of gravity of robot onto the support polygon formed from sole and wheel. Fall prevention functions are the strategies against the fall factors of real robot in the real environment that is a difficult situation to prevent fall by the only offline planning. Finally, this paper demonstrated the experiments using real bipedal robot to verify these effects. In particular, the experiments in this paper proved the remarkable outcome that stability margin maximization in the proposed approach can reduce power and energy consumption during locomotion. This lower power and energy consumption for the proposed locomotion is revealed from the experimental results as shown in Fig. 14 and it contributes to the improvement of locomotion efficiency of bipedal robot. This paper expanded locomotion style of leg-wheel robots without introducing complex mechanisms to legged robots, adding actuators, and modifying configuration of electrical system. On the other hand, there are some disadvantages as trade-offs such as (1) slow locomotion, (2) risk of backward falling due to passive wheels on the heel of robot, and (3) limits of robustness for floor conditions. Regarding (1), improvement of locomotion speed is expected and it is required to expand from static locomotion strategies to dynamic locomotion strategies. This can be feasible by replacing center of gravity trajectories with ZMP trajectories, which are planned by stability margin maximization proposed in this paper. Regarding (2), this paper showed the evidence that robot does not fall backward even when center of gravity of robot locates posterior to edge of feet soles in Fig. 4. However, the more terrible modeling error between robot model and real robot and the larger disturbance can affect the risk of backward falling. We consider that the switching between active wheels and passive wheels is required to brake robot and prevent falling against these external factors. Regarding (3), this paper demonstrated the locomotion experiment in Fig. 13 and the