Next Article in Journal
Effect of Reinforcement Ratio and Bond Characteristic on Flexural Behavior of Carbon Textile-Reinforced Concrete Panels
Previous Article in Journal
Advances in Two-Dimensional Magnetic Semiconductors via Substitutional Doping of Transition Metal Dichalcogenides
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermomechanical Simulation of Orthogonal Metal Cutting with PFEM and SPH Using a Temperature-Dependent Friction Coefficient: A Comparative Study

by
Juan Manuel Rodríguez Prieto
1,2,
Simon Larsson
2,* and
Mohamadreza Afrasiabi
3
1
Escuela de Ciencias Aplicadas e Ingeniería, Universidad EAFIT, Cra 49 n 7–sur–50, Medellín 050022, Colombia
2
Department of Engineering Sciences and Mathematics, Solid Mechanics, Luleå University of Technology, 971 87 Luleå, Sweden
3
Data-Driven & Computational Manufacturing Group, Inspire AG, 8005 Zürich, Switzerland
*
Author to whom correspondence should be addressed.
Materials 2023, 16(10), 3702; https://doi.org/10.3390/ma16103702
Submission received: 24 February 2023 / Revised: 16 March 2023 / Accepted: 12 April 2023 / Published: 12 May 2023

Abstract

:
In this work, we apply the Particle Finite Element Method (PFEM) and Smoothed Particle Hydrodynamics (SPH) to simulate the orthogonal cutting chip formation of two workpiece materials, i.e., AISI 1045 steel and Ti6Al4V titanium alloy. A modified Johnson–Cook constitutive model is used to model the plastic behavior of the two workpiece materials. No damage or strain softening is included in the model. The friction between the workpiece and the tool is modeled following Coulomb’s law with a temperature-dependent coefficient. The accuracy of PFEM and SPH in predicting thermomechanical loads at various cutting speeds and depths against the experimental data are compared. The results show that both numerical methods can predict the rake face temperature of AISI 1045 with errors less than 34%. For Ti6Al4V, however, the temperature prediction errors are significantly higher than those of the steel alloy. Errors in force prediction were in the range of 10% to 76% for both methods, which compare very well with those reported in the literature. This investigation infers that the Ti6Al4V behavior under machining conditions is difficult to model on the cutting scale irrespective of the choice of numerical method.

1. Introduction

Orthogonal cutting simulation remains a challenging task due to multiple computational and material complexities. Firstly, it involves the solution of differential equations that represent the physics of the problem in a domain that changes throughout the simulation, generating very large strains and deformations. Secondly, it needs a constitutive model that represents the behavior of the workpiece material in a wide range of strains, strain rates, and temperatures. Finally, (over) simplified friction models such as Coulomb’s law, which are mostly used in the numerical modeling framework, do not adequately represent the physics at the chip–tool interface, therefore, new friction models dependent on temperature and relative velocity are needed. The first challenge has been partially resolved with remeshing in finite element models or using completely different discretization approaches such as mesh-free or particle-based methods [1,2,3,4,5,6,7]. The second challenge is in continuous development because as the new materials appear, new constitutive models need to be developed; even so, physically based and dislocation density models have been developed that improve the prediction of the cutting forces in comparison with the Johnson–Cook constitutive model [5,8,9,10]. The last challenge appears to be less studied than the others in the numerical simulation of cutting, even though the use of better models for describing the friction behavior in metal cutting has been repeatedly concluded by many researchers in the field [11,12,13,14]. The contribution of this paper is thus directed mainly towards the third challenge in orthogonal cutting modeling. A sample of the expected results in the orthogonal cutting simulations is shown in Figure 1.
A coupled Eulerian–Lagrangian FE model as applied to predict the cutting force and chip morphology of DE 718 using a Johnson–Cook material model and a temperature dependent friction coefficient in [15], finding better predictions of cutting forces and chip shape compared to simulations with constant coefficient of friction. At the same time, a Lagrangian finite element model of metal cutting considering eight different friction coefficient models depending on the relative velocity v S , normal contact force F, and temperature T was developed using ABAQUS/Standard in [16]. The goal was to study the influence of each friction model on the temperature distribution in the workpiece. The research found a big difference in the rising cutting temperature and in the predicted chip shapes. Moreover, a recent review of the state of the friction in metal cutting identifies that future friction models must consider thermomechanical loads [17].
In addition to the well-known Finite Element Method (FEM), particle-based methods are another spatial discretization approach frequently used for metal cutting simulations. Among the recently developed particle methods are the Particle Finite Element (PFEM) and the Smoothed Particle Finite Element (SPH) methods. Gingold et al. [18] and Lucy [19], in 1977, introduced SPH into astrophysics applications. Since the original formulation, this method has been applied to various problems, including impact [20], shock-wave simulations [21], free-surface problems [22], multi-phase flows including thermal phenomena [23], laser-based manufacturing problems, including the modeling of material removal in [24,25], pulsed-laser ablation of aluminum in [26], a ball charge and its interaction with the mill structure [27], simulation of hydraulic jump using periodic open boundaries [28], modeling dam break evolution over a wet bed [29], prediction of wear in a dumper truck body [30], and the collapse of an axisymmetric granular column of dry potassium chloride (KCl) in granular form [31].
The PFEM was initially utilized for fluid dynamics and fluid—structure interaction [32,33] problems where nonlinearities and evolving domains are present. The PFEM combines the advantages of Lagrangian particle methods with the accuracy and robustness of mesh-based methods. The PFEM discretizes the physical domain with a mesh, on which the conservation equations are solved using a Lagrangian finite element approach. In the PFEM, the mesh nodes move according to the equation of motion, behaving similar to particles and transporting their momentum together with all their physical properties. At the same time, the mesh distortion problem is overcome by generating a new mesh with a Delaunay triangulation. Unlike standard finite element schemes, the PFEM keeps the nodes of the previous mesh fixed to avoid remapping from mesh to mesh. More details about the PFEM theory are given in Section 2.3.4. PFEM is known as a numerical tool for modeling fluid and solid under big configuration changes. Examples include the application in multi-fluid flows accounting for surface tension effects [34], granular flows [35], insertion problems in geotechnical engineering [36], landslide wave generation [37], industrial forming processes [38], and fresh cement paste [39].
Within the past 15 years, several researchers have used SPH to simulate orthogonal cutting both with commercial software [40,41,42] and in-house code [4,12,43,44,45,46]. The vast majority of these papers do do not consider heat transfer between bodies or consider that the friction coefficient is constant. A friction coefficient dependent on the temperature and relative velocity was proposed in [12,47]. Previous research shows that including the dependence of friction on temperature and relative velocity improves the prediction of process forces and temperatures.
More recent than SPH, PFEM has also been used to simulate orthogonal cutting in 2D [48,49,50] and 3D [51]. Through the early developments, the method has allowed for modeling continuous [48] and segmented chips [49], including the heat transfer between the chip and the tool with a constant friction coefficient. Although previous studies have demonstrated the capability of PFEM in simulating metal cutting processes, none of them account for temperature-dependent friction coefficients. In the present work, we fill this gap by presenting a PFEM-based orthogonal cutting simulation framework featuring a temperature-dependent friction model with the coefficient values identified in [11]. The prediction accuracy of PFEM are then compared with the SPH simulation results and experimental measurements taken from the literature. The thermomechanical model created for both PFEM and SPH contains: (1) heat generation due to the plastic and frictional work, (2) heat transfer from the workpiece to the tool, (3) heat conduction in the whole system, and (4) heat loss via convection and radiation over open surfaces. No damage or fracture is taken into account.
The remainder of the manuscript is organized as follows. Section 2 presents the simulation framework for metal cutting problems with SPH and FEM. As the core of this paper, Section 3 presents the numerical results alongside experimental validation by means of force and temperature measurements. Finally, the paper closes by summarizing the key findings of this research and discussing some opportunities for further improvements, given in Section 4.

2. Numerical Simulation Framework

Orthogonal cutting simulations involve complex thermal, mechanical, and material phenomena. This section provides an overview about the numerical modeling of machining with PFEM and SPH without describing their derivation procedures in detail. Figure 2 illustrates the setup of an orthogonal cutting model with the mechanical and thermal boundary conditions. More details about the PFEM and SPH model of metal cutting are given in Section 2.3.5 and Section 2.4.3.

2.1. Balance Equations

