This paper presents a new vector-field-based streamline smoothing method in the parametric space and a tool orientation optimization technique for five-axis machining of complex compound surfaces with torus-end cutters. Iso-planar tool path is widely used in the machining of various types of surfaces, especially for the compound surface with multiple patches, but the operations of intersecting the compound surface with a series of planes have depended considerably on the complicated optimization methods. Instead of intersecting the surface directly with planes, a novel and effective tool path smoothing method is presented, based on the iso-planar feed vector fields, for five-axis milling of a compound surface with torus-end cutters. The iso-planar feed vector field in the parametric domain is first constructed in the form of stream function that is used to generate the candidate streamlines for tool path generation. Then, a G^{1} blending algorithm is proposed to blend the vector fields within the adjacent parametric domains to ensure smooth transition of cross-border streamlines. Based on the smoothened streamlines in the parametric domains, pathlines along with their correspondent side sizes are selected as desirable tool paths. Concerning a high performance machining, detailed computational techniques to determine the tool axis orientation are also presented to ensure, at each cutter contact (CC) point, the torus-end cutter touches the part surface closely without gouging. Both the computational results and machined examples are demonstrated for verification and validation of the proposed methods.

## Introduction

Due to the enormous advantages in describing complex shapes, compound surfaces, which consist of multiple surface patches, have been commonly found in industrial products including consumer products, die/mold, aerospace and automotive parts, etc. Such parts are often machined on five-axis machine tools that bring some obvious benefits such as better surface finish and higher machining efficiency compared with three-axis machine tools [1,2]. When machining the surfaces with multiple patches as a whole, the iso-planar tool path method is the typical choice to eliminate the machining marks caused by cutter's entering and leaving any two adjacent patches. It is usually done by slicing the surface using a series of parallel planes. The principle of generating iso-planar tool path is simple, but the practical operations of intersecting the compound surface with a set of parallel planes are still dependent on considerably complicated optimization procedures, such as iterative intersecting line tracing [3]. On top of that, selection of correct path interval between adjacent tool paths is also challenging. Engineers typically end up selecting conservative path interval step sizes, which leads into unnecessarily long machining time [4]. Due to the challenge of computation efforts, most methods are based on three-axis machining using a simple ball-end mill cutter. However, for the situation of five-axis machining with torus-end cutters, once the iso-planar cutter contact (CC) paths are determined, sequential adjustments of tool orientations are performed to avoid local and global gouging. Such practice suffers from the lack of further consideration of resultant changes in the cutter's machining strip width (MSW) along the whole tool paths [5].

To address these problems, a novel feed vector field tool path generation method is proposed in this paper. The proposed method takes advantage of the iso-planar feed vector fields in the parametric domains for continuous and smooth tool motions and, at the same time avoids the complicated intersecting computations. Change of the MSW caused by tool orientation adjusting is also used as a feedback in the tool path generation to ensure tool path intervals between adjacent paths are optimized. In Sec. 2, literature review is provided on the related works of tool path generation and tool orientation selection for five-axis milling of compound surfaces.

## Literature Review

### Tool Path Generation for Compound Surface.

The current cutting strategies of a compound surface can be generally classified into three categories. The first feasible way is to approximate the compound surface by a set of triangular meshes, and then, the tool path generation methods, such as iso-planar methods [6,7] or the conformal mapping-based methods [8–11] are used to generate tool paths for milling the part surfaces. Although the tessellation of the compound surface simplifies considerably the process of tool path generation, it suffers from loss of accuracy and the degeneration of machining quality. Thus, in practical industrial applications, a compound surface can also be machined patch-by-patch [12,13]. Using this method, most of the existing tool path patterns, such as the iso-parameter tool path [14], the iso-scallop tool path [15,16], the spiral tool path [17], and the tool path along the directions of various preferred feed directions [18–23], can in principle be applied to the machining of each individual patch of the given compound surface. Although it is beneficial for maintaining the machining accuracy of each patch, a difficult issue inherently in patch-by-patch machining is to find a good solution to the linkage of tool paths among different patches. Otherwise, unappealing cutting marks may be left in the neighborhood region of shared patch boundary.

