Mechanics and dynamics of general milling cutters.: Part I: helical end mills

A variety of helical end mill geometry is used in the industry. Helical cylindrical, helical ball, taper helical ball, bull nosed and special purpose end mills are widely used in aerospace, automotive and die machining industry. While the geometry of each cutter may be different, the mechanics and dynamics of the milling process at each cutting edge point are common. This paper presents a generalized mathematical model of most helical end mills used in the industry. The end mill geometry is modeled by helical flutes wrapped around a parametric envelope. The coordinates of a cutting edge point along the parametric helical flute are mathematically expressed. The chip thickness at each cutting point is evaluated by using the true kinematics of milling including the structural vibrations of both cutter and workpiece. By integrating the process along each cutting edge, which is in contact with the workpiece, the cutting forces, vibrations, dimensional surface finish and chatter stability lobes for an arbitrary end mill can be predicted. The predicted and measured cutting forces, surface roughness and stability lobes for ball, helical tapered ball, and bull nosed end mills are provided to illustrate the viability of the proposed generalized end mill analysis.     

Effect of cutter runout on process geometry and forces in peripheral milling of curved surfaces with variable curvature

This paper presents a novel method for cutting force modeling related to peripheral milling of curved surfaces including the effect of cutter runout, which often changes the rotation radii of cutting points. Emphasis is put on how to efficiently incorporate the continuously changing workpiece geometry along the tool path into the calculation procedure of tool position, feed direction, instantaneous uncut chip thickness (IUCT) and entry/exit angles, which are required in the calculation of cutting force. Mathematical models are derived in detail to calculate these process parameters in occurrence of cutter runout. On the basis of developed models, some key techniques related to the prediction of the instantaneous cutting forces in peripheral milling of curved surfaces are suggested together with a whole simulation procedure. Experiments are performed to verify the predicted cutting forces; meanwhile, the efficiency of the proposed method is highlighted by a comparative study of the existing method taken from the literature.     

Cutting force prediction in ball end milling of sculptured surface with Z-level contouring tool path

This paper presents an approach to predict cutting force in 3-axis ball end milling of sculptured surface with Z-level contouring tool path. The variable feed turning angle is proposed to denote the angular position of feed direction within tool axis perpendicular plane. In order to precisely describe the variation of feed turning angle and cutter engagement, the whole process of sculptured surface milling is discretized at intervals of feed per tooth along tool path. Each segmented process is considered as a small steady-state cutting. For each segmented cutting, the feed turning angle is determined according to the position of its start/end points, and the cutter engagement is obtained using a new efficient Z-map method. Both the chip thickness model and cutting force model for steady-state machining are improved for involving the effect of varying feed turning angle and cutter engagement in sculptured surface machining. In validation experiment, a practical 3-axis ball end milling of sculptured surface with Z-level contouring tool path is operated. Comparisons of the predicted cutting forces and the measurements show the reliability of the proposed approach.     

球头铣刀铣削斜面的三维有限元仿真研究

在能源动力、汽车、航空航天、模具制造等关键零部件的加工过程中,球头铣刀因其特有的刀具几何结构,常作为零件加工的最终成型刀具。考虑到在球头铣刀立铣加工中不同的刀具与工件相对姿态会对切削过程产生不同的影响,本文研究切屑形成和不同走刀方式下切削过程中各物理量(切削力、切削温度等)的变化情况,结合有限元仿真技术在切削加工中的应用,建立硬质合金球头铣刀铣削斜面的有限元模型,模拟相同切削参数下,八种不同走刀方式的球头铣削过程,分析刀具切入切出工件时切屑的形成过程,探究切削力和切削温度的变化规律。仿真结果表明:不同的走刀方式,平均切削合力各不相同,同时切屑和工件的最大切削温度也出现较大差异,而斜坡上坡逆铣的走刀方式所对应的平均切削合力和最大切削温度均最优。     

Improving Cycle Time in Sculptured Surface Machining Through Force Modeling

In this paper, an enhanced mathematical model is presented for the prediction of cutting force system in ball end milling of sculptured surfaces. This force model is also used as the basis for off-line feed rate scheduling along the tool path in order to decrease the cycle time in sculptured surface machining. As an alternative for setting a constant feed rate all along the tool path in rough machining of sculptured surfaces, resultant cutting forces are aimed to be kept under a pre-set threshold value along the tool path by off-line scheduled piecewise variable feed rates. In this paper, it is shown that machining time, depending on complexity of sculptured surfaces, can be decreased significantly by scheduling feed rate along the tool path. The model is tested under various cutting conditions and some of the results are also presented and discussed in the paper.     

