i-manager's Journal on Mechanical Engineering

View PDF

Volume :6 No :1 Issue :-2016 Pages :24-33

Prediction of Cutting Forces in Turning of AISI-4140 Steel using FEA and Experimental Validation

M. Sivaramakrishnaiah *   P. Nanda Kumar **  G. Ranga Janardhana ***
* Assistant Professor, Department of Mechanical Engineering, Sri Venkateswara College of Engineering & Technology, Chittoor, A.P, India
** Professor, Department of Mechanical Engineering, NBKR institute of Science and Technology, Vidyanagar, Nellore, AP, India.
*** Professor, Department of Mechanical Engineering, Jawaharlal Nehru Technological University, Kakinada, AP, India.

Abstract

Turning is one of the most widely used machining techniques. Predictions of important process variables of turning such as, cutting forces and stress distributions play significant role in designing turning process parameter and optimizing cutting conditions. The present work is focused on investigating the effect of process parameters on feed, cutting, thrust forces and prediction of forces in turning of AISI 4140 steel based on FEA method. Analysis of 3D unsteady state forces in a metal cutting process is carried out by using a Deform 3D FEA code. The cutting speed, feed and depth of cut were used as the process parameters and the cutting forces are measured and predicted. This project covers a study on modeling and simulation of cutting forces in turning of AISI 4140 steel using coated and uncoated tungsten carbide cutting tool by finite element technique and finally with experimental results. The FEA is used to analyze cutting forces, stress and strain in turning. The simulation results are compared with those of the experimental results and found that both results are in close agreement with each other.

Keywords :

  • FEA,
  • Tungsten Carbide,
  • AISI 4140 Steel,
  • Tin Coating,
  • Turning

Introduction

The turning process used widely in industry, has countless applications. Traditionally, the process is used to reduce the diameter of a cylindrical work piece, or to change a work piece of non-circular cross-section. This is done by rotating the work piece about the machine's axis of spindle and removing the material with the cutting tool, fed in the perpendicular direction. For the past fifty years, metal cutting researchers, for simulation purposes, have developed many modeling techniques including analytical techniques, slip-line solutions, empirical approaches and finite element techniques. In recent years, the finite element analysis has particularly become the main tool for simulating metal cutting processes. Finite element analysis is widely used for calculating the stress, strain, strain-rate distributions in the primary, secondary and tertiary sub-cutting zones. In consequence, in the tool, chip and work piece, as well as cutting forces, elasto-plastic deformation (shear angles and chip thickness), chip formation and possibly its breaking can be determined faster than using costly and time consuming experiments. In this work, FEA based models are developed that are able to predict the effect of various process variables on the performance measures such as cutting forces, stresses and strains.

1. Methodology

1.1 Turning Process

Turning is a very important machining process in which a single-point cutting tool removes material from the surface of a rotating cylindrical work piece. The cutting tool is fed linearly in a direction parallel to the axis of rotation. Turning is carried out on a lathe that provides the power to turn the work piece at a given rotational speed and to feed the cutting tool at a specified rate and depth of cut. Therefore, the effect of five cutting parameters, i.e. cutting speed, feed rate, depth of cut, tool geometry (nose radius) and cutting condition on the cutting forces, stress and strains need to be studied in a turning operation.

Turning process is performed to modify shape, dimension, and the surface finish of a workpiece by removing the unwanted material in the form of chips. The theory of metal cutting and chip formation are complicated, where plasticity, thermodynamic and mathematical analysis are involved.

1.2 Finite Element Analysis of Turning Process

Finite Element Analysis is a most useful and accurate approach for the determination of field variables because of the advancements in computational and processing power of computers. Thus it is almost used for all the Computer Aided Design methodologies. Applications range from deformation and stress analysis to field analysis of heat flux, fluid flow, magnetic flux, seepage and other flow problems. In this method of analysis, a complex region defining a continuum is discredited into simple geometric shapes called finite elements. The present work is based on the application of finite element for determination of cutting forces in single point cutting tool for turning operation. Once the model is developed for determination of cutting forces for single point cutting tool, it can also be implemented for other multipoint processes like drilling, milling or grinding also.

