Research on Kinematic Calibration and Trajectory Tracking of the Dual-Robot Collaborative Grinding and Polishing System

1. Introduction

With the rapid development of modern manufacturing technology, three-dimensional curved surface workpieces have been increasingly and widely applied in such fields as automobiles, aerospace, additive manufacturing, and abrasive tool manufacturing. As a key link to achieve high-precision surface processing, the grinding and polishing technology has put forward higher requirements for production efficiency and product quality. For a long time, in the majority of situations, this type of work has relied on manual operation. However, owing to the adverse working environment, the high-intensity workload, and the escalating precision requirements, not only has it led to substantial labor costs, but quality problems have also frequently arisen. Consequently, it is an inescapable trend that automated [1,2,3] grinding technology is progressively replacing traditional manual work.

In the current domain of grinding and polishing, the application of robots predominantly centers around single-robot grinding [4,5]. Nevertheless, as industrial products are becoming more and more complex and diversified, and the forms of complex machined parts are constantly changing, the single-robot grinding system encounters a great many challenges in respect of the grinding coverage rate. Owing to the constraints on its motion and operation, it is arduous to comprehensively and efficiently cover all the sections of complex workpieces requiring grinding. The dual-robot collaborative technology offers a novel solution to this problem. In this technical solution, the system employs a dual-robot collaboration model. The master robot secures the workpiece to be processed via the flange, while the slave robot holds the grinding tool through the flange and executes the grinding and polishing operations. Specifically, the robot tasked with fixing the workpiece possesses the kinematic ability of six degrees of freedom. It is capable of dynamically adjusting its pose during the grinding operation in accordance with the demands of diverse grinding tasks and the intricate shapes of the workpieces, thereby attaining a higher level of precision and flexibility. This design is capable of ensuring that the grinding area is comprehensively covered, effectively resolving the technical problem that a single robot usually has difficulties in achieving full coverage when handling workpieces with complex shapes. Currently, in comparison with the research accomplishments in numerous fields both domestically and internationally, the studies regarding the collaborative motion of dual robots [6,7,8,9] and the calibration of dual robots remain relatively scarce. Qin et al. [10] proposed a calibration method combining dual quaternion closed-form solutions with an iterative Levenberg-Marquardt (LM) algorithm-based approach for simultaneously calibrating hand-eye, flange-tool, and inter-robot relationships in dual-robot collaborative systems. Wang et al. [11] addressed the kinematic calibration problem of robotic hand-eye systems by proposing a separable calibration framework based on dual quaternions. Through theoretical derivation and experimental validation, they conducted an in-depth analysis of the comparative advantages, limitations, and applicable conditions between simultaneous and separable calibration methods. Hua et al. [12] proposed a hand-eye calibration algorithm integrating convex relaxation-based global optimization with dual quaternions, which enhances robotic hand-eye coordination accuracy to meet the high-precision positioning requirements for switchgear operations.

In this study, the master-slave control approach is employed. The EPSON C4-A901S robot is designated as the master robot for gripping the workpiece, while the EPSON C4-A601S robot is selected as the slave robot tasked with the grinding operation. Kinematic constraints are first established based on the closed kinematic chain formed by the dual-robot system. The “handshake” method [13] and the seven-point method [14] are then employed to precisely define the spatial transformation between the base and tool coordinate systems of the dual robots, enabling system modeling in MATLAB 2022a. Subsequently, the master and slave robots are controlled to follow predefined trajectories, with comprehensive simulation analysis conducted to validate feasibility. Finally, an experimental platform integrating EPSON C4 series robots, a grinding module, a clamping unit, and a drive system is constructed. Coordinated motion experiments are performed on the dual-robot collaborative grinding system, verifying both the reliability of the proposed approach and the accuracy of the simulation model.

2. Kinematic Analysis of Dual-Robot Systems

In this study, an EPSON C4 series robotic system is employed, where the EPSON C4-A601S—selected for its high joint motion velocity and positioning accuracy—is equipped with a grinding module to function as the grinding robot, while the EPSON C4-A901 serves as the workpiece-holding robot. Figure 1 illustrates the coordinate system relationship within the dual-robot grinding and polishing system. In the figure, [R_i] denotes the robot base coordinate system, [E_i] represents the end-effector coordinate system of the robotic arm, and [T_i] indicates the tool coordinate system of the robotic arm, where i=1 corresponds to the EPSON C4-A901S robot and i=2 refers to the EPSON C4-A601S robot, while [w] signifies the workpiece coordinate system. When the grinding tool establishes contact with the workpiece, the dual-robot system forms a closed kinematic chain. Based on the master-slave control strategy for system decoupling, the constraint equations of the system can be derived as follows:

$${}^{R1}T_w={}^{R1}T_{E1}\cdot {}^{E1}T_{T1}\cdot {}^{T1}T_w$$

(1)

$${}^{R1}T_w={}^{R1}T_{R2}\cdot {}^{R2}T_{E2}\cdot {}^{E2}T_{T2}\cdot {}^{T1}T_w$$

(2)

Wherein: ^R₁T_w represents the transformation matrix from the workpiece coordinate system to the master robot’s base coordinate system. Upon specification of grinding task requirements and determination of the tool-path trajectory, this transformation matrix becomes fully defined and serves as a known parameter in the system; ^E1T_T1 and ^E2T_T2 denote the transformation matrices from the gripping tool to the master robot and from the grinding tool to the slave robot, respectively. While both matrices are structurally defined by the mechanical assembly, ^E2T_T2 must be calibrated owing to the unique properties of the grinding tool, resulting in a constant matrix post-calibration. The homogeneous transformation matrix ^R1T_R2 defines the spatial relationship between the base frames of both robotic systems, which is determined through calibration procedures and remains invariant during operation. The homogeneous transformation matrix ^T1T_w defines the spatial relationship between the workpiece frame and the gripping tool frame. This matrix is structurally determined by the mechanical assembly configuration, where the workpiece maintains a fixed pose relative to the gripping tool, resulting in a time-invariant transformation. The homogeneous transformation matrix ^T2T_w defines the dynamic spatial relationship between the workpiece frame and the grinding device frame. As the grinding tool must maintain persistent contact with the target point throughout the process, this matrix becomes path-dependent and exhibits time-varying characteristics [15].

2.1. Dual-Robot Base Frame Calibration

Kinematic analysis of the dual-robot system indicates that calibration of ^R1T_R2 is essential for establishing accurate inter-robot spatial relationships. Since this relative pose relationship serves as the fundamental basis for dual-robot collaborative operation, its precise measurement is imperative. The calibration of dual-robot base frames aims to determine the transformation matrix between coordinate system R1 and R2. This study employs a calibration method in which two EPSON robots are each equipped with a calibration needle at their end-effectors. The calibration procedure involves controlling the robots to bring their needle tips into contact at four spatially distinct positions, with the critical condition that these four contact points must not lie on the same plane. Based on the coordinate relationships obtained from these non-coplanar contact points, the calibration matrix between the two robot systems is then mathematically determined. This experimental method enables direct determination of the coordinate transformation relationship between the two robotic base frames without requiring auxiliary measurement tools.

The calibration procedure is conducted as follows:

(1)

Coordinate system definition and notation specification: The world coordinate system, base frames (R1 and R2 for Robot 1 and Robot 2 respectively), and tool frames (T1, T2) are explicitly defined. With calibration needles mounted on each robot’s end-effector, both manipulators are controlled to establish four non-coplanar contact points (P₁-P₄) within their shared workspace. The corresponding coordinates of each contact point are recorded in both robot base frames as ^R1P_k and ^R2P_k (k=1,2,3,4).

(2)

Rotation matrix computation: Assuming the world coordinate system coincides with Robot 1’s base frame, we obtain:

$${}^{bi}P={}^{R1}R_{R2}\cdot {}^{R2}P+{}^{R1}T_{R2}$$

(3)

Substituting the coordinates of the four contact points into the above equations, subtracting point coordinates and reformulating them in matrix form yields:

$$\left[\begin{array}{ccc}{}^{R1}x_1-{}^{R1}x_2& {}^{R1}x_1-{}^{R1}x_3& {}^{R1}x_1-{}^{R1}x_4\\ {}^{R1}y_1-{}^{R1}y_2& {}^{R1}y_1-{}^{R1}y_3& {}^{R1}y_1-{}^{R1}y_4\\ {}^{R1}z_1-{}^{R1}z_2& {}^{R1}z_1-{}^{R1}z_3& {}^{R1}z_1-{}^{R1}z_4\end{array}\right]={}^{R1}R_{R2}\cdot \left[\begin{array}{ccc}{}^{R2}x_1-{}^{R2}x_2& {}^{R2}x_1-{}^{R2}x_3& {}^{R2}x_1-{}^{R2}x_4\\ {}^{R2}y_1-{}^{R2}y_2& {}^{R2}y_1-{}^{R2}y_3& {}^{R2}y_1-{}^{R2}y_4\\ {}^{R2}z_1-{}^{R2}z_2& {}^{R2}z_1-{}^{R2}z_3& {}^{R2}z_1-{}^{R2}z_4\end{array}\right]$$

