Numerical simulation and prediction of cutting forces in five-axis milling processes with cutter run-out

Radial cutter runout is a common issue in milling processes and has a direct effect on milling stability due to variations of resulting chip load and forces. This paper presents a new method to effectively model and predict the instantaneous cutting forces in 5-axis milling processes with radial cutter runout based on tool motion analysis. First, the undeformed chip thickness model taking runout effect into account is established under continuous change of cutter axis orientation by means of the sweep traces of cutter edges. Second, the engaged cutting edge is determined and cutting coefficients are subsequently calibrated. Finally, the method of identifying runout parameters from the measured cutting forces is proposed, and mechanistic method is then applied to predict the cutting force. Since this method is completely based on the relative motion analysis of tool-part, it can reduce the prediction errors of cutting forces effectively and is suitable for generic rotation cutters. Several validation examples are given under different cutting conditions to prove its effectiveness and accuracy. The results reveal that the developed method can predict the cutter forces with a high accuracy and has the ability to be used in simulations and optimizations of five-axis machining.     

Phase width analysis of cutting forces considering bottom edge cutting and cutter runout calibration in flat end milling of titanium alloy

This paper is concerned with the combined cutting effects of both flank and bottom edges based on a systematic study of the cutting force in flat end milling of the titanium alloy. Besides the flank edge, the bottom edge of the cutter is also found to be an important factor influencing the cutting force distributions and can lead to uniform phase widths for non-zero cutting forces even under considerable cutter runout. One such phenomenon of uniform phase width induced by the bottom edge for the cutting force is deeply revealed. To do this, the models for characterizing the cutting force coefficients related to both edges are established based on the measured instantaneous cutting forces, and cutter runout is considered in the computation of process geometry parameters such as cutter/workpiece engagements and instantaneous uncut chip geometry parameters. Novel algorithms for identifying the cutter runout parameters and the bottom uncut chip width are also developed. Results definitely show that the flank cutting force coefficients can be treated as constants and that size effect obviously exists in the bottom cutting force coefficients that can be characterized by a power function of the bottom uncut chip width.The proposed model is validated through a comparative study with the existing model and experiments. From the outcomes of the current work, it is clearly seen that the prediction of cutting forces for titanium alloy can resort to the proposed model instead of traditional ones.     

Theoretical modelling and simulation of cutting forces in face milling with cutter runout

A new approach to theoretical modelling and simulation of cutting forces in face milling is presented. Based on a predictive machining theory, the action of a milling cutter is modeled as the simultaneous actions of a number of single-point cutting tools. The milling forces are predicted from the workpiece material properties, cutter parameters, tooth geometry, cutting conditions and types of milling. The properties of the workpiece material are considered as functions of strain, strain-rate and temperature in the cutting region. It takes into account the effect of the intermittent contact between each milling tooth and the workpiece on the temperature in the cutting region. It also takes into account the effect of cutter runout on the undeformed chip thickness. Milling experiments have been conducted to verify the proposed model. Good agreements between the experimental and simulated results are presented.     

Cutting force prediction of sculptured surface ball-end milling using Z-map

The cutting force in ball-end milling of sculptured surfaces is calculated. In sculptured surface machining, a simple method to determine the cutter contact area is necessary since cutting geometry is complicated and cutter contact area changes continuously. In this study, the cutter contact area is determined from the Z-map of the surface geometry and current cutter location. To determine cutting edge element engagement, the cutting edge elements are projected onto the cutter plane normal to the Z-axis and compared with the cutter contact area obtained from the Z-map. Cutting forces acting on the engaged cutting edge elements are calculated using an empirical method. Empirical cutting mechanism parameters are set as functions of cutting edge element position angle in order to consider the cutting action variation along the cutting edge. The relationship between undeformed chip geometry and the cutter feed inclination angle is also analyzed. The resultant cutting force is calculated by numerical integration of cutting forces acting on the engaged cutting edge elements. A series of experiments were performed to verify the proposed cutting force estimation model. It is shown that the proposed method predicts cutting force effectively for any geometry including sculptured surfaces with cusp marks and a hole.     

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.     

Milling force convolution modeling for identification of cutter axis offset

This paper discusses the application of a convolution integral force model to the identification of the geometry of cutter axis offset in milling operations. This analysis builds upon the basis of linear decomposition of elemental local cutting forces into a nominal component and an offset-induced component. The convolution of each elemental local cutting force component with the chip width density in the context of cutter angular position provides an integral expression for the total cutting forces. By virtue of the convolution integration property, the total cutting forces in the frequency domain can be derived as closed-form functions of the cutting pressure constants, various cutting conditions, as well as the cutter offset geometry. Subsequently, the magnitude and phase angle of cutter axis offset are shown to be algebraic and explicit functions of the Fourier series coefficients of cutting forces at the spindle frequency. Following the theoretical analysis, experimental study is discussed to illustrate the implementation procedure for offset identification, and frequency domain data are presented to verify the analytical results. Potential industrial applications of this work include the real-time monitoring of dynamic cutter runout and the in-process compensation for the loss of tolerance or finish using automatic controls based on the feedback information of offset magnitude and phase angle.     

Simplified and efficient calibration of a mechanistic cutting force model for ball-end milling

Accurate evaluation of the empirical coefficients of a mechanistic cutting force model is critical to the reliability of the predicted cutting forces. This paper presents a simplified and efficient method to determine the cutting force coefficients of a ball-end milling model. The unique feature of this new method is that only a single half-slot cut is to be performed to calibrate the empirical force coefficients that are valid over a wide range of cutting conditions. The instantaneous cutting forces are used with the established helical cutting edge profile on the ball-end mill. The half-slot calibration cut enables successive determination of the lumped discrete values of the varying cutting mechanics parameters along the cutter axis whereas the size effect parameters are determined from the known variation of undeformed chip thickness with cutter rotation. The effectiveness of the present method in determining the cutting force coefficients has been demonstrated experimentally with a series of verification test cuts.     

Sculpture surface machining: a generalized model of ball-end milling force system

A new mechanistic model is presented for the prediction of a cutting force system in ball-end milling of sculpture surfaces. The model has the ability to calculate the workpiece/cutter intersection domain automatically for a given cutter location (CL) file, cutter and workpiece geometries. Furthermore, an analytical approach is used to determine the instantaneous chip load (with and without runout) and cutting forces. In addition to predicting the cutting forces, the model also employs a Boolean approach for a given cutter, workpiece geometries, and CL file to determine the surface topography and scallop height variations along the workpiece surface which can be visualized in 3-D. The results of model validation experiments on machining Ti-6A1-4V are also reported. Comparisons of the predicted and measured forces as well as surface topography show good agreement.     

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.     

Systematic study on cutting force modelling methods for peripheral milling

This paper systematically studies the cutting force modelling methods in peripheral milling process in the presence of cutter runout. Emphasis is put on how to efficiently calibrate the cutting force coefficients and cutter runout. Mathematical derivations and implementation procedures are carried out based on the measured cutting force or its harmonics from Fourier transformation. Five methods are presented in detail. In the first three methods the cutting force coefficients are assumed to be constants whereas in the last two they are taken as functions of instantaneous uncut chip thickness. The first method and the fifth one are taken from literatures for comparison. The second, the third and the fourth methods are original contributions, which are carried out with optimization ideas. The second method proceeds using the first and Nkth harmonic forces as the source signal while the third and the fourth are derived based on the measured cutting forces and its first harmonics. The engagement of the cutter with the workpiece is considered in these three new calibration procedures without the requirement of a prior knowledge of the actual cutter runout. Comparisons among the calibrated results from different methods are made to study the limitations and consistency of the presented methods. Experiments are also conducted to show the prediction ability of all methods.     

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.     

Modelling and simulation of chatter in milling using a predictive force model

A predictive time domain chatter model is presented for the simulation and analysis of chatter in milling processes. The model is developed using a predictive milling force model, which represents the action of milling cutter by the simultaneous operations of a number of single-point cutting tools and predicts the milling forces from the fundamental workpiece material properties, tool geometry and cutting conditions. The instantaneous undeformed chip thickness is modelled to include the dynamic modulations caused by the tool vibrations so that the dynamic regeneration effect is taken into account. Runge–Kutta method is employed to solve the differential equations governing the dynamics of the milling system for accurate solutions. A Windows-based simulation system for chatter in milling is developed using the predictive model, which predicts chatter vibrations represented by the tool-work displacements and cutting force variations against cutter revolution in both numerical and graphic formats, from input of tool and workpiece material properties, cutter parameters, machine tool characteristics and cutting conditions. The system is verified with experimental results and good agreement is shown.     

Development of a virtual machining system, part 1: approximation of the size effect for cutting force prediction

In this three-part paper, components of a virtual machining system for evaluating and optimizing cutting performance in -axis NC machining are presented. Part 1 describes a new method of calculating cutting-condition-independent coefficient and its application to the prediction of cutting forces over a wide range of cutting conditions. The prediction of the surface form error and transient cutting simulations, described in Parts 2 and 3, respectively, can be effectively performed based on the cutting force model with the improved size effect model that is presented in Part 1.

The relationship between the instantaneous uncut chip thickness and the cutting coefficients is calculated by following the movement of the center position of the cutter, which varies with nominal feed, cutter deflection and runout. The salient feature of the presented method is that it determines the cutting-condition-independent coefficients using experimental data processed for one cutting condition. The direct application of instantaneous cutting coefficient with size effects provides more accurate predictions of the cutting forces. A systematic comparison of the predicted and measured cutting forces over a wide range of cutting conditions confirms the validity of the proposed mechanistic cutting force and size effect models.     

Feed rate scheduling model considering transverse rupture strength of a tool for 3D ball-end milling

This paper presents an analytical model of off-line feed rate scheduling to determine desired feed rates for 3D ball-end milling. Off-line feed rate scheduling is presented as the advanced technology to regulate cutting forces through change of feed per tooth, which directly affects variation of uncut chip thickness. In this paper, the uncut chip thickness is calculated by following the movement of the position of the cutter center, which is determined by runout and cutter deflection. Also, since the developed cutting force model uses the cutting-condition-independent coefficients, off-line feed rate scheduling can be effectively performed regardless of continuous change of cutting conditions. Transverse rupture strength of the tool is used to determine the reference cutting force at which resultant cutting forces are regulated through feed rate scheduling. Experiments validated that the presented feed rate scheduling model reduced machining time drastically and regulated cutting forces at the reference cutting force.     

New procedures for calibration of instantaneous cutting force coefficients and cutter runout parameters in peripheral milling

In this paper, new procedures are proposed to calibrate the instantaneous cutting force coefficients and the cutter runout parameters for peripheral milling. By combining with optimization algorithm, i.e., the Nelder–Mead simplex method, detailed calibration schemes are derived for a mechanistic cutting force model in which the cutting force coefficients are described as the exponential functions of the instantaneous uncut chip thickness. Three different cutter runout models are considered in the calculation of instantaneous uncut chip thickness. Only one or two tests are required to perform the calibration. Experimental verifications are also conducted to validate the proposed procedures, and the results show that they are efficient and reliable. To see the effect of different runout models on milling process, comparisons among the predicted results under a wide range of cutting parameters are made to study the consistency and limitations of different models. It is found that the radial cutter runout model is a recommendable one for cutting force modelling.     

The prediction of dimensional error for sculptured surface productions using the ball-end milling process. Part 1: Chip geometry analysis and cutting force prediction

The study of machining errors caused by tool deflection in the balkend milling process involves four issues, namely the chip geometry, the cutting force, the tool deflection and the deflection sensitivity of the surface geometry. In this paper, chip geometry and cutting force are investigated. The study on chip geometry includes the undeformed radial chip thickness, the chip engagement surface and the relationship between feed boundary and feed angle. For cutting force prediction, a rigid force model and a flexible force model are developed. Instantaneous cutting forces of a machining experiment for two 2D sculptured surfaces produced by the ball-end milling process are simulated using these force models and are verified by force measurements. This information is used in Part 2 of this paper, together with a tool deflection model and the deflection sensitivity of the surface geometry, to predict the machining errors of the machined sculptured surfaces.     

An investigation of cutting forces in machining with worn ball-end mill

Reliable tool wear monitoring technique is one of the important aspects for achieving an integrated and self-adjusting manufacturing system. In this study, an analytical model is proposed to estimate the cutting forces, the tool geometry, and the chip geometry in relation to the flank wear, when milling with a ball-end mill. Modeling is based on thermomechanical modelling of oblique cutting. The worn tool geometry is decomposed into a series of axial elementary cutting edges. At any active tooth element, the flank wear geometry is calculated and the chip formation is obtained from an oblique cutting process characterised by local undeformed chip section and local cutting angles. Coated carbide ball-end tool, and a titanium workpiece material have been considered in this paper. The results found by using developed models have shown good agreement with experimental results.     

Mechanistic model for the reaming process with emphasis on process faults

A mechanistic model is developed to predict the cutting forces for an arbitrary reamer geometry. Inputs to the model include: tool geometry, feed, speed, initial hole geometry, and process faults including parallel offset runout, spindle tilt, their respective locating angles, and tool/hole axis misalignment. Given these input parameters, the model predicts torque, thrust, and radial forces. The cutting edges of the reamer are divided into elements and the elemental forces are determined from a fundamental oblique cutting model. The model is calibrated over a range of feed, speed, and varying tool geometry. Model validation tests were conducted and model predictions match experimental data well. The effects of process faults on cutting forces are examined through model-based studies. It was found that parallel offset runout, spindle tilt, spindle tilt locating angle, and tool/hole axis misalignment have significant effects on the radial forces. These radial forces are shown to be correlated to hole quality.