(Microsoft Word - 93-106\335\321\307\323 \346\343\321\312\326\354) Al-Khwarizmi Engineering Journal Al-Khwarizmi Engineering Journal, Vol. 14, No.2, June, (2018) P.P. 93- 106 Creeping Gait Analysis and Simulation of a Quadruped Robot Firas A. Raheem* Murtadha Khudhair Flayyih** *,** Department of Control and System Engineering/ University of Technology *Email: 60124@uotechnology.edu.iq **Email: cse.60354@uotechnology.edu.iq (Received 22 May 2017; accepted 5 December 2017) https://doi.org/10.22153/kej.2018.12.004 Abstract A quadruped (four-legged) robot locomotion has the potential ability for using in different applications such as walking over soft and rough terrains and to grantee the mobility and flexibility. In general, quadruped robots have three main periodic gaits: creeping gait, running gait and galloping gait. The main problem of the quadruped robot during walking is the needing to be statically stable for slow gaits such as creeping gait. The statically stable walking as a condition depends on the stability margins that calculated particularly for this gait. In this paper, the creeping gait sequence analysis of each leg step during the swing and fixed phases has been carried out. The calculation of the minimum stability margins depends upon the forward and inverse kinematic models for each 3-DOF leg and depends on vertical geometrical projection during walking. Simulation and results verify the stability insurance after calculation the minimum margins which indicate clearly the robot COG (Center of Gravity) inside the supporting polygon resulted from the leg-tips. Keywords: Quadruped robot Creeping Gait, Stability Margins, Quadruped robot Simulation. 1. Introduction Leg-based walking systems utilization has turned out to be very familiar with the robotics field. These systems present better adaptability over wheeled systems, particularly when dealing with irregular territories [1]. The ability of legged robots in dealing with rough territories and obstacles which also have the planning technique of the standard walking gaits for a quadruped robot [2, 3]. Among legged robots, quadrupeds have the preferred of having a less complex structure while fulfilling the requirements for usage of a statically stable gait during walking, resulted from the stable support polygons. A lot of researches about the quadruped locomotion frequently have its origins in the investigation of natural biological walking gaits. These studies have prompted to the improvement of both statically stable walking gait and dynamic running gaits for legged robots [4]. The main point is the need to design and implement a robust and statically stable gait during the quadruped robot walking. The property of their gaits would be repeating the ability of step sequence and stability over irregular territories [5]. Step sequences for quadrupeds, and the criteria for stability verification on an inclined terrain, as well as the talking about moving between varieties in gait and omnidirectional statically stable walking [6]. Stability criteria have some elements such as COG (Centre of Gravity), Support Polygon and Stability Margin which is also talked about the gait transitions performed by the creeping gaits [7, 8]. The generation and the sequence of such leg motion in the quadruped robot is called gait [9]. There are many motions that quadruped robots doing, such as stair climbing with many performances [17], flying trotting [19], jumping [20], and walking/trotting over obstacles [18]. Generally, the gait of quadruped robots can be separated for two groups: the first one is the statically stable gait and the second one is the Firas A. Raheem Al-Khwarizmi Engineering Journal, Vol. 14, No. 2, P.P. 93- 106 (2018) 94 Dynamic running gaits, depending on the gait type [10]. Quadruped body must be stable and have the ability during any transition gait for moving from one position to another. The statically stable gait is simpler and slower contrasted with dynamic gaits. It indicates that the COG projection of quadruped robot is on the ground at instance time is inside support polygon during the walking. Among these properties, the quadruped robot meets the necessity of both force and moment adjustment during the walking gait and here the quadruped robot is said to be statically stable. Dynamic gaits such as the trotting gait [11] aim to the reality that if the COG projection of the quadruped robot on the ground is not inside the support polygon framed by the leg tips. In this case, the quadruped robot is in the unbalanced period and has a tendency of dumping. Then the quadruped robot in case of the dynamic running gaits needs to satisfy the dynamically stable requirement rather than static stability analysis. Most of the quadruped robots can be considered as statically stable, when there are three legs at least should be stayed in contact with the ground at the same instant time while another leg is swinging in the air. The most biological statically stable gait that matches the animals or insects is the creeping gait. In this paper, the problem of static stability analysis is achieved. This is sufficient for ensuring safety walking. During the movement of the quadruped robot, the motion will be stable if and only if the body center of gravity COG is either inside the supporting triangle in case of one leg in a transition phase or inside the supporting polygon in case all legs are on the ground. At the time period of the quadruped robot is walking the problem is how to produce and control sequence of putting and lifting all the legs extend at any time moment. The quadruped robot walks with stable steps if the stability margins have positive values, another case the robot will be unstable and the walking gait will fail. 2. Quadruped Robot Description To concentrate on the general guidelines of quadruped robot motion gaits, though, different models with some acknowledgment attributes is embraced. These attributes have three criteria [13] 1. The standard quadruped walking at constant speeds, at steady state movement, over an even supporting surface. In this case, the legs of the quadruped robot are repeated periodically during the locomotion. 2. The standard quadruped robots are symmetrical in both longitudinal and lateral directions. This property leads to the facts that the design is adapted to have a similar structure for every leg with less the mechanical complexity of the system and then less in cost . 3. In case of a regular shape quadruped robot, COG located at the robot geometrical center of the quadruped robot body, which is known as COG the center of gravity. The quadruped robot is shown in Figure (1) demonstrates the typical model of a quadruped robot which represented the above three criteria. In this research, the motion task assumption includes that the quadruped robot walks towards the X-axis only. This means that there is no change in the Y-axis stride. So the vertical projection of the stride for all legs is a line. These stride lines are situated in the longitudinal direction to fulfill the condition of forwarding walk. Fig. 1. Quadruped robot (top view) illustrate the length (2b) and the width (2c) and the COG point From Figure (1), the distance for all legs from the closest point of the stride line to the COG projection is (a), the distance for all legs from the farthest point of the stride line to the COG projection is (b)and (c) is the distance between two front stride line or two back stride line. The definition of stride line for every leg in X-axis is [13] Stride line (λ) = b – a …(1) When the quadruped robot body is moved the COG is moved, correspondingly followed by the Firas A. Raheem Al-Khwarizmi Engineering Journal, Vol. 14, No. 2, P.P. 93- 106 (2018) 95 four legs moving. According to the gait sequence of quadruped robot creeping gait, there is no necessity to move the robot body and legs at the same time [14]. The forward walking is towards the X-axis. In the first origin support pattern, the positions of the four legs at (x1, c) (x2,-c) (x3, c) (x4,-c), and the length of stride is λ. By this manner, the positions of every leg after one step movement are located at (x1+ λ, c) (x2+ λ, -c) (x3+ λ, c) (x4+ λ,-c). These locations should be matching the stride constraint as follows: � � �� � � � � �� � � �� � � � �� �� � � � �� � � �� � λ � � � � �� � λ � � �� � � � λ � �� �� � � � λ � �� Along these considerations, the largest stride of the reachable area for a quadruped robot can be calculated from equation (1). Most of the quadruped robots have the value of length equals to 2 b and the value of width equals to 2 c. henceforth, the ratio between the length and width is equal to c/b. In this paper, each leg contains three joints 3-DOF called (Coxa-joint, Femur-joint, and Tibia- joint). So the total number of DOFs is equal to 12. The overall the quadruped robot is showing in Figure (2). Where in Figure (2-a) shows the main part in the quadruped robot is represent a mimic of the natural anatomy of legged animals. From Figure (2-a), Leg1 is the right front (RF) side and Leg3 is the right rear (RR) side and Leg2 is the left front (LF) side and Leg4 is the left rear (LR) side. In this Figure, it is illustrated and guaranteed to be inside the supporting polygon. The forward direction of motion of the quadruped robot, in case of the forward direction of motion, is set to be with the X-axis direction and the lateral motion of the quadruped robot will be in Y-axis. The quadruped body has a symmetrical dimension in X and Y axes. This gave the robot better stability during the walking and making it statically stable. In Figure (2-b) shows each leg has three links and joints (Coxa, Femur, and Tibia). (a) (b) Fig. 2. Quadruped Robot Simulation Using Matlab (a) quadruped robot with four legs labeling (b) showing the three joints (Coxa, Femur, Tibia) in each leg. To grantee that the quadruped robot is at stable situation during its movement, the position of each leg plays a very important rule of the stability calculations. In order to ensure the static stability condition of the quadruped robot, the forward kinematics and inverse kinematics is needed priory before analysis. 3. Creeping Gait Analysis and Sequence Description The quadruped robot gaits can be classified into three types of main gaits depending also on the value of the duty factor ( �), where i=1,..4, like: crawl gait, running gait, and the galloping gait. The duty factor ( �) is the ratio between the time period of one leg spend in the air to the time period spend on the ground, according to the creeping gait sequence, creeping gait has duty factor equal to 0.75. The running gait has duty factor ranging between 0.5-0.75, and the galloping Firas A. Raheem Al-Khwarizmi Engineering Journal, Vol. 14, No. 2, P.P. 93- 106 (2018) 96 gait has a duty factor less than 0.5. These gaits that are utilized by quadruped animals have regularly been given names, which generally known as creeping gait, trotting gait and bounding gait [١6]. Some of these gaits are utilized by some mammals. For a very slow walking, animals such as elephants are using the creeping gait. The advantage of this gait, that it can be generally statically stable movement and it can be used at the range of low-speed motions. The creeping gait needs at least three legs on the ground at the translation phase of the other leg in the air. This condition is important for achieve a statically stable gait. The creeping gait used the mechanism that the back leg is touching-down the ground followed by the front leg in the same side leaving the ground. After the front leg is touching-down the ground, the back leg in the other side of the body will begin to leave the ground, and so on, this sequence is repeated continuously during the locomotion of the quadruped robot. There are six types of legs sequences for quadruped robot. In this paper, the leg swinging sequence is (RR→RF→ LR→LF) where (R: right, L: left, R: rear, F: front). Where RR is leg4, RF is leg2, LR is leg3, and LF is leg1. The main feature of this specific sequence is to have a safe walking that ensures the body robot of moving forward at the instance time [12]. There are two Cartesian coordinate reference frames in the quadruped robot locomotion: the first one is the ground frame, the second one is the coxa joint frame related to the robot body. In these coordinate frames, the front side of the body is set towards the X-axis, while the body left side direction is set towards the Y-axis. In this way, these two frames are important for the computations that ensure the COG location is inside the robot body projection during leg placing and lifting. [13] 4. Kinematics Modeling of a Quadruped Robot 4.1 Forward Kinematics Depends on the leg design, a quadruped robot relies on the configuration of each leg. Since all parts of walking are represented by the physical constraints of the leg. Each leg has three-revolute joints named in the kinematical chain (r1, r2, r3). These joints have been chosen for every leg mechanism. The basic idea is to model the robot leg structure as the biological leg structure for animals or insects. The forward kinematics is applied to drive the geometrical model for every leg related to the center of the robot body, and then find the position and orientation of the end- tip (xi, yi, zi) where i=1,2,3, and 4, for each leg, the leg structure appeared in Figure (3) Fig. 3. Coordinate frames for one leg of a quadruped robot. [9] The D-H (Denavit-Hartenberg) parameters of the leg are shown in Table (1). The leg's frame of this robot begins with the link zero which is the point on the robot body. The leg is jointed as link one is the coxa link, link two is the femur link and link three is the tibia link. Legs are set to be symmetrically around the coordinate axis in the direction of movement. In our case, the direction of motion is the X-axis. The general frame for the transformation matrix from a link i to a link (i-1) utilizing the D-H parameters is given in Equation (2.1) [9]: The definition of the transformation matrix is a multiple series of transformation along: 1. Translate d� along ���� axis. 2. Rotate �� about ���� axis. 3. Translate �� along ���� axis. 4. Rotate �� about ���� axis. So, for one link and joint the translations and rotations are illustrated in the transformation matrix as [21]: ����� = ������ ����������� ���������� ������������ ���������� ����������� �������00 �����0 ����� ��0 1 …(2) The link parameter table of one leg in the quadruped robot is Firas A. Raheem Al-Khwarizmi Engineering Journal, Vol. 14, No. 2, P.P. 93- 106 (2018) 97 Table 1, The Denavit-Hartenberg parameters table for one leg in our quadruped robot Link No. Link name !" (deg) #" (cm) $" (cm) %" (deg) 1 Coxa 90 �� �� �� 2 Femur 180 �� 0 �� 3 Tibia 0 � 0 � Where �� = (9 cm), ��=2.5 cm ��= 5 cm, and � = 9 cm. �� is the distance from Coxa joint to the ground. �� are the lengths of the leg links. Now, obtaining the overall transformation matrixes by making a product of three transformation matrixes as the following: �&'()*)+, = �&'()-,./0 ∗ �-,./02�*�) …(3) From Equation (2), the transformation matrixes for one quadruped’s leg are describing as below: ��3 = � �� 0 �� ������ 0 ��� ����00 10 0 ��0 1 …(4) ��� = � �� �� 0 ������ ��� 0 ����00 00 �1 0 0 1 ...(5) � � = � � �� 0 � � � � 0 � � 00 00 1 0 0 1 …(6) By multiplication the Equations (4) and (5) 4��3 ∗���5 yields: ��� 0 �� ������ 0 ��� ����00 10 0 ��0 1 � �� �� 0 ������ ��� 0 ����00 00 �1 0 0 1 = ����� ���� ��� ��4����6��5���� ���� �� ��4����6��5��0 ���0 0 ���� � ��0 1 …(7) In the same method and by multiplication Equations (7) by (6) 4��3 ∗ ���5 ∗ � � yielding equation (8): ���4��� 5 ����� ��� � ��4��� 5 � ��4����6��5��4��� 5 ��4��� 5 �� � ��4��� 5 � ��4����6��5��� 0 �4��� 50 0 � 4��� 5 � �2�2 � ��0 1 …(8) Where: �� = ����� and �� = ����� , for i=1,2, and 3. 4��� 5= cos 4�� � � 5 4��� 5= sin 4�� � � 5 From the Figure (3) and by using the equation (8), the coordinate position of each leg-tip can be calculated as the following: �� = �����4� 4������ 5 � 4�������6��55 (9) =� = �����4� 4������ 5 � 4�������6��55 (10) �� = � 4������ 5 � ������� � �� (11) Where: (��,=�,��) is the coordinates of the leg-tip, i=1…4. 4.2 Inverse kinematics The inverse kinematics method formulated the joint angles from a given position and orientation of the legs tip. These position and orientation are achieved from the forward kinematics. The objective is to find the leg three joint angles ��,�� and � , relating to the required of the leg tip position [15] as shown in Figure (4). Fig. 4. Front view showing one leg of the quadruped and its links (Coxa, Femur, Tibia) [15] To find the inverse kinematics in each leg joint, the method simplifies this problem from 3D problem view into two 2D problem view, after that, using the geometrical analysis to solve and find the solution for (��,�� and � ), as showing in Figure (5) Fig. 5. Quadruped top view showing %> From Figure above, �� can be calculated, as the following: Firas A. Raheem Al-Khwarizmi Engineering Journal, Vol. 14, No. 2, P.P. 93- 106 (2018) 98 (? = tan 4��5 Thus, �� = B����4(?5 …(12) To find the other two angels which are �� and � . These angles are in the same (Y-Z) plane. To find ��, divide this angle into ��and �� to simplify the problem. To find �� firstly working out on the distance L to get the angle �� as showing in the following Figure (6). Fig. 6. Quadruped robot one leg side view showing links (Coxa, Femur, Tibia) and %C and %D [15] From the Figure above calculate the followings: E = F�'--+,2� � 4E� � �5�G …(13) Where c is the Coxa link length, and: �� = �����4HIJJKLMN 5 …(14) From the Cosine rules, �� and � angles can be calculated. From Figure (6) shown above the three sides of a triangle are (Femur, Tibia, and L). To calculate �� : �� = O� � E� � 42 ∗ O ∗ E5cos 4��5 Thus: �� = ����� PQG6NG�RG�∗Q∗N S …(15) Where: F is the length of Femur link. T is the length of Tibia Link. Although, ��is equal to: �� = �� � �� �� = ����� PHIJJKLMN S � ����� PQG6NG�RG�∗Q∗N S …(16) Finally, � is calculated as: � = ����� PQG6RG�NG�∗Q∗R S …(17) 5. Stability Analysis for Quadruped Robot Creeping Gait 5.1 Static Stability Criterion With a specific goal for the quadruped robot to walk with statically stable movements, the COG projection on the ground should be inside the supporting polygon framed by the leg tips. This condition is considered as the main criteria for the statically stable walking. When this condition is completely satisfied, the quadruped robot won't overturn because of the gravity. Furthermore, this analysis can be made for the legs gait sequence as discussed to be (leg4→leg2→leg3→leg1). During the swing period of leg4, the area of support triangle is framed by the tips of the supporting legs (leg1, leg2, leg3), and it has a symbolic as Tr1. In the swing period of the leg2, the support triangle is shaped by the tips of the supporting legs (leg1, leg3, leg4). And this triangle has a symbolic as Tr2. At this moment the intersection of the Tr1 and Tr2 resulted in another triangle named as DST1. In which the COG projection must be inside this triangle to ensure the static stability when leg4 and leg2 are in swinging phase. By the same way, the edges of supporting triangles Tr3 and Tr4 resulted during the period of swing leg3 and leg1, respectively. The intersection of the Tr3 and Tr4 resulting another triangle named as DST2. By the same condition that the COG must be inside the triangle DST2 triangle when leg3 and leg1 are in swinging phase. Many series of DSTs (Double Supporting Triangles) will be shaped during the movement of the quadruped robot, thus, to ensure the statically stable condition the planned trajectory of the COG projection must be planned along the DST's progression sequence: DST1→DST2→DST3→DST4.... …(18) This grantee that the locomotion of the quadruped robot is statically stable [10]. The DST sequence of the quadruped robot motion is illustrated in Figure (7). Fig. 7. Series of Double Support Triangles on the ground frame [10]. Firas A. Raheem Al-Khwarizmi Engineering Journal, Vol. 14, No. 2, P.P. 93- 106 (2018) 99 The main advantage of making the walking locomotion of the quadruped periodically repeated is to achieve a constant velocity which leads to making the acceleration on the quadruped body is equal to zero, which reduces the disturbances on the swinging legs. 5.2 Stability Margin Analysis In case of quadruped robot creeping gait walking, and according to static stability condition, the stability margins are utilized the shortest distance between the COG projection of the quadruped robot to the boundary limits of the supporting pattern [13]. The following Figures explaining the mathematical model of a quadruped robot with the analysis of stability margins where the calculations of these stability margins as follows (a) (b) Fig. 8. Quadruped robot when leg4 is swinging (a) quadruped robot (b) Quadruped support triangle The first case, shown in Figure (8-a) when the quadruped robot leg4 swing in the air and legs (1, 2, and 3) on the ground. The supporting area will have three triangles, each triangle has a defined area (Area1, Area2, Area3) and the minimum distance from one of the supporting lines (L1, L2, L3) to the perpendicular COG projection on the ground which is denoted as (T1, T2, T3). Then the analysis of the stability margins at this period is showing in Figure (8-b) as: TUV�� = �� W 1 1 1�&'X �� � =&'X =� = Y Where (�&'X, =&'X) is the coordinate of the center of gravity of the quadruped robot on the ground. (��, =�) is the coordinate of leg1 position. (� , = ) is the coordinate of leg3 position. So, by expanding and simplify this matrix TUV�� can be calculated as the following: TUV�� = �� [\�� � �&'X]\= � =&'X] � \� � �&'X]\=� � =&'X]^ E� = _4�� � � 5� � 4=� � = 5� �� = 2 ∗ 4`0,)aNa 5 ...(19) By the same way, TUV�� and TUV� can be calculated: TUV�� = �� [\�� � �&'X]\=� � =&'X] �\�� � �&'X]\=� � =&'X]^ E� = _4�� � ��5� � 4=� � =�5� �� = 2 ∗ 4`0,)GNG 5 ...(20) TUV� = �� [\�� � �&'X]\= � =&'X] � \� � �&'X]\=� � =&'X] ̂ E� = _4�� � � 5� � 4=� � = 5� � = 2 ∗ 4`0,)bNb 5 …(21) Finally, the first stability margin �c� is the minimum of these three margins: �c� = c�� 4��, ��, � 5 …(22) According to the quadruped robot creeping gait sequence, the robot swinging leg 2 while the other legs (leg1, leg3, and leg4) are on the ground, �c� can be found as the following (a) Firas A. Raheem Al-Khwarizmi Engineering Journal, Vol. 14, No. 2, P.P. 93- 106 (2018) 100 (b) Fig. 9. Quadruped robot when leg2 is swinging (a) quadruped robot 3D view (b) Quadruped support triangle. The second case, showing in Figure (9-a) when the leg2 is swinging in the air and the other legs (1, 3, and 4) will be on the ground. Then the analysis of the stability margins at this period is showing in Figure (9-b) as: TUV�� = �� [\�� � �&'X]\= � =&'X] � \ � � �&'X]\=� � =&'X]^ E� = _4�� � � 5� � 4=� � = 5� �� = 2 ∗ 4`0,)aNa 5 …(23) By the same way, TUV�� and TUV� can be calculated: TUV�� = �� [\�� � �&'X]\= � =&'X] � \� � �&'X]\=� � =&'X]^ E� = _4�� � � 5� � 4=� � = 5� �� = 2 ∗ 4`0,)GNG 5 …(24) TUV� = �� [\� � �&'X]\= � =&'X] � \� � �&'X]\= � =&'X]^ E� = _4� � � 5� � 4= � = 5� � = 2 ∗ 4`0,)bNb 5 …(25) Finally, the second stability margin �c� is the minimum of these three margins: �c� = min 4��, ��, � 5 …(26) According to the quadruped robot creeping gait sequence, the robot swinging leg 3 and the other legs (leg1, leg2, and leg4) are on the ground, the �c can be found as the following (a) (b) Fig. 10. Quadruped robot when leg3 is swinging (a) quadruped robot 3D view (b) Quadruped support triangle. The Third case, showing in Figure (10-a) when the leg3 is swinging in the air and the other legs (1, 2, and 4) will be on the ground. Then the analysis of the stability margins at this period is showing in Figure (10-b) as: Area� = �� h4x� � xjkl 54y� � yjkl 5 � 4x� � xjkl 54y� � yjkl 5n L� = _4x� � x�5� � 4y� � y�5� T� = 2 ∗ 4qrstaua 5 …(27) By the same way, TUV�� and TUV� can be calculated: Area� = �� h4x� � xjkl 54y � yjkl 5 � 4x � xjkl 54y� � yjkl 5n L� = _4x� � x 5� � 4y� � y 5� T� = 2 ∗ 4qrstGuG 5 …(28) Area = �� h4x� � xjkl 54y � yjkl 5 � 4x � xjkl 54y� � yjkl 5n L = _4x� � x 5� � 4y� � y 5� T = 2 ∗ 4qrstbub 5 …(29) Firas A. Raheem Al-Khwarizmi Engineering Journal, Vol. 14, No. 2, P.P. 93- 106 (2018) 101 Finally, the Third stability margin �c is the minimum of these three margins: �c = min 4��, ��, � 5 …(30) According to the quadruped robot creeping gait sequence, the robot swinging leg 1 and the other legs (leg2, leg3, and leg4) are on the ground, the �c can be found as the following (a) (b) Fig. 11. Quadruped robot when leg3 is swinging (a) quadruped robot 3D view (b) Quadruped support triangle. The Fourth case, showing in Figure (11-a) when the leg1 is swinging in the air and the other legs (2, 3, and 4) are on the ground. Then the analysis of the stability margins at this period is showing in Figure (11-b) as: Area� = �� h4x� � xjkl 54y � yjkl 5 � 4x � xjkl 54y� � yjkl 5n L� = _4x� � x 5� � 4y� � y 5� T� = 2 ∗ 4qrstaua 5 …(31) By the same way, TUV�� and TUV� can be calculated: Area� = �� h4x� � xjkl 54y � yjkl 5 � 4x � xjkl 54y� � yjkl 5n L� = _4x� � x 5� � 4y� � y 5� T� = 2 ∗ 4qrstGuG 5 …(32) Area = �� h4x � xjkl 54y � yjkl 5 � 4x � xjkl 54y � yjkl 5n L = _4x � x 5� � 4y � y 5� T = 2 ∗ 4qrstbub 5 …(33) Finally, the Fourth stability margin �c is the minimum of these three margins: �c = min 4��, ��, � 5 … (34) 6. Simulation and Results This section is showing the simulation results of the analyzed stability margins which are necessary for the quadruped robot walking stability. According to the robot legs sequence, the results are achieved as the following (a) (b) Fig. 12. Showing (a) Leg 4 in swinging phase (b) The stability margin vw>(cm) when leg4 is swing. (a) 50 40 X 30 20 10 0 0 10 20 Y 30 40 Leg 4 swing Support triangle Firas A. Raheem Al-Khwarizmi Engineering Journal, Vol. 14, No. 2, P.P. 93- 106 (2018) 102 (b) Fig. 13. Showing (a) Leg 2 in swinging phase. (b) The stability margin vwC(cm) when leg2 is swing. (a) (b) Fig. 14. Showing (a) Leg 3 in swinging phase. (b) The stability margin vwD(cm) when leg3 is swing. (a) (b) Fig. 15. Showing (a) Leg 1 in swinging phase. (b) The stability margin vwx(cm) when leg1 is swing. From figures shown above, it can be seen that all stability margins (Sm�, Sm�, Sm , and Sm ) have positive values. These values are sufficient to satisfy a stable walking for the quadruped robot creeping gait as assumed in this paper. The results of the stable creeping gait sequence for all legs have been converted into joint angles (��,��,� ) for each leg by using inverse kinematics analysis. So, the results contain the changing of each leg three angles and the change of distance of each leg-tip along X-axis direction and Z-axis direction during the creeping gait walking of quadruped robot. The following Figures (12-15) illustrate the quadruped robot configurations and results. (a) (b) Fig. 16. When Leg 4 is swinging. (a) Leg 4 moved one step along X-axis. (b) Leg 4 moved one step along Z-axis. S m 2 ( c m ) 0 5 0 100 1 50 200 250 3 00 350 1 2 3 4 5 6 Iteration S m 3 ( c m ) sm3 Firas A. Raheem Al-Khwarizmi Engineering Journal, Vol. 14, No. 2, P.P. 93- 106 (2018) 103 (a) (b) Fig. 17. When Leg 2 is swinging. (a) Leg 2 moved one step along X-axis. (b) Leg 2 moved one step along Z-axis. (a) (b) Fig. 18: When Leg 3 is swinging. (a) Leg 3 moved one step along X-axis. (b) Leg 3 moved one step along Z-axis. (a) (b) Fig. 19. When Leg 1 is swinging. (a) Leg 1 moved one step along X-axis. (b) Leg 1 moved one step along Z-axis. From above Figures it can be seen the following: 1. In Figure (16): leg4 tip is swinging and moving from (-5 cm) to (-3.5 cm) along X-axis and (1.5 cm) along Z-axis. 2. In Figure (17): leg2 swinging and moving from (5 cm) to (7.4 cm) along X-axis, also the height of the lifting tip is (1.5 cm) along Z-axis. 3. In Figure (18): leg3 swinging and moving from (-5 cm) to (-2 cm) along X-axis, also the height of the lifting tip is (1.5 cm) along Z-axis. 4. In Figure (19): leg1 swinging and moving from (5 cm) to (8.5 cm) along X-axis, also the height of the lifting tip is (1.5 cm) along Z-axis. As discussed in this paper, the quadruped walking along X-axis only while the Y-axis is fixed, there is no changing in the Y-axis during the walking of the quadruped robot. Figure (20) showing this case. Firas A. Raheem Al-Khwarizmi Engineering Journal, Vol. 14, No. 2, P.P. 93- 106 (2018) 104 Fig. 20. Showing the changing of the leg tips along Y-axis during the quadruped walking. During the creeping gait sequence, the quadruped robot leg angles are changing. These angles are (��,�� and� ) of each leg. This changing is showing in Figure (21) (a) (b) (c) Fig. 21. changing of angles (a) Coxa-angle (b) Femur-angle (c) Tibia-angle From Figure (21) shown above it can be seen the following: 1. Figure (21-a) is the changing of leg1 angle��. This angle is varying between (15 degree) to (-15 degree). 2. Figure (21-b) is the changing of leg1 angle��. This angle is varying between (30 degree) to (65 degree). 3. Figure (21-c) is the changing of leg1 angle� . This angle is varying between (120 degree) to (140 degree). 7. Conclusion The main focus of this paper is how to analyze the quadruped robot creeping gait walking and then to derive, verify and prove that the robot is statically stable during walking. In order to satisfy the goal of this work, first, the full forward and inverse kinematics model has been derived and used for walking steps stability validation. Also, the interaction between creeping gait sequence made for this robot and the geometric modeling of robot legs-tip has been fully derived for finding all the static stability margins during walking. The results clearly prove that the stability margins have the positive values which guarantee to keep robot COG inside the supporting triangle during the swing phase of one leg while other legs are on the ground. Furthermore, a future work development is needed to analyze and enhance the quadruped robot walking on the hard and rough terrain. Notation COG Center of gravity �&'X Center of gravity in X-axis =&'X Center of gravity in Y-axis x1 leg1 tip in X-axis x2 leg2 tip in X-axis x3 leg3 tip in X-axis x4 leg4 tip in X-axis y1 leg1 tip in Y-axis y2 leg2 tip in Y-axis y3 leg3 tip in Y-axis y4 leg4 tip in Y-axis ��, ��, � distance from stability margins to the COG Greek letter z Stride line Firas A. Raheem Al-Khwarizmi Engineering Journal, Vol. 14, No. 2, P.P. 93- 106 (2018) 105 8. References [1] Kumar, T. Komal, R. Kumar, S. Mogily, D. Goel, N. Pareekutty, S. V. Shah, and K. Madhava Krishna. "Implementation of gaits for achieving omnidirectional walking in a quadruped robot."Proceedings of the 2015 Conference on Advances In Robotics - AIR '15(2015). [2] Rebula, J. R., Neuhaus, P. D., Bonnlander, B. V,.Johnson, M. J., Pratt, J. E. (2007). “A controller forthe littledog quadruped walking on rough terrain”.ICRA, 1467-1473. [3] Raibert, M., Blankespoor, K., Nelson, G., Playter, R.(2008). Bigdog, the rough-terrain quadruped robot. InProc, of the 17th World Congress, 17(1), 10822-10825. [4] Shimoyama, I., Miura, H., Horita, M., Hattori, Y. (1991). Analysis of walk for quadruped. IROS,1541-1544. [5] Fukuoka, Y., Kimura, H., Cohen, A. H. (2003). Adaptive dynamic walking of a quadruped robot onirregular terrain based on biological concepts. The Int.J. of Robotics Research, 22(3-4), 187-202. [6] Ma, S., Tomiyama, T., Wada, H. (2005). Omnidirectional static walking of a quadrupedrobot.Transactions on Robotics, 21(2), 152-161. [7] Lee, T. T., Shih, C. L. (1986). A study of the gaitcontrol of a quadruped walking vehicle. J. of Roboticsand Automation, 2(2), 61-69. [8] Kim, B. H. (2009). Centroid-based analysis ofQuadruped-robot Walking Balance. Int. Conf. onAdvanced Robotics, 1-6. [9] Firas A. Raheem, Hind Z. Khaleel, "Hexapod Robot Static Stability Enhancement using Genetic Algorithm Al-Khwarizmi Engineering Journal,Vol. 11, No. 4, P.P. 44- 59, 2015. [10] Chen, Binbin, Shaoyuan Li, and Dan Huang. "Quadruped robot crawl gait planning based on DST." Proceedings of the 33rd Chinese Control Conference (2014). [11] Liu, Fei. Research about "Gait Planning and Balance Control ofQuadruped Robot", Doctoral dissertation, University of Scienceand Technology of China, 2010. [12] Buchli, Jonas and Kalakrishnan, Mrinal and Mistry, Michaeland Pastor, Peter and Schaal, Stefan, Compliant quadruped locomotion over rough terrain, The 2009 IEEE/RSJ International Conference on Intelligent Robots and Systems: 814–820, 2009. [13] Liu, F., D. Wu, and K. Chen. "The Simplest Creeping Gait For A Quadruped Robot". Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science 227.1 (2012): 170-177. Web. [14] Pack DJ and Kak A. A simplified forward gait control fora quadruped walking robot. In: Proceedings of the IEEE/RSJ/GI international conference on intelligent robots and systems, Munich, Germany, 12–16 September 1994, .١٠١٢ pp.1011–1018. New York: IEEE. [15] Mănoiu-OlaruSorin, NiţulescuMircea, “MatlabSimulator for Gravitational Stability Analysis of a HexapodRobot,” The Romanian Review PrecisionMechanics, Optics &Mechatronics, 158 2011, No. 39 Mechanics, Optics &Mechatronics, 2011. [16] FreyrHardarson,"Stability analysis and synthesis ofstatically balanced walkingfor quadruped robots",Doctoral ThesisMechatronics LabDepartment of Machine DesignRoyal Institute of Technology, KTHS-100 44 Stockholm, Sweden,2002. [17] Y. Ito, S. Nozawa, J. Urata, T. Nakaoka, K. Kobayashi, Y. Nakanishi, K. Okada, and M. Inaba, “Development and verification of life- size humanoid with high-output actuation system,” in IEEE International Conference on Robotics and Automation(ICRA), 2014, pp. 3433–3438. [18] H. Khan, S. Kitano, M. Frigerio, M. Camurri, V. Barasuol, R. Featherstone, D. G. Caldwell, and C. Semini, “Development of the lightweight hydraulic quadruped robot - MiniHyQ,” in IEEE International Conference on Technologies for Practical Robot Applications (TEPRA), 2015. [19] A. Winkler, C. Mastalli, I. Havoutis, M. Focchi, D. G. Caldwell, and C. Semini, “Planning and execution of dynamic whole- body locomotion for a hydraulic quadruped on challenging terrain,” in IEEE International Conference on Robotics and Automation (ICRA), May 2015. [20] V. Barasuol, J. Buchli, C. Semini, M. Frigerio, E. R. D. Pieri, and D. G. Caldwell, “Areactive controller framework for quadrupedal locomotion on challenging terrain,” in IEEE International Conference on Robotics and Automation (ICRA), 2013. [21] Yang, Chenguang, et al. Advanced Technologies in Modern Robotic Applications. Springer Singapore, 2016. )2018( 93-106، صفحة 2د، العد14دجلة الخوارزمي الهندسية المجلمفراس عبد الرزاق رحيم 106 تحليل و محاكاة روبوت رباعي االرجل في حالة المشي **مرتضى خضير فليح *فراس عبد الرزاق رحيم الجامعة التكنلوجية /هندسة السيطرة والنظمقسم *,** uotechnology.edu.iq60124@البريد االلكتروني:* .uotechnology.edu.iq60354cse@البريد االلكتروني:** الخالصة على قدرته فضال عنعلى االرض المسطحة واالنتقال المشيروبوت رباعي االرجل له القابلية على القيام بمهمات عديدة منها ربوت هوالكوادروبيد على محاكاة الحيوانات او الحشرات رباعية االرجل من خالل تقليد الحركات التي تقوم بها هذه الحيوانات. إن المشي على االرض الوعرة. يقوم هذا الروبت القفز. ان حركة المشي التي تم تنفيذها في هذه الورقة هي والركض معالحركات التي يستخدمها هذا الروبوت في االنتقال هي: المشي والركض أشهرمن واحدة بعد االخرى خالل مدة زمنية معينة وبعد وصول الرجل الى المكان المناسب تبدأ رجل اخرى بالحركة. يستخدم هذا االيقاع في عبارة عن تحريك رجل االرجل سوفرفع احدى خالل عمليةانتقال الرجل المرفوعة. مدةصقة على االرض لتحقيق االستقرارية خالل لتم أرجلالحركة لكي نضمن بقاء ثالثة خالل هذه الحركات لذلك اثلث ناتج من اسقاط نقاط تثبيت بقية االرجل الثالث مع سنتر الروبوت على االرض.ان مهمة هذا الروبوت هو البقاء مستقريتكون م ) في هذه الورقة ايضا تم تحليل حركة الروبوت بالنسبة للمشي Stability Marginsيحتاج الى تحليل االستقرارية عن طريق حساب الـ ( ) ومن ثم حساب معادالت الكاينماتيك الخاصة للحصول على موقع االرجل على االرض والالزمة لتحقيق االستقرارية الثابتة Creeping Gaitطيء(الب للروبوت اثناء المشي. ومن بعدها تم الحصول على النتائج الالزمة لتحقيق هذه االستقرارية من خالل مواقع االرجل على االرض.