(4)

To ensure the solvability of ^R1R_R2, the selected four points must be non-coplanar (the determinant of the matrix on the right-hand side of the equation must be non-zero), from which the rotation matrix is derived as follows:

$${}^{R1}R_{R2}=\left[\begin{array}{ccc}{}^{R1}x_1-{}^{R1}x_2& {}^{R1}x_1-{}^{R1}x_3& {}^{R1}x_1-{}^{R1}x_4\\ {}^{R1}y_1-{}^{R1}y_2& {}^{R1}y_1-{}^{R1}y_3& {}^{R1}y_1-{}^{R1}y_4\\ {}^{R1}z_1-{}^{R1}z_2& {}^{R1}z_1-{}^{R1}z_3& {}^{R1}z_1-{}^{R1}z_4\end{array}\right]\cdot {\left[\begin{array}{ccc}{}_{R2}x_1-{}_{R2}x_2& {}_{R2}x_1-{}_{R2}x_3& {}_{R2}x_1-{}_{R2}x_4\\ {}_{R2}y_1-{}_{R2}y_2& {}_{R2}y_1-{}_{R2}y_3& {}_{R2}y_1-{}_{R2}y_4\\ {}_{R2}z_1-{}_{R2}z_2& {}_{R2}z_1-{}_{R2}z_3& {}_{R2}z_1-{}_{R2}z_4\end{array}\right]}^{-1}$$

(5)

(3)

Compute the translation matrix: Substituting the rotation matrix ^R1R_R2 into:

$${}^{R1}P_{R2}={}^{R1}P-{}^{R1}R_{R2}\cdot {}^{R2}P$$

(6)

(4)

Orthonormalization of the rotation matrix: Due to computational errors, the derived rotation matrix ^R1R_R2 may violate orthonormality constraints. To reconcile the discrepancies between matrices while preserving geometric validity, we optimize the solution using the Frobenius norm criterion. The specific formulation is given as follows:

$${‖A‖}_F={\left({\displaystyle \sum _{i=1}^m{\displaystyle \sum _{j=1}^n{\left|a_{ij}\right|}^2}}\right)}^{1/2},\mathrm{A}\in R^{m\times n}$$

(7)

Define the objective function as:

$$J={‖R-{}^{R1}R_{R2}‖}_F^2$$

(8)

Constraints:$R\cdot R^T=I$

The optimized orthonormalized rotation matrix ${}^{R1}\stackrel{\wedge }{R}_{R2}$ is obtained when the objective function J reaches its minimum.

$$J\left({}^{R1}\stackrel{\wedge }{R}_{R2}\right)=\underset{R\cdot R^T=I}{{\displaystyle \min }}{‖R-{}^{R1}R_{R2}‖}_F^2$$

(9)

Base coordinate system calibration experiment:

The calibration experimental platform consists of an RC700-A controller with host computer, two robotic manipulators, and a calibration probe, as illustrated in Figure 2. Following the experimental procedure, a dual-robot base frame calibration was conducted, with the coordinate data of key points for both the master and slave robots recorded in Table 1 and Table 2 respectively.

The rotation matrix ^R1R_R2 is derived from Equation (5).

$${}^{R1}R{}_{R2}=\left[\begin{array}{ccc}1.000& 0& 0.010\\ -0.004& 1.000& 0\\ -0.011& -0.012& 1.006\end{array}\right]$$

(10)

The translation matrix ^R1P_R2 is derived from Equation (6).

$${}^{R1}P{}_{R2}=\left[\begin{array}{c}-1085.811\\ 8.298\\ 7.362\end{array}\right]$$

(11)

After orthonormalization, the transformation matrix between the base frames of the two robots is obtained as:

$${}^{R1}T{}_{R2}=\left[\begin{array}{cccc}0.9999& 0.0039& 0.00002& -1085.811\\ 0& 0.9999& 0.0039& 8.298\\ -0.0099& -0.00003& 0.9999& 7.362\\ 0& 0& 0& 1\end{array}\right]$$

(12)

2.2. TCP Calibration

This paper adopts the seven-point method for Tool Center Point (TCP) calibration.

The calibration procedure is as follows:

(1)

Determine calibration points: A calibration needle is fixed on the optical platform, with its tip position designated as P. Control the robot to position the tool end-effector at point P with four distinct orientations, as illustrated in Figure 3, while recording its pose relative to the base coordinate frame.

(2)

Formulate the system of equations: Applying coordinate transformation principles yields the relevant equations. Since the tool end-effector reaches the same position four times, simultaneous elimination of point P coordinates gives matrix:

$$AX=B$$

(13)

Where matrix A is composed of differences between rotation matrices of the flange end coordinate system under different orientations, matrix B consists of corresponding positional differences, and X represents the tool center point (TCP) position in the end-effector coordinate system (^E2P_T).

(3)

Solve using least squares method: Obtain ^E2P_T2 by applying the least squares principle:

$$X={\left(A^{T}A\right)}^{-1}{A}^{T}B$$

(14)

(4)

Initial Position Adjustment and Recording: Move the robot near the origin point, adjust the joints to align the tool coordinate system’s X and Z axes parallel to the corresponding axes of the base coordinate system, then record the data parameters (for the 5th calibration point).

(5)

Motion Position Recording: Maintaining the current orientation, move from Point 5 along the X and Y axes of the base coordinate system by specified distances respectively, then record the parameters for Points 6 and 7 as illustrated in Figure 4.

(6)

Direction Vector Determination:

$$PX={}_{T2}{}^{R2}P{}_6-{}_{T2}{}^{R2}P_5$$

(15)

$$PZ={}_{T2}{}^{R2}P{}_7-{}_{T2}{}^{R2}P_5$$

(16)

(7)

Transformation Matrix Computation: Orthonormalize PX, PY, PZ and substitute to obtain:

$${}^{E2}R{}_{T2}={\left({}^{R2}R_{E2}\right)}^{-1}\cdot {}^{R2}R_{T2}$$

(17)

Combining the obtained ^E2P_T2, we derive ^E2T_T2:

$${}^{E2}T{}_{T2}=\left[\begin{array}{cc}{}^{E2}R_{T2}& {}^{E2}P_{T2}\\ 0& 1\end{array}\right]$$

(18)

TCP Calibration Experiment:

By operating the Epson C4-A601S robot through its controller, the tool end-effector was repeatedly positioned at the same target point in 4 distinct spatial orientations (as illustrated in Figure 5). The end-effector pose data for each repetition were recorded. The experimental dataset has been compiled in Table 3.

The position matrix ^E2P_T2 is obtained according to Equation (14):

$${}^{E2}P{}_{T2}=\left[\begin{array}{c}0.6442\\ -129.4600\\ -147.4710\end{array}\right]$$

(19)

Furthermore, since the Y-axis, Z-axis, and X-axis of the tool coordinate system are parallel to the X-axis, Y-axis, and Z-axis of the base coordinate system respectively, PX and PY represent the directional vectors of the Y-axis and Z-axis in the tool coordinate system. The robot was controlled to move 50mm along the X-axis and 50mm along the Y-axis of the base coordinate system, with the data shown in Table 4.

The matrices ^R2R_T2 and ^E2T_T2 are computed using (17) and (18).

$${}^{R2}R{}_{T2}=\left[\begin{array}{ccc}0& 1& 0\\ 0& 0& 1\\ 1& 0& 0\end{array}\right]$$

(20)

$${}^{\mathrm{E}2}T{}_{T2}=\left[\begin{array}{cccc}1& 0& 0& 0.6442\\ 0& 1& 0& -129.4600\\ 0& 0& 1& -147.4710\\ 0& 0& 0& 1\end{array}\right]$$

(21)

3. Master Robot Modeling and Verification

3.1. Robot Coordinate Frame

The EPSON C4-A901S robot is a serial articulated robot with six revolute joints. The modified Denavit-Hartenberg (D-H) parameter method was employed to establish the robot link coordinate system, as illustrated in Figure 6. Furthermore, the MDH parameters of the EPSON C4-A901S robot were obtained and are presented in Table 5. The EPSON C4-A601S robot differs in the MDH table only by the parameters a₃=250 and d₄=-250.

β denotes the minor angular deflection between nominally parallel axes; β_i represents the joint offset of the i-th joint; θ_i is the angle between two adjacent links; ai indicates the link length; α_i specifies the angle between two normal; and d_ᵢ defines the distance between two normal.

3.2. Simulation and Verification of Robot MATLAB Model

Based on the MDH parameters presented in Table 5, the robotic model of EPSON C4-A901S was simulated using the Link function from MATLAB Robotics Toolbox. The simulation results of the EPSON C4-A901S robot are obtained as shown in Figure 7.

The homogeneous transformation matrix is obtained by solving the simulation model using the fkine command:

$$T=\left[\begin{array}{cccc}0& -1& 0& 0\\ 0& 0& -1& 565\\ 1& 0& 0& 720\\ 0& 0& 0& 1\end{array}\right]$$

(22)

The homogeneous transformation matrix is obtained using the RC+ software provided with the EPSON robot:

$$T^{\text{'}}=\left[\begin{array}{cccc}0& -1& 0& 0\\ 0& 0& -1& 565\\ 1& 0& 0& 720\\ 0& 0& 0& 1\end{array}\right]$$

(23)

The MATLAB model demonstrates good agreement with the actual robot system.

Within the constrained workspace, Figure 8(a) shows the robot motion from initial position A (0, 0, 0, 0, 0, 0) to B (pi/2, pi/6, -pi/6, pi/3, -pi/4, 0) through teach pendant programming, while Figure 8(b) presents the results obtained via EPSON RC+ software control.

By comparing the recorded position and orientation data of the robotic arm in Figure 8, it can be concluded that the MATLAB-based robot model parameters exhibit high consistency with the actual robotic system. As shown in Figure 9, the joint motion of the EPSON C4-A901S robot model demonstrates excellent smoothness over time without abrupt changes. This indicates that the established simulation model accurately characterizes the kinematic behavior of the physical system.

3.3. MATLAB-Based Modeling and Simulation of the Master Robot

For ^E1T_T1, as illustrated in Figure 1, the gripper is connected to the end flange along the Z-axis via bolts, from which the transformation matrix T is derived through measurement.

$${}^{E1}T{}_{T1}=\left[\begin{array}{cccc}1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& 1& 230.5\\ 0& 0& 0& 1\end{array}\right]$$

(24)

Similarly, the transformation matrix^T1T_w is obtained.

$${}^{T1}T{}_w=\left[\begin{array}{cccc}1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& 1& 10\\ 0& 0& 0& 1\end{array}\right]$$

(25)

The master robot model was simulated using the tool function from the MATLAB Robotics Toolbox, as shown in Figure 10. In robotic kinematics and coordinate system definitions, the tool coordinate frame (TCF) is established as a local frame relative to the end-effector coordinate frame. Through the robot’s teach mode, the tool coordinate frame’s pose relative to the base frame is displayed, resulting in a change of the y-coordinate value from 565.0 to 805.5.

4. Collaborative Motion Simulation Analysis

A dual-robot grinding simulation platform was established using MATLAB Robotics Toolbox, with subsequent motion simulation analysis of the collaborative grinding system.

4.1. Establishment of the Dual-Robot Model

Based on the derived ^R1P_R2 and ^E2T_T2 transformations, a dual-robot grinding simulation model was established using MATLAB Robotics Toolbox, as illustrated in Figure 11.

4.2. Cooperative Motion Simulation of Dual Robots

Within constrained space, the cup’s edge surface is chosen for grinding; the actual workpiece is illustrated in Figure 12. Dual-arm trajectory control: The master robot performs linear spatial motion with fixed end-effector orientation, while the slave robot tracks a circular edge path on the cup within its restricted workspace.

The slave robot is positioned at the designated starting point for grinding, and its pose data is acquired to derive the transformation matrix of the workpiece relative to the base coordinate system as:

$${}^{R1}T{}_w=\left[\begin{array}{cccc}0& -0.01& -1& -343.223\\ 0& 1& -0.01& -115.673\\ 1& 0& 0& 578.234\\ 0& 0& 0& 1\end{array}\right]$$

(26)

Based on the calibration results, Equations (1-1) and (1-2) are used to determine the end-effector trajectories and corresponding joint angles for both robotic arms. The simulation results obtained from these equations are shown in Figure 13 and Figure 14.

As shown in Figure 13a and Figure 14a, the joint angles of all six joints for both robots demonstrate smooth variation characteristics over time, with no observable abrupt joint motion phenomena. Furthermore, comparative analysis of Figure 13b and Figure 14b demonstrates satisfactory position tracking accuracy for both robotic manipulators. Thus, the feasibility and effectiveness of the proposed dual-robot grinding motion planning model are successfully validated.

5. Coordinated Motion Experiment

As illustrated in Figure 15, the dual-robot coordinated motion hardware platform comprises the following components: two 6-DOF industrial robots (EPSON C4-A601s and EPSON C4-A901s), an upper computer, an air compressor, two RC700-A robot controllers, a pressure sensor, gripping tools, and grinding tools.

The system’s software architecture is implemented through EPSON’s native RC+ development platform, employing SPEL+ programming language to establish Ethernet-based communication with robot controllers for real-time robotic control.

During the experimental procedure, the host computer was initially interfaced with both EPSON C4-A901S controllers via Ethernet and USB communication protocols respectively. Subsequently, the SPEL+ programming language was employed to compile motion programs for both the master and slave robots. The compiled programs were then uploaded to the respective robot controllers. Upon receiving the code, each controller executed an internal processing sequence comprising: (1) path planning generation, (2) motion trajectory calculation, (3) conversion to motor control signals, and (4) actuation of robotic joints. Finally, joint angle data and positional coordinates during operation were recorded through the execution of Call instructions.

The critical aspects of dual-robot collaborative motion are illustrated in Figure 16. During dual-robot collaborative motion, the simulated reference trajectory and actual tracking trajectory of the master-slave robot end-effectors demonstrate high consistency, as shown in Figure 17(a). This phenomenon preliminarily verifies that the dual-robot collaborative simulation model exhibits high accuracy and reliability in simulating real-world motion scenarios. To quantitatively evaluate the positional tracking accuracy of dual-robot collaborative operation, an in-depth analysis was conducted on the motion errors of the grinding tool along the X, Y, and Z axes (as shown in Figure 17(b)). Precision measurements and data processing revealed maximum motion errors of the grinding tool along the X, Y, and Z axes to be 2.158 mm, 1.81 mm, and 0.855 mm respectively, all within the 2.5 mm tolerance threshold. The observed errors primarily originate from calibration inaccuracies, deviations between theoretical and actual MDH parameters of the robots, and mechanical structural errors. These results demonstrate that the dual-robot collaborative simulation model can effectively guide practical operations, ensuring precise tracking of desired trajectories during coordinated motion while meeting stringent performance requirements for positional accuracy in collaborative tasks. This provides robust technical assurance for efficient and stable dual-robot collaborative operations.

6. Conclusions

This section is not mandatory but can be added to the manuscript if the discussion is unusually long or complex.

This study addresses the high-precision positioning challenge in dual-robot collaborative grinding and polishing by proposing an innovative, simple yet efficient calibration method that integrates the seven-point method with handshake calibration. Using a calibration needle as the reference tool, the method achieves rapid and precise positioning of dual robots, breaking through the accuracy bottleneck of system coordination in complex working scenarios.

This study establishes a system of kinematic constraint equations for dual robots based on their closed-loop motion chain configuration. Based on the process parameters and technical requirements of the grinding task, this study employs trajectory planning algorithms to generate precise grinding paths. Under the master-slave control architecture, this study first plans the ideal posture trajectory for the master robot, then determines the slave robot’s pose parameters through inverse kinematics solving based on real-time updated collaborative motion constraint equations, thereby achieving efficient coordination of dual robots in grinding operations. During system calibration, the “handshake” calibration method is employed to align the base frames of dual robots. By establishing transformation relationships between the global coordinate system and individual robot coordinate systems, base frame installation errors are effectively eliminated. A seven-point calibration method is employed to achieve precise positioning of the Tool Center Point (TCP), where coordinated measurements via the robot’s end-effector ensure high-accuracy acquisition of TCP coordinates. The integration of these two calibration methods fulfills the high-precision operational requirements of dual-robot systems. In robotic modeling, an improved Denavit-Hartenberg (D-H) parameter method was employed to establish the kinematic model of the gripper robot (EPSON C4-A901S), with experimental validation confirming the model’s effectiveness and accuracy.

Based on the aforementioned calibration results, a dual-robot collaborative grinding simulation model was developed using the MATLAB Robotics Toolbox, followed by trajectory tracking experiments. Experimental results demonstrate that the grinding model achieves position tracking accuracy within 2.5mm, with comparative analysis against physical measurements confirming the validity and reliability of the dual-robot kinematic model. This research significantly enhances the positioning accuracy of dual-robot systems and improves their adaptive recognition capability for end-effector tools, thereby establishing a robust yet straightforward technical solution for high-precision collaborative grinding of complex workpieces.