The balance equations for a thermomechanically coupled modeling of machining in the Lagrangian representation follow the basic conservation laws (i.e., mass, momentum, and energy) and are given by:
d ρ d t = ρ · v ¯
d r ¯ d t = v ¯
ρ d v ¯ d t = · σ = + b ¯
ρ c p d T d t = k 2 T + Q s
where ρ is the density, v ¯ is the velocity vector, σ = is the Cauchy stress tensor, b ¯ is the body force per unit of volume, c p is the specific heat capacity, T is the temperature, k is the thermal conductivity, Q s is the heat source due to plasticity, and r ¯ is the position vector. The heat source due to friction is included as a boundary condition at the tool–chip interface.
The term Q s in Equations (4) and (14) can be calculated from
Q s = χ σ y ε ¯ ˙ p
where χ dimensionless parameter with range between 0 and 1, specifying the fraction of plastic work that is converted into heat, σ y is the yield stress of the material and ε ¯ ˙ p is the equivalent plastic strain rate.
The boundary conditions associated with Equations (2)–(4), (12) and (14) includes the heat loss that occurs via radiation and convection, the heat generation at the tool–chip interface due to friction force, imposed temperature and imposed displacement/velocity on fixed boundaries (see Figure 2). The heat loss due to convection and radiation can be expressed as:
q l = h c ( T s T ) + ϵ σ ( T s 4 T 4 )
in which h c is the heat transfer coefficient, T s is the temperature at the surface of the body, T is the background temperature, ϵ is the emissivity factor, and σ is the Stefan–Boltzmann constant. Heat generated due to friction can be written as:
q f = η f f r i c · v r e l
where η represents the fraction of frictional work converted into heat, f f r i c is the magnitude of friction force, and v r e l is the magnitude of the relative velocity of the chip over the tool.
In the following sections, besides an overview of material and friction modeling aspects of metal cutting, we remark on the different PDEs (but equivalent) solved by each of the numerical schemes to clarify that the displacement and pressure fields are the primary variables in PFEM, whereas the primary variable in SPH is velocity.

2.2. Material and Friction Modeling

This section describes the constitutive model for the workpiece as well as the friction model that will be used between the chip–tool interface.

2.2.1. Flow Stress Modeling

The Johnson–Cook (JC) model [52] is used to define the plastic response of the workpiece. The von Mises stress σ y is the product of components describing strain hardening ( f 1 ) as a function of plastic strain ε ¯ p , strain rate sensitivity ( f 2 ) as a function of strain rate ε ¯ ˙ p , and thermal softening ( f 3 ) as a function of temperature T:
σ y = [ A + B ( ε ¯ p ) n ] f 1 · 1 + C ln ε ¯ ˙ p ε 0 ˙ f 2 · 1 T T r T m T r m f 3
where A is the initial yield strength of the material at room temperature, B is the hardening modulus, n is the work-hardening exponent, C is the coefficient dependent on the strain rate, ε 0 is the reference strain rate, m is the thermal softening coefficient, T r is the room temperature, and T m is the melting temperature. The values of JC parameters (A, B, C, m, and  n) are taken from the numerical–experimental identification database of [11] and listed in Table 1.
The JC model predicts that the greater the deformations, the greater the stress. However, for steel materials starting at a strain of approximately ε ¯ p = 2 , the formation of new dislocations and the mutual extinction of existing dislocations remain in a state of equilibrium as deformation progresses. As a result, the yield stress remains constant for higher strain values, as was shown in [11]. To describe the behavior described above for AISI 1045 steel, the JC model of AISI 1045 will not increase for strains ε ¯ p for ε ¯ p > 2 as follows:
σ y ( ε ¯ p > 2 ) = [ A + B ( ε ¯ p = 2 ) n ] · 1 + C ln ε ¯ ˙ p ε 0 ˙ · 1 T T r T m T r m
The model defined above for AISI 1045 steel is referred to as the modified JC variant along this work.

2.2.2. Friction Modeling

The friction force f f r i c is approximated by Coulomb’s law at the interface between the chip and the tool:
f f r i c = μ ( T ) f c o n t
where f c o n t represents the normal contact force and μ ( T ) is a temperature-dependent friction coefficient. The dependency of the friction coefficient on temperature follows an expression similar to the one used to model thermal softening in the Johnson–Cook flow stress model, which reads:
μ ( T ) = μ 0 1 T T r T m T r q
where μ 0 is the friction coefficient at T = T r and q is a factor that describes how fast friction decreases as temperature increases. The parameters used in the numerical simulations were obtained from the literature [11] and are presented in Table 2.

2.3. Particle Finite Element Modeling

This section presents an overview of machining simulation using PFEM that describes the balance equations, the spatial discretization, the temporal discretization, the meshing procedure, and the PFEM simulation setup.

2.3.1. Balance Equations

In the PFEM formulation of metal cutting, the dynamic equilibrium, mass conservation, and energy balance in a Lagrangian frame consist of:
ρ 2 u t 2 = · σ = + b
p = κ ln ( J ) 3 κ α ( 1 ln ( J ) ) J ( T T 0 )
ρ c p D T D t = k 2 T + Q s
where u ¯ is the displacement vector, σ = is the Cauchy stress tensor, b is the body forces per unit of current volume, p is the pressure, κ is the Bulk modulus, α is the thermal expansion coefficient, J is the determinant of the deformation gradient, T 0 is the room temperature, c p is the specific isobaric heat capacity, T is the temperature k, the thermal conductivity and Q s is the heat source due to plasticity given by Equation (5). The coupling between Equations (12) and (13) is given by the decomposition of the Cauchy stress tensor σ = into isotropic and deviatoric parts:
σ = = p I = + s =
where I = is the second-order identity tensor and s = is the deviatoric stress tensor.
Equation (12) of the PFEM is mathematically equivalent to introduce the advection equation (Equation (2)) into the momentum equation (Equation (3)). Equation (1) is equivalent to Equation (13), the difference is due the solid approach followed in the PFEM to the modeling of the hydrostatic part of the Cauchy Stress tensor, meanwhile, SPH uses a fluid flow approach (Equation (1)). The reason to introduce the continuity equation and the decomposition of the stress tensor into its deviatoric and hydrostatic components is because we use a mixed displacement-pressure formulation with linear triangle finite elements here that handles the incompressibility constraint due to J 2 plasticity.

2.3.2. PFEM Spatial Discretization of Metal Cutting Equations

The semi-discretized dynamic equilibrium and mass conservation equation are given by:
M u ¯ ¨ + F ¯ int , m e c = F ¯ e x t , m e c
M p p ¯ F p , v o l + F p , s t a b = 0
where M is the mass matrix for displacement, u ¯ ¨ is the nodal acceleration vector, F ¯ int , m e c is the vector of internal forces, F ¯ e x t , m e c is the vector of external forces, M p is the mass matrix for pressure, F p , v o l corresponds to the spatial discretization of the right-side term of Equation (13), and F p , s t a b is the stabilization term to overcome volumetric locking phenomena and pressure oscillations. A stabilization technique ( F p , s t a b term) based on the polynomial pressure projection [53] is used in this work. More information about the spatial discretization of the dynamical equilibrium and mass conservation in PFEM is given in [5,48,49].
The thermal problem appears in the resolution of semi-discretized heat equation:
C T ¯ ˙ + F ¯ int , t h e r = F ¯ e x t , t h e r
where C is the heat capacity matrix, T ¯ is the nodal temperature, F ¯ int , t h e r is the vector of internal thermal forces, and F ¯ e x t , t h e r is the vector of external thermal forces. The above global force vectors are obtained as the assemblies of element vectors as it is standard in the finite elements method. In this work, the element force vectors are evaluated using Gaussian quadratures.

2.3.3. Temporal Discretization in PFEM

