Next Article in Journal
Research on Error Correction Technology in Underwater SINS/DVL Integrated Positioning and Navigation
Next Article in Special Issue
Robotic Intracellular Pressure Measurement Using Micropipette Electrode
Previous Article in Journal
Dynamic Resource Allocation for Network Slicing with Multi-Tenants in 5G Two-Tier Networks
Previous Article in Special Issue
Prehensile and Non-Prehensile Robotic Pick-and-Place of Objects in Clutter Using Deep Reinforcement Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Surface Tracking and Constant Force Control of a Grinding Robot

1
School of Mechanical Engineering, Yanshan University, Qinhuangdao 066000, China
2
School of Vehicle and Energy, Yanshan University, Qinhuangdao 066000, China
3
Suzhou Huichuan Control Technology Co., Ltd., Suzhou 215000, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(10), 4702; https://doi.org/10.3390/s23104702
Submission received: 17 April 2023 / Revised: 5 May 2023 / Accepted: 8 May 2023 / Published: 12 May 2023
(This article belongs to the Special Issue Sensors for Robots II)

Abstract

:
To improve the quality and efficiency of robot grinding, a design and a control algorithm for a robot used for grinding the surfaces of large, curved workpieces with unknown parameters, such as wind turbine blades, are proposed herein. Firstly, the structure and motion mode of the grinding robot are determined. Secondly, in order to solve the problem of complexity and poor adaptability of the algorithm in the grinding process, a force/position hybrid control strategy based on fuzzy PID is proposed which greatly improves the response speed and reduces the error of the static control strategy. Compared with normal PID, fuzzy PID has the advantages of variable parameters and strong adaptability; the hydraulic cylinder used to adjust the angle of the manipulator can control the speed offset within 0.27 rad/s, and the grinding process can be carried out directly without obtaining the specific model of the surface to be machined. Finally, the experiments are carried out, the grinding force and feed speed are maintained within the allowable error range of the expected value, and the results verify the feasibility and effectiveness of the position tracking and constant force control strategy in this paper. The surface roughness of the blade is maintained within Ra = 2~3 μm after grinding, which proves that the grinding quality meets the requirements of the best surface roughness required for the subsequent process.

1. Introduction

The development of the wind power industry is of great significance to environmental protection, and there has been rapid development, with increasing requirements for wind turbines in recent years [1]. Therefore, the manufacturing process of wind turbine blades and other components is gradually becoming more intelligent and automated [2], and the process of grinding wind turbine blades is essential. The roughness after grinding affects the blade’s subsequent effectiveness, and the grinding process itself determines the production cost and service life of the blade. Due to the irregular shape and large size of a wind turbine blade’s surface, artificial grinding, which is currently the most common method of treatment, is not merely inefficient, but also struggles to achieve uniformity, and polished glass fiber dust can cause environmental pollution and damage to workers’ health [3,4]. In consequence, the design of, and research into, wind turbine blade grinding robots has gradually attracted more interest in recent years. Some of them employ a gantry structure [5,6] equipped with a grinding head on each side, whereby the gantry moves along a guide rail to perform the grinding of the blades. This model achieves high grinding efficiency but requires much space. Others use mobile platforms equipped with multi-degree-of-freedom manipulators and grinding devices [7], which can effectively adapt to the site, but with more difficult manipulator control.
With the gradual increase in industrial automation in recent years, the requirements relating to the control and repetitive positioning accuracy of robots have increased, and the environmental aspects of robots need to be improved. The surfaces of wind turbine blades are complex, with curvature variation in the axial and circumferential directions; thus, the control of the grinding robot is particularly important. To ensure a suitable grinding effect, it is essential to maintain a constant grinding force and the real-time tracking of the blade surface.
Most previous studies have focused on force control, which can be divided into passive control and active control [8]. Guo et al. [9] proposed a new two-dimensional adjustable force mechanism that does not require additional sensors. However, because of the low accuracy and poor stability of this method of passive constant force control, active force control has become a research hotspot.
Active force control includes impedance control and force–position hybrid control. Shen et al. [10] proposed a fuzzy-based impedance control algorithm, with which the contact force error could be controlled within 2N for the external environment with unknown stiffness. Wahballa et al. [11] used an OSRDF impedance controller to offset the damping force and the contact force tracking error, which improved the system stability.
PID control is a common method applied to achieve active force control. The process is proportional, integral and differential, enabling quick responsiveness and high stability in control, and the robustness of the mode of PID control is strong. To resolve the issues of model construction and algorithm complexity that arise in the active flexibility control method, Yao et al. [12] proposed a fuzzy PID smooth grinding control algorithm based on a two-degree-of-freedom constant force grinding device driven by pneumatic force, wherein the position and force are decoupled and controlled separately, and verified that the control precision and response time of the fuzzy PID control algorithm were better than those of the conventional PID control method. In addition, Sun et al. [13] confirmed the effectiveness of the fuzzy PID control algorithm. Dai et al. [14] also presented a backstepping + PID control method to closely track the applied force.
As such, many scholars use PID as the control strategy in their research. Xu et al. [15] proposed the combination of an active strategy, based on PI/PD control, and a passive strategy, based on PID control, which was shown to improve the accuracy and efficiency of the controlled force and avoid over- and under-cutting. Zhang et al. [16] presented a method for optimizing the parameters of a robot PD constant grinding force controller using deep reinforcement learning DRL Rainbow, which can ensure that a constant force is applied during grinding by adjusting the grinding depth. Further, some studies have been carried out on the normal vector of the surface; in establishing the positioning of the grinding tool during processing, Zhao et al. [17] proposed an adaptive PD constant force controller and a normal vector search algorithm, which together ensure that constant force is applied during the grinding of workpieces with unknown shapes. Wang et al. [18] put forward a force–position hybrid control method based on a PD constant force control algorithm. Han et al. [19] presented a fuzzy gain scheduling PID controller, which could reduce the influence of dynamic characteristics on the control system.
As regards the influence of curvature variance, beginning with an improved material removal model, Li et al. [20] took the material to be removed from different positions on the object being machined as the research object, and constructed a hybrid force–position control algorithm to ensure the contour accuracy of workpieces with complex curved surfaces.
Some scholars have proposed other new algorithms to ensure constant force control. In order to resolve the shortcomings of traditional robots’ constant force output, Jia et al. [21] combined the PID control method with a particle swarm optimization algorithm (PSO) to realize the constant force output of a cylinder and drive the grinding device. Zhao et al. [22] developed a strategy for constant force grinding that reconstructs the surface of the workpiece based on the point cloud data obtained via visual measurements, without requiring an accurate CAD model of the workpiece. In the grinding process, noise in the machining environment will affect the constant force control of the robot. Dai et al. [23] proposed a new force control method using an extended state observer (ESO) that can reduce fluctuations in the grinding force in the normal direction. Due to the impact between the grinding tool and the surface being machined, there is constant fluctuation in the contact force, and so Xiao et al. [24] put forward an active disturbance rejection method. Here, in order to achieve constant force control during machining, an RBF neural network model and an iterative algorithm are combined. Zhang [25] presented a variable-gain iterative learning force–position hybrid control method with ideal contour force as the contact point processing idea, which could still ensure the control accuracy of contour force when the tracking error was not accurate enough. Li et al. [26] used the force–position hybrid control method to estimate and adjust the contact state between the robot and the unknown workpiece, which could provide a constant normal contact force. Zhang et al. [27] proposed a constant force grinding controller based on proximal strategy optimization, which improved the response ability of the system.
The roughness of the outer surface of the wind turbine blade before grinding is about Ra = 0.8~1.5 μm, and the requirements are met when this value after grinding is Ra = 1.5~8 μm. A roughness of Ra = 2.0~3.0 μm yields the best results in the subsequent painting of the blade.
In this paper, a wind-turbine-blade-grinding robot, along with its control strategy and related algorithms, are introduced. Given the size and complexity of the surfaces of wind turbine blades, this paper considers the surface tracking and constant force control of blades simultaneously, so as to enhance the surface quality of wind turbine blades after grinding. Firstly, the structure of the grinding robot is described, and the method for planning the motion of the robot is designed. Then, the strategies of surface tracking and constant force control are established; fuzzy PID control is applied to the power source of the manipulator, and gravity compensation for the grinding robot is conducted. Finally, the controlling algorithm constructed in this paper is applied to the operation of a grinding robot to validate the reliability of our method of surface tracking and constant force control; by comparison with artificial grinding, the feasibility of the effectiveness of the scheme is confirmed. The experimental results are discussed at the end of the paper.

2. Grinding Robot

2.1. Structure of the Grinding Robot

The grinding robot addressed in this paper is composed of three main parts: an omnidirectional motion chassis, a five-degree-of-freedom manipulator and a flexible grinding device, as shown in Figure 1. The omnidirectional motion chassis is equipped with a five-degree-of-freedom manipulator to enable the maneuverability of the grinding robot, and a flexible grinding device is installed at the end of the manipulator.
The structure of the pose adjustment system of the flexible grinding device is shown in Figure 2. Four small springs are arranged around the directional connecting block in the pose adjustment system to ensure that the flexible grinding device can move around the axial surface of the blade, and a micro displacement sensor is included under the large spring in the center, so the normal force of the grinding device on the blade is governed by measuring the compression distance of this spring.

2.2. Motion Planning and Control Scheme

Unlike the structure of a traditional manipulator, the wind-turbine-blade-grinding robot addressed in this paper is a type of vehicle-mounted combined manipulator. For simplicity, the control of the total motion of the grinding robot is separated into two steps. First is the position adjustment of the whole machine; the grinding robot is moved to the correct area for blade processing in this stage. Here, only the control of the moving chassis is considered, while the control of the manipulator is not considered. Second is the grinding stage, during which the mechanical arm maneuvers the front-end grinding device to carry out the grinding of the blade. At this stage, the moving chassis is stationary, and only the control of the motion and force of the mechanical arm are considered, so as to ensure the surface tracking of, and the control of force applied to, the blade. After completing a given grinding area, the above stages are repeated in the next area. In this paper, a longitudinal-spacing grinding method is adopted, and the planning of the overall motion of the grinding robot is illustrated in Figure 3.
In the first stage, only the position of the moving chassis needs to be adjusted, and must be kept at a distance of 1.5~2.5 m from the blade. Figure 4 shows the control method employed in this stage.
In the second stage of the grinding operation, due to the presence of different normal vectors at each grinding point on the outer surface of the wind turbine blade, each part of the grinding robot that can be manipulated needs to be coordinated and controlled to enable pose adjustment. Based on the Peterson equation, the empirical formula of the material removal rate during the grinding process [28,29] is
r = C A K A k t V b V w l w F A
where r is the instantaneous removal rate of material, C A is a constant determined via experimentation, K A is the resistance coefficient, k t is the durability coefficient of the grinding tool, V b is the cutting speed, V w is the feed rate, l w is the width of the grinding area and F A is the grinding force. The factors that most determine the impact of the grinding robot on the blade are the grinding force and the feed rate, so these two factors should be controlled in all following steps. The control strategy applied is shown in Figure 5, and this strategy enables real-time surface tracking and the constant control of the grinding force.
In summary, the control scheme applied to the grinding robot, in which the moving chassis, the manipulator and the front-end grinding device are hierarchically controlled, is shown in Figure 6. The commands controlling the chassis are sent via a remote controller, and each wheel group is controlled separately. Data on the distance and angle of movement are collected by the encoder and angle sensors and are sent back to the host computer. Information on the positioning of the manipulator and the grinding device is collected by multiple sensors and transmitted to the host computer. The two laser sensors installed on the upper and lower edges of the grinding device ensure that the tangent of the grinding roller and blade is perpendicular to the normal vector of the grinding point of the blade, and the displacement sensor included in the pose adjustment system described above controls the pressure of the grinding device on the blade. Cable sensor 1 on the lifting part of the manipulator, the angle sensor incorporated in the angle adjustment mechanism, and cable sensor 2 on the telescopic part jointly control the grinding force and speed of the robot.

3. Control Strategy of Grinding Robot

In the above analysis, the grinding robot is divided into two parts: the manipulator and the moving chassis, the latter of which can be controlled by the human operator. Therefore, control of the chassis is not addressed in detail in the subsequent discussion, wherein we consider the real-time tracking of the blade’s surface and the maintenance of consistent grinding via the control of the manipulator.

3.1. Surface Tracking Control Strategy

3.1.1. Surface Tracking of Blade Surface