Modeling 3D cutting process using finite element technique is an ongoing research activity due to significant cost savings and offers insights into the process which are not easily measured in experiments. In particular, heat transfer and the modeling of cutting process require careful consideration in any modeling activity. This paper presents approaches for modeling the turning process for AISI-4140 type of steel. In this study, Deform 3D Finite Element Analysis software is used to study the effects of cutting speed, feed rate, and depth of cut on cutting forces. The work-piece is modeled as Elastic-plastic material to take thermal, elastic, plastic effect. Work-piece is represented by a linear model with different length for each condition. The short material length was chosen to save computational time without compromising the model integrity as heat generation in machining in small areas around the cutting zone. The work-piece shape is constructed by the DEFORM machining module, and includes geometry created by a previous tool pass, including appropriate depth of cut and nose radius details. An unstructured tetrahedral finite element mesh was generated using DEFORM's automatic mesh generation system. Re-meshing parameters, including minimum element size, and parameters for adaptive mesh definition are set within the system. For these procedures of the simulations, a minimum element size of 0.25 of the feed rate was specified. The total number of elements are 100,000 to 130,000 depending upon work-piece size. The Lagrangian calculation embeds a computational mesh in the material domain and solves for the position of the mesh at discrete points in time. An incremental Lagrangian formulation with an implicit integration method designed for large deformation simulations is used to simulate the cutting process. The solver used was the sparse matrix with a direct integration method, because the conjugategradient offers an improved computational speed but less stability in convergence.

1.3 Tool and Work Material

The work piece considered is AISI 4140 Steel material and the tungsten carbide cutting tool is used for machining. Thermal and mechanical properties of tool and work piece material are shown in Table 1. AISI 4140 steel has 1% chromium - molybdenum and is medium hardenability, general purpose high tensile steel, generally supplied hardened and tempered in the tensile range of 850 - 1000 Mpa.

Table 1. Thermal and Mechanical Properties of the Work Piece and Tool Material

2. Literature Review

A Review of literature describes the study of the cutting forces of turning process with the help of finite element analysis and experimental results. Typical approaches for numerical modeling of metal cutting processes is Lagrangian and Eulerian techniques, as well as a combination of both called an arbitrary Lagrangian– Eulerian formulation [1, 2]. It may be noticed that, all these methods are mathematically equivalent. The major feature of Lagrangian formulation used in this study is that the mesh attached to the workpiece.

Also, the finite element analysis was performed by Johnson-Cook's constitutive equation with three different sets of material constants (found by the application of several methods) is implemented in the FE model to study the behavior of Ti6Al4V alloy during the machining process in conventional and high speed regimes [3]. Demand for higher productivity and good quality for machining parts has encouraged many researchers to study the effects of machining parameters using FEM simulation using either two or three dimension version [4].

Finite element method by Asmaa A. Kawi [5], is a successful technique to perform analysis to estimate cutting temperatures, a possibility of developing temperature forms adequately representing metal cutting temperature as a Polynomial model of third, fourth and fifth degree with time that gives steady state temperature and for the four alloy steels used and different operation conditions. All alloys have a severe increasing temperature with increasing feed rate, while it looks less sharp with increasing cutting speed. Also the ratio of the number of nodes have maximum temperature for any operating condition and any alloy used with respect to the total number of nodes is less than 1%.

Mofid Mahdi, and Liangchi Zhang [6], considered the chip breaking and developed a 2D cutting force model with the finite element method. The variation of the cutting force was investigated carefully against both the cutting conditions and the anisotropy of the material with the following development: (a) a constitutive model of a homogeneous anisotropic elastic material under plane deformation; (b) a failure model of the work material based on the Tsai Hill criterion; (c) a contact model of the mechanisms of the cutting process. A comparison with experimental measurements showed that the constitutive model leads to a reasonable prediction.

Tugrul Ozel, and Taylan Altan [7], used a methodology to determine simultaneously (a) the flow stress at high deformation rates that are encountered in the cutting zone, and (b) the friction at the chip–tool interface. This information is necessary to simulate high-speed machining using FEA based programs. A flow stress model based on process dependent parameters such as strain, strain-rate was used together with a friction model based on shear flow stress of the workpiece at the chip–tool interface. High-speed cutting experiments and process simulations were utilized to determine the unknown parameters in flow stress and friction models.

J.C. Outeiro, et al. [8], present the influence of cutting process parameters on machining performance and surface integrity generated during dry turning of Inconel 718 and austenitic stainless steel AISI 316L with coated and uncoated carbide tools. A three-dimensional Finite Element Model was also developed and the predicted results were compared with those measured.

Johnson-Cook [9] and Zerrili- Armstrong [10] models with the material constants obtained from experiments. These material models determine the workpiece material behavior in high strain, wide range of strain rates of cutting zone. The effect of the material model constants on the orthogonal cutting simulation accuracy is also investigated in the literature [11]. Tool-chip friction model is another important input parameter for the modeling of the cutting process. Many investigations are performed to study the influence of the friction modeling on metal cutting simulation results [12-14]. AISI 4140 is used extensively in most industry sectors for a wide range of applications such as Adapters, Arbors, Axle Shafts, Bolts, Crankshafts, Connection Rods, Chuck Bodies, Collets, Conveyor Pins & Rolls, Ejector Pins, Forks, and Gears.

3. Experimental Procedure

Turning is a popularly used machining process in which, a single point cutting tool removes unwanted material from the surface of a rotating workpiece. The conventional lathe machines play a major role in machining industry to enhance product quality as well as productivity.

3.1 Design of Experiments

The modeling of parameters for an AISI-4140 steel and tool material involved a large number of other factors, such as ways of holding the workpiece, the geometry of cutting tool, etc. However to facilitate the experimental data collection, only three dominant factors were considered in the planning of the experiment. The factors considered were depth of cut, feed rate, and cutting speed. The range of values of each factor was set at the three levels, namely low, medium, high, as shown in the Table 2. Based on L9 orthogonal array, a total of 9 experiments, each having a combination of different levels of factors as shown in Table 3, were carried out.

Table 2. Range of Values for Each Variable

Table 3. L9 Orthogonal Array

3.2 Experimentation

The experimentation was carried to measure the cutting forces at different set cutting parameters. Machining tests are carried out on AISI 4140 Steel specimen of 350 mm long and a 50 mm diameter is turned by using tungsten carbide tool (without coating and with Tin coating) on Kirloskar Turnmaster- 35. All are Geared conventional lathe. The tool insert used in the process is THN SNMG 08. The variable speed upto 850 rpm can be used on the machine. The workpieces are cleaned prior to the experiments by removing 0.3 mm thickness of the top surface in order to eliminate any surface defect and wobbling. The experimental setup of the dynamometer in conjunction with tool for the measurement of cutting forces in X,Y,Z directions is shown in Figure 1. From the setup, metal cutting is carried with automatic feed and depth of cut. At different setup cutting conditions such as speed, feed and depth of cut, the cutting forces are measured using dynamometer.

Figure 1. Experimental Setup of Measurement of Cutting Forces

4. Finite Element Modeling

The Deform-3D software is selected for the FEM simulation of the orthogonal cutting. The Finite Element Modeling (FEM) is considered as a very good technique to obtain the machining parameters such as cutting force, stresses, and others analysis. A typical finite element analysis in a software system requires the information procedures like, workpiece and tool is set to 30°C and simulations are carried out to reach steady state condition. The FEA is computed on Intel core i3, 2.17GHz, 4GB RAM,32 Bit computer system.

4.1 Material Modeling

Material and tool modeling by classifying (parts and properties), contact and failure laws analysis, meshing elements and boundary conditions. The effect of process parameters i.e. cutting speed, feed rate and depth of cut on the process response i.e. cutting forces, or stress, etc. can be studied. For turning of work piece, the dimensions and material data are entered in a pre-processor for creating the workpiece and tool and the data required like processing conditions of turning to analyze the simulation of forces. A simulation engine is analyzing the results by numerical calculations. Tool and workpiece material Modeling, a fully thermo-mechanically coupled implicit is considered and is shown in Figure 2. The materials and tools modeling are carried out by using Finite Element Modeling (FEM). The workpiece material dimension of 50 mm x 100 mm with machining parameters; cutting speed, v (m/min), cutting depth, d (mm), and feed rates, f (mm/tooth). Load tool holder and choose MTANR as tool holder from 'Import material from library' in workpiece material setup window.

Figure 2. FEA Preprocessing Data

4.2 Contact Analysis

In the finite element modeling, in addition to the cutting parameters, the following Johnson–Cook plasticity model was implemented. The material model employed is based on the Johnson Cook plasticity model. This constitutive relationship is commonly employed in modelling orthogonal cutting with use of FEM, since due to its being a strain rate dependent, it has a strong effect on the strain/stress relationship in the machining process. The Johnson–Cook parameters viz., d1 = 1.5, d2 = 3.44, d3 = - 2.12, d4 = 0.002, and d5 = 0.1 for AISI 4140 are taken from Deform-3D user manual. This model is suitable for problems where strain-rates vary over a large range (102- 106). The Johnson–Cook model is capable of predicting the plastic deformation as well as damage and failure.

4.3 Meshing and Boundary Conditions

The mesh generation starts by creating elements along the workpiece boundary and the cutting tool. The tool mesh with minimum element size of 0.25 mm is considered and shown in Figure 3 (a). The meshing elements and boundary conditions are used at the contact surface between the cutting tool edge and work piece. Total number of Tetrahedral elements generated are about 10,000 and 20,000 nodes for work piece and tool respectively as shown in Figure 3 (b). Enter the number of simulations as 200 steps and collect the forces data after completion of simulations on Deform-3D in the post processor group as shown in Figure 4.

Figure 3 (a). Work Piece Mesh Generation

Figure 3 (b). Tool Mesh Generation

Figure 4. FEM Simulation Result for Cutting Force

5. Results and Discussion

The cutting forces in X, Y & Z direction are determined from both simulation and experimental setup and are tabulated in Table 4 for uncoated cutting tool and Table 5 for coated cutting tool. Figures 5, 6 and 7 show the comparison of simulation and experimental results in turning used uncoated WC tool and Figures 8, 9, 10 show the comparison of simulation and experimental results in turning using TiN coated WC tool. The comparison shows that the forces using coated tool are lower than using the uncoated tool.

Table 4. Forces with uncoated Cutting Tool

Table 5. Forces with coated Cutting Tool

Figure 5. Comparison of Experimental and FEM Values of Cutting Force with Uncoated Tool

Figure 6. Comparison Experimental and FEM Values of Feed Force with Uncoated Tool

Figure 7. Comparison Experimental and FEM Values of Thrust Force with Uncoated Tool

Figure 8. Comparison Experimental and FEM Values of Thrust Force with Coated Tool

Figure 9. Comparison Experimental and FEM Values of Feed Force with Coated Tool

Figure 10. Comparision Experimental and FEM Values of Cutting Force with coated Tool

The FEM simulation results compared with experimental values are presented in Table 5. The graphs of cutting forces, feed forces and thrust forces are plotted between FEM values and experimental values respectively. Figure 8, 9 and 10 shows that, FEM values are very close to the experimental values and also they project similar trend. The FEM model is extensively useful. Figure 6 shows the feed forces of both experimental values are and FEM values are very close and their average values also very close. Both coated and uncoated cases of FEM feed force values are very close to the experimental values. Further, it is observed that with TIN coating the experimental forces are reduced in comparison with the uncoated tools. The same was reflected with simulation results in Figure 11.

Figure 11. Comparison of Resultant Forces of Simulation and Experimental Values

Conclusion

Finite Element Method using DEFORM 3D program is found to be a successful technique to perform trend analysis to estimate cutting forces in metal cutting with respect to various combinations of design variables (metal cutting speed, feed rate and depth of cut). The simulation results of cutting forces were agreeable to the experimental results and the errors are within 5%. Starting to rise rapidly, it gradually enters the mean steady state, where the cutting force began to decline, and ultimately remained stable. Increased cutting speeds (Vc) resulted in decreased cutting tool forces. The maximum forces of the workpiece during machining is reduced to a much lower value by using TiN coated WC tool.

References

[1]. M. H. El-Axir (2002). “A method of modeling residual stress distribution in turning for different material”. International Journal of Machine Tool & Manufacture, Vol. 42, pp. 1055-1063.
[2]. O. Pantal, J.-L. Bacaria, O. Dalverny, R. Rakotomalala, and S. Caperaa (2004). “2D and 3D numerical models of metal cutting with damage effects”. Comput. Methods Appl. Mech. Engg, Vol .193 .
[3]. Domenico Umbrello, (2008). “Finite element simulation of conventional and high speed machining of Ti6Al4V alloy”. Journal of Materials Processing Technology, Vol.196, pp.79–87 .
[4]. Hendri Yandra, Jaharah A. Ghani & Che Hassan Haron, (2010). “Effect of Rake angle on Stress, Strain and Temperature on the Edge of carbide cutting tool in Orthogonal Cutting using FEM Simulation”. ITB J. Eng. Sci., Vol. 42, No. 2 .
[5] Asmaa A. Kawi, (2011). “Temperature behavioure of some alloy steels in Turning Process under different operating conditions”. Al-Qadisiya Journal for Engineering Sciences, Vol. 4, No. 3 .
[6]. Mohid Mahdi, & Liangchi Zhang (2001). “A Finite element model for the orthogonal cutting of fiberreinforced composite materials”. Journal of Material Processing Technology , Vol .113 , pp. 373-377.
[7]. Tugrul Ozel, & Taylan Altan, (2000). “Determination of workpiece flow stress and friction at the chip–tool contact for high- speed cutting”, International Journal of Machine Tools & Manufacture, Vol.40, pp. 133-152.
[8]. J.C. Outeiro, J.C. Pina, R.M. Saoubi, F. Pusavec, and I.S. Jawahir, (2000). “Analysis of residual stresses induced by dry turning of difficult-to-machine materials”. CIRP Annals – Manufacturing Technology, Vol. 57, pp. 77-80 .
[9]. Baker, M. (2006). “Finite element simulation of highspeed cutting forces”. Journal of Materials Processing Technology ,Vol .176, pp.117– 126.
[10]. Johnson, G. R., & Cook, W. H. (1985). “Fracture characteristics of three metals subjected to various strains, strain rates, temperatures and pressures”. Engineering Fracture Mechanics , Vol .21, pp. 31-48.
[11]. Zerilli, F. J. & Armstrong R.W., (1987). “Dislocation mechanics- based constitutive relations for material dynamics calculations”. Journal of Applied Physics, Vol. 61 No. 5, pp.1816–1825.
[12]. Sima, M, & Ozel, T., (2010). “Modified materials constitutive models for serrated chip formation simulations and experimental validation in machining of titanium alloy Ti- 6Al-4V, International Journal of Machine Tools & Manufacture, Vol. 50, No.11, pp. 943-960.
[13]. Calamaz ,M, Coupard, D. & Girot, F. (2008). “A new Material Model for 2D numerical simulation of serrated clip formation when machining titaniumalloy Ti-6Al-4V”, International Journal of Machine Tools & Manufacturing, Vol .48, No.3-4, pp. 275-288.
[14]. Rhim, S-H. & Oh, S-I. (2006). “Predication of serrated clip formation in metal cutting process with new flow stress model for AISI1045 Steel”, Journal of Materials Processing Technology, Vol .171, No.3, pp. 417-422.