To improve the machining quality in patch boundary, Cho et al. [24] and Chen and Fu [25] proposed methods to generate cross-boundary tool path with *G*^{0} position continuity based on the topological relationships among patches, so that cutter's entry/exit marks were avoided when crossing the common boundary. However, if the connection of tool paths between adjacent patches is only *G*^{0} continuity, there are possibilities of unfavorable frequent acceleration and deceleration when the cutter is moving across the patch boundary, running a risk of worse dynamic performance [26] as well as higher energy consumption [27]. To overcome the previously mentioned difficulties, in our previous work presented in Ref. [28], a vector-field-based method was proposed to generate tool path with *G*^{1} continuity between adjacent surface patches, by which the good kinematic performance and better surface quality were achieved. The result presented in our earlier work in Ref. [28] indicates the proposed method works well for compound surface with *C*^{1} continuity.

Other than the above two kinds of machining strategies on handling compound surfaces, the compound surface can be also machined as a whole with the well-known iso-planar method [29,30]. However, current iso-planar methods are mainly focused on three-axis machining. Also, as aforementioned, the current practice of finding intersecting curves between the compound surface and the sequence of intersecting planes is pretty computation intensive since it involves solving the nonlinear optimization to determine the desired CC points. Despite of this, the iso-planar tool path methods can generate continuous paths across all the adjacent patches with at least *G*^{0} continuity. When processing the compound surface consisting of multiple surface patches, the iso-planar tool path methods do have obvious advantages in generating continuous tool paths across all surface patches without interruptions of cutter lifting during machining. To maintain this advantage and simultaneously avoid the surface/plane intersecting operations, a novel streamline-type tool path generation method is proposed in this paper. The proposed vector-based tool generation method can generate continuous tool paths with *G*^{1} continuity across the boundary of two adjacent patches based on the iso-planar feed vector field as discussed in Sec. 3.

### Tool Orientation Determination.

Tool orientation is usually determined so that local and global gouging can be simultaneously avoided. In this paper, global gouging refers to the rear gouging that is beyond the neighborhood region at, or next to the CC point. On the other hand, local gouging is focused on the vicinity of the CC point, ensuring no occurrence of overcut. A common method used in local gouging avoidance is the Sturz method [31], in which a large inclination angle is set artificially to avoid the local gouging and kept constant throughout the machining process. After that, some adaptive tool orientation methods were proposed in Refs. [32–34] by matching normal curvatures of the cutter's effective cutting profile with the part surface at a CC point. Further, by means of the Dupin indicatrices of the cutter's swept surface and the part surface, several new methods were developed to find the locally optimal tool orientation [35,36]. Lately, to increase the machining strip width, second-order approximation and third-order approximation methods were proposed in Refs. [37] and [38] for local tool orientation optimization. Theoretically, these local methods can produce a maximum nominal MSW without local gouging. However, since these methods ignore the possible global gouging, they might still suffer the risk of overcuts in nearby regions not immediately at, or next to the CC point.

To achieve the optimal tool orientation without global gouging, the multipoints machining method [39], rolling ball method [40], arc-intersection method [41], and penetration elimination method [42] have been developed in the last few years. These methods are focused on finding the smallest gouge-free tool inclination angle under the constraint of zero tilt angle. These results are satisfactory for global (rear) gouging avoidance but are still far from the optimal solutions in terms of the maximum MSW [43,44]. In Ref. [45], the part surface was approximated locally by a generic quadric approximation at the CC point, and then, the tool orientation angle was optimized to enlarge the MSW and to eliminate machining gouging at the same time. However, due to the inconsistency of the fitted quadric and the original surface, the results might be inaccurate in some extreme cases and still need further improvement. In Ref. [46], the tool orientations, which tend to maintain the hyper-osculating configurations when it is accessible, and smoothly switch to nearly two-contact configurations if it is not accessible due to some constraints, are generated by a global optimization procedure for five-axis machining of freeform surfaces using a flat end cylindrical tool. As a result, a good approximation quality of machining can be obtained at a relatively high computation cost. Recently, an improved rotary contact method (IRCM) [47] was proposed to adjust the tool axis by using a nested optimization of the two tool orientation angles until the cutter touches the part surface at two points. Although the method is simple and interesting, the apparent increase in computation complexity and demand of heavy computation efforts will become burdensome when trying to optimize the two orientation angles of the cutter at the same time.

In this paper, a new global tool positioning method is proposed by combining the advantages of curvature matching. Namely, the tool is initially positioned by using a local gouging avoidance method and then follows with a global gouge checking. Once the undesirable overcut occurs, the tilt angle is first optimized to decrease the penetration error. After that, adjustment of the inclination angle is performed for global gouging avoidance. This method is generally simple in implementation, and it can guarantee large MSW like other multipoint contact methods.

The remainder of this paper is arranged as follows: the vector fields fitted to the iso-planar feed directions, streamline formulation, and *G*^{1} blending of the vector fields are first presented in Sec. 3. The procedures of finding the optimal tool orientation at each CC point are discussed in detail in Sec. 4. Then, full length tool paths on the whole compound surface are generated in Sec. 5. The analytical results and experimental validations are presented in Sec. 6, followed by the concluding remarks in Sec. 7.

## Vector-Fields-Based Streamline Construction

### Vector Field Fitted to Iso-Planar Feed Directions in Parametric Domain.

*x*th patch of the compound surface has been uniformly discretized into a grid of points, and the total number of sample points is

*s*+ 1. The iso-planar feed direction $\delta ui,\delta vi,i=0,1\cdots s$ in the parametric plane satisfies

*i*th sample point, and $np$ is the unit normal of the intersecting planes. The iso-planar feed direction $\delta ui,\delta vi$ in the parametric domain can be obtained by solving Eq. (1) as follows:

where $Ei=Sui\xb7Sui,Fi=Sui\xb7Svi,andGi=Svi\xb7Svi$. For the convenience of vector field construction, the calculated vector $(\delta ui,\delta vi)$ by Eq. (2) needs to be normalized.

*k*-degree B-spline basis function defined over $U$, and $Nj,lv$ is a

*l*-degree B-spline basis function defined over $V$. $U$ and $V$ are both clamped and open knot vectors. $\Omega i,jx$ is the element of the control matrix $\Omega x$. Similar to the vector interpolation method presented in Ref. [48], in this paper, the streamlines are defined and derived as the level curves of a stream function. Taking the feed direction as the tangential direction of a streamline, the streamline function can be easily derived based on the fact that the gradient of a point on the stream function is orthogonal to its tangent on the level curve, as cast in the form as follows:

Using Eq. (5), tool paths on the *x*th patch can be generated by mapping the streamlines from the parametric domain to the physical space.

*G*^{1} Blending of Adjacent Vector Fields.

*k*th degree with the shared border at

*u*= 1 (

*w*= 0), and the corresponding knot vectors are given by

**,**

*U***and**

*V***,**

*W***. At a point on the shared border, the condition of**

*V**G*

^{1}blending streamline is given as follows:

**T**

^{1}and

**T**

^{2}are tangents of streamlines for $\phi 1(u,v)$ and $\phi 2(w,v)$, respectively, and $\lambda v$ is a ratio coefficient. For the convenience of processing, $\lambda v$ is assumed to be a constant $\lambda 0$. By simplifying Eq. (6) with the boundary condition, one can obtain the following equations:

where $w\xaf1,w\xaf1+k,u\xafm$, and $u\xafm+k$ are elements in the knot vectors ** W** and

**.**

*U**a*) and (7

*b*) valid, the original control coefficients of the control matrix $\Omega 1$ and $\Omega 2$ need to be adjusted under the principle of minimum deviation. That is, the deviation between the adjusted coefficients and the original coefficients has to be minimized. In the earlier method presented in Ref. [28], the relation expressed in Eq. (7

*a*) is processed by fixing the last row of $\Omega 1$ and adjusting only the first row of $\Omega 2$. In this paper, to handle the associated rows simultaneously, a new and more effective blending method is proposed. Let $\Omega \u0303i,j$ be the adjusted control coefficient of $\Omega i,j$. The deviation $\epsilon 1$ can be expressed as follows:

*a*) are determined, the control coefficients involved in the

*m*− 1th row of $\Omega 1$, and the second row of $\Omega 2$ can be further adjusted by using Eq. (7

*b*). The adjustment deviation error $\epsilon 2$ can be defined as follows:

where $a=w\xaf1+k\u2212w\xaf1/u\xafm+k\u2212u\xafm.$ Update the original control coefficients with the adjusted ones, two vector fields can be smoothly blended so that the smoothened streamlines with *G*^{1} continuity can be achieved at the shared border of adjacent parameter domains, as shown in Fig. 2(b).

## Determination of Adaptive Tool Orientation

In this paper, the tool orientation selection method starts first with a local method by matching the normal curvature of the cutter swept surface with that of a compound surface in all directions on the tangent plane, which can be also called “comprehensive curvature matching” [35]. After that, the overcut error of the tool positioning is evaluated to see whether the global gouging appears. If the cutter does penetrate the part surface, both the two rotational axis angles are adjusted to avoid the global gouging; otherwise, the current local tool orientation is selected as the optimal one. Details of the proposed tool orientation selection method are discussed in Secs. 4.1 and 4.2.

### Tool Orientation Via Curvature Matching.

**. The radii of**

*T***and the torus are**

*T**R*and

*r*, respectively. A local coordinate system is established at the CC point, in which the tool orientation is defined by the orientation angles $(\lambda ,\omega )$. $xL\u2009\u2009and$$zL$ axes of local coordinate system are defined as the unit tangent of the tool path and the surface normal

**at a CC point**

*n***; $yL$ is the cross product of $zL$ and $xL$. $\lambda $($0\u2264\lambda \u2264\pi /2)$ is the inclination angle measured counterclockwise around $yL$, and $\omega (\u2212\pi /2\u2264\omega \u2264\pi /2)$ is the tilt angle measured counterclockwise around $zL$. To determine the geometric relationships between the cutter and the part surface, two additional coordinate systems are required, which are the tool coordinate system (TCS) and the workpiece coordinate system (WCS), respectively. By defining $xT$ as the unit vector from the origin of TCS to the CC point, $zT$ as the unit tool axis vector and $yT=zT\xd7xT$, the swept surface $Et,\theta $ of cutting torus can be mathematically expressed in the WCS by the kinematic transformation as follows:**

*P*where $\kappa maxS$ and $\kappa minS$ are the principle curvatures of the part surface, $\beta $ specifies the maximum principle direction $fmax$ of the surface at a CC point, and $\tau $ is the surface torsion in terms of $\beta $.

*I*and

_{E}*II*are the first and second fundamental forms of $Et,\theta $. Detailed formulas of finding the first and second fundamental forms, namely,

_{E}*I*and

_{E}*II*are provided in the Appendix. Expanding Eq. (19) yields an explicit expression on $\kappa \alpha E$

_{E}From Eq. (22), it can be seen that a locally optimized tool orientation can be achieved by curvature matching in only one direction on the tangent plane, i.e., $\alpha =\pi /2$. However, due to complexity of the surface shape, the gained tool orientation via the curvature matching method alone cannot guarantee that it is gouge-free between the cutter and the neighborhood part surface except for the vicinity of CC point. Further, global gouging avoidance procedures are required as discussed in Sec. 4.2.

### Global Gouging Detection and Avoidance.

**in the WCS can be expressed as**

*T*The values of $ei|i=1,2,\u2026n$ reveal the geometrical relation clearly: a negative value implies the existence of cutter gouging, and a positive value indicates residual cusp material left on part surface after machining and zero indicates $Qi$ is a contact point. Therefore, the gouging status between the cutter and the surface can be quantified with the negative errors. The minimum negative value $\Delta =min(ei|ei\u22640)$ is called as the penetrating error (PE), which is a dominated factor in the adjustment of tool orientation.

**as $el$ and $er$, respectively. Generally, $el$ and $er$ are different under the initial tool orientation due to likely geometric asymmetry at the CC point as shown in Fig. 6. If $er>el$, the tilt angle should be adjusted anticlockwise, otherwise the cutter should be tilted clockwise. In the process of tilt angle rotating, the larger error will decrease while the smaller error will increase. The optimal tilt angle will lead to equal penetrating errors on both sides of the CC point, that is**

*T*Then, in this case, a *bisection searching routine* is used to find the optimal tilt angle $\omega opt$ with $\lambda 0$ fixed. To accelerate the searching speed, $el$ and $er$ are calculated quickly based on the shortest distances from the left and the right half of ** T** to the part surface and the shortest distance between two objects can be determined efficiently with Newton iteration method. From the error distribution curves shown in Fig. 7, one can see that the smallest PE increases dramatically after optimizing the tilt angle. In addition, the approximately symmetry error distribution (Fig. 7) can be also obtained, which is beneficial to superior machining quality.

where $GP$ denotes the point on cutter surface corresponding to the penetrating error $\Delta $.

In order to completely eliminate the possible gouging, the final inclination angle is required to make sure $g\u22650$. In other words, no points on the cutter will penetrate the part surface. To achieve this goal, we first retrieve the two points corresponding to $el$ and $er$ on the cutter. Denote the two points by GP1 and GP2, respectively. Then, GP1/GP2 will detach from the part surface as $\lambda $ increases. The optimal inclination angle $\lambda 1$ for GP1 arises when GP1 contacts the part surface. $\lambda 1$ can be computed with a bisection method, and the optimal inclination angle $\lambda 2$ for GP2 can be achieved with the same process. The adjusted inclination angle $\lambda opt$ is selected as the larger one from the calculated inclination angles ${\lambda 1,\lambda 2}$. Optimization of $\lambda $ will terminate if no gouging points are detected in the subsequent gouging check process. Otherwise, the above mentioned optimization algorithm on $\lambda $ will be performed iteratively until the obtained tool orientation is gouge-free.

The presented tool orientation method can guarantee at least two-point contact configuration and is similar in spirit to the multipoint contact method. Since the initial condition is properly adopted, difference between $\lambda opt$ and $\lambda 0$ is slight, so that the error distribution in $(\lambda opt,\omega opt)$ will be approximate to the desirable one in $\lambda 0,\omega opt,$ which denotes better surface quality after machining.

When machining the shared patch boundary, the locally optimal tool orientation may be influenced by geometric features of adjacent patches. For example, if a CC point is located at the common boundary, the principle curvatures at this point are different for the two adjacent patches. In this situation, a constant tilt angle is suggested, and the larger inclination angle computed by Eq. (22) is selected as the local tool orientation with respect to these two adjacent patches. After that, procedures of gouge checking and avoidance will be implemented to search the globally optimal tool orientation

## Tool Path Generation for Machining the Compound Surface

Streamlines in the vector field can be generated by solving Eq. (5) with a fourth-order Runge–Kutta solver. Once the initial point/seed point (SP) is determined, the Runge–Kutta integration can be implemented with a forward or backward step to trace a streamline. As long as a set of streamlines is generated with reasonably controlled intervals in the parameter domain, their corresponding pathlines on the part surface can be used as desired tool paths.

*h*, two points with the position error equal to

*h*can be found on both sides of a CC point under the optimal tool orientation. Denote their corresponding foot points on the surface with $Pl$ and $Pr$, the machining strip width

*L*is computed by the projection distance between $Pl$ and $Pr$ along the unit vector $yL$ as follows:

Denote surface parameters of $Pl$ and $Pr$ by $(ui,jl,vi,jl)$ and $(ui,jr,vi,jr)$, respectively, subscripts *i* and *j* denote the *j*th point on the *i*th streamline, *l* and *r* reveal their locations in terms of the feed direction at a CC point. To maximize the machining efficiency, the interval between two adjacent streamline paths should be as large as possible, meanwhile satisfying the machining tolerance. This is an iterative process as seen in Fig. 8. Assume that the streamlines are arranged in order from bottom to top in the parametric domain. The machined boundaries $Bil$, $Bir$ of the *i*th streamline and $Bi+1l$, $Bi+1r$ of the *i +* 1th streamline are formed by interpolating the points $(ui,jl,vi,jl)$, $(ui,jr,vi,jr)$, $(ui+1,jl,vi+1,jl)$, and $(ui+1,jr,vi+1,jr)$, respectively. To achieve the largest possible interval between the *i*th and *i +* 1th streamline, $Bil$ should lie on the upper side of $Bi+1r$ and be as close as possible to $Bi+1r$. This can be realized by adjusting the SP of the *i* + 1th streamline through trial and error. For a single patch, a set of feasible streamlines can be generated in sequence with this process until the whole parametric domain is covered.

For the streamline crossing the common boundary of adjacent patches, both the adjacent parametric domains and their associated stream functions are simultaneously taken into consideration in tool path generation. To search the SP satisfying the machining requirement of all domains crossed by the streamline, the SP for every parametric domain is first searched using the method discussed earlier; then, the intersection of a streamline with the common boundary is easily obtained based on the features of a streamline. Two intersections can be found for the boundary between two adjacent patches. The one that satisfies machining requirements of both patches is selected to be a potential SP. The final SP of the next cross-border streamline can be determined by recursively invoking the same process for every two adjacent patches until all the related domains are traversed.

**and the tool axis vector**

*O***, can be obtained with the CC point**

*I***and the assigned orientation angles ($\lambda ,\omega )$ as expressed in the following equation:**

*P*## Examples

The developed tool path generation method has been coded in C++ language and implemented on a 2.8 GHz personal computer. In the first test, a compound surface with two patches as illustrated in Fig. 9(a) is used to demonstrate the validity and efficiency of the proposed method via simulation. The compound surface is modeled in Unigraphics 9.0 by stitching two Bézier patches with *C*^{1} continuity. A torus-end cutter with *R* = 6 mm and *r* = 2 mm is adopted to machine the surface under the scallop height *h* = 0.05 mm. Feed directions of the sample points are chosen to be parallel to the *X*–*Z* plane. Then, the iso-planar feed vector fields are constructed and blended by using the method of Sec. 3. On basis of this, the streamline tool paths are generated by the proposed method and shown in Fig. 10(a). By comparison, the traditional iso-planar tool paths generated by slicing the compound surface using a series of *X*–*Z* planes are shown in Fig. 10(b). It can be seen that our tool paths are similar with iso-planar tool paths in topology. The computation times listed in Table 1 show that the proposed method is more efficient than the iso-planar method, and the computation efficiency is almost doubled.

For our proposed tool path and the iso-planar tool path, the tool orientation at a CC point is determined by the proposed tool orientation method presented in Sec. 4. Then, the milling simulation of these two kinds of tool paths is performed in Unigraphics 9.0. The machined tolerance maps of these two tool paths are shown in Figs. 11(a) and 11(b), respectively. One can see that there are no overcuts on the surface. Although the machining error exceeds the required scallop height in some fractional regions, it is small enough to be neglected in the practical machining. Table 2 lists the comparison of the proposed tool orientation method with the IRCM method at ten uniformly selected points. In this test, the offset value of IRCM is set to 0, and the initial tool orientation is selected as (0, 0). Considering that it is numerically intensive when optimizing the tilt angle with a constant step as in the IRCM, a bisection method is also used for the tilt angle optimization process just as ours. Both the methods produce almost the same optimal tool orientations. However, compared with our method, the IRCM is still lower in computation efficiency due to more iterative optimizations of the two orientation angles. The experimental results, listed in Table 2, demonstrate this point. This shows that our proposed method can not only generate gouging-free tool orientation but also have higher computation efficiency.

The second compound surface, consisting of four B-spline patches as shown in Fig. 9(b), is retrieved from the computer-aided design model of a mold. On this surface, for the need of test, the traditional iso-planar tool path method is performed patch-by-patch, and the intersecting planes are selected to be parallel to the diagonal line of each patch. The generated iso-planar tool paths are shown in Fig. 12(b). Through the sampled iso-planar feed directions, the proposed method is performed, and the generated vector-field-based tool paths are shown in Fig. 12(a). Time consumption of these two kinds of tool paths is also provided in Table 1, and similar conclusion can be drawn as in Test 1. Comparing the generated tool paths, as expected, it is observed from the magnified view shown in Fig. 12(c) that the discontinuity occurs at the shared boundary for the iso-planar tool paths, while the tool paths generated by the proposed method are very smooth and continuous. An actual machining experiment is also conducted on the proposed tool path pattern and the benchmarking one. The scallop height is set to be 0.05 mm, and a torus-end cutter of radius with *R* = 4 mm and *r* = 2 mm is utilized to mill the part. The machined results are shown in Fig. 13. It reveals that the surface machined by the proposed method has better surface finish than that of the iso-planar method at the shared borders. The machining precision of the surface profile is measured by a Prismo coordinate measuring machine (CMM) in our precision machining lab. Three cross sections at *X* = 55 mm, *Y* = 45 mm, and *X* = 10 mm (*X* = 55 mm and *Y* = 45 mm are shared patch boundaries) are measured, respectively, on the parts shown in Figs. 13(a) and 13(b). Figure 14 shows the differences between the theoretical profile and the measured profiles. In Fig. 14, CMM data1 denote the machined profile by the proposed method, and CMM data2 are for the iso-planar method. From the results, it is seen that the maximum undercut has been constrained below the specified machining tolerance and no overcut/gouging appears. Hence, the machined parts satisfy machining requirements and validate the effectiveness of the proposed method.

## Conclusions and Future Works

In this paper, the machining of compound surface has been realized with the proposed vector-field-based smoothened streamline tool path method. This presented method takes advantages of the iso-planar feed vector and avoids the complicated intersection operations of the model surface with the intersecting planes. Since it realizes the *G*^{1} blending of feed vector fields in the parametric domains of adjacent patches, the generated tool path can cross smoothly the shared border with *G*^{1} continuity compared with the traditional iso-planar tool path. The proposed tool path method is an attempt to schedule tool paths with *G*^{1} continuity for machining the whole compound surface. The experimental results have demonstrated its validity. In addition, an effective tool positioning method is also presented by eliminating penetration errors between the cutter and the part surface. The proposed method is easy to implement, and the computed tool orientations are advantageous in achieving reasonably good machined surface quality as demonstrated in both the simulation and the actual machining experiments. Compared with IRCM method, the computation time for the optimal tool orientation is also greatly reduced.

In the future work, the emphasis should be focused on the theoretical analysis of continuous tool path blending for the arbitrary-bounded adjacent patches. We are also interested in combining the specific machine tools' kinematic constraints and the tool orientation optimization for better machining performance.

## Funding Data

National Natural Science Foundation of China (51525501, 11290143, and 51621064).

SCP to Dalian University of Technology (JCKY2016212A506-0102).

National Science Foundation Grants to North Carolina State University (CMMI-1547105 and CMMI-1404916).