# Analysis of Stiffness and Energy Consumption of Nonlinear Elastic Joint Legged Robot

Authors: Dongliang Chen; Jindong Zhang; Xutao Weng; Yunjian Zhang; Zhonghui Shi
Journal: Applied Bionics and Biomechanics
 (2020)

---

## Abstract

In order to reduce the energy consumption of the legged robot in walking, this paper designs a kind of nonlinear elastic joint from the flexible variable-stiffness joint based on the mammal walking on the limb and optimizes the leg structure of the legged robot. The motor is rigidly connected to the articulated lever. When the lever is accelerated or decelerated, the elastic unit is introduced. The system can be considered as a special variable-rate elastic system. This paper will study it from theory and simulation experiments. Based on the dynamic analysis, a functional relationship between the output torque and the torsion spring stiffness and between the energy consumption and the torsion spring stiffness was established. By finding the extremum, the two optimum torsional spring stiffness that can minimize the required output average torque and the energy consumed during one cycle of motion were deduced. The results show that using this design in a reasonable position can effectively reduce the energy consumption of the system and can achieve up to a 50% reduction in energy consumption.


---

## Body

## 1. Introduction

Compared with the wheeled robot, the foot robot has great advantages in adaptability to the external complex environment and has great research value [1]. With the increase of the number of feet, the structure of the foot robot is more redundant, but it improves the fault-tolerant ability of the robot: the damage of one leg will not have a great impact on the robot, it can move smoothly with the cooperation of other legs. However, the foot robot needs to coordinate its limbs, and its motion efficiency is relatively low, and its energy loss is relatively large. In the process of movement, the friction loss between the body and the ground as well as between the mechanisms will also be dissipated by the heat generated by the driving winding due to the need to maintain a fixed joint angle in the process of standing still. Therefore, for the foot robot, how to reduce energy consumption to improve energy utilization is particularly important, which will bring great performance improvement to the foot robot.In recent years, the research on energy-saving optimization of foot robot can be classified into the following four categories: (1) study the energy conversion of bionic joint by analyzing the physiological structure of leg joint of foot animal [2]; (2) optimize the end trajectory curve under different conditions [3]; (3) reasonable gait planning and optimal gait generation parameters [4]; and (4) how to reasonably distribute the foot force [5]. A large number of experiments have shown that human beings and some creatures use their own physiological structure to complete energy transfer and transformation when they are in motion, which is a highly coupled function brought by evolution and the body [6]. For example, in the movement system composed of human bones, muscles, and tendons, the flexible components of muscles and tendons produce the effect similar to that of springs, which are regularly shaped in motion. Variable can reduce energy consumption and overall impact of the system, which can be approximately considered as a load spring inverted pendulum slip model [7, 8]. Inspired by this idea, this paper introduces the spring design idea to the joint motion of the foot robot for energy conversion. Nowadays, many scholars have studied the elastic joint of robot and discussed the structure and control method of the elastic joint, but most of them focus on the performance research of damping and reducing the impact force [9]; this paper will study the energy consumption brought by it.However, the flexibility of the robot joint has a significant impact on the positioning accuracy of the legs and the accuracy of the end of the trajectory when the robot moves [10], and the degree of freedom of its control system will increase by one level, resulting in its control more complex and easy to lose stability. Therefore, a large number of scholars have studied the impact on the control model after introducing the flexibility of the joint. Nowadays, there are many significant achievements in system modeling [11] and dynamic modeling [12] of robot with elastic joints. However, for the foot robot with high redundancy, such as the six-legged robot, the lack of accuracy of the foot end trajectory of one leg is not enough to cause the stability imbalance of the whole robot, but the resulting energy consumption is reduced, which can greatly improve the endurance of the foot robot. Secondly, most of the current research uses linear elastic elements and will participate in the whole movement. In this paper, a joint motor is used to drive the joint directly, and a torsion spring is installed at the joint which is in contact with the bar. By improving some elastic joints, a nonlinear elastic system is formed. Only when the joint motor decelerates to the stop and the initial speed is zero to the acceleration, the elastic joint will start to play an obvious role. In this way, it can help the joint motor decelerate in the deceleration process and accelerate in the acceleration process, so as to achieve an energy conversion function and reduce energy loss. In the following chapters, we will start with a simplified model of nonlinear elastic connections. Through dynamic analysis, the relationship between motor output torque and spring stiffness coefficient is found. Then, through the analysis of energy consumption, the value of joint stiffness coefficient is found. After that, we will design some simulation experiments to verify the effect of nonlinear elastic connection on the system. Finally, the relevant conclusions are summarized.


