Numerical Simulation of Metal Machining Process with Eulerian and Total Lagrangian SPH

04/04/2019 ∙ by Md Rushdie Ibne Islam, et al. ∙ 0

This paper presents numerical simulations of metal machining processes with Eulerian and Total Lagrangian Smoothed Particle Hydrodynamics (SPH). Being a mesh-free method, SPH can conveniently handle large deformation and material separation. However, the Eulerian SPH (ESPH) in which the kernel functions are computed based on the current particle positions suffers from the tensile instability. The Total Lagrangian SPH (TLSPH) is free of this instability as the kernel functions are calculated in the reference configurations. In this work, the metals are modelled using the Johnson-Cook constitutive model, which can capture strain hardening and thermal softening in metals. The processing/cutting tools are modelled as rigid bodies, while the metal-tool contact forces are considered using the standard SPH interaction and the particle-particle pinball contact in ESPH and TLSPH, respectively. The two methods are employed to model several cases with impact, pressing, and cutting; the results are compared with reference experimental and numerical results. It is found that both the two SPH methods can capture the salient phenomena in metal processing, e.g. strain localisation, large deformation, and material separation. However, the TLSPH approach provides a better simulation of strain localisation and chip morphology. This work shows that the TLSPH method has the potential to model the metal machining processes efficiently without any numerical instabilities.



There are no comments yet.


page 10

page 12

page 13

page 14

page 15

page 16

page 17

page 18

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

Metal machining processes such as cutting, pressing, forming, drilling and milling are one of the most important and widely-used procedures in the manufacturing industry. The optimisation of tools, specimens and processing flows are essential for achieving higher performance in manufacturing. However, the nature of the operation makes it a challenging task for numerical modelling. During processing, the specimen/workpiece undergoes finite deformation and the behaviour of the material changes based on the accumulated strains and strain rates. The machining procedure also involves heat generation; hence, the thermal softening behaviour needs to be considered. This entire process becomes even more complicated with the possibility of strain localisation, cracking, and material separation. Experimental and empirical observations (el2001effect, ; sievert2003simulation, ; tonshoff2005chip, ) are used to understand the metal cutting. Simple analytical models (merchant1945mechanics, ; lee1951theory, ) are also used. However, the cost and complexity of the experiment limit its practical application, whereas the analytical models are less reliable due to simplifications. Therefore, the numerical models are used as an alternative to model the metal machining process.

As the most widely-used numerical method, the finite element method (FEM) was employed to model orthogonal metal cutting in oh2004finite . However, the FEM has inherent difficulties in modelling problems involving finite deformations and material separations due to the utilisation of mesh. Mesh distortion and discontinuities require a special technique to handle. Several studies employed the Arbitrary Lagrangian-Eulerian (ALE) method (movahhedy2000simulation, ; pantale20042d, ) to treat the excessive element distortion and material separation, in which the Johnson-Cook material model (johnson1983constitutive, ) is used to account for the influence of plastic strain, strain rate and thermal softening. In ALE, the issues of element distortion, entanglement, and material separation are handled through continuous re-meshing. However, the re-meshing is not only computationally intensive but also difficult to perform in three-dimensional complex cases. Moreover, it is well-known that FEM suffers from mesh-dependency in modelling strain localisation, a dominant mechanism in metal processing. For instance, the influence of mesh size and orientation in metal cutting are investigated in (hortig2007simulation, ). It is found that the mesh influences the formation of shear bands, chips and the values of the computed cutting force. It is, therefore, challenging to obtain mesh-objective results with the ALE approach.

Lagrangian particle-based methods are naturally capable of capturing large material deformation due to the absence of meshes. The issues related to mesh, e.g. distortion, entanglement are not encountered. Therefore, some of the meshless methods are employed to model metal processing. For instance, the Material Point Method (MPM) is used to simulate the metal pressing/cutting process in ambati2012application . Another popular meshless method for metal cutting modelling is Smoothed Particle Hydrodynamics (SPH) (limido2007sph, ; villumsen2008simulation, ; limido2011metal, ; ruttimann2012simulation, ; zahedi2012application, ; eberhard2013simulation, ; olleak2015influence, ; olleak2015prediction, ). SPH is a truly meshless method as it does not require any background mesh. It was first developed for the astronomical problems (gingold1977smoothed, ; monaghan1985refined, ) and later gradually used for the simulation of fluid flows and solid deformations (monaghan1994simulating, ; libersky1993high, ). In SPH, the computational domain is discretised into Lagrangian particles carrying field variables and moving with the material. The meshless and Lagrangian nature of SPH makes it a perfect candidate for metal processing modelling. It is shown in previous works that shear bands and formation of chips are captured naturally without any requirement of special treatment such as nodal enrichment or re-meshing. SPH based simulations in LS-DYNA are performed (villumsen2008simulation, ) to predict the cutting forces, plastic strains and cutting planes with reasonable accuracy. A coupled FEM and SPH framework is also developed in (zahedi2012application, ).

