Design and transmission performance analysis of planthopper hip joint gear

To enhance the efficiency and load-carrying capacity of the jumping robot’s transmission system, the transmission characteristics of the planthopper hip gear were analyzed. The analysis revealed that the gear pair exhibited low transmission errors, high transmission efficiency, and significant torque capacity, making it suitable for high-speed, high-precision, and high-power density transmissions. A bionic design study of the jumping robot’s gear was conducted, drawing inspiration from the planthopper hip gear’s tooth shape and profile. A logarithmic spiral was assumed as the tooth profile based on its consistent pitch angle and curvature gradient characteristics. Parameters of the tooth profile of the bionic gear were calculated, and equations for the tooth surface were derived. A three-dimensional model of the bionic gear was constructed, and both its transmission efficiency and error were calculated. Using the finite element method, theoretical formulas for contact stress and bending stress in the bionic gear were deduced. The distribution characteristics of both bending and contact stress were analyzed, demonstrating that the bionic gear used in high-performance jumping robots possesses low transmission error, high transmission efficiency, and excellent contact strength, providing a theoretical foundation for enhancing their jumping capabilities.


Introduction
Jumping robots exhibit the remarkable capability to traverse obstacles tenfold or even several tens of times their size.This ability holds immense value and offers extensive application possibilities in interstellar exploration, terrain investigation, ore prospecting, and enemy detection.Most jumping robots store energy by propelling energy storage elements through a multi-stage variable-speed drive mechanism to enable their jumping motion. 1,2The pressure angle in the gear drive system of the jumping robot varies significantly, leading to substantial dynamic loads on the system. 3,4Gear pairs serve as fundamental components within the jumping robot transmission system and have a profound impact on transmission efficiency and load-carrying capacity, which, in turn, influence the jumping capacity of the robot.Therefore, the development of high-performance gears suitable for jumping robot transmission systems holds significant importance and offers extensive technical applications for enhancing the jumping abilities of these robots.
Jumping insects, after thousands of years of evolution, have excellent jumping abilities in natural environments.Based on biomimetic principles, the physiological structure of jumping insects is simulated to design mechanical structures with similar functions, which is an important way to improve the jumping performance of jumping robots.Zoologists Malcolm Burrows and Gregory Sutton, affiliated with the University of Cambridge, were the first to identify a gear-shaped biological structure in the planthopper, 5 situated in their hips (as depicted in Figure 1).This gear rapidly rotates and facilitates the synchronized movement of the planthopper's left and right legs during jumping.The planthopper, having evolved over an extensive period in its natural habitat, demonstrates a remarkable adaptation to its environment.The hip gear of the planthopper, a vital component in its jumping mechanism, exhibits characteristics of high synchronization accuracy and relative output torque. 5To meet the demands for high transmission efficiency and loadcarrying capacity in jumping robots, we have designed a new type of gear pair.This design draws inspiration from the structural features of the planthopper's hip gear pair, with the goal of achieving higher transmission efficiency, synchronization accuracy, and load capacity.A new type of gear is designed to meet the requirements of high-performance jumping robots for transmission efficiency and load capacity.
The tooth profile of the planthopper hip gear occurs naturally, while the logarithmic spiral is a naturally occurring curve known for its simplicity, structural stability, and uniform pitch angle.Incorporating the logarithmic spiral into the bionic design of the planthopper hip joint gears offers a fresh perspective.The logarithmic spiral's consistent pitch angle addresses issues related to high dynamic loads and unstable transmission, which stem from wide fluctuations in pressure angles within the gear transmission of jumping robots. 6urthermore, the curvature gradient characteristic of the logarithmic spiral tackle problems such as gear wear, power loss, and fatigue failure arising from variations in the curvature of gear tooth profiles. 7In recent years, logarithmic helix has gradually been applied to the tooth profile of bevel gear transmission, and the tooth profile of this type of bevel gear is still involute or circular. 8,9Yang 10 proposed that an imaginary rack cutter with the tooth shape of the planthopper could be used to generate a pair of gears.Multiple circular arcs were used to fit tooth profiles of the planthopper hip gear.
As a result, this paper introduces the hypothesis that the working tooth profile of the hip gear is a logarithmic spiral.Leveraging the tooth shape and tooth profile attributes of the planthopper hip gear, a novel bionic gear is developed, considering the mathematical properties of the logarithmic spiral.The meshing characteristics and transmission performance of the bionic gear are analyzed to enhance efficiency, thereby contributing to the improved efficiency and load capacity of the jumping robot's transmission system.

