Study on Dynamic Coupling Behavior of End-Meshing Harmonic Reducers

: To study the coupling mechanism and dynamic responses of an end-face-movable gear transmission system under complex excitation, a specific configuration of end-meshing, movable-gear reduction mechanism was used to achieve predetermined rigid-thrust-transmission and mismatched-gear-meshing functions, which solved the inherent defects of traditional harmonic gear mechanisms of thin-wall-flexible wheels that are easily damaged by fatigue. Considering the phenomenon of elastic deformation of live teeth that is accompanied by significant changes in meshing characteristics in the transmission process of an end-meshing harmonic reducer, the influences of dynamic meshing parameters, live tooth deformation, time-varying stiffness of tooth meshing, and time-varying backlash on nonlinear dynamic performance were explored, as well as the mechanisms of multi-parameter coupling effects on transmission performance. The nonlinear dynamics model of the end-meshing harmonic reducer was established to solve the chattering prediction problem. Finally, a comprehensive test bed for the transmission system of a harmonic reducer with a meshing type with an adjustable-characteristic end was built to verify the correctness of the theoretical model and provide the theoretical and technical basis for exploring the optimal parameter selection to address the passive vibration-suppression problem.


Introduction
For fields involving the precision operation of industrial robots and aerospace exploration activities, the important joint systems of its machines have significant requirements for the reliability and bearing capacity of harmonic gear transmissions [1].To mitigate the inherent shortcomings of the flexspline within traditional harmonic reducers, which are susceptible to fatigue damage, and to guarantee the precise operational state of the transmission system, an end-face engagement is employed to effectuate predetermined rigid-thrust transmission.This solution effectively reconciles the conflict between deformation and load.
The dynamic characteristics of harmonic gears are inherently complex due to the periodic changes in meshing parameters such as tooth deformation, meshing stiffness, and tooth backlash during the meshing process.Such changes cause obvious nonlinear coupling effects and mechanical jitter, causing continuous vibration in the transmission, which can have a detrimental effect on transmission performance.Noise and vibration not only affect the performance of machinery and equipment but can also affect the safety and stability of the system [2].Therefore, the study of the nonlinear dynamics of harmonic gear transmission systems can provide valuable engineering insights and a theoretical foundation for enhancing their performance.
Zhang et al. [3] first proposed the concept of an end-face harmonic gear drive with oscillating teeth, elucidated its transmission principle, and subjected its kinematics law attributes to comprehensive analysis.Subsequently, utilizing a virtual prototype model, they engineered and optimized the tooth profile to augment the transmission efficiency.The structural parameters pertinent to this drive system were calculated and optimized; however, the dynamic meshing behavior was not thoroughly explored.
Scholars [4][5][6][7][8][9][10] have established a nonlinear dynamic model of planetary gear sets based on considerations of the interaction of factors such as time-varying mesh stiffness and gear backlash.The influences of internal and external excitation on vibration and impact behavior was revealed, and measures to reduce vibration and impact were proposed.Zhu et al. [11] further extended application of the harmonic balance method to the nonlinear dynamic modeling of compound planetary gear sets.Through their research, they were able to study the nonlinear dynamic characteristics of gear sets and identify the influences of multiple factors, including dimensionless backlash, meshing stiffness, and error excitation amplitude, on frequency-response characteristics.Cui et al. [12][13][14] developed a coupled nonlinear dynamic model for compound planetary gears.This model was utilized to extract the natural frequency of the system and to evaluate the influences of moment of inertia, meshing stiffness, and other pertinent factors on its vibration response.Zheng et al. [15] proposed a translational-torsional coupled dynamics model of the RV reducer.This model was used to analyze the system's nonlinear time-varying behavior, yielding displacement responses of the various components under various conditions.Furthermore, the model was used to examine the sensitivity of the dynamic characteristics to both internal and external excitation factors.Tung et al. [16] developed a reduced-order, time-variant numerical model for the compound reducer that enabled the prediction of the robot's dynamic stiffness with enhanced accuracy.Hu et al. [17] developed a dynamic analysis model of the three-stage planetary transmission within a wind turbine reducer.This model involved utilization of the stiffness factor method to examine the mechanical properties of the interacting components.Subsequently, experimental validation was conducted to verify the reliability of the dynamic model.Yang et al. [18] and colleagues established a dynamic model of planetary gear trains utilizing the lumped parameter approach and conducted in-depth research on the vibration responses of planetary gear trains when subjected to both deterministic and random loads.In consideration of the effects of backlash and time-varying stiffness, Saeed Gheisari Hasnijeh et al. [19] have refined the conventional spur gear dynamics model by incorporating Gaussian white noise into the loading terms.Liu and Yang et al. [20,21] developed a nonlinear dynamic model of planetary gear transmission systems, which accounted for the effects of dynamic mesh parameters and stiffness.This model was subsequently subjected to experimental verification and analysis.
According to the literature review, current research on end-face harmonic gear drives predominantly involves kinematic analysis, with a limited focus on dynamic behavior.While numerous studies have been conducted on the system dynamics modeling of traditional harmonic and planetary gear reducers, it is evident that both internal and external excitations play significant roles in influencing the system's performance.Given the unique features of the end-face harmonic gear drive, the internal axial excitation has a significant effect on the dynamic behavior of the system, leading to significant buffeting phenomena that should not be overlooked.In contrast, traditional gear-system-dynamics modeling typically fails to analyze axial meshing parameters based on the radial meshing form.
For the reasons above, we constructed a prototype of the specific end-face harmonic gear drive of an oscillating-teeth reducer.According to previous analysis of the dynamic meshing-behavior law of the oscillating teeth, the influence of various dynamic meshing parameters, oscillating-teeth deformation, time-varying stiffness, time-varying backlash, and other coupling factors on the nonlinear dynamic performance were explored.Based on the lumped-mass method, the nonlinear dynamic model of the multi-tooth meshing system was developed.Then, the correctness of the theoretical model was verified through experiments.This study provides a theoretical basis for vibration-suppression analysis of the end-face harmonic gear transmission mechanism.

System Dynamic Incentive Analysis
In the process of gear transmission, the end-face harmonic gear drive of the oscillatingteeth reducer is subjected to a combined effect of internal and external excitations.Specifically, internal excitation primarily involves the elastic deformation of gear teeth under load, deformation of the support system derived from the assembly relationship, and the combined influence of time-varying meshing stiffness as a result of the intermittent engagement of teeth.External excitation, on the other hand, is primarily affected by the combined effects of power sources at the input end and the fluctuation of load and torque at the output end.As a result, the gear transmission system experiences vibration and noise due to the combined action of internal and external excitations.

