Sage Journals: Discover world-class research

Abstract

In the process of high-efficiency milling of titanium alloys, cutting vibration causes the frequent changes in the contact relationship between the flank face of the cutter tooth and the workpiece. The instantaneous friction speed, friction force and energy consumption of the milling cutter thus change continuously, resulting in the difficulty to precisely control the wear and life of milling cutter. This has become a bottleneck of restricting the further improvement of the cutting performance of high energy efficiency milling cutter. In this work, the model of instantaneous friction velocity and friction energy consumption of the flank face of the cutter tooth under vibration are developed. The distribution and evolution of instantaneous friction speed and friction energy consumption are studied. The proposed models are verified by high-efficiency milling experiments. The results show that the average relative error of the proposed model is 10.05%, by using this model, the influences of the cutting vibration on instantaneous friction energy consumption and friction wear boundaries could be effectively unveiled. The proposed model also could accurately describe the distribution and evolution of the above-mentioned parameters.

Keywords

High energy efficiency milling cutter wear friction speed friction force energy consumption

Introduction

With its excellent cutting performance and lightweight design, high energy efficiency milling cutter is widely used in high-speed and high-efficiency machining of large aircraft parts represented by titanium alloy structural parts. In the process of high energy efficiency milling cutter processing, affected by factors such as cutter tooth error, milling vibration, etc., the instantaneous friction boundary between the milling cutter tooth and the workpiece and the boundary of the wear area show nonlinear changes with contact angle.^1–3 The degree of friction and wear of milling cutter and the overall service life are difficult to accurately predict, which has become a bottleneck restricting the further improvement of performance of high energy efficiency milling cutter.^4–6

In the process of high-efficiency milling, under the action of cutting force and centrifugal force, high-frequency and small-scale random vibrations are generated between the milling cutter and the workpiece, which makes the contact angle between the flank face of the milling cutter and the workpiece dynamically change, resulting in uncertainty in the relative speed, mechanical state, enclosed boundary and enclosed area of the contact area between the flank face of the tool tooth and the workpiece.^7,8 Therefore, it is difficult to accurately predict the instantaneous friction speed, friction force and energy consumption rate of milling cutter tooth. The above results, under the influence of the milling cutter’s long-term cutting, can easily lead to the complexity of the friction boundary, the friction speed evolution and the distribution state of the flank face.

At present, the measurement and characterization of the flank face friction of milling cutter tooth mainly focus on the maximum wear width, ignoring the characteristic information of the entire wear boundary of the flank face. Wang et al.⁹ examined the microstructure and compositional characterization of worn tool inserts using SEM and EDS analysis. Oliveira et al.¹⁰ tested different tools wear under finishing, roughing, and moderate machining conditions. Wu et al.¹¹ studied the effect of aspect ratio on the machining process during milling of titanium alloys. Saketi et al.¹² studied the wear characteristics of the cutting tool in the milling process of two different quenched cold work tool steels, High-resolution field emission gun scanning electron microscope, energy spectrometer, Auger electron spectrometer and optical surface profiler were used to detect worn blades. The results show that due to friction chemical reaction, adhesive wear, slight abrasive wear, tooth face wear and crater wear were controlled by continuous wear. Vipindas and Mathew¹³ studied the wear behavior of TiAlN-coated solid tungsten carbide tools in the process of milling Ti-6Al-4V, the tool wear was studied from the aspects of tool face wear, adhesive wear, abrasive wear, blade radius increase, coating layering, etc. The cutting tool wear and face finish were characterized by a three-dimensional optical profiler. Zheng et al.¹⁴ conducted a high-speed turning test under dry conditions. In the turning test, the cutting tooth flank wear width was measured and revealed the anti-friction performance and wear mechanism of TiAlN/TiN coated tools on high-strength steel. With the change of sliding speed and load, the wear amount and wear rate of the coating showed an opposite trend. An et al.¹⁵ studied the changes of Ti6Al4V side milling tool wear, cutting torque and surface topography under four sustainable conditions. A theoretical model of the plane wear width VB was established, and the average prediction error was 15.87%. The above expression method of friction width cannot be used to accurately identify the specific geometric characteristics of the instantaneous friction boundary of the flank face.

At present, in terms of energy consumption model, Liu et al.¹⁶ conducted experiments on the impact of tool wear on milling energy consumption and process sustainability. He pointed out that changes in process parameters and tool wear play a dominant role in energy consumption. He found that the increase in tool wear was the main reason for the increase in energy consumption. Li et al.¹⁷ established an energy-saving optimization model for NC turning batch processing parameters that considers tool wear. This model considers the synergistic effect of tool wear status and process parameters on energy consumption. He pointed out that adjusting the process parameters according to the tool wear status can make the tool more in the low-wear area to process, which effectively reduced the total energy consumption. Lv et al.¹⁸ had established three models for material removal energy consumption in the turning process: specific energy model (SEM), cutting force model (CFM), and exponential function model (EFM). He pointed out that the source of energy consumption model coefficients was the main reason that affects the prediction accuracy. Sealy et al.¹⁹ simulated the energy consumption of precision hard milling, and they found that tool wear had the greatest impact on net cutting energy consumption. However, these studies mainly focus on qualitative analysis. There is still a lack of a clear expression of friction energy consumption changing with friction boundary.

This paper first proposed a method of identifying and characterizing the friction boundary of the flank face of cutter tooth and a method of selecting feature points in the friction area. Then, the dynamic cutting model of the milling cutter under the action of vibration was establish. The thermal-stress couple field simulation was carried out. According to the simulation effective stress distribution, a method to identify the characteristic points of the instantaneous friction boundary of the flank face was proposed. In order to solve the distribution function of friction energy consumption more accurately, a method for calculating the friction speed and friction force at the feature points was proposed. The proposed model is verified by cutting experiments.

Characterization of friction boundary of flank face

In order to characterize the characteristic and the friction boundary of the cutter tooth, the cutter tooth coordinate was established according to the geometric structures of the cutter tooth, and the milling cutter structure coordinate considering the error of the cutter tooth was established, as shown in Figure 1. The parameters explanation are shown in Table 1.

Figure 1.

High energy efficiency milling cutter and cutter tooth coordinate.

Table 1.

Parameters explanation of high energy efficiency milling cutter and cutter tooth coordinate.

Serial numbers	parameters	Parameters explanation
1	r_max	Maximum radius of gyration at the tip of the cutter tooth.
2	r_i	The radius of gyration of the tip point of the ith cutter tooth of milling cutter, i = 1, 2, 3.
3	θ_i	The angle between the ith cutter tooth and the (i + 1)th cutter tooth.
4	α ₀	Mounting angle of cutter tooth.
5	Δz_i	Axial error of ith cutter tooth and (i + 1)th cutter tooth.
6	Δr_i	Radial error of ith cutter tooth and (i + 1)th cutter tooth.

In Figure 1 o₀-x₀y₀z₀ is the milling cutter structure coordinate. The coordinate origin o₀ is the intersection point between the center axis of the milling cutter and the plane where the lowest point of the cutter tooth is located. The z₀ axis direction is the direction of the center axis of the milling cutter, and the y₀ axis is determined according to the principle of the Cartesian coordinate. o_i-x_iy_iz_i is the cutter tooth coordinate. The x_i axis direction of the cutter tooth coordinate is the straight line direction of the two end points of the cutter tooth installation plane. The straight line direction passing through the inner cutter tip point of the cutter tooth and perpendicular to the cutter tooth installation plane is the y_i axis direction. The intersection point is the coordinate origin o_i. The direction of the z_i axis is determined according to the principle of the Cartesian coordinate.

The instantaneous cutting behavior of the milling cutter and its cutter teeth directly affects the contact relationship between the flank face and the workpiece, thereby affecting the friction speed and friction force changes. In order to reveal the dynamic cutting process of the milling cutter under the action of vibration, the instantaneous contact relationship between the milling cutter and the workpiece at different cutting positions was characterized, as shown in Figure 2. The parameters explanation are shown in Table 2.

Figure 2.

Dynamic cutting process of high energy efficiency milling cutter under vibration.

Table 2.

Parameters explanation of dynamic cutting process of high energy efficiency milling cutter under vibration.

Serial numbers	Parameters	Parameters explanation
1	o-xyz	Workpiece coordinate. Point o is the intersection of the three surfaces of the bottom, back andright side of the workpiece. The x axis is parallel to the feed direction of the milling cutter, they axis is parallel to the width of the workpiece, and the z axis is parallel to the depth of theworkpiece.
2	o_v-x_vy_vz_v	Milling cutter cutting coordinate without vibration. The coordinate origin o_v is the intersectionpoint between the center axis of the milling cutter body and the lowest plane of the milling cutter,the x_v axis direction is the milling cutter feed speed direction, the y_v axis direction is the millingcutter cutting width direction, and the z_v axis direction is the milling cutter cutting depth direction.
3	o_v′-x_v′y_v′z_v′	The milling cutter cutting coordinate without vibration. The coordinate origin o_v′ and x_v′, y_v′, z_v′axis is the deviation of the coordinate origin o_v and x_v, y_v, z_v axis of the milling cutter cuttingcoordinate without vibration.
4	x_L	The length of the workpiece.
5	y_L	The width of the workpiece.
6	z_L	The height of the workpiece.
7	v_f	The feed speed of the milling cutter.
8	n	Rotation speed of milling cutter.
9	a_p	Cutting depth.
10	a_e	Cutting width.
11	t	Cutting time of milling cutter.
12	ψ_i(t)	When the cutting time is t, the angle between the x₀ axis of the milling cutter structure coordinateand the x_v′ axis of the milling cutter cutting coordinate with vibration.
13	φ_i(t)	When the cutting time is t, the angle between the x₀ axis of the milling cutter structure coordinateand the x_i axis of the cutter tooth coordinate.
14	θ_d	Contact angle. The angle between the connection line between the projection point of the millingcutter tip point on the lowest plane of the milling cutter and the origin of the milling cutter cuttingcoordinate and the x_v axis of the milling cutter cutting coordinate.
15	δ(t)	The attitude angle of the milling cutter.
16	δ′(t)	The direction angle of the milling cutter.
17	δ₁(t)	The projection angle of the milling cutter attitude angle δ(t) on the x_vo_vz_v plane.
18	δ₂(t)	The projection angle of the milling cutter attitude angle δ(t) on the y_vo_vz_v plane.
19	A_x(t)	The vibration displacement of the milling cutter in the x_v axis direction.
20	A_y(t)	The vibration displacement of the milling cutter in the y_v axis direction.
21	A_z(t)	The vibration displacement of the milling cutter in the z_v axis direction.
22	l	Milling cutter overhang.

According to the dynamic cutting process of the milling cutter under the action of vibration, the trajectory of the coordinate origin of the milling cutter coordinate was calculated, as shown in equation (1).

{\begin{matrix} x = v_{f} \cdot t - r_{max} + A_{x} (t) \\ y = y_{L} - a_{e} + A_{y} (t) \\ z = z_{L} - a_{p} + A_{z} (t) \end{matrix}

(1)

The deflection angle of the milling cutter caused by vibration can be expressed by the attitude angle and the direction angle of the milling cutter. The solution of the attitude angle can be solved by first calculating its two projection angles, such as equations (2) and (3). Then the projection angles were synthesized and solved, as shown in equation (4). The solution of the direction angle is shown in equation (5).

δ_{1} (t) = \arctan (\frac{A_{x} (t)}{l - A_{z} (t)})

(2)

δ_{2} (t) = \arctan (\frac{A_{y} (t)}{l - A_{z} (t)})

(3)

δ (t) = \arctan \frac{\sqrt{{A^{2}}_{x} (t) + {A^{2}}_{y} (t)}}{l - A_{z} (t)}

(4)

δ^{'} (t) = {\begin{matrix} \arcsin \frac{A_{x} (t)}{\sqrt{A_{x}^{2} (t) + A_{y}^{2} (t)}}, A_{x} (t) \geq 0 \\ \arcsin \frac{A_{x} (t)}{\sqrt{A_{x}^{2} (t) + A_{y}^{2} (t)}} + π, A_{x} (t) < 0 \end{matrix}

(5)

The matrix transformation relationship between the coordinate in Figure 2 can be used to solve the trajectory equation of any point on the cutter tooth, as shown in equation (6).

{[\begin{matrix} x y z 1 \end{matrix}]}^{T} = A_{3} A_{2} B_{3} B_{4} B_{2} A_{1} B_{1} {[\begin{matrix} x_{i} & y_{i} & z_{i} & 1 \end{matrix}]}^{T}

(6)

In the equation, each matrix is specifically expressed as follows:

\begin{matrix} A_{1} = [\begin{matrix} 1 & 0 & 0 & r_{i} \sin φ_{i} \\ 0 & 1 & 0 & - r_{i} \cos φ_{i} \\ 0 & 0 & 1 & - Δ zi \\ 0 & 0 & 0 & 1 \end{matrix}] \\ B_{1} = [\begin{matrix} \cos φ_{i} & - \sin φ_{i} & 0 & 0 \\ \sin φ_{i} & \cos φ_{i} & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{matrix}] \\ A_{2} = [\begin{matrix} 1 & 0 & 0 & A_{x} (t) \\ 0 & 1 & 0 & A_{y} (t) \\ 0 & 0 & 1 & A_{z} (t) \\ 0 & 0 & 0 & 1 \end{matrix}] B 2 = [\begin{matrix} \cos ψ_{i} (t) & \sin ψ_{i} (t) & 0 & 0 \\ - \sin ψ_{i} (t) & \cos ψ_{i} (t) & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{matrix}] \\ B 2 = [\begin{matrix} \cos ψ_{i} (t) & \sin ψ_{i} (t) & 0 & 0 \\ - \sin ψ_{i} (t) & \cos ψ_{i} (t) & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{matrix}] \end{matrix}

\begin{matrix} B 3 = [\begin{matrix} 1 & 0 & 0 & 0 \\ 0 & \cos δ 1 (t) & - \sin δ 1 (t) & 0 \\ 0 & \sin δ 1 (t) & \cos δ 1 (t) & 0 \\ 0 & 0 & 0 & 1 \end{matrix}] \\ B 4 = [\begin{matrix} \cos δ 2 (t) & 0 & \sin δ 2 (t) & 0 \\ 0 & 1 & 0 & 0 \\ - \sin δ 2 (t) & 0 & \cos δ 2 (t) & 0 \\ 0 & 0 & 0 & 1 \end{matrix}] \\ A_{3} = [\begin{matrix} 1 & 0 & 0 & v_{f} \cdot t - r_{max} \\ 0 & 1 & 0 & y_{L} - a_{e} \\ 0 & 0 & 1 & z_{L} - a_{p} \\ 0 & 0 & 0 & 1 \end{matrix}] \end{matrix}

Equations (1)–(6) can be used to calculate the motion trajectory and attitude angle of the milling cutter and cutter tooth under the action of vibration, which reflects the irregular changes in the position and posture of the cutter tooth during the machining process. This irregular change directly leads to the uncertainty of the contact relationship between the flank face of the cutter tooth and the workpiece, which will cause many inflection points on the boundary of the friction area, and these inflection points constitute the friction boundary.

In order to identify the friction boundary and the inflection points, the basis for identification was proposed. The basis is that there is a clear difference in color between the friction and non-friction areas. The inflection points with obvious color gaps were identified. These points were connected to form the friction boundary, as shown in Figure 3.

Figure 3.

Characterization of the friction boundary and selection of its feature points.

The speed and force of points at different positions in the friction area are different. In order to comprehensively evaluate the distribution of speed and force in the entire friction area, a method for selecting feature points in the friction area of the flank face was proposed. The sections parallel to the plane y_io_iz_i were taken at equal intervals on the flank face, and the interval size was Δx_i. The intersection of any section x = x_ie and the friction upper boundary is denoted as e₁, and the intersection of any section with the friction lower boundary is denoted as e₂. e₁ and e₂ are the feature points of the friction upper boundary and the friction lower boundary. e is the point between points e₁ and e₂ on the section. The method for selecting feature points is also shown in Figure 3, and the parameters explanation are shown in Table 3.

Table 3.

Parameters explanation of the characterization of the friction boundary and selection of its feature points.

Serial numbers	Parameters	Parameters explanation
1	x_io_iy_i	Two-dimensional projection plane of cutter tooth coordinate o_i-x_iy_iz_i.
2	x_i0	The x_i axis value corresponding to the left edge of the friction area of the flank face ofthe cutter tooth.
3	x_i1	The x_i axis value corresponding to the right edge of the friction area of the flank face of thecutter tooth.
4	x_ie	Any x_i axis value between x_i₀ and x_i₁.
5	e₀	The intersection of an arbitrary section parallel to the plane y_io_iz_i and the original boundary ofthe cutting edge.
6	e₁	The intersection of an arbitrary section parallel to the plane y_io_iz_i and the friction upper boundary.
7	e₂	The intersection of an arbitrary section parallel to the plane y_io_iz_i and the friction lower boundary.
8	e	Any point between e₁ and e₂ on the flank face.
9	y_ie	The coordinate value of the y_i axis corresponding to the point e.
10	Δx_i	The distance between any two adjacent sections.

In Figure 3, y_i= y₀(x_i) (i = 1,2,3) is the original boundary of the cutting edge of the cutter tooth i. y_i = y_ri(x_i,t) (i = 1,2,3) is the friction upper boundary of the flank face of cutter tooth i when the cutting time of the milling cutter is t. y_i = y_hi(x_i,t) (i = 1,2,3) is the friction lower boundary of the flank face of cutter tooth i when the cutting time of the milling cutter is t.

The characterization method of the friction boundary and the selection method of the feature points provide a basis for the identification of the instantaneous friction boundary, and provide a basis for the construction of the friction energy consumption distribution function.

Identification of instantaneous friction boundary on flank face

The instantaneous friction boundary in the cutting process cannot be measured in real time through experiments. In order to obtain the instantaneous friction boundary at a certain moment, the stress of the cutter tooth during the cutting process was simulated.

The thermal coupling field simulation of the high efficiency milling cutter adopts the finite element simulation software Deform-3D. The milling cutter cutting process is simulated using its Machining [Cutting] module. UG was used to build tool and workpiece models and import them into the simulation module. The Johnson-Cook constitutive parameters of TC4 titanium alloy material are shown in Table 4.

Table 4.

Constitutive parameters of TC4 titanium alloy J-C.

Material	A (MPa)	B (MPa)	n	c	m
TC4	860	683	0.47	0.035	1.0

In the Table 4, A, B, n, c, m are yield stress intensity, strain hardening constant, strain hardening exponent, strain rate hardening parameter, and temperature strain rate sensitivity, respectively.

The boundary conditions of the simulation are shown in Table 5.

Table 5.

Boundary conditions of finite element simulation.

Toolsmaterial	Coatingmaterial	Thermal conductivityλ(w/m°C)	Rotating speedn (r/min)	feed rate v_f(mm/min)	Cutting widtha_e (mm)	Cutting deptha_p (mm)
WC	Ti-CN	45	1225	500	16	0.5

The number of meshes divided by the cutter teeth in the simulation process is 500,000. In order to make the simulation results more accurate, when the workpiece is meshed with absolute size, the convergence effect is best when the step distance does not exceed 1/3 of the minimum element size of the workpiece grid element.

The simulated stroke is 5 m. The effective stress distribution cloud diagram obtained by simulation is shown in Figure 4.

Figure 4.

Effective stress distribution cloud diagram of the flank face.

According to the distribution of effective stress in Figure 4, a method for judging the instantaneous friction boundary of the flank face was proposed.

The friction upper boundary was formed by the wear of the cutting edge. From the perspective of material mechanics, the wear of the cutting edge of the milling cutter is due to the cutter tooth material reaching its yield strength σ_s. Therefore, the yield strength can be used as the basis for judging the feature points of the instantaneous friction upper boundary of the flank face of the cutter tooth. If the effective stress of a point reaches the yield strength, this point is the feature point e₁ of instantaneous friction upper boundary. In order to find these feature points, the parallel sections were taken on the flank face. The sections x₁ = 6.36, x₁ = 7.36, and x₁ = 8.36 were taken as examples, the effective stress distribution along these sections were constructed, as shown in Figure 5.

Figure 5.

Effective stress distribution of the flank face.

In Figure 5, the coordinates of the feature points e₁ on the sections x₁ = 6.36, x₁ = 7.36, and x₁ = 8.36 were identified. More sections were taken, and the feature point e₁ on each section was judged by the method in Figure 5. When all the feature points were connected, the instantaneous friction upper boundary of the flank face could be constructed, as shown in Figure 6.

Figure 6.

Instantaneous friction upper boundary of flank face.

The instantaneous friction upper boundary in Figure 6 was identified by yield strength σ_s. Because the instantaneous friction lower boundary may not reach the yield strength, yield strength σ_s cannot be used to determine the instantaneous friction lower boundary.

Because the instantaneous friction lower boundary is the boundary between the friction area and the non-friction area of the flank face, to identify the feature point e₂ of it, it is necessary to analyze the stress of the friction area and the non-friction area. During the cutting process, stress concentration will be formed in the friction area, which leads to a large difference in the effective stress between the friction area and the non-friction area of the flank face. According to this stress distribution, as long as the position where the effective stress on the flank face changes the most, the feature point of the instantaneous friction lower boundary can be found. The position where the effective stress changes the most can be judged by solving the effective stress change rate.

The three sections in Figure 5 were taken as examples. The effective stress change rate at these sections was calculated, as shown in Figure 7.

Figure 7.

Effective stress change rate of flank face.

In Figure 7, the coordinates of the feature points e₂ on the sections x₁ = 6.36, x₁ = 7.36, and x₁ = 8.36 were identified. More sections were taken, and the feature point e₂ on each section was judged by the method in Figure 7. When all the feature points were connected, the instantaneous friction lower boundary of the flank face could be constructed, as shown in Figure 8.

Figure 8.

Instantaneous friction lower boundary of flank face.

The instantaneous friction lower boundary in Figure 8 and the instantaneous friction lower boundary in Figure 6 together constitute the boundary of the instantaneous friction area. This boundary defines the range of instantaneous friction energy consumption that needs to be calculated. The identification of the instantaneous friction boundary of the flank face of the cutter tooth is the basis for constructing the distribution function of the instantaneous friction energy consumption.

Instantaneous friction energy consumption distribution function on flank face

Determination of friction speed

In order to solve the distribution function of the instantaneous friction energy consumption of the flank face more accurately, the friction speed and friction force were calculated. The calculation of friction speed can be obtained by partial derivative of the cutter tooth trajectory with respect to the cutting time.

The equation (7) of the cutter tooth trajectory could be calculated by matrix operation on the equation (6).

{\begin{matrix} x = x (t) \\ y = y (t) \\ z = z (t) \end{matrix}

(7)

The equation (8) of the movement speed of any point on the cutter tooth in the directions of the three coordinate axes x, y, and z was calculated:

{\begin{matrix} v_{x} (t) = \frac{\partial x (t)}{\partial t} \\ v_{y} (t) = \frac{\partial y (t)}{\partial t} \\ v_{z} (t) = \frac{\partial z (t)}{\partial t} \end{matrix}

(8)

In Figure 9, G(x(t), y(t), z(t)) is the flank face equation of the cutter tooth. I(x(t), y(t), z(t)) is the transition face equation of the workpiece. A_α is the flank face. θ_vm(t) is the space angle projected by the flank face movement speed and the movement speed in the common tangent plane when the cutting time is t. v_m is the friction speed. −v_m and v_m are equal in size and opposite in direction. $\overset{⇀}{m}$ is the unit direction vector of the friction speed direction.

Figure 9.

Method for solving the friction speed of feature points in the friction area.

The direction of the friction speed on the flank face is determined by the contact relationship between the flank face and the workpiece. The direction of the friction speed should be the opposite direction of the projection of the movement speed on the common tangent plane. The common tangent plane could be solved by separately obtaining the partial derivatives of the two contact surfaces and then simultaneously solving the equations, as shown in equation (9).

{\begin{matrix} \frac{\partial G (x (t), y (t), z (t))}{\partial t} = 0 \\ \frac{\partial I (x (t), y (t), z (t))}{\partial t} = 0 \end{matrix}

(9)

The unit direction vector projected by the movement speed on the common tangent plane can be expressed as equation (10) in the workpiece coordinate.

\overset{⇀}{m} = (m_{x}, m_{y}, m_{z})

(10)

Where m_x, m_y, and m_z are the unit coordinate values of $\overset{⇀}{m}$ in the directions of the x, y, and z axes.

The movement speed of the feature point of the flank face was also expressed in the form of a vector in the workpiece coordinate:

{\overset{⇀}{v}}_{v} (t) = (v_{x} (t), v_{y} (t), v_{z} (t))

(11)

Where v_x(t), v_y(t), and v_z(t) are the velocity components of the movement speed v_v(t) at the feature points of the flank face in the directions of the x, y, and z axes.

The friction speed v_m(t) was calculated by equation (12).

v_{m} (t) = v_{v} (t) \cos (π - θ_{vm} (t))

(12)

Where θ_vm(t) can be calculated by equation (13).

\begin{matrix} θ_{vm} (t) = \arccos \\ (\frac{m_{x} \cdot v_{x} (t) + m_{y} \cdot v_{y} (t) + m_{z} \cdot v_{z} (t)}{\sqrt{m_{x}^{2} + m_{y}^{2} + m_{z}^{2}} \cdot \sqrt{v_{x}^{2} (t) + v_{y}^{2} (t) + v_{z}^{2} (t)}}) \end{matrix}

(13)

The friction speed of the feature point of the flank face was expressed as a vector in the workpiece coordinate:

{\overset{⇀}{v}}_{m} (t) = (v_{m x} (t), v_{m y} (t), v_{m z} (t))

(14)

Where v_mx(t), v_my(t), and v_mz(t) are the velocity components of the friction speed v_m(t) at the feature points of the flank face in the directions of the x, y, and z axes.

In order to calculate the friction energy consumption in the same reference coordinate, the friction speed was transformed into the cutter tooth coordinate through coordinate transformation, as shown in equation (15).

\begin{matrix} {[\begin{matrix} x_{i} & y_{i} & z_{i} & 1 \end{matrix}]}^{T} = B_{1}^{- 1} A_{1}^{- 1} B_{2}^{- 1} B_{4}^{- 1} B_{3}^{- 1} A_{2}^{- 1} A_{3}^{- 1} \\ {[\begin{matrix} v_{mx} (t) & v_{my} (t) & v_{mz} (t) & 1 \end{matrix}]}^{T} \end{matrix}

(15)

The vector of the friction speed of the feature point in the tooth coordinate was expressed as follows:

{\overset{⇀}{v}}_{m i} (x_{i}, y_{i}, z_{i}, t) = (v_{m x_{i}} (t), v_{m y_{i}} (t), v_{m z_{i}} (t))

(16)

Determination of friction force and friction energy consumption

In order to calculate the magnitude of the effective stress component in the friction speed direction, a decomposition model of the effective stress was established, as shown in Figure 10, the parameters explanation are listed in Table 6. The grid divided in the simulation process was a tetrahedron, so the vector sum of the stress along the three edges should be the same as the effective stress of the feature point, which converted the effective stress at the feature point into three directional stresses.

Figure 10.

Effective stress decomposition model on the flank face.

Table 6.

Parameters explanations of effective stress decomposition model on the flank face.

Serial numbers	Parameters	Parameters explanation
1	f ₁, f₂, f₃	The vector sum of the effective stress along the three sides of the tetrahedral element should bethe effective stress value at the intersection of the three sides.
2	θ_m1	The space angle between the effective stress f₁ and the friction speed v_mi in the cutter toothcoordinate.
3	θ_m2	The space angle between the effective stress f₂ and the friction speed v_mi in the cutter toothcoordinate.
4	θ_m3	The space angle between the effective stress f₃ and the friction speed v_mi in the cutter toothcoordinate.
5	v_mxi	The component of the friction speed v_mi in the cutter tooth coordinate in the direction ofthe x_i axis.
6	v_myi	The component of the friction speed v_mi in the cutter tooth coordinate in the directionof the y_i axis.
7	v_mzi	The component of the friction speed v_mi in the cutter tooth coordinate in the directionof the z_i axis.

According to the effective stress distribution of the flank face, the effective stresses f₁, f₂, and f₃ were expressed in vector in the cutter tooth coordinate, as shown in equation (17).

\begin{matrix} {\overset{⇀}{f}}_{1} = (x_{f 1}, y_{f 1}, z_{f 1}) \\ {\overset{⇀}{f}}_{2} = (x_{f 2}, y_{f 2}, z_{f 2}) \\ {\overset{⇀}{f}}_{3} = (x_{f 3}, y_{f 3}, z_{f 3}) \end{matrix}

(17)

In the equation, x_f1, y_f₁, and z_f₁ are the force components of the effective stress f₁ in the x_i, y_i and z_i axis directions, x_f₂, y_f₂, and z_f₂ are the force components of the effective stress f₂ in the x_i, y_i and z_i axis directions, x_f₃, y_f₃ and z_f₃ are the force component of the effective stress f₃ in the x_i, y_i, and z_i axis directions.

The three stresses were projected in the direction of the friction speed, and the calculation method was shown in equation (18).

f_{qv} (x_{i}, y_{i}, z_{i}, t) = f_{q} \cos θ_{mq} (t) (q = 1, 2, 3)

(18)

In the equation, f_q(q = 1,2,3) is the effective stress of the three sides of the tetrahedron. f_qv(x_i,y_i,z_i,t) is the component stress of f_q in the direction of friction speed. θ_mq(t) is the angle between the three effective stress vectors and the friction speed vector, and its magnitude was calculated by equation (19).

\begin{matrix} θ_{mq} (t) = \arccos \\ \frac{x_{fq} \cdot v_{m x_{i}} (t) + x_{fq} \cdot v_{m y_{i}} (t) + x_{fq} \cdot v_{m z_{i}} (t)}{\sqrt{x_{fq}^{2} + y_{fq}^{2} + z_{fq}^{2}} \cdot \sqrt{v_{m x_{i}}^{2} (t) + v_{m y_{i}}^{2} (t) + v_{m z_{i}}^{2} (t)}} \end{matrix}

(19)

The effective stress vector sum in the direction of the friction speed is the friction force, as shown in equation (20).

f_{m} (x_{i}, y_{i}, z_{i}, t) = \sum_{q = 1}^{3} f_{qv} (x_{i}, y_{i}, z_{i}, t)

(20)

According to equations (16) and (20), the instantaneous friction energy consumption distribution function of the flank face was constructed, as shown in equation (21).

P_{m} (x_{i}, y_{i}, z_{i}, t) = f_{m} (x_{i}, y_{i}, z_{i}, t) \cdot {v_{m}}_{i} (x_{i}, y_{i}, z_{i}, t)

(21)

The distribution function of instantaneous friction energy consumption on the flank face of the cutter tooth is constructed by equations (7)–(21). It can be used to calculate the instantaneous friction energy consumption at any point in the friction area of the flank face, and obtain the instantaneous friction energy consumption distribution state of the entire flank face.

Distribution and evolution of friction energy consumption on flank face

Distribution of friction energy consumption

Due to the direct contact between the flank face of the milling cutter and the machined surface, the magnitude and distribution of friction energy consumption on the flank face of the milling cutter directly affect the friction and wear state of the flank face and the quality of the machined surface. Therefore, this paper mainly studies the distribution and evolution of friction energy consumption on the flank face of milling cutter teeth.

In order to calculate the distribution and evolution of the friction energy consumption of the flank face, the instantaneous friction energy consumption is calculated according to the proposed distribution function of the instantaneous friction energy of the flank face. The distribution of the instantaneous friction energy consumption of the feature points along the y_i axis of the cutter tooth coordinate was constructed. The change of the instantaneous friction energy consumption of the feature points with the contact angle θ_d was constructed.

Three sections x₁= 4.36, x₁ = 5.36, and x₁ = 6.36 were taken on cutter tooth 1, and the distribution of friction speed and friction force along the y_i axis of the cutter tooth coordinate were drawn, as shown in Figures 11 and 12.

Figure 11.

Distribution of instantaneous friction speed on the flank face.

Figure 12.

Distribution of instantaneous friction force on the flank face.

Figure 11 shows the distribution of the instantaneous friction speed at the three sections on the flank face of the cutter tooth 1. The shaded area in the figure is the non-friction area on the flank face. In this area, the tooth does not contact the workpiece. So the friction speed in this area is zero. The boundary between the non-friction area and the friction area of the flank face is the friction lower boundary. Due to the influence of vibration during the cutting process, the contact relationship between the flank face of the cutter tooth and the workpiece is constantly changing. Therefore, the friction lower boundary changes irregularly at all times. And the feature point e₂ of the friction lower boundary at different sections is not at the same coordinate position. The friction upper boundary is formed by the wear of the original cutting edge. Since the occurrence of wear is also irregular, the feature point e₁ of the friction upper boundary at different sections is not at the same coordinate position.

In the friction area, the friction speed on the same section will not change much. The reason is that the friction speed of the feature point is a velocity component of the movement velocity. The changes of movement speed play a decisive role in the friction speed. A major component of the movement speed is the cutting speed produced by the milling cutter speed. The cutting speed is affected by the radius of gyration. The radius of gyration of the feature points on the same section is different, but the difference is very small. Therefore, the friction speed change trend on the same section is not obvious. The friction speed of each section is enlarged and expressed separately, and its changing trend can be clearly seen. The relative positional relationship between the cutter teeth and the cutter body determines that the radius of gyration of the feature points on the section gradually increases along the positive direction of the y₁ axis. Therefore, the friction speed on the same section increases along the positive direction of the y₁ axis. The friction speed on different sections varies greatly. This is also mainly determined by the radius of gyration. Since the radial position of different sections is very different, the influence of the radius of gyration on its speed is much greater. The section x₁ = 4.36 is the closest to the center of rotation, so the friction speed is the smallest. The section x₁ = 6.36 is the farthest from the center of rotation, so its friction speed is the largest.

The friction force distribution interval of the three sections in Figure 12 is consistent with the friction speed. The friction force only exists in the friction area. The distribution of friction force in different sections is basically the same. The friction force of different sections increases along the positive direction of the y₁ axis. During the cutting process, the friction area of the flank face of the cutter tooth interacts with the workpiece. However, because the area close to the friction upper boundary participates in the material removal of the workpiece, it will withstand greater interaction forces. This leads to greater friction force in the area close to the friction upper boundary. The friction force at the feature point near the friction upper boundary has a larger variation. The reason is that the part close to the friction upper boundary has a large interaction force with the workpiece when it participates in cutting, causing friction and friction changes to increase.

The instantaneous friction energy consumption of the entire friction area was calculated. The flank face was used as a reference surface with zero energy consumption. The instantaneous friction energy consumption surface of the flank face was fitted. The sections were taken on the friction energy consumption surface in the same way. The instantaneous distribution diagram of frictional energy consumption at the section along the y₁ axis is obtained, as shown in Figure 13.

Figure 13.

Distribution of instantaneous friction energy consumption on the flank face.

It can be seen from the fitted friction energy consumption distribution surface that the friction energy consumption of the friction area increases along the positive direction of the y₁ axis. The frictional energy consumption value reaches the maximum at the friction upper boundary. And the frictional energy consumption is the smallest at the friction lower boundary. This distribution of friction energy consumption is the same as the trend of friction force. The reason is that due to the fact that the friction upper boundary participates in the removal of workpiece material and bears greater stress. The friction energy consumption in the friction area is distributed in the direction of the x₁ axis with a middle height and a low side. Due to the irregular changes in the curved structure of the flank face and the contact relationship between the milling cutter and the workpiece, which makes the middle part of the friction area more likely to be squeezed. The middle part of the friction area is more stressed than the two sides, and the friction energy consumption is also greater.

Evolution of friction energy consumption

The friction speed, friction force and friction energy consumption of the feature points of the flank face of the cutter tooth 1 were calculated at different contact angles. The method of selecting feature points had been given in Part 1. This part selected three feature points on the same section as an example to calculate. They are the intersection point e₁ between the section x₁ = 5.36 and the friction upper boundary, the intersection point e₂ between the section x₁ = 5.36 and the friction lower boundary, and the point e between e₁ and e₂. The change of friction speed and friction force with contact angle at selected feature points are shown in Figures 14 and 15.

Figure 14.

Variation of friction speed of flank face.

Figure 15.

Variation of friction force of flank face.

On the whole, the friction speed changes periodically during the cutting process. In Figure 14, one cycle of friction speed is T. One cycle is divided into two parts, T₁ and T₂, by whether the friction speed is zero. The reason for the periodic change of friction speed is related to the structural parameters of the milling cutter itself and the cutting processing parameters. In the cutting process, each revolution of the milling cutter is a cycle T. Since the cutting width selected in the experiment is the same as the radius of gyration of the milling cutter, in a cycle T, the cutter starts to cut into the workpiece when the contact angle is 0, and cuts out the workpiece when the contact angle is π/2. The cutter tooth has friction speed on the flank face only when it participates in the cutting process. The part T₁ in Figure 14 is the friction speed at the feature points of the flank face when the tooth 1 participates in cutting.

It can be seen from the right part of Figure 14 that the friction speed of the feature points fluctuate irregularly with the contact angle. This volatility is mainly affected by two factors. On the one hand, the vibration speed generated by the vibration has an impact on the speed of the cutter tooth movement, which in turn affects the change of the friction speed; on the other hand, the instantaneous friction boundary at different contact angles is different, resulting in the position coordinates of the feature point in the cutter tooth coordinate. When the change occurs, the radius of gyration of the feature point also changes, so the change in the movement speed of the feature point causes the friction speed to eventually change. In general, these two aspects indirectly affect the friction speed by influencing the movement speed.

In Figure 15, the feature point friction also changes periodically with the contact angle. The reason is the same as that of the periodic change of the friction speed in Figure 14, which is determined by the rotation period of the milling cutter and the period of a single tooth participating in the cutting. The friction force of the feature point also fluctuates irregularly with the contact angle. However, the reason for the irregular fluctuation is different from the friction speed. The main reason is that the instantaneous friction boundary change causes the change of the effective stress distribution on the flank face. At the same time, the change of the contact relationship between the milling cutter and the workpiece also affects the effective stress transmission effect on the flank face. Eventually, the friction force produced irregular changes.

The instantaneous friction energy consumption at multiple contact angles was calculated and fitted to the instantaneous friction energy distribution surface. Sections and feature points were taken on the energy consumption surface using the same method. The curves of frictional energy consumption at feature points with contact angle was plotted, as shown in Figure 16.

Figure 16.

Variation of friction energy consumption of flank face.

It can be seen from Figure 16 that the areas and positions of friction energy consumption on the flank face at different contact angles are different. This is mainly determined by the change of the instantaneous friction boundary with the contact angle. The position of the maximum friction energy consumption of the entire flank face is also different. This shows that the most severe position of friction and wear at different contact angles has changed. The contact relationship between the flank face and the workpiece is constantly changing, resulting in uneven force distribution on the flank face, which changes the distribution of friction energy consumption on the entire flank face. The law of friction energy consumption at the feature point with the contact angle is similar to the law of friction in Figure 15. From Figure 14, it can be seen that the friction speed fluctuates only in a small range with the change of the contact angle, so the change law of friction energy consumption is more affected by the friction force. The change process of friction energy consumption is more susceptible to friction force.

Experimental verification

In order to verify the instantaneous friction energy consumption model, the cutting experiments of the high energy efficiency milling cutter were carried out. The cutting experiments were carried out on a three-axis milling machining center (VDL-1000E).

The energy efficient milling cutter used in the experiment is an indexable three-tooth equal-pitch milling cutter with a diameter of 32 mm produced by Walter, as shown in Figure 17. The titanium alloy material of model TC4 is used. The length, width and height of the titanium alloy workpiece are 250 × 100 × 20 mm.

Figure 17.

Energy efficient milling cutter.

The grade of titanium alloy used in this study is TC4, and its chemical composition is shown in Table 7.

Table 7.

Chemical composition of Ti-6Al-4V.

Element	Al	V	Fe	O	Si	C	N	H	Others	Ti
Contents (%)	5.5–6.8	3.5–4.5	0.3	0.2	0.15	0.1	0.05	0.01	0.5	Margin

During the experiments, the charge output PCB three-axis high-frequency acceleration sensor was used to collect the vibration acceleration signal at the workpiece surface to be processed. After the experiments, the DH5922 transient signal test and analysis system was used to process the vibration signal. The experiment setup is shown in Figure 18, and the specific experimental parameters are shown in Table 8.

Figure 18.

Experiment setup.

Table 8.

Experiment parameters.

Cutting stroke	Cutting parameters				Axial error (mm)			Radial error (mm)
L (m)	N (r/min)	v_f (mm/min)	a_p (mm)	a_e (mm)	Δz_i ¹	Δz_i ²	Δz_i ³	Δr_i ¹	Δr_i ²	Δr_i ³
5	1225	500	0.5	16	0	0.006	0.004	0	−0.009	−0.011

After the experiments, the friction boundary of the flank face of the milling cutter tooth was identified, and the boundary of the friction boundary were obtained, as shown in Figure 19.

Figure 19.

Experimental friction boundary.

Since the friction boundary measured by the experiments is the cumulative result of all instantaneous friction boundaries in the entire cutting stroke, when calculating the friction energy consumption boundary, multiple instantaneous friction energy consumption boundaries in the same cutting stroke should be calculated. These instantaneous friction energy consumption boundaries were accumulated, and the maximum outer contour of the boundaries were taken as the final friction energy consumption boundary. The final cumulative friction energy consumption boundary is comparable to the experimentally measured boundary. The instantaneous friction energy consumption surface at multiple moments were accumulated. The cumulative friction energy consumption boundary was identified, as shown in Figure 20.

Figure 20.

The calculated boundary of friction energy consumption.

In order to compare the calculated and experimental friction boundary more intuitively, the boundary and inflection points were reconstructed in the cutter tooth coordinate, as shown in Figure 21.

Figure 21.

Friction energy consumption boundary and friction boundary.

It can be seen from Figure 21 that the friction upper boundary obtained by calculating the friction energy consumption is more similar to the measured boundary, and the calculated friction lower boundary was also relatively similar to the measured boundary, but the boundary had a higher degree of similarity. The reason for the deviation of the inflection point between the calculated and measured friction boundary mainly comes from two aspects. On the one hand, it was the influence of the simulation error on the solution result. On the other hand, the inflection points of the friction boundary were affected by the unstable vibration during the experiment.

The curves alone cannot quantitatively characterize the degree of similarity between the solution and the change characteristics of the measured boundary. The method of relative error calculation needed to be used to further verify the effectiveness of the model and method.

The relative error is defined as:

W = | (y_{p} - y_{1}) / y_{1} | \cdot 100 %

(22)

Where W is the relative error, y_p is the coordinate value of the calculated energy consumption boundary, and y₁ is the coordinate value of the measured friction boundary.

By equation (22), the maximum relative error and the minimum relative error of the calculated and measured boundary coordinate values were calculated, and the calculation results are shown in Table 9.

Table 9.

The calculation results of the relative error.

Relative error	Upper boundary (%)	Lower boundary (%)
W_max	5.71	9.71
W_min	8.93	15.85

From the results in Table 9, it can be seen that the maximum relative error between the calculated friction boundary and the actual measurement is 15.85%, the minimum relative error is 5.71%, and the average relative error is 10.05%. It shows that the calculation results of the friction energy consumption calculation model are in good agreement with the experimental results. The proposed model can reveal the evolution characteristics of the instantaneous friction energy consumption of the flank. It can realize the correct calculation of the instantaneous friction energy consumption of the flank face of the milling cutter tooth.

Conclusions

The dynamic cutting process model of the milling cutter under vibration was established. The motion trajectory and attitude angle of the milling cutter and cutter tooth were calculated in detail. The identification and characterization method of the friction boundary of the flank face were proposed. The selection method of the characteristic points in the friction area of the flank face was given.

The thermal-stress couple field simulation of the cutting process was carried out. And based on the simulated effective stress, a method to identify the feature points of the instantaneous friction boundary was proposed. The results show that the geometry and position of the friction boundary under the action of vibration are constantly changing, and the time-varying characteristics of the friction and wear behavior of the tool-workpiece interface are revealed.

A method for solving the friction speed of feature points in the friction area was proposed. The decomposition model of effective stress was established. The calculation method of friction force was given based on the decomposition model. The distribution function of instantaneous friction energy consumption on the flank face was established. The results show that the friction speed and force in the friction area of the flank face have volatility, and the value and direction of them are different with the change of position, which leads to that the distribution function of friction energy consumption is a multidimensional nonlinear function.

The identification methods of instantaneous friction speed, friction force and friction energy consumption distribution were proposed. A method to identify the evolution of friction speed, friction force and friction energy consumption was proposed. The high-efficiency milling machining experiments were carried out to verify the accuracy of the proposed model. The results show that the distribution of instantaneous friction energy consumption on the flank is high in the middle and low on both sides, and the friction energy consumption is the largest near the upper boundary. The frictional energy consumption varies with the contact angle in the same way as the frictional force. The average relative error between the calculated friction boundary and the actual measurement is 10.05%. It shows that the proposed calculation model of friction energy consumption is accurate.

Footnotes

Handling Editor: Chenhui Liang

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work are supported by National Natural Science Foundation of China (51875145) and Natural Science Foundation of Heilongjiang Province of China (ZD2020E008).

ORCID iD

Peiyi Zhao

References

Totis

Albertelli

Torta

, et al. Upgraded stability analysis of milling operations by means of advanced modeling of tooling system bending. Int J Mach Tools Manuf 2017; 113: 19–34.

García-Ordás

Alegre-Gutiérrez

González-Castro

, et al. Combining shape and contour features to improve tool wear monitoring in milling processes. Int J Prod Res 2018; 56: 3901–3913.

Xun

Ziming

Shenliang

, et al. Influence of turning tool wear on the surface integrity and anti-fatigue behavior of Ti1023. Adv Mech Eng 2021; 13: 2021.

Huang

Yang

Zhang

. Modeling and parameter optimization of NC milling process for energy consumption. China Mech Eng 2016; 27: 9. (in Chinese).

Kuram

Ozcelik

Bayramoglu

, et al. Optimization of cutting fluids and cutting parameters during end milling by using D-optimal design of experiments. J Clean Prod 2013; 42: 159–166.

Zhang

Wear optimization of titanium alloy cutter milled with microtextured ball-end milling cutter. Adv Mech Eng 2020; 12: 168781401989210.

Liu

Jiang

Gao

, et al. Investigation of the effect of vibration amplitude on the surface integrity in high-speed rotary ultrasonic elliptical machining for side milling of Ti-6Al-4V. J Mech Eng 2019; 55: 215.

Lei

Qin

Zhite

, et al. Cutting force analysis of the fixed grain tool in rotary ultrasonic machining based on tool-workpiece contact rate. Machinery 2015; 42: 1–5–+15.

Wang

Guo

Microstructural analysis of failure progression for coated carbide tools during high-speed milling of Ti-6Al-4V. Wear 2020; 456–457: 203356.

10.

Oliveira

da Silva

Baldin

, et al. Effect of tool wear on the surface integrity of Inconel 718 in face milling with cemented carbide tools. Wear 2021; 476: 203752.

11.

Chen

Bian

, et al. Investigations on Burr formation mechanisms in micro milling of high-aspect-ratio titanium alloy ti-6al-4 v structures. Int J Mech Sci 2020; 185: 105884.

12.

Saketi

Sveen

Gunnarsson

, et al. Wear of a high cBN content PCBN cutting tool during hard milling of powder metallurgy cold work tool steels. Wear 2015; 332–333: 752–761.

13.

Vipindas

Mathew

Wear behavior of TiAlN coated WC tool during micro end milling of Ti-6Al-4V and analysis of surface roughness. Wear 2019; 424–425: 165–182.

14.

Zheng

Zhao

Cheng

, et al. Frictional and wear performance of TiAlN/TiN coated tool against high-strength steel. Ceram Int 2018; 44: 6878–6885.

15.

Cai

Zou

, et al. Tool wear and machined surface characteristics in side milling Ti6Al4V under dry and supercritical CO₂ with MQL conditions. Tribol Int 2020; 151: 106511.

16.

Liu

Guo

Sealy

, et al. Energy consumption and process sustainability of hard milling with tool wear progression. J Mater Process Technol 2016; 229: 305–312.

17.

Xiao

A cutting parameter energy-saving optimization method for CNC turning batch processing considering tool wear. J Mech Eng 2021; 57: 217. (in Chinese)

18.

Tang

, et al. An investigation into methods for predicting material removal energy consumption in turning. J Clean Prod 2018; 193: 128–139.

19.

Sealy

Liu

Zhang

, et al. Energy consumption and modeling in precision hard milling. J Clean Prod 2016; 135: 1591–1601.

Instantaneous friction energy consumption and its evolution of high energy efficiency milling cutter

Abstract

Keywords

Introduction

Characterization of friction boundary of flank face

Identification of instantaneous friction boundary on flank face

Instantaneous friction energy consumption distribution function on flank face

Determination of friction speed

Determination of friction force and friction energy consumption

Distribution and evolution of friction energy consumption on flank face

Distribution of friction energy consumption

Evolution of friction energy consumption

Experimental verification

Conclusions

Footnotes

Declaration of conflicting interests

Funding

ORCID iD

References