Interference-free tool path generation for the robotic polishing of blisks

In the robotic polishing process, the tool and the blisk interfere easily because of the narrow operation space and seriously twisted curved surfaces. Algorithms are proposed to detect and avoid collisions with high efficiency and accuracy. First, the curved surface of the blade is discretized into a set of points, and the collision detection between the tool and the blisk is converted into the calculation of distances between the tool and points on blade surface. Then, the tool axial vector is adjusted with the minimum rotation angle to avoid collision, which reduces the impact on surface profile accuracy after changing the tool postures. The machining quality is finally guaranteed by con-trolling the material removal depth of the polishing process. The proposed method realizes the collision detection and interference avoidance of the blisk polishing effectively, while it also ensures the surface quality of workpiece when adjusting the tool posture. Simulation and experiments are carried out to verify the feasibility and advantages of the proposed method.


Introduction
As the core part of aero engine, the surface quality of the blisk affects the service life and efficiency of the engine directly. Therefore, it is necessary to grind and polish the blisk after milling to improve the surface quality. Due to the complex geometric structure of the blisk, and its characteristics of distortion, deep and narrow flow channel [1], the tool interferes with the blisk easily during polishing, which makes the trajectory design of the blisk polishing compli-cated and challenging. The research of non-interference polishing technology for blisk is significant to improve the processing quality and promote the development of advanced manufacturing technology for aviation equipment [2,3].
Generally, the interferences between the tool and blisk are divided into two types: local interference and global in-terference [4]. Local interference is the interference between the cutting edge and the tool contact point on the ma-chined surface. It occurs when the tool radius is larger than the curvature radius at the contact point of the curved surface. When local interference occurs, the tool removes excessive materials near the contact point, which directly affects the quality of the machined surface. Global interference is the collision between the tool handle and the machining surface or other surfaces in the machining environment, which causes irreversible damage to the blisk and tool. A lot of scholars have conducted in-depth research on local interference, such as Ahmet [5], Lin [6], Fard [7], etc. Many methods based on differential geometry and curvature 2 analysis have been successfully applied in detection and correction of local interference for complex surface processing. For global interference, finding the best pose where tools do not collide with the workpiece is a process of repeated inspections and corrections, which requires a large amount of calculation [8]. Therefore, the handling of global interference is much more complicated and difficult than local interference. This article focuses on the research of global interference of the blisk polishing.
In existing studies, different collision detection and avoidance algorithms have been proposed for processing interference problems of curved surface machining. Wang et al. [9] proposed a method of constructing collision-free regions without interference checking. This method searches the critical points between the holder and surface and calculates the corresponding critical tool vectors. Then, the critical tool vectors are mapped to a two-dimensional space and the collision-free regions are constructed. Ilushin et al. [10] proposed a collision detection method that uses space subdivision techniques and ray-tracing algorithms to derive a highly accurate polygon/surfacetool intersection algorithm. Ahmet et al. [5] proposed a method for determining the detection area and detection point based on the parameter area, and the interference is checked by calculating the maximum inclination of the connection between the tangential contact point and the detection point. Although this method is technically feasible, the calculation is large. Huang et al. [11] constructed a hierarchical OBB bounding box for the tool, which used the separation axis theory to identify the interference point, and realized the collision avoidance by adjusting the parameters of the tool and the tool bounding box vector. However, it does not consider the effect of changing tool posture on the surface quality of the workpiece.
From the above, it can be seen that existing methods of global interference detection and collision avoidance for curved surface do have contribution to check and avoid collision in order to achieve continuous and safe machining. However, very few methods consider the influence of the tool posture changes on the surface quality of the blisk. This paper proposes a method for generating a non-interference tool path by changing the tool axial vector to a safe position with the minimum rotation angle in the plane where the maximum interference is located, which ensures the machining quality by controlling the material removal depth of the polishing.

Collision detection method
Path planning methods for complex curved surfaces include isoparametric curve method [12], section plane method [13,14], constant scallop-height method [15] and projection method [16], etc. The collision detection and avoidance methods in this article are not limited to the initial path planning method, i.e., different path planning method can be selected according to needs. This paper uses the isoparametric method to generate the initial tool path of the blade. Initial tool pos-

Discretization of blade surface
The collision detection between the tool and the blisk is an essential problem of judging whether the two sets of curved surfaces intersect. However, the calculation of collision between the surfaces is less effective and more difficult. The polishing tool is a standard cylinder, if the surface is discretized into a point set, the collision detection between the tool and the blisk is transferred to calculate the distance from the discrete point to the tool axis, which will improve the calculation efficiency.
When processing a certain blade surface, the tool will only interfere with two adjacent surfaces of the same flow channel, so only the two surfaces need to be discretized. B-spline surface is adopted to construct the blade surface, which is defined as

Simplification of discrete points
In theory, the higher resolution the surfaces are discretized, the higher accuracy the collision detection can reach. However, as the discrete resolution increases, the number of discrete points will increase. If collision detection is performed on all discrete points, the calculation load will be heavy. In order to reduce the calculation burden, it is necessary to exclude discrete points that are impossible to interfere with the tool, in order to reduce the number of discrete points for collision detection.
As shown in Fig.1, the local workpiece coordinate system is established on the polishing blade. The tool handle is projected to the surfaces of two blades along the negative and positive directions of the Z axis in the workpiece coordinate system. Only discrete points located inside the tool handle projection may interfere with the tool.
As shown in Fig 3, the tool is projected to the surface of the polishing blade along the negative direction of the Z axis.
whereU is the projected unit vector of the tool axis vectorT on the XOY plane of the workpiece coordinate system, and K is the unit vector along the positive Z axis of the workpiece coordinate system.
In order to judge whether the discrete points are in the projection of the tool handle, temporary coordinate system WCS  is created whose origin is the tool position point and the coordinate axis directions are the vectors of B , U , K respectively. The transformation matrix from the workpiece coordinate system to the temporary coordinate system is 0 0 5 The coordinates ,, of the discrete point in the temporary coordinate system is calculated as TT If the position of the discrete point The point i F ′ fall inside the projected section of the tool handle, where 1 L is the width of the polishing head, 2 L is the length of the tool handle, and 2 R is the radius of the tool handle. Only the discrete points fall inside the projected section of the tool handle may cause collision. For the points that are not in the projected section of the tool handle, collision could not happen, and they are eliminated in the following collision detection task in order to reduce the calculation burden.

Collision detection
By comparing the shortest distance from the discrete point to the tool axis with the radius of the tool handle, it is judged whether the tool interferes with the machined surface of the blisk. The detection process repeats until all discrete points are traversed.
where is the angle between By noticing as the safe distance between the tool and the blisk, if 6 L=R 2 +ds (9) where ds is the safe distance between the tool and the blisk.
the tool and the blade collide at this tool position, where is the radius of the tool handle. As shown in Fig 5,

Collision avoidance when the tool interferes with the polishing blade
At the tool location point P where collision occurs between the tool and the polishing blade, all discrete points in the projected section of the tool handle are traversed to calculate the shortest distance a L between the tool axis and the discrete points that have interference with the tool. By noting the coordinates of the discrete point corresponding to the shortest distance as a F , the projected coordinate a F  on the tool axis is whereT is unit vector of the tool axis.
where ds is the safe distance between the tool and the blisk.
The tool can be adjusted to a safe orientation with the minimum rotation angle when rotating the tool in the plane M. The minimum rotation angle of the tool axis is calculated by the following formula.
θ =arcsin The change of the tool posture has the least impact for the blade surface machining quality when rotating the tool shaft with the minimum rotation angle .The new interference-free tool axis vector ( ) Changing the tool axis vector will change the bonding state of the tool's grinding head and the blade surface, which causes excessive polishing and affect the surface quality of the workpiece. Therefore, further processing is required to meet the machining accuracy requirement. Whether the machining accuracy meets the requirement is checked by calculating the material removal depth h of the blade surface after rotating the tool. 22 1 where 1 L is the width of the polishing head, 1 R is the radius of the polishing tool head, h is the ma- 8 terial removal depth of the blade's surface after machining, and  is the error constraint value.
If the material removal depth does not meet the requirement of machining accuracy, the tool needs to move a certain distance along the normal vector of the tool contact point to satisfy the requirement. After that, the tool axis vector remains unchanged, and the new tool position point P ′ is calculated by the following formula.
Pꞌ=P+h×N (15) where P ′ is the final interference-free tool position, ( )

Collision avoidance treatment when the tool interferes with the adjacent blade
If interference occurs between the tool and the adjacent blade, the tool rotation angle is calculated according to Eq. (11) and Eq. (12), and then the tool is lifted according to Eq. (14) and Eq. (15). The lifting direction is opposite to the rotating direction relative to the polishing blade surface. The tool may interfere with the blisk after being lifted. Therefore, it is necessary to perform iterative calculation. The collision avoidance processing operation should be performed iteratively until non-interference occurs between the tool and the adjacent blade after rotating and lifting the tool. The distance between the discrete point a F and the tool axis is 9 L j =R 2 +j×ds (16) where j is the number of collision avoidance cycles, ds is the safe distance between the tool and the blisk, and 2 R is the radius of the tool handle. The rotation angle of the tool is The follow-up processing is the same as Section 3.1, therefore details of the mathematical calculations are not repeated here.

Simulation
Visual Studio is used to design the simulation environment of the blisk interference-free tool path generation based on the secondary development of NX 12.0. The experiment is divided into three parts, namely initial tool path generation, collision detection and tool posture correction.
Through the secondary development of NX 12.0, the key information such as the normal vector and curvature of the blade surface is extracted, and the tool path is planned on the surface using the isoparametric method to generate the initial tool path as shown in Fig 8 [17]. First, interference detection is performed on the initial tool path by using the proposed method in section 2. As shown in Fig.9, two interference situations involving the interference with the polishing blade and the adjacent blade are detected. If the interference is detected, the minimum adjusting angle of the tool is calculated to adjust the 10 tool to a safe posture by use the proposed method in section 3. Comparisons are performed with the collision avoidance method proposed by Cai [18], which rotate the tool along the normal direction of the interference point to avoid collision. In the tests, there are 22 tool positions that have interference between tool and the polishing blade on one of the tool paths. The rotation angles of the two methods for collision avoidance are shown in Table 1.
Through the comparison, it can be seen that the rotation angles obtained by the proposed method are smaller, which reduces the impact on blade's surface processing quality after changing the tool postures.
After adjusting postures of the tool, the bonding state between the surface of the polishing tool head and the surface of the blade is changed, which causes excessive polishing of the machined surface and affect the processing quality. Taking the above collision points as an example, the material removal depth comparison after rotating the tool posture is shown in Fig 10. It can be clearly seen from Fig.10 that rotating the tool posture will cause over-polish of the blade surface. The maximum material removal depth can reach 0.308mm, which causes over polishing and affects the processing quality of the blisk seriously. After lifting the tool, the over polishing is eliminated effectively. The distance between the polishing head and the blade surface is controlled within 0mm-0.02mm. The force control mechanism can control the grinding head to fit the blade surface for polishing effectively. Finally, an interference-free tool path that meets the processing requirements is generated. In addition, the tool posture before and after collision avoidance treatment is shown in Fig 11.   The rotation angle, the material removal depth before and after treatment during collision avoidance treatment are shown in Table 2. It can be seen that the material removal depths of generated tool pose satisfy the accuracy requirement after adjusting the tool orientation and position. The 12 force control mechanism can control the grinding head to fit the blade surface during the polishing process.

Experiment verification
Comau robot is used to perform interference-free polishing experiments of the blisk in this paper. As shown in Fig 13, Comau robot holds the force control mechanism to polish the blisk. Force control polishing is realized by real-time feedback controlling the positive pressure through the force control mechanism. In experiments, the polishing tool path is generated by using the proposed non-interference tool path generation method, and actual polishing tests are conducted on the blisk. As shown in Fig.14, the interference-free polishing of the blisk is effectively realized in the actual machining process.
The blade surface roughness is analyzed to further verify the machining quality of the proposed method. The surface roughness of the blade before and after polishing is shown in Fig.15. It can be seen that the surface roughness of the blade before polishing is 1.406um and the surface roughness of the blade after polishing is 0.358um. The effect of polishing is obvious. 13 (a)Polishing of the blade convex (b)Polishing of the blade back In summary, it can be seen that the interference-free tool path generation method of the blisk polishing in this paper does not only achieves collision detection and collision avoidance, but also ensures the machining quality of workpiece surface.

Conclusions
In order to address the problem of the interference in the polishing process of the blisk，this paper proposes an interference-free tool path generation method of the blisk polishing. The main conclusions are as follows: (1) The distance between the discrete points on the surface of the blade and the axis of the tool is used to determine whether the tool and the blisk have collision, which improves the efficiency of collision detection.
(2) If the collision is detected, the tool axial vector is adjusted through the minimum rotation angle to avoid collision, which reduces the impact on blade's surface quality after changing the tool postures.
(3) The tool posture after collision avoidance is further adjusted in order to control the material removal depth within the preset limitations, which ensures the quality of workpiece surface.
Simulation and experiments verify that the proposed method realizes the collision detection and avoidance of the blisk polishing effectively. The surface quality of the workpiece is also guaranteed after adjusting the tool posture to avoid interference.