Optimal Design of the Modular Joint Drive Train for Enhancing Cobot Load Capacity and Dynamic Performance

Automation advancements prompts the extensive integration of collaborative robot (cobot) across a range of industries. Compared to the commonly used design approach of increasing the payload-to-weight ratio of cobot to enhance load capacity, equal attention should be paid to the dynamic response characteristics of cobot during the design process to make the cobot more flexible. In this paper, a new method for designing the drive train parameters of cobot is proposed. Firstly, based on the analysis of factors influencing the load capacity and dynamic response characteristics, design criteria for both aspects are established for cobot with all optimization design criteria normalized within the design domain. Secondly, with the cobot in the horizontal pose, the motor design scheme is discretized and it takes the joint motor diameter and gearbox speed ratio as optimization design variables. Finally, all the discrete values of the optimization objectives are obtained through the enumeration method and the Pareto front is used to select the optimal solution through multi-objective optimization. Base on the cobot design method proposed in this paper, a six-axis cobot is designed and compared with the commercial cobot. The result shows that the load capacity of the designed cobot in this paper reaches 8.4 kg, surpassing the 5 kg load capacity commercial cobot which is used as a benchmark. The minimum resonance frequency of the joints is 42.70 Hz.


Introduction
The collaborative robot (cobot) represents a novel category of robots specifically designed to physically interact with humans in shared environments, eliminating the traditional barriers or protective cages found in conventional robotics systems [1][2][3].Possessing characteristics such as lightweight design, enhanced safety features, and simplified mechanisms for seamless human interaction, cobot find widespread applications in processing, assembly, surgery, and rehabilitation.Comprising lightweight linkages and modular joints, these components serve as the driving force for the cobot's movements.Modular joints, as the driving components of cobot, have been extensively studied by numerous scholars due to their ability to achieve the unique performance of cobot through their drive train design.
Unlike the high load capacity characteristic of industrial robots, collaborative robots have a relatively lower load capacity, which limits their scope of application.Reducing the weight of the cobot itself or increasing the load-to-self-mass ratio is an effective method to enhance the cobot's load capacity.In this regard, many scholars have conducted relevant work.An optimization objective that includes all the mass of motors and reducers is proposed by Zhou et al. [4][5][6], and the motors and reducers are selected by minimizing the objective.Yin et al. research the hybrid structure cobot links and the drive train regression design to achieve the lightweight of the cobot [7][8][9].With the development of cobot, more performance is concerned.A multi-objective optimization function is proposed to reduce the mass and expand the workspace of the cobot [10].The workspace and the mass of the cobot, including motors, gearboxes, and links, are combined with the scalarization method.Different from Ref. [10], Kawaharazuka et al. replace the workspace with manipulability, increasing flexibility while reducing the cobot mass [11].The motor's mass, the reducer's mass, and the reducer's cost are added to the objective function to design the drive train with the scalarization method of multi-objective [12].Padilla-Garcia et al. combine the total mass, tracking error, and the consumed energy to construct the optimal index, which is used to obtain comprehensive performance by configuring different motors [13].The tracking error, manipulability, and consumed energy are optimized by selecting the motors and reducers [14].
Recently, there has been an increased focus on dynamic characteristics, given that a high dynamic response is not only effective in suppressing vibrations but also crucial for enhancing cobot compliance.This enhancement makes cobot safer for interaction with the environment and facilitates easier direct teaching and cooperation [15].Du et al. conducted optimization on manipulator geometrical parameters and reducer selection to maximize fundamental frequency and payload-to-self-mass ratio [16].The design objective is a two-link system, where frequency improvement is considered by accounting for only the reducer's stiffness and the manipulator's weight.
The servo control variables of robots are generally the torque, position, or velocity of the motors.The dynamic performance of each joint is affected not only by output inertia and reducer stiffness but also by motor axis inertia.Hence, it is imperative to adopt a more comprehensive design approach for the joint drive train to enhance cobot dynamic performance.Additionally, the cobot's load capacity, a critical indicator, must be taken into account in the design process-a factor frequently overlooked by other designers.
In this study, optimization functions for both the dynamic response performance and load capacity of cobot are concurrently developed.The Pareto front is employed to select the optimal joint drive chain design, ensuring the overall performance of the cobot is maximized.To validate the effectiveness of the design method, a case study is conducted using a commercial collaborative cobot with a 5-kilogram payload as the benchmark.
The results indicate that this design method can simultaneously optimize the dynamic response performance and load capacity of cobot, achieving optimal performance.This paper is organized as follows: Section 2 analyzes the drive train structure, encompassing the cobot and the modular joint.The dynamics of the joint are also examined in this section.Section 3 covers the motor's design process, used for calculating the motor mass and rotor inertia.The optimization objective and design process are presented in Section 5, and a design example is implemented to illustrate the validity of the design method.Finally, Section 5 provides the conclusion of the work.

Cobot Structure Analysis
Figure 1 displays the commonly used types of cobot: sixaxis and seven-axis.In the seven-axis cobot, it is typical for the three joints near the base to form a spherical pair by inter-secting axes at a point, which is also observed in the terminal three joints, as depicted in Figure 1(a).
The six-axis cobot can generally be classified into two categories, without considering the bias between the joints.In the first category, joints 3 and 4 form the elbow of the cobot, with the axes of these two joints being vertical, as illustrated in Figure 1(b).The second category features one joint forming the elbow, while the fourth joint is positioned near the end effector, as shown in Figure 1(c).
The seven-axis cobot stands out for its high flexibility, making it suitable for obstacle avoidance.Conversely, the six-axis cobot exhibits higher stiffness and precision, making it more suitable for industrial applications [17].The fundamental units consist of modular joints and lightweight thin-walled links, with the modular joint serving as the actuator that significantly influences the dynamic performance of the cobot [18][19][20].Due to the typical structure and greater applicability of six-axis Figure 1 The configuration of cobot: (a) Configuration of the seven-axis cobot, (b) Six-axis cobot with two parallel axis configuration, (c) Six-axis cobot with three parallel axis configuration collaborative robots, the analysis and design are conducted focusing on the joints of the six axes.

Modular Joint Structure Analysis
Figure 2 depicts a typical modular joint, serving as the design example in this study.The components are conventional, with the exception of the torque sensor, enhancing torque measurement accuracy.Typically, a harmonic reducer is employed, capitalizing on its features of a substantial reduction ratio and low mass to enhance joint integration [21].
The torque sensor is positioned on the reducer's output port, and the motor rotor is connected to the reducer's wave-generator through a hollow shaft.A brake near the motor can secure the hollow shaft to bring the motor to a stop.The encoder gauges the position of the motor or the output flange of the joint.Moreover, two encoders are configured to measure position information from both the motor and output port.The driver oversees motor movement, transmitting information through the bus to the controller.

Modular Joint Dynamic Analysis
The cobot is constructed with modular joints and lightweight links in a serial configuration.The output torque of the joint can be computed using Eqs.( 1) and ( 2).This torque comprises inertial torque, Coriolis torque, gravity torque, damping torque, and external torque.The damping torque is the aggregate of viscous friction torque and Coulomb friction torque [22].
In the equations, τ ∈ R n represents the vector of joint output torque, M(q) ∈ R n×n is the symmetric and pos- itive-definite inertia matrix, C(q, q) q ∈ R n denotes the centripetal and Coriolis vector, g(q) ∈ R n is the gravity vector, τ f ∈ R n signifies the damping vector, τ ext ∈ R n is the torque vector generated by the external forces and (1) τ = M(q) q + C(q, q) q + g(q) + τ f + τ ext , (2) moments, D ∈ R n is the viscous friction torque coeffi- cient vector, and µ ∈ R n is the Coulomb friction torque coefficient vector.q ∈ R n , q ∈ R n , q ∈ R n represent position, velocity and acceleration vector of the joints.Superscript 'n' denotes the number of joints.The external forces and speeds can be transformed into joint space using the inverse Jacobian matrix in Eqs. ( 3) and (4): In the equations, J ∈ R n×6 represents the Jacobian matrix, F ∈ R 6×1 is the vector of external force and moment, and ẋ ∈ R 6×1 is the vector of endpoint speed.
The load and speed of the joint are crucial factors in selecting the harmonic reducer, ensuring compliance with the constraints specified in Eqs. ( 5) to (7).τ repre- sents the reducer output torque, and τ peak is the maxi- mum output torque of the reducer.q denotes the joint speed, and qpeak is the maximum output speed of the reducer.τ rrms stands for the maximum mean torque of the reducer, and the joint mean torque τ jrms must be less than this value.
The joint motor shaft rotational inertia J m,k is expressed in Eq. ( 8), where J rotor,k denotes the motor rotor rota- tional inertia, and J parts,k is the rotational inertia of vari- ous components, encompassing the wave-generator, CNC-machined aluminum alloy shaft, bearing rotation, and encoder rotation.Given the low density of the rotating component of the brake, its rotational inertia can be considered negligible.The variable subscript 'k' refers to the k-th joint starting from the cobot base, and the subsequent descriptions of variables follow the same convention.
The transmission mechanism of the harmonic reducer involves the wave-generator pressing against the gear of the flexspline on the circular spline.The flexspline should be designed to be thin and soft compared to the gearbox.This design imparts softness to the harmonic reducer, and its flexibility should be taken into account in the dynamic modeling and control of the modular joint.

