An Effect of a Carbon-Containing Additive in the Structure of a Friction Material on Temperature of the Wet Clutch Disc

This paper consists of two parts. The first one contains a description and methodology of the composite material used as friction material in clutches. Four variants of such material, differing in the type of carbon additive (the elemental graphite, pencil graphite and foundry coke powder of various fractions) were considered. Thermal conductivity, thermal diffusivity as well as the specific heat all materials were determined experimentally. On the inertial IM-58 stand, a simulation of the braking process of the friction pair consisting of a steel disc with friction material and a counterpart in the form of a homogeneous steel disc was carried out. On this basis, averaged coefficients of friction, unchanging in the entire sliding process, were found for the four friction pairs. The experimental data obtained in the first stage were used in the second stage to develop two (2D and 3D) numerical models of the friction heating process of the friction pairs under consideration. For four variants of the friction material, a comparative spatial-temporal temperature analysis was performed using both models. It was found that a simplified axisymmetric (2D) model can be used to estimate the maximum temperature with high accuracy. The lowest maximum temperature (115.6 °C) obtained for the same total friction work was achieved on the friction surface of the material with the addition of GP-1.


Introduction
The friction materials are widely used in friction pairs of vehicles, tractors, motorcycles, machine tools and others. Such are the hydro mechanical gearboxes, the oil-cooled brakes, the clutches and etc. [1][2][3]. From the reliable operation of such units and mechanisms depends on the safety of both service personnel and surrounding people [4].
The main advantages of friction materials are high efficiency in torque transmission, low wear [5], smoothness and transparency of operation [4], absence of seizure and formation of burns [6], effective dissipation of thermal energy generated during a slip [7], possibility to minimize the force of switching on and switching off the friction device due to creation of the optimum kinematic scheme.
Known friction materials containing fibers [8], are based on a polymer matrix [9], however, for heavy-loaded friction systems, the most common are powder sintered friction materials (PSFM) [10]. The PSFMs are produced by powder metallurgy and have a complex composite structure that combines metal and non-metal components. As the basis of the friction material, iron and copper have become the most widespread. An iron is used in friction pairs without lubrication. For friction pairs with lubrication, copper-based materials are the most common. Copper allows for efficient material processing, heat dissipation, and low wear [11,12]. However, pure copper does not provide the required values of the coefficient of friction and wear intensity. To change these properties, graphite powder but also thermal stresses at a specific point on the working surface and inside the friction material were developed [25].
However, the greatest advances in temperature modeling of friction systems such as disc brakes and clutches stem from the use of numerical methods. Without listing all publications on this subject in detail, one should refer the reader to relevant reviews [26][27][28][29]. The studies that were not included in these articles cover the 2D FE computational model, implemented in ANSYS software [30]. It was assumed that both working elements of the clutch are heated separately with heat flux proportional to the friction power density. The ratio of heat distribution in the form of the Charron's formula was used as the proportionality coefficient. Single and multiple clutch operation modes were considered.
The subject of temperature is closely related to other aspects of clutch operation, in particular, such as strength of materials and mechanical properties. The analysis of the temperature field and its gradients was performed during a single engaging of the carbon fabric wet clutch in order to predict damage to the composite material [31]. The latter were defined as the increase in surface roughness and the occurrence of cracks as observed on scanning electron microscopy (SEM) micrographs.
A review and evaluation of materials used for friction elements of clutches was presented in relation to specific design solutions developed from the end of the 19th century [32]. The design aspects of the complex hybrid mechanical material system of a woven fiber yarn (glass fiber with aromatic polyamide, copper, and poly-acrylic-nitrile) reinforced friction material was performed with the intention of creating separated component groups (matrix and fiber groups). It was noted that the obtained data can be used as input parameters for the development of thermomechanical models of clutches.
The purpose of this article was: (1) development of effective 2D and 3D computational models with the use FEM to determine the temperature of the friction clutch; (2) carrying out a comparative analysis of the temperature fields obtained by means of both models with the same friction work done; (3) investigating the effect of a carbon-containing additive in the structure of the new four friction materials on the temperature mode of the clutch. (4) When developing the computational models of frictional heating of the clutch, the authors used their previous experience, gained earlier in the modeling of the temperature mode of disc brakes [33,34].

Friction Materials
The basis of the composite friction material was bronze obtained from copper powder and 10% tin (BrO 10 ), which contained 20% of the additive in the form of elemental graphite GE-1, pencil graphite GP-1 and foundry coke powder of various fractions C-1 and C-2.
The charge of the material was obtained by mixing the initial copper powders of the PMS-1 grade with an average particle size of 80 µm (Figure 1a), tin of the PO-1 grade with an average particle size of 30 µm (Figure 1b), graphite of the element grade GE-1 with an average particle size of 100 µm (Figure 1c), pencil graphite GP-1 with an average particle size of 8 µm (Figure 1d), foundry coke powder with a size of less than 60 µm (C-1) ( Figure 1e) and 160-200 µm (C-2). Mixing of the initial powders was carried out in a mixer of the "drunken" barrel type for 50 min.
It should be noted that the tribotechnical properties of the friction material are largely influenced by the amount and size of the used powder additives, which determine the contact area during friction. In addition, structure has an impact. Graphite of any grade has a layered crystalline structure, while coke is amorphous.
The used coke powder additive is obtained by crushing large pieces of coke, which is used as fuel, the cost of which is not high. Graphites of grades GE-1 and GP-1 are natural, the cost of which is 20-30% higher than the cost of foundry coke.

Experiment
The essence of the test method was to simulate the braking process on the inertial stand IM-58 and record the dependence of the friction moment on the velocity and time

Experiment
The essence of the test method was to simulate the braking process on the inertial stand IM-58 and record the dependence of the friction moment on the velocity and time of braking. The scheme of the IM-58 stand is shown in Figure 2. The tests were carried out with friction by the ends of a rotating counterpart 12 and a stationary test specimen 24, which are mounted in the mandrels 11, 22 and 23 ( Figure 2b). The axes of their rotation are located horizontally and mounted on bearings 4, 7, 10, 20. The rotating mandrel 11 is fixed on the spindle 25, which rotates from the electric motor 28 through the pulleys 1 and 2 and the V-belt transmission 3. The stationary mandrel 21 is fixed on the piston rod 15, mounted on two radial bearings and having the ability to rotate around its axis. The mandrel 21 and the rod 15 are held against rotation by a strain gauge beam. The load on the samples was created by a gear oil pump 18 powered from the oil tank 19 through the filter 17. The working oil pressure was regulated by the throttle and recorded by the pressure gauge 16. Fixation of readings of number of turns was carried out through the intermediate disk 6 and the speed sensor 27. The inertial masses 8 necessary to obtain the required kinetic energy are attached to the spindle flange 25 using bolts 9 and 26. After reaching the selected rotation velocity, the spindle 25 with the inertial masses 8 was disconnected from the electric motor using an overrunning cam clutch 5 and the stationary sample 24 was pressed with the required force against the rotating counterpart 12. With the electric motor turned off, the unique energy accumulator was the rotating inertial masses. Under the action of frictional forces between the counterpart and the test sample, the drive shaft along with the inertial masses gradually stopped, transmitting the torque to the amplifier 14 and personal computer 13.
The appearance of the test specimens is shown in Figure 3. The friction disc was a steel base in thickness 2.5 mm, made of 65G steel on which a layer of friction material with a thickness 0.5 mm was attached. The outer and inner diameters of the friction pad were 95 mm and 65 mm, respectively. The counterpart samples used for research were made of 65G steel (0.62-0.7% C, 0.17-0.37% Si, ≤0.9-1.2% Mn, ≤0.035% P, ≤0.035% S, ≤0.25% Cr, ≤0.2% Cu) with the thickness of 10.0 mm, the Brinell hardness of 180-200 MPa and the roughness of the surface of friction Ra = 0.7-0.8 μm. Hydraulic oil of grade A (TU 38.1011282-89) intended for operation in torque converters and automatic transmissions was used as a coolant. Oil was supplied from the inner part of the disc at a rate of 1.01 min −1 .
The results of the influence of the type of carbon-containing additive in the composition of the friction material on the thermal properties and the coefficient of friction are shown in Table 1. The coefficient of friction was recorded after 300 cycles of each test, i.e., in a steady state wear mode. The tests were carried out with friction by the ends of a rotating counterpart 12 and a stationary test specimen 24, which are mounted in the mandrels 11, 22 and 23 ( Figure 2b). The axes of their rotation are located horizontally and mounted on bearings 4, 7, 10, 20. The rotating mandrel 11 is fixed on the spindle 25, which rotates from the electric motor 28 through the pulleys 1 and 2 and the V-belt transmission 3. The stationary mandrel 21 is fixed on the piston rod 15, mounted on two radial bearings and having the ability to rotate around its axis. The mandrel 21 and the rod 15 are held against rotation by a strain gauge beam. The load on the samples was created by a gear oil pump 18 powered from the oil tank 19 through the filter 17. The working oil pressure was regulated by the throttle and recorded by the pressure gauge 16. Fixation of readings of number of turns was carried out through the intermediate disk 6 and the speed sensor 27. The inertial masses 8 necessary to obtain the required kinetic energy are attached to the spindle flange 25 using bolts 9 and 26. After reaching the selected rotation velocity, the spindle 25 with the inertial masses 8 was disconnected from the electric motor using an overrunning cam clutch 5 and the stationary sample 24 was pressed with the required force against the rotating counterpart 12. With the electric motor turned off, the unique energy accumulator was the rotating inertial masses. Under the action of frictional forces between the counterpart and the test sample, the drive shaft along with the inertial masses gradually stopped, transmitting the torque to the amplifier 14 and personal computer 13.
The appearance of the test specimens is shown in Figure 3. The friction disc was a steel base in thickness 2.5 mm, made of 65G steel on which a layer of friction material with a thickness 0.5 mm was attached. The outer and inner diameters of the friction pad were 95 mm and 65 mm, respectively. The counterpart samples used for research were made of 65G steel (0.62-0.7% C, 0.17-0.37% Si, ≤0.9-1.2% Mn, ≤0.035% P, ≤0.035% S, ≤0.25% Cr, ≤0.2% Cu) with the thickness of 10.0 mm, the Brinell hardness of 180-200 MPa and the roughness of the surface of friction Ra = 0.7-0.8 µm. Hydraulic oil of grade A (TU 38.1011282-89) intended for operation in torque converters and automatic transmissions was used as a coolant. Oil was supplied from the inner part of the disc at a rate of 1.01 min −1 .
The results of the influence of the type of carbon-containing additive in the composition of the friction material on the thermal properties and the coefficient of friction are shown in Table 1. The coefficient of friction was recorded after 300 cycles of each test, i.e., in a steady state wear mode.   The obtained data revealed that the highest value of the thermal conductivity of be due to the fact that the area of the metal contact of the formed tin bronze, in the case of GE-1, is significantly larger than that of GP-1 ( Figure 4a). The fine powder of graphite GP-1 in the process of mixing is able to be located both between copper particles and inside them, leading to a decrease in the area of the metal contact ( Figure 4b).
(a) (b)   The obtained data revealed that the highest value of the thermal conductivity of 44 W m −1 K −1 was obtained for friction material GE-1, containing coarse graphite. In the case of fine graphite GK-1, the thermal conductivity was only 28 W m −1 K −1 . This may be due to the fact that the area of the metal contact of the formed tin bronze, in the case of GE-1, is significantly larger than that of GP-1 ( Figure 4a). The fine powder of graphite GP-1 in the process of mixing is able to be located both between copper particles and inside them, leading to a decrease in the area of the metal contact ( Figure 4b).  The obtained data revealed that the highest value of the thermal conductivity of be due to the fact that the area of the metal contact of the formed tin bronze, in the case of GE-1, is significantly larger than that of GP-1 ( Figure 4a). The fine powder of graphite GP-1 in the process of mixing is able to be located both between copper particles and inside them, leading to a decrease in the area of the metal contact ( Figure 4b).
(a) (b)  An even greater decrease in the thermal conductivity of the friction materials C-1 and C-2 containing foundry coke powder may be due to the significantly lower value of the thermal conductivity of coke in comparison with graphite. So in the monographs [35,36] it is shown that the thermal conductivity of coke powder of various fractions is 0.249-0.542 W m −1 K −1 , while after graphitization is equal to 0.58-1.34 W m −1 K −1 . For example, the thermal conductivity of graphite determined by the spatial orientation of the layers can be 233 W m −1 K −1 , while pure copper is equal to 400 W m −1 K −1 [37].
The highest value 0.06 of the friction coefficient was obtained for the friction material C-2 with the addition of foundry coke powder, while for graphite's GE-1 and GP-1 0.45 and 0.35, respectively. This may be due to both the structure of the crystal lattice and spatial orientation. Thus, graphite's GE-1 and GP-1 have a hexagonal crystal lattice, while coke is amorphous.

Numerical Model
Numerical simulation of the friction heating process clutch during a single deceleration from the initial speed to the stop was performed on the inertial stand of the IM-58 ( Figure 2). The analyzed clutch consists of two working elements: a steel disc with friction material applied to one face, and a steel counterpart ( Figure 3). The calculations were performed for four (no. 1, 2, 3, 4) variants of the friction material with thermophysical properties and friction coefficients from Table 1. The values of the other input parameters are shown in Table 2. The given values of contact pressure, velocity, moments of inertia are characterized as parameters of heavily loaded friction units, clutches, mechanisms of automotive vehicles, special purpose vehicles, devices, machine tools. The most important dimensions of both elements are shown in Figure 5.
and C-2 containing foundry coke powder may be due to the significantly lower valu the thermal conductivity of coke in comparison with graphite. So in the monogra [35,36] it is shown that the thermal conductivity of coke powder of various fraction 0.249-0.542 W m −2 K −1 , while after graphitization is equal to 0.58-1.34 W m −2 K −1 . For ample, the thermal conductivity of graphite determined by the spatial orientation of layers can be 233 W m −2 K −1 , while pure copper is equal to 400 W m −2 K −1 [37]. The highest value 0.06 of the friction coefficient was obtained for the friction m rial C-2 with the addition of foundry coke powder, while for graphite's GE-1 and G 0.45 and 0.35, respectively. This may be due to both the structure of the crystal lattice spatial orientation. Thus, graphite's GE-1 and GP-1 have a hexagonal crystal lattice, w coke is amorphous.

Numerical Model
Numerical simulation of the friction heating process clutch during a single dece ation from the initial speed to the stop was performed on the inertial stand of the IM (Figure 2). The analyzed clutch consists of two working elements: a steel disc with tion material applied to one face, and a steel counterpart (Figure 3). The calculations w performed for four (no. 1, 2, 3, 4) variants of the friction material with thermophy properties and friction coefficients from Table 1. The values of the other input parame are shown in Table 2. The given values of contact pressure, velocity, moments of ine are characterized as parameters of heavily loaded friction units, clutches, mechanism automotive vehicles, special purpose vehicles, devices, machine tools. Wm K 600 The most important dimensions of both elements are shown in Figure 5.  Two numerical models were developed using FEM: axisymmetric (2D) and spatial (3D). The differences in these models were in the material of the lining and the disc connected to it, which had fasteners on the outer cylindrical surface limiting its rotation. In the case of the lining from the 3D model, cuts (grooves) were taken into. In the 2D model, these cuts were not included due to axial symmetry. Thus, nominal contact surfaces in 2D and 3D models were different. The final FE (quadratic Lagrange elements) meshes created for both models are shown in Figure 6, and the corresponding numbers of elements are included in Table 3. Two numerical models were developed using FEM: axisymmetric (2D) and spatial (3D). The differences in these models were in the material of the lining and the disc connected to it, which had fasteners on the outer cylindrical surface limiting its rotation. In the case of the lining from the 3D model, cuts (grooves) were taken into. In the 2D model, these cuts were not included due to axial symmetry. Thus, nominal contact surfaces in 2D and 3D models were different. The final FE (quadratic Lagrange elements) meshes created for both models are shown in Figure 6, and the corresponding numbers of elements are included in Table 3.
and the clamping force P of the clutch elements during braking was:  Taking into account the data contained in Table 2, the nominal contact surface area A a in the 2D model, determined by the values of the inner r p and outer R p radius of the disc containing the friction material, is equal to: and the clamping force P of the clutch elements during braking was: Because the friction surface area in the 3D model was smaller than in the 2D model, the contact pressure was higher. On the basis of the dimensions (Figure 5), an appropriate spatial CAD geometric model of two elements of the friction pair was developed and the friction surface area was determined A (3D) a = 3759 · 10 −6 m 2 . Then, taking into account the Formula (2), the contact pressure for the 3D model was determined: With the constant contact force P and the friction coefficients unchanged in the braking process f i (Table 1), the reduction of the angular velocity ω in time t from the initial value ω 0 to zero in the moments of stopping t = t s,i is linear: where the equivalent radius r eq was calculated from the formula [33]: Taking into account the Formulas (3) and (6) and the values of the coefficients of friction from the Formula (5) contained in Table 1, the next stopping times were determined for each of the four variants of the friction material: t s,1 = 7.02 s, t s,2 = 5.46 s, t s,3 = 4.91 s and t s,4 = 4.09 s. Due to the constant contact force of the discs P (2) in both calculation models, also the stopping times for 2D and 3D models with the use of the corresponding friction materials were the same. The change in the process of braking power Q i and work of friction W i , i = 1, 2, 3, 4: are shown in Figure 7. Based on the Formulas (8) and (9), it was established that in the moments of stopping t = t s,i , the friction work for all four calculation variants was the same and amounted to: W s,i ≡ W i (t s,i ) = 0.5W 0,i = 19.43 kJ, i = 1, 2, 3, 4. This value is equal to the initial kinetic energy of the system, calculated from the formula: The simulated process concerned frictional heating during braking in time 0 ≤ t ≤ t s,i , which made it possible to stop the rotating masses at a given moment of inertia, and then cooling the friction system during t s,i ≤ t ≤ 90 s, i = 1, 2, 3, 4. During the slip of the clutch elements resulting in generation of heat in the time interval 0 ≤ t ≤ t s,i , in the contact area, the thermal contact of friction was perfect, i.e., the sum of heat fluxes directed perpendicularly from the contact surface of the friction material and the steel disc to the inside was equal to the power density of the frictional forces, and the temperature values of the opposite surfaces were the same. At t > t s,i , i = 1, 2, 3, 4 the friction surfaces of the elements were adiabatic. Throughout the analyzed process (0 ≤ t ≤ 90 s), convection cooling on the friction-free surfaces with the constant heat transfer coefficient took place. The exception was the two outer frontal surfaces of the sliding elements, which were thermally insulated (h = 0) both during heating and cooling ( Figure 6).  (Figure 7), special tools, name ODEs (ordinary differential equations) and DAEs (differential-algebraic equati Mathematics module were incorporated.

Numerical Analysis
The numerical simulation concerned the comparative analysis of the tem fields in four friction pairs (no . 1, 2, 3, 4) obtained on the basis of the axisymm and spatial (3D) numerical models.
The temperature changes on the friction surface of the disc are shown in In accordance with the experimental tests carried out, the relative rotational movement of the clutch components was simulated. The calculations were conducted in COMSOL Multiphysics ® software by using the Heat Transfer Module [38]. Additionally, in order to determine the friction work done W i (Figure 7), special tools, namely Global ODEs (ordinary differential equations) and DAEs (differential-algebraic equations) from Mathematics module were incorporated.

Numerical Analysis
The numerical simulation concerned the comparative analysis of the temperature fields in four friction pairs (no. 1, 2, 3, 4) obtained on the basis of the axisymmetric (2D) and spatial (3D) numerical models.
The temperature changes on the friction surface of the disc are shown in Figure 8. The assumptions made in the analysis regarding the constant contact pressure and the linearly decreasing angular velocity of the rotating elements lead to a typical change in the temperature of the contact surface of the friction clutch elements-its value increases at the beginning of the process to the maximum value within approximately 2-4 s and then decreases slightly until it stops (Figure 8a). The differences in the time courses of temperature on the friction surface, obtained by means of two FE models during braking, are insignificant and occur only at the stage of cooling, after stopping (Figure 8b). The difference of the maximum temperature values obtained at the braking stage with the use of both calculation models for all friction materials does not exceed 2% (Table 4). In the cooling phase after stopping, the temperature of the friction surface of the disc drops gradually.
The temperature values at a fixed point in time obtained by means of the 2D or 3D model for all four friction materials differ negligibly. It is noticeable that the temperature found when the axisymmetric model is each time higher than that obtained with the use of the spatial model. gradually. The temperature values at a fixed point in time obtained by means of the 2D or 3D model for all four friction materials differ negligibly. It is noticeable that the temperature found when the axisymmetric model is each time higher than that obtained with the use of the spatial model.  The temperature distributions on the friction surface in the radial direction are presented in Figure 9. In the case of the axisymmetric model, when the temperature does not depend on the angular variable, the choice of the radial straight line along which the calculations were performed was arbitrary, while when using the 3D model, a line along the axis x was chosen (Figure 6b). In the contact area { ≤ ≤ 30 mm 47.5 mm r , = 0 z }, marked by vertical solid lines in Figure 9, the temperature distribution is shaped by a directly proportional linear dependence of the friction power density on the radial variable r (Figure 9a). The linear temperature rise in the contact area along the radius is the best seen in the initial braking stage at = 0.5 s t . Outside the area of heating, a clear drop in temperature is visible. The lack of heating after stopping results in a gradual leveling of temperature (Figure 9b). The highest temperature is achieved in the area of heating in the vicinity of the equivalent radius, which, according to Formula (6), is equal to = 39.4 mm   The temperature distributions on the friction surface in the radial direction are presented in Figure 9. In the case of the axisymmetric model, when the temperature does not depend on the angular variable, the choice of the radial straight line along which the calculations were performed was arbitrary, while when using the 3D model, a line along the axis x was chosen (Figure 6b). In the contact area {30 mm ≤ r ≤ 47.5 mm, z = 0}, marked by vertical solid lines in Figure 9, the temperature distribution is shaped by a directly proportional linear dependence of the friction power density on the radial variable r (Figure 9a). The linear temperature rise in the contact area along the radius is the best seen in the initial braking stage at t = 0.5 s. Outside the area of heating, a clear drop in temperature is visible. The lack of heating after stopping results in a gradual leveling of temperature (Figure 9b). The highest temperature is achieved in the area of heating in the vicinity of the equivalent radius, which, according to Formula (6), is equal to r eq = 39.4 mm. The difference between the maximum temperature values obtained with the 2D and 3D models is about 3 • C. Over time, the differences in temperature values obtained with both models increase and at the stopping time they are about 5 • C. At the time of the stop, the maximum temperature of 117.5 • C (2D model) was reached for the material no. 4 on the radius r = 39.5 mm close to the equivalent radius r eq = 39.4 mm (Figure 9a). At the end of the cooling process of the clutch friction elements, the above-mentioned temperature differences are more noticeable, but they relate to a relatively much lower temperature than during the braking stage. At the time t = 90 s, the most heated material was that denoted no. 2, and the least-material no. 3, but the difference between the temperatures of these materials was insignificant (≈ 0.9 • C).
above-mentioned temperature differences are more noticeable, but they relate to a relatively much lower temperature than during the braking stage. At the time = 90 s t , the most heated material was that denoted no. 2, and the least-material no. 3, but the difference between the temperatures of these materials was insignificant ( ≈°0.9 C ).
Temperature changes in the axial z direction for a stationary disc with the friction material ( − ≤ ≤ 3 mm 0 z ) and a rotating steel disc (counterpart ≤ ≤ 0 10mm z ) are shown in Figure 10. In the stop times = , s i t t , = 1, 2, 3, 4 i the calculation results obtained on the basis of 2D and 3D models agree well (Figure 10a). The temperature of the stationary steel disc decreases with the distance from the contact surface = 0 z . The slight increase in temperature at a deflection from the contact surface in the stationary disc is associated with the known lagging effect-the decrease in temperature inside the body is slower than on its heated surface. At the end of the process = end t t of unforced convection cooling of the clutch after stopping, the temperature distribution along the thickness of the friction pair elements is even (Figure 10b). However, as in the radial direction (Figure 9b), the differences between the temperature values obtained with the two calculation models are noticeable ( ≈°4 C ), but the temperature itself at that moment is not high. The temperature distributions on the friction surface obtained by means of the 2D and 3D models at the time moments max,i t , = 1,2,3,4 i ; when the maximum values of temperature max T were reached (Table 4) are shown in Figure 11. They show the differences between the results obtained with the use of two calculation models on one hand, and on the other they show the explicitly temperature effect of using different friction materials. For the same total friction work, the temperature of the friction surface increases with the braking time reduction ( , = .     Temperature changes in the axial z direction for a stationary disc with the friction material (−3 mm ≤ z ≤ 0) and a rotating steel disc (counterpart 0 ≤ z ≤ 10 mm) are shown in Figure 10. In the stop times t = t s,i , i = 1, 2, 3, 4 the calculation results obtained on the basis of 2D and 3D models agree well (Figure 10a). The temperature of the stationary steel disc decreases with the distance from the contact surface z = 0. The slight increase in temperature at a deflection from the contact surface in the stationary disc is associated with the known lagging effect-the decrease in temperature inside the body is slower than on its heated surface. At the end of the process t = t end of unforced convection cooling of the clutch after stopping, the temperature distribution along the thickness of the friction pair elements is even (Figure 10b). However, as in the radial direction (Figure 9b), the differences between the temperature values obtained with the two calculation models are noticeable (≈ 4 • C), but the temperature itself at that moment is not high.
The temperature distributions on the friction surface obtained by means of the 2D and 3D models at the time moments t max,i , i = 1, 2, 3, 4; when the maximum values of temperature T max were reached (Table 4) are shown in Figure 11. They show the differences between the results obtained with the use of two calculation models on one hand, and on the other they show the explicitly temperature effect of using different friction materials. For the same total friction work, the temperature of the friction surface increases with the braking time reduction (t s,1 = 7.02 s, t s,2 = 5.46 s, t s,3 = 4.91 s, t s,4 = 4.09 s). The location of the area with the increased temperature in the vicinity of the equivalent radius is visible r eq (6).
The isotherms in the cross-section of the clutch in the plane rz at the moments of reaching the maximum temperature in each of the four friction materials are shown in Figure 12. The highest value of the maximum temperature T max = 145.1 • C was observed in the case of the material no. 4, and the lowest T max = 115.6 • C was determined when using the material no. 1. The maximum temperature was reached on the radii of 40.4 mm, 40.7 mm, 40.75 mm, 40.85 mm, for the material no. 1, 2, 3, 4 respectively. It can be seen that the disc with the friction material is heated in the entire thickness of 3 mm, while the counterpart is heated only to about half of its thickness (≈5 mm).  The isotherms in the cross-section of the clutch in the plane rz at the moments of reaching the maximum temperature in each of the four friction materials are shown in Figure 12.   The isotherms in the cross-section of the clutch in the plane rz at the moments of reaching the maximum temperature in each of the four friction materials are shown in Figure 12.

Conclusions
This paper presents the results of experimental studies on the thermal and physical properties and the coefficient of friction as well as numerical simulations of the temperature mode of the new composite friction material used in the clutches. The data obtained

Conclusions
This paper presents the results of experimental studies on the thermal and physical properties and the coefficient of friction as well as numerical simulations of the temperature mode of the new composite friction material used in the clutches. The data obtained for the four variants (no. 1, 2, 3, 4) of such material, differing in the type of carbon additive, were analyzed. The axisymmetric (2D) and the spatial (3D) numerical models of the wet clutch in order to determine the transient temperature fields of the friction discs during a single engagement and after the stop were developed. Based on the obtained results, it was found that: (1) the highest thermal conductivity and specific heat has the material no. 2, while the lowest values of the quantities has the material no. 3; (2) the highest coefficient of friction appears for the steel disc combined with the friction material no. 4, and the lowest for steel disc and friction material no. 1; (3) the estimation of the maximum clutch temperature can be carried out with sufficient accuracy using a 2D model, which allows for the reduction of labor losses at the preparatory stage and computational time. The single simulation case carried out on the workstation with CPU Intel ® Xeon ® E5-2698 v4 @ 2.20GHz; RAM 64 GB (DDR4) lasted approximately 40 s, and 2200 s, when using 2D and 3D models, respectively. On the other hand, the determination of the temperature field in the elements of the friction clutch is better to carry out with the use of the 3D model; (4) At the same total friction work done during a single clutch engagement, the lowest temperature was achieved when using friction material no. 1, and the highest for material no. 4. Also, the use of these materials resulted in the longest and shortest periods of frictional sliding, respectively.
It should be noted that in the proposed numerical models with the use of FEM, the averaged, experimentally obtained, thermophysical properties of the composite material were used. The development of the computational models that take into account the structure of the composite and the properties of individual components belongs to the future. Subsequent studies will aimed at using the developed 2D and 3D models for assessing the effect of additives of powders (ceramic, metallic, intermetallic compounds), their particle size, distribution and the amount in the composition of the friction material on the thermal effects during friction. One of the important aspects would be to include realistic properties of a coolant. This will make it possible to assess the efficiency of their use without expensive studies on the full scale inertial test bench.

Conflicts of Interest:
The authors declare no conflict of interest.