Twin-Tool Orientation Synchronous Smoothing Algorithm of Pinch Milling in Nine-Axis Machine Tools

Pinch milling is a new technique for slender and long blade machining, which can simultaneously improve the machining quality and efficiency. However, two-cutter orientation planning is a major challenge due to the irregular blade surfaces and the structural constraints of nine-axis machine tools. In this paper, a method of twin-tool smoothing orientation determination is proposed for a thin-walled blade with pinch milling. Considering the processing status of the two cutters and workpiece, the feasible domain of the twin-tool axis vector and its characterization method are defined. At the same time, an evaluation algorithm of global and local optimization is proposed, and a smoothing algorithm is explored within the feasible domain along the two tool paths. Finally, a set of smoothly aligned tool orientations are generated, and the overall smoothness is nearly globally optimized. A preliminary simulation verification of the proposed algorithm is conducted on a turbine blade model and the planning tool orientation is found to be stable, smooth, and well formed, which avoids collision interference and ultimately improves the machining accuracy of the blade with difficult-to-machine materials.


Introduction
Pinch milling uses two cutters that can simultaneously mill the opposing surfaces of two blades from top to top, which can significantly improve machining efficiency and quality.Turbine blades are typically elongated components, distinguished by thin walls, asymmetry, and twisted free surfaces.A twin-tool orientation, formed by three rotary axes, plays a crucial role in machining complex parts like a turbine blade.Each tool orientation requires continuous and smooth transitions, but it is very difficult to simultaneously plan these tool orientations.Meanwhile, inappropriately defined tool orientations in such applications can cause fatal collisions and damage the in-processing part, and unsmooth segments of tool orientation may lead to unwanted vibrations due to sudden fluctuations in tool movements and, thus, limit the full production capabilities of multiaxis machine tools and inevitably impact machining precision.
In order to obtain smooth tool orientation, many studies have been conducted to achieve optimal tool orientation.For example, the tool-sweep surface-based method [1], curvature matching method [2][3][4], multipoint machining method [5][6][7], and simultaneous optimization with feed direction [8,9]; all share the goal of selecting a noninterference and collision-free tool orientation.Meanwhile, Mi et al. [10] presented a feasible C-space computation algorithm for triangular mesh models to plan smooth tool orientation along a tool path.Chen et al. [11] presented a reference plane to generate a set of smoothly aligned tool orientations along a tool path.Gong et al. [12] proposed a method to find the optimal tool orientation based on a ruled surface; the optimization objective function was set to minimize the vibration of rotary axes.Hong et al. [13] described a practical method for tool orientation generation for ball cutters commonly used in complex-surface-finish machining.Dong et al. [14] proposed a multiscale tool orientation generation method that considers both the machining strip width and roughness scales.Sun et al. [15] presented the prediction of an automatic tool axis orientation algorithm that avoids chatter while improving productivity.
In recent years, researchers' attention has also been directed to smoothing tool orientation, and researchers have proposed many typical methods, such as the shortest-path linking method [16], interpolation method [17,18], and optimization-based method [19][20][21].Yuan et al. [22] proposed a method for generating smooth tool orientation by considering the relationship between tool orientation and strip width for five-axis machining with flat-end cutters.Dong et al. [23] proposed a multiscale tool orientation smoothness method that considers both machining strip width and roughness scales.On this basis, a novel method based on the best curvature matching was used to generate smoothing tool orientation [24].Yana et al. [25] proposed a mathematical framework to generate smoothing tool axis variation even on partial surfaces lacking G2 and/or G1 continuities.Wang et al. [26] constructed a selection strategy for the smoothing tool axis from the discrete domain of feasible orientations.All the technology deeply improved the smoothness of tool orientation vector sequences.
The above theoretical studies are dedicated to planning smoothing tool axis vectors in single-tool machining.But pinch milling has two cutters, and any nonsmooth tool axis direction can seriously affect the cutting process.In early stages, the twin-tool milling path is planned [27], the opposite cutting contact points and paths are planned, the twintool orientation is characterized, and the initial dual-tool axis vectors are also planned.However, in actual machining, it was found that the angle between the cutters changed discontinuously and unevenly, which limited the cutting speed of the machine tool and caused obvious vibration during cutting processes.Meanwhile, the envelope surface is jointly constructed by the tool and workpiece at the cutting contact, which plays a crucial role in the cutting process.The state of the tool orientation has a significant impact on the formation of the envelope surface and the machining quality of the workpiece surface.It is a key item in the twin-tool cutting process.Therefore, smoothness orientation planning for pinch milling is an enormous challenge.
In this paper, the twin-tool orientation synchronous smoothing algorithm is proposed for thin-wall blades used in pinch milling.Considering the two cutters' structural layouts, the twin-tool orientations and coupling relationship are characterized; a method is also proposed for defining the feasible region of the cutters.In order to guarantee the smoothness of the tool orientations, the tool axis vector is parameterized and tool posture curve of global optimization is formulized, and the evaluation algorithm of local optimization is investigated.Finally, the twin-tool orientation selection and optimization scheme is explored within the feasible domain along the two tool paths.Furthermore, the smoothness tool orientation is planned for a typical turbine blade.

The Method of Pinch Milling
Pinch milling is a method that involves simultaneously milling both sides of an irregular blade profile from top to top, using two tools concurrently.The milling process is shown in Figure 1; two tools are oppositely and simultaneously assigned to mill the dorsal (convex) and basin (concave) surfaces, and the directions along blade height are the feed directions for the two tools.Meanwhile, it considers the shape and size characteristics of the leading and trailing edges in order to avoid collisions between the two cutters, which are milled along the length direction with either one of the twin tools.
the leading and trailing edges in order to avoid collisions between the two cutters, which are milled along the length direction with either one of the twin tools.At the top-to-top contact point, the tool axes orientations of the two cutters should influence each other, and the rotation angle of the workpiece simultaneously affects the cutting angle of the tool.

Twin-Tool Orientations Description
In pinch milling, in addition to the translations along the length direction of the blade, two cutters can be rotated with the workpiece axis to adapt to complex curved surfaces.As shown in Figure 2, the local coordinate system OL-XL-YL-ZL is defined at the cutter contact (CC) point C. The XL-axis is defined along the instantaneous cutting direction, the ZL-axis is defined along the direction of local surface normal, and the YL-axis is defined by the XLand ZL-axes with the right-hand-rule.The tool orientation is determined with an inclination angle αL around the YL-axis and a tilt angle λL around the ZL-axis.The tool orientation for the other tool is also determined by the same rule.In the local coordinate system OL-XL-YL-ZL, the unit vectors f, n, and I are defined along the direction of the XL-axis, ZL-axis, and tool-axis, respectively.Then, the inclination angle αL and the tilt angle λL can be derived as follows: where At the top-to-top contact point, the tool axes orientations of the two cutters should influence each other, and the rotation angle of the workpiece simultaneously affects the cutting angle of the tool.

Twin-Tool Orientations Description
In pinch milling, in addition to the translations along the length direction of the blade, two cutters can be rotated with the workpiece axis to adapt to complex curved surfaces.As shown in Figure 2, the local coordinate system O L -X L -Y L -Z L is defined at the cutter contact (CC) point C. The X L -axis is defined along the instantaneous cutting direction, the Z L -axis is defined along the direction of local surface normal, and the Y L -axis is defined by the X L -and Z L -axes with the right-hand-rule.The tool orientation is determined with an inclination angle α L around the Y L -axis and a tilt angle λ L around the Z L -axis.The tool orientation for the other tool is also determined by the same rule.
Materials 2024, 17, x FOR PEER REVIEW 3 of 14 the leading and trailing edges in order to avoid collisions between the two cutters, which are milled along the length direction with either one of the twin tools.At the top-to-top contact point, the tool axes orientations of the two cutters should influence each other, and the rotation angle of the workpiece simultaneously affects the cutting angle of the tool.

Twin-Tool Orientations Description
In pinch milling, in addition to the translations along the length direction of the blade, two cutters can be rotated with the workpiece axis to adapt to complex curved surfaces.As shown in Figure 2, the local coordinate system OL-XL-YL-ZL is defined at the cutter contact (CC) point C. The XL-axis is defined along the instantaneous cutting direction, the ZL-axis is defined along the direction of local surface normal, and the YL-axis is defined by the XLand ZL-axes with the right-hand-rule.The tool orientation is determined with an inclination angle αL around the YL-axis and a tilt angle λL around the ZL-axis.The tool orientation for the other tool is also determined by the same rule.In the local coordinate system OL-XL-YL-ZL, the unit vectors f, n, and I are defined along the direction of the XL-axis, ZL-axis, and tool-axis, respectively.Then, the inclination angle αL and the tilt angle λL can be derived as follows: In the local coordinate system O L -X L -Y L -Z L , the unit vectors f , n, and I are defined along the direction of the X L -axis, Z L -axis, and tool-axis, respectively.Then, the inclination angle α L and the tilt angle λ L can be derived as follows: where Considering structural constraints of pinch milling, the two tools can only horizontally swing and parallel each other.As Figure 3, in the workpiece coordinate system O W -X W -Y W -Z W , the points C 1 and C 2 are top-to-top cutting contacts, a 1 and a 2 are supposed to be the tool orientations at the opposite cutting points, respectively, n 1 and n 2 are supposed to be the normal vectors at the point on the dorsal and basin surfaces respectively, then m 1 and m 2 are the projection vectors of n 1 and n 2 onto the plane Y W -O W -Z W .

Twin-Tool Orientations Coupling Relationship
Local coordinate systems OL1XL1YL1ZL1 and OL2XL2YL2ZL2 are established at the point C1 and C2 according to the right-hand rule; f1 is the unit vector of tool feed direction at the point C1, and f2 is the unit vector of the tool feed direction at the point C2.If the tool axis Twin-tool cutting path Assuming normal vector n 1 is [x 1 , y 1 , z 1 ] and normal vector n 2 is [x 2 , y 2 , z 2 ], then [0, y 1 , z 1 ] is expressed as the vectors m 1, and [0, y 2 , z 2 ] is expressed as the vectors m 2 .The angle between the vectors m 1 and m 2 is denoted as φ, which can be expressed as The vectors b 1 and b 2 are the projection vectors of a 1 and a 2 onto the plane Y W -O W -Z W , respectively.Suppose that the angle between m 1 and b 1 is φ 1 , and the angle between m 2 and b 2 is φ 2 , then the relationship among φ 1 , φ, and φ 2 can be expressed as follows: Suppose that the angle between a 1 and b 1 is γ 1 , which is the rotation angle from a 1 to b 1 ; Similarly, the angle between a 2 and b 2 is γ 2 , which is the rotation angle from a 2 to b 2 .Then translating the vectors into the same plane X W -O W -Z W , a 1 and a 2 can be calculated and expressed as where

Twin-Tool Orientations Coupling Relationship
Local coordinate systems O L1 X L1 Y L1 Z L1 and O L2 X L2 Y L2 Z L2 are established at the point C 1 and C 2 according to the right-hand rule; f 1 is the unit vector of tool feed direction at the point C 1 , and f 2 is the unit vector of the tool feed direction at the point C 2 .If the tool axis vector a 1 and a 2 are all unit vectors, the inclination angle α L1 and tilt angle λ L1 of the cutter at the dorsal surface can be derived as follows: where Similarly, the inclination angle α L2 and tilt angle λ L2 of the cutter at the basin surface can be derived as follows: ) where It will be seen that the rotation angle φ 1 simultaneously determines the inclination angles α L1 and α L2 , as well as the tilt angles λ L1 and λ L2 of the tool axis vectors, resulting in mutual influence between the two tool axis vectors without a linear relationship.Meanwhile, the sizes of α L1 and λ L1 of the tool orientation a 1 are directly affected by the parameter of the angle γ 1 , and the sizes of α L2 and λ L2 of the tool orientation a 2 are directly affected by the parameter of the angle γ 2 .This phenomenon indicates that there is a mutually dependent coupling relationship between the two tool axis vectors, meaning that when the parameters of one tool axis vector change, it triggers corresponding changes in the parameters of the other tool axis vector.This interrelated state directly influences the machining conditions during the twin-tool cutting process.Therefore, adjustments to the parameters of either tool axis vector require careful consideration of their impact on the overall machining state.

Identification of Twin-Tool Orientation Accessible Region
Due to the structural constraints of the pinch milling, the two tool movement spaces are all imposed limitations, and we thoroughly delved into the mutual influence relationship.Based on the parameters of φ 1 , γ 1 , and γ 2 , we can solve the preliminary feasible space of all tool axis vectors, which all satisfy the structural constraints.This specific space is defined as the feasible domain of Q(α L1 , λ L1 , α L2 , λ L2 ) and can be specifically described as Under the premise of satisfying the structural constraints on pinch milling, due to multiple factors such as tool size, blade, and fixture shapes, the tool axis vector cannot cover all ranges within the feasible domain Q.In addition, more constraints are considered to avoid collisions and interferences.On the top-to-top cutting contact, the parameters are calculated based on tool data and local geometric data of the curved surface, such as the distance d e between the tool and the machine tool or workpiece, the effective cutting radius R e of the tool, the radius of curvature ρ k at the contact point, and the distance d s between the tool bottom and the workpiece.Using these data, we can determine a collision-free or interference-free tool axis vector space.In the twin-tool milling process, any of these tools must be ensured to avoid collisions and interferences.Therefore, the feasible domain is further restricted and marked as R 1 (α L1 , λ L1 , α L2 , λ L2 ) and can be specifically described as Meanwhile, the geometry of the cutting surface is affected by the angle of the tool orientation, which manifests as excessive residual heights, as well as overcutting and undercutting phenomena, thereby affecting the machining quality of the workpiece surface.Thus, under the condition of the allowable scallop height h max and chord error δ tmax requirements, the reachable motion space of the tool axis vector can be determined, which is defined as the feasible region R 2 (α L1 , λ L1 , α L2 , λ L2 ) under geometric error constraints, and is expressed as where h r is the actual scallop height and δ tr is the actual chord error corresponding to the two-tool axis vectors.
The true feasible region is actually the intersection of all feasible regions considered under the constraints, thus forming a closed area reachable by the twin-tool axis vectors, as shown in Figure 4.This area is defined as the feasible region Ω(α L1 , λ L1 , α L2 , λ L2 ) for the posture of dual tools, and is represented as free or interference-free tool axis vector space.In the twin-tool milling process, any of these tools must be ensured to avoid collisions and interferences.Therefore, the feasible domain is further restricted and marked as R1(αL1, λL1, αL2, λL2) and can be specifically described as (10) Meanwhile, the geometry of the cutting surface is affected by the angle of the tool orientation, which manifests as excessive residual heights, as well as overcutting and undercutting phenomena, thereby affecting the machining quality of the workpiece surface.Thus, under the condition of the allowable scallop height hmax and chord error δtmax requirements, the reachable motion space of the tool axis vector can be determined, which is defined as the feasible region R2(αL1, λL1, αL2, λL2) under geometric error constraints, and is expressed as (11) where hr is the actual scallop height and δtr is the actual chord error corresponding to the two-tool axis vectors.
The true feasible region is actually the intersection of all feasible regions considered under the constraints, thus forming a closed area reachable by the twin-tool axis vectors, as shown in Figure 4.This area is defined as the feasible region Ω(αL1, λL1, αL2, λL2) for the posture of dual tools, and is represented as Assuming that the feasible region of a1 is Ω1(αL1, λL1), and the feasible region of a2 is Ω2(αL2, λL2), the final feasible region is the union of the feasible regions of the two tool axis vectors.The equivalent expression Ω is represented as ( , , , )= ( , ) ( , ) In addition, to clearly describe the positions of the dual tool axis vectors in the workpiece coordinate system OWXWYWZW, a Gaussian sphere of unit tool axis vectors is constructed at the origin OW.Meanwhile, the unit vector a1 and a2 are translated to the origin OW, respectively, so that a1 and a2 are equivalent to a particle on the Gaussian sphere, and the corresponding feasible regions Ω1 and Ω2 represent the ranges of motion of the particles on the Gaussian sphere surface, as shown in Figure 5.However, the inclination angle αL and tilt angle λL only describe the tool's posture in the local coordinate system at the contact point, and its tool axis vector is actually a directional vector in the workpiece coordinate system, which can be represented as (i, j, k).Therefore, the feasible region Ω of twin-tool orientation can be re-expressed as Assuming that the feasible region of a 1 is Ω 1 (α L1 , λ L1 ), and the feasible region of a 2 is Ω 2 (α L2 , λ L2 ), the final feasible region is the union of the feasible regions of the two tool axis vectors.The equivalent expression Ω is represented as In addition, to clearly describe the positions of the dual tool axis vectors in the workpiece coordinate system O W X W Y W Z W , a Gaussian sphere of unit tool axis vectors is constructed at the origin O W .Meanwhile, the unit vector a 1 and a 2 are translated to the origin O W , respectively, so that a 1 and a 2 are equivalent to a particle on the Gaussian sphere, and the corresponding feasible regions Ω 1 and Ω 2 represent the ranges of motion of the particles on the Gaussian sphere surface, as shown in Figure 5.However, the inclination angle α L and tilt angle λ L only describe the tool's posture in the local coordinate system at the contact point, and its tool axis vector is actually a directional vector in the workpiece coordinate system, which can be represented as (i, j, k).Therefore, the feasible region Ω of twin-tool orientation can be re-expressed as Materials 2024, 17, x FOR PEER REVIEW 7 of 14 , , ) Figure 5. Feasible region in the workpiece coordinate system.

Global Smoothness of Twin-Tool Orientation
In multiaxis machining, it is expected that the tool axis vectors change continuously and uniformly along the contact points to achieve global smoothing of the tool axis vectors.However, the tool axis vectors formed during trajectory planning may have excessive or uneven angular changes, as shown in Figure 6, and the concept of "smoothness" is inherently complicated and often inconsistent.For unit tool axis vectors, rotation angle of tool orientation is equivalent to the trajectory of a particle moving on the spherical surface.The smoothing is transformed into a curve fitting problem at the particle, which means finding a smooth curve that ensures all tool axis vectors lie on this curve.Assuming the tool axis vector at the cutting contact point is defined as ack, and the values are [ai,ck, aj,ck, ak,ck], where the subscripts i, j, k represent the component values in each axis of the coordinate system, and the subscript ck indicates the position of the cutting contact point), and the corresponding control vertex is ai, then the parametric B-spline curve can be defined by the basis functions Ni,q(s), control points ai, and degree n with the following form: The basis functions Ni,q(s) are functions of the geometric parameter s and knot vector S = [s0, s1, …, sn+1], and are defined as follows:

Global Smoothness of Twin-Tool Orientation
In multiaxis machining, it is expected that the tool axis vectors change continuously and uniformly along the contact points to achieve global smoothing of the tool axis vectors.However, the tool axis vectors formed during trajectory planning may have excessive or uneven angular changes, as shown in Figure 6, and the concept of "smoothness" is inherently complicated and often inconsistent.For unit tool axis vectors, rotation angle of tool orientation is equivalent to the trajectory of a particle moving on the spherical surface.The smoothing is transformed into a curve fitting problem at the particle, which means finding a smooth curve that ensures all tool axis vectors lie on this curve.
Materials 2024, 17, x FOR PEER REVIEW 7 of 14 Figure 5. Feasible region in the workpiece coordinate system.

Global Smoothness of Twin-Tool Orientation
In multiaxis machining, it is expected that the tool axis vectors change continuously and uniformly along the contact points to achieve global smoothing of the tool axis vectors.However, the tool axis vectors formed during trajectory planning may have excessive or uneven angular changes, as shown in Figure 6, and the concept of "smoothness" is inherently complicated and often inconsistent.For unit tool axis vectors, rotation angle of tool orientation is equivalent to the trajectory of a particle moving on the spherical surface.The smoothing is transformed into a curve fitting problem at the particle, which means finding a smooth curve that ensures all tool axis vectors lie on this curve.Assuming the tool axis vector at the cutting contact point is defined as ack, and the values are [ai,ck, aj,ck, ak,ck], where the subscripts i, j, k represent the component values in each axis of the coordinate system, and the subscript ck indicates the position of the cutting contact point), and the corresponding control vertex is ai, then the parametric B-spline curve can be defined by the basis functions Ni,q(s), control points ai, and degree n with the following form: The basis functions Ni,q(s) are functions of the geometric parameter s and knot vector S = [s0, s1, …, sn+1], and are defined as follows:

Initial tool orientation
Cutting path Assuming the tool axis vector at the cutting contact point is defined as a ck , and the values are [a i , ck , a j,ck , a k , ck ], where the subscripts i, j, k represent the component values in each axis of the coordinate system, and the subscript ck indicates the position of the cutting contact point), and the corresponding control vertex is a i , then the parametric B-spline curve can be defined by the basis functions N i,q (s), control points a i , and degree n with the following form:

Anticipant tool orientation
The basis functions N i,q (s) are functions of the geometric parameter s and knot vector S = [s 0 , s 1 , . .., s n+1 ], and are defined as follows: Meanwhile, in order to achieve the analytical solution of control points, the number of basis functions and control points in Equation ( 14) is set equal to the number of tool axis vector a k , which allows for the linear system to be solved for the control points.For this, the knot vector S is solved based upon the angles between the tool axis vector a k .Thus, the parameter values s k are assigned to each tool axis vector a k , which can be expressed as Using the assigned parameter values s k , the knot vector is solved as follows: Based on the defined knot vector S and the assigned parameter values s k , the control point a i can be solved through the line system of equations and represented as Once the knot vector S and control points a i have been determined, the B-spline curve P c (s) of the tool axis can be obtained from Equation (14), which is called the tool posture curve P c (s), and describes the trajectory of a particle on the spherical surface, representing the angle variation of the tool axis vector at the contact point.Therefore, new tool axis vectors are interpolated the contact point along the tool path; all tool axis vectors are located on the tool posture curve, and the global smoothing tool axis vector can be obtained.

Local Smoothness of Twin-Tool Orientation
All cutting contact points are all expected to have tool axis vectors that achieve global smoothness variation.However, the optimized cutting axis vector may not be in the feasible region, resulting in the uncertainty of the local cutting contact, so it is necessary to further solve the local optimization problem of the cutting axis vector.For the discrete data, the variation of the tool axis vector can be described as local angular velocity, which can be approximately represented by mean angular velocity in each interval between two neighboring tool orientations.Assuming that the tool center has a constant feed rate, the average angular velocity at unit speed of the tool axis vector between the ith interval can be expressed as where θ i is the angle between tool axis vectors a i and a i+1 , and ℓ(c i c i+1 ) stands for the arc length between the two cutting contact points.
In real situations, usually, two consecutive cutting contact points are close, and θ i is small, so the approximate relationship between angles and arc length can be expressed as Materials 2024, 17, 2977 9 of 14 Therefore, the angular speed ω i are re-expressed as The smoothness of angular velocity can be achieved by minimizing the harmonic mean of ω i or the almost constant ω i along the tool path to ensure the optimal tool axis vector.Then, the objective function for the optimization problem can be formulated as The optimal tool axis vector is sought to achieve uniform and smooth variation of the tool axis at adjacent consecutive contact points within the feasible region.

Twin-Tool Orientation Selection and Optimization Schemes
As long as the optimization algorithm is clearly defined, the tool posture determination process can be carried out.For the pinch milling, both the orientation of two tools are required to have continuous and smooth changes along the cutting path.Therefore, the selection and optimization process is established to solve the optimal tool axis vector; the flowchart of twin-tool smoothness orientation planning is shown in Figure 7.If a 1 c1k belongs to Ω and a 1 c2k does not belong to Ω, then the objective function H2 is represented as Contact point data file Initial Cutting contact points p c1 、p c2 and tool vector a c1K 、a c2k P c2 (s) For the two cutters on the dorsal and basin surfaces of the blade, all reachable regions Q(α L1 , λ L1 , α L2 , λ L2 ) of the twin-tool axis vector are determined based on the structural constraints.And the feasible domain Ω(α L1 , λ L1 , α L2 , λ L2 ) is calculated based on the cutting contact points p c1 and p c2 on the cutting path and the tool axis vectors a c1k and a c2k .At the same time, the tool posture curves P c1 (s) and P c2 (s) are determined corresponding to the primitive tool axis, and the global optimal vectors a 1 c1k and a 1 c2k are solved.If a 1 c1k and a 1 c2k both belong to Ω, simultaneous optimization vectors on both sides can be achieved, which can be expressed as a 2 c1k and a 2 c2k .However, it is difficult to achieve simultaneous optimization in the actual calculation process.For a 1 c1k and a 1 c2k do not belong to Ω, it is necessary to ensure that the change in angular velocity for each of the cutters is minimized.The objective function H 1 is represented as If a 1 c1k belongs to Ω and a 1 c2k does not belong to Ω, then the objective function H 2 is represented as If a 1 c1k does not belong to Ω and a 1 c2k belongs to Ω, then the objective function H 3 is represented as In the actual optimization process, a 3 c1k and a 3 c2k are defined as the local optimal vectors.If a 1 c1k belongs to Ω, a 3 c1k is the same as a 1 c1k .If a 1 c2k belongs to Ω, a 3 c2k is the same as a 1 c2k , otherwise, they are not the same.Therefore, the final tool axis vector a ′ ck can be expressed as Based on the above algorithm, all the smoothness orientations can be calculated for the top-to-top contact point, and are all optimal tool axis vectors.

Twin-Tool Orientation Smoothness Planning
The implementation of the presented twin-tool orientation optimization algorithm was successfully carried out using C++.A turbine blade is used to test the algorithm, as shown in Figure 8, and the length and rotational diameter of the blade are 633 and 165 mm, respectively.It can be seen from this CAD model that the blade body is thinwall, asymmetric, and torque-shaped; planning the path and tool axis vector for two opposing milling cutters is a complex and challenging task, and, hence, a good test example for the proposed algorithm.The tool path is assumed to be generated from the upper stream procedures.
was successfully carried out using C++.A turbine blade is used to test the algorithm, as shown in Figure 8, and the length and rotational diameter of the blade are 633 and 165 mm, respectively.It can be seen from this CAD model that the blade body is thin-wall, asymmetric, and torque-shaped; planning the path and tool axis vector for two opposing milling cutters is a complex and challenging task, and, hence, a good test example for the proposed algorithm.The tool path is assumed to be generated from the upper stream procedures.Two of the same torus-shaped milling cutters with the radius R of 16 mm and the corner radius r of 6 mm are employed to synchronously mill the oppose surface of the blade, respectively.For the opposing cutting contact points and initial twin-tool axis vectors generated during the path planning of blade profiles, there is no interference during the machining process of the blade surface, and it also satisfies the structural constraints of the two cutters layout.However, there is an obvious nonsmooth phenomenon among the tool axis vectors, and the tool axis vectors are not in a unified sequence.
Using this optimization method, the feasible domain of the two cutters axis vectors are calculated at each cutting contact point, and the tool posture curve is constructed; then, the optimal twin-tool axis vectors are planned through global and local algorithms.The results are shown in Figure 8.By comparing and analyzing the optimization results, it is found that the cutting contact point positions on the path remain unchanged, and the cutting parameters before and after optimization are macroscopically consistent.However, the optimized tool axis vector has visually visible smoothness and changes uniformly along the path, while avoiding interference during the machining process.
The average angular velocities between the initial and optimal tool axis vectors were calculated separately, and the results are shown in Figure 9.The variations in angular velocities are not entirely consistent on both surfaces of the blade, but the significant angle Two of the same torus-shaped milling cutters with the radius R of 16 mm and the corner radius r of 6 mm are employed to synchronously mill the oppose surface of the blade, respectively.For the opposing cutting contact points and initial twin-tool axis vectors generated during the path planning of blade profiles, there is no interference during the machining process of the blade surface, and it also satisfies the structural constraints of the two cutters layout.However, there is an obvious nonsmooth phenomenon among the tool axis vectors, and the tool axis vectors are not in a unified sequence.
Using this optimization method, the feasible domain of the two cutters axis vectors are calculated at each cutting contact point, and the tool posture curve is constructed; then, the optimal twin-tool axis vectors are planned through global and local algorithms.The results are shown in Figure 8.By comparing and analyzing the optimization results, it is found that the cutting contact point positions on the path remain unchanged, and the cutting parameters before and after optimization are macroscopically consistent.However, the optimized tool axis vector has visually visible smoothness and changes uniformly along the path, while avoiding interference during the machining process.
The average angular velocities between the initial and optimal tool axis vectors were calculated separately, and the results are shown in Figure 9.The variations in angular velocities are not entirely consistent on both surfaces of the blade, but the significant angle variations between the initial tool axis vectors have been corrected, i.e., the maximum value of the initial tool orientation on the basin surface is 0.5302, and the final is 0.4646.Moreover, the amplitude of variation in average angular velocity after optimization is relatively reduced, and the variations of adjacent tool axes are decreased.It can be seen that the algorithm given obtains continuously changing tool axis vectors with relatively uniform amplitudes of variation, and the planning twin-tool orientation is nearly smoothness.
value of the initial tool orientation on the basin surface is 0.5302, and the final is 0.4646.Moreover, the amplitude of variation in average angular velocity after optimization is relatively reduced, and the variations of adjacent tool axes are decreased.It can be seen that the algorithm given obtains continuously changing tool axis vectors with relatively uniform amplitudes of variation, and the planning twin-tool orientation is nearly smoothness.
(a) Angular velocity on dorsal surface of blade (b) Angular velocity on basin surface of blade

Pinch Milling Experimental
The nine-axis machine tool for pinch milling was used to validate the proposed tool orientation optimization algorithm.The two cutters are arranged symmetrically with respect to the workpiece rotation, which are both the same torus-shaped milling cutters.Any cutter has an independent rotational speed, which is set to 2200.The material of the blade is 1Cr12Ni2W1Mo1V, which is a difficult-to-machine material.The material of the two cutters is hard alloy.Using the planned optimal tool axis vector, the pinch milling process is implemented for the blade, as shown in Figure 10.It is shown that the twin-tool orientation continuously and uniformly varies along the tool path, the milling process is smoothest, the cutting speed is increased by 46% for the initial process, the noise during the cutting process is noticeably reduced, and cutting stability is achieved along the entire path.The machined surface is smooth and flawless, with no visible interference marks, which is feels delicate and smooth to the touch with low roughness.Meanwhile, the product's size and shape closely match the specialized inspection mold, and the machining accuracy fully meets the preset requirements.Therefore, the algorithm is very beneficial for ensuring good machining quality of the blade surface.Additionally, it can better leverage the performance of the machine tool to enhance processing efficiency.

Pinch Milling Experimental
The nine-axis machine tool for pinch milling was used to validate the proposed tool orientation optimization algorithm.The two cutters are arranged symmetrically with respect to the workpiece rotation, which are both the same torus-shaped milling cutters.Any cutter has an independent rotational speed, which is set to 2200.The material of the blade is 1Cr12Ni2W1Mo1V, which is a difficult-to-machine material.The material of the two cutters is hard alloy.Using the planned optimal tool axis vector, the pinch milling process is implemented for the blade, as shown in Figure 10.It is shown that the twin-tool orientation continuously and uniformly varies along the tool path, the milling process is smoothest, the cutting speed is increased by 46% for the initial process, the noise during the cutting process is noticeably reduced, and cutting stability is achieved along the entire path.The machined surface is smooth and flawless, with no visible interference marks, which is feels delicate and smooth to the touch with low roughness.Meanwhile, the product's size and shape closely match the specialized inspection mold, and the machining accuracy fully meets the preset requirements.Therefore, the algorithm is very beneficial for ensuring good machining quality of the blade surface.Additionally, it can better leverage the performance of the machine tool to enhance processing efficiency.
value of the initial tool orientation on the basin surface is 0.5302, and the final is 0.4646.Moreover, the amplitude of variation in average angular velocity after optimization is relatively reduced, and the variations of adjacent tool axes are decreased.It can be seen that the algorithm given obtains continuously changing tool axis vectors with relatively uniform amplitudes of variation, and the planning twin-tool orientation is nearly smoothness.
(a) Angular velocity on dorsal surface of blade (b) Angular velocity on basin surface of blade

Pinch Milling Experimental
The nine-axis machine tool for pinch milling was used to validate the proposed tool orientation optimization algorithm.The two cutters are arranged symmetrically with respect to the workpiece rotation, which are both the same torus-shaped milling cutters.Any cutter has an independent rotational speed, which is set to 2200.The material of the blade is 1Cr12Ni2W1Mo1V, which is a difficult-to-machine material.The material of the two cutters is hard alloy.Using the planned optimal tool axis vector, the pinch milling process is implemented for the blade, as shown in Figure 10.It is shown that the twin-tool orientation continuously and uniformly varies along the tool path, the milling process is smoothest, the cutting speed is increased by 46% for the initial process, the noise during the cutting process is noticeably reduced, and cutting stability is achieved along the entire path.The machined surface is smooth and flawless, with no visible interference marks, which is feels delicate and smooth to the touch with low roughness.Meanwhile, the product's size and shape closely match the specialized inspection mold, and the machining accuracy fully meets the preset requirements.Therefore, the algorithm is very beneficial for ensuring good machining quality of the blade surface.Additionally, it can better leverage the performance of the machine tool to enhance processing efficiency.

Conclusions
In this research, the twin-tool orientation smoothness planning method is proposed for thin-wall blade with pinch milling.Based on the two cutters structural layout, we clarify the coupling relationship of twin-tool posture, and a calculation method of the feasible region is defined considering multiple constraints.Then, the tool posture curve is constructed for tool orientation, and the global and local optimal planning algorithms are studied within the feasible domain.Furthermore, taking the tool axis vector on the tool posture curve as the optimal value and aiming for the minimum average angular velocity change, the twin-tool synchronization smoothing planning method is proposed for pinch milling.
The optimal twin-tool axis vector are successfully planned for a typical turbine blade using the proposed method, and are nearly smooth, and the collision-free requirement is guaranteed.By simulation and verification experiments, it is shown that the tool axis vectors of two cutters change uniformly along the path, and the pinch milling processing is very smooth.Therefore, efficient and high-precision processing of pinch milling can be achieved, ultimately improving the machining accuracy of the blade with difficult-tomachine materials.

Figure 1 .
Figure 1.The method of pinch milling.

Figure 2 .
Figure 2. Definition of tool orientation angles.

Figure 1 .
Figure 1.The method of pinch milling.

Figure 1 .
Figure 1.The method of pinch milling.

Figure 2 .
Figure 2. Definition of tool orientation angles.

Figure 2 .
Figure 2. Definition of tool orientation angles.

Figure 3 .
Figure 3. Twin-tool orientations constraint relationship of the CC point.

Figure 3 .
Figure 3. Twin-tool orientations constraint relationship of the CC point.

Figure 5 .
Figure 5. Feasible region in the workpiece coordinate system.

Figure 7 .
Figure 7. Flowchart of twin-tool orientation for pinch milling.
(a) the tool orientation on dorsal surface of blade (b) the tool orientation on basin surface of blade

Figure 8 .
Figure 8.The final twin-tool orientation for pinch milling.

Figure 8 .
Figure 8.The final twin-tool orientation for pinch milling.

Figure 9 .
Figure 9.Comparison of angular velocity by initial and finial tool orientation.

Figure 9 .
Figure 9.Comparison of angular velocity by initial and finial tool orientation.

Figure 9 .
Figure 9.Comparison of angular velocity by initial and finial tool orientation.