Multi-Tooth Parameter Analysis of an End-Face Gear Drive
The principle behind an end-face harmonic gear drive of oscillating teeth [3] is shown in Figure 1a.When the end faces of the oscillating teeth in contact with the cam enter the cam lift profile, the oscillating teeth enter the meshing.The specific configuration of the end-face harmonic reducer is shown in Figure 1b.and other coupling factors on the nonlinear dynamic performance were explored.Based on the lumped-mass method, the nonlinear dynamic model of the multi-tooth meshing system was developed.Then, the correctness of the theoretical model was verified through experiments.This study provides a theoretical basis for vibration-suppression analysis of the end-face harmonic gear transmission mechanism.

System Dynamic Incentive Analysis
In the process of gear transmission, the end-face harmonic gear drive of the oscillating-teeth reducer is subjected to a combined effect of internal and external excitations.Specifically, internal excitation primarily involves the elastic deformation of gear teeth under load, deformation of the support system derived from the assembly relationship, and the combined influence of time-varying meshing stiffness as a result of the intermittent engagement of teeth.External excitation, on the other hand, is primarily affected by the combined effects of power sources at the input end and the fluctuation of load and torque at the output end.As a result, the gear transmission system experiences vibration and noise due to the combined action of internal and external excitations.

Multi-Tooth Parameter Analysis of an End-Face Gear Drive
The principle behind an end-face harmonic gear drive of oscillating teeth [3] is shown in Figure 1a.When the end faces of the oscillating teeth in contact with the cam enter the cam lift profile, the oscillating teeth enter the meshing.The specific configuration of the end-face harmonic reducer is shown in Figure 1b.The end-face cam rotates to push the oscillating teeth to move axially.The oscillating teeth always keep contact with the end face of the cam under the action of the spring, rise along the cam profile of the lift, and gradually engage with the gear teeth of the circular spline, which are forced to complete the rotation by the reaction force of the circular spline, and the power is output through the grooved pulley.This structure has both wheelbase and backlash adjustment functions.
To ensure optimally effective meshing between the oscillating teeth and the circular spline teeth, the tooth tops are appropriately shaped and trimmed to prevent potential interference during the meshing process.Due to the large number of oscillating teeth, analyzing each one individually would lead to redundant work, increasing the workload and analytical complexity.The method of tooth extraction allows for a clearer observation of the engagement conditions of individual teeth, thereby providing more precise data The end-face cam rotates to push the oscillating teeth to move axially.The oscillating teeth always keep contact with the end face of the cam under the action of the spring, rise along the cam profile of the lift, and gradually engage with the gear teeth of the circular spline, which are forced to complete the rotation by the reaction force of the circular spline, and the power is output through the grooved pulley.This structure has both wheelbase and backlash adjustment functions.
To ensure optimally effective meshing between the oscillating teeth and the circular spline teeth, the tooth tops are appropriately shaped and trimmed to prevent potential interference during the meshing process.Due to the large number of oscillating teeth, analyzing each one individually would lead to redundant work, increasing the workload and analytical complexity.The method of tooth extraction allows for a clearer observation of the engagement conditions of individual teeth, thereby providing more precise data and information.This facilitates a deeper understanding of the performance of the gear system.Furthermore, in the event of faults or anomalies within the gear system, tooth extraction enables timely detection of potential damage or wear, thus facilitating effective fault diagnosis.For the reasons above, the tooth extraction method is chosen, with each oscillating tooth being assigned three teeth to ensure optimal meshing strength and minimal motion vibration.The theoretical total number of oscillating teeth can be calculated by Equation (1): In the equation, N is the number of oscillating tooth blocks; Z A is the actual number of oscillating teeth; and Z V is the number of teeth removed from the oscillating gear.
The difference between the number of circular spline teeth Z E and the theoretical total number of oscillating teeth Z O represents the number of dislocated teeth required for each rotation of the cam, while the number of cam waves represents the number of reciprocating movements of the oscillating teeth in a cycle.To ensure the correct meshing of the gear teeth, the two values must be guaranteed to be equal.In this design, it is a single-stage harmonic drive, and the transmission ratio can be expressed as in Equation (2).
The cam-profile surfaces, which consist of two symmetrical helical surfaces, are integral to the determination of the axial displacement of the oscillating teeth.In each cycle, the displacement of the oscillating teeth situated on the right-hand and left-hand helical surfaces of the cam can be expressed as in Equation ( 3), where h is the cam lift, n is the cam revolution speed, and t is the cam running time.

Analysis of Meshing Stiffness between Circular Spline and Oscillating Teeth
An oscillating tooth's force condition following its entry into the meshing phase is given as depicted in Figure 2.
( ) In the equation,  is the number of oscillating tooth blocks;   is the actu of oscillating teeth; and   is the number of teeth removed from the oscillatin The difference between the number of circular spline teeth   and the total number of oscillating teeth   represents the number of dislocated teet for each rotation of the cam, while the number of cam waves represents the reciprocating movements of the oscillating teeth in a cycle.To ensure the correc of the gear teeth, the two values must be guaranteed to be equal.In this des single-stage harmonic drive, and the transmission ratio can be expressed as in (2).
The cam-profile surfaces, which consist of two symmetrical helical surface gral to the determination of the axial displacement of the oscillating teeth.In the displacement of the oscillating teeth situated on the right-hand and left-ha surfaces of the cam can be expressed as in Equation ( 3 where ℎ is the cam lift,  is the cam revolution speed, and  is the cam runn

Analysis of Meshing Stiffness between Circular Spline and Oscillating Teeth
An oscillating tooth's force condition following its entry into the meshin given as depicted in Figure 2. The meshing area and load point of the oscillating teeth are directly infl the displacement of the oscillating teeth.Given that the reaction force of the circ on the oscillating teeth is assumed to be an even load, the force generated b The meshing area and load point of the oscillating teeth are directly influenced by the displacement of the oscillating teeth.Given that the reaction force of the circular spline on the oscillating teeth is assumed to be an even load, the force generated by multiple teeth is equivalent to the concentrated load F of the center gear tooth.Based on the principles of force and torque equilibrium and input torque, the force equation of the oscillating teeth can be determined using Equation (4).
In the formula, L S is the length of the groove; L W is the distance between the force point at the bottom of the oscillating tooth and the bottom of the groove; and L M is the distance between the bottom of the oscillating tooth and the top of the tooth.
At the same time, the number of meshing teeth changes alternately with the cam angle.When the input torque is constant, the load on the oscillating teeth changes accordingly.The displacement of each oscillating tooth is denoted as z 1 ~z7 , respectively.The meshing relationship between it and the gear teeth of the circular spline is shown in Figure 3. teeth is equivalent to the concentrated load F of the center gear tooth.Based on the principles of force and torque equilibrium and input torque, the force equation of the oscillating teeth can be determined using Equation (4).
( ) In the formula,   is the length of the groove;   is the distance between the force point at the bottom of the oscillating tooth and the bottom of the groove; and   is the distance between the bottom of the oscillating tooth and the top of the tooth.
At the same time, the number of meshing teeth changes alternately with the cam angle.When the input torque is constant, the load on the oscillating teeth changes accordingly.The displacement of each oscillating tooth is denoted as  1 ~7 , respectively.The meshing relationship between it and the gear teeth of the circular spline is shown in Figure 3.The horizontal line in the figure represents the displacement required for the oscillating teeth to enter the meshing state.When the displacement of the oscillating teeth is higher than the horizontal line, the oscillating teeth enter the meshing.In the contrary situation, oscillating teeth are disengaged.That is, the number of meshing teeth at the intersection of the horizontal line and the fold line changes.
The oscillating gear teeth can be simplified to a tapered beam on the flexible body of the oscillating teeth, and the load  is equivalent to the right end of the microelement, which is decomposed into force   along the x-axis direction, force   along the y-axis direction, and the equivalent bending  generated by , as shown in Figure 4.
The effective length of the tapered beam is ℎ, i.e., the distance between the base point  of the oscillating tooth and the top of the tooth.The gear tooth is divided into a series of rectangular micro-elements along the x-axis direction from the bottom of the tooth to the load point.Each microelement is represented by the symbol , and its width is represented by the symbol   ; F is the normal load on the oscillating tooth;  is the angle between the load and the y-axis direction;   is the distance between the microelement along the x-axis direction and the load point;   is the distance between the load point along the x-axis direction and the tooth top;  is the distance between the load point and the x-axis; and   is the half-tooth width of the tooth root.The horizontal line in the figure represents the displacement required for the oscillating teeth to enter the meshing state.When the displacement of the oscillating teeth is higher than the horizontal line, the oscillating teeth enter the meshing.In the contrary situation, oscillating teeth are disengaged.That is, the number of meshing teeth at the intersection of the horizontal line and the fold line changes.
The oscillating gear teeth can be simplified to a tapered beam on the flexible body of the oscillating teeth, and the load F is equivalent to the right end of the microelement, which is decomposed into force F x along the x-axis direction, force F y along the y-axis direction, and the equivalent bending M generated by F, as shown in Figure 4.
The effective length of the tapered beam is h, i.e., the distance between the base point N of the oscillating tooth and the top of the tooth.The gear tooth is divided into a series of rectangular micro-elements along the x-axis direction from the bottom of the tooth to the load point.Each microelement is represented by the symbol i, and its width is represented by the symbol T i ; F is the normal load on the oscillating tooth; β is the angle between the load and the y-axis direction; L i is the distance between the microelement along the x-axis direction and the load point; L j is the distance between the load point along the x-axis direction and the tooth top; Y is the distance between the load point and the x-axis; and Y N is the half-tooth width of the tooth root.The calculation of deformation can be divided into three parts: the bending deformation, shear deformation, and axial compression deformation of the oscillating tooth body; the additional deformation caused by the elasticity of the oscillating tooth; and the local contact deformation of the tooth at the meshing point of the oscillating tooth.
(1) Calculation of Deformation of the Oscillating Tooth Body Under load, the oscillating gear teeth produce axial compression deformation, shear deformation, and bending deformation along the equivalent deformation of the load direction.The deformation of a single microelement can be calculated and superimposed.It is assumed that the left end of each microelement  is fixed, and that the portion connected to the right end of the microelement is regarded as a rigid body.The amounts of compression deformation, shear deformation, and bending deformation of the microelement can be obtained using Equation (5).
( ) In Equation ( 5),   ,   , and   are the meshing-point deformations caused by the compression, shear, and bending of the oscillating tooth body, and   ,   , and   are the deformations caused by the compression, shear, and bending of the microelement, respectively;  is the Poisson's ratio of the material; and   is the equivalent elastic modulus of the tooth.According to Cornell's analysis, the ratio of the tooth width   to the tooth thickness   is   /  < 5, so it is a narrow tooth.At this time, the value of   is the elastic modulus of the material;   is the cross-sectional area of the gear teeth;  1 is the deflection under the action of   ;  1 is the angle under the action of   ;  2 is the deflection under the action of ; and  2 is the angle under the action of , respectively, as shown in Equation (6).The calculation of deformation can be divided into three parts: the bending deformation, shear deformation, and axial compression deformation of the oscillating tooth body; the additional deformation caused by the elasticity of the oscillating tooth; and the local contact deformation of the tooth at the meshing point of the oscillating tooth.
(1) Calculation of Deformation of the Oscillating Tooth Body Under load, the oscillating gear teeth produce axial compression deformation, shear deformation, and bending deformation along the equivalent deformation of the load direction.The deformation of a single microelement can be calculated and superimposed.It is assumed that the left end of each microelement i is fixed, and that the portion connected to the right end of the microelement is regarded as a rigid body.The amounts of compression deformation, shear deformation, and bending deformation of the microelement can be obtained using Equation (5).
In Equation ( 5), δ a , δ s , and δ b are the meshing-point deformations caused by the compression, shear, and bending of the oscillating tooth body, and δ ai , δ si , and δ bi are the deformations caused by the compression, shear, and bending of the microelement, respectively; ν is the Poisson's ratio of the material; and E e is the equivalent elastic modulus of the tooth.According to Cornell's analysis, the ratio of the tooth width b i to the tooth thickness s p is b i /s p < 5, so it is a narrow tooth.At this time, the value of E e is the elastic modulus of the material; A i is the cross-sectional area of the gear teeth; ω 1 is the deflection under the action of F y ; θ 1 is the angle under the action of F y ; ω 2 is the deflection under the action of M; and θ 2 is the angle under the action of M, respectively, as shown in Equation (6).
In Equation ( 6), I i is the moment of inertia of the microelement.Substituting Equation (6) into Equation ( 5), results can be obtained using Equation (7).
(2) Calculation of Oscillating Tooth Body Deformation For narrow teeth, we deal with the problem of plane stress with Equation (8).
In Equation ( 8), E is the elastic modulus of the material; H f is the tooth thickness at the tooth root N; and L f is the equivalent arm of force. (

3) Calculation of Local Contact Deformation of Oscillating Teeth
The contact deformation of the meshing point of the oscillating tooth surface is caused by the contact and compression deformation of the gear meshing line, which can be expressed by Equation (9).
The total deformation of the meshing points of the movable teeth can be obtained by adding up each deformation as in Equation (10).
When the normal load F is constant and when the oscillating tooth displacement z > h t − h, the oscillating tooth and the circular spline tooth are meshing, as in Equation (11).

) Comprehensive Meshing Stiffness Calculation for Oscillating teeth
A pair of intermeshing gear teeth can be regarded as a pair of springs in series.The compressive stiffness, shear stiffness, bending stiffness, deformation stiffness, and Hertz contact stiffness of oscillating gear teeth are expressed as Equations ( 12)-( 16).
In summary, the meshing stiffness of a single pair of gears at the load point can be obtained using Equation (17).
In the formula, the meshing stiffness of the driving gear and the driven gear, namely the oscillating gear teeth and the circular spline teeth, at the load point is a function of the position of the meshing point.
In the equation, K 1 and K 2 represent the mesh stiffness of the driving and driven gears, respectively.They are functions of the meshing point position, specifically the contact point between the oscillating teeth and the circular spline teeth.
Setting φ W as the cam rotation angle, through the self-compiled software (Mathematica 12.1) for example analysis and drawing curves, the stiffness of each part of the first oscillating tooth can be obtained as shown in Figures 5-9.