In conventional SPH, particle positions are updated after each computational step. The kernel functions are computed based on the updated particle positions; therefore, the traditional SPH (ESPH) kernel function is termed as Eulerian kernel function in the sense that particles can move in and out of its influence domain. Although widely used, the Eulerian kernel function leads to the well-known tensile instability swegle1995smoothed , i.e. a local clustering of particles forming unphysical numerical fracture. A popular correction is the artificial pressure/stress (monaghan2000sph, ; gray2001sph, ), which can effectively alleviate the tensile instability. However, this approach requires a proper tuning of the parameter to avoid any unphysical effect in the simulation. This tuning process is case-dependent thus highly undesirable. Studies reveal that tensile instability is associated with the use of Eulerian kernel function belytschko2000unified . It is found that this instability can be avoided if the reference particle position is used to calculate the kernel functions (belytschko2002stability, ; rabczuk2004stable, ). The reference configuration-based kernel function is called Lagrangian kernel, and the corresponding SPH formulation is termed Total Lagrangian SPH (TLSPH). The TLSPH is entirely free of tensile instability (vignjevic2006sph, ; belytschko2000unified, ; bonet2002alternative, ); therefore, it is an ideal option for metal processing modelling.

The present work aims to simulate the metal pressing/cutting processes in the conventional/ Eulerian SPH (ESPH) and TLSPH frameworks and provide a comparison between these two methods. In our study, the metal materials are modelled using the Johnson-Cook model. The processing tools are considered as rigid bodies in contact with the specimens. In the ESPH and TLSPH, the tool-metal contact forces are modelled using the standard SPH interaction and the particle-particle contact algorithm, respectively. The presented methods are first validated using a Taylor impact case. Then metal pressing and cutting problems are modelled. The results from the two formulations are compared. Numerical results from the literature are also considered as a reference to evaluate the simulations.

The remaining part of the paper is organised as follows. The formulations of the ESPH and TLSPH are introduced in section 2. The validation of the two SPH frameworks and the numerical simulations of the metal pressing and cutting are presented in section 3. Some conclusions are drawn in section 4.

2 Computational Frameworks

In this section, the ESPH and TLSPH are introduced. In SPH, the computational domain is represented by Lagrangian particles (Libersky93, ; liu2003smoothed, ; liu2006restoring, ; liu2010smoothed, ). The particles interact with each other through a kernel function. In ESPH, the kernel function is calculated based on the current configuration, which is updated in every computational step. On the other hand, the Lagrangian kernel, which is based on the undeformed reference configuration, is used in the TLSPH.

The governing equations in the current liu2010smoothed and reference de2013total configurations are

Mass conservation:


Momentum conservation:


where is the density,

is the velocity vector in the current configuration

. The subscript indicates the values are computed at the reference configuration . The current and reference descriptions are related through a mapping as . is the Jacobian computed based on the deformation gradient matrix as .

is the Cauchy stress tensor,

the first Piola Kirchhoff stress tensor. and denote the divergence/gradient operator, and the outer product between two vectors, respectively.

The deformation gradient is defined as


where is the displacement calculated as . represents the identity tensor. Similarly, the rate of deformation gradient is


where is the velocity.

2.1 Discretised forms in Eulerian SPH (ESPH)

In ESPH, any field variable and its derivative are approximated based on the current configuration as follows


where is the field variable value at the -th particle, is the vector from particle to particle . is the kernel function defined in the current configuration. represents the volume of the -th particle in the current configuration. The following cubic B-spline function is used as kernel function in this work.


where is the normalisation coefficient taken as in 2D and in 3D, is the smoothing length and is the normalised distance associated with the particle pair.

The governing Eqs. (1) and (3) in the current configuration can be discretised chakraborty2013pseudo as,


where is the kernel gradient evaluated at the current position of particle . is the artificial viscosity used to stabilise the computation in the presence of shock in field variables. The Monaghan type artificial viscosity MonaghanGingold83 with controlling parameters , is used.




with the relative velocity, the coefficient used to prevent singularity when two particles are too close. is the sound speed in the medium with Young’s modulus. A bar above a quantity indicates that it is the averaged value over the two particles.

In the ESPH, the tensile instability is common in the simulations of solids. The artificial pressure monaghan2000sph ; gray2001sph is used to treat the tensile instability. The form of artificial pressure writes


where is the pressure. is a case-dependent tuning parameter in the magnitude of 0.1. It is zero if and . is a constant taken as , is the initial particle spacing.

2.2 Discretised forms in Total Lagrangian SPH (TLSPH)

In the TLSPH description, a function and its derivative are approximated in the reference configuration as


where is the kernel gradient, and represents the distance vector and the volume of the -th particle in the reference configuration.

In the TLSPH formulation, the mass conservation equation is trivially satisfied; therefore, it is not needed to be solved numerically. The discretised form of the momentum conservation Eq. (4) writes


where is the first Piola Korchhoff stress. is the artificial viscosity. The TLSPH is free from any tensile instability due to the use of reference configuration in the computation. Therefore, there is no need of any numerical stabilization term.

Additionally, the discretised forms of the deformation gradient and its rate are


2.3 Constitutive equations

The Cauchy stress tensor can be expressed in terms of deviatoric stress and hydrostatic pressure as . The deviatoric stress rate is calculated from the material frame invariant Jaumann stress rate as


where is the shear modulus, is the strain rate tensor and is the spin tensor. is the trace of the strain rate tensor . The strain rate and spin tensors are calculated from the velocity gradient tensor as follows


In ESPH, the velocity gradient tensor is obtained as


whereas in TLSPH, the velocity gradient tensor is computed as .

In this work, the hydrostatic pressure is obtained from a linear equation of state


where bulk modulus with the poison’s ratio.

2.3.1 Viscoplastic material model

Metals are modelled using the Johnson-Cook johnson1983constitutive (J-C) model, which can consider salient metal behaviours such as plastic hardening, rate dependency, and thermal softening. In the J-C model, the von Mises yield stress is calculated as


where is the initial yield stress of the material; , , and are the hardening parameters. is the accumulated effective plastic strain. denotes the dimensionless plastic strain rate defined as


where and are the current effective plastic strain rate and the reference strain rate. The homologous temperature is defined as


where , and are the current temperature, the room temperature, and the melting temperature of the metal, respectively. It is assumed that the temperature increases due to some fraction of the plastic work being converted into heat.


where , and are the accumulated plastic work density, the specific heat density at constant pressure, and the Taylor-Quinney empirical constant. denotes the percentage of plastic work under the adiabatic condition that is dissipated as heat.

The von Mises yield criterion is used to account for metal plasticity, where is the yield stress obtained using Eq. (25), and is the second invariant of the deviatoric stress tensor. The Wilkins criterion is used for return mapping when the trial elastic stress state is beyond the yield surface, with and the corrected deviatoric stress tensor. Finally, the following equations are used to compute the increment of plastic strain, the increment of effective plastic strain, and the accumulated plastic work density


2.4 Contact Force

The interaction between different bodies in ESPH can be modelled naturally with the conventional SPH particle interaction. However, the contact force cannot be modelled using the SPH particle interaction in TLSPH due to the use of the reference configuration. Consequently, an explicit algorithm for contact force is essential for the TLSPH for modelling multibody interaction. The pin-ball contact campbell2000contact is applied in this work. Here, the contact algorithm is discussed briefly. For more details, the readers may refer to islam2019total .

In the pin-ball contact, the particles are assumed to have a predefined finite influence zone, as shown in Figure 1. The influence zone is represented by a sphere and a circle in 3D and 2D, respectively. The radius of the circle or the sphere is , where and are the smoothing length and a constant factor, respectively. For one particle, the contact force is computed if another particle from a different body comes within its influence radius .

Figure 1: Penetration between particles for contact force

The contact force is a function of material parameters, penetration depth, and rate of penetration. First, the magnitude of penetration of the particle pair is computed as


where and are the influence radius of -th and -th particle, respectively. indicates material penetration, which results in the contact force belytschko1993splitting




where , , and are the time step, the rate of penetration, and the scaling factor generally chosen through numerical experiment. The momentum Eq. (17) in the TLSPH is modified by including the contact force as


Note that the contact force is only employed in simulations using TLSPH.

3 Numerical simulation

In this section, the ESPH and TLSPH methods are first validated using a Taylor impact test batra2008ssph . The numerical simulations of metal pressing and cutting are then performed. The predictor-corrector time integration is used to integrate the discretised governing equations. The Courant–Friedrichs–Lewy (CFL) condition is applied to determine the stable time step. In all the simulations, the tool-metal contact is assumed to be frictionless. In all simulations with the ESPH, the artificial pressure coefficient is taken as 0.3. In TLSPH, the reference configuration is updated to the current configuration when the deviation between the current and reference frames is considerable leroch2016smooth . The particle volumes and the material coordinates are updated as and respectively. The updates preserve the stress tensors. The updates of the reference frames are essential for problems involving a large plastic flow of materials to capture the changes in particle interaction leroch2016smooth .

3.1 Taylor impact test

A low strength 4340 steel cylinder of 37.97 mm length and 7.595 mm diameter strikes a rigid wall with an initial velocity of 181.0 m/s batra2008ssph as shown in Figure 2. The material and computational parameters are given in Table 1 and 2, respectively.

Figure 2: Setup for the taylor impact test
Mechanical properties J-C strength parameters J-C thermal parameters
kg/m = 792 MPa = 293 K
GPa = 510 MPa = 1293 K
GPa = 0.26 = 1.03
= 0.014 J/(kgK)
Table 1: Johnson-Cook (J-C) Material parameters of low strength 4340 steel
Inter-particle Smoothing Time step Artificial viscosity
spacing () length () () parameters ()
0.5 mm 0.75 mm s (1.5, 1.5)
Table 2: Computational parameters for Taylor impact test

The numerical results of the final length of the cylinder after impact and the final diameter at the mushrooming end are compared with the experimental results house1995estimation and the numerical results obtained with LS-DYNA batra2008ssph . The final diameters of the mushrooming end are found to be 9.80 mm and 9.83 mm with the ESPH and TLSPH, respectively. The experimental value is 9.5 mm, and the computed value from LS-DYNA is 10.6 mm. The predictions of the final diameter of the cylinder at the mushrooming end are very close to the experimental observation. The current predictions deviate from the experimental observation by 3.16 and 3.47 with the ESPH and TLSPH, respectively. However, the deviation of the computed diameter by LS-DYNA is 11.58 to the experimental result. The final lengths of the deformed cylinder are also calculated and are found to be 34.19 mm and 34.18 mm with the ESPH and TLSPH, respectively. The experimental and computed values with LS-DYNA are found to be 34.6 mm and 34.4 mm. The present predictions are close to experimental and numerical observation. The differences between our results and the experimental results are 1.18% and 1.21% for the ESPH and TLSPH, respectively.

The change in the length of the cylinder is also compared with the LS-DYNA results in Figure 3. The present predictions are in good agreement with the LS-DYNA results. The contours of the accumulated effective plastic strain and temperature are shown in Figure 4 and 5 respectively. It can be readily observed that the distributions of the plastic strain and temperature are captured quite well and are similar to the LS-DYNA results. However, the contour plot from the TLSPH provides a better agreement in the distribution of plastic strain and temperature with the LS-DYNA. Furthermore, the particle distributions in the TLSPH results are more regular, whereas, in the results obtained from the ESPH, the irregular particle distribution caused by tensile instability can be observed even with the artificial pressure.

Figure 3: Comparison of the change in length of the cylinder over time
Figure 4: Comparison of accumulated effective plastic strain with the LS-DYNA results from batra2008ssph
Figure 5: Comparison of temperature (K) contour with the LS-DYNA results from batra2008ssph

3.2 Metal pressing

In this section, the pressing of a metal block of AISI 4340 steel by a rigid tool is modelled, as shown in Figure 6. The tool is pressed from the free right side, while the left side of the metal block is fixed. The workpiece is of 1 mm 0.5 mm 0.1 mm, and the pressing tool is of 0.1 mm 0.7 mm 0.1 mm. The pressing tool moves with a constant velocity of 200 m/s. A three-dimensional simulation is performed, using 184,624 particles for the workpiece and 25,199 particles for the pressing tool. The response of the workpiece is modelled by the Johnson-Cook material model and the pressing tool is modelled as a rigid body. The constitutive parameters are listed in Table 3. The other relevant computational parameters are given in Table 4.

Figure 6: Setup for the metal pressing process
Mechanical properties J-C strength parameters J-C thermal parameters
kg/m = 792 MPa = 298 K
GPa = 510 MPa = 1573 K
= 0.26 = 1.03
= 0.014 J/(kgK)
Table 3: Johnson-Cook (J-C) Material parameters for AISI4340 steel workpiece
Inter-particle Smoothing Time step Artificial viscosity
spacing () length () () parameters ()
0.00667 mm 0.01 mm s (0.5, 0.5)
Table 4: Computational parameters for metal pressing

The numerical results from the ESPH and TLSPH are compared with those from the Material Point Method (MPM) and Finite Element Method (FEM) simulations ambati2012application . The final shape and distribution of equivalent plastic strain are shown in Figure 7, with results from all the four approaches. It is observed that the results from the four simulations have similar final shape and distribution of equivalent plastic strain, indicating the both the ESPH and TLSPH can correctly capture the large deformation and strain localisation in the pressing. Compared with MPM, the final shapes obtained from SPH simulations are closer to the FEM results, as the MPM results show large deformation along the secondary group of shear bands, leading to somewhat different final shape. Also, the equivalent plastic strains in the MPM results are diffusive and lack distinct shear bands. On the other hand, the results from the two SPH methods show clear shear bands, like those from FEM. However, the ESPH results show slightly more diffused distribution of equivalent plastic strain. This is due to the application of the updated particle position in the computation.

(a) FEM
(b) MPM
(c) ESPH
Figure 7: Distribution of plastic strain accumulation for metal pressing with FEM ambati2012application , MPM ambati2012application , ESPH and TLSPH

The accumulated equivalent plastic strain across the metal block is plotted in Figure 8. The two methods give very similar results except for the most left part. The figure shows a consistent trend that the right part of the metal block is subjected to more plastic deformation. The increase in temperature is shown in Figure 9. It is observed that the area and magnitude of temperature change are consistent with the equivalent plastic strain. This increase in temperature can further reduce the yield stress in shear bands, making the sling of material in these regions more easily.

Figure 8: The accumulated equivalent plastic strain across the block.
(a) ESPH
Figure 9: Distribution of temperature for metal pressing with ESPH and TLSPH

To analyse the effect of initial particle spacing, four simulations with different resolution ( mm, 0.0133 mm, 0.01 mm, 0.0067 mm) are performed. In the four simulations, the ratio between the smoothing length and initial particle spacing is kept constant as 1.5. The results are compared with the results obtained from the particle spacing 0.0067 mm in Figure 10. It can be observed that the plastic strain bands become more distinct as the particle resolution is refined. Moreover, results from TLSPH show more regular particle distribution than those from ESPH.

(a) ESPH
Figure 10: Distribution of plastic strain accumulation for metal pressing with the ESPH and TLSPH framework at different particle resolution ()

3.3 Metal cutting

Formation of shear bands and accumulation of plastic shear strain play an essential role in metal cutting processes. As the cut materials often undergo considerable plastic deformation, high strain localisation, and material separation, mesh-based methods such as FEM have difficulties because of mesh distortion and entanglement. The discontinuities in metal cutting processes are also problematic for mesh-based methods. As SPH is meshless, it is free from mesh deficiencies and naturally capable of simulating large deformation. Moreover, material separations can be simulated naturally in SPH.

Figure 11: Setup for the metal cutting process

The setup for the orthogonal metal cutting process is shown in Figure 11. The workpiece (3 mm 1 mm) is made of AISI 4340 steel. The JC material model is used to simulate the material behaviour. The JC parameters are shown in Table 3. In the metal cutting process, a rigid cutting tool (0.5 mm 1 mm) moves with a constant velocity of 50 m/s ambati2012application . In the present work, it is further assumed that there is no friction between the workpiece and the cutting tool. Other computational parameters are listed in Table 5. The pinball contact model is used for the computation of the contact forces in the TLSPH.

Inter-particle Smoothing Time-step Artificial viscosity
spacing () length () () parameters ()
0.01 mm 0.018 mm s (0.5, 0.5)
Table 5: Computational parameters for metal cutting

3.3.1 Results with feeding length of 0.3 mm

The plastic strain and temperature distribution are important output variables in the metal cutting process. The results of accumulated equivalent plastic strain from the ESPH and TLSPH are shown in Figure 12 and 13, respectively. It can be observed that both approaches produce continuous chip morphologies in the metal cutting process. As the tool cuts further into the workpiece, plastic strain accumulates in front of the toe of the tool, and then the strain localisation propagates towards the surface in an inclination approximately equal to 45. This numerical observation is inconsistent with the plastic theory. Once a shear band is fully developed, the cut material above the newly-formed shear band slides along it, leading to a new chip. With this mechanism, results from both the ESPH and TLSPH show a number of chips separated by shear bands. The increase of temperature is also shown in Figure 14 and 15. The increase in temperature is consistent with the distribution of equivalent plastic strain.

Figure 12: Generation of chips and accumulation of effective plastic strain in metal cutting process using ESPH with 0.3 mm feeding length
Figure 13: Generation of chips and accumulation of effective plastic strain in metal cutting process using TLSPH with 0.3 mm feeding length
Figure 14: Distribution of temperature in metal cutting process using ESPH with 0.3 mm feeding length
Figure 15: Distribution of temperature in metal cutting process using TLSPH with 0.3 mm feeding length

The formation of chips in the metal cutting process is compared with results from MPM and FEM in Figure 16. Numerical results of chip formation and distribution of shear bands from ESPH and TLSPH are in good agreement with the MPM and FEM results. However, the formation of chips and the deformation of the cut stripe captured using TLSPH are closer to the MPM and FEM results. The reason may lie in the formulation of TLSPH in which the number of instabilities is less than the ESPH. Moreover, the bending of the cut stripe in the TLSPH is more than the ESPH simulation, which is closer to the results from MPM and FEM. Furthermore, the shear bands are more distinct in results from the TLSPH. More importantly, the simulations using TLSPH does not need to tune the coefficient of artificial pressure like in ESPH.

Figure 16: Comparison of metal cutting process with ESPH, TLSPH, MPM ambati2012application , FEM ambati2012application

3.3.2 Effect of feeding length on shear band

The shear bands and the accumulation of plastic strains play an important role in the metal cutting process. The increase of the feeding length increases the cutting force in the simulation. This, in turn, increases the amount of plastic deformation and heat dissipation in the workpiece. The depth of feeding is increased to 0.5 mm from 0.3 mm used in the previous metal cutting simulation. This influences the metal cutting process and changes the chip morphology. The deformation of the workpiece with 0.5 mm feeding length are shown in Figure 17 and 18 by the ESPH and TLSPH respectively. It is observed that the increase in the cutting forces leads to the increment of the size of the chips, i.e. the number of shear bands decreases. This is due to the increase of the localisation of material plastic strains and leads to larger chips between shear bands as shown in Figure 17 and 18. The computed results are compared with the results from MPM ambati2012application . It can be observed in Figure 19 that the present results are in good agreement with the results from ambati2012application .

Figure 17: Generation of chips and accumulation of effective plastic strain in metal cutting process using ESPH with 0.5 mm feeding length
Figure 18: Generation of chips and accumulation of effective plastic strain in metal cutting process using TLSPH with 0.5 mm feeding length
Figure 19: Comparison of metal cutting process with ESPH, TLSPH, MPMambati2012application

4 Conclusion

The numerical predictions of metal machining process such as metal pressing and cutting are a complex problem. These processes involve localisation of plastic strain, shear band formation, material separation and temperature increment leading to material softening. The mesh-based methods such as FEM are not a suitable candidate due to mesh distortion, entanglement etc. In this paper, the ESPH and TLSPH frameworks are used to model the metal machining process. Due to the particle nature, SPH is inherently capable of capturing the large plastic deformation of materials. However, the ESPH suffers from tensile instability leading to numerical fracture of material. The TLSPH is free from this due to the computation of kernel functions in the reference configuration. The ESPH and TLSPH frameworks are first verified using the Taylor impact test. The computed final length, diameter and the change in length over time are in good agreement with the experimental and numerical observations.

The ESPH and TLSPH are used to simulate the metal pressing and cutting process. The accumulated effective plastic strains and the distribution of temperature are compared with the numerical result available in the literature. The presented results are consistent with the observations from the literature. It is further observed that the TLSPH provides a more realistic chip morphology, deformation of the workpiece and generation of shear bands. The change in feeding length influences the material deformation. With the increase in the feeding length, the cutting force increases. This leads to a reduction in the number of shear bands and an increase in the width between two adjacent shear bands. In the current work, the effect of friction in the metal cutting process is neglected. However, the influence of frictional force and material failure may have impacts on metal machine processing. In future, more advanced contact algorithm will be considered for the frictional contact.


This work has received funding from the European Unions Horizon 2020 research and innovations program under grant agreement No. 778627 and from FFG the Austrian Funding Agency under the project No. 865963.


  • (1) T. El-Wardany, M. Elbestawi, Effect of material models on the accuracy of highspeed machining simulation, Scientific Fundamentals of High-Speed Cutting, Hanser (2001) 77–91.
  • (2) R. Sievert, A. Hamann, H. Noack, P. Löwe, K. Singh, G. Künecke, R. Clos, U. Schreppel, P. Veit, E. Uhlmann, et al., Simulation of chip formation with damage during high-speed cutting, Tech. Mech 23 (2003) 216–233.
  • (3) H. Tönshoff, B. Denkena, R. B. Amor, A. Ostendorf, J. Stein, C. Hollmann, A. Kuhlmann, Chip formation and temperature development at high cutting speeds, Hochgeschwindigkeitsspanen, Wiley-vch (2005) 1–40.
  • (4) M. E. Merchant, Mechanics of the metal cutting process. i. orthogonal cutting and a type 2 chip, Journal of applied physics 16 (5) (1945) 267–275.
  • (5) E. Lee, The theory of plasticity applied to a problem of machining, ASME J. Appl. Mech. 18 (1951) 405.
  • (6) J.-D. Oh, J. C. Aurich, A finite-element-analysis of orthogonal metal cutting processes, in: AIP Conference Proceedings, Vol. 712, AIP, 2004, pp. 1390–1395.
  • (7) M. Movahhedy, M. Gadala, Y. Altintas, Simulation of the orthogonal metal cutting process using an arbitrary lagrangian–eulerian finite-element method, Journal of materials processing technology 103 (2) (2000) 267–275.
  • (8) O. Pantalé, J.-L. Bacaria, O. Dalverny, R. Rakotomalala, S. Caperaa, 2d and 3d numerical models of metal cutting with damage effects, Computer methods in applied mechanics and engineering 193 (39-41) (2004) 4383–4399.
  • (9) G. R. Johnson, W. H. Cook, A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures, in: Proceedings of the 7th International Symposium on Ballistics, Vol. 21, The Hague, The Netherlands, 1983, pp. 541–547.
  • (10) C. Hortig, B. Svendsen, Simulation of chip formation during high-speed cutting, Journal of Materials Processing Technology 186 (1-3) (2007) 66–76.
  • (11) R. Ambati, X. Pan, H. Yuan, X. Zhang, Application of material point methods for cutting process simulations, Computational Materials Science 57 (2012) 102–110.
  • (12) J. Limido, C. Espinosa, M. Salaün, J.-L. Lacome, Sph method applied to high speed cutting modelling, International journal of mechanical sciences 49 (7) (2007) 898–908.
  • (13) M. F. Villumsen, T. G. Fauerholdt, Simulation of metal cutting using smooth particle hydrodynamics, LS-DYNA Anwenderforum, C-III (2008) 17.
  • (14) J. Limido, C. Espinosa, M. Salaun, C. Mabru, R. Chieragatti, J.-L. Lacome, Metal cutting modelling sph approach, International journal of machining and machinability of materials 9 (3-4) (2011) 177–196.
  • (15) N. Rüttimann, Simulation of metal cutting processes using meshfree methods, Ph.D. thesis, ETH Zurich (2012).
  • (16) A. Zahedi, S. Li, A. Roy, V. Babitsky, V. V. Silberschmidt, Application of smooth-particle hydrodynamics in metal machining, in: Journal of Physics: Conference Series, Vol. 382, IOP Publishing, 2012, p. 012017.
  • (17) P. Eberhard, T. Gaugele, Simulation of cutting processes using mesh-free lagrangian particle methods, Computational Mechanics 51 (3) (2013) 261–278.
  • (18) A. A. Olleak, M. N. Nasr, H. A. El-Hofy, The influence of johnson-cook parameters on sph modeling of orthogonal cutting of aisi 316l.
  • (19) A. A. Olleak, H. A. El-Hofy, Prediction of cutting forces in high speed machining of ti6al4v using sph method, in: ASME 2015 International Manufacturing Science and Engineering Conference, American Society of Mechanical Engineers, 2015, pp. V001T02A018–V001T02A018.
  • (20) R. A. Gingold, J. J. Monaghan, Smoothed particle hydrodynamics: theory and application to non-spherical stars, Monthly Notices of the Royal Astronomical Society 181 (3) (1977) 375–389.
  • (21) J. J. Monaghan, J. C. Lattanzio, A refined particle method for astrophysical problems, Astronomy and astrophysics 149 (1985) 135–143.
  • (22) J. J. Monaghan, Simulating free surface flows with sph, Journal of computational physics 110 (2) (1994) 399–406.
  • (23) L. D. Libersky, A. G. Petschek, T. C. Carney, J. R. Hipp, F. A. Allahdadi, High strain lagrangian hydrodynamics: a three-dimensional sph code for dynamic material response, Journal of computational physics 109 (1) (1993) 67–75.
  • (24) J. Swegle, D. Hicks, S. Attaway, Smoothed particle hydrodynamics stability analysis, Journal of computational physics 116 (1) (1995) 123–134.
  • (25) J. J. Monaghan, SPH without a tensile instability, Journal of Computational Physics 159 (2) (2000) 290–311.
  • (26) J. P. Gray, J. J. Monaghan, R. P. Swift, SPH elastic dynamics, Computer Methods in Applied Mechanics and Engineering 190 (49) (2001) 6641–6662.
  • (27) T. Belytschko, Y. Guo, W. Kam Liu, S. Ping Xiao, A unified stability analysis of meshless particle methods, International Journal for Numerical Methods in Engineering 48 (9) (2000) 1359–1400.
  • (28) T. Belytschko, S. Xiao, Stability analysis of particle methods with corrected derivatives, Computers & Mathematics with Applications 43 (3-5) (2002) 329–350.
  • (29) T. Rabczuk, T. Belytschko, S. Xiao, Stable particle methods based on lagrangian kernels, Computer methods in applied mechanics and engineering 193 (12-14) (2004) 1035–1063.
  • (30) R. Vignjevic, J. R. Reveles, J. Campbell, Sph in a total lagrangian formalism, CMC-Tech Science Press- 4 (3) (2006) 181.
  • (31) J. Bonet, S. Kulasegaram, Alternative total lagrangian formulations for corrected smooth particle hydrodynamics (csph) methods in large strain dynamic problems, Revue Européenne des Éléments Finis 11 (7-8) (2002) 893–912.
  • (32) L. D. Libersky, A. G. Petschek, T. C. Carney, J. R. Hipp, F. A. Allahdadi, High strain lagrangian hydrodynamics: A three-dimensional SPH code for dynamic material response, Journal of Computational Physics 109 (1993) 67–75.
  • (33) G.-R. Liu, M. B. Liu, Smoothed particle hydrodynamics: a meshfree particle method, World Scientific, 2003.
  • (34) M. B. Liu, G. R. Liu, Restoring particle consistency in smoothed particle hydrodynamics, Applied Numerical Mathematics 56 (1) (2006) 19–36.
  • (35) M. B. Liu, G. R. Liu, Smoothed particle hydrodynamics SPH: an overview and recent developments, Archives of Computational Methods in Engineering 17 (1) (2010) 25–76.
  • (36) T. De Vuyst, R. Vignjevic, Total lagrangian sph modelling of necking and fracture in electromagnetically driven rings, International Journal of Fracture 180 (1) (2013) 53–70.
  • (37) S. Chakraborty, A. Shaw, A pseudo-spring based fracture model for SPH simulation of impact dynamics, International Journal of Impact Engineering 58 (2013) 84–95.
  • (38) J. J. Monaghan, R. A. Gingold, Shock simulation by the particle method SPH, Journal of Computational Physics 52 (2) (1983) 374–389.
  • (39) J. Campbell, R. Vignjevic, L. Libersky, A contact algorithm for smoothed particle hydrodynamics, Computer methods in applied mechanics and engineering 184 (1) (2000) 49–65.
  • (40) M. R. I. Islam, C. Peng, A total lagrangian sph method for modelling damage and failure in solids, arXiv preprint arXiv:1901.08342.
  • (41) T. Belytschko, I. Yeh, The splitting pinball method for contact-impact problems, Computer methods in applied mechanics and engineering 105 (3) (1993) 375–393.
  • (42) R. C. Batra, G. M. Zhang, SSPH basis functions for meshless methods, and comparison of solutions with strong and weak formulations, Computational Mechanics 41 (4) (2008) 527–545.
  • (43) S. Leroch, M. Varga, S. Eder, A. Vernes, M. R. Ripoll, G. Ganzenmüller, Smooth particle hydrodynamics simulation of damage induced by a spherical indenter scratching a viscoplastic material, International Journal of Solids and Structures 81 (2016) 188–202.
  • (44)

    J. W. House, J. C. Lewis, P. P. Gillis, L. Wilson, Estimation of flow stress under high rate plastic deformation, International journal of impact engineering 16 (2) (1995) 189–200.