雕刻表面球形铣削的分段变进给率加工

雕刻表面球形铣削加工中,根据刀头的受力模型,可以计算整个刀具路径上刀头的受力.刀具路径上加工深度变化时,若采用恒进给率,则刀具受力是变化的,为了安全保险起见,往往按最大受力来选择较低的恒进给率,整个加工效率很低.采用变进给率加工,即加工深度小时提高进给率,使刀具在整个路径上受力均匀.本文给出了路径上力的计算仿真方法,以及分段进给率的计算方法,可以指导实际的加工过程,这样既保证了加工过程的安全,又提高了生产率.     

Five-axis milling mechanics for complex free form surfaces

Accurate and fast prediction of machining forces is important in high performance cutting of free form surfaces that are commonly used in aerospace, automotive, biomedical and die/mold industries. This paper presents a novel and generalized approach for prediction of cutting forces in five-axis machining of parts with complex free form surfaces. Engagement simulations between cutter and part are performed precisely along the tool path by a recently developed boundary representation method. Moreover, mathematical model for five-axis milling mechanics is developed for any given solid model of parts with complex free form surfaces. Theoretical simulations and experimental validations show that cutting forces are predicted fast and precisely for five-axis machining of complex free form surfaces.     

Chip volume prediction using a numerical control verification model

This paper presents a prediction model of chip volume per tooth in milling operations using a numerical control (NC) verification model that can describe actual cutting. This is important from the viewpoint of cutting force calculation and prediction of the chip load on a cutter tooth. In this paper, the cutter swept volumes along the path were calculated from the information obtained from NC codes as well as the feed-rate. The chip volume per tooth was kinetically modeled for the machining state of free-formed surfaces with flat, ball and round endmills, respectively. To investigate the effectiveness of the proposed model, simulation studies were conducted. Following this simulation, the chip volume per tooth precisely mimicked real machining.     

Closed form formulation of cutting forces for ball and flat end mills

In this paper, a set of basis functions for the calculation of cutting forces in milling are introduced. It is shown that the cutting force at any tool position can be determined by the linear combination of the force basis functions. The method is based on the projection of the chip load area onto the reference coordinate planes. Due to the analytical integration of the cutting forces along the cutter edge, the developed closed form equations provide a fast means of calculating the cutting forces. The validity of the method is experimentally verified for both flat and ball end mills.     

Estimation of cutter deflection and form error in ball-end milling processes

This paper presents a method to analyze the 3-dimensional form error of a ball-end milled surface due to the elastic compliance of the cutting tool. In order to estimate the form error in various cutting modes, the cutting force and the cutter deflection models including the effect of the surface inclination were established. The cutting forces were calculated by using the cutter contact area determined from the Z-map of the surface geometry and the current cutter location. The tool deflection responding to the cutting force was then calculated by considering the cutter and the holder stiffness. The cutter was modeled as a cantilever beam consisting of the shank and the flute. The stiffness of the holder was measured experimentally. Various experimental works have been performed to verify the validity of the proposed model. It is shown that the proposed method is capable of accurate prediction of cutting forces and the surface form error.     

Process geometry modeling with cutter runout for milling of curved surfaces

Prediction of cutting forces and machined surface error in peripheral milling of curved geometries is non-trivial due to varying workpiece curvature along tool path. The complexity in this case, arises due to continuously changing process geometry as workpiece curvature varies along tool path. In the presence of cutter runout, the situation is further complicated owing to changing radii of cutting points. The present work attempts to model process geometry in machining of curved geometries and in the presence of cutter runout. A mathematical model computing process geometry parameters which include cutter/workpiece engagements and instantaneous uncut chip thickness in the presence of cutter runout is presented. The developed model is more realistic as it accounts for interaction of cutting tooth trajectory with that of preceding teeth trajectories in computing process geometry. Computer simulation studies carried for this purpose has shown that it is essential to account for teeth trajectory interactions for accurate prediction of process geometry parameters. This aspect is further confirmed with machining experiments, which were conducted to validate this aspect. From the outcomes of present work, it is clearly seen that the computation of process geometry during machining of curved geometries and in presence of cutter runout is not straightforward and requires a systematic approach as presented in this paper.     

Evaluation of cutter orientations in high-speed ball end milling of cantilever-shaped thin plate

This work was carried out in order to evaluate machinability in high-speed ball end milling of thin cantilever-shaped aluminum plate, and the effect of tool path, tool angle on geometric accuracy and surface integrity. The experiments were performed at tilt angles 0°, 15°, 45° of workpieces using four cutter orientations (horizontal outward, horizontal inward, vertical outward, vertical inward). Workpiece deflection, surface profile and machined surface have been measured. The results show that when machining cantilever-shaped thin plate, the best cutting strategy is vertical outward at a workpiece tilt angle of 45°.     

Effect of direction of parameterization on cutting forces and surface error in machining curved geometries

The present paper investigates the effect of two variables, namely direction of parameterization and cutter diameter on process geometry, cutting forces, and surface error in peripheral milling of curved geometries. In machining of curved geometries where the curvature varies continuously along tool path, the process geometry variables, namely feed per tooth, engagement angle, and maximum undeformed chip thickness too vary along tool path. These variations will be different when a given geometry is machined from different parametric directions and with different cutter diameters. This difference in process geometry variations result in changed cutting forces and surface error along machined path. This aspect has been studied for variable curvature geometries by machining from both parametric directions and using cutters of different diameter. The computer simulation studies carried out show considerable amount of shift in the location of peak cutting forces with the change in cutting direction and cutter diameter, particularly in concave regions of workpiece geometry. A new parameter γ that relates the instantaneous curvature of workpiece with cutter radius is defined. The larger value of γ is an indicator of greater shift in the location of peak forces from the point of maximum curvature on the workpiece. The simulation results are validated by carrying out machining experiments with curved workpiece geometry and are found to be in good agreement.     

An improved process simulation system for ball-end milling of sculptured surfaces

A simulation system is developed in this paper, which deals with the geometry and mechanics of machining with ball-end milling cutters. The geometry of the workpiece, the cutter, and the cutter/workpiece engagement is modeled using a geometric simulation system. This module uses a commercial solid modeler (ACIS) as a geometric engine and automatically extracts the critical geometric information required for the physical simulation system. To calculate the instantaneous cutting forces, a new mechanistic force model is developed. This force model takes into account the variations of the cutting coefficients along the cutting edge, and considers the variations of the rake angle and the chip flow direction on the rake face. The calibration of the developed model is performed for half-immersion ball-end milling operation. The applicability of the developed system is verified experimentally for various up-hill angles. It is shown that as the up-hill angle increases, the ball-nose tip engagement decreases which in turn significantly affects the magnitude of the resultant forces. Also, lower cutting forces and powers are experienced if cutting with the vicinity of the tool tip is avoided.     

Dynamic response analysis in end milling using pre-twisted beam finite elements

This paper develops an analytical model for estimating the dynamic responses in end milling, i.e. dynamic milling cutter deflections and cutting forces, by using the finite-element method along with an adequate end milling-cutting force model. The whole cutting system includes the spindle, the bearings and the cutter. The spindle is modelled structurally with the Timoshenko-beam element, the milling cutter with the pre-twisted Timoshenko-beam element due to its special geometry, and the bearings with lumped springs and dampers. Because the damping matrix in the resulting finite-element equation of motion for the whole cutting system is not one of proportional damping due to the presence of bearing damping, the state-vector approach and the convolution integral is used to find the solution of the equation of motion. To assure the accuracy of prediction of the dynamic response, the associated cutting force model should be sufficiently precise. Since the dynamic cutting force is proportional to the chip thickness, a quite accurate alogorithm for the calculation of the variation of the chip thickness due to geometry, run-out and spindle-tool viration is developed. A number of dynamic cutting forces and tool deflections obtained from the present model for various cutting conditions are compared with the experimental and analytical results available in the literature, good agreement being demonstrated for these comparisons. The present model is useful, therefore, for the prediction of end milling instability. Also, the tool deflections obtained using the pre-twisted beam element are found to be smaller than those obtained using the straight beam element without pre-twist angle. Hence neglecting the pre-twist angle in the structural model of the milling cutter may overestimate the tool deflections.     

Increasing productivity in sculpture surface machining via off-line piecewise variable feedrate scheduling based on the force system model

Sculpture surface machining is a critical process commonly used in various industries such as the automobile, aerospace, die/mold industries. Since there is a lack of scientific tools in practical process planning stages, feedrates for CNC machining are selected based on the trial errors and previous experiences. In the selections of the process parameters, production-planning engineers are conservative in order to avoid undesirable results such as chipping, cutter breakage or over-cut due to excessive cutter deflection. Currently, commonly used CAD/CAM programs use only the geometric and volumetric analysis, but not the physics of the processes, and rely on experience based cutting tool database and users’ inputs for selection of the process parameters such as feed and speed. Usually, the feeds and cutting speeds are set individual constant values all along the roughing, semi-finishing, and finishing processes. Being too conservative and setting feedrate constant all along the tool path in machining of sculpture surfaces can be quite costly for the manufacturers. However, a force model based on the physics of the cutting process will be greatly beneficial for varying the feedrate piecewise along the tool path.The model presented here is the first stage in order to integrate the physics of the ball-end milling process into the selection of the feeds during the sculpture surface machining. Therefore, in this paper, an enhanced mathematical model is presented for the prediction of cutting force system in ball end milling of sculpture surfaces. This physical force model is used for selecting varying and ‘appropriate’ feed values along the tool path in order to decrease the cycle time in sculpture surface machining. The model is tested under various machining conditions, and some of the results are also presented in the paper.     

Efficient calibration of instantaneous cutting force coefficients and runout parameters for general end mills

This paper aims at developing a unified approach to identify the cutting force coefficients together with the cutter runout parameters for general end mills such as cylindrical, ball, bull nose ones, etc. The cutting forces that are modeled using the instantaneous cutting force coefficients are analyzed and separated into two terms: a nominal component independent of the runout and a perturbation component induced by the runout. The nominal component enables the calibration of the instantaneous cutting force coefficients whereas the runout parameters are determined from the perturbation component. The validity of the present method is demonstrated with simulation and experimented data.     

Milling error prediction and compensation in machining of low-rigidity parts

The paper reports on a new integrated methodology for modelling and prediction of surface errors caused by deflection during machining of low-rigidity components. The proposed approach is based on identifying and modelling key processing characteristics that influence part deflection, predicting the workpiece deflection through an adaptive flexible theoretical force-FEA deflection model and providing an input for downstream decision making on error compensation. A new analytical flexible force model suitable for static machining error prediction of low-rigidity components is proposed. The model is based on an extended perfect plastic layer model integrated with a FE model for prediction of part deflection. At each computational step, the flexible force is calculated by taking into account the changes of the immersion angles of the engaged teeth. The material removal process at any infinitesimal segment of the milling cutter teeth is considered as oblique cutting, for which the cutting force is calculated using an orthogonal–oblique transformation. This study aims to increase the understanding of the causes of poor geometric accuracy by considering the impact of the machining forces on the deflection of thin-wall structures. The reported work is a part of an ongoing research for developing an adaptive machining planning environment for surface error modelling and prediction and selection of process and tool path parameters for rapid machining of complex low-rigidity high-accuracy parts.     

Mechanics and dynamics of thread milling process

This paper presents the mechanics and dynamics of thread milling operations. The tool follows a helical path around the wall of the pre-machined hole in thread milling, which has varying tool-part engagement and cut area during one threading cycle. The variation of cut area that reflects the kinematics of threading as well as structural vibrations is modeled along the helical, threading path. The mechanics of the process are first experimentally proven, followed by the formulation of dynamic thread milling which is periodic in threading cycle, in a semi-discrete time domain. The stability of the operation is predicted as a function of spindle speed, axial depth of cut, cutter path and tool geometry. The mechanics and stability models are experimentally proven in opening M16×2 threads with a five-fluted helical tool on a Steel AISI1045 workpiece.     

Prediction of cutting forces in milling of circular corner profiles

This paper proposes an approach to predict the cutting forces in peripheral milling of circular corner profiles in which varying radial depth of cut is encountered. The geometric relationship between an end mill and the corner profile is investigated and a mathematical model is presented to describe the different phases of the cutter/workpiece contact. The milling process for circular corner is discretized into a series of steady-state cutting processes, each with different radial depth of cut determined by the instantaneous position of the end mill relative to the workpiece. A time domain analytical model of cutting forces for the steady-state machining conditions is introduced to each segmented process for the cutting force prediction. The predicted cutting forces can be calculated in terms of tool/workpiece geometry, cutting parameters and workpirece material property, as well as the relative position of the tool to workpiece. Experiments are conducted and the measured forces are compared to the predictions for the verification of the proposed method.