Next Article in Journal
Tribological Performance and Rheological Properties of Engine Oil with Graphene Nano-Additives
Next Article in Special Issue
Numerical Wear Modeling in the Mixed and Boundary Lubrication Regime
Previous Article in Journal
Super-Hard DLC Coatings as an Alternative to Polycrystalline Diamond for Cutting Tools: Predictive Analysis of Aluminium Alloy Surface Quality
Previous Article in Special Issue
A Closer Look at the Contact Conditions of a Block-on-Flat Wear Experiment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Micro-Wear Simulation of Braking Interface Based on Particle Discrete Element Modeling

College of Mechanical Engineering, Dalian Jiaotong University, Dalian 116028, China
*
Author to whom correspondence should be addressed.
Submission received: 28 March 2022 / Revised: 8 June 2022 / Accepted: 24 June 2022 / Published: 26 June 2022
(This article belongs to the Special Issue Advances in Wear Predictive Models)

Abstract

:
For material fracture and severe wear in braking conditions, the discrete element method (DEM) is used to simulate the wear process of the braking interface explicitly. Based on the central difference method, particle motion equations are established considering the influence of elemental damping on particle contact. Combined with the Particle Flow Code (PFC) software, a DEM wear model of the braking interface is established using the parallel bond modeling method. The braking wear process is simulated, and the material damage process is investigated. The simulation results demonstrate that with the increase of the initial braking load and the initial braking speed, the wear depth increased by 24.75% and 16.22%, respectively. The increase in the number of detached particles leads to an increasing trend of fracture force chains, which increases the thickness of the flowing particle layer, revealing the micro-wear mechanism of the braking interface.

1. Introduction

When a high-power disc brake is subjected to high-speed and heavy-load braking conditions, the wear phenomenon of the braking interface is severe. The brake disc and pad surface could have varied damage. Therefore, there is no doubt that the wear characteristics of the braking interface should be studied.
Under braking conditions, wear and fracture as a universal phenomenon affect braking efficiency. The form of material wear is based on the detachment of wear debris, including blocks and particles. The continuous mechanical approaches usually simulate the mechanical behavior of the material. The finite element method (FEM) is difficult to use for the real-time observation of the wear process and estimating the wear amount [1]. However, simulation of non-continuous phenomena such as multi-fracturing is not well-adapted to a continuous description. In this case, the discrete element method (DEM) is an excellent alternative to consider discontinuities. DEM could simulate the entire friction process more realistically, which has obvious advantages in studying material structure destruction and describes a wear phenomenon of the braking interface. Many researchers have shown interest in this approach for wear and fracture simulation.
In the early years, DEM was proposed by Cundall et al. to separate a non-continuum into rigid elements, which solved the motion equations of each element by iterative methods and then obtained the motion morphology of the non-continuum [2]. Euser et al. found that DEM could simulate inter-particle connection bond breaking and crack formation, which makes DEM applications more and more widely used in cutting and processing [3].
The braking interface could produce the third-body and other wear debris during the braking process. In practice, the classical DEM modeling approach is a challenge considering that the flows of the third-body and the wear debris depend on the mechanical properties of the material. Consequently, Renouf et al. considered the mechanical properties of the third-body and its thermal and physicochemical properties, proposed an extended DEM, and applied it to the simulation of third-body flows [4]. It was found to be an effective numerical method for calculating engineering problems involving particle materials. However, the representation of irregular particles using DEM is a challenging issue. Irazábal et al. presented a new insight in applying one of those simplifications known as rolling friction, which avoided excessive rotation when irregularly shaped materials were simulated as spherical particles [5].
Similarly, to overcome excessive rolling of such an oversimplified shape of the particle, Aboul et al. considered contact rolling resistance. Meanwhile, the influences of the elastic and the plastic parameters were discussed. The results showed that the plastic macroscopic behavior of the particle assembly only depended on the plastic parameters, which mainly depended on the plastic rolling moment reflecting the shape of the particle [6].
Particle–particle interactions play an important role in particle flows. Huang et al. found that the flows could have some strong nonlinear characteristics and be different from normal materials, either solids or liquids [7]. Huang et al. used a developed 2D coupled method to analyze the mechanical response of rocks. Compared with the theoretical approach, the displacement by numerical simulation corresponds with an analytical equation for cases with different ground stress ratios [8]. These suggest that DEM could be used to model a continuum, including concrete or rocks. However, it requires a significant calibration effort, even to capture just the linear elastic behavior of a continuous model via the classical force-displacement relationships at the contact interfaces between particles. Celigueta et al. proposed a new way for computing the contact force between discrete particles. The newly proposed force considered the surroundings of the contact, not just the contact itself, which provided an accurate approximation to an elastic continuum and avoided calibration of DEM parameters for the purely linear elastic range [9].
Fracture of the force chain indicates that the material structure is destroyed. Vandenende et al. investigated the effects of inter-particle friction and a distributed grain size on macroscopic strain rates. They found that the discrete element model could reproduce realistic compaction behavior. The strain rates produced by the model corresponded with uniaxial compaction experiments [10]. Furthermore, the mechanical response of a particle system is influenced by bulk material properties and various factors. Desu et al. investigated the influences of the factors as mentioned above on the macroscopic stress-strain response, which helped understand macroscopic damage in terms of microscopic factors [11].
With the continuous development of virtual simulation technology, the application range of DEM has become large. DEM could be used to simulate wear prediction. Zha et al. introduced different constitutive laws for contacts to represent the stiffness and strength characteristics of the metal panels and the core materials. The approach could predict detailed local damage and failure modes, an advantage over the conventional numerical methods [12]. Fillot et al. presented a fundamental approach to studying the wear process by considering the detachment of particles. This numerical model was constructed to visualize and accurately measure the influence of the interfacial layer on the wear process [13]. In addition, Chaudry et al. found that the maximum tensile stress-based criterion could be used to address computational aspects of crushing/comminution of particle materials [14]. Damien et al. defined the materials by Young’s modulus, Poisson’s ratio, and density. DEM could be used to fit the static and dynamic mechanical behavior of the material [15]. Phan et al. demonstrated that it could study abrasive roll wear [16]. Boemer et al. suggested a generic procedure that predicted the wear evaluation of the lining surfaces [17].
To ensure the accuracy and stability of the DEM wear prediction, one of the essential tasks is to perform sensitivity analysis at the single-particle level. Chen et al. conducted the pin-on-disc wear tests modeled to predict the sliding wear by individual iron ore particles. The reference demonstrated that accurate and stable results of sliding wear volume could be achieved. Increasing particle density or radius caused more sliding wear at a steady state. By contrast, particle Poisson’s ratio, particle shear modulus, geometry mesh size, rotating speed, coefficient of restitution, and time step had no impact on wear volume [18].
Meanwhile, Hossein et al. researched the abrasive wear of iron ore pellets by numerical simulation using DEM and an analytical model. The variable parameters of the analysis model, including average normal contact force and sliding distance, were predicted from DEM simulation results [19].
Moreover, Ke et al. investigated the influences of different normal loads on the seismic moment, macroscopic friction coefficient, kinetic energy, gouge layer thickness, and recurrence time between slips. The simulations revealed the influences of normal loads on the dynamics of sheared particle fault gouges and demonstrated the capabilities of DEM for studying the stick-slip dynamic behavior of the particle fault system [20]. In summary, it can be seen that DEM has certain advantages for the wear study of the material.
This paper mainly focuses on the wear process of the braking interface to carry out DEM wear simulation. Combined with elemental damping, the dynamic relaxation method is used to establish the contact relationship between discrete particles. The periodic boundary simulates the dynamic motion of the brake disc and pad. Considering the high-speed and heavy-load conditions, the DEM wear model of the braking interface is established to analyze the dynamic wear process. Based on the simulation results, the velocity, displacement, and force chains of discrete particles can describe the surface wear and fracture process. The effect of initial braking load and initial braking speed on surface wear is discussed. Meanwhile, the simulation results are verified by pin-on-disc wear tests.

2. Materials and Methods

2.1. Particle Motion Analysis

DEM is a method of discrete and numerical simulation of continuous media. Figure 1 is the flow chart of the iterative solution within the DEM wear simulation, whose details are as follows:
(1)
The DEM wear model should be established by adding geometry particles and selecting the contact model.
(2)
The contact particles can be identified, and the contact force is calculated.
(3)
The force-displacement criterion determines the magnitude of the particle force so that particles automatically contact and separate according to the wear criterion of the material, and the wear of the material particles is realized.
(4)
According to Newton’s second law, acceleration and displacement can be calculated.
(5)
The DEM wear model is remapped to different arid points.
(6)
All particle messages can be updated, and the new potential contacts of the new particle system can be checked.
(7)
The computer performs an iterative solution until the dynamic process of the DEM iteration simulation is completed.
To analyze the metal deformation mechanism using DEM, the metal material model needs to be discretized. The motion model of the metal material is discretized into the motion model of the particle unit. During the DEM wear simulation process, the collision force between particles could cause each particle to translate and rotate. According to Newton’s second law, the external force produces acceleration, so particles gain motion speed and produce displacement. The change of position causes the contact relationship between particles to change. The new displacement increment results in a new external force and particles generate a new displacement and contact relationship until the system is stable. In the force model of a single degree of freedom system, the differential equation of motion is usually used to describe [21]. Based on the dynamic relaxation method, the contact characteristics of the discrete particle element are considered:
{ m u ¨ t + η u ˙ t + λ u t = F t I φ ¨ t + η φ ˙ t + λ φ t = M t
where, Ft is the external force, η is the spring damping coefficient, λ is the elastic coefficient, I is the rotational inertia, u is the displacement vector, φ is the rotation angle vector, and Mt is the resultant moment.
To solve the motion process of a single particle, the motion equation is established, the target solution is sought, and the motion state is analyzed. Figure 2 shows the motion state of a single particle from position Ot to position Ott in the period time Δt, and the motion equation under the action of the combined external force can be expressed:
v t + Δ t = v t + a t Δ t
where, vt is the velocity vector at time t, vtt is the velocity vector at time tt, and at is the acceleration vector at time t.
The motion equation of the particle is determined by the form of translational motion and rotational motion, and the central difference method is used to convert the problem of solving the differential equation system into the issue of solving the algebraic equation system [22].
(1)
Calculation of particle translational motion
The acceleration and velocity in Equation (1) are expressed by the displacement vector u, as shown in the following equation:
{ u ˙ t = u t + Δ t u t Δ t 2 Δ t u ¨ t = u t + Δ t 2 u t + u t Δ t ( Δ t ) 2
Substituting Equation (1) into Equation (3), the displacement vector utt at time t + Δt can be obtained:
u t + Δ t = ( Δ t ) 2 F t + ( η 2 Δ t m ) u t Δ t + [ 2 m k ( Δ t ) 2 ] u t m + η 2 Δ t
The displacement increment Δu can be obtained:
Δ u = ( Δ t ) 2 m + η 2 Δ t F t + η 2 Δ t m m + η 2 Δ t u t Δ t + m η 2 Δ t k ( Δ t ) 2 m + η 2 Δ t u t
Substituting Equation (4) into Equation (3), the acceleration vector u ¨ t and velocity vector u ˙ t can be obtained.
(2)
Calculation of particle rotational motion
The angular acceleration and angular velocity in Equation (1) are expressed by angle vector φ, as shown in the following equation:
{ φ ˙ t = φ t + Δ t φ t Δ t 2 Δ t φ ¨ t = φ t + Δ t 2 φ t + φ t + Δ t ( Δ t ) 2
Substituting Equation (1) into Equation (6), the angular vector φtt at time tt can be obtained:
φ t + Δ t = ( Δ t ) 2 M t + ( η 2 Δ t I ) φ t Δ t + [ 2 I k ( Δ t ) 2 ] φ t I + η 2 Δ t
The angular increment Δφ can be obtained:
Δ φ = ( Δ t ) 2 I + η 2 Δ t M t + η 2 Δ t I I + η 2 Δ t φ t Δ t + I η 2 Δ t k ( Δ t ) 2 I + η 2 Δ t φ t
Substitute Equation (7) into Equation (6), the angular acceleration vector φ ¨ t and angular velocity vector φ ˙ t can be obtained.
The particle motion equation is composed of Equation (3) to Equation (8). Furthermore, the DEM wear simulation can calculate the motion equation of the particle to determine the motion form and realize the discrete elemental motion simulation.

2.2. Wear Criterion

The essence of wear is when part of the material is separated from the matrix material. From the perspective of the DEM, the wear process can be explained as the bond of material particles breaking, causing particles or particle clusters to detach from the matrix material and form wear debris. A parallel bond model describes the connection between particles, as shown in Figure 3. The force and moment generated by particles are related to the maximum normal stress and tangential stress under the boundary conditions [23]. Stiffness coefficients and damping coefficients are added to the contact position to simulate collision and squeeze state of particles. If the value of any stress is greater than the corresponding bond strength, it means that the state of the parallel bond is broken. The parallel bond model has five parameter definitions: the normal stiffness λn, the tangential stiffness λs, the normal damping ηn, the tangential damping ηs, and the bond radius R.
The resultant force FP in the parallel bond model can be decomposed into two parts, the normal force Fn and the tangential force Fs:
F P = F n + F s
The resultant moment MP of the contact cross-section can be decomposed into two parts: the normal moment Mn and the tangential moment Ms:
M P = M n + M s
The maximum normal stress σmax and maximum tangential stress τmax at the contact cross-section are:
{ σ max = | F n | A + | M s | I R τ max = | F s | A + 2 | M n | I R
If the maximum stress value of the particles is greater than the compressive strength (σmaxσc) or shear strength (τmaxτc), the parallel bond of the connected particles is broken, and the values of the contact force and moment is zero. The separation phenomenon occurs between particles, and particles are separated from the matrix material. That is, the process can be judged as material wear.

3. DEM Wear Simulation of Braking Interface

3.1. Parallel Bond Model

To establish the DEM wear model of the braking interface, the model needs to be simplified appropriately. According to the braking mechanism and wear theory, the wear area of the brake disc and pad is mainly concentrated on the contact surface. Therefore, the wear behavior between the brake disc and pad can be analyzed by constructing a DEM wear model of the braking interface. The following critical assumptions made in particle-based DEM wear simulations can be stated. Elemental particles are rigid, have mass and rotational inertia, and can be analytically described.
(1)
The motion of each particle is independent, including the translational motion and the rotational motion.
(2)
The contact between particles occurs over a tiny area, and each contact involves only two particles.
(3)
Particles can overlap slightly at the contact points, which can transmit tensile and compressive forces in the normal contact direction and a tangential force orthogonal to the normal contact force.
(4)
Considering the separation distance between two particles, the tensile force can be calculated. When the tensile force exceeds the maximum tension force at the contact point, particles could move away from each other, and the contact is no longer considered in calculating the contact force.
The discrete element properties of the braking interface need to be determined by the parallel bond model method, whose details are as follows:
(1)
Discrete particles are generated in the specified area based on the micro-parameters of the known model, and each particle has its mechanical properties by assigning mechanical parameters.
(2)
Pressure is applied to particles to generate isotropic force. To ensure the compactness of the particle model, floating particles with less than three pairs of contact pairs need to be deleted.
(3)
Parallel bonds are generated, wall constraints are removed, and a parallel bond model with material properties can be obtained, as shown in Figure 4.
(4)
The FISH function is used to release the stress of the particle model, which can avoid the influence of the internal stress generated in the modeling process on the simulation.
The DEM model requires assumptions regarding aggregate geometry and grain size distribution to solve the contact stress. The force chain represents the contact mode of particles, and the types of force chain can be divided into tensile force chain and compressive force chain. Furthermore, the force chain thickness can describe the force chain strength, as shown in Figure 5. In the entire DEM system, the strong force chain is responsible for transferring large external force, supporting most of the weight and the external load of the particle system, but the number of strong force chains is small. The weak force chain is responsible for transferring small external force, but the number of weak force chains is huge. Meanwhile, the weak force chain has an auxiliary effect on the stability of the strong force chain and participates in reconstructing the strong force chain after a fracture. Therefore, the fracture of the force chain can characterize the failure of the material structure.
The micro-parameters of the DEM model materials cannot be directly obtained through tests, so the correct micro-parameter settings affect the accuracy of the DEM wear simulation. To make the established DEM wear model have mechanical properties that match the actual material, the macro-mechanical model is used to adjust the micro-parameters of the DEM wear model. When the numerical values of the simulation results match the actual material mechanics test values, the micro-parameter settings in the model are considered reasonable, and the following simulation work can be carried out.
The macroscopic material properties of the friction pair are shown in Table 1. The mechanical properties of materials can be calibrated by a uniaxial compression test [24]. Uniaxial compression tests are performed based on the micro-parameters, as shown in Figure 6. (1) Model establishment: The height value of particle model is 0.1 mm, the width value of particle model is 0.05 mm, and the particle radius is 3.6 × 10−3 mm. The upper wall and lower wall move to the center at a speed of 0.1 m/s, and a small amount of damping 0.1 is applied. The test is run until stress falls below 70% of the peak, and the stress-strain curve of the particle model is recorded during the uniaxial compression test. (2) Micro-parameter calibration: Based on the stress-strain curve results of the uniaxial compression test, the Young’s modulus and Poisson’s ratio are calculated with the stress peak as the calibration point:
{ E = Δ σ y Δ ε y ν = Δ ε x Δ ε y
where E is Young’s modulus, ν is Poisson’s ratio. σy is the stress, and εx and εy are the strains in the x and y directions, respectively.
The calibration tests of the DEM model are performed using the DEM micro-parameters. When the mechanical properties of the DEM model are close to the values in Table 1 and within the allowable error range of Table 2, it is considered that the DEM micro-parameters can describe the macroscopic mechanical properties of the brake disc and pad materials. The calibration tests show that the relative error of Young’s modulus on the brake disc is 1.47%, the relative error of Poisson’s ratio on the brake disc is 2.58%, the relative error of Young’s modulus on the brake pad is 1.67%, and the relative error of Poisson’s ratio on the brake pad is 2.00%. Therefore, the revised results for the micro-parameters follow the calibrated material properties and the micro-parameters of the discrete element model can be obtained, as shown in Table 3.

3.2. Establishment of DEM Model

Based on the braking mechanism and wear theory, it is known that the wear area of the brake disc and pad mainly concentrates on the braking interface. Therefore, the wear behavior is analyzed by constructing a DEM wear model of the braking interface. The wear phenomenon is the cumulative result of a long time. When the wear process enters the stable wear stage, the wear trends of the braking interface are unchanged. In the process of DEM analysis, the braking interface wear is simplified. When establishing the DEM wear model of the braking interface, the PFC software’s unique periodic boundary is used to simplify the continuous contact model into a periodic flow model. Using periodic boundaries can simulate huge particle sets by considering only selected subdomains (called periodic cells).
In the periodic space, discrete particles and particle clusters constitute the main components of the material composition, which would tend to move outside the boundaries due to particle motion, as shown in Figure 7. The simulation region is then ‘remapped’ so that particles re-enter the simulation region at the corresponding location on the opposite periodic face. The criterion for this remapping considers the particle centroid coordinates. When particle B and particle cluster A–C are detected to move outside the boundary wall, they are introduced in the boundary wall at the same elevation as the new particle B’ and particle cluster A’–C’.
The local wear interface of the brake disc and pad is selected. During the rotation of the brake disc, the change of radian in the local area does not change significantly. Therefore, it is assumed that the local contact between the brake disc and pad is linear relative motion, and the DEM wear model of the braking interface is established, as shown in Figure 8.
The upper particle area of the model is the brake pad, and the lower particle area is the brake disc. The model boundary conditions are applied by adding constraints on the particle boundary wall, mainly consisting of three parts: the pressure wall, the support wall, and the periodic domain wall. The pressure wall of the brake pad provides the load for the discrete model and ensures the contact state of the brake disc and pad by applying the nodal pressure in the -y-direction. The pressure wall of the brake disc provides support for the brake disc material and restricts the degree of freedom in the y-direction. The support wall supports the brake pad and restricts the degree of freedom in the x-direction. The periodic domain wall sets the particle circulation area. A velocity vector in the -x-direction of particles is applied to simulate the motion of the brake disc. The speed loaded in the model is the linear speed of the brake disc.
In addition, changes in temperature can cause expansion and compression of materials. The thermal stress increment associated with free expansion corresponding to temperature increment ΔT has the form:
Δ σ i j = E α t Δ T δ i j
where Δσij is the thermal stress increment, αt is the coefficient of liner thermal expansion, and δij is the Kronecker delta [25]. Furthermore, thermal stress affects the value of contact stress.
Therefore, heat sources need to be added to the DEM wear model. The temperature curve is extracted from the FEM simulation results in the literature [26], as shown in Figure 9. Because the discrete element research scale is small and the temperature change of the overall model is little, the unified heat sources are set as the boundary. Heat transfers only when a parallel bond connection forms between the particles; when the bond is broken, heat cannot be transferred.
High-power braking conditions can be simulated based on the control variable method. During the loading process of DEM wear simulation boundary conditions, initial braking load and initial braking speed need to be converted into pressure and relative speed. In addition, based on PFC discrete element software simulation, the DEM analysis program is written by particle flow code, and the wear criterion is calculated to analyze the microscopic wear process of the braking interface. During the discrete element simulation, the braking simulation conditions are shown in Table 4. It is assumed that the braking process is stable and the braking conditions are unchanged. The braking time is set to 8 s.

3.3. Pin-on-Disc Wear Test

A pin-on-disc wear testing machine was used for pin-on-disc wear tests. The friction disc sample adopted Q345B, and the pin sample adopted copper-based powder metallurgy; the material properties are shown in Table 1. The friction load was 60 N, and the friction speed was 1200 r/min, as shown in Figure 10. Before tests, the friction surface of the sample was polished to make the surface roughness Ra = 0.3 μm. During testing, the room temperature was 20 °C, and the relative humidity was 25%~50%. After tests, the friction surface was cleaned with acetone solution, and the solution was filtered and dried to obtain wear debris. An optical microscope measured the size of the wear debris, and the friction surface morphology was observed by an optical three-dimensional surface profile measuring instrument. Then the friction state and the wear degree of the friction pair were judged.

3.4. DEM Simulation Analysis

Based on braking condition 5 in Table 2, the single braking DEM wear simulation is carried out. According to the wear criterion, when the maximum stress value exceeds its ultimate strength, the material would be damaged and destroyed, which leads to surface wear. In the DEM wear simulation results, the wear debris exists in the form of particles and particle clusters. The vertical displacement of the particle indicates the wear depth of the surface. The range of the wear layer includes the damaged surface of the brake disc and pad, as shown in Figure 11. Figure 11a is the particle displacement cloud diagram in the initial simulation. The maximum wear depth is 3.66 × 10−4 mm. The braking process begins, and the wear area concentrates on the contact surface of the brake disc and pad. The range of the wear layer is small, and the wear depth is low. Figure 11b is the particle displacement cloud diagram in the end simulation. The maximum wear depth of the surface is 7.69 × 10−2 mm. The range of the wear layer becomes large, and the wear depth is high. In addition, the detached material forms the wear debris in the wear layer and moves away from the surface.
According to the damage range, the wear range ratio of the brake disc to the brake pad is 1:3. However, only a small part of the brake disc has material wear, and the wear area only occurs on the brake pad surface. The brake pad surface is severely worn, and the wear depth is large. This is because the difference in the properties of friction materials causes the actual wear of the brake pad to be greater than the wear of the brake disc. The mechanical properties of the structural alloy steel used in the brake disc are high. The powder metallurgy used in the brake pads is prone to material damage, and the brake pad is an easy wear part. Therefore, during a single braking process, the wear of the brake disc is often ignored, and the wear of the brake pad is mainly used to analyze the surface wear law.
To study surface wear during braking, the maximum displacement of particles under different braking time steps is extracted, as shown in Figure 12. Friction and wear, cumulative damage, and local fracture of the surface are observed. Figure 12a shows the contact state of the brake disc and pad. Because the brake disc and pad have the surface characteristics of pits and asperities, the interface contacts each other and scratches in the initial braking, as shown in Figure 12b. Figure 12c further illustrates that the interaction of the asperity with the rough surface aggravates the surface wear. This is because the strength of the brake disc material is higher than that of the brake pad, and the brake pad is the first to suffer material damage. Figure 12d shows asperities are destroyed, and the material falling off the braking surface turns into wear debris. A part of the wear debris is filled in pits of the rough surface, and the other part leaves the wear layer of the surface with the brake disc movement. When moving in the wear layer, the area of the asperity is stacked with wear debris, as shown in Figure 12e. Currently, the extrusion phenomenon increases the surface wear. When the wear debris and asperity are destroyed, the surface forms pits and the range of the wear layer is further expanded, as shown in Figure 12f.
During the braking process of the disc brake, the surface material undergoes elastoplastic deformation. In severe wear cases, the material undergoes shear damage, and the material is separated from the brake disc and pad for transfer, causing surface wear. The velocity vector diagram describes the motion trend of metal materials, and the wear process is analyzed, as shown in Figure 13.
When the brake disc and pad contact, the brake pad moves downward under the initial braking load. The rotation motion of the brake disc causes its moving speed direction to be horizontal, as shown in Figure 13a. During the relative motion of the brake disc and pad, the surface material is squeezed. When the tangential force and the normal force exceed the allowable value of the material, shear fracture occurs in part of the material, which leaves the surface and moves in the braking interface, as shown in Figure 13b,c. However, the brake disc and pad are in a state of squeezing each other, and the detached material accumulates, forming blocks or clusters of wear debris, as shown in Figure 13d.
Furthermore, a large amount of wear debris forms a material transfer layer at the braking interface. Through the action of the shear force and the extrusion force, the wear debris accumulated in the material transfer layer is destroyed again, forming discrete wear debris with a smaller volume, as shown in Figure 13e. Discrete wear debris leaves the surface, and the empty positions on the surface are supplemented by new discrete materials. The material transfer layer is always in a dynamic partial detaching and re-repairing process, as shown in Figure 13f. However, the damage to the surface cannot be repaired, leading to the continuous expansion of the wear layer, so the surface of the brake disc and pad experiences material wear.

4. Discussion

In order to compare the feasibility of the DEM wear model, the change trend of the average wear depth is compared with the FEM simulation results, and the change trend of the surface wear morphology is compared with the pin-on-disc test results as well.

4.1. Comparative Analysis of Average Wear Depth

Based on the FEM wear simulation results obtained under the same simulation conditions [26], the cubic B-spline curves of the average wear depth of the FEM and DEM simulation results are drawn by means of data fitting, as shown in Figure 14. Within the same time step, the growth trends of the wear curves of the FEM and DEM simulation results are basically the same. It can be seen that the average wear depth and time are positively correlated. Due to the uneven wear plane of the material matrix, there will be slight errors.

4.2. Comparative Analysis of Test Morphology

The DEM wear simulation results show that the wear form of the material on the friction pair surface changes during the friction process. Therefore, the results of the DEM analysis are verified by the results of the pin-on-disc wear tests. As shown in Figure 15, the shape and size of the wear debris present different characteristics, which affect the change of the surface topography. Figure 15a shows that when the friction time is 2 min, the size of the wear debris generated is less than 10 μm. In the initial wear test, the structure of the sample surface is rough, and the asperity on the friction surface is mutually squeezed. When the stress exceeds the yield limit, the friction material is damaged, and the asperity is destroyed to produce small wear debris. The material composition is the same as that of the friction pair. As shown in Figure 15b, when the soft structure of the material is worn away, hard particles scratch the surface, resulting in apparent scratches in local areas of the surface.
As the test progresses, the size of the wear debris increases, as shown in Figure 15a. When the friction time is 5 min, the size of the wear debris is 30 μm, and the wear debris exists in two forms, including spherical wear debris and flake wear debris. Spherical wear debris is due to the mechanical mixing action of the transfer material of the friction pair, and the wear debris is accessible on the braking interface. The accumulation of a large amount of small wear debris forms flake wear debris, and the contact area between the flake wear debris and the friction surface is large, causing surface damage, as shown in Figure 15b. Under the action of the normal force, the meshing effect of metal materials is strengthened, and the local damage to the braking interface gradually expands to large-area surface damage. Under the action of the shear force, the material transfer occurs on the friction surface. The friction material falls off the surface structure in the form of massive wear debris, and the number of pits and scratches formed on the sample surface increases. Therefore, test results verify the change process of the wear form during the DEM wear simulation process.

5. Effect of Braking Parameters on Brake Disc Wear

5.1. Effect of Initial Braking Load on Brake Disc Wear

In the process of DEM wear simulation, the change in the initial braking load affects the force state of the material. The particle displacement results of the surface in simulation braking conditions 1, 2, and 5 are extracted, respectively, as shown in Figure 16. The influence mechanism of the initial braking load on surface wear is analyzed. From the simulation results, the material failure trend is distributed in layers, the wear debris leaves the surface during the braking process, and the interspace is generated. However, the interspace variation is minimally affected by macroscopic surface contact. When the relative position of the brake disc and pad remains unchanged, the increase in surface clearance will change the surface wear state.
Figure 16a is a cloud diagram of particle displacement when the initial braking load is 13 kN. The wear phenomenon occurs on the surface of the brake disc and pad, leading to material damage and destruction. The detached material leaves the surface as wear debris, resulting in an interspace in the braking interface. When the initial braking load is 17 kN, the wear becomes severe due to the increase of the load, and the wear debris of the surface clearance increases, as shown in Figure 16b. The reason is that when the rough surfaces are squeezed against each other, the large load makes the asperity easier to destroy, and the wear debris that does not leave the surface accumulates. When the initial braking load is 21 kN, the wear layer of the surface becomes significant, as shown in Figure 16c. There is an apparent fault phenomenon between the wear debris and the surface, and the interspace becomes large, which further proves that the increase of the load would aggravate surface wear. Similarly, it shows that surface wear results from the interaction between the rough surfaces and is affected by the motion of the wear debris on the surface.
Initial braking load can directly affect the number of particle force chain breaks, which cause material particles to fall off, as shown in Figure 17. In the DEM wear simulation, the force chain distribution trend of the particle indicates the severity of the material damage. Black compressive force chains are more, and red tensile force chains are less. During the braking process, the organizational damage of the material is mainly affected by the squeezing force.
When the initial braking load is 13 kN, strong force chains of the primary materials of the brake disc and pad are densely distributed, but strong force chains of the wear area are sparsely distributed, as shown in Figure 17a. It shows that the load of the primary body material is within the allowable range, but the material structure is unstable. There is a small volume of loose wear debris, and the wear phenomenon only occurs on the surface. When the initial braking load is 17 kN, the number of strong force chains increases, but the number of strong force chains between the brake disc and the main body of the brake pad is significantly reduced. The increased initial braking load causes some strong force chains to break with maximum stress values exceeding their strength limits.
Meanwhile, decreased strong force chains reduce the stability of the particle system, resulting in a sparser distribution trend of force chains, as shown in Figure 17b. The material structure is damaged, and wear debris increases. Under the action of a large load, the strength of the surface material is weakened, causing material damage to occur more easily, and aggravating surface wear. When the initial braking load is 21 kN, the number of strong force chains between braking interfaces is small. The strong force chain structure of the brake disc and brake pad body close to the surface is loose, and the strong force chain structure far away from the surface is dense, which shows that the damage range of the particle system increases, as shown in Figure 17c. The surface damage extends to the inside of the material, and the surface wear continues to expand with the material damage.
In addition, Figure 16 shows the changes in the number of force chains of particles and clearly distinguishes the range of force chains between the brake disc and the main body of the brake pad. The strength limit of the brake disc material is high, but the strength limit of the brake pad material is low. Therefore, in the braking process, the brake pad material undergoes shear failure first under the action of the squeezing force. The surface wear of the brake pad is greater than that of the brake disc.
To quantitatively analyze surface wear, the wear depth of the surface and the number of force chains of the material organization are extracted, as shown in Figure 18. It can be seen from Figure 18a that the increase in the initial braking load causes surface wear to increase. Braking load affects the formation of the surface material transfer layer. With the increase in the initial braking load, the number of wear debris in the material transfer layer increases to a steady state. It means that the initial braking load is large within a certain range, which causes serious damage to the surface material, and more wear debris is generated on the surface. By counting the number of force chains of particles, it can be found that as the initial braking load increases, the number of force chains of particles shows a downward trend, as shown in Figure 18b. The large load causes damage and destruction of the material, resulting in a large number of force chains breaking, which reduces the structural strength of the surface material. This is because the initial braking load is large, and the normal force on the surface is large. The material is severely squeezed, and the material reaches the limit of compressive strength early. Meanwhile, the damage degree of material tissue is severe.
In addition, when the initial braking load rises from 13 kN to 17 kN, the wear depth increases by 13.14%, and the number of force chains decreases by 23.37%. When the initial braking load rises from 17 kN to 21 kN, the wear depth increases by 11.61%, and the number of force chains decreases by 5.96%. It can be seen that the initial braking load rises uniformly, and the increase in the wear depth is similar. However, the decrease rate of the break number of force chains is large when the initial braking load is 17 kN, indicating that the material structure is unstable and wear is prone to occur. Therefore, in the actual braking process, it is recommended to select an initial braking load of 13 kN for mechanical braking under the premise of ensuring braking demand.

5.2. Effect of Initial Braking Speed on Brake Disc Wear

The initial braking speed determines the braking speed of the brake disc during braking. The simulation braking conditions 3, 4, and 5 are extracted, and the influence mechanism of the initial braking speed on surface wear is analyzed, as shown in Figure 19.
It can be seen from Figure 19a that when the initial braking speed is 800 r/min, the particle distribution is not uniform. Still, the damage of the particle system is concentrated on the surface. When the initial braking speed is 1000 r/min, the damage range of the particle system is enlarged, and the number of discrete particles in the wear layer increases, as shown in Figure 19b. The initial braking speed increases and the mutual impact of rough surfaces is strengthened. The wear debris impacts the surface under the rotation influence of the brake disc, causing damage to the asperity on the surface. When the initial braking speed is 1200 r/min, it can be found that surface wear is even. The motion speed of particles is accelerated at a high rotation speed, and surface particles fall off by impact force easily, as shown in Figure 19c. As the initial braking speed increases, the motion speed of particles in the wear layer accelerates, and the particles do not have enough time to fill the interspace, so faults occur during wear.
The initial braking speed affects the motion state of particles and changes the shear force generated by the relative motion. The relative motion speed of the brake disc and pad is large, and the surface material is subjected to a large shear force. Under high-speed braking conditions, the initial braking speed increases the mutual extrusion of the braking surface and the collision between the wear debris and the surface, causing local deformation and plastic fracture on the surface. In addition, a higher relative motion speed causes the wear debris to leave the braking interface faster, and thereby, surface wear is accelerated.
When the initial braking load is the same, the force chain cloud diagram can reveal the influence of the initial braking speed on the surface damage degree, as shown in Figure 20. When the initial braking speed is 800 r/min, the strong force chain structure of the brake disc and the main body of the brake pad is relatively loose, and the weak force chain assists the strong force chain supporting material structure shown in Figure 20a. When the initial braking speed is 1000 r/min, the number of strong force chains increases significantly, and the structure of the force chain between the brake disc and the main body of the brake pad is compact. The reason is that the initial braking speed is large, the energy carried by particles is high, and strong force chains are needed to carry the load. When the maximum stress value of the force chain exceeds the strength limit of the strong force chain, strong force chains will break, and parts of weak force chains are transformed into strong force chains as a supplement to maintain the material structure, as shown in Figure 20b.
Under the same load, the increase of the initial braking speed increases the speed difference between the relative motion of the brake disc and pad. A large amount of kinetic energy is converted into internal energy, leading to changes in material properties and accelerating the surface wear. When the initial braking speed is 1200 r/min, it can be seen from Figure 20c that the force chain distribution in the wear area is sparse, and the strong force chain distribution between the brake disc and the main body of the brake pad is dense. The reason is that the initial braking speed increases, and the compressive and tensile forces of particles increase, accelerating the formation of strong force chains and making the particle system more stable. It shows that increasing the initial braking speed increases the shear force of the material, making it easier for the shear stress of the material to reach the shear strength, resulting in material failure.
The wear depth and the number of force chains at different initial braking speeds are extracted, as shown in Figure 21. Increased initial braking speed aggravates the wear of the braking interface, as shown in Figure 21a. The growth of the initial braking speed increases the friction force and impact force, increasing the mechanical engagement of the rough surfaces. Meanwhile, the process of converting kinetic energy into internal energy releases a large amount of heat, which reduces the hardness and strength of the material. The action of shear force increases the damage range of the surface. As shown in Figure 21b, the initial braking speed is high and the particle shear stress is large, causing the fracture of strong force chains in the particle system. However, due to the decrease in the number of weak force chains, the structural stability of the particle system decreases, and damage continues to occur. It shows that increased initial braking speed increases the squeezing force of the material, reducing the surface material strength. Furthermore, the interaction with the braking interface becomes intense, and the damage range of the material increases, which eventually leads to severe wear of the braking interface.
Moreover, when the initial braking speed rises from 800 r/min to 1000 r/min, the wear depth increases by 7.45%, and the number of force chains decreases by 15.58%. When the initial braking speed rises from 1000 r/min to 1200 r/min, the wear depth increases by 8.77%, and the number of force chains decreases by 3.86%. It can be seen that the wear depth rises with increasing initial braking speed. However, when the initial braking speed is 1000 r/min, the number of force chain fractures decreases greatly, indicating that increased initial braking speed would increase the damage of the material structure and aggravate the surface wear. This phenomenon can affect the contact state of the surface and reduce the braking capacity of the disc brake. Therefore, in the actual braking process, selecting an initial braking speed of 800 r/min for mechanical braking under the premise of ensuring braking demand is recommended.

6. Conclusions

Based on the braking conditions of the high-power disc brake, the braking interface of the DEM model was established and analyzed, which showed a dynamic process of surface wear. The wear depth and force chains of particles were examined with numerical simulations. It can be concluded that braking parameters play an important role in braking wear. According to these investigations, a list of conclusions can be drawn as follows:
(1)
A wear phenomenon of the braking interface is caused by the continuous destruction of the material, and the wear value is positively correlated with time.
(2)
Through the comparison and analysis of simulation and tests, it can be seen that the surface material transfer changes dynamically during the braking process. The wear debris exists in the form of particle clusters or discrete particles. A part of the wear debris detached with the action of the friction pair and the other part is embedded in the surface.
(3)
A large initial braking load would decrease the surface material strength and the surface interspace would increase. Therefore, the increased motion range of the wear debris results in a fracture between the wear debris and the surface. As the initial braking load rises from 13 kN to 21 kN, the wear depth increases by 24.75%.
(4)
Under high-speed working conditions, the large centrifugal force makes the wear debris easily leave the surface. The violent squeezing action causes the shear stress to increase, and local deformation of the surface and plastic fracture of the material occur similarly. As the initial braking speed rises from 800 r/min to 1200 r/min, the wear depth increases by 16.22%.

Author Contributions

Conceptualization, Z.S. and Q.H.; methodology, Z.S. and Q.H.; software, Z.S. and Q.H.; validation, Q.H. and Y.L.; formal analysis, Q.H. and Y.L.; investigation, Q.H. and F.M.; resources, Q.H. and Y.L.; data curation, Q.H. and Y.L.; writing—original draft preparation, Z.S. and Q.H.; writing—review and editing, Q.H. and Y.L.; visualization, Y.L. and J.Y.; supervision, Z.S. and Y.L.; project administration, Y.L. and S.Z.; funding acquisition, Z.S. and J.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 51675075), Scientific Research Project of the Education Department of Liaoning Province (No. LJKZ0479), and Dalian Science and Technology Innovation Fund (No. 2021JJ12GX009).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schmidt, A.A.; Schmidt, T.; Grabherr, O.; Bartel, D. Transient wear simulation based on three-dimensional finite element analysis for a dry running tilted shaft-bushing bearing. Wear 2018, 408–409, 171–179. [Google Scholar] [CrossRef]
  2. Cundall, P.; Marti, J.; Beresford, P.; Last, N.; Asgian, M. Computer modeling of jointed rock masses. In Technical Report; Army Engineer Waterways Experiment Station: Vicksburg, MS, USA, 1978. [Google Scholar]
  3. Euser, B.; Rougier, E.; Lei, Z.; Knight, E.E.; Frash, L.P.; Carey, J.W.; Viswanathan, H.; Munjiza, A. Simulation of fracture coalescence in granite via the combined finite-discrete element method. Rock Mech. Rock Eng. 2019, 52, 3213–3227. [Google Scholar] [CrossRef] [Green Version]
  4. Renouf, M.; Cao, H.; Nhu, V. Multiphysical modeling of third-body rheology. Tribol. Int. 2010, 44, 417–425. [Google Scholar] [CrossRef] [Green Version]
  5. Irazábal, J.; Salazar, F.; Oñate, E. Numerical modelling of granular materials with spherical discrete particles and the bounded rolling friction model. Application to railway ballast. Comput. Geotech. 2017, 85, 220–229. [Google Scholar] [CrossRef]
  6. Hosn, R.A.; Sibille, L.; Benahmed, N.; Chareyre, B. Discrete numerical modeling of loose soil with spherical particles and interparticle rolling friction. Granul. Matter 2017, 19, 4. [Google Scholar] [CrossRef] [Green Version]
  7. Huang, Y.J.; Nydal, O.J.; Ge, C.; Yao, B. An introduction to discrete element method: A meso-scale mechanism analysis of granular flow. J. Dispers. Sci. Technol. 2015, 36, 1370–1377. [Google Scholar] [CrossRef]
  8. Huang, F.; Wang, Y.; Wen, Y.; Lin, Z.; Zhu, H. The deformation and failure analysis of rock mass around tunnel by coupling finite difference method and discrete element method. Indian Geotech. J. 2019, 49, 421–436. [Google Scholar] [CrossRef]
  9. Celigueta, M.A.; Latorre, S.; Arrufat, F.; Oñate, E. Accurate modelling of the elastic behavior of a continuum with the Discrete Element Method. Comput. Mech. 2017, 60, 1460–1480. [Google Scholar] [CrossRef]
  10. Ende, M.P.A.V.D.; Marketos, G.; Niemeijer, A.R.; Spiers, C.J. Investigating compaction by intergranular pressure solution using the discrete element method. J. Geophys. Res. Solid Earth 2018, 123, 107–124. [Google Scholar] [PubMed] [Green Version]
  11. Desu, R.K.; Gan, Y.; Kamlah, M.; Annabattula, R.K. Mechanics of binary crushable granular assembly through discrete element method. Nucl. Mater. Energy 2016, 9, 237–241. [Google Scholar] [CrossRef] [Green Version]
  12. Zha, X.; Wan, C.; Fan, Y.; Ye, J. Discrete element modeling of metal skinned sandwich composite panel subjected to uniform load. Comput. Mater. Sci. 2013, 69, 73–80. [Google Scholar] [CrossRef]
  13. Fillot, N.; Iordanoff, I.; Berthier, Y. Modelling third body flows with a Discrete Element Method a tool for understanding wear with adhesive particles. Tribol. Int. 2007, 40, 973–981. [Google Scholar] [CrossRef]
  14. Chaudry, M.; Wriggers, P. On the computational aspects of comminution in discrete element method. Comput. Part. Mech. 2017, 5, 175–189. [Google Scholar] [CrossRef]
  15. André, D.; Iordanoff, I.; Charles, J.-L.; Néauport, J. Discrete element method to simulate continuous material by using the cohesive beam model. Comput. Methods Appl. Mech. Eng. 2012, 213, 113–125. [Google Scholar] [CrossRef] [Green Version]
  16. Phan, H.; Tieu, A.; Zhu, H.; Kosasih, B.; Zhu, Q.; Grima, A.; Ta, T. A study of abrasive wear on high speed steel surface in hot rolling by Discrete Element Method. Tribol. Int. 2017, 110, 66–76. [Google Scholar] [CrossRef]
  17. Boemer, D.; Ponthot, J. A generic wear prediction procedure based on the discrete element method for ball mill liners in the cement industry. Miner. Eng. 2017, 109, 55–79. [Google Scholar] [CrossRef]
  18. Chen, G.; Schott, D.L.; Lodewijks, G. Sensitivity analysis of DEM prediction for sliding wear by single iron ore particle. Eng. Comput. 2017, 34, 2031–2053. [Google Scholar] [CrossRef] [Green Version]
  19. Hossein, M.; Masoumi, A.; Mehdi, V. Analysis of abrasive damage of iron ore pellets. Powder Technol. 2018, 331, 20–27. [Google Scholar]
  20. Gao, K.; Euser, B.J.; Rougier, E.; Guyer, R.A.; Lei, Z.; Knight, E.E.; Carmeliet, J.; Johnson, P.A. Modeling of stick-slip behavior in sheared granular fault gouge using the combined finite-discrete element method. J. Geophys. Res. Solid Earth 2018, 123, 5774–5792. [Google Scholar] [CrossRef] [Green Version]
  21. Thomson, W.; Marie, D. Theory of Vibration with Applications, 5th ed.; Tsinghua University Press: Beijing, China, 2005. [Google Scholar]
  22. Rapaport, D. The Art of Molecular Dynamics Simulation; Cambridge University Press: New York, NY, USA, 2004. [Google Scholar]
  23. Gang, W.; Rui, W.; Meng-meng, W.; Peng-fei, W.; Yu-Ming, Z. Simulation of conventional triaxial test on coal under hydro-mechanical coupling by particle flow code. Rock Soil Mech. 2016, 37, 537–546. [Google Scholar]
  24. Wang, S.; Wang, S.; Xu, Q.; Chen, H. Crushed model and uniaxial compression analysis of random packed ceramic pebble bed by DEM. Fusion Eng. Des. 2018, 128, 53–57. [Google Scholar] [CrossRef]
  25. Brandt, F.T.; Das, A.; Frenkel, J.; Taylor, J.C. Kronecker delta energy terms in thermal field theory. Phys. Rev. D Part. Fields 2002, 65, 085008. [Google Scholar] [CrossRef] [Green Version]
  26. Zhang, S.; Hao, Q.; Liu, Y.; Jin, L.; Ma, F.; Sha, Z.; Yang, D. Simulation study on friction and wear law of brake pad in high-power disc brake. Math. Probl. Eng. 2019, 2019, 6250694. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flow chart of iterative solution within DEM wear simulation.
Figure 1. Flow chart of iterative solution within DEM wear simulation.
Lubricants 10 00136 g001
Figure 2. Illustration of particle motion.
Figure 2. Illustration of particle motion.
Lubricants 10 00136 g002
Figure 3. Illustration of particle-particle contact force geometry.
Figure 3. Illustration of particle-particle contact force geometry.
Lubricants 10 00136 g003
Figure 4. Generation process of parallel bond model.
Figure 4. Generation process of parallel bond model.
Lubricants 10 00136 g004
Figure 5. Composition of force chain within DEM.
Figure 5. Composition of force chain within DEM.
Lubricants 10 00136 g005
Figure 6. Uniaxial compression tests.
Figure 6. Uniaxial compression tests.
Lubricants 10 00136 g006
Figure 7. Consideration of boundaries in periodic elements.
Figure 7. Consideration of boundaries in periodic elements.
Lubricants 10 00136 g007
Figure 8. Diagram of braking interface model.
Figure 8. Diagram of braking interface model.
Lubricants 10 00136 g008
Figure 9. Temperature curve of braking interface.
Figure 9. Temperature curve of braking interface.
Lubricants 10 00136 g009
Figure 10. Pin-on-disc wear test.
Figure 10. Pin-on-disc wear test.
Lubricants 10 00136 g010
Figure 11. Displacement of braking interface particles. (a) Braking time, 0.1 s; (b) Braking time, 8 s.
Figure 11. Displacement of braking interface particles. (a) Braking time, 0.1 s; (b) Braking time, 8 s.
Lubricants 10 00136 g011
Figure 12. Displacement of particles at different times. (a) Braking time, 0.1 s (average wear depth is 3.66 × 10−4 mm); (b) braking time, 0.5 s (average wear depth is 5.42 × 10−3 mm); (c) braking time, 1.0 s (average wear depth is 1.16 × 10−2 mm); (d) braking time, 2.0 s (average wear depth is 2.15 × 10−2 mm); (e) braking time, 4.0 s (average wear depth is 4.15 × 10−2 mm); (f) braking time, 8.0 s (average wear depth is 7.69 × 10-2 mm).
Figure 12. Displacement of particles at different times. (a) Braking time, 0.1 s (average wear depth is 3.66 × 10−4 mm); (b) braking time, 0.5 s (average wear depth is 5.42 × 10−3 mm); (c) braking time, 1.0 s (average wear depth is 1.16 × 10−2 mm); (d) braking time, 2.0 s (average wear depth is 2.15 × 10−2 mm); (e) braking time, 4.0 s (average wear depth is 4.15 × 10−2 mm); (f) braking time, 8.0 s (average wear depth is 7.69 × 10-2 mm).
Lubricants 10 00136 g012
Figure 13. Velocity vector change of braking interface particles. (a) Contact of interfacial material; (b) relative sliding of interface; (c) elastic deformation of material; (d) formation of wear debris; (e) repair of material transfer layer; (f) formation of transferred particle layer.
Figure 13. Velocity vector change of braking interface particles. (a) Contact of interfacial material; (b) relative sliding of interface; (c) elastic deformation of material; (d) formation of wear debris; (e) repair of material transfer layer; (f) formation of transferred particle layer.
Lubricants 10 00136 g013
Figure 14. Comparison curves of the average wear depth between the FEM and DEM simulation.
Figure 14. Comparison curves of the average wear depth between the FEM and DEM simulation.
Lubricants 10 00136 g014
Figure 15. Results of pin-on-disc wear test. (a) Observation of surface wear debris by electron microscope; (b) observation of surface topography by optical microscopy.
Figure 15. Results of pin-on-disc wear test. (a) Observation of surface wear debris by electron microscope; (b) observation of surface topography by optical microscopy.
Lubricants 10 00136 g015
Figure 16. Distribution of particle displacement with different initial braking loads. (a) Initial braking load is 13 kN; (b) initial braking load is 17 kN; (c) initial braking load is 21 kN.
Figure 16. Distribution of particle displacement with different initial braking loads. (a) Initial braking load is 13 kN; (b) initial braking load is 17 kN; (c) initial braking load is 21 kN.
Lubricants 10 00136 g016
Figure 17. Distribution of strong force chains with different initial braking loads. (a) Initial braking load is 13 kN; (b) initial braking load is 17 kN; (c) initial braking load is 21 kN.
Figure 17. Distribution of strong force chains with different initial braking loads. (a) Initial braking load is 13 kN; (b) initial braking load is 17 kN; (c) initial braking load is 21 kN.
Lubricants 10 00136 g017
Figure 18. Effect of initial braking loads on braking interface wear. (a) Wear depth; (b) force chain.
Figure 18. Effect of initial braking loads on braking interface wear. (a) Wear depth; (b) force chain.
Lubricants 10 00136 g018
Figure 19. Distribution of particle displacement with different initial braking speeds. (a) Initial braking speed is 800 r/min; (b) initial braking speed is 1000 r/min; (c) initial braking speed is 1200 r/min.
Figure 19. Distribution of particle displacement with different initial braking speeds. (a) Initial braking speed is 800 r/min; (b) initial braking speed is 1000 r/min; (c) initial braking speed is 1200 r/min.
Lubricants 10 00136 g019
Figure 20. Distribution of strong force chains with different initial braking speeds. (a) Initial braking speed is 800 r/min; (b) initial braking speed is 1000 r/min; (c) initial braking speed is 1200 r/min.
Figure 20. Distribution of strong force chains with different initial braking speeds. (a) Initial braking speed is 800 r/min; (b) initial braking speed is 1000 r/min; (c) initial braking speed is 1200 r/min.
Lubricants 10 00136 g020
Figure 21. Effect of initial braking speeds on braking interface wear. (a) Wear depth; (b) force chain.
Figure 21. Effect of initial braking speeds on braking interface wear. (a) Wear depth; (b) force chain.
Lubricants 10 00136 g021
Table 1. Material parameters of brake disc and pad.
Table 1. Material parameters of brake disc and pad.
Material PropertiesBrake Disc
(Q345B)
Brake Pad
(Copper-Based Powder
Metallurgy)
Density/(kg/m3)78505250
Young’s modulus/(GPa)201178
Poisson’s ratio0.3020.294
Table 2. Mechanical properties of brake disc and pad.
Table 2. Mechanical properties of brake disc and pad.
Mechanical PropertiesBrake DiscBrake Pad
Calibrated ResultsRelative ErrorCalibrated ResultsRelative Error
Young’s modulus204 GPa1.47%180 GPa1.67%
Poisson’s ratio0.32.58%0.312.00%
Table 3. DEM micro-parameters of brake disc and pad.
Table 3. DEM micro-parameters of brake disc and pad.
Micro-Parameters of PBM ModelBrake DiscBrake Pad
Specify the effective Young’s modulus/(MPa)9.35 × 1048.45 × 104
Normal-to-shear stiffness ratio1.751.75
Tensile strength/(MPa)440245
Cohesion/(MPa)440245
Linear expansion coefficient9.5 × 10−61.11 × 10−5
Normal critical ratio0.50.5
Friction coefficient0.50.5
Radius/(mm)3.6 × 10−33.6 × 10−3
Bond gap/(mm)5 × 10−65 × 10−6
Table 4. Simulation braking conditions.
Table 4. Simulation braking conditions.
Braking
Conditions
Initial Braking
Load/(kN)
Simulation Pressure/(MPa)Initial Braking
Speed/(r/min)
Simulation Relative Speed/(m/s)
1130.33120047.12
2170.43120047.12
3210.5380031.41
4210.53100039.27
5210.53120047.12
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sha, Z.; Hao, Q.; Yin, J.; Ma, F.; Liu, Y.; Zhang, S. Micro-Wear Simulation of Braking Interface Based on Particle Discrete Element Modeling. Lubricants 2022, 10, 136. https://0-doi-org.brum.beds.ac.uk/10.3390/lubricants10070136

AMA Style

Sha Z, Hao Q, Yin J, Ma F, Liu Y, Zhang S. Micro-Wear Simulation of Braking Interface Based on Particle Discrete Element Modeling. Lubricants. 2022; 10(7):136. https://0-doi-org.brum.beds.ac.uk/10.3390/lubricants10070136

Chicago/Turabian Style

Sha, Zhihua, Qiang Hao, Jian Yin, Fujian Ma, Yu Liu, and Shengfang Zhang. 2022. "Micro-Wear Simulation of Braking Interface Based on Particle Discrete Element Modeling" Lubricants 10, no. 7: 136. https://0-doi-org.brum.beds.ac.uk/10.3390/lubricants10070136

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