Due to the changing curvature and irregular shape of wind turbine blades, it is very important for a grinding robot to track its position on a blade’s surface during grinding, and so it is necessary to establish a reasonable control strategy [30]. Two PANFEE L2S laser ranging sensors are included above and below the front-end polishing device to track the blade surface, as shown in Figure 7a.
Here, the parameters L 1 and L 2 represent the distance between the grinding device and the surface of the wind turbine blade as measured by the laser sensors above and below the grinding device, respectively; D represents the distance between the two laser sensors, and θ represents the angle of offset between the attitude of the grinding device and the deviation of the blade, the calculation formula of which is
  θ =   arctan ( L 2 L 1 D )
After calculating θ , the pitch angle of the grinding device can be adjusted using fuzzy PID control until θ reaches zero and the parameters L 1 and L 2 are equal, which ensures that the grinding device is parallel to the tangent line of the blade’s surface at the point of grinding, and this process is carried out in real time during grinding.
In addition, when the two relevant degrees of freedom—of the three degrees of freedom of lifting, stretching and angle adjustment in the manipulator—are determined, the end of the grinding roller can be moved along any trajectory. Therefore, in order to optimize the control algorithm, the position surface tracking of the robot is planned according to the three limit conditions in the grinding process, as shown in Figure 7b.
Firstly, when the grinding area is the trailing edge of the blade, the grinding surface is almost an inclined plane and the curvature change is comparatively small. Therefore, combining the two degrees of freedom of lifting and stretching, the robot can realize the grinding work at the trailing edge of the blade.
Secondly, when the grinding area is the middle shell, the curvature of the blade begins to change significantly, and the effect of the telescopic manipulator is small, so the degrees of freedom of lifting and angle adjustment are used at this time.
Finally, when the grinding area is the leading edge of the blade, the surface shape is convex, and it is also the bottom end of a single grinding process. At this time, the lifting degree of freedom no longer applies; in consequence, combining the two degrees of freedom of stretching and angle adjustment, the grinding action can satisfy the surface of the leading edge of the blade.
In this way, the tracking of blade surface in a single grinding process can be realized.

3.1.2. Adjustment of the Robot’s Position

To ensure the safe operation of the grinding robot and improve the overall performance, the grinding robot must keep a safe vertical distance from the blade (set at about 1.5~2.5 m in this paper), which is measured using laser sensor 3, shown in Figure 8. In order to prevent the lift height of the robot’s grinding device from exceeding its stroke distance, the top edge of the blade should be monitored to prevent the grinding device from sticking and causing failure. In this paper, two laser sensors are included at different heights on the lifting manipulator. Laser sensor 3 is set horizontally to detect the horizontal distance between the grinding robot and the wind turbine blade. Laser sensor 4 is attached at the top of the grinding robot to detect the distance between the robot and the highest point of the blade.
Here, d 1 is the vertical distance between the grinding robot and the real-time position of the blade, β is the angle of laser sensor 4, d 2 is the distance measured by laser sensor 2, and d x is the distance as roughly calculated from angle β and d 1 , the relevant formula for which is
  d x = d 1 sin β
d x as calculated by Equation (3) is compared with the distance d 2 measured by laser sensor 4. If d 2 is much larger than d x , the grinding robot has reached the highest edge of the wind turbine blade; when d 2 gives the distance from laser sensor 2 to the ground, then the lift manipulator is commanded to stop moving and rise no further.
Figure 9 shows the control process of ensuring the grinding robot’s safe positioning. Firstly, it is determined whether the vertical distance d1 between the robot and the blade, as measured by laser sensor 3, is within the safe distance of 1.5~2.5 m. If it is within this range, the chassis stops moving. Then the distance d 2 is compared with d x , and if d 2 d x , then the lift manipulator stops rising; otherwise, it commences the grinding action. If d 1 1.5   m , the robot is alerted and will stop moving forward. If d 1 2.5   m , the robot will continue to move closer to the wind turbine blade until it reaches a safe distance.

3.2. Constant Force Grinding Control Strategy

According to the above analysis, the strategy used for controlling the grinding force and feed speed of the device is the key to realizing the constant control and compliance of the force of grinding. There are five degrees of freedom between the manipulator and the grinding device of the robot. The synchronous control strategy is illustrated in Figure 10, the ideal output of which is the hydraulic cylinder that controls the driving manipulator and the servo motor that drives the grinding device. In order to improve the response speed of the hydraulic system, fuzzy PID is used to control each hydraulic cylinder in the manipulator.
In the control strategy, the input variable used to control the grinding device is its tangential angle from the surface being processed, and the grinding speed and direction of the grinding force are the final outputs. The input variables used to control the grinding force are the blade parameters and the preset grinding force, and the final output is the grinding speed and the grinding force. The control of the input variable related to the manipulator is the final output of the former two inputs; the controller outputs the running speed of each part of the manipulator in order to maintain the grinding force and grinding speed required. Using the methods described above, the relevant control strategies for adjusting the grinding speed and the direction of grinding force are adjusted. Figure 11 shows the analysis method used to control the grinding speed and the grinding force.
Here, v H and v g are the lifting speed of the lifting arm and the telescopic speed of the telescopic arm, respectively, and v α is the speed of the hydraulic cylinder used for angle adjustment, which is not equal to the actual angular rotational velocity of the manipulator. Their corresponding relationship is
  d α d t L 4 = v α sin γ
where d α d t is the angular velocity of the opening and closing of the mechanical arm and represents the actual speed of the grinding device; L 3 and L 4 , respectively, represent the distances between the installed position of the hydraulic cylinder used for angle adjustment and the rotational center of the manipulator being adjusted; α is the angle between the lifting and telescopic manipulators; and γ is the angle between the telescopic manipulator and the hydraulic cylinder used to adjust the angle. α can be measured in real time by the angle sensor, and γ can be obtained via the following relationship:
γ   =   arctan ( L 3 sin α L 4 L 3 cos α )
Combining Equations (4) and (5), to meet the required opening and closing speed of the manipulator, the speed provided by the hydraulic cylinder for angle adjustment must be
v α = d α d t · L 4 sin [ arctan ( L 3 sin α L 4 L 3 cos α )   ]
The grinding device adjusts the pressure applied on the blade by changing the displacement l, which is measured by the micro displacement sensor installed in the pose adjustment device. The resultant vector of the differentiated feed speed v F and the stretching speed of the grinding device must both be adjusted in real time to ensure the correct grinding force, as shown by the following equation:
v F i n t e g r a t i o n = v F + v = d l d t + v
The mathematical relationship between the coordinate systems of the grinding robot and those of each joint is essentially nonlinear, and can be expressed by the following nonlinear vector value function:
x ( t ) = f [ q ( t ) ]
So, the operating speed ultimately required is given by
x ˙ = [ x 1 ˙ , x 2 ˙ , x m ˙ ] T = v i n t e g r a t i o n
To provide the speed ultimately required, the speed required from each joint is
q ˙ = [ q ˙ 1 , q ˙ 2 , q ˙ n ] T
The relationship between the required speed of the end and the speed of each link can be obtained by applying a kinematic calculation to the robot, which is shown in the following expression:
x ˙ = J q ˙
where J is the Jacobi matrix of the robot, the specific value of which can be derived from the kinematics of the grinding robot, as shown in the following expression:
J T ( q ) = [ J T 1 ( q )   J T 2 ( q )   J T 3 ( q )   J T 4 ( q )   J T 5 ( q ) ]
The elements in the formula are
J T 1 ( q ) = [ 0 ( d 2 d 4 c θ 3 ) ( c θ 5 s θ 3 c θ 3 s θ 5 ) + ( d 4 s θ 3 + a 2 ) ( s θ 3 s θ 5 c θ 3 c θ 5 ) d 4 s θ 3 + a 2 c θ 5 s θ 3 c θ 3 s θ 5 s θ 3 s θ 5 c θ 3 c θ 5 0 ] J T 2 ( q ) = [ 0 0 1 0 0 0 ] J T 3 ( q ) = [ 0 d 4 s θ 5 0 s θ 5 c θ 5 0 ] J T 4 ( q ) = [ 0 0 1 0 0 0 ] J T 5 ( q ) = [ 0 0 0 0 0 1 ]
The above-designed control strategy ensures that the correct speed is provided by each connecting rod at any given time in order to meet the grinding requirements, as
q ( t ) = J 1 ( q ) q ˙ ( t )
In order to resolve the problem whereby the inverse Jacobian matrix of the grinding robot does not exist because of the redundancy of the mechanical arm, the following expression can be used [31]:
C = 1 2 q ˙ T A q ˙ + λ T [ x ˙ J ( q ) q ˙ ]
where C is an m × m symmetric positive definite matrix, and λ is the cost judgement. By minimizing C , the speed of each link and λ, we can obtain the following expression:
q ˙ ( t ) = A 1 J T ( q ) λ
q ˙ ( t ) = J ( q ) q ˙ ( t )
By combining the above two expressions, λ can be derived as
λ = [ J ( q ) A 1 J T ( q ) ] 1 x ˙ ( t )
By combining Equations (16) and (18), the desired speed x ( t ) is given for the end of the grinding device, and the speed of each joint is obtained according to Equation (18):
q ˙ ( t ) = A 1 J T ( q ) [ J ( q ) A 1 J T ( q ) ] 1 x ˙ ( t )
where A is an identity matrix. Equation (19) can be simplified to Equation (14), which relates the grinding blade’s position and time. The position of the end of the grinding device, which is related to the sample taken at any given time, is a fixed value, and the sampling period is 200 ms.
Due to the redundancy of the grinding robot, there are multiple solutions to determining the expected speed of the robot end at any time. In order to ensure optimal control of the robot, the change in speed of the connecting rod must be as small as possible in any sampling period. Using the sensor and a calculation of the mathematical relationship, the angle between the desired direction of travel at the end and the direction of each connecting rod at any given time can be obtained, as shown in Figure 12.
Give its three-part structure, the robot can be moving in two directions at any given time. For example, the lifting system can be raising or lowering at any time. In order to ensure its compliance to control, the angle between its direction of motion and the desired direction of the end must be as small as possible [32]. The speeds required of each connecting rod are calculated by Equation (19). On this basis, the required driving force, or torque, which is also the force provided by the hydraulic cylinder at any time, can be obtained by the following expression:
f 4 = D 44 q ¨ 4 ( t ) + D 45 q ¨ 5 ( t ) + D 433 q ˙ 3 ( t ) q ˙ 4 ( t )
where D i j k ( q ) is the inertial matrix of the grinding robot.

3.3. Control Strategy of Mechanical Arm Power System

In maintaining constant control over the force of the grinding robot, the position of the manipulator will change at any given time, so the PID parameters of the hydraulic control system need to be dynamically adjusted in real time. Based on the necessity of real-time surface tracking of the turbine blade and the application of constant force during grinding, the fuzzy PID control principle used to control the hydraulic system is shown in Figure 13.
Through the analysis of the grinding control system, it is finally determined that the inputs of the fuzzy PID control are the deviation and deviation increment of the combined vector of the expected grinding feed speed and the expected grinding force and the combined vector of the two vectors at any time, which is used as the input to control the adjustment of the fuzzy parameters Kp, Ki and Kd in real time.

3.3.1. Fuzzy Input Value

In order to achieve fuzzification, it is necessary to fuzzify the calculated values. In this design, the fuzzy subset is {NB, NM, NS, ZO, PS, PM, PB}, the domain is {−6, −5, −4, −3, −2, −1, 0, 1, 2, 3, 4, 5, 6}. There is an interval for the measured force and velocity signals, which is made as same as the PID adjustment time. Deviation and deviation change rate are shown as:
f ( e ) = 6 · e V m a x V m i n
f ( e ) = 6 · e c 2 ( V m a x V m i n )
where e is the deviation and ec is the deviation increment. The quantization of input error and error change rate can be obtained by using the above equation. The functional relationship of the linear membership function is shown in Figure 14.

3.3.2. Establish Fuzzy Rule Table

After modeling the input, Kp, Ki and Kd in the PID control are adjusted. The specific fuzzy rules are shown in Table 1. According to the fuzzy rule table and the deviation and its change rate calculated by the value collected by the sensor and the expected value, the fuzzy subsets corresponding to ∆Kp, ∆Ki and ∆Kd are obtained.

3.3.3. Demoulding Processing

Since Kp, Ki and Kd are adjusted in the actual fuzzy control, it can be quantified by using the center of gravity calculation according to the triangular membership function, as shown in Equation (23).
v 0 = i = 0 n M i F i i = 0 n M i
where M is membership and F is fuzzy quantization value. The sum of membership is 1, so the denominator can be omitted, so that the calculation of each object is actually a matrix operation:
K = [ M e 1   M e 2 ] [ F a   F b F e   F d ] [ M e c 1   M e c 2 ] T
Coefficients can also be introduced to enlarge and reduce the variation of Kp, Ki and Kd after the increment is obtained.
K ( n ) = K ( n 1 ) + Δ K · λ
where λ is set to 1 in this control strategy and Δ K is the calculated value. The hydraulic cylinder driving the lifting manipulator is taken as an example for this experimental analysis. The performance indices of the hydraulic cylinder with fuzzy PID control and without fuzzy PID control are compared, as shown in Figure 15, and their specific indices are shown in Table 2.
It can be seen that, when using fuzzy PID control, the hydraulic cylinder shows a small degree of overshoot and fast response, and the errors arising from the designed static control strategy can be greatly reduced.

3.4. Analysis of Gravity Compensation

The grinding device of the robot described in this paper is attached to the end of the manipulator by a fixed connection, and a schematic diagram of the robot’s coordinates of movement is shown in Figure 16. In addition to the parameters relevant to the control strategy, the factors influencing the stability of the control strategy also include the contribution of the grinding device’s own gravity to the force applied.
The overall force includes the contact force F E , the gravity of the grinding device G T and the inertial force of acceleration during the movement of the manipulator F I , so the force measured by the sensor cannot be directly used in the algorithm. The influence of inertial force can be ignored because the acceleration of the manipulator is negligible, and gravity compensation can be incorporated to eliminate the influence of the grinding device’s gravity. Figure 17 shows the analysis of the force measured by the sensor in the grinding device, and the relation is shown in Equation (26).
F S = F E + G T + F I
This can be simplified as
F S = F E + G T
Using the rotation transformation matrix of the robot, the force of the gravity of the grinding device as measured by the sensor can be expressed as
G T S = R E S R O E G T O
where G T S is the gravitational force vector of the front-end grinding device in the coordinate system given by the displacement sensor, R E S is the rotation matrix derived from the system of axes of the terminal part of the robotic arm in relation to the sensor’s coordinate system and R O E is the rotation matrix of the frame of the axes of the grinding device in relation to the basal coordinate system, the relative expression of which is
O E R = c θ 1 c θ 3 c θ 5 + s θ 3 s θ 5 c θ 1 s θ 3 s θ 5 c θ 3 s θ 5 s θ 1 s θ 1 s θ 3 c θ 5 + s θ 3 s θ 5 s θ 1 s θ 3 c θ 5 c θ 3 s θ 5 c θ 1 s θ 3 c θ 5 c θ 3 s θ 5 s θ 3 s θ 5 c θ 3 s θ 5 0
where θ i (i = 1~5) is the angle of the joints between mechanical arms I to V. There is a fixed deflection θ between the fixed system of axes of the grinding device and the system of axes of the sensor used to measure the pressure, as shown in Figure 18.
The rotational transformation of the system of axes of the sensor relative to the flange’s coordinate system is
R E S = [ c θ s θ 0 s θ c θ 0 0 0 1 ]
Combining the above four formulas, the following expression can be obtained:
[ F x F y F z ] = c θ s θ 0 s θ c θ 0 0 0 1 · c θ 1 c θ 3 c θ 5 + s θ 3 s θ 5 c θ 1 s θ 3 s θ 5 c θ 3 s θ 5 s θ 1 s θ 1 s θ 3 c θ 5 + s θ 3 s θ 5 s θ 1 s θ 3 c θ 5 c θ 3 s θ 5 c θ 1 s θ 3 c θ 5 c θ 3 s θ 5 s θ 3 s θ 5 c θ 3 s θ 5 0 · 0 0 G T
Simplifying Equation (31), the following expression can be given:
c o s θ = F x 2 + c o s 2 ( θ 3 + θ 5 ) c o s 2 ( θ 3 + θ 5 ) + s i n 2 ( θ 3 θ 5 )
Multiple sets of data pertaining to the installation deflection angle θ can be derived from information on the angles of the robot’s different poses, and an average value can be obtained. Because the micro displacement sensor is used to collect force information, the force that is measured is only unidirectional, and its calculation formula is
F s E = F G T c o s ( θ 3 θ 5 ) c o s ( θ + θ 1 ) / c o s 2 θ 5 t a n 2 θ 5 + ( c o s 2 θ 3 + s i n   θ 3 s i n   θ 1 ) + c o s   θ 3 t a n   θ 5 ( s i n   θ 3 s i n   θ 1 )
where F s E is the grinding force as a component of the coordinate system measured by the pressure sensor, and F is the grinding force calculated by the displacement sensor. In the process of grinding a given area, the rotational angle of the grinding robot is always ±90°, so Equation (34) can be simplified as
F s E = F G T c o s ( θ 3 θ 5 ) s i n θ / c o s 2 θ 5 t a n 2 θ 5 + c o s 2 θ 3 + c o s   θ 3 t a n   θ 5 s i n   θ 3
The operational flow of gravity compensation in the grinding robot is as follows:
  • The grinding device performs a profiling movement on the blade’s surface when it is not in contact with the blade and the grinding force is zero, after which time the measurements of the pressure sensor and the manipulator in each attitude are collected;
  • The force and attitude information collected are substituted into Equation (33) to obtain the installation deflection angle;
  • The above two steps are repeated to collect multiple different sets of information and thus obtain multiple installation deflection angles, which can be used to calculate the average installation deflection angle θ ;
  • The deflection angle obtained in the previous step is substituted into Equation (34) to determine whether the result is 0; if it is not 0, the above steps are repeated until the contact force obtained is 0.
In this way, gravity compensation can be realized; the influence of the gravitational force of the grinding device can be eliminated using the above formula, and the real contact force can thus be obtained, which can then be applied in the constant force control algorithm.

4. Experiments

To verify the feasibility of our control strategy, the surface tracking and constant force control of the grinding robot have been experimentally verified in a grinding operation. The experimental results of robotic grinding are here compared with the artificial results, and this verifies the feasibility of the designed control strategy, according to which the grinding robot can meet the requirements of a given processing site.

4.1. Control Algorithm Comparison Experiment

Taking the driving hydraulic cylinder at the opening and closing manipulator as an example, the experimental analysis is carried out. The performance indices of hydraulic cylinders using fuzzy PID control and normal PID control when the length of the telescopic arm is 218.5 mm are compared as shown in Figure 19. The maximum overshoot is 0.27 rad/s.
The reason for the large difference between the two response curves is that the opening and closing angle corresponds to the expansion length of the hydraulic cylinder. Normal PID has different control effects when the hydraulic cylinder telescopic length is different. However, the three parameters Kp, Ki and Kd of fuzzy PID are variable, which can adapt to the different lengths of hydraulic cylinder. The experiment shows that fuzzy PID control is more adaptable and superior than normal PID control.

4.2. Surface Tracking Experiment

This surface tracking experiment uses time intervals of 200 ms. The coordinate sizes of the tangent point of the grinding device and the blade relative to the base coordinates of the grinding robot are calculated using data given by the five laser sensors installed on the manipulator, enabling us to monitor the actual trajectory of the blade.
In a given grinding area, the grinding robot does not change direction on the x-axis, and this value is always 0. The formulae for calculating the coordinates of the tangent point of the robot and the wind turbine blade in the y-direction and z-direction are
y   = d 4 sin ( π θ 3 180 ° ) + 600 cos ( π θ 3 + θ 5 180 ° 180 ° )
z   = d 2 d 4 cos ( π θ 3 180 ° ) + 600 sin ( π θ 3 + θ 5 180 ° 180 ° )
where d 2 is the position of the lifting manipulator, d 4 is the telescopic length of the telescopic manipulator, θ 3 is the opening and closing angle of the mechanical arms, and θ 5 is the pitch angle of the polishing device. The calculated trajectory of the cutting point coordinates in a grinding area is shown in Figure 20, and they are shown to be very close to the contours of the processing area, which verifies the utility of this control strategy for surface tracking during grinding.

4.3. Grinding Force Tracking Experiment

The expected values of the grinding force and feed speed are set to 1500 N and 150 mm/s, respectively, and these are compared with data collected for a grinding area at sampling periods of 200 ms. Other factors affecting the grinding are set to constant values, such as the revolution speed of the grinding device, which is 1500 r/min. Figure 21 shows the grinding contact force and feed speed, collected from top to bottom during the grinding of a processing area. The grinding contact force is maintained at about 1500 N, and the grinding feed speed remains basically stable at about 150 mm/s after contact, which verifies the utility of our strategy.

4.4. Field Processing Experiment

Actual processing experiments of the grinding robot also needed to be carried out to verify that it could meet the needs of field processing. In the artificial grinding area and the robot grinding area, two square areas with side lengths of 500 mm were randomly selected. Within these areas, detection points were set every 100 mm along the axial direction of the grinding device and the grinding feed direction to measure the roughness of the machined surface, and the measured data are shown in Figure 22. The average roughness yielded by artificial grinding was 5.0934 μm. The average roughness yielded by robot grinding was 2.4828 μm, and all measurement points yielded data that met the surface roughness requirements for an optimal subsequent painting process.

5. Discussion

In this paper, we introduced a robot that can be used for grinding the surface of wind turbine blades, and we proposed related control and motion planning algorithms; we found that tracking the blade’s surface and employing the constant force control algorithm improved the grinding effect. Fuzzy PID control of the power source, focusing on the hydraulic cylinder of the manipulator, was carried out, and this greatly improved the response speed and reduced errors in the static control strategy. In order to eliminate the impact of the weight of the grinding device on the information collected regarding grinding force, a gravity compensation method for the grinding device was designed.
The experiments show that fuzzy PID has excellent adaptability and can meet the design requirements. The effect of tracking the blade’s surface during the grinding process is very impressive, enabling the average grinding force and grinding feed speed to be maintained at near the predicted values. In the actual processing experiment, the surface qualities of blades after robot grinding and artificial grinding were compared. The mean roughness of the blade after grinding by robot was Ra = 2.4828 μm, which is within the optimal range of Ra = 2.0~3.0 μm, and this makes the quality of the blade surface benefit the subsequent painting process.
Future work should study the joint motion of the chassis and the manipulator of the grinding robot. In addition to the grinding force and feed speed, other factors that influence the grinding effect, such as the grinding speed and the amount of sandpaper in the grinding device, can also be further studied.

Author Contributions

Conceptualization, X.S.; methodology, X.S. and S.F.; software, S.F.; validation, M.L., Y.D. and S.F.; formal analysis, Y.D.; investigation, M.L.; resources, X.S.; data curation, Y.D.; writing—original draft preparation, M.L. and S.F.; writing—review and editing, X.S., M.L. and Y.D.; visualization, S.F.; supervision, X.S.; project administration, X.S.; funding acquisition, X.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dai, J.C.; Yang, X.; Wen, L. Development of wind power industry in China: A comprehensive assessment. Renew. Sustain. Energy Rev. 2018, 97, 156–164. [Google Scholar] [CrossRef]
  2. Liu, Y.D.; Wang, J. The SWOT Analysis and Countermeasure Research on the Development of Wind Power Industry in China. In Proceedings of the IOP Conference Series: Earth and Environmental Science, Chengdu, China, 9–11 July 2021. [Google Scholar]
  3. Chen, G.W.; Huang, X.M.; Zhang, L.A.; Wang, J.H.; Liu, W.S. Study on sand grinding technology of wind turbine blade surface. Acta Energ. Sol. Sin. 2021, 42, 302–307. [Google Scholar]
  4. Yan, C.; Chen, X.L.; Li, G.L.; Li, C.L. Design and development of an efficient and adaptive grinding head system for a new type of wind turbine blade automatic grinding robot. Compos. Sci. Eng. 2021, 12, 114–119. [Google Scholar]
  5. Xu, Z.L.; Lu, S.; Yang, J.; Feng, Y.H.; Shen, C.T. A wheel-type in-pipe robot for grinding weld beads. Adv. Manuf. 2017, 5, 182–190. [Google Scholar] [CrossRef]
  6. Ma, W.C. Design and Research on Constant Force Control Device for Robot Grinding. Master’s Thesis, Huazhong University of Science & Technology, Wuhan, China, May 2019. [Google Scholar]
  7. Xiao, M. Research on Constant Force Control Methods in Robot Grinding Process. Ph.D. Thesis, South China University of Technology, Guangzhou, China, May 2020. [Google Scholar]
  8. Zhang, T.; Yu, Y.; Zou, Y.B. An Adaptive Sliding-Mode Iterative Constant-force Control Method for Robotic Belt Grinding Based on a One-Dimensional Force Sensor. Sensors 2019, 19, 1635. [Google Scholar] [CrossRef] [PubMed]
  9. Kuo, Y.L.; Huang, S.Y.; Lan, C.C. Sensorless force control of automated grinding/deburring using an adjustable force regulation mechanism. IEEE Int. Conf. Robot. Autom. 2019, 2019, 9489–9495. [Google Scholar]
  10. Shen, Y.; Lu, Y.; Zhuang, C. A fuzzy-based impedance control for force tracking in unknown environment. J. Mech. Sci. Technol. 2022, 36, 5231–5242. [Google Scholar] [CrossRef]
  11. Wahballa, H.; Duan, J.; Dai, Z. Constant force tracking using online stiffness and reverse damping force of variable impedance controller for robotic polishing. Int. J. Adv. Manuf. Technol. 2022, 121, 5855–5872. [Google Scholar] [CrossRef]
  12. Yao, M.Y.; Zheng, C.W.; Wang, Z.H.; Qian, J.Z.; Xi, Q.; Kuang, S.L. Research on the method of skateboard edge grinding by combination of robot and pneumatic constant force grinding device. In Proceedings of the Fourth International Conference on Mechanical, Electric and Industrial Engineering (MEIE2021), Kunming, China, 22–24 May 2021. [Google Scholar]
  13. Sun, M.J.; Guo, K.; Sun, J. Research on robot grinding force control method. Intell. Robot. Appl. 2021, 13014, 821–829. [Google Scholar]
  14. Dai, S.J.; Li, S.N.; Ji, W.B.; Sun, Z.L.; Zhao, Y.F. Force tracking control of grinding end effector based on backstepping + PID. Ind. Robot. Int. J. 2021, 49, 34–46. [Google Scholar] [CrossRef]
  15. Xu, X.H.; Chen, W. Hybrid active/passive force control strategy for grinding marks suppression and profile accuracy enhancement in robotic belt grinding of turbine blade. Robot. Comput. Integr. Manuf. 2021, 67, 102047. [Google Scholar] [CrossRef]
  16. Zhang, T.; Yuan, C.; Zou, Y.B. Online Optimization Method of Controller Parameters for Robot Constant Force Grinding Based on Deep Reinforcement Learning Rainbow. J. Intell. Robot. Syst. 2022, 105, 85. [Google Scholar] [CrossRef]
  17. Zhao, W.; Xiao, J.L.; Liu, S.J.; Dou, S.X.; Liu, H.T. Robotic direct grinding for unknown workpiece contour based on adaptive constant force control and human–robot collaboration. Ind. Robot. 2022, 50, 376–384. [Google Scholar] [CrossRef]
  18. Wang, G.; Deng, Y.; Zhou, H.; Yue, X. PD-adaptive variable impedance constant force control of macro-mini robot for compliant grinding and polishing. Int. J. Adv. Manuf. Technol. 2023, 124, 2149–2170. [Google Scholar] [CrossRef]
  19. Han, J.L.; Shan, X.L.; Liu, H.T.; Xiao, J.L.; Huang, T. Fuzzy gain scheduling PID control of a hybrid robot based on dynamic characteristics. Mech. Mach. Theory 2023, 184, 105283. [Google Scholar] [CrossRef]
  20. Li, D.W.; Yang, J.X.; Zhao, H.; Ding, H. Contact force plan and control of robotic grinding towards ensuring contour accuracy of curved surfaces. Int. J. Mech. Sci. 2022, 227, 107449. [Google Scholar] [CrossRef]
  21. Jia, W.D.; Jiang, Z.F.; Dai, Y. Constant Force Control Method of Grinding Device. Lect. Notes Comput. Sci. 2022, 13457, 378–387. [Google Scholar]
  22. Zhao, X.W.; Lu, H.; Yu, W.F.; Tao, B.; Ding, H. Vision-based Mobile Robotic Grinding for Large-scale Workpiece and Its Accuracy Analysis. IEEE/ASME Trans. Mechatron. 2022, 28, 895–906. [Google Scholar] [CrossRef]
  23. Dai, S.J.; Zhang, W.H.; Ji, W.B.; Zhao, Y.F.; Zheng, H.W.; Mu, J.H.; Li, P.W.; Deng, R.Q. Research on constant force grinding control of aero-engine blades based on extended state observer. Ind. Robot 2022, 49, 1077–1088. [Google Scholar] [CrossRef]
  24. Xiao, M.; Zhang, T.; Zou, Y.B.; Chen, S.Y. Robotic constant force grinding control based on grinding model and iterative algorithm. Ind. Robot 2020, 48, 270–279. [Google Scholar] [CrossRef]
  25. Zhang, X.X. Hybrid Coutour Force/Position Control of Industrial Robot Based on Variable-Gain Iterative Learning. Ph.D. Thesis, China University of Geosciences, Wuhan, China, May 2022. [Google Scholar]
  26. Li, J.; Guan, Y.; Chen, H.; Wang, B.; Zhang, T. Robotic Polishing of Unknown-Model Workpieces With Constant Normal Contact Force Control. IEEE/ASME Trans. Mechatron. 2023, 28, 1093–1103. [Google Scholar] [CrossRef]
  27. Zhang, T.; Yuan, C.; Zou, Y.B. Research on the algorithm of constant force grinding controller based on reinforcement learning PPO. Int. J. Adv. Manuf. Technol. 2023, 1–14. [Google Scholar] [CrossRef]
  28. Peterson, M.B. Design considerations for effective wear control. In Wear Control Handbook; American Society of Mechanical Engineering: New York, NY, USA, 1980; pp. 413–474. [Google Scholar]
  29. Tsai, M.J.; Huang, J.F.; Kao, W.L. Robotic polishing of precision molds with uniform material removal control. Int. J. Mach. Tools Manuf. 2009, 49, 885–895. [Google Scholar] [CrossRef]
  30. Lai, H.S.; Wu, L.; Chen, X.D.; Yang, Z.Q. Research on the force control of the robot and surface tracking with unknown parameters. Manuf. Technol. Mach. Tool 2018, 6, 21–25. [Google Scholar]
  31. Xiong, M.Q. Research on Key Techniques of Robot Force/Position Control Grinding and Polishing. Master’s Thesis, Changchun University of Science and Technology, Changchun, China, June 2020. [Google Scholar]
  32. Zhang, T.; Xiao, M.; Zou, Y.B.; Xiao, J.D. Research on robot constant force control of surface tracking based on reinforcement learning. J. Zhejiang Univ. Eng. Sci. 2019, 53, 1865–1873. [Google Scholar] [CrossRef]
Figure 1. Grinding robot.
Figure 1. Grinding robot.
Sensors 23 04702 g001
Figure 2. Pose adjustment system.
Figure 2. Pose adjustment system.
Sensors 23 04702 g002
Figure 3. Motion planning of the grinding robot.
Figure 3. Motion planning of the grinding robot.
Sensors 23 04702 g003
Figure 4. Strategy of controlling the grinding robot in the position-adjustment stage.
Figure 4. Strategy of controlling the grinding robot in the position-adjustment stage.
Sensors 23 04702 g004
Figure 5. Control strategy of grinding robot applied in the grinding stage.
Figure 5. Control strategy of grinding robot applied in the grinding stage.
Sensors 23 04702 g005
Figure 6. Scheme of control unit.
Figure 6. Scheme of control unit.
Sensors 23 04702 g006
Figure 7. Surface tracking diagram. (a) Surface tracking method; (b) Different grinding areas of blade.
Figure 7. Surface tracking diagram. (a) Surface tracking method; (b) Different grinding areas of blade.
Sensors 23 04702 g007
Figure 8. Detection of the blade’s highest edge.
Figure 8. Detection of the blade’s highest edge.
Sensors 23 04702 g008
Figure 9. Flowchart showing the safe position control of the grinding robot.
Figure 9. Flowchart showing the safe position control of the grinding robot.
Sensors 23 04702 g009
Figure 10. Five-axis synchronous control strategy diagram.
Figure 10. Five-axis synchronous control strategy diagram.
Sensors 23 04702 g010
Figure 11. Diagram of the control strategy for grinding force and speed.
Figure 11. Diagram of the control strategy for grinding force and speed.
Sensors 23 04702 g011
Figure 12. Instantaneous velocity relationship of contact force control strategy.
Figure 12. Instantaneous velocity relationship of contact force control strategy.
Sensors 23 04702 g012
Figure 13. Fuzzy PID control principle.
Figure 13. Fuzzy PID control principle.
Sensors 23 04702 g013
Figure 14. Membership function of fuzzy controller.
Figure 14. Membership function of fuzzy controller.
Sensors 23 04702 g014
Figure 15. The performance indices. (a) with fuzzy PID control and (b) without fuzzy PID control.
Figure 15. The performance indices. (a) with fuzzy PID control and (b) without fuzzy PID control.
Sensors 23 04702 g015
Figure 16. Diagram of relative positional coordinates.
Figure 16. Diagram of relative positional coordinates.
Sensors 23 04702 g016
Figure 17. Analysis diagram of force measured by sensor.
Figure 17. Analysis diagram of force measured by sensor.
Sensors 23 04702 g017
Figure 18. The grinding device and the sensor’s coordinate systems.
Figure 18. The grinding device and the sensor’s coordinate systems.
Sensors 23 04702 g018
Figure 19. Curve of velocity response.
Figure 19. Curve of velocity response.
Sensors 23 04702 g019
Figure 20. Trajectory of grinding cutting point.
Figure 20. Trajectory of grinding cutting point.
Sensors 23 04702 g020
Figure 21. Feedback values of grinding force and grinding feed speed.
Figure 21. Feedback values of grinding force and grinding feed speed.
Sensors 23 04702 g021
Figure 22. Grinding surface roughness.
Figure 22. Grinding surface roughness.
Sensors 23 04702 g022
Table 1. Fuzzy control rules table.
Table 1. Fuzzy control rules table.
Kp Ki Kdec
NBNMNSZOPSPMPB
eNBPB NB PSPB NB NSPM NM NBPM NM NBPS NS NBZO ZO NMZO ZO PS
NMPB NB PSPB NB NSPM NM NBPS NS NMPS NS NMZO ZO NSNS ZO PS
NSPM NB ZOPM NM NSPM NS NMPS NS NMZO ZO NSNS PS NSNS PS ZO
ZOPM NM ZOPM NM NSPS NS NSZO ZONSNS PS NSNM PM NSNS PS ZO
PSPS NM ZOPS NS ZOZO ZO ZONS PS ZONS PS ZONM PM ZONB PB ZO
PMPS ZO PBZO ZO NSNS PS PSNM PS PSNM PM PSNM PB PSNB PB PB
PBZO ZO PBZO ZO PMNM PS PMNM PM PMNM PM PSNB PB PSNB PB PB
Table 2. Comparison of specific indicators of hydraulic cylinders with different control methods.
Table 2. Comparison of specific indicators of hydraulic cylinders with different control methods.
Control ModeMaximum OvershootRise TimePeak TimeSteady State ErrorAccommodation Time
fuzzy PID control3.58%700 ms1500 ms±0.5 cm/s3100 ms
proportional servo valve control5.58%1500 ms2700 ms±2 cm/s5300 ms
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shi, X.; Li, M.; Dong, Y.; Feng, S. Research on Surface Tracking and Constant Force Control of a Grinding Robot. Sensors 2023, 23, 4702. https://doi.org/10.3390/s23104702

AMA Style

Shi X, Li M, Dong Y, Feng S. Research on Surface Tracking and Constant Force Control of a Grinding Robot. Sensors. 2023; 23(10):4702. https://doi.org/10.3390/s23104702

Chicago/Turabian Style

Shi, Xiaohua, Mingyang Li, Yuehu Dong, and Shangyu Feng. 2023. "Research on Surface Tracking and Constant Force Control of a Grinding Robot" Sensors 23, no. 10: 4702. https://doi.org/10.3390/s23104702

APA Style

Shi, X., Li, M., Dong, Y., & Feng, S. (2023). Research on Surface Tracking and Constant Force Control of a Grinding Robot. Sensors, 23(10), 4702. https://doi.org/10.3390/s23104702

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop