Flat-end tool orientation based on rotation-minimizing frame

Well-designed tool orientation is crucial for the quality execution of five-axis machining with flat-end tools. Tool orientation can be specified in terms of tool inclination and tool tilt angles. For a given surface path, these aspects need to be specified for each cutter contact point along the path, because poor tool orientation choice can cause large axial acceleration of the machine tool, leading to inferior quality of the machined surface. The rotation-minimizing frame aims to reduce unnecessary frame rotation during movement. This concept has been adopted to develop tool orientation methods that allow variation in the inclination angle, the tilt angle, and both the inclination and tilt angles. The intention is to reduce unnecessary rotation of the tool frame as the tool follows a specified path. Evaluation was conducted in the context of five-axis flat-end tool machining. Based on these methods, tool orientation was planned along surface paths of a torus, sphere, and dome. Changes in tool orientation were always smooth. From the perspective of reducing tool orientation changes and axial acceleration, it was demonstrated that simultaneous variation of the inclination and tilt angles based on rotation minimization provided the best results.


Introduction
Tool path planning for five-axis computer numerical control (CNC) machining has attracted significant interest in recent years. Constraints to the planning process arise from gouging avoidance, kinematic limits, singularity avoidance, cutting precision, and efficiency requirements. However, complete solutions to path planning have yet to be identified.
An attempt was made to divide the part surface into convex and non-convex regions, and to determine the optimal cutter orientation for local and rear gouging avoidance, by considering the offset of the machined surface and the offset cylinder of the cutter [1]. Maximization of the machining quality was achieved by identifying gouging-free tool paths, analyzing hyper-osculating circles that have third-order contact with the target surface, and formulating two-contact configurations as systems of algebraic constraints [2]. The concept of forbidden circles was proposed to avoid singularities in five-axis machining [3]. The problem was considered two-dimensional, as a three-dimensional ''singular cone'' associated with tool orientation was projected onto a two-dimensional space. Efforts were made to compensate for geometric errors to improve machining precision based on a generalized inverse kinematic model of the five-axis machine tool [4,5]. According to variations in cutting force and stability considerations along the tool path, efficient surface milling parameters were identified for tool path generation, the aim being to shorten the milling time [6].
Typical path planning involves generating an initial tool path, including the identification of cutter contact (CC) points and cutter positions, and construction of adjacent tool paths and CC points with consideration of factors such as material removal rate, scallop height, gouging, and machine dynamics. It is crucial to identify the tool orientation at each of the CC points. The focus of this paper is to address tool orientation for five-axis machining with flatend tools.
Tool orientation is an important factor that affects the surface quality of a machined part [7]. Existing research in tool orientation for five-axis machining has focused mainly on gouging avoidance, rotary axes movement optimization, and cutting efficiency.
To avoid local gouging, rear gouging, and global tool collision, search methods based on configuration space (C space) [8][9][10] were proposed to determine the tool orientation, in which the motion plan of the object was transformed into the motion plan of a point in the C space. Ezair and Elber [11] further utilized conservative computation of accessibility for a finite surface and orientation elements to ensure the correctness of collision-free orientation.
For chatter-free tool orientation, structural dynamics of the tool and work piece were analyzed in terms of cutter workpiece engagement coordinates [12]. An iterative search for gouging-free tool axis orientation was performed using the Nyquist criterion.
To maximize material removal and avoid local gouging, tool orientation was planned along the tool path to maintain a constant scallop height [13,14]. Decisions were made based on the local surface curvatures and analyses of the effective cutting radius [15,16]. In actuality, there is a trade-off between increasing material removal rate and reducing scallop height between adjacent paths [17].
Li et al. [18] presented a strategy for tool orientation planning, which divided a complex surface into regions based on both the Gaussian and mean curvatures, to reduce the rate of change in the tool orientation. Castagnetti et al. [19] proposed the concept of domain-of-admissible-toolorientation to minimize the movement of each rotary axis. Ho et al. [20] adopted a circular arc representation of tool path segments to smoothen tool orientation changes between adjacent CC points. Jun et al. [21] formalized a weight criterion in the C space, which consisted of the square root of variation of the tool tilt and inclination angles between two consecutive CC points, for smoothing the tool orientation.
Hu and Tang [22] focused on improving the dynamic performance of the machine tool and reducing the maximum angular acceleration at the rotary axes. An optimal work piece setup was first identified through a heuristic search. Subsequently, the tilt and yaw angles were found to optimize the tool orientation. An accurate cutter/work piece engagement model was developed by Zhang et al. [23] for optimizing the tool inclination and tilt angles for efficient and accurate machining. Xu [24] proposed an interpolation scheme for online calculation of tool orientation at each interpolated point to improve machining efficiency. Li et al. [25] suggested a two-stage process to reduce the cutter posture change between two adjacent cutter contact points. This was achieved by recognizing the range of cutter posture angles for collision avoidance and estimating the direction for smoothest posture change. Lin et al. [26] developed a method to change the tool orientation to avoid five-axis singularities and improve the texture of the machined surface.
The above literature review shows the following: (i) tool orientation can be represented by the tilt angle and inclination angle, (ii) tool orientation algorithms mainly involve finding the best angles based on certain objectives while satisfying certain requirements, and (iii) these algorithms should aim at stable and non-excessive rotary axes movement.
Farouki and Li [27] and Farouki et al. [28] used a rotation-minimizing frame (RMF) to reduce variation in tool orientation in five-axis machining with ball-end tools. The tool inclination angle was fixed to maintain constant cutting speed, and the component of the tool vector on the tangent plane was specified as the parallel transport vector. This resulted in smooth overall tool orientation variation, and the total tool rotation about the surface normal was smaller than that in Sturz's method. However, Han [29] pointed out that tool axis movement was not always small for an RMF with a parallel transport vector. Alternate RMFs for various applications were reviewed by Farouki [30]. It appears that methods favoring a small and stable tool axis speed improve efficiency and ensure robust machine tool performance.
In general five-axis machining with flat-end tools, both tool inclination and tilt angles may be varied during machining. Previous work has demonstrated the benefits of using an RMF to promote smooth changes in the tool axis. However, the cited works only considered changing the tool inclination angle, while the tilt angle remained constant. For cutting with ball-end tools and corner radius tools, a constant tilt angle implies a constant cutting speed. For cutting with flat-end tools, it may be advantageous to vary the tilt angle as well. This paper intends to address the use of an RMF to adjust the tool inclination angle and/or the tool tilt angle. Basic concepts pertaining to RMFs are reviewed in the next section. Section 3 presents methods for tool orientation for RMFs. Comparison of these methods, which includes tool paths on a torus, sphere, and dome, is conducted in Section 4. This is followed by a simulation of tool orientation changes on a free-form surface. The final section concludes the paper.

A review of RMFs
Due to its superior rotation motion properties, RMFs have drawn much attention for diversified applications such as rigid body movement planning [31][32][33][34][35] and construction of RMFs on multi-curves and surfaces [36][37][38][39][40]. In this section, basic concepts related to RMFs are reviewed, and the elements of RMFs relevant to this study are introduced.
For an orthogonal frame Þmoving along a curve C n ð Þ in three-dimensional space, the movement is said to be rotation-minimizing with respect to a certain reference direction if the angular velocity of F has no component along that direction.
For convenience, when the reference direction is f 1 , it is named the rotary vector, and f 2 and f 3 are the basic vectors. The angular velocity x of F can then be expressed as where r ¼ C 0 n ð Þ j j is the rate of change of arc length associated with the curve.
The derivative of the basic vectors with respect to n become Basic vectors f 2 and f 3 do not rotate instantaneously about f 1 , and their derivatives must be parallel to f 1 as For an orthogonal frame moving along a three-dimensional space curve, with one of the frame vectors (i.e., f 1 ) aligned to certain reference direction, the frame is still free to rotate about that reference direction. According to the concept behind the RMF, f 2 and f 3 are selected so that they do not rotate instantaneously about f 1 . In other words, there is no ''unnecessary rotation'' of the frame about the rotary vector as the frame moves along the space curve.

RMF-based tool orientation
Þ on the surface. The surface is smooth and oS ou Â oS ov does not vanish on the surface. Let n be the unit surface normal, t the unit curve tangent that points to the feed direction, and k ¼ t Â n the unit bi-normal at the CC point P on the path. From computational geometry, where j g , j n , and s g are the geodesic curvature, normal curvature, and geodesic torsion, respectively. Vector T is the unit tool vector. Inclination angle a is the angle between n and T. Vector m is the unit vector aligned to the projection of T on the tangent plane at P. Tilt angle b is the angle between t and m. Vectors T and m can be written as In five-axis machining, tool orientation must be prescribed as it follows path C to machine surface S. Tool orientation can be treated as selecting a and b along the path. In this section, RMFs are constructed according to various settings of the inclination and tilt angles. These include the following: (i) varying the tilt angle and maintaining a fixed inclination angle, (ii) varying the inclination angle and maintaining a fixed tilt angle, and (iii) varying both the tilt and inclination angles.

Varying the tilt angle
In this setting, tilt angle b is allowed to vary while inclination angle a is fixed. This can be considered as a direct application of Farouki's method [27]. As shown in Fig. 1, when a is fixed, the component of T along the normal direction is also fixed. When b is allowed to vary, the component of T on the tangent plane may vary. RMF n; h; m ð Þis selected with n as the rotary vector and m and h as basic vectors, where h ¼ m Â n. This choice of the RMF reduces variation in the basic vectors as rotary vector n varies along the tool path. For instance, the basic vectors do not rotate about n and their instantaneous change is parallel to n. This rotation method is named RM I. From Eqs. (4) and (6), the derivative of m with respect to n along the curve can be written as As t and k are orthogonal to n, and m 0 Â n ¼ 0 according to Eq. (3), components of m 0 along the directions of t and k in Eq. (7) must vanish. Obviously, sin b and cos b in these components cannot vanish simultaneously. Hence, the remaining option for satisfying m 0 Â n ¼ 0 is and Eq. (7) can be simplified as where b 0 is an integration constant.
Computation of the tool orientation T requires selection of the fixed tool inclination angle a and the initial tilt angle b 0 . Then, T can be calculated using Eqs. (5) and (6), and b n ð Þ can be evaluated numerically using Eq. (10) as the tool moves along the tool path.
Furthermore, from Eq. (5), the derivative of T is which can be written as Using Eqs. (4) and (9) and rotation transformations, t ¼ À sin bh þ cos bm, k ¼ cos bh þ sin bm, and w ¼ h Â T ¼ À cos am þ sin an.

Varying the inclination angle
In this setting, inclination angle a is allowed to vary while tilt angle b is fixed. As shown in Fig. 1, h is normal to the m-n plane and T is fixed on that plane. If b is fixed, so is h. When a varies, the plane, together with T, rotates about h. From Eqs. (4) and (6) and using rotation transformation, the derivatives of n and m can respectively be written as It is noted that the term Àrj g h is present in Eq. (14) but not in Eq. (7). This is because b is fixed in RM II but is allowed to vary as b 0 ¼ rj g in RM I. From Eq. (5), the derivative of T can firstly be written as and then rearranged using Eqs. (13) and (14) as As m and n are orthogonal to h, and T 0 Â h ¼ 0 according to Eq. (3), the components of T 0 along directions of m and n in Eq. (16) must vanish. Obviously, sin a and cos a in these components cannot vanish simultaneously. Hence, the remaining option for satisfying Equation (16) can then be simplified as Integrating Eq. (17) leads to where a 0 is an integration constant. Computation of tool orientation T requires selection of fixed tool tilt angle b and initial inclination angle a 0 . Then, T can be calculated using Eqs. (5) and (6), and a n ð Þ can be evaluated numerically using Eq. (19) as the tool moves along the path.

Varying the tilt and inclination angles (dual RMF)
To adapt to the local surface geometry more successfully, it may be preferable for the tool to be able to rotate about the normal direction, and to vary its inclination as it goes to the next tool contact point along the path. Thus, in a more general setting, this implies that both the inclination angle a and tilt angle b are allowed to vary to allow more flexibility for orienting the tool. Change in the tool orientation involves a two-sequential rotation of the tool axis. An RMF is employed so that variation of the basic vectors is reduced in each rotation by suppressing the rotation of the basic vectors about the rotary vector. Thus, this sequential rotation method is named dual RMF (DRM), which combines the primary RMF (PRMF) with the secondary RMF (SRMF).
In PRMF, the tool is set to rotate about the normal axis. Thus, n is selected as the rotary vector, with h and m as the basic vectors. Tilt angle b is allowed to vary while inclination angle a is fixed. In other words, the PRMF is the same as RM I, and Eq. (8) is applicable to the determination of b.
In SRMF, the tool is set to rotate about h, one of the basic vectors in PRMF. Thus, h is selected as the rotary vector, with T and w as the basic vectors. Inclination angle a is allowed to vary while tilt angle b is fixed.
From Eq. (5), the derivative of T with respect to n can be written as Here, n does not depend on b but m does. The derivation of n 0 in Eq. (13) does not involve any assumption with regard to b, and it can be used to express n 0 in Eq. (20). Conversely, the derivation of m 0 in Eq. (14) involves assumptions with regard to b and it cannot be used. On the other hand, as the determination of b in PRMF and that in RM I are exactly the same, Eq. (9) can be used to express m 0 . Therefore, Eq. (20) may be modified to T 0 ¼ Àa 0 sin an þ cos a Á r j n sin b À s g cos b À Á h À À j n cos b þ s g sin b À Á m Á þ a 0 cos am þ sin a Á r j n cos b þ s g sin b À Á n ¼ a 0 À r j n cos b þ s g sin b À Á À Á cos am À sin an ð Þ þ r cos a j n sin b À s g cos b À Á h: As the components of T 0 along the directions of m and n in Eq. (21) must vanish, the same restriction on T 0 needs to be imposed here (as in Eq. (17)), and T 0 can then be further simplified as Computation of tool orientation T requires selection of initial inclination angle a 0 and initial tilt angle b 0 . Then, a n ð Þ and b n ð Þ can be evaluated numerically using Eqs. (19) and (10) respectively along the tool path, and T can be calculated using Eqs. (5) and (6) (see Fig. 2).

Kinematics
Tool contact location and tool orientation information (together with tool and cutting parameters) are translated to axis command for the execution of five-axis machining [26]. The translation is conducted according to the kinematic chain of a machine. There are several common kinematic chains in five-axis machine tools [23][24][25]. For an A-B type five-axis machine (i.e., an orientable-spindle machine), rotary axes A and B are parallel to the translation axes X and Y.  Such a machine model (see Fig. 3) is adopted in this paper for kinematics calculations. The tool axis is along the Z-axis. x t ; y t ; z t ½ , x m ; y m ; z m ½ and x p ; y p ; z p Â Ã represent the tool coordinate system, machine coordinate system (MCS), and part coordinate system (PCS), respectively. Rotary axes A and B are parallel to translation axes x t and y t . The kinematic equation can be established in MCS as where T rans represents the transport matrix and R ot is the rotation matrix; x 0 ; y 0 ; z 0 ð Þis the initial tool contact location, and L T ½ denotes the tool contact location and orientation in MCS. From Eq. (23), the relationship between the tool axis vector and rotary angles can be obtained as Then, the rotary angles can be solved as In practical machining, X, Y, Z, A, and B are generated after tool path generation and tool orientation determination. The focus of the current work is tool orientation determination with an RMF. After the tool path on the part surface has been generated, the tool orientation is determined (tilt and inclination angles). Then, rotary angles A and B can be calculated using Eq. (25). It should be mentioned that a singularity may arise during the calculation of B. Singularity avoidance is an important issue, but is not considered in the current work.

Constraints on inclination and tilt angles
For five-axis machine with a flat-end tool (see Fig. 2), the effective cutter radius r eff is expressed as where r is the tool radius, a ¼ r sin a cos b, and h ¼ arctan tan a tan b ð Þ . The effective cutter radius must not be larger than the radius of curvature at a CC point. In general, a should not be large enough to promote material removal. Figure 4 shows an undesirable situation of rear gouging in which the rear of the tool interferes with the part surface.
This can happen when a is too small as the tool passes a concave part of a surface. Thus, to avoid rear gouging, a should not be too small. Further, b should be small so that h can be small, allowing r eff to be large.
These considerations should be accounted for in the selection of the initial inclination and tilt angles. Other issues, such as collision avoidance and material removal efficiency, are also important for tool path and orientation planning; however, these aspects are beyond the scope of this paper.

Evaluation
The focus of the evaluation is to reduce the variation of the tool orientation. Small variation is preferred as it implies that the orientation of the tool need not change rapidly along a tool path, and less demand for rapid response is exerted on the rotary axes of the machine.
Variation in the tool orientation can be evaluated using the rate of change of the unit tool vector T 0 j j, which can be readily computed using Eqs. (12), (18), and (22)   can be computed to reflect demands on the rotary axes for changes in the angular velocity. Here, L i and L iþ1 correspond to the arc lengths between these path points.

Comparison results
Three sets of comparisons are reported. To facilitate comparison, the first and second sets are based on the paths used by Farouki and Li [27]. Tool orientation variation and rotary axes movement are computed along paths on these surfaces. Results for RM I, RM II, and DRM are indexed as 1, 2, and 3, respectively, in Figs. 5-17.

Torus
The torus surface in Cartesian space is defined as with R ¼ 200 mm, d ¼ 100 mm, and u; v 2 0; 2p ½ , and the tool path to be used is defined in the u; v ð Þ parametric space as with n 2 0; 1 ½ . Grid size Dn ¼ 0:01 is used in the experiments.Tool vector T along the path is depicted in Fig. 5 for the three methods. Initial angles a 0 ¼ b 0 ¼ 5 are found to give reasonable tool orientation and are adopted.
For RM I with a fixed at 5 , the tool vector is close to the surface normal n. Notably, the actual tool location also depends on the inclination angle b, which varies along the path. Thus, the tool vector points almost horizontally at the start. Subsequently, the tool vector gradually changes together with the surface normal, and points almost vertically upwards near the end.
For RM II with b fixed at 5 , the tool vector simply stays close to the t-n plane. Initially, it also points almost horizontally, then it gradually tilts away from n as the tool moves along the path. For DRM, the tool vector falls on the m-n plane. It has the same initial orientation as in RM I and RM II. Then, it progressively inclines away from the t-n plane, and also tilts away from n along the path.
A comparison of the rate of change of tool orientation T 0 j j for the three methods is given in Fig. 6. For RM I, T 0 j j starts near 2 and then goes above 3 before coming down. For RM II, T 0 j j starts near 1.2 and decreases to zero, and then goes up to about 1.4. It appears that the tool tilt direction reverses, leading to the intermediate zero T 0 j j. T 0 j j also starts near 1.2 for DRM. However, its trend appears rather steady and T 0 j j stays close to 1 throughout the path. Thus, RM I has the largest T 0 j j. The average T 0 j j values are Flat-end tool orientation based on rotation-minimizing frame 263 similar for RM II and DRM, but DRM has more steady T 0 j j values and a smaller peak value than RM II.
Analytical representations of T 0 j j for the three methods are given in Eqs. (12), (18), and (22). According to Eq. (22), T 0 j j for DRM contains the term r cos a j n sin b À s g cos b À Á h only. Thus, the change direction is orthogonal to T. As a is typically small, the rate of change depends mainly on s g if b is small, and it also depends on j n if b is not.
Besides the above term, Eq. (12) for RM I contains the term r j n cos b þ s g sin b À Á w, which is orthogonal to that term. This suggests that RM I uses a larger T 0 j j than DRM. The change direction is also orthogonal to T as w is orthogonal to T. For RM II, the extra term in Eq. (18) is Àrj g sin ah. This indicates that the change direction is the same as that for DRM. While T 0 j j depends on j n and s g in DRM, it also depends on j g in RM II.
Rotary angles A and B and the corresponding acceleration terms S RA and S RB are shown in Figs. 7-10. The rotary angles start in the same manner for all the three methods because the tool orientation is the same initially (see Fig. 5). Subsequently, the angles assume quite different values as different methods are used to change the inclination and tilt angles. The acceleration terms also vary in different ways for the three methods (see Figs. 8 and 10). However, it is observed that DRM uses smaller acceleration than the other methods. This may be attributed to the facts that both the inclination and tilt angles are allowed to vary in DRM and that rotation minimization is adopted for Flat-end tool orientation based on rotation-minimizing frame 265 those changes. The term S RAB is explained in Fig. 11. This term reflects the combined acceleration of the individual rotary axes and also indicates that DRM uses the smallest acceleration. It can be seen that S RAB is the smallest for DRM, and S RAB values for RM II are slightly smaller than those for RM I.

Sphere
The sphere surface is defined as with R ¼ 200 mm, u 2 0; 2p ½ , and v 2 À 1 2 p; 1 2 p Â Ã . Tool path C n ð Þ is as defined in Eq. (29). The same grid size is adopted, and initial a 0 ¼ 45 and b 0 ¼ 0 . For the three methods, the tool vector along the path is given in Fig. 12. Again, although the initial tool orientation is the same (i.e., aligned to the surface normal), subsequent tool orientation is quite different because the inclination and tilt angles are changed using different ways in the three methods.
The rate of change of tool orientation T 0 j j is given in Fig. 13. The value of T 0 j j for RM I is large in the beginning and drops to a value close to that of DRM in the end. Conversely, the value of T 0 j j for RM II starts small but grows rapidly from mid-course onward. It seems that DRM has the smallest overall T 0 j j value. Moreover, T 0 j j is always smooth for DRM. Figure 14 clearly shows that DRM has a relatively steady axial acceleration and the smallest S RAB . Both RM I and RM II use similar acceleration all the way until about two-thirds of the path. Then, a surge in acceleration occurs for RM II so as to keep up with the planned tool orientation. This may be attributed to the local geometry, in combination with the configuration of the machine tool, which might not favor orienting the tool with a fixed tilt angle in that vicinity. The surge is not obvious from the tool orientation (see Fig. 12) and its rate of change (see Fig. 13). The value of S RAB is small for both DRM and RM I. However, that of T 0 j j is much larger for RM I.

Dome
The surface is defined as with R ¼ 100 mm and u; v 2 À1; 1 ½ . The tool path is C n ð Þ ¼ n; n ð Þ with n 2 0; 1 ½ . Initial a 0 and b 0 are 45 and 5 , respectively. The tool path and orientation are given in Fig. 15. The path starts at the top of the dome and ends at one of the base corners. For RM I, the tool axis maintains an angle of 45 to the surface normal along the path. RM II allows changes in the inclination angle. The tool axis gradually approaches the normal direction as the tool moves away from the dome center. In this instance, DRM has the same tool orientation as RM II. This happens when the Darboux frame is rotation-minimizing with respect to the surface normal n, and DRM and RM II have the same initial b 0 . For the current tool path, it seems that the key to reduce tool orientation changes is to allow the tool inclination angle to vary, which is the case in RM II and DRM but not in RM I. The rate of change of the tool orientation T 0 j j and the acceleration term S RAB are plotted in Figs. 16 and 17. In all cases, the values of T 0 j j and S RAB are relatively large in the beginning and small toward the end. The curves for RM II and DRM are exactly the same. The peak T 0 j j and peak S RAB values for RM I are a few times higher than those for RM II and DRM. They demonstrate that for some surface paths, varying a according to rotation minimization is critical to reducing the tool orientation changes and lightening the burdens on the rotary axes.

Further simulations
Free-form surfacing is common in five-axis machining. For machining free-form surfaces, the surfaces are typically modelled as spline surfaces. Here, part of the surface is modelled as a bicubic B-spline surface S u; v ð Þ with 0 u; v 1. The tool path for illustration is generated on the surface using u n ð Þ ¼ v n ð Þ ¼ n with n 2 0; 1 ½ . Tool orientation is rendered along path locations.
The initial tool orientation setup is a ¼ 20 and b ¼ 10 . Besides the three proposed RMF tool orientation methods, Sturz's method (where a is constant, which is also set at 20 , and b ¼ 0 as the tool axis is tilted toward the path tangent) has been adopted in the simulations. Tool orientation along the path is illustrated in Fig. 18 at locations of fixed increments of n, with the path starting from the left of the surface.
These simulations help elucidate how tool orientation can be adapted using these methods. It seems that RM I and Sturz's method exhibit similar tool orientation changes. This is because a is fixed at the same value for both methods. While b ¼ 0 for Sturz's method, b probably varied over a small range for RM I along this particular path. RM II and DRM also seem to have similar tool orientation changes, but they are smaller than those of RM I and Sturz's method, which is confirmed by checking the total tool orientation changes along the path (see Table 1). It is also confirmed that, among the four methods, DRM commanded the least amount of tool orientation changes during the tool path execution.

Conclusions
The concept of RMF was explored for tool orientation planning in the context of five-axis machining. Three RMF methods were presented: RM I (which allows varying the inclination angle), RM II (which allows varying the tilt angle), and DRM (which allows varying both the inclination and tilt angles). These methods specify how tool orientation can be adjusted along the tool path without unnecessary rotation of the tool frame. The derivative of the tool orientation was also derived for these methods.
These methods were applied to flat-end tool machining. Tool orientation was computed for surface paths on a torus, sphere, and dome. Changes in the tool orientation were always smooth for these methods, thus supporting the use of RMF for tool orientation planning. Overall, DRM had the smallest rate of change in tool orientation as well as the smallest acceleration along the rotary axes. The values were sometimes smaller for RM I, and at other times, smaller for RM II. This depends on the local geometry and machine tool kinematics. It is reasonable to expect that to reduce the tool orientation changes and axial acceleration, sometimes, it is more critical to vary the inclination angle, and at other times, the tilt angle. DRM outperforms the other methods because both angles are allowed to vary. Implementation of these tool orientation methods is rather straightforward as the computation only involves the local geometry of the surface and path. Notwithstanding issues such as collision avoidance, gouging, and cutting quality, which cannot be dealt with in this paper, RMF is an attractive tool for orientation planning. Further, DRM, is particularly useful in five-axis machining. The proposed methods can be extended in a future study to plan stable and smooth tool orientation within the constraints arising from these issues.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://crea tivecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.  Flat-end tool orientation based on rotation-minimizing frame 267