1.275
In summary, the meshing stiffness of a single pair of gears at the load poin obtained using Equation (17).
In the formula, the meshing stiffness of the driving gear and the driven gear, the oscillating gear teeth and the circular spline teeth, at the load point is a functio position of the meshing point.
In the equation,  1 and  2 represent the mesh stiffness of the driving and gears, respectively.They are functions of the meshing point position, specifically tact point between the oscillating teeth and the circular spline teeth.
Setting   as the cam rotation angle, through the self-compiled software matica 12.1) for example analysis and drawing curves, the stiffness of each part of oscillating tooth can be obtained as shown in Figures 5-9.The comprehensive stiffness is shown in Figure 10.The comprehensive stiffness is shown in Figure 10.The comprehensive stiffness is shown in Figure 10.The remaining oscillating teeth have the same laws and only differ in phas ences.It can be obtained from the curve that the stiffness increases gradually w oscillating tooth gradually entering the meshing, and that it decreases gradually oscillating tooth gradually withdrawing from the meshing; the bending stiffnes stiffness, and base stiffness are positively correlated with load angle, while comp stiffness is the opposite.The comprehensive stiffness of the oscillating tooth mes improved by increasing .

Contact Stiffness Analysis
A universal ball bearing is used between the oscillating teeth and the cam contact point expands into a contact surface under the action of load N. The con face is projected onto the vertical surface of the contact normal, as shown in Figur the contact area between the universal ball and the circular spline, the contact distributed in the semi-ellipsoid, as shown in Figure 11b.The remaining oscillating teeth have the same laws and only differ in phase differences.It can be obtained from the curve that the stiffness increases gradually with the oscillating tooth gradually entering the meshing, and that it decreases gradually with the oscillating tooth gradually withdrawing from the meshing; the bending stiffness, shear stiffness, and base stiffness are positively correlated with load angle, while compression stiffness is the opposite.The comprehensive stiffness of the oscillating tooth mesh can be improved by increasing β.

Contact Stiffness Analysis
A universal ball bearing is used between the oscillating teeth and the cam, and its contact point expands into a contact surface under the action of load N. The contact surface is projected onto the vertical surface of the contact normal, as shown in Figure 11a.In the contact area between the universal ball and the circular spline, the contact stress is distributed in the semi-ellipsoid, as shown in Figure 11b.The contact between the cam and the universal ball is deformed as in Equation ( 18).
In Equation (18), N is the normal pressure at the contact point between the cam and the universal ball, and its direction is perpendicular to the cam profile; ′ is the sum of the equivalent elastic modulus;  is the main curvature; ma is the long half-axis coefficient of the contact ellipse; and () is the first type of complete elliptic integral related to the eccentricity of the ellipse e.The contact stiffness between the cam and the oscillating tooth can be obtained using Equation ( 19): The contact between the cam and the universal ball is deformed as in Equation (18).
In Equation (18), N is the normal pressure at the contact point between the cam and the universal ball, and its direction is perpendicular to the cam profile; E′ is the sum of the equivalent elastic modulus; Σρ is the main curvature; ma is the long half-axis coefficient of the contact ellipse; and K(e) is the first type of complete elliptic integral related to the eccentricity of the ellipse e.The contact stiffness between the cam and the oscillating tooth can be obtained using Equation (19): For the contact problem between the grooved pulley and the oscillating teeth, because the mass of the grooved pulley is significantly larger than the mass of the oscillating teeth the grooved pulley can be regarded as a rigid body.The oscillating tooth is squeezed by the groove pulley during the meshing process, and the expression of the compressive stiffness between the groove pulley and the oscillating tooth can be obtained using Equation ( 20): In Equation ( 20), b 3 is the width of the oscillating tooth rod; E 3 is the elastic modulus of the oscillating tooth material; and A 3 is the contact area between the oscillating tooth rod and the groove pulley.

