Next Article in Journal
Some Finite Difference Methods to Model Biofilm Growth and Decay: Classical and Non-Standard
Previous Article in Journal
Correlation Effects in Trimeric Acylphloroglucinols
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical and Numerical Solutions for the Thermal Problem in a Friction Clutch System

by
Laith A. Sabri
1,
Katarzyna Topczewska
2,
Muhsin Jaber Jweeg
3,
Oday I. Abdullah
4,5,6,* and
Azher M. Abed
7
1
Department of Mechatronics, Al-Khwarizmi College of Engineering, Baghdad 10003, Iraq
2
Faculty of Mechanical Engineering, Bialystok University of Technology, 15-351 Bialystok, Poland
3
College of Technical Engineering, Al-Farahidi University, Baghdad 10005, Iraq
4
Department of Energy Engineering, College of Engineering, University of Baghdad, Baghdad 10003, Iraq
5
Department of Mechanics, Al-Farabi Kazakh National University, Almaty 050040, Kazakhstan
6
System Technologies and Engineering Design Methodology, Hamburg University of Technology, 21073 Hamburg, Germany
7
Air Conditioning and Refrigeration Techniques Engineering Department, Al-Mustaqbal University College, Hillah 51001, Iraq
*
Author to whom correspondence should be addressed.
Submission received: 9 October 2021 / Revised: 6 November 2021 / Accepted: 11 November 2021 / Published: 15 November 2021
(This article belongs to the Section Computational Engineering)

Abstract

:
The dry friction clutch is an important part in vehicles, which has more than one function, but the most important function is to connect and disconnect the engine (driving part) with driven parts. This work presents a developed numerical solution applying a finite element technique in order to obtain results with high precision. A new three-dimensional model of a single-disc clutch operating in dry conditions was built from scratch. As the new model represents the real friction clutch including all details, the complexity in the geometry of the clutch is considered one of the difficulties that the researchers faced using the numerical solution. The thermal behaviour of the friction clutch during the slip phase was studied. Meanwhile, in the second part of this work, the transient thermal equations were derived from scratch to find the analytical solution for the thermal problem of a clutch disc in order to verify the numerical results. It was found, after comparison of the numerical results with analytical results, that the results of the numerical model are very accurate and the difference between them does not exceed 1%.

1. Introduction

Nowadays, owing to the increase of demands for high-quality industrial products, along with paramount consideration for lowering cost and weight, induce many attempts in the design and development of elastic systems. Friction clutches constitute vital components in many applications such as vehicles, machines of productions, etc. A popular known application of a clutch in an automotive vehicle is used to connect and disconnect the engine with the gear box. An accurate and reliable stress analysis of the clutch system parts is one of the major concerns to the mechanical design engineers to ensure reliability and long-life operation. However, an accurate computation of thermal stresses and contact pressure is very difficult to obtain due to many complexities involved in this process such as the geometry of the clutch system, variation of material properties, coefficient of friction, surface roughness, etc. The high temperatures will appear on the contacting surfaces of the clutch’s elements, which lead to generate of high thermal stresses. These thermal stresses are responsible for several disadvantages such as surface cracks, permanent distortions, etc. Finally, in some cases these disadvantages may lead the surfaces of contact to fail before the expected lifetime. Figure 1 illustrates the main parts of the dry clutch system (clutch disc, pressure plate, and flywheel) that are responsible to transfer the power from the driving part to the driven part.
After studying the available literature that related to the numerical solutions of the thermal problem of dry friction clutches through the sliding/heating phase, it was found that the vast majority of researchers used two-dimensional numerical models (axisymmetric models) to simulate the thermal problem of the clutch system to obtain the thermal behaviour and performance during the thermal stage [1,2,3,4,5,6,7,8,9]. Meanwhile, there are a very limited number of studies that developed three-dimensional numerical models, where most of these models failed to represent all the complexities of the geometry of a frictional clutch and they were approximated models. However, they obtained acceptable results compared with the analytical and experimental solutions.
On the other hand, a limited number of researchers resorted to an analytical solution to find the temperatures and thermal stresses of the automotive friction clutches and brakes based on some assumptions in an attempt to reduce the complications of the thermal problem of the clutch. One of these hypotheses is assuming that the variation in temperature is in one dimension only and neglecting many other important factors such as convection, surface roughness, distribution of contact pressure, etc. [10,11,12,13,14,15,16,17,18,19,20,21]. The purpose of all these hypotheses is to obtain an analytical solution by avoiding the solutions for 2D or 3D models.
Chao et al. [22] applied the finite element method to study the thermal behaviour of the dry dual clutch under different load conditions. A 3D model was created using CATIA software and then the ABAQUS software was used for the numerical analysis. Their results presented the temperature and the thermal stresses during the sliding time. They applied the thermal load based on the basic approach, that adopted the heat partition factor. Its mean was to find the solution of each element of the clutch system separately. Using the same approach, Bao et al. [23] investigated the effect of grooves in the frictional clutch disc in order to study the change in the thermal behaviour due to the shape of the groove. Three types of grooves were used, which were a two-way parallel groove, a three-way parallel groove, and a waffle groove. They found that the the waffle groove shape is the most optimal one because of the minimum level of heat generated.
Zhang et al. [24] presented a 2D numerical model of the wet clutch system. It was included in the numerical model the effect of the viscous shear, asperity friction and convection effect. The Runge–Kutta method was applied to find the solution of the frictional thermal model. They investigated the influence of some significant factors such as material parameters, applied pressure, initial angular velocity, friction lining permeability, and surface roughness. They found that the surface roughness and the material properties have a great effect on the surface temperature.
The researchers were faced with challenges of cost and complexities of the clutch system when they resorted to the experimental work. Furthermore, the time issues to design, perform and validate a new test rig. Therefore, it was concluded that the most appropriate approach among the solutions which are currently available based on the reasons mentioned above to find the solution of the frictional thermal problem of automotive clutches and brakes is the finite element method. During the last years, the researchers overcame great challenges to build an advanced mathematical model, taking into account many important factors that interact during realistic operating conditions.
The main objective of this work is to present a new developed three-dimensional FE model of a dry friction clutch and adopt it to compute the distributions of temperature during the slipping period (single engagement). This model is considered the most accurate model that is available until now, based on the available literature. Where all shape complexities were taken into consideration for all elements of the clutch system (Figure 1, Figure 2 and Figure 3). Furthermore, it was derived a new mathematical model to the clutch system from scratch to obtain the thermal analytical solution in order to simulate the thermal problem during the sliding time. The comparisons were made between the results of two solutions, where the maximum difference between them is not exceeded 1%.

2. Formulation of Thermal Problem

During the heating phase of the clutch’s working cycle, slipping will occur between driving and driven parts. As result, high temperature will appear on the contacting surfaces of the flywheel, the clutch disc and the pressure plate. Where, it was supposed that all amounts of heat energy were all generated from all kinetic energy, which means no losses. Thus, the frictional heat generated at any time between the surfaces of contacting elements is [9],
q t ( r , t ) = μ   p   r ω s
where, ts is the slipping time, ωs is the angular sliding velocity (rad/s), p is the contact pressure (N/m2), μ is coefficient of friction and r is the disc radius (m). Assuming that the angular slip rate decreases as following from,
ω s ( t ) = ω o ( 1 1 t s )
ωo: the preliminary angular sliding velocity (ts = 0). Therefore, the final equation for the rate of heat generated between the contacting elements is,
q c ( r , t ) = μ   p   r   ω o ( 1 t t s )
All the parameters of the clutch’s elements have been indexed with (c) for the friction material, (cu) for the axial cushion, (p) for the pressure plate and (f) for the flywheel. The high frictional heat will be generated as a result of the sliding velocity with the existing of applied pressure between the contacting surfaces, where the frictional material has high roughness. The heat will dissipate between the elements of the clutch system by conduction and to the working environmental by convection. Hence, that the slipping time is too short (less than 1 s), therefore it can be ignored the influence of radiation in the numerical analysis with no significant effect on the accuracy of results.
It can be considered that the creating of the 3D geometry of the main parts of the clutch system that including in the numerical solution is the fundamental stage. Solidworks 2020 was used to create the model of the clutch system based on the realistic dimensions. After that, a numerical approach was developed to simulate the thermal problem of the clutch system using Ansys (workbench) 2021 in order to obtain the complete solution of the problem. It was assumed the temperature at the initial time of sliding was,
T ( r , θ , z , 0 ) = T i
Next step, the heat conduction equation was applied for the main elements of the frictional clutch as shown in Figure 2,
2 T r 2 + 1 r T r + 1 r 2 2 T θ 2 + 2 T z 2 = 1 α T t
where, α is the thermal diffusivity (α = k/ρ c). r, θ and z are the cylindrical coordinates. It was taken into consideration the convection effect on the surfaces that exposed to the working environment during the sliding time. Figure 3 illustrates the thermal load that generates between the surfaces of contacting elements.

3. Finite Element Analysis

The main steps of the numerical approach (FE) will present in this section, including the hypotheses required to obtain the numerical solution. The thermal problem in the friction clutch system should be represented as a transient problem due to the large variation of temperature during a very short time. The solution of the thermal problem of the clutch system based on the finite element method can be written as [25],
[ C ] { T } + [ K ] { T } = { F }
where, [C] and [K] are the matrices of specific heat and conductivity. {T}, {} and {} are the vectors of the temperatures of nodes, derivative of temperature with respect to time and thermal load, respectively. One of the significant factors affect the accuracy for finite element results is the selection mesh of the FE model. Where, it was made the Mesh Independence Study was conducted in order to obtain the most suitable mesh to obtain highly accurate results of temperatures. The Crank–Nicolson technique was applied to the scheme. It was assumed that friction clutch works under dry conditions and the heating phase is very short. The boundary and initial conditions of the elements for dry friction clutch system (single disc) are given as follows (Figure 2 and Figure 3):
  • Axial cushion: the condition at the outer disc radius,
    K c u T r | r = r 0 c u = h [ T ( r o c u , z , t ) T a ] ,                         0 θ 2 π ,   ( t c u / 2 ) z ( t c u / 2 ) ,   t 0
    and at the inner disc radius is,
    K c u T r | r = r i c u = 0 ;               0 θ 2 π ,   ( t c u / 2 ) z ( t c u / 2 ) ,   t 0
    where Ta is the ambient surrounding temperature and h is the convection heat transfer coefficient.
  • Friction clutch disc: the condition at the inner disc radius; there is a convection,
    K c T r | r = r i c = h [ T ( r i c , z , t ) T a ] ,         0 θ 2 π ,   [ ( t c u / 2 ) + t c ] z [ ( t c u / 2 ) + t c ] ,   t 0
    and at the outer disc radius is,
    K c T r | r = r o c = h [ T ( r o c , z , t ) T a ] ,                         0 θ 2 π ,   [ ( t c u / 2 ) + t c ] z ( [ t c u / 2 ) + t c ] ,   t 0
    the frictional heat generated on the contact surfaces of clutch disc is,
    K c T r | z = t c + ( t c u / 2 ) = q c ( r , t ) ,           r i c r r o c ,     0 θ 2 π ,   0 t t s
    K c T r | z = [ t c + ( t c u / 2 ) ] = q c ( r , t ) ,           r i c r r o c ,       0 θ 2 π ,   0 t t s
  • Flywheel: the conditions on the inner and outer disc radii are,
    K f T r | r = r i f = h [ T ( r i f , z , t ) T a ] ,         0 θ 2 π ,   [ ( t c u / 2 ) + t c ] z [ t f + ( t c u / 2 ) + t c ] ,   t 0
    K f T r | r = r o f = h [ T ( r o f , z , t ) T a ] ,         0 θ 2 π ,   [ ( t c u / 2 ) + t c ] z [ t f + ( t c u / 2 ) + t c ] ,   t 0
    and at the back side of flywheel as,
    K f T r | z = [ t f + ( t c u / 2 ) + t c ] = h [ T ( r , z , t ) T a ] ,           r i f r r o f ,     0 θ 2 π ,   0 t t s
  • Pressure plate: the conditions at the inner and outer radii are:
    K p T r | r = r i p = h [ T ( r i p , z , t ) T a ] ,         0 θ 2 π ,   [ ( t c u / 2 ) + t c ] z [ t p + ( t c u / 2 ) + t c ] ,   t 0
    K p T r | r = r o p = h [ T ( r o p , z , t ) T a ] ,         0 θ 2 π ,   [ ( t c u / 2 ) + t c ] z [ t p + ( t c u / 2 ) + t c ] ,   t 0
  • The initial temperature is,
    T ( r , θ , z , 0 ) = T i
There are two theories for designing the friction clutch system which are the uniform pressure and uniform wear. In this numerical analysis, the uniform wear theory was adopted, which is considered the fundamental step to compute the magnitude of frictional heat generated. The materials that are selected in the numerical analysis are isotropic and homogeneous. All details of the inputs of the numerical simulation (materials, working conditions and dimensions) are listed in Table 1. The selected type of the frictional material is sintered metal (dry), where the range of maximum pressure is (2–2.75 MPa), the maximum allowable temperature for continuous and instantaneous cases are 622 and 822 K, and the maximum allowable velocity is 1097.28 m/s [26].
In order to find the optimal mesh of the friction clutch system, it was achieved the mesh independence test, as shown in Figure 4, where the selected finite element model of the clutch system has 268,675 elements (1,164,367 nodes). Figure 5 shows the Flowchart of the developed numerical approach (FE) to find the solution of transient thermal problem. A 3D Finite Element model (Figure 6a) of a complete frictional clutch system was built as the first step, where it could be used to study unsymmetric load, wear and hotspot. In this analysis, the objective is to validate the results of the new numerical model with analytical solutions to be reliable for any future investigations. According to the uniform wear assumption, the frictional thermal load is uniform over the contacting area. Owing to this assumption, a part/sector of the complete FE model (1/8 of complete model) will be sufficient to simulate the thermal problem in the dry friction clutch (Figure 6b). The selected thermal element is SOLID90 with twenty nodes and the DOF is the temperature. The reason to select this kind of element is the complex geometry of the clutch system that has many boundaries with curves.
The main development in this analysis is to extend the numerical solution from 2D axisymmetric problem [1,11,20,27,28] to 3D solution that included all the complexities of the friction clutch system. The new solution can be applied for more complex problems that occur in the clutch system such as hotspot, wear, flash temperature, the effect of roughness, etc.

4. Analytical Model

The thermal analysis conduction by means of the finite element method was compared with results that were obtained from the one-dimensional analytical model. In order to apply this model, the following assumptions were introduced:
1. The frictional heat is totally absorbed in the normal direction to the contact surface by the friction pair elements and they are only significantly heated in the near distance from the friction zone.
2. Wear and convective heat exchange on the lateral surfaces of the clutch elements are negligible.
3. Thermal sensitivity of friction materials and the instability of friction coefficient under temperature variations are negligible due to relatively low maximum temperature level.
4. The frictional heat is uniformly distributed over the friction surface area.
5. The perfect thermal contact conditions were adopted on the contact surfaces.
These assumptions are reasonable only for the short-duration, heavy-loaded friction processes [12,13,14,15,16]. Owing to the symmetry of the problem, it was assumed that the friction processes on the flywheel–clutch disc and clutch disc–pressure plate contact surfaces are the same. Based on the above assumptions, the one-dimensional, boundary–value problem of heat conduction formulated for a semispace–semispace tribosystem can be written as follows [13,14]:
2 T ( z , t ) z 2 = 1 α c T ( z , t ) t ,   z > 0 ,   0 < t t s ,
2 T ( z , t ) z 2 = 1 α f T ( z , t ) t ,   z < 0 ,   0 < t t s ,
K f T ( z , t ) z | z = 0 K c T ( z , t ) z | z = 0 + = q ( t ) ,   0 < t t s ,
T ( 0 , t ) = T ( 0 + , t ) T ( t ) ,   0 < t t s ,
T ( z , t ) 0 ,   | z | ,   0 < t t s ,
T ( z , 0 ) = T a ,   | z | < .
All values referring to the clutch disc and flywheel (or pressure plate) will, respectively, have subscripts c and f. Taking into account the Equations (1)–(3). The specific friction power q ( t ) , that presented in the boundary condition Equation (21), has the following form:
q ( t ) = q 0 q * ( t ) = μ   p   r i   ω o   q * ( t ) = μ   p   r i   ω o ( 1 t / t s ) ,     0 t t s ,
where q 0 —the specific friction power (nominal value), r i —inner radius of the friction surface, and the rest parameters are included in the Table 1. The above problem of heat conduction can be written in the dimensionless form:
2 T * ( ζ , τ ) ζ 2 = 1 α * T ( ζ , τ ) τ ,   ζ > 0 ,   0 < τ τ s ,
2 T * ( ζ , τ ) ζ 2 = T ( ζ , τ ) τ ,   ζ < 0 ,   0 < τ τ s ,
T * ( ζ , τ ) ζ | ζ = 0 K * T * ( ζ , τ ) ζ | ζ = 0 + = q * ( τ ) ,   0 < τ τ s ,
T * ( 0 , τ ) = T * ( 0 + , τ ) T * ( τ ) ,   0 < τ τ s ,
T * ( ζ , τ ) 0 ,   | ζ | ,   0 < τ τ s ,
T * ( ζ , 0 ) = 0 ,   | ζ | < .
where,
T * = T T a T 0 ,   T 0 = q 0 a f K f ,   ζ = z a f ,   τ = α f t a f 2 ,   τ s = α f t s a f 2 ,   α * = α f α c ,   K * = K f K c ,
a f = 3 α f t s —effective depth of heat penetration in the flywheel (or pressure plate). Based on Duhamel’s theorem, the exact solution to the initial boundary problem of heat conduction Equations (26)–(32) can be found from the following relation:
T * ( ζ , τ ) = 0 τ q ( s ) τ T ( 0 ) ( ζ , τ s ) d s ,   | ζ | < ,   0 τ τ s ,
where T ( 0 ) is the known formula obtained from solution to the problem Equations (26)–(32) for constant friction power q ( 0 ) ( τ ) = 1 , 0 τ τ s :
T ( 0 ) ( ζ , τ ) = 2 ( 1 f c ) τ   ierfc   ( n ζ / 2 τ ) ,   0 τ τ s ,
ierfc ( x ) = π 1 / 2 e x 2 x [ 1 erf ( x ) ] ,   n = {                   1                       for       ζ 0 , α * 1               for   ζ 0 .
where erf ( x ) is Gauss error function and f c is already defined in Equation (4). Substituting the derivative of the solution Equation (24) and the specific friction power Equation (25) into the Equation (33), it will achieve,
T * ( ζ , τ ) = ( 1 f c ) π 1 / 2 0 τ ( 1 s τ s ) ( τ s ) 1 / 2 e ( n ζ / 2 τ ) 2 d s ,   | ζ | < ,   0 τ τ s ,
On the friction surface for ζ = 0 the Equation (36) takes the form:
T * ( τ ) = ( 1 f c ) π 1 / 2 0 τ ( 1 s / τ s ) ( τ s ) 1 / 2 d s ,   0 τ τ s .
After integration of Equation (37), it was achieved the solution to the initial boundary problem Equations (26)–(32), in the following form:
T * ( τ ) = 2 ( 1 f c ) [ 1 2 τ / ( 3 τ s ) ] τ / π ,   0 τ τ s .
which is convergent to the results achieved by other researchers [12,13,14].

5. Results and Discussions

This research paper consists of two parts, in the first part, a new 3D FE model of the dry friction clutch was built from scratch in order to enhance the numerical solution of the thermal problem during the slipping period. Meanwhile, in the second part, a new analytical solution was obtained based on the methodical model that was derived from scratch to obtain the variation of temperature and used it to validate the new numerical model.
The numerical and analytical results of variations for maximum temperature that occurred in the contacting surfaces during the complete slipping period are illustrated in Figure 7. Generally, it can be seen that the numerical results are very accurate and the maximum difference with analytical results is not exceeded 1%. According to the obtained results, it can be seen three important points in the curve that represents the thermal behaviour of the clutch system. Firstly, the initial point (Ti = 300 K) that occurred at ts = 0; secondly the peak point where the maximum temperature occurred (Tmax = 443 K at ts = 0.2 s). While, the last point is the end point that represents the end of the slipping period (Tf = 384 K at ts = 0.4 s). It can be observed that the peak point occurred at approximately the middle of the slipping period. Where, the maximum increasing in the temperature that occurred during the slipping period was 120.8 K. The temperatures that were generated during the sliding period were a function of many parameters, such as torque variation during the engagement period, surface roughness, level of cooling, thermal properties of the contacting materials, friction characteristics, sliding speed, etc.
Figure 8 shows the distribution of temperatures of the main elements of the clutch system (friction clutch disc, flywheel and pressure plate) at different times of the slipping period. It is observable that the distribution of surface temperature was approximately uniform, that corresponding to the applied frictional thermal load (uniform distribution). The results presented the thermal behaviour of the used frictional clutch disc, where the uniform wear rate on the contact surfaces was assumed. Based on this theory, the frictional heat generation is applied uniformly on the sliding surfaces. Additionally, the thermal behaviour of the new clutch system is essential, while the pressure will be uniform over the contacting area. This will lead to focusing the highest fractional generated at the zone of the outer disc of the clutch. This important point will be investigated deeply in the subsequent investigation.

6. Conclusions and Remarks

This research is an essential achievement in the direction of obtaining the optimal numerical and analytical simulations of the frictional clutch system. A new three-dimensional FE model was built using a developed numerical approach after verifying the reliability of the temperature results of the clutch system during start-up based on the comparison made with new enhanced analytical solutions. The main goal was to reach the highest accuracy of the numerical results, so this presented approach is the basis for designing the friction clutches.
The most important point that has been proven is the reliability of the developed numerical model after the obtained results were compared with the analytical solution of temperatures during the entire slip period. Extreme temperatures are the most important results that must be focused on, as the level of this temperature must be within the permissible temperature of the selected materials. Whereas the maximum surface temperatures occurred approximately at the mid-slide time for all cases. In addition, when assuming that the rate of wear is uniform over the clutch surfaces, the temperature distribution is uniform over the surfaces of the clutch elements, as explained in the results. The amounts of frictional heat that entered to the steel elements (flywheel and pressure plate, 92.93% of total amount of heat) were much greater than the frictional heat that entered to the friction disc (only 7.07% of the total amount of heat) due to the poor thermal properties of the friction materials. Furthermore, it was found that the highest temperatures appeared on the contact surfaces, while if it moved towards the thickness, it was found that the effect of heat generated approached zero near the middle of the thickness of the frictional material.
There are some researchers that developed finite element models of a clutch disc to study the thermal buckling, in order to find the critical value of the buckling temperature [29]. On the other hand, some research papers focused on the mathematical modelling in order to determine the behaviour of the automotive brakes and clutches during the thermal processes, taking into consideration the unequal distribution of the heat flux and thermal conductivity. They also compared the results with experimental work in order to verify the developed mathematical models [30,31].
A future stage in this research project is to find the numerical solutions to more complex problems of the friction clutch system such as the effects of surface roughness, thermal stresses, unsymmetric thermal load, hotspots, etc.

Author Contributions

L.A.S.: software and review; K.T.: analytical methodology and writing—review and editing; M.J.J.: writing—review and editing; O.I.A.: methodology and writing—review and editing; A.M.A.: investigation and visualization. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The study did not report any data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abdullah, O.I.; Schlattmann, J. Effect of band contact on the temperature distribution for dry friction clutch. Tribol. Ind. 2013, 35, 317–329. [Google Scholar]
  2. Al-Shabibi, A.M. Transient behavior of initial perturbation in multidisk clutch system. Tribol. Trans. 2014, 57, 1164–1171. [Google Scholar] [CrossRef]
  3. Al-Shabibi, A.M. Solution of heat conduction problem in automotive clutch and brake systems. Heat Transf. Summer Conf. 2008, 48470, 321–326. [Google Scholar]
  4. Senatore, A.; D’Agostino, V.; Di Giuda, R.; Petrone, V. Experimental investigation and neural network prediction of brakes and clutch material frictional behaviour considering the sliding acceleration influence. Tribol. Int. 2011, 44, 1199–1207. [Google Scholar] [CrossRef]
  5. Pisaturo, M.; Senatore, A. Simulation of engagement control in automotive dry-clutch and temperature field analysis through finite element model. Appl. Therm. Eng. 2016, 93, 958–966. [Google Scholar] [CrossRef]
  6. Pica, G.; Cervone, C.; Senatore, A.; Lupo, M.; Vasca, F. Dry dual clutch torque model with temperature and slip velocity effects. Intell. Ind. Syst. 2016, 2, 133–147. [Google Scholar] [CrossRef] [Green Version]
  7. Abdullah, O.I.; Schlattmann, J.; Majeed, M.H.; Sabri, L.A. The distribution of frictional heat generated between the contacting surfaces of the friction clutch system. Int. J. Interact. Des. Manuf. 2019, 13, 487–498. [Google Scholar] [CrossRef]
  8. Abdullah, O.I.; Schlattmann, J.; Jobair, H.; Beliardouh, N.E.; Kaleli, H. Thermal stress analysis of dry friction clutches. Ind. Lubr. Tribol. 2018, 72, 189–194. [Google Scholar] [CrossRef]
  9. Abdullah, O.I.; Schlattmann, J.; Lytkin, M. Effect of surface roughness on the thermoelastic behaviour of friction clutches. FME Trans. 2015, 43, 241–248. [Google Scholar]
  10. Kuciej, M.; Topczewska, K. Determination of the maximum temperature of the friction pad in a one-disk clutch. Mater. Sci. 2019, 54, 819–826. [Google Scholar] [CrossRef]
  11. Abdullah, O.I.; Schlattmann, J.; Majeed, M.H.; Sabri, L.A. The Temperatures Distributions of a Single-Disc Clutches Using Heat Partitioning and Total Heat Generated Approaches. Case Stud. Therm. Eng. 2018, 11, 43–54. [Google Scholar] [CrossRef]
  12. Yevtushenko, A.; Kuciej, M.; Topczewska, K. Effect of the temporal profile of the friction power on temperature of a pad-disc brake system. J. Theor. Appl. Mech. 2019, 57, 461–473. [Google Scholar] [CrossRef]
  13. Topczewska, K.; Schlattmann, J.; Abdullah, O.I. Temperature and thermal stresses distributions in a dry friction clutch. J. Theor. Appl. Mech. 2020, 58, 351–360. [Google Scholar] [CrossRef]
  14. Olesiak, Y.; Pyryev, A. Yevtushenko: Determination of temperature and wear during braking. Wear 1997, 210, 120–126. [Google Scholar] [CrossRef]
  15. Yevtushenko, A.; Kuciej, M.; Topczewska, K. Some theoretical model for determining the temperature field of a multi-disk brake. Adv. Mech. Eng. 2020, 12, 1–15. [Google Scholar] [CrossRef]
  16. Yevtushenko, A.; Kuciej, M.; Topczewska, K. Frictional Heating during Braking of the C/C Composite Disc. Materials 2020, 13, 2691. [Google Scholar] [CrossRef]
  17. Yevtushenko, A.; Kuciej, M.; Topczewska, K. Analytical model for investigation of the effect of friction power on temperature in the disk brake. Adv. Mech. Eng. 2017, 9. [Google Scholar] [CrossRef]
  18. Stevens, K.; Tirovic, M. Heat dissipation from a stationary brake disc, Part 1: Analytical modelling and experimental investigations. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2018, 232, 1707–1733. [Google Scholar] [CrossRef] [Green Version]
  19. Abdullah, O.I.; Abd Al-Sahb, W.; Al-Shabibi, A. Finite element analysis of transient thermoelastic behavior in multi-disc clutches. In SAE 2015 World Congress and Exhibition; SAE International: Pittsburgh, PA, USA, 2015. [Google Scholar]
  20. Yevtushenko, A.; Grzes, P. Mutual influence of the sliding velocity and temperature in frictional heating of the thermally nonlinear disc brake. Int. J. Therm. Sci. 2016, 102, 254–262. [Google Scholar] [CrossRef]
  21. Al-Shabibi, A.M. Solution of a Non-Linear Heat Conduction Problem in Automotive Clutch and Brake Systems. In Proceedings of the ASME 2009 Heat Transfer Summer Conference collocated with the InterPACK09 and 3rd Energy Sustainability Conferences, San Francisco, CA, USA, 19–23 July 2009; American Society of Mechanical Engineers: New York, NY, USA, 2009; pp. 415–419. [Google Scholar]
  22. Xie, C.; Wu, G.; He, S.; Zhang, Y.; Tian, J. Finite Element Analysis of Thermal Load Characteristics of Dry Dual Clutch. J. Phys. Conf. Ser. 2021, 1748, 052047. [Google Scholar] [CrossRef]
  23. Bao, H.; Kong, W.; Hou, X.; Zhu, R. Analysis on temperature field of friction pair of aviation friction clutch based on different groove shapes of friction disk. J. Mech. Sci. Technol. 2021, 35, 3735–3742. [Google Scholar] [CrossRef]
  24. Zhang, Z.; Zou, L.; Liu, H.; Chen, Y.; Zhang, B. Effects of Operating and Material Parameters on the Thermal Characteristics of a Wet Clutch. Available online: https://0-journals-sagepub-com.brum.beds.ac.uk/doi/10.1177/16878140211034101 (accessed on 9 October 2021).
  25. Reddy, J.N.; Gartling, D.K. The Finite Element Method in Heat Transfer and Fluid Dynamics; CRC Press: Boca Raton, FL, USA, 2010. [Google Scholar]
  26. Shigley, E.J.; Mischke, R.C. Standard Handbook of Machine Design, 2nd ed.; McGraw-Hill: New York, NY, USA, 1996. [Google Scholar]
  27. Al-Shabibi, A.M.; Barber, J.R. Transient Solution of 2-D TEI Problem Using Reduced Order Models. SAE Trans. 2001, 110, 1244–1248. [Google Scholar] [CrossRef]
  28. Al-Shabibi, A.M. Solution for transient heat conduction problem with low Biot number. In Proceedings of the 9th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics, Paceville, Malta, 16–18 July 2012. [Google Scholar]
  29. Zhao, J.; Chen, Z.; Yang, H.; Yi, Y.-B. Finite element analysis of thermal buckling in automotive clutch plates. J. Therm. Stress. 2016, 39, 77–89. [Google Scholar] [CrossRef]
  30. Awrejcewicz, J.; Grzelczyk, D. Modeling and Analysis of Thermal Processes in Mechanical Friction Clutch—Numerical and Experimental Investigations. Int. J. Struct. Stab. Dyn. 2013, 13, 1340004. [Google Scholar] [CrossRef]
  31. Awrejcewicz, J.; Grzelczyk, D. Modeling and Analytical/Numerical Analysis of Wear Processes in a Mechanical Friction Clutch. Int. J. Bifurc. Chaos 2011, 21, 2861–2869. [Google Scholar] [CrossRef]
Figure 1. Main elements of clutch system (single disc).
Figure 1. Main elements of clutch system (single disc).
Computation 09 00122 g001
Figure 2. Three-dimensional model of frictional clutch system.
Figure 2. Three-dimensional model of frictional clutch system.
Computation 09 00122 g002
Figure 3. Boundary and load conditions of the friction clutch system.
Figure 3. Boundary and load conditions of the friction clutch system.
Computation 09 00122 g003
Figure 4. The mesh independence test of the FE model of clutch system.
Figure 4. The mesh independence test of the FE model of clutch system.
Computation 09 00122 g004
Figure 5. Flowchart of developed numerical approach (FE) to find solution of transient thermal problem.
Figure 5. Flowchart of developed numerical approach (FE) to find solution of transient thermal problem.
Computation 09 00122 g005
Figure 6. 3D FE model of frictional clutch. (a) Complete FE model; (b) FE model (sector).
Figure 6. 3D FE model of frictional clutch. (a) Complete FE model; (b) FE model (sector).
Computation 09 00122 g006
Figure 7. Highest surface temperature of clutch disc through the slipping time. Dotted–dashed line presents FE results and solid line presents analytical results (32), (38).
Figure 7. Highest surface temperature of clutch disc through the slipping time. Dotted–dashed line presents FE results and solid line presents analytical results (32), (38).
Computation 09 00122 g007
Figure 8. Temperature [K] distributions of disc clutch, flywheel and pressure plate at selected time.
Figure 8. Temperature [K] distributions of disc clutch, flywheel and pressure plate at selected time.
Computation 09 00122 g008
Table 1. The material properties and working parameters of clutch system [5].
Table 1. The material properties and working parameters of clutch system [5].
ParaameterValue
Inner radius of clutch disc and axial cushion (ri-c & ri-cu), [m]0.06
Outer radius clutch disc and axial cushion (ro-c & ro-cu), [m]0.0792
Frictional facing thickness of tcu, [m]0.004
Inner radius of pressure plate (ri-p), [m]0.06
Outer radius of pressure plate (ro-p), [m]0.091
Pressure plate thickness, [m]0.01
Inner radius of flywheel (ri-f), [m]0.0485
Outer radius of flywheel (ro-f), [m]0.097
Flywheel thickness of tf, [m]0.0194
Applied pressure pa, [MPa]1
Coefficent of friction (μ)0.3
Angular sliding velocity (ωo), [rad/s]295
Number of friction surfaces, n2
Density of friction material (ρc), [kg/m3]1000
Density of axial cushion, flywheel and pressure plate (ρcu, ρf & ρp), [kg/m3]7200
Specific heat of friction material (cc), [J/kg K]1400
Specific heat of axial cushion, flywheel and pressure plate (ccu, cf & cp), [J/kg K]450
Thermal condutivity of of friction material (Kc), [W/mK]0.75
Thermal condutivity of of axial cushion, flywheel and pressure plate (Kcu, Kf & Kp), [W/mK]56
Slipping time (ts), [s]0.4
Heat transfer coefficient (h), (W/m2K)40.89
Initial temperature (Ta), [K]300
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sabri, L.A.; Topczewska, K.; Jweeg, M.J.; Abdullah, O.I.; Abed, A.M. Analytical and Numerical Solutions for the Thermal Problem in a Friction Clutch System. Computation 2021, 9, 122. https://0-doi-org.brum.beds.ac.uk/10.3390/computation9110122

AMA Style

Sabri LA, Topczewska K, Jweeg MJ, Abdullah OI, Abed AM. Analytical and Numerical Solutions for the Thermal Problem in a Friction Clutch System. Computation. 2021; 9(11):122. https://0-doi-org.brum.beds.ac.uk/10.3390/computation9110122

Chicago/Turabian Style

Sabri, Laith A., Katarzyna Topczewska, Muhsin Jaber Jweeg, Oday I. Abdullah, and Azher M. Abed. 2021. "Analytical and Numerical Solutions for the Thermal Problem in a Friction Clutch System" Computation 9, no. 11: 122. https://0-doi-org.brum.beds.ac.uk/10.3390/computation9110122

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