Figure 2 Section view of a typical modular joint
Illustrated in Figure 3, τ m,k represents the motor torque applied by the stator to the rotor.b m,k is the damping coefficient of the motor shaft, influenced by the motor bearing and the harmonic reducer wave generator.K k denotes the stiffness of the joint, which is the reducer stiffness if the joint lacks a torque sensor.b l,k is the damp- ing coefficient on the joint output side, and J l,k is the inertia of the joint output side.
The mathematical expressions are presented in Eqs. ( 9) to (12).Here θ k represents the position of the motor, while θk and θk denote the acceleration and velocity of the motor, respectively.τ mo,k is the effective torque transmit- ted by the reducer from the motor, and N k is the reduc- tion ratio of the reducer τ l,k and τ ext,k represent the torque applied to the link by the reducer and the external force, respectively.q k signifies the position of the link, while qk and qk are the acceleration and velocity of the link, respectively.
Since the joint control primarily targets the motor, the transfer function from torque to motor velocity is derived, as presented in Eqs.(13) to (15).In this context, external forces act as disturbances, and the damp torque is negligible, allowing it to be disregarded.
The two-inertia system exhibits resonance frequency and anti-resonance frequency, depicted in Eqs. ( 16) and ( 17), especially when the damping coefficient is small.( 9) The resonance frequency indicates that the joint experiences intense vibration, while the anti-resonance frequency implies that the motor cannot effectively execute instructions from the driver.Consequently, the joint may become uncontrollable at these specific frequencies.
The transfer function exhibits zeros and poles at ω ARF ,k and ω NTF ,k .Given that a higher resonance frequency allows for greater control bandwidth, thereby enhancing dynamic performance [23], one of the design objectives for the cobot is to increase the resonance frequency.

Motor Design Process
Considering that the motor's mass impacts the inertia of other joints and the cobot's lifting capacity, and the shaft inertia influences the resonance frequency, the motor design process is thoroughly examined to identify the optimal design scheme.This section primarily focuses on the motor design process concerning mass and rotor inertia.

Constraints of Torque and Speed
The motor's rated torque is contingent on the reducer, with a larger reduction ratio requiring less torque.Eqs. ( 18) and ( 19) represent the peak load and the rated load, respectively, where η k is the reducer efficiency at 55 ℃ and 2000 r/min rated input speed.
The motor speed must adhere to the constraint outlined in Eq. ( 20) to ensure that the joint speed is adequate.

Constrains of Power Supply
The joint functions as a servo system configured with a driver, which propels the motor in accordance with commands from the controller.Since the driver operates in space vector pulse width modulation (SVPWM) mode [24,25] the motor voltage must adhere to Eq. ( 16)  (21).Here, V rated represents the bus voltage, the V lost is the lost voltage on the driver, and u ph and u phmax denote the phase voltage and the maximum phase voltage, respectively.
The current at the operating point must conform to Eq. ( 22), where i q,k represents the cross axis, and I rated,k is the root mean square (RMS) of the motor's rated current.

Calculation of Motor Mass and Rotor Inertia
Motor design encompasses various parameters, including both dimensional and electromagnetic aspects.In this paper, the design work primarily focuses on the motor's dimensions to calculate mass and inertia.Illustrated in Figure 4, the motor consists of the stator steel, armature, rotor, magnets, and yokes, all assembled together.The rotor's inner and outer diameters are denoted by d 1,k and d 2,k respectively, while the stator's inner and outer diameters are D 1,k and D 2,k .The param- eters t 1,k , t 2,k , t 3,k , t 4,k , t 5,k , t 6,k , they are teeth height, pole boot height, teeth width, slot depth, the yoke width of the stator, and the width of the notch, respectively.The magnet steel width is b mg,k , and its thickness is l m,k .The length of the air gap is denoted by l g,k , and the stator's thickness is l k , representing the depth in the vertical straight direction on the paper.The motor shaft is designed to be hollow to facilitate routing and reduce mass.The variable in motor design is the motor length, with other parameters changing correspondingly based on the motor length.(21) Eq. ( 23) can be used to calculate the width of the magnet b mg,k , where t s,k is the polar arc coefficient, and p k repre- sents the number of polar pairs.The electromagnetic parameters are computed as follows: where k l,k is the magnetic leakage factor, k c,k is the Carter factor, t se,k is the equivalent air gap length, B g,k is the air gap magnetic density, φ l,k and κ k are the intermediate variables, z k is the teeth number, B dw,k and B stl,k are mag- netic flux density of the stator steel and the rotor yoke, L dw,k and L stl,k are the length of the magnet path, µ rdw,k , µ rstl,k , and µ r are the relative permeability of stator steel, rotor yoke steel and the magnet steel.The length of the motor can be calculated using Eq. ( 30), and the area of the slot is given by Eq. (32): where W rated,k is the rated power of the motor, n rated,k is the motor rated speed, N ph,k is the phase turns number, k w,k is the winding factor, S slot,k is the single slot area, and (23)  Here, m mag,k represents the mass of a single mag- netic steel, m steelr,k is the rotor mass, m steel,k is the sta- tor steel mass, m copp,k is the mass of all the copper wire, and m all,k is the total motor mass.ρ steel and ρ copp are the density of the steel and the copper.d k is the diameter of the copper wire.The motor rotor inertia is then presented below:

Cobot Optimization Strategy
To enhance the cobot's load capacity and increase flexibility in human-cobot interactions, this study focuses on designing the modular joint drive train to optimize the comprehensive performance of load capacity and dynamics.Constructing a comprehensive performance evaluation index poses a challenge, given the varying load capacity and dynamic parameters with changes in cobot pose.The attributes of these two indicators are inconsistent, leading to controversy in constructing the evaluation index.
To address this challenge, the proposed design method in this work adopts a multi-objective optimization approach.Two optimization indices are constructed independently, and the Pareto method is employed to determine the optimal design parameters. (33) (37) Regarding the construction of optimization design indices, we formulate the indices when the robot is in a horizontal pose, which is shown in Figure 5, where the output force at the end effector of the robot is minimized.The optimization indices are outlined below.

Optimal Objective
Given the inter-joint load effects and the influence of the control algorithm on end effector dynamic characteristics, this study utilizes the minimum value of all joint resonance frequencies in the horizontal pose to formulate the dynamic performance function, as depicted in Eq. (41) Here, ω represents the vector of joint resonance frequencies which is shown in Eq. (40), and min(ω) is the minimum value within the vector.ω NTR,n is the resonance frequency of the nth joint, and n is the number of the joints.Additionally, ω max and ω min represent the maxi- mum and minimum values, respectively, of all min(ω) set values, serving as the boundaries of the value range in the design domain.As the resonance frequency contains information of anti-resonance frequency, and the dynamic performance of the robot is related to the inertia of the motor axis, this work utilizes resonance frequency to construct optimization indices.In this work, the end effector mass is also a crucial factor for optimization, calculated as the nominal load in the industry divided by gravitational acceleration.Therefore, Eq. (41) normalizes the resonance frequency, where a larger f d (x) implies superior dynamic performance for a cobot.Here, x rep- resents the vector of design variables, and ω is the func- tion of the vector x.
(40 The load capacity performance function, f l (x) is pre- sented in Eq. ( 42), where F uplift represents the uplift force at the end effector, also being a function of vector x.F max and F min denote the maximum and minimum values of F uplift within the design domain, respectively.Equation (42) normalizes F uplift , and a higher value of f l (x) indi- cates superior load capacity for a cobot.
Ensuring not only load capacity optimization but also the maximization of joint output torque is crucial.To achieve this, the joint motor's rated torque is calculated based on the reducer's maximum average load, and the rated speed is set accordingly.Employing a high torque density design approach for the motors reduces joint mass, thereby enhancing the cobot's load capacity [26,27].Consequently, for identical output torque and speed, the motor design can be determined by the air gap diameter [28] as the high torque density design approach of motor makes that only one variable is allowed between length and diameter.
In this study, the outer diameter of the motor rotor is selected as the optimal design variable.Hence, vector x is a function of each joint motor rotor outer diameter and the joint reduction ratios, as illustrated in Eq. (43):

Constraints
All the formulas from Eqs. (1) to (39) represent constraints, encompassing joint dynamics and the motor design process.Additional constraints are outlined in Eq. ( 44), where [d 2,k min , d 2,k max ] defines the range for the kth joint motor rotor outer diameter.{N 1,k , N 2,k , • • • , N y,k } signifies the reduction ratio vector of the kth joint, where y is the number of reduction ratios.M(q) (k,k) denotes the inertia of the cobot converted to the kth joint.K k repre- sents the stiffness of the kth joint, constituting the kth parameter of the stiffness vector K (k) .n is the total number of joints. (42) (44)

Optimization Strategy
The enumeration method is employed for optimizing the objective function, mitigating the risk of being trapped in local optima and simplifying the analysis of the design process.Given the computational intensity of motor design, discretizing the motor design scheme proves beneficial to enhance the optimization design speed.The specific optimization steps are as follows: Step 1: Establish the configuration and dimensions of the links.Additionally, identify the type of harmonic reducer and input its parameters, encompassing mass, wave generator rotational inertia, sets of reduction ratios, stiffness, rated load, and speed.
Step 2: Preliminarily determine the size of the joint structure.Design the motor considering the constraints of joint size, the rated load, and speed of the reducer, establishing the design range for the motor.Discretize the design range for the motor rotor's outer diameter.
Step 3: Employ the enumeration method to compute ω max and ω min , F min , and F max adhering to the constraints specified in Eq. (44).
Step 4: Calculate f d (x) and f l (x) using the enumeration method, and ascertain the optimal design scheme by analyzing the Pareto front.

Preliminary Design
The optimization design method is employed to design a six-axis cobot, as illustrated in Figure 5.The configuration of the six-axis cobot is depicted in Figure 1(b), with the fourth joint positioned close to the third joint to enhance the cobot's load capacity.Considering factors such as product appearance design, manufacturing, and comparison with commercial cobot, this cobot is equipped with only three types of joints.Specifically, joint 1, joint 3, and joint 5 are identical, while joint 2, joint 4, and joint 6 share similarities, with slight offsets introduced in joint 6 to simplify the machining process.Therefore, in the worst-case scenario pose, the optimization design of joint 2, joint 3, and joint 5 becomes critical in the cobot's design.
Building on the aforementioned analysis, the preliminary design and 3D modeling of the cobot are conducted, with a benchmark against a 5 kg load capacity industrial cobot.The density and dimensional parameters closely approximate those of the real prototype, ensuring that the physical parameters are accurate enough for the design optimization.
Joint 1 and joint 2 are fitted with the HSG-25 reducer model, while joint 3 and joint 4 utilize the HSG-20 reducer model.Joint 5 and joint 6 are equipped with the HSG-14 reducer model.
The initial parameters are presented in Table 1.Here, L 1 signifies the length of Link 1, while L 2 represents the length of Link 2. Additionally, L 3 indicates the distance from the axis of joint 5 to the center of mass of the end effector.L g1 and L g2 correspond to the distances from the axes of joints 1 and 3 to the centers of mass of Links 1 and 2, respectively.M is the mass of the end effector, serving as a reference for the 5 kg load capacity of a commercial cobot.Furthermore, M 1 and M 2 denote the mass values of Link 1 and Link 2, respectively.g is the gravitational acceleration.Lastly, M j,1 represents the mass value of the kth joint component, excluding the joint rotation shaft components, while J wg,k signifies the rotational inertia of the wave generator of the kth joint harmonic reducer.

Setup of the Boundary
The initial joint design establishes the parameters for motor design, as expressed in Eq. ( 45).These constraints en-compass the range of the motor rotor's inner diameter, the motor stator's outer diameter, the motor length, the reducer's rated torque, the motor's rated velocity, and the high-torque -density design approach for the motor.
The viable motor design space is investigated by adhering to the constraints outlined in Eq. ( 45), and the (45) resultant feasible design options are illustrated in Figure 6.The design of these motors emphasizes high torque density, entailing the imposition of various constraints throughout the design procedure.Key factors influencing the motor specifications include the rated torque, rated speed, and distinctive dimensions.
In this investigation, the chosen dimension variable is the motor rotor outer diameter ( d 2,k ), with the rated torque and rated speed being determined by the rated parameters of the reducer.Figure 6(a), (b) and (c) illustrate the motor design spaces for joint 2, corresponding to reducer reduction ratios of 80, 100, and 120.Likewise, (d), (e) and (f ) display the motor design spaces for joint 3, featuring reducer reduction ratios of 50, 80, and 100, while (g), (h) and (i) showcase the motor design spaces for joint 5, with reducer reduction ratios of 50, 80, and 100. Figure 6 graphically represents the motor mass and the rotational inertia of the motor rotor.
The findings from Figure 6 reveal a positive correlation between the motor rotor's rotational inertia and outer diameter, while motor length and mass decrease with an increase in rotor outer diameter when the output torque and speed are constants.Constraints for motor diameter are outlined in Eq. ( 46), defining the limitations for optimal functions.Specifically, the range for d 2,2 is from 44 mm to 58 mm, for d 2,3 it's from 32 mm to 42 mm, and for d 2,5 it's from 24 mm to 32 mm.Here, N 2 , N 3 , and N 5 denote the reducer's reduction ratios.
In this study, the stiffness value of the harmonic reducer is chosen from the middle section of the time varying stiffness, aligning with real industrial applications.The motor design process and the preliminary cobot design also serve as key components of the boundary conditions for the optimized design. (46)

Optimization Process
The discretization of the motor design feasible solution can significantly enhance calculation speed.To achieve this, the variable d 2,k is discretized with engineering precision.The discrete set for d 2,2 includes {44, 46, 48, 50, 52, 54, 56, 58}, d 2,3 includes {32, 34, 36, 38, 40, 42}, and d 2,5 includes {24, 26, 28, 30}.Consequently, there are a total of 5184 design solutions to be considered, accounting for the reduction ratio.The enumeration method is employed to identify the optimal solution, aiding in avoiding issues related to local optimality.The steps for this process are as follows.
Algorithm Enumeration

Optimal Solution
The calculations indicate that ω min is 27.55 Hz and ω max is 42.78 Hz occurring at the 4969th and the 204th steps, respectively.F min is 75.17N and F max is 82.56 N appear- ing at multiple steps, indicating the presence of several optimal solutions.All the data of f d (x) and f l (x) are plot- ted in Figures 7 and 8.The data presented in Figures 7 and 8 illustrate the presence of local optimal solutions for the function f d (x) and the periodicity of the function f l (x) .This observa- tion suggests that the enumeration method is suitable for finding the global optimal solution for the drivetrain design.The maximum values of functions f d (x) and f l (x) occur before the 216th step in the process, with a focus on sections A and B during the analysis.
In section B, the data shows a growth trend in the form of bands.However, in section A, the data exhibits an increasing trend, except for the last band, which shows a downward trend.Consequently, the maximum values of the functions f d (x) and f l (x) do not coincide at the same  In Figures 9 and 10, C1 and C2 correspond to the 204th step of the optimization process, representing the optimal values of bands A1 and B1.Similarly, D1 and D2 correspond to the 216th step, signifying the optimal values of bands A2 and B2.However, the condition that C1 is larger than C2 and D1 is less than D2 introduces uncertainty regarding the global optimal value.To address this, the Pareto front is plotted in Figure 11 using data from A1, A2, B1 and B2 to obtain the global optimal value.
The horizontal axis represents f d (x) , while the verti- cal axis represents f l (x) in Figure 11.The ideal solution is represented by point I, indicating that functions f d (x) and f l (x) simultaneously achieve their optimum values.Therefore, the distance from the function values to point I(1,1) serves as a metric to evaluate the optimality of each design scheme, with the optimum design achieved when this distance is minimized.Intuitively, the optimum is reached at either C(C1, C2) or D(D1, D2), and the distance is calculated using Eqs.( 47) and (48): As �f d < �f l , it indicates that point D is the global optimum, which is reached at the 216th step, and the corresponding design results are shown in Table 2.

Optimization Result Analysis
The design results indicate that the payload capacity of the cobot has reached the global maximum, and the optimization function for the dynamic performance of the cobot is also close to its maximum value.During the optimization process, due to the discrete nature of the speed ratio selection of the deceleration mechanism and the use of an enumeration method for optimization, the overall presentation shows discrete band characteristics, and there are local optimal solutions.In Figure 7, the three regions E, F, and G correspond to the speed ratios of joint 2 in the optimization process, namely 80, 100, and 120.The starting segment of region F is higher than the ending segment of region E, and similarly, the starting segment of region G is higher than the ending segment of region F. The starting segment of region E is identified as the globally optimal segment.The optimized results show that the reducer's reduction ratio selected for joint 2 is the largest in the set.The corresponding motor resulting from optimization has a relatively large length-to-diameter ratio, larger mass, and small rotor rotational inertia.For joints 3 and 5, the reducers' reduction ratio is the largest in the set.The motor selected through optimization has a small mass and a larger shaft rotational inertia.The motor design results are indicated by the red stars in Figure 6.
Figure 12 displays a distribution plot of min(ω), demonstrating a similar partition trend to the function f d (x) .However, its vertical axis represents the actual values of min(ω).The maximum value of min(ω) is 42.78 Hz, and the minimum value is 27.55 Hz.
The resonance frequency of each joint is crucial for the overall dynamic characteristics of cobot.Therefore, the distribution of resonance frequencies for the joint 2, 3, 5 of the cobot is studied in the design state during the optimization process, as shown in Figures 13,14,15.Figure 13 illustrates that the maximum resonance frequency of joint 2 is 42.78 Hz at the 216th step, while the minimum resonance frequency is 27.55 Hz at the 4969th step.Moving to Figure 14, it presents the data of joint 3 resonance frequency, with the maximum value reaching 78.09Hz at the first point of each band and the minimum value 42.60 Hz at the 205th step.The 216th step value is 42.70 Hz which is the optimization result.Similarly, Figure 15 provides data for joint 5 resonance frequency, where the maximum value is 84.28 Hz at the first point of each band and the minimum value is 64.24Hz at the last point of each band.Consequently, the last band in section H of Figure 12 represents data for joint 3, while the remaining data corresponds to joint 2.
Analyzing diagrams 13, 14, and 15, we can observe that the resonance frequency of joint 3 is almost entirely higher than that of joint 2. Additionally, for joint 5, its resonance frequency is mostly higher than that of joint 3. Therefore, for cobot, joints closer to the base have the most significant impact on the overall dynamic performance of the cobot.The optimal joint drive train scheme is achieved at step 216th, where the resonance frequency of joint 2 is at its maximum, while the resonance frequencies of joints 3 and 5 are at their minimum.Figure 16 is the value of the cobot lifting force.It can be observed that the cobot's maximum lifting force is 82.56 N, and the minimum lifting force is 75.17N. The maximum value is achieved at the 216th step of the optimization process.The maximum value is 10% greater than the minimum value.
The optimized results indicate that successfully improving the cobot's load capacity is achieved by reducing the mass of joints 3 to 6.This reduction enables the joints to exert greater external force while simultaneously increasing the resonance frequencies of joints 2 and 1, thereby enhancing the cobot's dynamic response performance.
Considering that the resonance frequency distribution band of joint 2 is the lowest among joints 2, 3 and 5, the optimization design process selects a scheme for joint 2 with the maximum resonance frequency.However, this choice involves a trade-off, as it sacrifices the optimal resonance frequencies of joints 3 and 5 to enhance the cobot's load capacity.Based on the optimized design results, a cobot prototype is produced and shown in Figure 17.
To effectively demonstrate the effect of the optimized design, this study compares its results with the parameters of commercial cobot, as listed in Table 3.
Due to the unavailability of dynamic characteristic parameters for commercial cobot, the comparison in this study is limited to load capacity, self-weight and working radius parameters, which are used to illustrate the cobot load performance.Given that the design case in this study aligns with commercial 5 kg payload cobot, a comparison is conducted with UR5e [29] and AUBO-i5 [30].From Table 3, it can be seen that the design case has the maximum payload, and the largest working radius.Therefore, the cobot's load capacity in the design case is superior, indicating the effectiveness of the design.

Conclusions
The lightweight structural design results in a relatively lower load capacity for the cobot.Simultaneously, enhancing the dynamic response characteristics of cobot is crucial for improving human-robot collaboration performance.Therefore, this study focuses on optimizing the overall performance of the cobot in terms of load capacity and dynamic response characteristics through the design of the joint transmission system.The conclusions are summarized as follows: (1) Using joint resonance frequency and cobot load capacity to construct optimization indicators, and applying the Pareto frontier to optimize the joint transmission system, can achieve the optimal comprehensive performance of cobot's dynamic response and load capacity.(2) The design process reveals that joints near the base have the lowest resonance frequency design domain.Therefore, it is necessary to concentrate on  designing the drive train for these joints to enhance the cobot's dynamic response characteristics.The design results indicate that joints near the base are best suited for the adoption of harmonic reducers with a small reduction ratio and motors with low rotor inertia, which helps improve the cobot's dynamic response characteristics.(3) For joints far from the base, their drive train need to adopt harmonic reducers with a large reduction ratio and lightweight motors to enhance the cobot's load capacity.(4) Compared to commercial cobot, the cobot in this design case has higher load capacity.
Based on the analysis above, the method proposed in this work can enable cobot to achieve optimal integration of load and dynamic response characteristics.This method is versatile, applicable not only to six-axis collaborative robots but also suitable for designing seven-axis collaborative robots and other electromechanical devices

Figure 3
Figure 3 Dynamics model diagram of modular joint

Figure 4
Figure 4 Motor dimension design illustration

Figure 5
Figure 5 Six-axis cobot of the optimization case

Figure 6
Figure 6 Motor design domains of the joint 2, 3 and 5: (a) The design spaces of joint 2 with N=80, (b) The design spaces of joint 2 with N=100, (c) The design spaces of joint 2 with N=120, (d) The design spaces of joint 3 with N=50, (e) The design spaces of joint 3 with N=80, (f) The design spaces of joint 3 with N=100, (g) The design spaces of joint 5 with N=50, (h) The design spaces of joint 5 with N=80, (i) The design spaces of joint 5 with N=100

Figure 7 Figure 8
Figure 7 Values distribution of the objective f d (x) in the optimization process .

Figure 9 Figure 10 Figure 11
Figure 9The last two bands of the section A

Figure 12 Figure 13 Figure 14 Figure 15
Figure 12 Values distribution of min(ω) in the optimization process

Figure 16 Figure 17
Figure 16 Values distribution of F uplift in the optimization process

Table 1
Parameters of the six-axis cobot

Table 2
Result of the optimization

Table 3
Parameter Comparison table of relevant benchmark commercial robots