## 2. Dynamic Modeling and Analysis

The single-leg structure of the legged robot model and the prototype are shown in Figure1. The robot can switch arbitrarily between a hexapod robot and a spherical robot. Each leg has three degrees of freedom and three joint motors. This paper focuses on the second joint, that is to say, only the leg joints are subjected to a simulation experiment with nonlinear elastic joints. For each revolute joint, a similar dynamic model will be generated when the torsion spring of the elastic joint is assembled in a suitable position. They have the same properties and only differ in numerical values.Figure 1
(a) Legged robot single-leg model; (b) a hexapod robot prototype.
(a)(b)The single-leg simplified model with nonlinear elastic joints is shown in Figure2. Since this paper only does theoretical analysis and simulation, in order to simplify the model, the ideal position is taken during assembly without actual assembly. Some elastic joints are horizontal and symmetrical with the axis of the revolute joint.Figure 2
Single-leg simplified model with nonlinear elastic joints.Through dynamic modeling, the relationship between joint torque and joint motion parameters of a nonlinear elastic joint robot can be deduced. At the same time, the relationship between energy consumption and joint movement can also be deduced. This paper uses the Lagrangian method to derive the dynamic equations of the elastic joint model. The elastic part of the torsion spring has no moment of inertia and is elastically deformed only in one direction, ignoring the damping of the elastic unit.Assume that the initial momentum of the system is zero. In the theoretical derivation and simulation experiments, this paper concludes that the trajectory is accurately tracked, and there is no error such as offset and phase difference. Figure3 shows the simplified nonlinear elastic joint mechanism. The parameters of the mechanism diagram are shown in Table 1; the torsion spring that mainly acts in the elastic joint is defined as the main torsion spring k1, and its stiffness is far greater than that of the secondary torsion spring k2 that only reduces the impact, that is, k1 is much larger than k2 and k2≈0, and the barycentric coordinate of the lever is set as x,y.Figure 3
Simplified nonlinear elastic joint mechanism.Table 1
Nonlinear elastic joint mechanism diagram parameters and meaning.
ParameterDefinitionθJoint ground angleθkJoint and torsion spring contact anglek1The main torsion spring stiffnessk2The secondary torsion spring stiffnessLBody lengthlHalf body lengthMOutput torqueIThe moment of inertia of bodymBody massThe joint space is divided into two parts: an elastic joint with a very small stiffness and an elastic joint with a large stiffness. In order to simplify the equation, it is necessary to make the following assumption. In the collision process, the collision is stable and the time is not counted. The collision is surface contact. The position and direction of the lever and the contact surface do not change during the collision. The inertia of each component does not change, and the deformation does not occur during the collision. Lagrangian equation is still established during collision [13].As shown in Figure3, we separately model the dynamics of different stiffness. In an ideal state, assuming that the rigid joint driving torque is Mθ, the torque of the torsion spring and the joint rotation angle in the elastic system are, respectively, the functions of the relationship Tk1θ and Tk2θ. When the lever has not yet reached the torsion spring with the large stiffness, the required torque for the joint can be given by
(1)Tq1θ=Mθ+Tk2θ.The required torque after the lever comes into contact with the large stiffness torsion spring can be given by(2)Tq2θ=Mθ+Tk1θ+Tk2θ.After the trajectory curve is determined, the joint output torque curve in each period is also determined, that is, the relationship betweenTq1,Tq2 and time is a determined curve. The sum of the output torque of the joint and the torque generated by the elastic joint is a certain value at a certain point of time. At this time, the greater the torque generated by the elastic joint, the smaller the required output torque of the joint.However, in order to make the dynamic model more accurate, the two parts of the joint space are now integrated into a dynamic analysis. The D-H method [14] is used to establish the coordinate system of the joint. Only the second joint is taken as the research object.The kinetic equation of the robot is represented as follows:(3)E=12mẋ2+ẏ2+12lθ̇2,where x,y can be given by
(4)x=lsinθ,y=lcosθ.Equation (4) can be rewritten as
(5)E=12ml2θ̇2+12Iθ̇2.Gravity energy can be given by(6)P1=mgy=−mglcosθ.There is almost no elastic potential energy in the system before the torsion spring with a large stiffness is touched, and the stiffness coefficient can be regarded as linear after the main torsion spring is contacted. So it can be considered as a special nonlinear variable stiffness elastic model in one cycle. The elastic potential energy can be given by(7)P2=12kθαΔθ2,where α is the angle conversion factor when the torsion spring and the rotation center do not coincide with each other, and α keeps the value of the rotation angle of the torsion spring and the rotation angle of the joint united. Δθ is the rotation angle of the torsion spring with respect to the joint and can be approximated as the joint rotation angle θ, because the lever and the motor are rigidly connected. kθ is a function of the nonlinear torsion spring for the joint rotation, which can be expressed in this paper by a piecewise function:
(8)k=0,θ−θk≤0,k1,θ−θk>0.The total potential energy can be given by(9)P=P1+P2=−mglcosθ+12kθαΔθ2.The Lagrangian operator can be given by(10)L=E−P=12ml2θ̇2+12Iθ̇2+mglcosθ−12kθαΔθ2.The Lagrangian kinetic equation can be given by(11)∂L∂θ=−mglsinθ±k̇θαΔθ,(12)ddt∂L∂θ̇=ml2θ¨+Iθ¨.According to the Lagrange equation, the joint torque can be obtained as(13)M=ddt∂L∂θ̇−∂L∂θ=ml2θ¨+Iθ¨+mglsinθ∓k̇θαΔθ.Assume that the joint output rotation function isft, then the output torque after contact with the elastic joint can be given by
(14)Mk=ml2f¨t+If¨t+mglsinft−kftαft−θk.Equation (14) has the same conclusion as equations (1) and (2). Equation (14) shows that before the end of the stance phase and the swing phase, the robot needs a joint motor to decelerate. With the participation of the elastic joint, due to the compression of the torsion spring, the torsion spring will produce an increasingly larger countertorque. It can effectively reduce the output torque of the joint motor under the same combined torque. At the beginning of a new swing phase or stance phase, the joint motor output torque is required to provide sufficient acceleration to the joint. In the model with elastic joints, the torsion spring is compressed when the previous movement decelerates. So the torsion spring and the joint motor collectively output torque in the acceleration phase. When the required torque is the same, the output torque of the joint motor can be greatly reduced. Therefore, theoretically, the design of the elastic joint can reduce the output peak torque of the joint motor, which can effectively reduce the energy consumption of the robot during walking.In order to obtain the best stiffness function of the torsion spring with the purpose of minimizing the output torque, differentiate equation (14) of time and let its value be zero:
(15)ḟk=dMkdt=0.kMθ can be obtained from equation (15). If f¨kMθ≥0, then kMθ (N∗m/rad) is the function of the torsion spring stiffness that minimizes the output torque under this input function.For the robot’s trajectory design, many scholars have done a lot of researches [15, 16], such as composite cycloids, modified elliptical trajectories, and classic multi-item fitting curves. For the joint tip trajectory planning with continuous smooth acceleration and the appropriate initial conditions for position and velocity, a simpler sinusoidal trajectory will be used, i.e.,
(16)ft=Asinωt+ϕ+B.Equation (14) can be rewritten as
(17)Mk=−ml2+IAω2sinωt+ϕ+mglsinAsinωt+ϕ+B−kAsinωt+ϕ+BαAsinωt+ϕ+B−θk.Differentiate equation (17) of time and ignore the effect of gravity. When t∈arcsinθk−B/A−ϕ/ω,π/2ω, equation (17) can be rewritten as
(18)ḟk=dMkdt=−ml2+IAω3cosωt+ϕ−Ak2αωcosωt+ϕ.Let equation (18) be equal to zero and get the result as follows:
(19)k2=ω2ml2+Iα.Iff¨k2>0, then kM=k2 is the optimal stiffness that will make the output torque the smallest torsion spring.For different trajectories, the specific positions of the main torsion springs for joint contact need different analysis, and it is necessary to ensure that the torsion spring contacts the deceleration when the joint driving torque is greatly reduced, and a symmetrical acceleration torque is required. When using different driving torque curves for acceleration and deceleration, the situation will become more complicated and changeable.


## 3. Energy Consumption Analysis

Without considering other factors such as friction loss, the energy consumed by a single joint model with a nonlinear elastic joint in one exercise cycle can be expressed as(20)E=∫0Tpdt=∫0Tτωdt=∫0TMkθ̇dt,where p is the joint drive power, τ is the joint drive torque, ω is the joint rotation angular velocity, and θ is the joint rotation angle. From equation (20), it can be seen that the energy consumption is related to the integral of the output torque and the angular velocity at this moment, and the output torque and the joint angular velocity are both functions of time. Equations (16) and (17) can be expressed as
(21)E=∫0Tml2+If¨tḟt+ḟtmglsinft−ḟtkftαft−θkdt.If we want to obtain a stiffnesskE that minimizes E, then
(22)gkft=E.Letġk=0 to get kE. If g¨kE>0, then kE is the demand.Assume that the energy consumption ratioε is an evaluation standard for energy consumption [17]. This assessment of energy consumption is based on the metabolic effects of the organism [18, 19]. The energy consumption ratio means that the unit of weight requires a unit of energy to dissipate the distance. The smaller the value is, the higher the energy consumption can be obtained.
(23)ε=EmcgS,where mc represents the total mass of the robot and S is the standard of the robot’s moving distance.The rotational movement of the joint is the main movement, so make some changes to equation (23) to conform to the simulation movement in this paper. Let the moving standard S=θ and the total weight be replaced with the moment of inertia about the rotation axis, equation (23) can be rewritten as
(24)ε=∫0TMkθ̇dtIθ.


## 4. Simulation Experiments

### 4.1. Experiment Environment

Using Adams software for simulation experiments, the experimental model was modeled in Creo software and imported into Adams, as shown in Figure4. After removing the components that are not related to the movement, add the parameters of friction and force of the restraint and the motion parameter. Table 2 is the setting of the frictional force parameter of the joint rotation; Table 3 lists the parameter setting of the contact force between the lever and the torsion spring; Table 4 is the setting of the secondary torsion spring physical parameters; Table 5 lists the parameters of the main torsion spring physical. The joint adopts a general point driving method. After inputting the endpoint trajectory, the required torque value of the joint is sampled. In order to keep the motor from no-load spinning, the end of the joint is loaded with an equivalent mass load. The trajectory function is shown in
(25)ft=0.3sin2t,where t is time.Figure 4
Simulation model of robotic leg with nonlinear elastic joint.Table 2
Friction parameters of rotary motion.
ParameterValueMu static0.5Mu dynamic0.8Friction arm10Bending reaction arm10Pin radius10Static transition velocity0.1Max static friction deformation0.01Friction torque preload0Table 3
Parameter of contact force.
ParameterValueStiffness1.00E+05Force exponent2.2Damping10Penetration depth0.1Friction forceCoulombStatic coefficient0.3Dynamic coefficient0.1Static translation Vel100Friction translation Vel1000Table 4
The physical parameters of main torsion spring.
Spring parameterValueSpring stiffnessk1 (N∗mm/deg)400Damping coefficient (N∗mm∗s/deg)1Preload (N)0Table 5
The physical parameters of secondary torsion spring.
Spring parameterValueSpring stiffnessk1 (N∗mm/deg)0.01Damping coefficient (N∗mm∗s/deg)0.0001Preload (N)0The model with a nonlinear elastic joint was set as the experimental group, and the model without the elastic joint was the control group. The physical dimensions and simulation parameters of the two models were basically the same, as shown in Table6.Table 6
Parameters of two simulation models.
ParameterValueLeg lengthL84.000 mmCenter of mass to axis of rotationl47.141 mmEquivalent body massm3.611E-002 kgEquivalent moment of inertiaI107.059kg∗mm2Initial joint angleθ0 deg


### 4.2. Simulation Experiment Results

The simulation experiments found that under the same quality and external conditions, the output torque required by the experimental group model compared to the control group model was greatly reduced. The simulation results of the required output torque at the joint of the control group model are shown in Figure5, and the simulation results of the required output torque at the joint of the experimental group are shown in Figure 6.Figure 5
Required output torque of joint without nonlinear elastic joint.Figure 6
Output torque required at torque springK=400N∗mm/deg.We changed the stiffness of the torsion spring and carried out the simulation experiment, as shown in Figure7, which were the required output torque of the joint when the stiffness of the torsion spring was different.Figure 7
Required output torque for joints with different stiffness.It can be clearly seen from Figure6 that due to the limitation of the function kθ in equation (13), the joint torque is a normal output before contacting the main torsion spring. In Figures 5 and 6, there is a coincident torque output value under different stiffness. Figure 7 show the torsion spring has a stiffness ranging from 100N∗mm/deg to 1200N∗mm/deg. The required output torque of the joint first slowly decreases to a certain value and then begins to rebound. This rebound phenomenon is most obvious when the joint gets close to the limit of rotation. The torsion spring in this case has an optimum stiffness that minimizes the required output torque. The trajectory parameters and the model physical parameters are brought into equation (19), where θk=7.32deg and α=1.05, resulting in an optimum stiffness kM=713.547N∗mm/deg.The abovekM and nearby values were used as the stiffness of the simulation model. The results (kM) are shown in Figure 8. The experimental results showed that the simulation experiment was not exactly the same as the theoretical value and there was a certain error.Figure 8
The required output torque of the joint around the theoretical stiffnesskM.After several simulations aroundk=713N∗mm/deg, the simulation results showed that the required output torque was slightly compared with other stiffness at the same time around k=740N∗mm/deg. The following numerical calculation will verify its correctness.Matlab was used to process the discrete torque derived output torqueM from the Adams. After using the trapez function to integrate its data over a period, and then dividing it by the total time in one cycle, the average output torque under different stiffness per unit time was obtained. As shown in Figure 9, the simulation result around the theoretical derivation kM is shown in Figure 10.Figure 9
Average required output torque for different stiffness in one cycle.Figure 10
Average required output torque in one cycle aroundkM.Figure10 shows that the average moment around k=740N∗mm/deg is the lowest, and the maximum error between the simulated value and the theoretical value is about 3.7%.The trajectory curve was brought into equation (21). So that it was integrated over the interval θk to π/4, resulting in the torsional spring stiffness kE=784.717N∗mm/deg that optimized the energy. The simulation results of the required joint output torque under the stiffness around kE are shown in Figure 11.Figure 11
The required output torque of the joint under the theoretical stiffness aroundkE.Then, Matlab is used to multiply the derived angular velocityθ̇ and the output torque Mk discrete data at the sampling point, and the trapez function is used to integrate the data over a period of time. After bringing its value into equation (24), it got the energy ratio ε at different stiffness, as shown in Figure 12. The simulation result around the theoretical derivation value kE is shown in Figure 13.Figure 12
Energy consumption ratio under different stiffness.Figure 13
Energy consumption ratio under the condition of the stiffness aroundkE.Figure13 shows that the average moment around k=800N∗mm/deg is the lowest, and the maximum error between the simulated value and the theoretical value is about 1.95%.According to the data of simulation experiments, the difference between the energy ratio in the case of nonlinear elastic joints and the energy ratio in the case of the optimal energy-saving torsion spring stiffness was first calculated. Then, calculate the ratio of the difference and the energy consumption when no elastic joint was added. It can be concluded that energy consumption of up to 52.81% can be reduced in the case of an optimal energy-saving stiffness.


### 4.3. Discussion of the Experimental Results

Chen et al. [20] use a torsional spring to participate in the movement of the quadrupedal robot joint, assisting the hydraulic drive to significantly reduce the hydraulic output during walking, thereby reducing the energy consumption. Zhou and Fu [21] propose a method based on biomechanics and bionic control strategy for the compliant quadruped robots. It mainly includes mechanical compliance element and control compliance element that can reduce the contact force both on hip and knee joints. Finally, the experiment results proved to effectively optimize joint torque and prevent damage to the robot. Ma et al. [22], inspired by small jumping animals, designed a leapfrog robot named “Grillo.” Its passive forelimbs can cushion the landing impact and can be converted into elastic potential energy and released during jumping to increase the peak power output of the motor. The above papers and the theoretical analysis of this paper and the simulation experiments indicate that adding elastic elements at certain positions of the robot joint can indeed reduce the energy consumption of the robot during walking, but may cause instability of the output torque. It may cause the decline of the stability and precision of endpoint of the robot, so it is necessary to weigh the advantages and disadvantages brought about before the introduction of the elastic joint.


## 4.1. Experiment Environment

Using Adams software for simulation experiments, the experimental model was modeled in Creo software and imported into Adams, as shown in Figure4. After removing the components that are not related to the movement, add the parameters of friction and force of the restraint and the motion parameter. Table 2 is the setting of the frictional force parameter of the joint rotation; Table 3 lists the parameter setting of the contact force between the lever and the torsion spring; Table 4 is the setting of the secondary torsion spring physical parameters; Table 5 lists the parameters of the main torsion spring physical. The joint adopts a general point driving method. After inputting the endpoint trajectory, the required torque value of the joint is sampled. In order to keep the motor from no-load spinning, the end of the joint is loaded with an equivalent mass load. The trajectory function is shown in
(25)ft=0.3sin2t,where t is time.Figure 4
Simulation model of robotic leg with nonlinear elastic joint.Table 2
Friction parameters of rotary motion.
ParameterValueMu static0.5Mu dynamic0.8Friction arm10Bending reaction arm10Pin radius10Static transition velocity0.1Max static friction deformation0.01Friction torque preload0Table 3
Parameter of contact force.
ParameterValueStiffness1.00E+05Force exponent2.2Damping10Penetration depth0.1Friction forceCoulombStatic coefficient0.3Dynamic coefficient0.1Static translation Vel100Friction translation Vel1000Table 4
The physical parameters of main torsion spring.
Spring parameterValueSpring stiffnessk1 (N∗mm/deg)400Damping coefficient (N∗mm∗s/deg)1Preload (N)0Table 5
The physical parameters of secondary torsion spring.
Spring parameterValueSpring stiffnessk1 (N∗mm/deg)0.01Damping coefficient (N∗mm∗s/deg)0.0001Preload (N)0The model with a nonlinear elastic joint was set as the experimental group, and the model without the elastic joint was the control group. The physical dimensions and simulation parameters of the two models were basically the same, as shown in Table6.Table 6
Parameters of two simulation models.
ParameterValueLeg lengthL84.000 mmCenter of mass to axis of rotationl47.141 mmEquivalent body massm3.611E-002 kgEquivalent moment of inertiaI107.059kg∗mm2Initial joint angleθ0 deg


## 4.2. Simulation Experiment Results

The simulation experiments found that under the same quality and external conditions, the output torque required by the experimental group model compared to the control group model was greatly reduced. The simulation results of the required output torque at the joint of the control group model are shown in Figure5, and the simulation results of the required output torque at the joint of the experimental group are shown in Figure 6.Figure 5
Required output torque of joint without nonlinear elastic joint.Figure 6
Output torque required at torque springK=400N∗mm/deg.We changed the stiffness of the torsion spring and carried out the simulation experiment, as shown in Figure7, which were the required output torque of the joint when the stiffness of the torsion spring was different.Figure 7
Required output torque for joints with different stiffness.It can be clearly seen from Figure6 that due to the limitation of the function kθ in equation (13), the joint torque is a normal output before contacting the main torsion spring. In Figures 5 and 6, there is a coincident torque output value under different stiffness. Figure 7 show the torsion spring has a stiffness ranging from 100N∗mm/deg to 1200N∗mm/deg. The required output torque of the joint first slowly decreases to a certain value and then begins to rebound. This rebound phenomenon is most obvious when the joint gets close to the limit of rotation. The torsion spring in this case has an optimum stiffness that minimizes the required output torque. The trajectory parameters and the model physical parameters are brought into equation (19), where θk=7.32deg and α=1.05, resulting in an optimum stiffness kM=713.547N∗mm/deg.The abovekM and nearby values were used as the stiffness of the simulation model. The results (kM) are shown in Figure 8. The experimental results showed that the simulation experiment was not exactly the same as the theoretical value and there was a certain error.Figure 8
The required output torque of the joint around the theoretical stiffnesskM.After several simulations aroundk=713N∗mm/deg, the simulation results showed that the required output torque was slightly compared with other stiffness at the same time around k=740N∗mm/deg. The following numerical calculation will verify its correctness.Matlab was used to process the discrete torque derived output torqueM from the Adams. After using the trapez function to integrate its data over a period, and then dividing it by the total time in one cycle, the average output torque under different stiffness per unit time was obtained. As shown in Figure 9, the simulation result around the theoretical derivation kM is shown in Figure 10.Figure 9
Average required output torque for different stiffness in one cycle.Figure 10
Average required output torque in one cycle aroundkM.Figure10 shows that the average moment around k=740N∗mm/deg is the lowest, and the maximum error between the simulated value and the theoretical value is about 3.7%.The trajectory curve was brought into equation (21). So that it was integrated over the interval θk to π/4, resulting in the torsional spring stiffness kE=784.717N∗mm/deg that optimized the energy. The simulation results of the required joint output torque under the stiffness around kE are shown in Figure 11.Figure 11
The required output torque of the joint under the theoretical stiffness aroundkE.Then, Matlab is used to multiply the derived angular velocityθ̇ and the output torque Mk discrete data at the sampling point, and the trapez function is used to integrate the data over a period of time. After bringing its value into equation (24), it got the energy ratio ε at different stiffness, as shown in Figure 12. The simulation result around the theoretical derivation value kE is shown in Figure 13.Figure 12
Energy consumption ratio under different stiffness.Figure 13
Energy consumption ratio under the condition of the stiffness aroundkE.Figure13 shows that the average moment around k=800N∗mm/deg is the lowest, and the maximum error between the simulated value and the theoretical value is about 1.95%.According to the data of simulation experiments, the difference between the energy ratio in the case of nonlinear elastic joints and the energy ratio in the case of the optimal energy-saving torsion spring stiffness was first calculated. Then, calculate the ratio of the difference and the energy consumption when no elastic joint was added. It can be concluded that energy consumption of up to 52.81% can be reduced in the case of an optimal energy-saving stiffness.


## 4.3. Discussion of the Experimental Results

Chen et al. [20] use a torsional spring to participate in the movement of the quadrupedal robot joint, assisting the hydraulic drive to significantly reduce the hydraulic output during walking, thereby reducing the energy consumption. Zhou and Fu [21] propose a method based on biomechanics and bionic control strategy for the compliant quadruped robots. It mainly includes mechanical compliance element and control compliance element that can reduce the contact force both on hip and knee joints. Finally, the experiment results proved to effectively optimize joint torque and prevent damage to the robot. Ma et al. [22], inspired by small jumping animals, designed a leapfrog robot named “Grillo.” Its passive forelimbs can cushion the landing impact and can be converted into elastic potential energy and released during jumping to increase the peak power output of the motor. The above papers and the theoretical analysis of this paper and the simulation experiments indicate that adding elastic elements at certain positions of the robot joint can indeed reduce the energy consumption of the robot during walking, but may cause instability of the output torque. It may cause the decline of the stability and precision of endpoint of the robot, so it is necessary to weigh the advantages and disadvantages brought about before the introduction of the elastic joint.


## 5. Conclusions

(1)
This paper designs a nonlinear elastic joint and deduces it theoretically. It is proved that this design can effectively reduce the energy consumption of the legged robot when walking. The simplification of the model may bring about errors between the theoretical analysis and the simulation experiment. Even if there is an error, it can reduce about 50% of the energy consumption compared with the full-rigid joint robot without considering external friction and gravity(2)
This design will reduce the order of the robot’s control equation to a certain extent and increase the difficulty of control. However, for a more redundant robot, the reduction in energy consumption may increase the overall performance of the robot. This paper will be validated in later work. However, the position of the elastic joint in the overall mechanism is not flexible enough. If the joint rotation trajectory curve does not coincide with it, it may cause a reaction. That is to say, when the output torque of the joint motor accelerates the joint, the joint mechanism may come into contact with the torsion spring in advance, which causes the joint motor to drive to provide extra torque. And in the real world, there are constants of the acceleration of gravity. The upper and lower elastic joints cannot be assembled symmetrically. It is even more difficult to fit the trajectory to the elastic joints(3)
The future work of this paper is to give the appropriate elastic joint assembly position and the optimal stiffness for the more general end trajectory and try to design a universal adjustable nonlinear elastic joint



---

*Source: 8894399-2020-07-14.xml*