Analysis of transmission characteristics of planthopper hip gear
During the planthopper jumping process, the leg and tibia segments with the most significant influence on jumping are considered in the planthopper's jumping foot, while other limb segments are simplified.The weight of the jumping foot is not factored into the analysis, as it typically constitutes a small proportion of the insect's total body weight (generally less than 10%). 11s depicted in Figure 2, the planthopper's jumping structure is simplified to comprise the hip gears H 1 and H 2 , with O representing the gear contact point and M S as the center of gravity for the planthopper's body.H 1 B and H 1 C denote the leg segments, while AB and AC represent the tibia segments.The upward jumping force, F y , is generated at the point of contact between A and the ground.
The planthopper hip gears are always in a meshing state during the jumping process.Each joint angle is shown in Figure 3(a).Referring to the literature, 11 the leg-ground angle a during the planthopper jumping is obtained, and a polynomial function is used for fitting.The planthopper hip gear speed is acquired by the derivation of the leg-ground angle a, as depicted in  The maximum speed of the planthopper hip gear is 12000 r/min, indicating that the gear is used in high-speed transmission. 12ccording to the literature, 5 the maximum output torque T out of the planthopper hip gear drive is 0.7 mN.mm.The average weight M b of the planthopper is 6.6 mg.Its average body length L b is 4.1 mm, and the diameter D c of the planthopper hip gear is about 0.4 mm.To calculate the maximum output power P of the planthopper hip gear, the weight M of the planthopper hip gear pairs can be calculated in proportion to their length, as depicted in equation (1).
The power density of the planthopper hip gear drive is 1.3743 kW/kg, which is used in the high-power density drive. 13he transmission error u er of the gear refers to the relative difference between the actual angle u out of the output shaft and the theoretical angle u in =i ð Þ when the input shaft rotates at any angle, 14 as depicted in equation (9).
Where i is the transmission ratio.
Referring to the literature, 5 it is known that the whole jumping time of the planthopper is 78.4 ms, and the time error is 21 ms.The jumping acceleration time is 1.8 ms.t er is the error of the jumping acceleration time, as equation (3).
Where t in is the jumping acceleration time; t all is the whole jumping time; t aer is the time error.The driving gear rotates 90°during acceleration (u in ).And the average angular velocity v in of the driving gear is obtained as equation (4).
The planthopper hip gear pair ensures precise synchronization, i = 1.The average angular velocity v in of the driving gear is equal to the average angular velocity v out of the driven gear, that is, v in = v out .The transmission error u er can be obtained as equation (5).
The transmission error time t er of the planthopper hip gear is 0.5 ms, and the average angular velocity v in is 8333.3r/min.The transmission error of the planthopper hip gear is 36" at an angle of 360°angle, which is a high-precision transmission. 15,16ssume that the input power of the planthopper hip gear is P in and the output power during the jump  acceleration is P out .The angular displacement of the driving gear (H 1 ) is S in , and that of the driven gear (H 2 ) is S out .Because of the symmetry of the planthopper's jumping foot, it is known that F 1 = F 2 .Since the radii of the planthopper hip gears are equal, the arc length is S in = S out .The transmission efficiency of the hip gear is given by equation (6).
The transmission efficiency of the planthopper hip gear transmission is 99.97% as depicted in equation ( 6), which is a high-efficiency transmission device. 17