The PFEM uses unconditionally stable implicit staggered thermomechanical algorithms, including inertial effects to integrate conservation equations presented in Section 2.3.1. It is important to briefly mention that SPH, unlike PFEM, uses an explicit time integration scheme (more information in Section 2.4.2).
In the PFEM, the thermomechanical problem is solved with an implicit staggered scheme that splits it in two sub-problems: the mechanical and the thermal part. At first, the dynamic equilibrium and mass conservation are solved at fixed temperature in a coupled way and, then, the energy balance is solved for a fixed domain. The coupling terms are included in the temperature-dependent constitutive law (see Section 2.2.1) and in the heat equation solution in the updated configuration.
Using the trapezoidal rule approximation, the time discretization of Equations (16) and (17) are written as
M n + 1 u ¯ ¨ n + 1 + F ¯ n + 1 int , m e c = F ¯ n + 1 e x t , m e c
M n + 1 p p ¯ n + 1 F n + 1 p , v o l + F n + 1 p , s t a b = 0
where M n + 1 is the mass matrix for displacement evaluated at time n + 1 , u ¯ ¨ n + 1 is the nodal acceleration vector at time n + 1 , F ¯ n + 1 int , m e c is the vector of internal forces at time n + 1 , F ¯ n + 1 e x t , m e c is the vector of external forces at time n + 1 , M n + 1 p is the mass matrix for pressure at time n + 1 , F p , v o l corresponds to the spatial discretization of the right side term of Equation (13) at time n + 1 , F n + 1 p , s t a b is the stabilization term at time n + 1 .
The velocity and acceleration at time n + 1 reads as
u ¯ ˙ n + 1 = u ¯ ˙ n + 2 Δ t ( u ¯ n + 1 u ¯ n )
a ¯ ¨ n + 1 = a ¯ ¨ n 4 Δ t u ¯ ˙ n + 4 Δ t 2 ( u ¯ n + 1 u ¯ n )
where u ¯ ˙ represents the vector of nodal velocities at time n + 1 . Equations (19) and (20) are solved in a coupled way using a Newton–Raphson scheme [50].
The thermal equation is solved at each time step by an implicit Backward Euler scheme numerical scheme as follows:
C ( T ¯ n + 1 T ¯ n Δ t ) + F ¯ n + 1 i n t , t h e r m = F ¯ n + 1 e x t , t h e r m
where T ¯ n + 1 is the nodal temperature at time n + 1 , T ¯ n is the nodal temperature at time n, Δ t is the time interval, F ¯ n + 1 i n t , t h e r m is the vector of internal thermal forces evaluated at time n + 1 , and v n + 1 e x t , t h e r m is the vector of external thermal forces evaluated at time n + 1 . The above equation is solved with the Newton–Raphson method. Further derivation and theoretical details about the Newton–Raphson method applied to the solution of the thermal problem can be found in [50].

2.3.4. PFEM Meshing

For completeness, the basics steps of the PFEM modeling framework are outlined (Algorithm 1). Given a set of particles that represent the workpiece at the beginning of the simulation.
Algorithm 1: Particle Finite Element Method Algorithm in Metal Cutting Simulations.
s t e p 1
while  s t e p n s t e p  do
1.
    Generate a Delaunay triangulation using the particles as nodes.
2.
    Identify the external and internal boundaries of the computational domain using the a l p h a -shape (see Figure 3 and Figure 4).
3.
    Implicitly solve the Lagrangian form of the mechanical and mass problem.
4.
    Update particles position, velocity, acceleration, pressure, and temperature.
5.
    Set tool incremental displacement = tool incremental displacement + time step tool displacement.
    if remeshing is needed then
        Go to step 1.
    else
        Go to step 3.
    end if
     n s t e p n s t e p + 1
end while
The Delaunay triangulation of a set of points that represent a piece is shown in Figure 3, and the identification of the external/internal boundary using the a l p h a -shape is shown in Figure 3. It is important to remark that sometimes the a l p h a -shape algorithm does not match perfectly the real internal and external boundaries, as it is shown in Figure 3, where the α shape method cannot perfectly recover the shape of the letter M. Some improvements to the predicted shape of the body can be obtained by changing the value of the alpha parameter as it is shown, for example, in [54]. However, if α is too large, some distorted or really large elements can be incorporated into the mesh. On the other hand, if α is too small, many elements are removed creating nonphysical holes. Figure 3 shows an example where the choice of α can perfectly predict the point-set geometry. The α applied to the prediction of the chip shape in the numerical modeling of metal cutting processes is shown in Figure 4. Figure 4 shows that the α shape is able to predict the chip shape with high accuracy.
Figure 3. Generation of Delaunay triangulation using particles as nodes and identification of the external/internal boundary using the α -shape algorithm for two different geometries.
Figure 3. Generation of Delaunay triangulation using particles as nodes and identification of the external/internal boundary using the α -shape algorithm for two different geometries.
Materials 16 03702 g003
Figure 4. The α -shape method applied to a triangulated metal cutting geometry.
Figure 4. The α -shape method applied to a triangulated metal cutting geometry.
Materials 16 03702 g004

2.3.5. PFEM Simulation Setup

Two-dimensional simulations of orthogonal cutting of Ti6Al4V titanium and AISI 1045 steel alloys using PFEM were carried out with the material properties, friction parameters, and process parameters given in Table A1, Table A2, and Table A3. Within the PFEM modeling framework, the workpiece was modeled as a thermo elastic–plastic body, whereby the tool is modeled as an elastic body with a cutting radius of r = 10 μ m. Figure 5 illustrates the initial setup and boundary conditions. A constant temperature of T = 293 K was imposed on the fixed surfaces (see Figure 5) as the Dirichlet boundary condition. Further, a Neumann boundary condition was imposed on the free surfaces accounting for convection and radiation heat losses (see Figure 2). The PFEM simulation runs until the tool has reached 2/3 of the workpiece length. This choice of the cutting length is made to ensure that the results are not affected by the fixed boundary condition (see the left boundary of the workpiece in Figure 2). Moreover, with the previous cutting length in all the numerical simulations a steady-state value of cutting forces is reached.
With respect to the piece geometry, the total height was chosen to be 3 times the cutting depth, and the length of the piece was set to be 3 mm for all the numerical simulations. For the discretization of the workpiece, a graded distance of particles was chosen (see Figure 5), which provides smaller distance between particles in areas with high gradients of the internal variables (typically the primary and the secondary shear zone) and where the curvature of the surface change along the simulation in order to increase the accuracy of the results. A minimum of 10, 15, and 20 particles were used along the uncut chip thickness for the cases in which t u = 0.1, 0.15, and 0.20 mm, respectively. The grading distance between particles used in the PFEM offers an advantage compared to other particle methods where the distance between particles is more or less constant along the simulation. Further, the same distance between nodes near the tool tip size was used to mesh the tool using linear triangle finite elements for solving displacement, elastic stresses, and temperatures. The numerical models was solved using a in-house Matlab PFEM code developed by the authors. Preliminary version of the PFEM software have been reported in [48,50,55].

2.4. Smoothed Particle Hydrodynamics Modeling

The Smoothed Particle Hydrodynamics (SPH) method, introduced in 1977 by Monaghan and Gingold [18] and Lucy [19], is a mesh-free Lagrangian discretization technique suitable for simulating solid and fluid flows. In SPH, the computational domain is discretized only by a set of points called particles without using a background mesh. In this section, the SPH approximations required for cutting simulations are described in brief.

2.4.1. SPH Spatial Discretization of Metal Cutting Equations

The SPH method approximates a field quantity f using a smoothed interpolation over the particles. This allows for a particle-based approximation of the partial differential equations and is well suited for large deformations problems such as metal cutting. Using the kernel approximation and subsequent particle approximation, the basic function interpolation over the bounded domain Ω can be converted into a finite, discretized form:
f i = Ω f ( x ) δ ( x x ) d x Ω f ( x ) W ( x x , h ) d x kernel approx . j = 1 N f j W i j V j particle approx .
where δ represents the Dirac delta function, W ( x x , h ) is the SPH kernel function (denoted as W i j = W ( x i x j , h ) ) with a smoothing length h, N is the total number of neighboring particles within the support domain of particle i, and V j = m j / ρ j is the integration weight.
As a Lagrangian mesh-free method particularly suitable for large deformations and large strains problems, SPH is used alongside PFEM for chip formation simulations. Previous studies such as [4,40,46] have proven the efficiency of SPH in machining modeling and also shown its potential for parallel computing [12]. A brief overview of the fundamental equations required for solving Equations (1) and (2) with SPH follows.
According to [56,57], the SPH approximates the gradient f and the Laplacian 2 f for a particle i in the following way:
f i ρ i j f i ρ i 2 + f j ρ j 2 W i j m j
2 f i j 2 f i j | r i j | e ¯ i j · W i j V j
where W i j = W h ( r ¯ i r ¯ j , h ) is a kernel function with h as a smoothing length and e ¯ i j = r ¯ i j / | r i j | is a unitary vector. At the same time, the notation ( ) i j ( ) j ( ) j is used for brevity. The Wendland quintic kernel function W is used in this paper is, which was originally developed by [58]. In [59], it was stated that previous kernel is the optimal choice for numerical stability of SPH.
The SPH approximations presented in Equations (24)–(26) are applied to the balance equations given in Section 2.1, obtaining, as a result, the semi-discretized form of the equations as:
d ρ i d t ρ i j v ¯ i j · W i j V j
d v ¯ i d t j σ = i ρ i 2 + σ = j ρ j 2 + Π i j I = + Λ = i j stabilizers · W i j m j + b ¯ i m i
d T i d t a i j 2 T i j | r i j | e ¯ i j · W i j V j + S i
where a i = k / ( ρ c p ) is the thermal diffusivity of particle i, S i = Q s / ( ρ c p ) its heat source term due to plasticity, is the artificial viscosity, and Λ = is the artificial stress tensor. The last two terms are used for numerical stability (see their original derivation in [60,61]). A deep explanation stabilizers and their tuning parameters for SPH cutting models can be found in [12,45]. The parameters used in this work are the same as used in the mentioned references.