Stiffness Analysis of Support System
Assuming that the radial stiffness of the bearing is isotropic since the input shaft is supported by a single bearing, the radial support stiffness of the cam is the radial support stiffness of the bearing at the input end, and its value can be calculated by Equation (21).
In Equation ( 21), F b1 is the radial load on the bearing at the input end; δ b11 is the radial elastic displacement of the bearing at the input end; δ b12 is the contact deformation between the outer ring of the bearing at the input end and the box hole; and δ b13 is the contact deformation between the inner ring of the bearing at the input end and the shaft diameter.The three expressions are expressed by Equations ( 22)-(24), respectively, where β b1 is the elastic displacement coefficient of the bearing at the input end, which is found from the standard according to the relative clearance g b1 /δ 01 ; δ 01 is the radial elastic displacement when the clearance in the bearing at the input end is zero; g b1 is the clearance or preload in the bearing at the input end; ∆ 1 is the fit clearance in the diameter direction between the outer ring of the bearing at the input end and the inner hole of the frame at the input end; H b11 is the elastic coefficient of the bearing at the input end; F b1 is the radial load; H b12 is the deformation coefficient; b b1 is the width of the bearing ring at the input end; and d b1 is the inner diameter of the bearing at the input end.
Similarly, the radial stiffness of the bearing at the output end can be obtained by Equation (25).
Because a pair of bearings are installed on the output shaft, the comprehensive radial support stiffness of the output terminal can be obtained by the method of spring parallel connection as shown in Equation (26).
The two shaft structures are both stepped shafts, and the torque provided by the motor or the load is received during the movement.The segmentation and torque are shown in Figure 12a,b.For the input shaft, the convex is simplified to a cylinder of equal width, ignoring the thin-walled support part between its inner diameter and the hub, and with processing in parallel with the matching shaft section, the total torsion angle of the input shaft can be obtained by Equation ( 27).
In Equation ( 27),  1 is the torsion angle of the i-th shaft section of the input shaft;  1 is the input torque; G is the shear modulus of the input shaft material;  11 ,  12 , and  13 are the lengths of each shaft section of the input shaft;  11 ,  12 , and  13 are the shaft diameters of each shaft section of the input shaft; and  14 and  15 are the inner and outer diameters of the cam, respectively.For the input shaft, the convex is simplified to a cylinder of equal width, ignoring the thin-walled support part between its inner diameter and the hub, and with processing in parallel with the matching shaft section, the total torsion angle of the input shaft can be obtained by Equation ( 27).
In Equation ( 27), φ 1i is the torsion angle of the i-th shaft section of the input shaft; T 1 is the input torque; G is the shear modulus of the input shaft material; l 11 , l 12 , and l 13 are the lengths of each shaft section of the input shaft; d 11 , d 12 , and d 13 are the shaft diameters of each shaft section of the input shaft; and d 14 and d 15 are the inner and outer diameters of the cam, respectively.
For the output shaft, the sheave is simplified into a cylinder of equal width, and its matching shaft section is regarded as a solid shaft as a whole, and the total torsion angle of the output shaft can be obtained by Equation (28).
In Equation (28), T 2 is the load torque; G is the shear modulus of the output shaft material; l 2i is the length of the i-th shaft section of the output shaft; and d 2i is the shaft diameter of the i-th shaft section of the output shaft.

Other Excitation Parameters
The mass and moment of inertia of each component are shown in Table 1.

System Dynamics Model
The following assumptions are made in the analysis: (1) During the motion process of the transmission system, all transmission components remain on the same axis.(2) Components such as the box are considered rigid bodies, and the contact between oscillating teeth and end-face cams, as well as that between the rigid circular spline and grooved pulley, is modeled as a spring-damping system.(3) Additional vibration phenomena arising from factors such as assembly and transmission errors among components are not considered.(4) The mass properties of each oscillating tooth are the same.
The relevant dynamic parameters that are set are shown in Table 2. Due to the dynamic behavior of each oscillating tooth being entirely identical, there is a phase difference in kinematics.Therefore, only the case of a single oscillating tooth needs to be considered.
The relevant space coordinate systems are established as follows: O − xyz is the fixed space coordinate system built on the rigid circular spline.O − x 1 y z 1 is the following space coordinate system rotating with the cam.O − x 2 y 2 z 2 is the following space coordinate system rotating with the grooved pulley.O − x 3 y 3 z 3 the following space coordinate system rotating with the i-th oscillating tooth.
The dynamic model of the transmission system is illustrated in Figure 13.The relevant parameters are listed in Table 3.
Table 3.The parameters relevant to Figure 13.Due to the dynamic behavior of each oscillating tooth being entirely identical, there is a phase difference in kinematics.Therefore, only the case of a single oscillating tooth needs to be considered.
The relevant space coordinate systems are established as follows:  −  is the fixed space coordinate system built on the rigid circular spline. −  1  1  1 is the following space coordinate system rotating with the cam. −  2  2  2 is the following space coordinate system rotating with the grooved pulley. −  3  3  3 the following space coordinate system rotating with the i-th oscillating tooth.
The dynamic model of the transmission system is illustrated in Figure 13.The relevant parameters are listed in Table 3.
Table 3.The parameters relevant to Figure 13.Taking the direction of the meshing point between the cam and the oscillating tooth toward the oscillating tooth as the positive direction, the projection of the cam displacement relative to the oscillating tooth along the direction of the mesh line can be obtained using Equation (29).

Symbol
Taking the direction of the meshing point between the grooved pulley and the oscillating tooth toward the oscillating tooth as the positive direction, the projection of the grooved pulley displacement relative to the oscillating tooth along the direction of the mesh line can be obtained using Equation (30).
Taking the direction of the meshing point between the rigid circular spline and the oscillating tooth toward the oscillating tooth as the positive direction, the projection of the rigid circular spline displacement relative to the oscillating tooth along the direction of the mesh line can be obtained using Equation (31).

System Dynamics Differential Equations
According to the system dynamics model, dynamics differential equations of transmission components can be obtained using Equations ( 32), ( 34) and (36).
For the cam, where T 1 is the input torque; x 13i , y 13i , and z 13i are the projections of δ 13i in the cam following coordinate system O − x 1 y 1 z 1 along the x 1 , y 1 , z 1 directions, respectively, which can be expressed as Equation (33).
For the grooved pulley, where T 2 is the output torque; x 23i , y 23i , and z 23i are the projections of δ 23i in the grooved pulley following coordinate system O − x 2 y 2 z 2 along the x 2 , y 2 , and z 2 directions, respectively, which can be expressed as Equation (35).
For the (i)-th oscillating tooth: where x 31i , y 31i , and z 31i are the projections of δ 31i in the (i)-th oscillating tooth following coordinate system O − x 3 y 3 z 3 along the coordinate axis; x 32i , y 32i , and z 32i are the projections of δ 32i in the (i)-th oscillating tooth following coordinate system O − x 3 y 3 z 3 along the coordinate axis; and x 34i , y 34i , and z 34i are the projections of δ 34i in the (i)-th oscillating tooth following coordinate system O − x 3 y 3 z 3 along the coordinate axis, respectively.The variables above can be expressed as Equation (37).
Variables f (x 34i ), f (y 34i ), and f (z 34i ) are the nonlinear functions of tooth side clearance along the x-axis, y-axis, and z-axis directions, respectively.Assuming that the tooth side clearance is 2b, the displacement of the oscillating tooth along the tooth surface in the normal direction is x n .The variables above can be expressed as (Equations ( 38), (39) and (40)), respectively.
From the equations above, the dynamics differential equations of the end-cam input mechanism, oscillating teeth, and grooved pulley output mechanism have been obtained.Organizing and sequencing the equations of each component and the overall dynamic differential equations of the system in matrix form is performed using Equation (41), where M is the system generalized quality matrix; X is the system generalized displacement matrix; C is the system generalized damping matrix; K is the system generalized rigidity matrix; and T is the system's externally-excited array matrix.

System Vibration Response Analysis
The nonlinear dynamics differential equations can be solved using the Runger-Kutta method.Setting speed and load of the cam in three different states and setting damping as Rayleigh damping, the transient vibration displacement response of the cam along each independent coordinate is obtained as shown in Figures 14-17 From the equations above, the dynamics differential equations of the end-cam input mechanism, oscillating teeth, and grooved pulley output mechanism have been obtained.Organizing and sequencing the equations of each component and the overall dynamic differential equations of the system in matrix form is performed using Equation (41), where  is the system generalized quality matrix;  is the system generalized displacement matrix;  is the system generalized damping matrix;  is the system generalized rigidity matrix; and  is the system's externally-excited array matrix.

System Vibration Response Analysis
The nonlinear dynamics differential equations can be solved using the Runger-Kutta method.Setting speed and load of the cam in three different states and setting damping as Rayleigh damping, the transient vibration displacement response of the cam along each independent coordinate is obtained as shown in Figures 14,15    From the equations above, the dynamics differential equations of the end-cam input mechanism, oscillating teeth, and grooved pulley output mechanism have been obtained.Organizing and sequencing the equations of each component and the overall dynamic differential equations of the system in matrix form is performed using Equation (41), where  is the system generalized quality matrix;  is the system generalized displacement matrix;  is the system generalized damping matrix;  is the system generalized rigidity matrix; and  is the system's externally-excited array matrix.

System Vibration Response Analysis
The nonlinear dynamics differential equations can be solved using the Runger-Kutta method.Setting speed and load of the cam in three different states and setting damping as Rayleigh damping, the transient vibration displacement response of the cam along each independent coordinate is obtained as shown in Figures 14,15     In Figures 14-17, the results show that the cam's transverse-displacement response is slightly larger than the radial displacement response.The reasons for this are that the cam-profile-surface lift angle and return angle are along the transverse direction, there is a fixed angle of contact between the cam and the oscillating teeth, and most of the vibration is transmitted to the transverse direction of the cam.
The transient vibration displacement response of the grooved pulley along each independent coordinate is shown in Figures 18,19    In Figures 14-17, the results show that the cam's transverse-displacement response is slightly larger than the radial displacement response.The reasons for this are that the cam-profile-surface lift angle and return angle are along the transverse direction, there is a fixed angle of contact between the cam and the oscillating teeth, and most of the vibration is transmitted to the transverse direction of the cam.
The transient vibration displacement response of the grooved pulley along each independent coordinate is shown in Figures 18,19   In Figures 14-17, the results show that the cam's transverse-displacement response is slightly larger than the radial displacement response.The reasons for this are that the cam-profile-surface lift angle and return angle are along the transverse direction, there is a fixed angle of contact between the cam and the oscillating teeth, and most of the vibration is transmitted to the transverse direction of the cam.
The transient vibration displacement response of the grooved pulley along each independent coordinate is shown in Figures 18-21, respectively.In Figures 14-17, the results show that the cam's transverse-displacement response is slightly larger than the radial displacement response.The reasons for this are that the cam-profile-surface lift angle and return angle are along the transverse direction, there is a fixed angle of contact between the cam and the oscillating teeth, and most of the vibration is transmitted to the transverse direction of the cam.
The transient vibration displacement response of the grooved pulley along each independent coordinate is shown in Figures 18,19      In Figures 18-21, the results show that the grooved pulley-transverse-displacement response is slightly less than those of the radial displacement response and axial displacement response.The reason for this is that while the oscillating teeth are meshing with the teeth of the rigid circular spline, there exists a central angle on the tooth side pointing to the center of rotation and a tooth angle on the tooth surface along the direction of the axis, which causes generation of the vibration of the oscillating gear to be more in the radial direction of the grooved pulley and to be transmitted to the grooved pulley.
The vibration patterns between pairs of dissimilar oscillating teeth are the same and there is only a phase difference.Taking the first pair of oscillating teeth as an example,   In Figures 18-21, the results show that the grooved pulley-transverse-displacement response is slightly less than those of the radial displacement response and axial displacement response.The reason for this is that while the oscillating teeth are meshing with the teeth of the rigid circular spline, there exists a central angle on the tooth side pointing to the center of rotation and a tooth angle on the tooth surface along the direction of the axis which causes generation of the vibration of the oscillating gear to be more in the radia direction of the grooved pulley and to be transmitted to the grooved pulley.
The vibration patterns between pairs of dissimilar oscillating teeth are the same and there is only a phase difference.Taking the first pair of oscillating teeth as an example   In Figures 18-21, the results show that the grooved pulley-transverse-displacement response is slightly less than those of the radial displacement response and axial displacement response.The reason for this is that while the oscillating teeth are meshing with the teeth of the rigid circular spline, there exists a central angle on the tooth side pointing to the center of rotation and a tooth angle on the tooth surface along the direction of the axis, which causes generation of the vibration of the oscillating gear to be more in the radial direction of the grooved pulley and to be transmitted to the grooved pulley.
The vibration patterns between pairs of dissimilar oscillating teeth are the same and there is only a phase difference.Taking the first pair of oscillating teeth as an example, In Figures 18-21, the results show that the grooved pulley-transverse-displacement response is slightly less than those of the radial displacement response and axial displacement response.The reason for this is that while the oscillating teeth are meshing with the teeth of the rigid circular spline, there exists a central angle on the tooth side pointing to the center of rotation and a tooth angle on the tooth surface along the direction of the axis, which causes generation of the vibration of the oscillating gear to be more in the radial direction of the grooved pulley and to be transmitted to the grooved pulley.In Figures 22-24, the results show that the transverse displacement of the oscillating teeth is the largest.The reason for this is that the setting of the angle between the oscillat ing teeth and the teeth of the rigid circular spline results in the line of action of the load being more directed in the transverse direction of the oscillating teeth.Therefore, the os cillating gear teeth are subjected to greater deformation and vibration responses th    In Figures 22-24, the results show that the transverse displacement of the oscillating teeth is the largest.The reason for this is that the setting of the angle between the oscillat ing teeth and the teeth of the rigid circular spline results in the line of action of the load being more directed in the transverse direction of the oscillating teeth.Therefore, the os cillating gear teeth are subjected to greater deformation and vibration responses along th    In Figures 22-24, the results show that the transverse displacement of the oscillatin teeth is the largest.The reason for this is that the setting of the angle between the oscillat ing teeth and the teeth of the rigid circular spline results in the line of action of the load being more directed in the transverse direction of the oscillating teeth.Therefore, the os cillating gear teeth are subjected to greater deformation and vibration responses along th transverse direction on entering meshing.In Figures 22-24, the results show that the transverse displacement of the oscillating teeth is the largest.The reason for this is that the setting of the angle between the oscillating teeth and the teeth of the rigid circular spline results in the line of action of the load being more directed in the transverse direction of the oscillating teeth.Therefore, the oscillating gear teeth are subjected to greater deformation and vibration responses along the transverse direction on entering meshing.
The cam, the grooved pulley, and the oscillating teeth all generate regular periodic vibrations in the transverse, radial, and axial directions at their equilibrium positions.When the oscillating teeth enter into and disengage from engagement, the vibration response changes abruptly.With an increase in the input speed of the single-wave cam and the load borne by the output end, the vibration amplitudes of the oscillating teeth in all directions increase gradually.

Dynamic Behavioral Validation
To validate the results of the nonlinear dynamics study of the end-face harmonic transmission system, we built the test bench shown in Figure 25.220V AC power is stepped down by a transformer to supply power to the stepper motor, torque sensor, and magnetic particle brake.The inputs control the stepper motors for power input through the driver, and the development board used is powered by 5V DC.The torque sensor is mounted on the input side of the end-meshing harmonic reducer to measure speed and torque, which are then read by the dynamic torsion-triple-display meter.The encoder is mounted to the output side to measure the real-time rotational speed of the output shaft, and the signal conversion is accomplished by the debugger.Magnetic particle brakes are mounted to the end to provide load and are adjusted in size by the load controller.The resulting sensing signals are transmitted to the PC for data processing.
Mathematics 2024, 12, x FOR PEER REVIEW 21 The cam, the grooved pulley, and the oscillating teeth all generate regular peri vibrations in the transverse, radial, and axial directions at their equilibrium posit When the oscillating teeth enter into and disengage from engagement, the vibration sponse changes abruptly.With an increase in the input speed of the single-wave cam the load borne by the output end, the vibration amplitudes of the oscillating teeth i directions increase gradually.

Dynamic Behavioral Validation
To validate the results of the nonlinear dynamics study of the end-face harm transmission system, we built the test bench shown in Figure 25.220V AC pow stepped down by a transformer to supply power to the stepper motor, torque sensor, magnetic particle brake.The inputs control the stepper motors for power input thro the driver, and the development board used is powered by 5V DC.The torque sens mounted on the input side of the end-meshing harmonic reducer to measure speed torque, which are then read by the dynamic torsion-triple-display meter.The encod mounted to the output side to measure the real-time rotational speed of the output s and the signal conversion is accomplished by the debugger.Magnetic particle brakes mounted to the end to provide load and are adjusted in size by the load controller.resulting sensing signals are transmitted to the PC for data processing.In Figures 26-30, the results show that the theoretical and experimental curves a general agreement in terms of trend and are numerically closer, which verifies the cor ness of the theoretical model.In addition to the influence of the load at the output end the rotational speed at the input end, the increase in the center distance makes the sup position of the rod part of the oscillating teeth gear system change, which leads to reduction in the meshing rigidity of the of the gear teeth and the increase in the vibra of the transmission system.According to the analysis of the meshing stiffness abov can be seen that the increase in pressure angle makes the meshing stiffness increase, w can reduce the vibration of the system.
The main reasons for the error between the experimental and simulation data fo vibration response are as follows: (1) Errors in the machining of components; (2) There is clearance between the rod part of the oscillating teeth gear system and moving part of the grooved pulley; (3) The box and bench support are not absolutely rigid; (4) The sensors are subject to a certain amount of deviation during the measurem   In Figures 26-30, the results show that the theoretical and experimental curves a general agreement in terms of trend and are numerically closer, which verifies the cor ness of the theoretical model.In addition to the influence of the load at the output end the rotational speed at the input end, the increase in the center distance makes the sup position of the rod part of the oscillating teeth gear system change, which leads to reduction in the meshing rigidity of the of the gear teeth and the increase in the vibra of the transmission system.According to the analysis of the meshing stiffness abov can be seen that the increase in pressure angle makes the meshing stiffness increase, w can reduce the vibration of the system.
The main reasons for the error between the experimental and simulation data fo vibration response are as follows: (1) Errors in the machining of components; (2) There is clearance between the rod part of the oscillating teeth gear system and moving part of the grooved pulley; (3) The box and bench support are not absolutely rigid; (4) The sensors are subject to a certain amount of deviation during the measurem process from the vibration of the test bench.As is shown in Figure 26, Case 1 is for the testing under 88 r/min, 0.65 N•m at the standard center distance.
As is shown in Figure 27, Case 2 is for the testing under 88 r/min, 0.98 N•m at the standard center distance.
As is shown in Figure 28, Case 3 is for the testing under 220 r/min, 0.65 N•m at the standard center distance.
As is shown in Figure 29, Case 4 is for the testing under 88 r/min, 0.65 N•m after reduction of the center distance by 0.5 mm.
As is shown in Figure 30, Case 5 is for the testing under 88 r/min, 0.65 N•m after reduction of the center distance by 1 mm.
In Figures 26-30, the results show that the theoretical and experimental curves are in general agreement in terms of trend and are numerically closer, which verifies the correctness of the theoretical model.In addition to the influence of the load at the output end and the rotational speed at the input end, the increase in the center distance makes the support position of the rod part of the oscillating teeth gear system change, which leads to the reduction in the meshing rigidity of the of the gear teeth and the increase in the vibration of the transmission system.According to the analysis of the meshing stiffness above, it can

Figure 3 .
Figure 3.The meshing relationship between the oscillating teeth and the circular spline gear teeth.

Figure 3 .
Figure 3.The meshing relationship between the oscillating teeth and the circular spline gear teeth.

Figure 4 .
Figure 4.The force model of the oscillating teeth.

Figure 4 .
Figure 4.The force model of the oscillating teeth.

Figure 10 .
Figure 10.Comprehensive stiffness of the oscillating tooth.

Figure 10 .
Figure 10.Comprehensive stiffness of the oscillating tooth.

Figure 11 .
Figure 11.(a) Contact between oscillating teeth and circular spline.(b) Contact stress between oscillating teeth and circular spline.

Figure 11 .
Figure 11.(a) Contact between oscillating teeth and circular spline.(b) Contact stress between oscillating teeth and circular spline.
of bearings are installed on the output shaft, the comprehensive radial support stiffness of the output terminal can be obtained by the method of spring parallel connection as shown in Equation (26).

Figure 12 .
Figure 12.(a) Segmentation of the input.(b) Segmentation of the output.
of the cam on the O − xyz due to vibration y 1 the longitudinal displacement of the cam on the O − xyz due to vibration z 1 the axial displacement of the cam on the O − xyz due to vibration φ 1 the torsion angle of the cam on the O − xyz due to vibration x 2 the transverse displacement of the grooved pulley on the O − x 1 y 1 z 1 due to vibration y 2 the longitudinal displacement of the grooved pulley on the O − x 1 y 1 z 1 due to vibration z 2 the axial displacement of the grooved pulley on the O − x 1 y 1 z 1 due to vibration Table 3. Cont.Symbol Explanation φ 2 the torsion angle of the grooved pulley on the O − x 1 y 1 z 1 due to vibration x 3i the transverse displacement of the (i)-th oscillating tooth on the O − x 3 y 3 z 3 due to vibration y 3i the longitudinal displacement of the (i)-th oscillating tooth on the O − x 3 y 3 z 3 due to vibration z 3i the axial displacement of the (i)-th oscillating tooth on the O − x 3 y 3 z 3 due to vibration θ W the spiral rise angle of the cam profile θ E the central angle corresponding to the rigid circular spline half-tooth r 1 the radius of the cam r 2 the turning radius of the grooved pulley φ W the cam rotation angle φ G the grooved pulley rotation angle τ i the relative rotation angle between O − x 3 y 3 z 3 and O − x 1 y 1 z 1 Mathematics 2024, 12, x FOR PEER REVIEW 14 of 25  2 the torsional stiffness of the grooved pulley support systems  2 the radial damping of the grooved pulley support systems  2 the axial damping of the grooved pulley support systems  2 the torsional damping of the grooved pulley support systems

Figure 13 .
Figure 13.The diagram of the transmission system dynamic model.

Explanation 𝑥 1
the transverse displacement of the cam on the  −  due to vibration  1 the longitudinal displacement of the cam on the  −  due to vibration  1 the axial displacement of the cam on the  −  due to vibration  1 the torsion angle of the cam on the  −  due to vibration  2 the transverse displacement of the grooved pulley on the  −  1  1  1 due to vibration  2 the longitudinal displacement of the grooved pulley on the  −  1  1  1 due to vibration  2 the axial displacement of the grooved pulley on the  −  1  1  1 due to vibration  2 the torsion angle of the grooved pulley on the  −  1  1  1 due to vibration  3 the transverse displacement of the (i)-th oscillating tooth on the  −  3  3  3 due to vibration  3 the longitudinal displacement of the (i)-th oscillating tooth on the  −  3  3  3 due to vibration  3 the axial displacement of the (i)-th oscillating tooth on the  −  3  3  3 due to vibration   the spiral rise angle of the cam profile   the central angle corresponding to the rigid circular spline half-tooth  1 the turning radius of the cam  2 the turning radius of the grooved pulley

Figure 13 .
Figure 13.The diagram of the transmission system dynamic model.
The vibration patterns between pairs of dissimilar oscillating teeth are the same and there is only a phase difference.Taking the first pair of oscillating teeth as an example, their transient vibration displacement responses along each independent coordinate are shown in Figures22-24, respectively.Mathematics 2024, 12, x FOR PEER REVIEW 20 of 2 their transient vibration displacement responses along each independent coordinate ar shown in Figures 22, 23 and 24, respectively.

Figure 25 .
Figure 25.The diagram of the test bench.

Figure 25 .
Figure 25.The diagram of the test bench.Measurement at standard center distance was performed under the following conditions: 88 r/min, 0.65 N•m; 88 r/min, 0.98 N•m; and 220 r/min, 0.65 N•m.After reduction of the center distance by 0.5 mm and 1 mm, respectively, the data under 88 r/min, 0.65 N•m was measured.The conditions above were recorded as Working Conditions (1-5).After data processing, the comparison of results between the output shaft-circumferential torsionalvibration response and the simulation results was performed as shown in Figures 26-30.

Figure 26 .
Figure 26.Comparison of simulation and experimental results for Case 1.

Figure 27 .
Figure 27.Comparison of simulation and experimental results for Case 2.

Figure 28 .
Figure 28.Comparison of simulation and experimental results for Case 3.

Figure 26 .
Figure 26.Comparison of simulation and experimental results for Case 1.

Figure 26 .
Figure 26.Comparison of simulation and experimental results for Case 1.

Figure 27 .
Figure 27.Comparison of simulation and experimental results for Case 2.

Figure 28 .
Figure 28.Comparison of simulation and experimental results for Case 3.

Figure 27 .
Figure 27.Comparison of simulation and experimental results for Case 2.

Figure 26 .
Figure 26.Comparison of simulation and experimental results for Case 1.

Figure 27 .
Figure 27.Comparison of simulation and experimental results for Case 2.

Figure 28 .
Figure 28.Comparison of simulation and experimental results for Case 3.

Figure 28 .
Figure 28.Comparison of simulation and experimental results for Case 3.

Figure 29 .
Figure 29.Comparison of simulation and experimental results for Case 4.

Figure 30 .
Figure 30.Comparison of simulation and experimental results for Case 5.

Figure 29 .
Figure 29.Comparison of simulation and experimental results for Case 4.

Figure 29 .
Figure 29.Comparison of simulation and experimental results for Case 4.

Figure 30 .
Figure 30.Comparison of simulation and experimental results for Case 5.

Figure 30 .
Figure 30.Comparison of simulation and experimental results for Case 5.

Table 1 .
The mass and moment of inertia of each component.

Table 2 .
The relevant dynamic parameters to be set.