Analysis of tooth profile characteristics of planthopper hip gear
Mathematical properties of logarithmic spirals.Logarithmic spirals are some of the most amazing curves in nature, ranging from tornado and star orbits to nautilus and spiral fungal structures, as depicted in Figure 4. Conch shells grow naturally with the growth of organisms.And they behave with great strength and toughness.The spider's web, winding around the center, achieves infinite extension outwards in a logarithmic spiral shape.These phenomena show that logarithmic spirals are the products of natural selection.
The logarithmic spiral has the following properties: (1) the pitch angles of each point on the same logarithmic spiral are equal everywhere; (2) a logarithmic spiral is an equidistant curve; (3) its evolute and vertical traces are still logarithmic spirals; (4) the radius of gradient curvature is an increasing function of the polar angle.
Extraction and fitting of planthopper hip gear tooth.Considering the actual tooth profile of the planthopper hip gear, the edge detection method 17 is used to identify and extract its working tooth profile.Inspired by nature, the working tooth profile of planthopper hip gear is assumed to be a logarithmic helical.
(1) Binarization processing The original images of the gear teeth of the planthopper hip gear are acquired as depicted in Figure 5.
The adaptive thresholding binarization method 18 is used to binarize Figure 5.If the grey value of the pixel is greater than the adaptive threshold 128, the pixel is displayed as white; if it is less than the adaptive threshold 128, the pixel is displayed as black.The custom denoising method 18 is used to further suppress the noise  in the background of the binary image.As depicted in Figure 6(a)-(d), the binarized images of the planthopper hip gear with different denoising thresholds are obtained.From Figure 6(c), it can be seen that the edge contour of the planthopper hip gear is complete, with no noise points near the contour, which is selected for further tooth profile extraction.
(2) Extraction and segmentation of the planthopper hip gear profile The extracted marginal tooth profile of Figure 6(c) is depicted in Figure 7(a), which is not smooth.The marginal tooth profile is then smoothed 19 to obtain the smooth planthopper hip tooth profile as depicted in Figure 7(b).
Figure 8 shows the conceptual diagram of the tooth profile comparison between the planthopper hip gear and the involute gear. 5According to the real size of the planthopper hip gear, the tooth height of the planthopper hip gear is 30 mm.The dedendum of the planthopper hip gear is 16.67 mm and the addendum is 13.33 mm, referring to the distribution ratio of the dedendum and addendum of the standard involute gear.The tooth profile of the planthopper hip gear is segmented based on the involute gear meshing principle. 20Curves AB and EF are the root transition tooth profile, the curves BC and DE are the working tooth profile and curve CD is the top transition curve.The nodes of the planthopper gear are the points of intersection k1 and k2 of the nodal line with the working tooth profiles BC and DE.
According to the above method of segmentation and node determination, the real extracted tooth profile of the planthopper hip gear is segmented as depicted in Figure 9(a).To study the characteristics of the working tooth profile BC and DE of the planthopper hip gear, the working tooth profile BC and DE are extracted separately as depicted in Figure 9(b), and the working tooth profile BC and DE are referred to as the convex tooth profile and the concave tooth profile, respectively.
(3) Fitting the logarithmic spiral function to the working tooth profile of the planthopper hip gear.
The parametric equation of the logarithmic spiral is shown as equation ( 7) 21 : where r 0 is the starting circle radius; t is a constant; and u is the polar angle.
According to the principle of least squares, 22 it is transformed into a minimum value of the sum of the error squares, which is solved step by step using an iterative method.The coordinate transformation of the parametric equation of the logarithmic spiral in the Oxyz plane is given by equation ( 8).
where a d is the variation in the horizontal coordinate of the logarithmic spiral and b d is the variation in the vertical coordinate of the logarithmic spiral.
The extracted convex tooth profile BC is translated so that the end point B of the convex tooth profile coincides with the coordinate origin.Similarly, the extracted concave tooth profile DE is translated so that the end point D of the tooth profile coincides with the coordinate origin.The processed convex tooth profile BC and the concave tooth profile DE are fitted separately as depicted in Figure 10(a) and (b).The fitted logarithmic spiral parameters are depicted in Table 1, and the errors are depicted in Table 2.
The determination factor of the fitted convex and concave tooth profiles are both greater than 0.99, and the maximum relative error is within 0.25%, indicating that the logarithmic spiral fits the convex and concave tooth profiles with high accuracy.It also verifies the assumption that the working tooth profile of the planthopper hip gear is a logarithmic spiral.

Research on the bionic gear tooth profile design
Construction of logarithmic spiral gear tooth shape of bionic gear.The characteristics of the planthopper hip gear are: (i) the tooth profile is a concave-convex asymmetric tooth profile; (ii) the teeth of the meshing gears are concave-convex meshing, which achieves bidirectional transmission.According to the principle of gear meshing, 23 the conditions for meshing concave-convex meshing gear pairs are: (i) the meshing tooth profile must be a smooth curve; (ii) two tooth profiles mesh at one point, and the common normal of the contact point must pass through the node; (iii) the gear is a helical gear.
The bionic tooth profile of the planthopper hip gear is adjusted to comply with the meshing principle.The planthopper hip gear profile is depicted is Figure 11(a) by referring to the pitch position of the involute gear  profile.The center point of the convex tooth profile is O 1 , the pressure angle is u 0 , the center point of the concave tooth profile is O 3 , the pressure angle is u 3 , and the tooth top height and tooth root height of the gear are h a , h f , respectively.According to the principle of gear meshing, 20 the pressure angles on both sides of the tooth profile are guaranteed to be equal in size, that is, u 0 = u 3 , as depicted in Figure 11(b).The contact points k1 and k2 on both sides of the tooth profile are on the nodal line to satisfy the requirements of the common normal through the node.And the teeth are helical.To ensure that the tooth profile is simple and easy to machine, the adapted bionic tooth profile is further simplified, as depicted in Figure 11(c).The upper transition curve CD is a part of the upper circle of the tooth.The working tooth profiles BC and DE are logarithmic spirals, and smooth circular arcs are adopted for the root transition curves AB and EF, 24 as depicted in Figure 11(c).According to the pure rolling contact characteristics of the planthopper hip gear, the contact points on both sides of the tooth profile are limited to the gear nodal line in the design of the bionic gear to achieve theoretical pure rolling.
Design of logarithmic spiral tooth profile parameters for bionic gears.S s1 (P 1 2x s1 y s1 ) (depicted in Figure 12) is the coordinate system established with the center point of   the gear pitch thickness P 1 as the origin, the center line of the pitch tooth thickness as the y s1 axis, and the pitch line as the x s1 axis.S s2 (P 2 2x s2 y s2 ) is the coordinate system established where the center point of the gear tooth thickness is the origin P 2 , the tooth thickness center line is the y s2 axis and the pitch line is the x s2 axis.The angle range of the convex tooth profile of the logarithmic helix is (u 1 , u 2 ).Its pressure angle is u 0 .The center point is O 1 , which has the coordinates (w 1 , 2f 1 ) in the coordinate system S s1 (P 1 2x s1 y s1 ).The angular range of the concave tooth profile of the logarithmic spiral on the other side of the same gear tooth is (u 4 , u 5 ).Its pressure angle is u 3 .The center point is O 3 with coordinates (w1, 2f1) in the coordinate system S s2 (P 2 2x s2 y s2 ) with coordinates (w 2 , f 2 ).The pressure angles of the logarithmic spiral convex and concave tooth profiles are equal in size, that is, u 0 = u 3 .The vertical distance from the root transition point to the nodal line of the logarithmic spiral convex and concave tooth profiles are h g1 and h g2 .Its tooth top height and tooth root height are h a and h f .The tooth thickness of the gear is s.The side clearance is j, and the tooth slot width is s + j.The parameter formulas of the bionic gear tooth profile are derived as depicted in Table 3.
3D modeling of the bionic gear.The angle between the center line of the adjacent teeth of the planthopper gear is about 8.5°.If the teeth are distributed over the whole circle, the number of teeth is about 42.Therefore, the number of teeth in the bionic gear pair is 42.The size of the planthopper gear is enlarged by 100 times.The height of the bionic gear tooth is 3 mm.Referring to the pure rolling single arc gear, 24 the dedendum and addendum are 1.8 and 1.2 mm respectively, as depicted in Table 4.
According to the fitting results in Table 1, the starting circle radius of the convex tooth profile is r 1 = 4.2 mm, and the logarithmic helix coefficient is t l = 22.5 3 10 23 after the gear is amplified 100 times.To ensure that the gear does not interfere, the starting circle radius of the concave tooth profile is r 2 = 8.4 mm, u 1 = 6°and u 5 = 37°.Structural parameters of logarithmic helix tooth profile are depicted in Table 5.
The normal meshing tooth profile of the bionic gear is drawn and projected onto the end face to obtain the end meshing tooth profile of the bionic gear.The end meshing tooth profile is copied onto the other end face to create the tooth helix, 25 as depicted in Figure 13(a).The end face tooth profile is swept along the respective tooth helix to obtain the meshing tooth model of the bionic gear, as depicted in Figure 13(b).

Analysis of transmission efficiency and error of bionic gear pair
A geometric model of the bionic gear pair is imported into the Dynamics software.The rotating gear pair and the contact pair of meshing tooth surfaces are set up in a gravity-free working environment in ADAMS software.The contact stiffness coefficients are obtained from Hertzian theory, 16 as depicted in equation ( 9).
The contact stiffness of the bionic gear is K = 4.2246 3 10 5 N/mm.To prevent the shock caused by excessive starting torque, the input speed is changed from 0 to 1000 rpm in 0;0.2 s, and the speed is kept constant in 0.2;10 s with a load of 1 3 10 5 N.mm.The simulation time is 10 s with a step size of 0.003 s.

Concave profile working section
Concave tooth profile process angle Analysis of transmission efficiency of bionic gear pair.The input and output torque of the bionic gear pair are measured to calculate its transmission efficiency.Lowpass filtering is used to reduce the useless noise information in the input torque, as depicted in Figure 14(a).The output torque is the added load depicted in Figure 14(b).The calculated transmission efficiency is depicted in Figure 14(c)-(d).The transmission efficiency of the bionic gear pair is more than 97.5%.It indicates that this designed bionic gear pair is one of the high-efficiency transmissions. 16alysis of bionic gear transmission error.When the driving gear rotating at a constant speed, the displacement of driven gear is measured.The transmission error of the bionic gear is depicted in Figure 15 based on equation (2).The transmission error is within 25'' when the bionic gear rotates 360°.It indicates that the designed bionic gear pair is of high precision.

Contact strength analysis of bionic gear
According to Hertz contact theory, the contact stress of the bionic gear tooth surface is deduced to predict its contact strength.
(1) Theoretical calculation formula of maximum contact stress of bionic gear tooth surface When the moving point P changes with the curve's coordinate parameters u and v, its trajectory is a spatial surface.As depicted in Figure 16(a).The point P of the hypersurface has an infinite number of tangents.The  plane formed by all the tangents is called the tangent plane of point P on the hypersurface.The normal to point P on the surface is a line perpendicular to the tangent plane of the surface point, as depicted in Figure 16(b).
Referring to the theory of differential geometry, 26 the maximum value of normal curvature is the principal curvature of the tooth profile direction in the gear tooth surface, set as K 1 .The minimum value of normal curvature is the principal curvature of the tooth direction, set as K 2 .The two main directions in which the main curvature is located are perpendicular to each other. 27If the parametric equation of the surface is known, the mean curvature H of a certain point on the gear tooth surface can be obtained.The tooth surface equation of the bionic gear is uniformly expressed according to the indirect conjugation principle, as depicted in equation (10).
Assuming that the average curvature of the convex tooth surface of the bionic gear is positive, the pressure angle of the gear contact point is u 0 , and the pitch radius is R 1 .The average curvature of the contact point of the convex tooth surface is H (1) according to the differential geometry theory.Assuming that the average curvature of the concave tooth surface of the bionic gear is negative, the pressure angle of the gear contact point is u 0 = p + u 0 ,  and the pitch circle radius is R 2 .The mean curvature of the concave tooth surface contact point is H (2) .The relative mean curvature of the contact point of the bionic gear tooth surface is H n = H (1) + H (2) .It is known that the minimum of the relative principal curvature of the bionic gear tooth surface is zero and the maximum of that is 2H n .The radius of the relative principal curvature is the reciprocal of the maximum value of the relative principal curvature.Therefore, the relative principal curvature radius of the bionic gear meshing point can be obtained from equation (11).
According to the literature, 27 the basic formula for calculating the maximum contact stress of circular arc gear tooth surface by Hertz contact theory is as equation (12).
Where F t is the circumferential force on the pitch circle of the end face, u 0 is the pressure angle, L is the length of the contact curve, r n is the relative radius of the main curvature, E 1 and E 2 are the Young's Modulus of the driving and driven gear material respectively, v 1 and v 2 are the Poisson's ratio of the driving and driven gear material respectively.The maximum contact stress of the bionic gear tooth surface is deduced by combining equations ( 11) and ( 12), as depicted in equation (13).
The contact curve length L of the bionic gear is unknown, and the mathematical model of the contact curve length L is derived according to the finite element simulation results.
(2) Mathematical model of gear contact curve length based on finite element analysis Referring to [28], the bionic gear pair model is imported into ANSYS (the finite element analysis software).The material is 45Cr, Young's modulus is 209 GPa, and Poisson's ratio is 0.3.Hexahedral elements are used.The mesh of the contact tooth surface is encrypted separately, depicted in Figure 17.Fixed support is used to constrain the driven gear, and all freedom degrees of the driven gear are constrained.Considering the relative sliding between tooth surfaces, the friction contact type is selected.The friction coefficient between tooth surfaces is 0.15.The input torque of the driving gear is T 1 = 1 3 10 5 N.mm, depicted in Figure 18.
When the mesh size of the contact tooth surface is changed from 0.15 to 0.1 mm, the maximum contact stress of the tooth surface changes from 763.02 to 763.11 MPa in Table 6 with a change rate of 0.09%.Therefore, the contact tooth surface mesh is individually encrypted to 0.15 mm to ensure the accuracy of contact stress simulation.
Referring to [29], the contact stress at the center of the tooth width of the bionic gear is analyzed.The contact stress nephograms of the driving and driven gears are depicted in Figure 19(a) and (b) when contacting at the center of the tooth width, respectively.The contact area is elliptical.Stress concentration occurs at the geometric mutation (tooth root).The maximum contact stress of 762.95 MPa occurs on the concave tooth surface of the driven gear in the contact area, indicating that the simulation results reflect the mechanical characteristics of the working bionic gear pair.After the tooth surface gradually runs in, the curvature radius of the contact arc between the convex and concave tooth surfaces tends to be equal, achieving the ideal situation of instantaneous contact along the tooth height line, which reduces the contact stress concentration.
(3) Contact curve length L of bionic gear The gear contact curve length is mainly related to the spiral angle and module.Therefore, the control variable method is used to change the spiral angle and module of the bionic gear while keeping other parameters unchanged.A multi-group gear pair model is established, as depicted in Table 8.When the driven gear is fixed and the input torque is 1 3 10 5 N.mm of the driving gear, the maximum contact stress of the tooth surface is s H in equation (13).And the contact curve length L in equation ( 13) is derived.s H and L are depicted in Table 7.
The contact curve length L varying with the spiral angle and modulus is fitted by a polynomial as equation (14) shown.
The formula for calculating the maximum contact stress of the bionic gear tooth surface is obtained by combining equations ( 13) and ( 14).15) is deduced based on the geometric relationship.
The bending stress of the convex tooth root is greater than that of the concave tooth root, 31 so the maximum bending stress at the convex tooth root is selected for analysis.As depicted in Figure 21, the line passing through the center O and G points of the driving gear is the x-axis, the center of the gear circle is the coordinate origin, and the line perpendicular to the line OG is the y-axis.A G B G is the extension line of the line where the convex curve of the tooth profile transition arc curvature is located.The B G point is on the pitch circle of the gear.a is the included angle between OB G and the xaxis.The coordinate of A G is (x t , y t ).Since the included angle between the line m-m and the x-axis is small enough to be negligible, the line m-m and the x-axis are assumed to coincide.If y t 6 ¼ 0, the position of A G on the transition arc of the convex tooth profile is: Where a = (H cotg)/R, when y t = 0, x t = R2h f .As depicted in Figure 22(a), OO Z is the center line of the tooth thickness of the gear teeth.A G Z G is the tangent line of the A G point.g a is the included angle between the tooth thickness center line OO Z and A G Z G .
u a is the selected angle between OO Z and the x-axis.C G point is the center of the transition arc curve.p-p is the line through the B G point, tangent to the pitch circle.nn is the line where the direction of curvature is through the A G point.A perpendicular line is drawn through the B G point, perpendicular to the n-n line.Draw a perpendicular line through point C G , intersecting at point D G .The point of intersection point E G of the line OD G and the line n-n is the curvature center of the tooth root transition curve based on the Euler-Sauvard theorem. 31s depicted in Figure 22(b), S t is the vertical distance from A G point to OO Z , the center line of the tooth thickness.h F is the distance between two intersection points, A G normal and the line OO Z and k normal and the line OO Z .The load application point k coordinates are (x f , y f ).
Equation ( 17) is deduced based on the geometric relationship.
Where, S t = ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi x t 2 + y t 2 p sin u a .The load point coordinates of the convex tooth profile are: As depicted in Figure 22(b), h F is deduced as equation (19).
Where, S f = ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi x f 2 + y f 2 p sin u a .H a is deduced as equation (20). 30a = À cos g a r is the curvature radius of the stress point on the transition curve of the tooth profile root.If r = r g1 , the equation of the tooth profile bending stress s 0 is deduced as equation (21).
The maximum bending stress of the bionic gear tooth root occurs at the junction of the transition arc and the tooth root, which is the point A G ' depicted in Figure 20.Therefore, the calculation parameters of the maximum bending stress at the root of the bionic gear are depicted in Table 8.
The designed bionic gear is a kind of helical gear with an asymmetric tooth shape, which is quite different from an involute spur gear.Therefore, the correction value s C of tooth shape stress is added in equation (21).The mathematical model of the maximum bending stress of the bionic gear tooth root is deduced as equation (22) shown.
The mathematical expression of s C .Reference [31] shows that the spiral angle and modulus are important parameters of s C in the design of a single circular arc gear.The control variable method is adopted to select 20 groups of bionic gear pair models in Table 9 for finite element analysis of bending stress in ANSYS software.
When the driven gear is fixed and the input torque is 1 3 10 5 N.mm of the driving gear, the maximum simulation value of tooth root bending stress s F is subtracted from the theoretical stress value s 0 to obtain s C , as depicted in Table 9.The tooth shape stress correction value s C varying law with the spiral angle and modulus is polynomial fitted as equation (23).Fu et al.
The formula for calculating the maximum bending stress of the bionic gear tooth root is obtained by combining equations ( 22) and (23).

Strength comparison between bionic gear and conventional gear
Strength comparison between bionic gear and involute helical gear.The teeth number, tooth width, helix angle, pressure angle, and modulus of the selected involute helical gear match those of the bionic gear in section ''Research on the bionic gear tooth profile design.''For a more intuitive comparison of the strength of the involute helical gear and the bionic gear at different angular positions, the maximum contact stress and bending stress of the tooth root at different angles in the width direction are depicted in Figure 23(a) and (b), respectively.It is depicted in Figure 23 as follows.
(1) The involute helical gear is convex and convex profile meshing near the end face, while the bionic gear is convex and concave profile contact.The relative main curvature radius is different, so the maximum contact stress of the bionic gear tooth surface is small.The involute helical gear has a long contact line near the center of the tooth width, and the contact stress is distributed along the contact line.The bionic gear pair meshes at a point with a much smaller contact area than the line contact area.Therefore, the maximum contact stress of the involute helical gear tooth surface is smaller at the center of the tooth width.(2) Since the tooth height and root thickness of the bionic gear are smaller than those of the involute helical gear, the maximum bending stress of the bionic gear root is larger than that of the involute helical gear.
Strength comparison between bionic gear and pure rolling single arc gear.The teeth number, tooth thickness, tooth width, helix angle, pressure angle, and modulus of the pure rolling single arc gear are the same as those of the bionic gear in section ''Research on the bionic gear tooth profile design.''To avoid the tooth profile interference of the single arc gear, the working profile radius of the convex gear is 4.5 mm, 24 and the working profile radius of the concave gear is 12 mm.For a more intuitive comparison of the strength of the single arc gear and the bionic gear at different angle positions, the maximum contact stress and bending stress of tooth root at different angles in the width direction are depicted in Figure 24(a) and (b), respectively.
It is depicted in Figure 23 as follows.
(1) The difference of the tooth profile curvature radius of the meshing bionic gear pair is smaller than that of the single arc gear pair to avoid the meshing interference, so the maximum contact stress of the bionic gear tooth surface is smaller than that of the single arc gear.(2) The thickness of the bionic gear root is larger than that of the single arc gear, so the maximum bending stress of the bionic gear root is smaller than that of the single arc gear.

Conclusions
Based on the requirements of the transmission efficiency and load capacity of the jumping robot transmission system, this paper analyzes the transmission characteristics of the planthopper hip joint gear.A bionic gear is designed based on the tooth shape and profile characteristics of the planthopper hip joint gear.The logarithmic spiral is taken as the working tooth profile of the bionic gear.The design method of the bionic gear is proposed.The transmission efficiency and error of the gear are analyzed.In addition, the theoretical formulas for the contact strength and bending strength of the bionic gear are deduced.The results show that the bionic gear is more suitable for the transmission system of high-performance jumping robots.(1) The transmission efficiency of bionic gears is 97.5%.The bionic gear rotates 360°with a transmission error of 25".Its characteristics of small transmission error and high transmission efficiency meet the requirements of the jumping robot transmission system.
(2) The maximum contact stress of the bionic gear tooth surface decreases by 29.8% compared with the involute helical gear.The relative sliding coefficient of the meshing pair is close to 0, which reduces the tooth surface damage and improves the load capacity.The tooth height of the bionic gear is more than twice as small as that of the involute gear, and the design of the teeth and tooth slots is more compact, which is helpful for the lightweight of the jumping robot.(3) Compared with the single arc gear, the maximum contact stress of the gear tooth surface decreases by 3.5%, and the maximum bending stress of the tooth root decreases by 2.1%.Meanwhile, the bionic gear requires only one tool during processing, which improves the process versatility and reduces the manufacturing cost.

Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Figure 3 (
Figure 3(b).The maximum speed of the planthopper hip gear is 12000 r/min, indicating that the gear is used in high-speed transmission.12According to the literature,5 the maximum output torque T out of the planthopper hip gear drive is 0.7 mN.mm.The average weight M b of the planthopper is 6.6 mg.Its average body length L b is 4.1 mm, and the diameter D c of the planthopper hip gear is about 0.4 mm.To calculate the maximum output power P of the planthopper hip gear, the weight M of the planthopper hip gear pairs can be calculated in proportion to their length, as depicted in equation(1).

Figure 3 .
Figure 3. Planthopper jumping foot model and its speed: (a) Jumping foot model and (b) The angular speed of the planthopper hip gear.

Figure 7 .
Figure 7. Smoothing of the tooth profile: (a) Before smoothing and (b) After smoothing.

Figure 9 .
Figure 9. Extraction of the working tooth profile of the planthopper hip gear: (a) Segmented tooth profile and (b) The working tooth profile.

Figure 10 .
Figure 10.The fitted logarithmic spiral of the concave and convex tooth profiles: (a) The fitted convex tooth profile and (b) The fitted concave tooth profile.

Figure 11 .
Figure 11.Schematic diagram of the adjusted bionic tooth profile: (a) The planthopper hip gear profile, (b) The pressure angle on both sides and (c) The simplified tooth profile.

Figure 13 .
Figure 13.Meshing teeth model of bionic gear in the end face: (a)The tooth helix and (b)The meshing tooth model.

Figure 14 .
Figure 14.Transmission efficiency of bionic gear pair: (a) The input torque, (b) The output torque, (c) The transmission efficiency and (d) The transmission efficiency from 5.5 to 6s.

Figure 15 .
Figure 15.Transmission error of bionic gear pair.

Figure 16 .
Figure 16.Surfaces and their tangent planes and normals: (a) The spatial surface and (b) The tangent planes and normals.

Figure 20 .
Figure 20.Partial parameters of bionic gear tooth groove.

Figure 19 .
Figure 19.Equivalent stress nephogram of bionic gear: (a) The driving gear and (b) The driven gear.

Figure 22 .
Figure 22.Convex tooth profile root analysis: (a) Meshing state and (b) Non meshing state.

Figure 21 .
Figure 21.Tooth root position of bionic gear.

Figure 23 .
Figure 23.Comparison of contact and bending stress between bionic gear and involute helical gear at different angles in width direction: (a)The maximum contact stress and (b) The maximum bending stress.

Table 1 .
Parameters of the logarithmic spiral fitting equation.

Table 2 .
Fitting error of the working tooth profile and logarithmic spiral.

Table 3 .
Tooth profile and transition arc calculation formulae.

Table 4 .
Basic parameters of the bionic gear.

Table 5 .
Structural parameters of logarithmic helix tooth profile.
Radius of the starting circle of the logarithmic helix of the convex and concave tooth profile r 1 =r 2 (mm) 4.2/8.4Logarithmic spiral coefficient t l 22.5 3 10 23 Convex and concave tooth profile process angle u 1 =u 4 (°) 6 14.7789 Convex and concave tooth profile termination angle u 2 =u 5 (°) 46.6539 37 Logarithmic spiral tooth profile center lateral offset w 1 =w 2 (mm) 1.3912 5.0024 Longitudinal offset of logarithmic helical tooth profile center f 1 =f 2 (mm) 1.6479 3.2958

Table 6 .
Mesh size and corresponding tooth surface maximum contact stress.

Table 7 .
Different s H and L under different gear parameters.

Table 8 .
Calculation parameters of bionic gear tooth root bending stress.

Table 9 .
Bending stress analysis under different gear parameters.