2.4.2. Temporal Discretization

Here, a second-order leapfrog is used as an explicit time integration scheme for evolving the system of SPH particles. The maximum time step in this method is calculated a priori using a Courant–Friedrichs–Lewy (CFL) criterion through:
Δ t CFL × min Δ t t , Δ t m Δ t m a x
where the time step constraint of thermal Δ t t and mechanical Δ t m problem are included as [62,63]:
Δ t t = ρ c p h 2 k
Δ t m = h ( c 0 + | v ¯ | m a x )
in which c 0 is the sound speed, and the CFL coefficient is taken between 0.25 and 0.40.

2.4.3. SPH Simulation Setup

The thermomechanical SPH model of the orthogonal metal cutting test (see Figure 2) was built by discretizing the tool and workpiece with particles of uniform size. For solving the governing mechanical PDEs, we assumed the cutting tool to be a perfectly rigid body in contact with a deformable workpiece. As a result, only the workpiece was discretized by SPH particles for the mechanical parts of the simulations, while SPH interactions were inactive on the tool side. On the contrary, the governing thermal equations were solved for the entire geometry. Figure 6 illustrates the initial step of the present SPH simulation and how the thermomechanical boundary conditions are imposed on the respective surfaces.
In the SPH model setup, the workpiece dimensions were considered to be sufficiently large such that unwanted boundary effects are alleviated. To that end, we defined a 3 mm long workpiece with a height of at least 3× higher than the chip thickness in all test cases. The smoothing length factor was h = 1.3 Δ x for all thermal and mechanical interactions. It is worth mentioning that the solver generated an identical particle spacing for discretizing the cutting tool and workpiece to ensure numerical consistency.
Furthermore, the same distance between particles was used to discretize the tool to solve the heat equation. Properties of the workpiece materials are tabulated in Table A1 and Table A2. A reduced heat capacity approach [64] was considered, in order to decrease the computing time for reaching a steady-state tool temperature. The tool is modeled as a rigid body with a radius of r = 10 μ m.
In order to enable fast simulations within manageable runtimes, we utilized the GPU-accelerated version of the open-source metal cutting code mfree-iwf, Available on GitHub: https://github.com/iwf-inspire/mfree_iwf-ul_cut_gpu (accessed 10 January 2023). This efficient SPH solver was developed from [4,12,45,47,65,66] and can simulate 2D cutting and chip formation processes in the order of several minutes, depending on the process settings, desired resolution, and, of course, the hardware specifications. In this work, we ran all SPH parallel simulations on an nVidia™ Quadro GP100 graphics card, which is suited for scientific computing with double-precision calculations.

3. Results and Discussion

3.1. Prediction of Forces

As a first validation of the numerical simulations, forces predicted by the PFEM and SPH are compared against the experimental measurements provided by Afrasiabi et al. [11]. The comparison is made for the cutting forces F c and the passive forces F p . For each material, a combination of three cutting speeds v c and three cutting depths t u was carried out. In total, for two numerical methods, two materials and nine combinations of cutting velocity/uncut chip thickness, 36 numerical simulations were performed.
Figure 7 and Figure 8 show the force predictions for AISI 1045 and Ti6Al4V using PFEM and SPH. The results show that, as expected, the cutting forces increases when the uncut chip thickness increases. Further, Figure 7 indicates that the predicted cutting force decrease as the cutting velocity increase. Meanwhile, for SPH and Ti6Al4V, increasing the cutting speed has a significant effect on the predicted forces in most cases (see Figure 8). At the same time, the SPH method reveals a negligible influence of the cutting speed on the predicted cutting force for AISI 1045, the same situation can be observed in the case of using PFEM for Ti6Al4V.
Figure 9 shows that the highest percentage error for AISI 1045 in the predicted forces using the PFEM takes places at a cutting velocity of 180 m/min and cutting depth of 0.1 mm with a relative error of 33.37% for cutting forces F c and 64.96% for the passive forces F p . Figure 9 shows that the error for AISI 1045 in the predicted passives forces F p is around two times the error in the cutting forces F c . The big errors in the prediction of F p can be attributed to the in-exactitude of friction modeling at the tool–chip interface, especially at high contact pressures and very high temperatures. Furthermore, the large error can explained by the big deviation between the calculated and the measured contact length between chip and tool. The errors in the predictions of the cutting forces F c and F p for Ti6Al4V (see Figure 10) using the PFEM are similar to the ones obtained for AISI 1045 being greater for passive forces than for cutting forces.
According to the Figure 9, for the SPH simulation of AISI 1045, the maximum error occurs at a cutting velocity of 120 m/min and a cutting depth of 0.1 mm for F c (45.9%) and for F p it occurs at a cutting speed of 60 m/min and a cutting depth of 0.2 mm with 38.6%. Figure 10 also shows the error in the cutting force F c and passive force F p using the SPH for Ti6Al4V. F c indicates the highest error at a cutting speed of 20 m/min and a cutting depth of 0.15 mm with 63.4%, meanwhile, F p reaches a highest error of 53.7%. More details about the predicted forces is given in Appendix B.
The above results show that for both materials, SPH predicts the passive forces better than PFEM. On the other hand, the PFEM outperforms SPH in terms of predicting the cutting forces. Although both numerical methods implement the same constitutive and frictional model, the differences in force predictions can be explained with the difference between the algorithms that are used to identify the boundary particles for imposing the boundary conditions in each of the numerical methods. That is, PFEM uses the α -shape [67] whereas SPH constructs a free-surface detection algorithm. In future work, a study of the influence of the free-surface detection algorithm on the predicted processes needs to be carried out.
The interpretation of the errors obtained in the prediction of the forces must take into account the strong dependence on the magnitude of the forces in the choice of the parameters of the JC model (see [55] for more details), because only changing the JC parameters can generate a significant alteration of errors. It is important to clarify that the objective of the work was to analyze the advantages and disadvantages of the numerical modeling approaches presented and not to carry out a sensitivity analysis to the input parameters (JC parameters and friction law). Follow-up developments could consider identifying the JC parameters from the cutting experiments instead of the standard material tests such as SHPB in order to improve the prediction accuracy, because the extreme strains, strain rates, and temperatures encountered in metal cutting can be reflected in the identified JC material parameters in that way. Similarly, the predictions for Ti6Al4V could be improved including some form of strain-softening [49] to obtain serrated chip formation, because the JC model alone is not able to capture chip serration [68].

3.2. Prediction of Thermal Loads

Numerical simulations were also validated with experimentally measured temperatures taken from [11]. In the following figures, T R F refers to rake face temperature and T F C S refers to free chip surface temperature. Figure 11 and Figure 12 show that, generally, predictions of the rake temperatures T R F are significantly more accurate for AISI 1045 than for Ti6Al4V in both numerical methods. It can be seen that SPH performs more accurate in predicting the free chip surface temperature T F C S of Ti6Al4V whereas PFEM predicts T F C S for AISI 1045 better than SPH. The predicted temperatures used for calculating the errors shown in the previous figures are provided in Appendix B.
Although both numerical methods implement the same constitutive and frictional model, the differences in temperature predictions are due to the fact that the methods use different algorithms to identify the particles at the boundary. A different boundary particle identification implies that the contact area between the chip and the tool is different, therefore the heat flux would be different.
The suggestions to improve the prediction of the forces made in the previous section will also help to reduce the errors in the prediction of the temperature due to the fact that coupled thermomechanical phenomena occur in the metal cutting.

3.3. Prediction of Chip Shapes

Figure 13 shows a comparison of the predicted chips by both numerical methods for Ti6Al4V at a cutting speed of 40 m/s and a uncut chip thickness of 0.1 mm. The deformed chip thickness for the PFEM is larger compared to the SPH. The deformed chip thickness for the PFEM is larger compared to the SPH; therefore, the chip predicted by the SPH bends faster. This can be explained with the temperature predicted by the PFEM at the chip–tool interface, which, being higher than the SPH, generates higher coefficients of friction and therefore makes it more difficult for the chip to slide over the tool in PFEM simulations.

4. Conclusions

This paper investigated the numerical modeling of orthogonal cutting of two workpiece materials (i.e., AISI 1045 steel and Ti6Al4V) using the PFEM and the (SPH). The purpose of this work was to compare the capabilities of these two methods, which are known to be efficient for simulating large-deformation problems, in predicting forces, temperatures, and chip shapes in orthogonal cutting of two materials with different mechanical and thermal properties.
We presented the first PFEM-SPH comparison of metal cutting simulations, validated by force and temperature measurements, representing an advancement in numerical models of machining operations. For AISI 1045, the present numerical simulation framework predicted the thermal and mechanical loads of the cutting experiments with an average error below 40%. For Ti6Al4V, the prediction errors of thermomechanical loads for both numerical methods were significantly higher. This inaccuracy suggests that friction and plasticity models are far too simple to realistically describe the complex behavior at high temperatures, strains, and strain rates encountered in cutting.
The methodological difference between the two methods stems from their different time integration algorithms, i.e., implicit PFEM and explicit SPH. This has a direct impact on the calculation of contact and friction forces, thus justifying the significant discrepancy observed in some numerical cutting experiments. No concluding remark about the computational effort and runtime can be made as the two numerical methods were written in different languages (C++/CUDA vs. Matlab) and implemented on dissimilar computing platforms (GPU vs. CPU). A potential benefit of PFEM over SPH, nevertheless, appears to be the non-uniform distribution of particles enabled by the PFEM mesh structure, thereby minimizing the number of degrees of freedom required for modeling the cutting process (see Figure 5 and Figure 6).
Future development will focus on incorporating more sophisticated material models that describe the cutting process more realistically. It is also possible to speed up the PFEM calculations using GPU and other parallel computing techniques. Finally, a free-surface detection algorithm similar to SPH can be implemented into the PFEM code in order to bring the performance analysis to a more comparable state.

Author Contributions

Conceptualization, J.M.R.P., S.L. and M.A.; methodology, J.M.R.P., S.L. and M.A.; software, J.M.R.P. and M.A.; validation, J.M.R.P., S.L. and M.A.; formal analysis, J.M.R.P., S.L. and M.A.; investigation, J.M.R.P., S.L. and M.A.; writing—original draft preparation, J.M.R.P., S.L. and M.A.; writing—review and editing, J.M.R.P., S.L. and M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The last author (M.A.) would like to show his gratitude to Matthias Röthlin, Hagen Klippel, Konrad Wegener, and Markus Bambach for their invaluable contributions to developing the iMFREE software tool at ETH Zurich.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Thermomechanical Properties

Appendix A.1. Workpiece

Table A1. Material properties of AISI 1045 and Ti6Al4V used for simulations, taken from [11].
Table A1. Material properties of AISI 1045 and Ti6Al4V used for simulations, taken from [11].
PropertySymbolUnitAISI 1045Ti6Al4V
Density ρ kg / m 3 78704430
Young’s modulusEGPa205110
Poisson ratio ν 0.290.342
Heat conductivitykW/(m K)50.77.44
Specific heat capacity c p J/(kg K)486553
Melting temperature T m K17731877

Appendix A.2. Tool

Table A2. Material properties of the carbide cutting tool, taken from [11].
Table A2. Material properties of the carbide cutting tool, taken from [11].
PropertySymbolUnitValue
Density ρ k g / m 3 15.700
Young’s modulusEGPa600
Poisson ratio ν 0.2
Heat conductivitykW/(m K)59
Specific heat capacity c p J/(kg K)15

Appendix A.3. Tool–Chip Contact Zone

Table A3. Friction properties at the tool–chip interface, taken from [11].
Table A3. Friction properties at the tool–chip interface, taken from [11].
PropertySymbolUnitValue
Reference temperature T r K293
Coeff. of plastic work into heatc0.9
Coeff. of frictional work into heath1

Appendix B. Predicted Forces and Temperatures

Table A4. Predicted forces for AISI 1045 for nine simulations obtained from the combination of three cutting velocities and three uncut chip thicknesses using PFEM.
Table A4. Predicted forces for AISI 1045 for nine simulations obtained from the combination of three cutting velocities and three uncut chip thicknesses using PFEM.
Case v c (m/min) t u (mm) γ θ f r ( μ m) F c  (N) F p  (N) F p / F c
1600.101010499.73311.350.62
2 0.15 689.69427.60.61
3 0.2 843.01492.790.58
41200.101010407.6210.910.51
5 0.15 590.68323.790.54
6 0.2 788.58470.890.59
71800.101010350.48152.760.43
8 0.15 517.79224.60.43
9 0.2 696.91322.720.46
Table A5. Predicted forces for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using PFEM.
Table A5. Predicted forces for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using PFEM.
Case v c (m/min) t u (mm) γ θ f r ( μ m) F c  (N) F p  (N) F p / F c
10200.101010383.24164.620.42
11 0.15 535.51223.840.41
12 0.2 624.17228.180.36
13400.101010352.62169.890.48
14 0.15 482.86213.440.44
15 0.2 623.92263.250.42
16600.101010337155.60.46
17 0.15 464.26189.020.40
18 0.2 608240.140.39
Table A6. Predicted forces for AISI 1045 for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH.
Table A6. Predicted forces for AISI 1045 for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH.
Case v c (m/min) t u (mm) γ θ f r ( μ m) F c  (N) F p  (N) F p / F c
1600.101010703.54323.850.45
2 0.15 967.55418.720.43
3 0.2 1191.91492.560.41
41200.101010717.60342.500.48
5 0.15 958.99422.120.44
6 0.2 1217.72498.550.41
71800.101010661.76285.420.43
8 0.15 965.00424.380.44
9 0.2 1222.50507.840.42
Table A7. Predicted forces for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH.
Table A7. Predicted forces for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH.
Case v c (m/min) t u (mm) γ θ f r ( μ m) F c  (N) F p  (N) F p / F c
10200.101010678.79250.460.37
11 0.15 1027.53398.650.39
12 0.2 1319.81494.790.37
13400.101010550.84156.250.28
14 0.15 883.82255.020.29
15 0.2 1155.34391.250.34
16600.101010528.64119.940.23
17 0.15 800.95179.740.22
18 0.2 1011.86185.780.18
Table A8. Percent error of force prediction for AISI 1045 for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using PFEM. Mean error for F c = 21.57%, Mean error for F p = 43.81%.
Table A8. Percent error of force prediction for AISI 1045 for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using PFEM. Mean error for F c = 21.57%, Mean error for F p = 43.81%.
Case F c  (N) F p  (N) F p / F c
119.66%36.33%20.75%
217.20%32.87%18.92%
317.51%38.55%25.51%
417.15%47.79%36.98%
519.64%41.87%27.67%
616.46%29.72%15.87%
733.37%64.96%47.42%
830.59%57.62%38.95%
922.57%44.64%28.51%
Table A9. Percent error of force prediction for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using PFEM. Mean error for F c = 18.77%, Mean error for F p = 36.47%.
Table A9. Percent error of force prediction for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using PFEM. Mean error for F c = 18.77%, Mean error for F p = 36.47%.
Case F c  (N) F p  (N) F p / F c
1035.81%45.67%15.37%
1114.86%37.12%26.15%
1216.55%32.69%19.34%
1318.75%33.64%18.32%
1419.79%36.66%21.04%
1514.30%28.27%16.30%
1618.20%39.92%26.55%
1720.09%43.74%29.60%
1810.59%30.60%22.38%
Table A10. Percent error of force prediction for AISI 1045 for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH. Average error F c = 26.92%, Average Error F p = 26.55%.
Table A10. Percent error of force prediction for AISI 1045 for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH. Average error F c = 26.92%, Average Error F p = 26.55%.
Case F c  (N) F p  (N) F p / F c
113.1%33.8%41.4%
216.2%34.3%43.4%
316.6%38.6%47.3%
445.9%15.2%41.9%
530.5%24.2%41.9%
629.0%25.6%42.3%
725.8%34.5%48.0%
829.4%19.9%38.1%
935.8%12.9%35.9%
Table A11. Percent error of force prediction for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH. Mean error for F c = 44.54%, Mean error for F p = 32.41%.
Table A11. Percent error of force prediction for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH. Mean error for F c = 44.54%, Mean error for F p = 32.41%.
Case F c  (N) F p  (N) F p / F c
1013.7%17.3%27.3%
1163.4%12.0%31.5%
1276.4%46.0%17.3%
1326.9%39.0%51.9%
1446.8%24.3%48.5%
1558.7%6.6%32.8%
1628.3%53.7%63.9%
1737.9%46.5%61.2%
1848.8%46.3%63.9%
Table A12. Predicted temperatures for AISI 1045 for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using PFEM.
Table A12. Predicted temperatures for AISI 1045 for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using PFEM.
Case v c (m/min) t u (mm) γ θ f r ( μ m) T RF  (K) T FCS  (K)
1600.101010850643
2 0.15 863646
3 0.2 880660
41200.1010101008668
5 0.15 1017667
6 0.2 1025680
71800.1010101050710
8 0.15 1080712
9 0.2 1120714
Table A13. Predicted temperatures for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using PFEM.
Table A13. Predicted temperatures for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using PFEM.
Case v c (m/min) t u (mm) γ θ f r ( μ m) T RF  (K) T FCS  (K)
10200.101010900880
11 0.15 953926
12 0.2 1120940
13400.1010101150955
14 0.15 1250950
15 0.2 1260961
16600.1010101280975
17 0.15 1350983
18 0.2 1400990
Table A14. Percent error in temperature prediction for AISI 1045 for 9 simulation obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using PFEM. Mean error for T R F = 14.36%, Mean error for T FCS = 3.61%.
Table A14. Percent error in temperature prediction for AISI 1045 for 9 simulation obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using PFEM. Mean error for T R F = 14.36%, Mean error for T FCS = 3.61%.
Case T RF  (K) T FCS  (K)
127.25%0.62%
218.22%1.10%
320.22%2.64%
426.79%0.45%
513.89%1.52%
63.02%4.78%
713.51%4.72%
84.75%8.21%
91.63%8.51%
Table A15. Percent error in temperature prediction for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using PFEM. Mean error for T R F = 40.16%, Mean error for T F C S = 38.26%.
Table A15. Percent error in temperature prediction for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using PFEM. Mean error for T R F = 40.16%, Mean error for T F C S = 38.26%.
Case T RF  (K) T FCS  (K)
1028.39%23.77%
1138.12%34.79%
1245.27%40.93%
1346.87%38.81%
1443.18%42.86%
1540.00%43.86%
1639.43%40.90%
1739.90%42.46%
1840.28%35.99%
Table A16. Predicted temperatures for AISI 1045 for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH.
Table A16. Predicted temperatures for AISI 1045 for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH.
Case v c (m/min) t u (mm) γ θ f r ( μ m) T RF  (K) T FCS  (K)
1600.101010671.50522.08
2 0.15 659.18509.08
3 0.2 618.31476.13
41200.101010743.82534.79
5 0.15 728.02513.68
6 0.2 676.39487.28
71800.101010796.50532.21
8 0.15 798.56546.06
9 0.2 725.65493.09
Table A17. Predicted temperatures for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH.
Table A17. Predicted temperatures for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH.
Case v c (m/min) t u (mm) γ θ f r ( μ m) T RF  (K) T FCS  (K)
10200.1010101289.11801.86
11 0.15 1376.49774.27
12 0.2 1334.20747.94
13400.1010101316.03777.58
14 0.15 1361.25767.15
15 0.2 1399.15772.01
16600.1010101303.87792.86
17 0.15 1350.74787.08
18 0.2 1334.11762.14
Table A18. Percent error in temperature prediction for AISI 1045 for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH. Mean error for T R F = 16.55%, Mean error for T F C S = 21.80%.
Table A18. Percent error in temperature prediction for AISI 1045 for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH. Mean error for T R F = 16.55%, Mean error for T F C S = 21.80%.
Case T RF  (K) T FCS  (K)
10.5%19.3%
29.7%20.3%
311.3%26.0%
46.4%20.3%
518.5%21.8%
632.0%24.9%
713.9%21.5%
822.5%17.0%
934.2%25.1%
Table A19. Percent error in temperature prediction for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH. Mean error for T R F = 61.28%, Mean error for T F C S = 12.77%.
Table A19. Percent error in temperature prediction for Ti6Al4V for 9 simulations obtained from the combination of 3 cutting velocities and 3 uncut chip thicknesses using SPH. Mean error for T R F = 61.28%, Mean error for T F C S = 12.77%.
Case T RF  (K) T FCS  (K)
1083.9%12.8%
1199.5%12.7%
1273.0%12.1%
1368.1%13.0%
1455.9%15.4%
1555.5%15.6%
1642.0%14.6%
1740.0%14.1%
1833.7%4.7%

References

  1. Rodríguez, J.; Carbonell, J.; Jonsen, P. Numerical methods for the modelling of chip formation. Arch. Comput. Methods Eng. 2020, 27, 387–412. [Google Scholar] [CrossRef]
  2. Arrazola, P.; Özel, T.; Umbrello, D.; Davies, M.; Jawahir, I. Recent advances in modelling of metal machining processes. CIRP Ann.-Manuf. Technol. 2013, 62, 695–718. [Google Scholar] [CrossRef]
  3. Sadeghifar, M.; Sedaghati, R.; Jomaa, W.; Songmene, V. A comprehensive review of finite element modeling of orthogonal machining process: Chip formation and surface integrity predictions. Int. J. Adv. Manuf. Technol. 2018, 96, 3747–3791. [Google Scholar] [CrossRef]
  4. Afrasiabi, M.; Röthlin, M.; Klippel, H.; Wegener, K. Meshfree simulation of metal cutting: An updated Lagrangian approach with dynamic refinement. Int. J. Mech. Sci. 2019, 160, 451–466. [Google Scholar] [CrossRef]
  5. Rodriguez, J.; Jonsén, P.; Svoboda, A. Simulation of metal cutting using the particle finite-element method and a physically based plasticity model. Comput. Part. Mech. 2017, 4, 35–51. [Google Scholar] [CrossRef]
  6. Afrasiabi, M. Thermomechanical Simulation of Manufacturing Processes Using GPU-Accelerated Particle Methods. Ph.D. Thesis, ETH Zurich, Zurich, Switzerland, 2020. [Google Scholar] [CrossRef]
  7. Korkmaz, M.E.; Gupta, M.K. A State of the Art on Simulation and Modelling Methods in Machining: Future Prospects and Challenges. Arch. Comput. Methods Eng. 2022, 30, 1–29. [Google Scholar] [CrossRef]
  8. Rodriguez, J.M.; Jonsén, P.; Svoboda, A. Dislocation Density Based Material Model Applied in PFEM-simulation of Metal Cutting. Procedia CIRP 2017, 58, 193–197. [Google Scholar] [CrossRef]
  9. Wedberg, D.; Svoboda, A.; Lindgren, L.E. Modelling high strain rate phenomena in metal cutting simulation. Model. Simul. Mater. Sci. Eng. 2012, 20, 085006. [Google Scholar] [CrossRef]
  10. Kalhori, V.; Wedberg, D.; Lindgren, L.E. Simulation of mechanical cutting using a physical based material model. Int. J. Mater. Form. 2010, 3, 511–514. [Google Scholar] [CrossRef]
  11. Afrasiabi, M.; Saelzer, J.; Berger, S.; Iovkov, I.; Klippel, H.; Röthlin, M.; Zabel, A.; Biermann, D.; Wegener, K. A Numerical-Experimental Study on Orthogonal Cutting of AISI 1045 Steel and Ti6Al4V Alloy: SPH and FEM Modeling with Newly Identified Friction Coefficients. Metals 2021, 11, 1683. [Google Scholar] [CrossRef]
  12. Afrasiabi, M.; Meier, L.; Röthlin, M.; Klippel, H.; Wegener, K. GPU-accelerated meshfree simulations for parameter identification of a friction model in metal machining. Int. J. Mech. Sci. 2020, 176, 105571. [Google Scholar] [CrossRef]
  13. Arrazola, P.J.; Ozel, T. Investigations on the effects of friction modeling in finite element simulation of machining. Int. J. Mech. Sci. 2010, 52, 31–42. [Google Scholar] [CrossRef]
  14. Arrazola, P.; Ugarte, D.; Dominguez, X. A new approach for the friction identification during machining through the use of finite element modeling. Int. J. Mach. Tools Manuf. 2008, 48, 173–183. [Google Scholar] [CrossRef]
  15. Klocke, F.; Döbbeler, B.; Peng, B.; Schneider, S. Tool-based inverse determination of material model of Direct Aged Alloy 718 for FEM cutting simulation. Procedia CIRP 2018, 77, 54–57. [Google Scholar] [CrossRef]
  16. Zanger, F.; Bollig, P.; Schulze, V. Simulative investigations on different friction coefficient models. Procedia CIRP 2017, 58, 140–145. [Google Scholar] [CrossRef]
  17. Liang, X.; Liu, Z.; Wang, B.; Chunjin, W.; Cheung, B. Friction behaviors in the metal cutting process: State of the art and future perspectives. Int. J. Extrem. Manuf. 2022, 5, 012002. [Google Scholar] [CrossRef]
  18. Gingold, R.A.; Monaghan, J.J. Smoothed particle hydrodynamics: Theory and application to non-spherical stars. Mon. Not. R. Astron. Soc. 1977, 181, 375–389. [Google Scholar] [CrossRef]
  19. Lucy, L.B. A numerical approach to the testing of the fission hypothesis. Astron. J. 1977, 82, 1013–1024. [Google Scholar] [CrossRef]
  20. Monaghan, J.; Kos, A.; Issa, N. Fluid motion generated by impact. J. Waterw. Port Coast. Ocean. Eng. 2003, 129, 250–259. [Google Scholar] [CrossRef]
  21. Afrasiabi, M.; Mohammadi, S. Analysis of bubble pulsations of underwater explosions by the smoothed particle hydrodynamics method. In Proceedings of the ECCOMAS International Conference on Particle Based Methods, Barcelona, Spain, 25–27 November 2009. [Google Scholar]
  22. Marrone, S.; Colagrossi, A.; Le Touzé, D.; Graziani, G. Fast free-surface detection and level-set function definition in SPH solvers. J. Comput. Phys. 2010, 229, 3652–3663. [Google Scholar] [CrossRef]
  23. Afrasiabi, M.; Roethlin, M.; Wegener, K. Thermal simulation in multiphase incompressible flows using coupled meshfree and particle level set methods. Comput. Methods Appl. Mech. Eng. 2018, 336, 667–694. [Google Scholar] [CrossRef]
  24. Afrasiabi, M.; Chatzi, E.; Wegener, K. A Particle Strength Exchange Method for Metal Removal in Laser Drilling. Procedia CIRP 2018, 72, 1548–1553. [Google Scholar] [CrossRef]
  25. Afrasiabi, M.; Wegener, K. 3D thermal simulation of a laser drilling process with meshfree methods. J. Manuf. Mater. Process. 2020, 4, 58. [Google Scholar] [CrossRef]
  26. Alshaer, A.; Rogers, B.; Li, L. Smoothed Particle Hydrodynamics (SPH) modelling of transient heat transfer in pulsed laser ablation of Al and associated free-surface problems. Comput. Mater. Sci. 2017, 127, 161–179. [Google Scholar] [CrossRef]
  27. Jonsén, P.; Pålsson, B.I.; Häggblad, H.Å. A novel method for full-body modelling of grinding charges in tumbling mills. Minerals Eng. 2012, 33, 2–12. [Google Scholar] [CrossRef]
  28. Jonsson, P.; Andreasson, P.; Hellström, J.I.; Jonsén, P.; Staffan Lundström, T. Smoothed Particle Hydrodynamic simulation of hydraulic jump using periodic open boundaries. Appl. Math. Model. 2016, 40, 8391–8405. [Google Scholar] [CrossRef]
  29. Jonsson, P.; Jonsén, P.; Andreasson, P.; Lundström, T.S.; Hellström, J.G.I. Modelling dam break evolution over a wet bed with smoothed particle hydrodynamics: A parameter study. Engineering 2015, 7, 248. [Google Scholar] [CrossRef]
  30. Forsström, D.; Lindbäck, T.; Jonsén, P. Prediction of wear in dumper truck body by coupling SPH-FEM. Tribol. Mater. Surf. Interfaces 2014, 8, 111–115. [Google Scholar] [CrossRef]
  31. Larsson, S.; Gustafsson, G.; Häggblad, H.Å.; Jonsén, P. Experimental and numerical study of potassium chloride flow using smoothed particle hydrodynamics. Minerals Eng. 2018, 116, 88–100. [Google Scholar] [CrossRef]
  32. Idelsohn, S.R.; Oñate, E.; Pin, F.D. The particle finite element method: A powerful tool to solve incompressible flows with free-surfaces and breaking waves. Int. J. Numer. Methods Eng. 2004, 61, 964–989. [Google Scholar] [CrossRef]
  33. Idelsohn, S.R.; Oñate, E.; Del Pin, F.; Calvo, N. Fluid–structure interaction using the particle finite element method. Comput. Methods Appl. Mech. Eng. 2006, 195, 2100–2123. [Google Scholar] [CrossRef]
  34. Mier-Torrecilla, M.d.; Idelsohn, S.; Oñate, E. Advances in the simulation of multi-fluid flows with the particle finite element method. Application to bubble dynamics. Int. J. Numer. Methods Fluids 2011, 67, 1516–1539. [Google Scholar] [CrossRef]
  35. Larsson, S.; Rodríguez Prieto, J.M.; Gustafsson, G.; Häggblad, H.Å.; Jonsén, P. The particle finite element method for transient granular material flow: Modelling and validation. Comput. Part. Mech. 2021, 8, 135–155. [Google Scholar] [CrossRef]
  36. Monforte, L.; Arroyo, M.; Carbonell, J.M.; Gens, A. Numerical simulation of undrained insertion problems in geotechnical engineering with the Particle Finite Element Method (PFEM). Comput. Geotech. 2017, 82, 144–156. [Google Scholar] [CrossRef]
  37. Mulligan, R.P.; Franci, A.; Celigueta, M.; Take, W.A. Simulations of landslide wave generation and propagation using the particle finite element method. J. Geophys. Res. Ocean. 2020, 125, e2019JC015873. [Google Scholar] [CrossRef]
  38. Oñate, E.; Franci, A.; Carbonell, J.M. A particle finite element method for analysis of industrial forming processes. Comput. Mech. 2014, 54, 85–107. [Google Scholar] [CrossRef]
  39. Reinold, J.; Meschke, G. Particle finite element simulation of fresh cement paste–inspired by additive manufacturing techniques. PAMM 2019, 19, e201900198. [Google Scholar] [CrossRef]
  40. Limido, J.; Espinosa, C.; Salaün, M.; Lacome, J.L. SPH method applied to high speed cutting modelling. Int. J. Mech. Sci. 2007, 49, 898–908. [Google Scholar] [CrossRef]
  41. Villumsen, M.F.; Fauerholdt, T.G. Simulation of Metal Cutting Using Smooth Particle Hydrodynamics. LS-DYNA Anwenderforum, C-III. 2008, p. 17. Available online: https://www.dynamore.it/de/download/papers/forum08/dokumente/C-III-02.pdf (accessed on 10 January 2023).
  42. Demiral, M. 1365. Smoothed particle hydrodynamics modeling of vibro-assisted turning of Ti alloy: Influence of vibration parameters. J. Vibroeng. 2014, 16, 2685–2695. [Google Scholar]
  43. Eberhard, P.; Gaugele, T. Simulation of cutting processes using mesh-free Lagrangian particle methods. Comput. Mech. 2013, 51, 261–278. [Google Scholar] [CrossRef]
  44. Spreng, F.; Eberhard, P. Machining process simulations with smoothed particle hydrodynamics. Procedia CIRP 2015, 31, 94–99. [Google Scholar] [CrossRef]
  45. Röthlin, M.; Klippel, H.; Afrasiabi, M.; Wegener, K. Metal cutting simulations using smoothed particle hydrodynamics on the GPU. Int. J. Adv. Manuf. Technol. 2019, 102, 3445–3457. [Google Scholar] [CrossRef]
  46. Roethlin, M.; Klippel, H.; Afrasiabi, M.; Wegener, K. Meshless single grain cutting simulations on the GPU. Int. J. Mechatron. Manuf. Syst. 2019, 12, 272–297. [Google Scholar] [CrossRef]
  47. Afrasiabi, M.; Klippel, H.; Roethlin, M.; Wegener, K. Smoothed Particle Hydrodynamics Simulation of Orthogonal Cutting with Enhanced Thermal Modeling. Appl. Sci. 2021, 11, 1020. [Google Scholar] [CrossRef]
  48. Rodriguez, J.M.; Carbonell, J.; Cante, J.; Oliver, J. Continuous chip formation in metal cutting processes using the Particle Finite Element Method (PFEM). Int. J. Solids Struct. 2017, 120, 81–102. [Google Scholar] [CrossRef]
  49. Prieto, J.R.; Carbonell, J.M.; Cante, J.; Oliver, J.; Jonsén, P. Generation of segmental chips in metal cutting modeled with the PFEM. Comput. Mech. 2017, 61, 639–655. [Google Scholar]
  50. Rodriguez, J.; Larsson, S.; Carbonell, J.; Jonsén, P. Implicit or explicit time integration schemes in the PFEM modeling of metal cutting processes. Comput. Part. Mech. 2022, 9, 709–733. [Google Scholar] [CrossRef]
  51. Carbonell, J.M.; Rodríguez, J.; Oñate, E. Modelling 3D metal cutting problems with the particle finite element method. Comput. Mech. 2020, 66, 603–624. [Google Scholar] [CrossRef]
  52. Johnson, G.R. A constitutive model and data for materials subjected to large strains, high strain rates, and high temperatures. In Proceedings of the Seventh International Symposium on Ballistics, The Hague, The Netherlands, 19–21 April 1983; pp. 541–547. [Google Scholar]
  53. Dohrmann, C.R.; Bochev, P.B. A stabilized finite element method for the Stokes problem based on polynomial pressure projections. Int. J. Numer. Methods Fluids 2004, 46, 183–201. [Google Scholar] [CrossRef]
  54. Cremonesi, M.; Franci, A.; Idelsohn, S.; Oñate, E. A state of the art review of the particle finite element method (PFEM). Arch. Comput. Methods Eng. 2020, 27, 1709–1735. [Google Scholar] [CrossRef]
  55. Rodríguez, J.M.; Larsson, S.; Carbonell, J.M.; Jonsén, P. Dislocation density based flow stress model applied to the PFEM simulation of orthogonal cutting processes of Ti-6Al-4V. Materials 2020, 13, 1979. [Google Scholar] [CrossRef]
  56. Brookshaw, L. A method of calculating radiative heat diffusion in particle simulations. In Proceedings of the Astronomical Society of Australia; Cambridge University Press: Cambridge, UK, 1985; Volume 6, pp. 207–210. [Google Scholar]
  57. Monaghan, J.J. Smoothed particle hydrodynamics. Annu. Rev. Astron. Astrophys. 1992, 30, 543–574. [Google Scholar] [CrossRef]
  58. Wendland, H. Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree. Adv. Comput. Math. 1995, 4, 389–396. [Google Scholar] [CrossRef]
  59. Dehnen, W.; Aly, H. Improving convergence in smoothed particle hydrodynamics simulations without pairing instability. Mon. Not. R. Astron. Soc. 2012, 425, 1068–1082. [Google Scholar] [CrossRef]
  60. Monaghan, J.; Gingold, R. Shock simulation by the particle method SPH. J. Comput. Phys. 1983, 52, 374–389. [Google Scholar] [CrossRef]
  61. Gray, J.; Monaghan, J.; Swift, R. SPH elastic dynamics. Comput. Methods Appl. Mech. Eng. 2001, 190, 6641–6662. [Google Scholar] [CrossRef]
  62. Cleary, P.W.; Monaghan, J.J. Conduction modelling using smoothed particle hydrodynamics. J. Comput. Phys. 1999, 148, 227–264. [Google Scholar] [CrossRef]
  63. Monaghan, J.J. Smoothed particle hydrodynamics. Rep. Prog. Phys. 2005, 68, 1703. [Google Scholar] [CrossRef]
  64. Childs, T.H.; Arrazola, P.J.; Aristimuno, P.; Garay, A.; Sacristan, I. Ti6Al4V metal cutting chip formation experiments and modelling over a wide range of cutting speeds. J. Mater. Process. Technol. 2018, 255, 898–913. [Google Scholar] [CrossRef]
  65. Afrasiabi, M.; Roethlin, M.; Wegener, K. Contemporary Meshfree Methods for Three Dimensional Heat Conduction Problems. Arch. Comput. Methods Eng. 2020, 27, 1413–1447. [Google Scholar] [CrossRef]
  66. Afrasiabi, M.; Klippel, H.; Röthlin, M.; Wegener, K. An improved thermal model for SPH metal cutting simulations on GPU. Appl. Math. Model. 2021, 100, 728–750. [Google Scholar] [CrossRef]
  67. Calvo, N.; Idelsohn, S.R.; Oñate, E. The extended Delaunay tessellation. Eng. Comput. 2003, 20, 583–600. [Google Scholar] [CrossRef]
  68. Childs, T. Adiabatic Shearing in Metal Machining. In CIRP Encyclopedia of Production Engineering; Laperrière, L., Reinhart, G., Eds.; Springer: Berlin/Heidelberg, Germany, 2014; pp. 27–33. [Google Scholar] [CrossRef]
Figure 1. Orthogonal metal cutting process and its numerical modeling output at the scale of the tool. The left-most image created by Florian Schott, downloaded from https://commons.wikimedia.org/wiki/File:SchlichtenDrehen.jpg (accessed on 3 January 2023) under CC BY-SA 4.0.
Figure 1. Orthogonal metal cutting process and its numerical modeling output at the scale of the tool. The left-most image created by Florian Schott, downloaded from https://commons.wikimedia.org/wiki/File:SchlichtenDrehen.jpg (accessed on 3 January 2023) under CC BY-SA 4.0.
Materials 16 03702 g001
Figure 2. Mechanical and thermal boundary conditions for the PFEM and SPH machining simulations.
Figure 2. Mechanical and thermal boundary conditions for the PFEM and SPH machining simulations.
Materials 16 03702 g002
Figure 5. Configuration of the thermomechanical PFEM model: Initial geometry and boundary conditions.
Figure 5. Configuration of the thermomechanical PFEM model: Initial geometry and boundary conditions.
Materials 16 03702 g005
Figure 6. Configuration of the thermomechanical SPH model: Initial geometry and boundary conditions.
Figure 6. Configuration of the thermomechanical SPH model: Initial geometry and boundary conditions.
Materials 16 03702 g006
Figure 7. Force prediction of PFEM and SPH for AISI 1045.
Figure 7. Force prediction of PFEM and SPH for AISI 1045.
Materials 16 03702 g007
Figure 8. Force prediction of PFEM and SPH for Ti6Al4V.
Figure 8. Force prediction of PFEM and SPH for Ti6Al4V.
Materials 16 03702 g008
Figure 9. Force errors of SPH and PFEM for AISI 1045.
Figure 9. Force errors of SPH and PFEM for AISI 1045.
Materials 16 03702 g009
Figure 10. Force errors of SPH and PFEM for Ti6Al4V.
Figure 10. Force errors of SPH and PFEM for Ti6Al4V.
Materials 16 03702 g010
Figure 11. Temperature errors of SPH and PFEM for AISI 1045.
Figure 11. Temperature errors of SPH and PFEM for AISI 1045.
Materials 16 03702 g011
Figure 12. Temperature errors of SPH and PFEM for Ti6Al4V.
Figure 12. Temperature errors of SPH and PFEM for Ti6Al4V.
Materials 16 03702 g012
Figure 13. Predicted chip shapes using PFEM and SPH: Ti6Al4V, cutting speed of 40 m/s, and a uncut chip thickness of 0.1 mm.
Figure 13. Predicted chip shapes using PFEM and SPH: Ti6Al4V, cutting speed of 40 m/s, and a uncut chip thickness of 0.1 mm.
Materials 16 03702 g013
Table 1. JC parameters for AISI 1045 and Ti6Al4V, taken from [11].
Table 1. JC parameters for AISI 1045 and Ti6Al4V, taken from [11].
PropertyAISI 1045Ti6Al4V
A (MPa)288940
B (MPa)695965
C0.03400.0387
m1.35580.9835
n0.28350.8372
ε 0 ˙ (1/s)0.0040.004
Table 2. Parameters of temperature friction law for AISI 1045 and Ti6Al4V, taken from [11].
Table 2. Parameters of temperature friction law for AISI 1045 and Ti6Al4V, taken from [11].
Material μ 0 q
AISI 10450.812.09
Ti6Al4V0.515.76
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rodríguez Prieto, J.M.; Larsson, S.; Afrasiabi, M. Thermomechanical Simulation of Orthogonal Metal Cutting with PFEM and SPH Using a Temperature-Dependent Friction Coefficient: A Comparative Study. Materials 2023, 16, 3702. https://doi.org/10.3390/ma16103702

AMA Style

Rodríguez Prieto JM, Larsson S, Afrasiabi M. Thermomechanical Simulation of Orthogonal Metal Cutting with PFEM and SPH Using a Temperature-Dependent Friction Coefficient: A Comparative Study. Materials. 2023; 16(10):3702. https://doi.org/10.3390/ma16103702

Chicago/Turabian Style

Rodríguez Prieto, Juan Manuel, Simon Larsson, and Mohamadreza Afrasiabi. 2023. "Thermomechanical Simulation of Orthogonal Metal Cutting with PFEM and SPH Using a Temperature-Dependent Friction Coefficient: A Comparative Study" Materials 16, no. 10: 3702. https://doi.org/10.3390/ma16103702

APA Style

Rodríguez Prieto, J. M., Larsson, S., & Afrasiabi, M. (2023). Thermomechanical Simulation of Orthogonal Metal Cutting with PFEM and SPH Using a Temperature-Dependent Friction Coefficient: A Comparative Study. Materials, 16(10), 3702. https://doi.org/10.3390/ma16103702

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop