Adaptive Impedance Control for Force Tracking in Manipulators Based on Fractional-Order PID

: Force tracking control in robot arms has been widely used in many industrial applications, particularly in tasks involving end effectors and environmental contact, such as grinding, polishing, and other similar operations. However, these environments are not always precisely known. In order to address the force tracking control problem in unknown environments, this paper proposes a fractional-order PID adaptive impedance control strategy based on traditional impedance control. The unknown environmental information is estimated online using the adaptive impedance control algorithm, and the estimated parameters are used to generate reference trajectories to reduce force tracking errors. Fractional-order PID control is then introduced into the system to improve the control performance of the system model, and the theoretical proof of strategy stability is conducted. Finally, a comparison of four strategies was conducted through simulations: traditional impedance control, adaptive hybrid impedance control, adaptive variable impedance control, and the fractional-order PID impedance control proposed in this paper. The simulation results demonstrate that the strategy proposed in this paper exhibits robustness, virtually eliminates overshoot, and enhances response speed. In contrast, both adaptive hybrid impedance control and adaptive variable impedance control exhibit approximately 30% to 45% overshoot during interactions with the environment. Furthermore, in terms of force tracking error, the proposed strategy in this paper outperforms the above two strategies by approximately 29% to 60%, achieving excellent force tracking control performance.


Introduction
As we look to the future, it seems likely that robots will assume a range of new roles in people's daily lives [1].They may replace human labor in some industries, provide valuable services, or even offer specialized training to individuals [2].Some people, when compared to communicating with humans in daily life, even prefer to use robots as communication partners in many roles and situations [3].These developments promise to revolutionize the way we live and work, giving rise to new possibilities and transforming our daily routines in unexpected ways [4].In recent years, robots have become an increasingly common sight in a variety of human environments, including nursing homes, hospitals, shopping malls, and the education industry [5,6].In the field of robotics, there has long been a focus on designing and controlling robots that can interact with their environment, particularly in cases where the objects of interaction are human [7].
As robots become increasingly integrated into open working environments, it is inevitable that robots will come into contact and interact with humans [8], and ensuring safe interactions between these machines and humans is becoming a key area of focus in robotics research.When humans provide motion correction, fine-tuning control, and situational guidance to working robots, the proximity required for these interactions can pose a significant risk of injury.This challenge has driven the development of new approaches to robot control, with an emphasis on creating security rules that enable safe and controlled physical interactions between robots, humans, and the environment [9,10].Therefore, to guarantee the safety of the surrounding environment and humans around robots, it is necessary to study the interaction mechanisms that enable robots to adapt to the surrounding environment [11,12].In addition, there is a need to balance the compliance and accuracy of the robot.
Impedance control was first proposed by Hogan [13], and later Kazeroni developed and implemented an automated impedance control method for the precise deburring of industrial robots to replace the time-consuming and expensive manual deburring approach [14].However, impedance control still faces many challenges: the nonlinear and time-varying characteristics of the robot's own dynamics and kinematics make it insufficient to rely solely on impedance control to solve practical problems [15].Therefore, in practical applications, impedance control is often combined with other control methods, such as adaptive control.
In Oladayo's study [16], an adaptive force tracking impedance controller was developed, and the controller's parameters were adjusted using a genetic algorithm (GA) under force-torque constraints.The controller adjusts the contact force by adapting the target stiffness in response to changes in the stiffness of the contact environment, with the aim of achieving adaptive control.In a study by Wei et al. [17], a force-based impedance control system was designed, and three self-calibrating functions for impedance parameters were constructed to improve the design of the impedance control system to achieve adaptive control.Chen et al. [18] optimized the PD control parameters by using particle swarm optimization and adaptive iterative algorithms based on the robot's machining state.K. Wakita et al. [19] optimized a robot's position control and reference trajectory by estimating human movement direction and using a Kalman filter, and they also achieved adaptive control for robots.
However, some adaptive impedance control methods, while showing promising performance compared to traditional impedance control, still have many issues.Zhou H et al. [20] proposed an adaptive control-based hybrid impedance control strategy for achieving compliant force control in industrial robots.Thunyajarern et al. [21] presented an adaptive force controller for six degrees of freedom robots and addressed the environmental identification problem by using particle swarm optimization.However, the optimization algorithm was performed offline and did not yield a satisfactory adaptation to the environment.Zhou Y et al. [22] proposed a variable impedance control algorithm that achieved superior force tracking control when compared to traditional impedance control.However, the end effector exhibited significant overshoot when in contact with the environment, potentially leading to equipment damage.Cao et al. [23] combined PID controllers with impedance control to enhance control performance during contact.However, traditional PID controllers have limited degrees of freedom and are prone to integral saturation.Li et al. [24] incorporated PID control into the design of a model reference adaptive impedance controller to improve impedance relationships.Nevertheless, significant overshoot still exists during the contact process, and a stability analysis of a system after introducing PID has not been provided.Zhao et al. [25] proposed an improved fractionalorder PID impedance control method based on traditional impedance control, achieving good force tracking control results.However, their proposed method lacks theoretical derivation of its validity and stability.
In order to address the issues and limitations observed in the aforementioned control algorithms, such as overshooting during robotic contact with the environment, offline algorithm optimization, and the lack of corresponding stability theory derivation, this paper proposes a fractional-order PID adaptive impedance control strategy and provides theoretical stability and validity derivations.This strategy allows for the real-time dynamic estimation of environmental information in uncertain and changing environments.It updates the reference trajectory in real-time based on estimated parameters to achieve adaptive control.Additionally, the flexibility of fractional-order PID further enhances force tracking control performance in complex environments, minimizing overshoot and maintaining stable force tracking.This ensures equipment safety and smooth task execution in applications with stringent machining precision requirements within complex environments.
The main contributions of this paper are as follows: a fractional-order PID adaptive impedance control strategy was proposed to enhance the performance of robotic interactions with the environment and the subsequent force tracking; a mathematical derivation of the rationality and stability of the proposed strategy; and the validation of its superiority through comparative simulations.
The remaining parts of the paper are organized as follows.Section 2 describes the impedance control system studied in this paper.Section 3 analyzes the stability of the impedance model and presents the design of the adaptive control with an online estimation of environmental parameters.Section 4 proposes the design of the fractional-order PID impedance controller and provides the theoretical proof of stability.Section 5 sets up simulations to validate the proposed strategy.Finally, Section 6 provides the conclusion of the paper.

System Description 2.1. Force Interaction Model
Generally, the process of a robot coming into contact with an unknown environment can be described as transitioning from free-space motion to contact and then to gradually stabilizing the contact forces.The contact process is shown in Figure 1.During the contact process, the dynamic characteristics of the robot and the behavior of the contact force vary at different stages.Figure 2 shows the variation of contact force during the contact process.
As shown in Figure 2, the contact process of the robot can be divided into three stages: Stage I: linear motion in free space when approaching the environment.
Stage II: a nonlinear region with significant overshoot and intense collision when contact occurs.
Stage III: an approximately linear region where the contact force gradually stabilizes after contact.
This paper focuses on Stage II and Stage III after contact, as Stage I is a basic ability that the robot should possess.

Impedance Control Model
The core idea of impedance control is to transform almost all physical systems into a "mass-spring-damper" system, as shown in Figure 3, and the force situation of this system can be expressed as In the equation, E = X d − X m , where X d is the desired position, and X m is the actual position.M r is the desired inertial matrix, B r is the desired damping matrix, and K r is the desired stiffness matrix.This paper employs position-based impedance control (sometimes referred to as admittance control) [26] to manipulate industrial robots.This control structure consists of an internal position control loop and an external force control loop, which equates the controller to an admittance system: input force; output position and equates the robot to an impedance system: input position; output force.Overall, it maps from force to position.The schematic diagram is shown in Figure 4.Where X r is the reference trajectory corresponding to the expected force, X d denotes the trajectory information sent by the impedance model to the robot, and X represents the trajectory information of the end effector of the robot's actual movement.In the position control mode, X d is usually equal to X, and the trajectory error caused by contact force can be expressed as E = X − X r = X d − X r .F r is the desired tracking force, F e is the actual contact force, and the error between the two is ∆F = F e − F r .
Usually, the error force ∆F is represented as a linear second-order system, as shown in Equation ( 1), which can be obtained by substituting F and E into it.
where M r is the desired inertia matrix, B r is the desired damping matrix, and K r is the desired stiffness matrix.The impedance model is a second-order transfer function: where the force error ∆F = F e − F r and the position error E = X r − X d , from which the trajectory information sent to the robot can be obtained as X d = X r + E.

Stability Analysis of Impedance Control Models
Usually, the model parameters, M r , B r , and K r , of the impedance controller and the environmental dynamics parameters, K e , are positive definite diagonal matrices, with the forces in each direction decoupled.Therefore, for the convenience of subsequent research analysis in this paper, the applied force is decomposed and studied.It is assumed that the contact force only acts in one direction, which is the direction of the force.Lowercase letters x, x r , x d , m r , b r , k r , f r , f e , x e , k e , and b e are used to represent variables and coefficients in a single direction, respectively, instead of the uppercase letters X, X r , X d , M r , B r , K r , F r , F e , X e , K e , and B e .Thus, the impedance relationship in Equation ( 2) can be expressed in a single direction: We assume the environmental dynamics model to be as is shown in Equation ( 5), where x e , k e , and b e represent the environmental position, environmental stiffness, and environmental damping, respectively.
Additionally, when the robot makes contact with the environment, its acceleration and movement speed are very small, and the contact force between the robot and the environment can be simplified to Taking the Laplace transform of Equation ( 4) yields By transforming and further simplifying Equation (7), we obtain where As mentioned in Section 2.2, it is generally assumed that there is no tracking error in the position control mode of the robot, that is, x = x d .By incorporating the desired force into Equation ( 9), we obtain Substituting Equation ( 8) into Equation (10) yields The steady-state tracking error in impedance control can be expressed as follows: From Equation ( 13), we can see that setting ( f r − k e (x r − x e )) = 0 or k r = 0 can make the steady-state tracking error ∆ f ss converge to zero as the system approaches a stable state.However, setting the desired stiffness k r to zero directly may reduce the force steady-state error to zero, but this would result in the system losing its rigid response to external forces.This would make the system unable to effectively resist external disturbances and may cause the system to become excessively compliant, slowing down its response speed and thereby affecting the system's dynamic performance.Therefore, the reference trajectory x r is given by From Equation ( 14), we can see that the reference trajectory x r is a dynamic function composed of the environmental position x e , environmental stiffness k e , and desired force f r .However, this formula requires accurate environmental information, which is often inaccurate or unknown in the actual operation of the robot.As a result, traditional impedance control struggles to achieve precise force control during contact.

Adaptive Control Based on Environmental Information Estimation
From the previous subsection, it is known that the actual working environment is typically time-varying and unknown [27].In such cases, the real-time online estimation of environmental parameters is required.Let ke and xe represent the estimated values of k e and x e , respectively.Then, Equation ( 14) can be expressed as When substituting the estimated values of the environmental parameters into Equation ( 6), and by taking the difference with Equation ( 6), we can obtain the estimated value of the contact force fe and the estimation error of the contact force: When simplifying Equation (17), we obtain where, φ k = ke − k e , φ x = ke xe − k e x e , φ = [φ k φ x ] T .Therefore, it is possible to set the control law to an online estimate of the environmental parameters ke and xe such that as t approaches infinity, fe approaches f e , and f e approaches f r , achieving the desired force tracking objective.The design of the Lyapunov function is shown below: where, Γ = [r 1 , r 2 ]I 2×2 is a positive definite matrix, r 1 and r 2 are positive real constants, and φ is represented as By combining Equations ( 19) and (20), and taking the derivative of Equation ( 19), we obtain Because V is positive definite, it is evident from Equation ( 21) that it possesses the semi-negative definite property, ensuring the stability of the controller.When substituting the positive definite symmetric matrix Γ into Equation (20), we obtain According to the analysis, simplifying Equation (22) yields For convenience, we define η 1 = 1/r 1 and η 2 = 1/r 2 .Then, the adaptive control algorithm based on environmental information can be represented as follows: From Equations ( 24)-( 26), it can be observed that the designed adaptive control algorithm estimates the environmental information parameters x e and k e in real time based on the robot's end-effector position and actual contact force; it then computes the desired trajectory for the robot.Utilizing this control approach ensures that fe converges to f e , achieving the desired force tracking objective.

Control Method Design
In traditional impedance controllers, the end effector usually exhibits significant overshoot upon contact with the environment, and it has a slow response to the error in contact force [28].In order to improve the performance of the manipulator when interacting with the environment, this paper combines the advantages of fractional-order PID and proposes a novel fractional-order PID adaptive impedance control method.The fractionalorder integration and differentiation help prevent integral saturation, reduce contact force overshoot, enhance system robustness, and further improve control performance.
Recently, fractional calculus has been applied in many fields, such as artificial neural networks [29], computer vision [30], and the extension of PID controllers [31].Due to the differentiation and integration orders no longer being limited to integers but extended to any real number, it can provide additional flexibility to PID controllers, allowing them to be extended to fractional-order PID (FOPID) controllers [32].This enables it to better accommodate complex systems with time-varying and nonlinear characteristics.The equation for fractional calculus is as follows: where D is the fundamental operator of calculus, a and t are the boundaries of this operator, and q denotes that this operation is not restricted to integer orders.Currently, various mathematical representations of fractional calculus can be found in the literature, with the most common ones being Caputo, Riemann-Liouville (RL), and Grunwald-Letnikov (GL) [33][34][35].
Equation (28) provides the definition of RL: Equation ( 29) provides the definition of Caputo: where (n − 1 < q < n), and Γ(q) is represented by the following equation.
Similarly, Equation (31) provides the definition of GL: Unlike integer-order functions, fractional-order functions cannot obtain their exact solutions, and approximate methods are usually used for calculations.Among them, the pole distribution-based approach proposed by Oustaloup has shown good approximation effectiveness in achieving the desired accuracy [36].Therefore, this provides an implementation method for fractional-order PID.The approximate fractional filter within the frequency range (w b , w h ) can be represented as where α represents the order of the differentiator, w k , w k , and K, respectively, denote its zeros, poles, and gain.
where ω u = √ ω h /ω b and the control effort of the PID controller is expressed as Its transfer function C(s) can be represented as When introducing fractional order into PID control, the transfer function of fractional order PID can be represented as where, µ and λ represent the differentiation and integration orders of the controller, respectively.The expression for introducing fractional-order PID impedance control is as follows: where ∆x = x d − x r , ∆ f = f r − f e , k p , k i , and k d are the proportional, integral, and derivative parameters of the fractional-order PID controller, respectively.By selecting appropriate fractional-order PID parameters, Equation (37) will achieve a faster response, more stable force tracking control, and lower force impact when compared to Equation (4).Converting Equation ( 37) into the frequency domain yields When combining Equation ( 12) with the final value theorem, we can obtain the expression for the steady-state error of the system: From Equation ( 39), we can clearly see that after introducing the fractional-order PID controller, the expression in the numerator representing the reference trajectory is not affected; only the denominator is modified, and the addition to the denominator can reduce the steady-state error.Therefore, we can conclude that the introduction of the fractionalorder PID is reasonable, as it not only reduces the steady-state error and improves force tracking control performance but also coexists with the adaptive control strategy that estimates the environment parameters online.
Based on the above, we can obtain the principle block diagram of the adaptive impedance control strategy based on fractional-order PID, as shown in Figure 5.

Stability Analysis
The simplified model for the interaction between the manipulator and the environment is shown in Figure 6.In the diagram, the fractional-order impedance control model is represented as υ(s) = k p + k i s λ +k d s µ m r s 2 +b r s+k r , and among them, 0 < λ, µ ≤ 1.It is typically assumed that there is no error in the robot's position control, which means R(s) = 1.The environmental interaction model is represented as K(s) = k e .According to the superposition theorem, the closed-loop output of the control system (s) is given by Equation (40), as shown below: Lemma 1 ([37]).A commensurate order system described by following rational transfer function G(s) = Q(ω)/P(ω) where ω = s α , α ∈ R + and (0 < α < 2) are stable if |arg(ω i )| > απ/2 with ω i is the i-th root of P (ω) = 0.
When λ = µ = 1, the integer-order PID control system output is as shown in Equation (41): Obtain the characteristic equation of the integer-order system from Equation (41) to derive its stability conditions: According to the Routh-Hurwitz criterion (here, for ease of derivation, we set m r = 1), we can obtain the stability conditions for Equation (42) as follows: According to Lemma 1, by comparing the stability conditions and boundaries of the two control systems, we can find that when 0 < α < 1,|arg(ω i )| > π/2 > απ/2.It implies that the stability region of the fractional-order system is larger than that of the integer-order system.When the integer-order system is stable, its stability conditions (Equation (43)) will definitely apply to the fractional-order system [38], ensuring the stability of the proposed strategy.

Model of Manipulator
Robot joints typically exhibit highly nonlinear and strongly coupled dynamics [39].In order to mitigate the adverse effects of these conditions on the dynamic control of robots, this paper validates the proposed adaptive control strategy on a two-link manipulator model, as shown in Figure 7.In this model, m 1 and m 2 represent the weights of link 1 and link 2, while l 1 and l 2 represent their respective lengths.The paper assumes that only the x-axis direction is affected by forces.
The dynamic model of the robotic arm is represented by Equation (44), M 0 (q) q + C 0 (q, q) q + G 0 (q where M 0 ∈ R 2×2 is the inertia matrix of the arm itself, q, q, and q ∈ R 2 represent the joint angles, velocities, and accelerations of the robotic arm, respectively.C 0 q, q ∈ R 2 represents the Coriolis and centrifugal forces, G 0 ∈ R 2 is the gravity term, and τ, τ d ∈ R 2 denotes the control torque and external disturbances, respectively.

Algorithm Simulation
In order to validate the effectiveness of the proposed strategy, we conducted comparative tests on a two-degrees-of-freedom robotic arm model in two aspects (robustness test and force tracking test) against four control algorithms: traditional impedance control (CIC), adaptive hybrid impedance control (AIC) [20], adaptive variable impedance control (VIC) [22], and the fractional-order PID adaptive impedance (FOAIC) control proposed in this paper.The common parameters of these four algorithms remained consistent across various simulation environments.

Control Robustness Simulation
For the four control algorithms, we conducted robustness simulations in three types of perturbed environments: sudden changes in environmental stiffness, environmental position, and desired force.In the selection of abrupt parameter values, we take into consideration factors such as mathematical theory and the limitations of the simulated experimental manipulator and their impact.The parameters are randomly chosen within suitable ranges.The specific perturbation parameters and timing are presented in Table 1.During each robustness test involving a particular perturbation, the remaining parameters were kept constant.The desired force for stiffness and position perturbations was set at f r = 100 N. The simulation results of the robustness tests are illustrated in Figure 8.As shown in Figure 8a, when there is a sudden change in environmental stiffness, all four algorithms demonstrate good force tracking control.However, FOAIC performs the best in terms of overshoot and settling speed during contact with the environment and when the environmental stiffness changes.Similarly, as depicted in Figure 8b,c, during the position and desired force perturbations, CIC exhibits the slowest settling speed.AIC and VIC both exhibit overshoot during contact with the environment and when conditions change.On the other hand, FOAIC maintains the fastest response speed while significantly reducing overshoot.The simulations conducted under these three perturbation scenarios collectively illustrate the notable robustness of FOAIC in force tracking.

Force Tracking Control Simulation
In reality, when manipulators engage in tasks involving interactions with the environment, the contact with the environment is almost never flat, and it is likely to be accompanied by time-varying desired forces.Similarly, for the four control algorithms, simulations were conducted in four different environments (plane, inclined plane, sinusoidal surface, and complex surface) under three conditions for each (constant force, variable force tracking, and position tracking).Additionally, when it comes to setting simulated parameter values for different environmental conditions, we considered various combinations and factors.These factors include the theoretical considerations in Section 4, prior literature settings, and the constraints imposed by the manipulator used in the simulation experiments.Furthermore, for different environments and conditions, we incrementally introduced time-varying and increased complexity into the simulated parameter values.This approach enhances the persuasiveness of the simulation results and provides deeper insights for future research in this area.
Assuming that the robot end effector is in contact with a plane, the environment stiffness is set to k e = 5000 N/m.The environment position is set to x e = 0.3 m, and the robot impedance parameters are configured as m r = 0.15, b r = 300, and k r = 80.The estimated initial values of the environment parameters are x e (0) = 0.295 m, k e (0) = 4500 N/m, η 1 = 75, and η 2 = 13.The parameters for the FOAIC control are k p = 0.85, k i = 0.1, k d = 0.3, λ = 0.5, and µ = 0.8.The simulation time is set to 5 s with a step size of 0.0002.Additionally, the basic parameters for the subsequent force tracking simulations are kept consistent with the above-mentioned parameters.For constant force tracking, the desired force is set to f r = 12 N, while for variable force tracking, the desired force is set to f r = 12 + 3sin(π/3t).The simulation results are illustrated in Figure 9. From Figure 9a, it can be observed that during the flat force tracking with ẋe = ẍe = 0, when the parameter k r of CIC is set to 0, CIC achieves excellent force tracking, confirming the derivation of steady-state tracking error in Section 3.1.According to the simulation results, on the plane environment, all four algorithms exhibit good force tracking performance with minimal overshoot, yet FOAIC demonstrates the fastest response speed.
Assuming that the robot end effector is in contact with an inclined plane, the environment position is set as x e = 0.3 + 0.05t, ẋe = 0, ẍe = 0.The FOAIC control parameters are set as k p = 0.95, k i = 0.15, k d = 0.3, λ = 0.5, µ = 0.8.The desired force fr for constant force tracking and variable force tracking is set as 12 N and 12 + 3sin(π/3t)N, respectively.The simulation results are presented in Figure 10.As shown in Figure 10a, during the force tracking control on the inclined plane, CIC exhibits a force tracking error of about 2N.On the other hand, AIC, VIC, and FOAIC maintain good force tracking performance.However, AIC and VIC have overshoots of 34.5% and 25.6%, respectively, while FOAIC demonstrates minimal overshoot.Larger overshoots can result in impacts when the robot interacts with the environment, increasing the difficulty of actual robot force control and potentially causing damage to the robot, leading to unnecessary losses.FOAIC significantly reduces overshoot while achieving a shorter time to reach a steady state, showcasing excellent force tracking performance.
Assuming that the robot end effector is in contact with a sinusoidal surface, where x e = 0.3 + 0.01sin(2πt/7), ẋe = 0, and ẍe = 0.The FOAIC control parameters are set as k p = 0.9, k i = 0.23, k d = 0.25, λ = 0.8, and µ = 0.8.The desired force f r for constant force tracking and variable force tracking are set as 12 N and 12 + 3sin(πt/3)N, respectively.The simulation results for force tracking on the sinusoidal surface are presented in Figure 11.As shown in Figure 11a, during the force tracking control on the sinusoidal surface, CIC exhibits poor force tracking performance, being unable to track both the desired force and the environmental position.On the other hand, the other three algorithms still manage to maintain good force tracking control, albeit with minor steady-state errors.AIC and VIC, however, exhibit force overshoots of 38.3% and 33.4%, respectively, while FOAIC demonstrates almost no overshoot.Furthermore, for a comprehensive evaluation of the force tracking control performance of AIC, VIC, and FOAIC, quantitative and visual comparisons were conducted using statistical indicators, including mean absolute error (MAE), root mean square error (RMSE), and mean square error (MSE), as indicated in Equation (45).These indicators were employed for analyzing the reliability of dynamic force tracking (1 s-5 s).The comparison results are presented in Table 2. Based on the specific quantified data in Table 2, it is evident that in terms of MAE, MSE, and RMSE metrics, the force tracking control performance of FOAIC surpasses that of AIC and VIC.Taking MSE as an example, during the sinusoidal surface force tracking control, FOAIC achieves a 26% improvement in force tracking accuracy compared to AIC and a 45% improvement compared to VIC.This demonstrates its exceptional force tracking control performance.
Lastly, assuming that the robot end effector is in contact with a complex surface.The environmental position is set as x e = 0.3 + (0.015t)sin(πt/2)e (−1/2t) , ẋe = 0, and ẍe = 0.The FOAIC control parameters are configured as k p = 0.96, k i = 0.3, k d = 0.3, λ = 0.6, µ = 0.8.The desired values for constant force tracking and variable force tracking are set at 12 N and 12 + 3sin(πt/3)N, respectively.The simulation results are depicted in Figure 12.As shown in Figure 12a, CIC still fails to achieve effective force and position tracking control, exhibiting no force steady-state error.Additionally, the overshoot of AIC and VIC rises to 43.6% and 39.2%, respectively.However, FOAIC remains unaffected, maintaining negligible overshoot and a shorter settling time.Yet, due to the rapid changes in the environment, compared to the sinusoidal surface, AIC, VIC, and FOAIC display more noticeable force steady-state tracking errors.Similarly, based on Equation (45), a comprehensive quantitative and visual comparison of force tracking errors for these three algorithms is conducted, and the results are presented in Table 3.According to the data in Table 3, similarly, taking MSE as a quantifiable perspective, during the force tracking control on the complex surface, FOAIC achieves an improvement of 40.3% in force tracking accuracy compared to AIC, and an improvement of 58.6% compared to VIC.This indicates that in more complex environments, FOAIC demonstrates greater in force tracking accuracy, highlighting its advantage.Therefore, even in uncertain and complex environments, it maintains excellent force tracking control performance.
In summary, based on the simulation results of the four scenarios mentioned above, it can be observed that all four algorithms maintain satisfactory force tracking performance in a flat environment.Although CIC exhibits force steady-state error in the inclined plane environment, it still manages to maintain force tracking.AIC and VIC are able to perform force tracking stably but show overshoot, with longer adjustment times compared to FOAIC.In the sinusoidal and complex surface environments, CIC struggles to achieve stable force tracking, and both AIC and VIC experience increased overshoot.FOAIC, on the other hand, nearly eliminates overshoot, exhibits faster adjustment times, and, according to MSE quantification data, demonstrates superior force tracking accuracy compared to AIC and VIC, with even greater improvement in the complex surface environment.Hence, the proposed FOAIC is better suited for robot interactions in complex and dynamically changing environments, offering shorter adjustment times, minimal overshoot, and stable and reliable force tracking performance.

Conclusions
With the increasing number of robot tasks involving interaction with unknown environments, ensuring stability, speed, and steady force tracking during the interaction process has become increasingly important.This paper proposes a fractional-order PID adaptive impedance control strategy to enhance the stability and force tracking performance of robot-environment interactions and provides the theoretical stability and validity derivations.The algorithms CIC, AIC, VIC, and FOAIC were compared and validated through two aspects of simulation-robustness testing and force tracking testing.The results indicate that when unexpected changes in the parameters occur, such as desired force, FOAIC maintains excellent force tracking performance while nearly eliminating overshoot and exhibiting fast adjustment times, showcasing robustness.Regarding force tracking control in different environments, CIC is only suitable for plane environments.In sinusoidal and complex surface environments, AIC demonstrates an overshoot of 38.3% and 43.6%, respectively, while VIC shows an overshoot of 33.4% and 39.2%.FOAIC, on the other hand, has almost no overshoot.According to the MSE analysis, the force tracking accuracy of FOAIC surpasses that of AIC and VIC by 29% and 45%, respectively, in the sinusoidal environment, and by 40.3% and 58.6%, respectively, in the complex surface environment, showcasing outstanding force tracking performance.
Moreover, the proposed strategy demonstrates strong scalability, as it can adapt to various interactive tasks by estimating environmental information to update the robot's reference trajectory online.Therefore, this control strategy introduced in this paper can be applied to a range of force control tasks involving interactions with the environment, such as deburring, polishing, assembly, machining tasks in deformable environments, and production tasks based on robot force control.

Figure 1 .
Figure 1.The contact process between the robot and the environment.

Figure 2 .
Figure 2. The process of variation for contact force.
the transfer function, as shown in Equation(3).It represents the desired impedance model of the robotic arm, which converts the force error during the contact process into a correction for the reference trajectory.Taking the Laplace transform of the environmental dynamics model shown in Equation (6) yields f e (s) = k e x(s) − x e (s) , x > x e

Figure 6 .
Figure 6.Model of robotic interaction with the environment.

Table 1 .Figure 8 .
Figure 8. Robustness testing simulations of the algorithm.(a) Shift in environmental stiffness; (b) shift in desired force; (c) shift in environmental position.

Figure 9 .
Figure 9. Simulation of force tracking control on the plane.(a) Comparison of constant force tracking control; (b) comparison of variable force tracking control; (c) comparison of constant force position tracking control.

Figure 10 .
Figure 10.Simulation of force tracking control on the inclined plane.(a) Comparison of constant force tracking control; (b) comparison of variable force tracking control; (c) comparison of constant force position tracking control.

Figure 11 .
Figure 11.Simulation of force tracking control on the sinusoidal surface.(a) Comparison of constant force tracking control; (b) comparison of variable force tracking control; (c) comparison of constant force position tracking control.

Figure 12 .
Figure 12.Simulation of force tracking control on the complex surface.(a) Comparison of constant force tracking control; (b) comparison of variable force tracking control; (c) comparison of constant force position tracking control.