Abstract
The frontend accessory drive belt drive system is a critical component in the vehicle engine. To avoid thermal deterioration under static state operating conditions, the thermal distribution for the belt drive system at each condition must be determined in an efficient manner. Due to the numerical approach is not feasible to address this concern because of its high computational cost, this paper proposes a reliable and efficient novel analytical thermal model to achieve this goal. This work develops the stateoftheart heat transfer ordinary differential equations (ODEs) describing the thermal flow and heat dissipations on the complex structures of pulleys. Then it integrates these ODEs with heat transfer governing equations of the belt and heat exchanges to establish an innovative system of equations that can be solved within a few seconds to provide temperature plots. Moreover, experiments were conducted on a dynamometer to verify the accuracy of the proposed model under a wide range of conditions. The results indicate that the measured temperatures are in good agreement with the corresponding analytical results. Owing to its efficiency, the proposed model can be integrated with other mechanical characterizations of the belt drive system in terms of design, optimization, and thermal fatigue analyses.
Introduction
A frontend accessory drive (FEAD), as shown in Fig. 1, is a critical belt drive system in an automotive engine. It consists of a driving pulley (DR pulley), a belt, and a number of driven pulleys (DN pulleys). An engine crankshaft drives the DR pulley, inputs energy into the belt drive system, and transmits energy through the belt to DN pulleys. Then, these pulleys drive other parts of the engine, such as a water pump or an alternator, to maintain engine operation. Currently, the DN pulley material, which is a kind of fiberreinforced plastic (FRP) with low thermal conductivity^{1}, is prone to deterioration and blistering under high temperature in the underthehood environment and a large amount of heat generation within the FEAD system during engine operation. Eventually, this deterioration causes pulley structural failure. To overcome this problem, this study developed a thermal analytical model that can predict realtime temperatures for target locations on the belt and pulleys to prevent material thermal failure under engine operating conditions. Unlike computationally expensive methods such as computational fluid dynamics^{2}, this model functions as a fast and economical tool to provide blackbox analysis for subsequent research investigations. For example, this analytical model can rapidly provide temperature variations under numerous operating conditions, enabling material scientists to address the thermal fatigue of FRP and thus improve its reliability.
The first step of thermal analysis for a belt drive system is to confirm the heat generation flux and its locations. A number of studies related to the power loss of the belt drive system have focused on improving the efficiency of the belt drive system. However, the power loss can also provide the heat generation within the system in this case. One prevalent theory for Vrib rubber belt drive systems considers five forms of energy loss based on Vbelt movement mechanisms, namely, belt bending, stretching, shear, radial compression, and slip^{4}. Manin et al.^{5} extended this method to a multipulley serpentine belt transmission system. Silva et al.^{6,7} considered that primary power losses are caused by belt hysteresis, which can be classified as the bending, stretching, shear, flank, and radial compression of belt rubber, and applied this theory to a polyV belt drive. Bertini et al.^{8} developed an analytical power loss model by considering belt sliding, hysteresis, and engagement/disengagement frictional losses. Chen et al.^{9} experimentally investigated the power loss of a rubber Vbelt based on various speed ratios, belt tensions, rotational speeds, external loads, and diameters of a pulley. Various factors may influence the belt drive mechanism and its power loss; one such factor is belt slip. Balta et al.^{10} used the response surface method to establish a relationship between belt slip behavior and belt drive parameters under various operating conditions. The transverse dynamic hysteretic damping characteristics of a serpentine belt were estimated by utilizing the variable stiffness and damping model^{11}. Kim et al.^{12} developed a belt force equation based on the classical Euler formula to predict the normal and tangential belt force distribution of a flat belt drive. Qin et al.^{13} developed an analytical model for the friction damping of round clamp band joints, which can predict the energy dissipation caused by the friction contact between joint components. Shen et al.^{14} used a dynamic finite element model (FEM) to simulate belt contact deformations and belt–pulley interaction, which can potentially be used to predict transient power loss. Other studies have focused on the power loss of another type of belt drive system, i.e., continuously variable transmission (CVT). Julió et al.^{15} presented a model to predict the transmission ratio time response for the rubber belt CVT under various operating conditions. This model can be used for CVT design and optimization. Zhu et al.^{16} studied the torque and speed loss of the rubber Vbelt CVT based on its mechanism. Several methods have been developed to improve the efficiency of this type of CVT.
The second step of thermal analysis is to use the information of power loss or heat generation to perform heat transfer analysis. Wurm et al.^{17,18} simulated the heat transfer analysis for a CVT system using CFD under critical load cases. They used a novel method to simulate rotational symmetric temperature profiles for nonrotating pulleys in the numerical program. Gerbert^{19} proposed a static analytical thermal model for a simple flat twopulleyonebelt drive system. This model established a set of ODE equations for the heat conduction within the components and the heat exchange of belt–pulley contacts, combined with power loss as the heat source, to calculate the temperature distributions of the belt and pulleys. Other studies have focused on heat analysis and thermal influence for only belts. Abe et al.^{20} developed a thermal analysis model for a typical timing belt. The analytical model considered the bending deformation of a timing belt cord as the primary source of heat generation and adopted a tirebased thermal analysis method to calculate heat generation and belt temperature. Merghache et al.^{21} proposed a thermal model for heat transfer of the AT10 timing belt. Belt–pulley heat flux was calculated using a mathematical model, and it was verified by the results of numerical simulations. Babak et al.^{22} investigated the transient heat transfer and stress wave propagation in nanocomposite sandwich plates with carbon nanotubes. Wu et al. performed experiments to investigate local and average heat transfer characteristics on the surface of the test wheel by using naphthalene sublimation technique^{23}. Krane et al. developed an analytical model to calculate the temperature distribution and heat transfer capacity of a lightweight belttype radiator for power plants^{24}. This model analyzed two thermal conductions, i.e., from the drum to the belt and from the thermal dissipation of the belt to the environment, to predict temperatures under steadystate operation. Some experiments focused on the turbulence properties caused by the rotation of the solid hollow circular disc with rectangular wings insert in a circular tube^{25}. Zhang et al.^{26} investigated the effects of rotating angular speed and pin configuration on the temperature maps and convective heat transfer characteristics on a rotating disk by use of numerical method. Kayhani et al. provided the analytical solutions for the heat transfer in 2dimensional cylindrical composite laminate^{27}. Merghache et al.^{28} experimentally investigated various tooth forms of synchronous belts and their impact on belt temperature. McPhee and Johnson^{29} focused on the heat transfer and dissipation for a brake disc. Belhocine and Bouchetara presented the thermal–mechanical coupled analysis of the dry contact between the brake disk and pads for the vehicle braking system^{30}. While Talati and Jalalifar focused on the surface temperatures and heat partition on the pad–disk contact surface^{31}. A numerical–experimental approach to calculate surface temperatures and flow partition coefficient in padondisc contact application under tribological solicitations^{32}, which is similar to the belt–pulley contact in this study. Detailed transient temperature calculations are presented for the thermal design of a high pressure compressor rotor of an aeroengine^{33}. Song et al.^{34} developed an FEM to investigate the influence of different ambient temperatures on belt shear deformation, tension and velocity variation, and contact stress distribution during operation. Chen et al.^{35} examined the relationship between rubber belt friction and the presence or absence of an interfacial ice film under different temperatures and belt operating conditions. The heat conduction model for a twolayer cylinder and its analytical algorithm^{36} also provide useful knowledge for the thermal calculation of the pulley in this case.
However, to the best of our knowledge, the aforementioned studies have provided accurate but not efficient temperature prediction for a belt drive system. Moreover, only few works have focused on the thermal analysis of a multipulley belt drive system. Toward this end, this paper presents an innovative thermal model that uses an analytical algorithm to calculate thermal distribution and simplify it as a system of equations that can be solved using low computing resources. The model can provide temperature prediction within a few seconds for the optimization of part structure in the design stage; this is considerably superior to any numerical analysis method. Additionally, the above analytical algorithm considers complex pulley structures to ensure the accuracy of results. This model also can be used to inversely determine the optimal operating conditions to avoid failure and enhance the system lifecycle. Similar research also performed in various applications. Liu and Ma developed a novel 3D levelset topology optimization approach to provide the best shape and topology for parts^{37}. Ning et al. adopted chip formation model and iterative gradient search method to calculate the Johnson–Cook model constants for ultrafinegrained titanium^{38}.
The remainder of this paper is organized as follows: “Research method” section describes the procedures for establishing the analytical algorithms for the thermal model to calculate the heat dissipation of pulley fins, belt–pulley interaction, and the final temperature distributions for multipulley systems. “Experiment” section discusses the setup and procedure of the belt drive system experiment. “Results and discussion” section presents and analyzes the experimental and analytical results to validate the proposed model. The analytical results are in good agreement with the experimental temperature data. Finally, “Concluding remarks” section concludes the paper.
Research method
Preliminaries
One preliminary for establishing a thermal model is to identify each pulley and contact surface. In an arbitrary multipulley belt drive system, the total number of pulleys is defined as N. Each pulley is assigned a value (1,2,…, N) in a counterclockwise direction starting from the DR pulley in the belt layout. Further, the total number of belt–pulley contact surfaces is also N, and the nth belt–pulley contact surface denotes the contact surface between the belt and the nth pulley in the system. An example of a twopulleyonebelt system is shown in Fig. 2. The left DR pulley is No. 1 and the right DN pulley is No. 2, under the condition of the left pulley driving the right pulley. The belt–pulley contact surfaces are highlighted as orange surfaces in Fig. 2.
Another preliminary for the thermal model is to calculate the heat flux generated within the belt drive system. There are five forms of energy losses^{4} according to movement mechanisms, as shown in Fig. 3. These power losses are converted into heat according to the conservation of energy. To simplify the thermal analysis of the belt drive system, these losses are subsequently classified into two types of heat sources according to their generation locations, namely, belt internal power loss P_{h}^{39} inside the belt and contact surface power loss P_{fn} at the nth belt–pulley contact surface. Figure 2 shows that P_{f1} and P_{f2} are located at the belt–pulley contact areas and P_{h} is located inside the belt. In this study, these power losses are calculated based on the power loss algorithm developed by Gerbert, which has been used for decades.
Overall thermal calculation
Two assumptions are made in this thermal analysis. The first is that the temperature of the belt is uniform because of the high velocity and small thickness (5–10 mm) of the belt. The second is that the thermal distribution on the pulley is tangentially homogeneous and radially symmetric owing to the highspeed spinning of the pulley and symmetric flow conditions^{40,41,42,43}. The temperature gradient is along the pulley radial direction from the outer radius to the inner radius of the pulley.
Figure 4 shows the general calculation procedure for the thermal model of an arbitrary multipulley belt drive system in this study. The first stage is to determine the heat flux and its locations within the system. Power losses P_{h} and P_{fn} are transformed into heat within the belt system. When the system runs under a constant operating condition and reaches thermal equilibrium, all generated heat is dissipated to the environment through the belt surfaces and each pulley surface.
where the unknown value Φ_{b} is the heat dissipation flux per second from the belt surfaces to the ambient environment and the unknown value Φ_{pn} is the heat dissipation flux per second from the nth pulley surface to the ambient environment. Moreover, there is no heat source inside any pulley component (Fig. 2). A portion of P_{fn} flows into the pulley and dissipates to the environment from the pulley surfaces. Therefore, a new coefficient, ξ_{n}, is introduced, which represents the percentage of P_{fn} that flows into the nth pulley and dissipates to the environment. Thus,
Combining Eqs. (1) and (2) gives
Equation (3) implies that all generated heat that does not flow into the pulley is dissipated to the environment through the belt surfaces.
In the second stage, the model investigates the thermal behavior via three types of thermal analysis and establishes a system of equations. A function is defined for each type of analysis to illustrate the overall calculation procedure; these functions are explained later. First, function Φ_{b} = f(T_{b}) is established for the belt internal thermal analysis to describe the relationship between belt surface temperature T_{b} and the heat dissipation flux per second, Φ_{b}. This calculation is highly related to the following parameters: belt geometries matrix D_{b} (surface areas, lengths of edges and etc.), heat transfer coefficient between the belt and pulley h_{exn} and ambient temperature T_{a}. The following equation is obtained by substituting this function into Eq. (3):
Second, for the belt–pulley heat exchange analysis, the temperature at the belt surface, T_{b}, may be different from that at the outer radius of the nth pulley, T_{pn}, causing pulley–belt heat flux exchange and disturbing the distribution of P_{fn} on the pulley, Φ_{pn}. This calculation is highly related to the following parameters: contact surface area between the belt and the nth pulley A_{exn}, thermal conductivity for the belt Λ_{b}, heat transfer coefficient with surrounding air h_{b}. The following function can be used to describe the influence of the thermal gradient (T_{b} − T_{pn}) on the heat allocation of P_{fn}:
Substituting Φ_{pn} with ξ_{n}P_{fn} in Eq. (2) gives
Third, there is a relationship between Φ_{pn} and T_{pn} in the pulley internal thermal analysis. This can be explained as follows: the higher the pulley surface temperature, the faster is the heat dissipation to the environment from the pulley surfaces. This calculation is highly related to the following parameters: pulley geometries D_{pn} (surface areas, lengths of edges and etc.), thermal conductivity for nth pulley Λ_{pn}, heat transfer coefficient between the nth pulley and air h_{pn}, and T_{a}. This relationship is simply defined through a function as
Substituting Φ_{pn} with ξ_{n}P_{fn} in Eq. (2) gives
In the final stage, the model calculates the temperatures within the belt drive system. Equations (4), (6), and (8) form a system of equations for obtaining belt temperature T_{b} and the outer surface temperature of each pulley, T_{pn}. Consider the twopulleyonebelt system as an example. There are two pulleys and one belt in the system, and the equation set is
There are five unknown values in these five equations, namely, T_{p1}, T_{p2}, T_{b}, ξ_{1}, and ξ_{2}; hence, there is only one solution. When a third pulley is added to this system, there will be additional versions of Eqs. (6) and (8) in Eq. (9), which ensures that the two additional unknown values, T_{p3} and ξ_{3}, have only one solution. Extra pulleys can be added to the system in the same manner. By utilizing this technique, this model can calculate the temperature for an arbitrary belt drive system.
However, functions f(T_{b}), g(T_{b} − T_{pn}, T_{pn}), and k(T_{pn}) must be established to make this algorithm functional. The details of establishing these functions are presented below.
Belt inner calculation
The belt dissipates heat to the ambient environment during its operation. There is one assumption that the temperature inside the belt is assumed to be uniform because it has a simple structure with a high lengthtothickness ratio. Consequently, this study uses T_{b} to represent the surface and inner temperature of the belt. Consequently, the relationship between the heat dissipation rate from the belt surfaces, Φ_{b}, and belt temperature T_{b} is given by
where A_{b} is the total belt surface area exposed to the environment, h_{b} is the heat transfer coefficient of the belt rubber, and T_{a} is ambient temperature. Combining Eq. (10) with Eq. (3) gives
Individual pulley–belt heat exchange calculation
The heat analysis at the belt–pulley engaged surfaces is complex. The nth belt–pulley engaged surface involves not only heat flux exchanges between the belt and pulley but also the distribution of P_{fn} to the belt and pulley. The heat exchange flux per second at the nth belt–pulley engaged surface (Φ_{exn}) from the belt to the nth pulley is related to the pulley–belt temperature difference (T_{b} − T_{pn}) as follows:
where A_{exn} is the area of the nth belt–pulley engaged surface and h_{exn} is the thermal contact conductance coefficient.
Moreover, Φ_{exn} influences the heat allocation of P_{fn} at the belt–pulley contact surfaces. When T_{pn} = T_{b}, the pulley acquires half of the frictional heat and ξ_{n} = 0.5. There is no heat exchange between the belt and pulley in this situation. When T_{pn} < T_{b}, heat flows from the belt to the pulley and the pulley acquires additional Φ_{exn}, making ξ_{n} > 0.5. Hence, this phenomenon can be expressed as
Combining Eqs. (12) and (13) yields the following equation, which provides the details of function g(T_{b} − T_{pn}, P_{fn}) in Eq. (5).
Individual Pulley inner thermal calculation
A general industrial pulley can be simplified by suppressing the irrelevant chamfer and fillet features. Its cross section is shown in Fig. 5a. It consists of top and bottom flanges and a middle web section attached by fins on two sides to improve structural stiffness under high belt hub load.
One assumption is made in the thermal calculation for the pulley structure. Owing to the pulley rotates at highspeed spinning during the operation, the temperature profiles are tangentially homogeneous and radially symmetrical distribution^{40,42,43}. Therefore, the assumption is that the temperatures on two flanges and the middle web are uniform in the axial direction of a pulley. On the other hand, this assumption does not apply to the temperature profile on the fins because of their dramatic surface heat dissipations and large temperature decrease in the pulley’s axial direction. This assumption is only applicable for the wellstructured crosssections such as I, L or C shapes. However, these crosssections account for 95% of the current industrial pulleys used in engine belt drive, making this model widely used in the industrial world.
The pulley can be divided into three layers, as shown in Fig. 5b. Layer 1 is the top flange, layer 2 is the web with two fin sections, and layer 3 is the bottom flange. In this structure, a belt and a shaft with a uniform temperature on their surfaces are in contact with the top and bottom flanges of the pulley, making the temperatures on both flanges uniform along the pulley axial direction as well. Moreover, the web in layer 2 is extremely thin; hence, the temperature gradient in the pulley axial direction is neglected. As a result, the temperature in the hatched areas is considered uniform along the pulley axial direction.
All three layers have rectangular cross sections highlighted as three different hatched areas in Fig. 5b. Consider one rectangular cross section in the ith layer as an example. The outer and inner radii are defined as R_{x} and R_{y}, respectively. Figure 6 shows the thermal behavior in this Sect. ^{44}.
The balance of heat gives
Function T(r) describes the temperature distribution along the pulley radial direction, W_{i} is the half width of the ith layer, Λ_{p} is the thermal conductivity of the pulley material, and h_{p} is the heat transfer coefficient between the pulley surface and ambient environment. Equation (15) changes to an ODE as follows:
However, in this study, as layer 2 consists of fins on both sides, Eq. (16) has a limitation and the thermal model must consider the heat dissipation on the fin surfaces. For a differential increase, dr, shown in Fig. 7, a certain amount of heat flows into the fin through the web–fin connecting area (green area in Fig. 7) and then dissipates to the environment through the fin surfaces. This area has an amplified heat dissipation rate owing to the fin, while the remaining areas (purple area in Fig. 3) have the natural heat dissipation rate, h_{p}, of the pulley material. Therefore, this study introduces a new coefficient, η_{f}, to indicate that the heat dissipation rate in the dashed area increases to η_{f}·h_{p}. When the number of fins connected to the web on one side is N_{f} in the cylinder body of differential dr, the dissipation rate on one side of the body is [N_{f}⋅W_{f}⋅dr⋅η_{f}⋅h_{p} + (2πr⋅dr − N_{f}·W_{f}·dr)·h_{p}]. However, the original dissipation rate is 2πr⋅dr⋅h_{p} if there is no fin attached. Therefore, h_{p} on the web side surfaces is amplified by a factor of [1 + N_{f}⋅W_{f}⋅(η_{f} − 1)/2πr]. Equation (16) is modified to
The equation is applicable to the layers with and without fins. When there is no attached fin, such as in layers 1 or 3, η_{f} is equal to 1 and Eq. (15) reverts to Eq. (14). When fins are attached, such as in layer 2, η_{f} > 1 and Eq. (15) considers the extra heat dissipation on the fin surfaces. To solve Eq. (15), the key factor is to determine the value of η_{f} and whether η_{f} varies along pulley radius r.
The fin connected to the cylinder body of differential dr can be considered as shown in Fig. 8. The temperature on the plate decreases from the web end to the edge end because of heat dissipation. The general heat dissipation equation for this flat plate is^{45}
where A_{fc} is the connected area (green area in Fig. 7) between the fin and web, G(0) is the temperature at A_{fc}, which is equal to T(r), and G(x) is the temperature distribution function along the pulley axial direction. The analytical solution of Eq. (18)^{45} provides the temperature distribution from 0 to fin length L_{f}:
The total heat dissipation for each fin is the summation of the heat dissipation flux per second on the fin end surfaces, Φ_{fs} (dark blue surface in Fig. 8), and the heat dissipation flux per second on the two side surfaces of the fin, Φ_{fe} (yellow surface in Fig. 8):
Therefore, the fin coefficient can be calculated according to its definition:
The abovementioned calculation implies that η_{f} does not change under various web surface temperatures, G(0), when the system is in thermal equilibrium; hence, it is not influenced by radius r or T(r). The fin dimensions (W_{f}, L_{f}) and heat dissipation flux per second (h_{p}) are the only parameters required for calculating η_{f}. Thus, for each pulley, η_{f} is a constant in Eq. (17).
As η_{f} is independent of r, Eq. (17) can be regarded as Kummer’s equation^{46}. Its analytical solution, T(r), provides the temperature distribution of the ith layer as follows:
where constant coefficients C_{1} and C_{2} can be derived from the thermal analysis boundary conditions. The temperature at the inner diameter, T(R_{y}), is not infinity and then C_{2} = 0. Therefore, C_{1} can be calculated as
The temperature distribution for one layer can be obtained as
Assume that T_{pn} or the temperature at the outer radius of layer 1 is known. The temperature distribution of a pulley can be calculated using Eq. (24) three times from layer 1 to layer 3. The connection of two layers shares the same temperature, and hence, T(R_{y}) in one layer is equal to T(R_{x}) in the next layer.
Subsequently, the heat dissipation rate at the two side surfaces in the ith layer is related to the integration of surface temperature distribution.
In function K(r,i), input parameter i determines the value of W_{i}. The total heat dissipated per second at the side surfaces of all three layers (red surfaces in Fig. 5b), Φ_{ps}, is calculated layer by layer. It is given by
where R_{o} and R_{i} are the outer and inner radii of the investigated pulley, respectively, and R_{1} and R_{2} are the inner radii of layers 1 and 2, respectively, as shown in Fig. 5b. Furthermore, this model includes the heat dissipated from other surfaces, Φ_{pb} (blue surfaces marked as 1, 2, and 3 in Fig. 5b), which have a uniform temperature:
where A_{1}, A_{2}, and A_{3} are the surfaces areas identified by balloons 1, 2, and 3 in Fig. 5b, and h_{s} is the heat transfer coefficient between the pulley surface and shaft material. Therefore, the total heat dissipation rate of the investigated pulley is
It is impossible to calculate Φ_{pn} without the value of T_{pn}. However, based on the above equations, Φ_{pn} is directly proportional to T_{pn} when the belt drive system is under thermal equilibrium. A new coefficient, λ_{pn}, is introduced for the nth pulley to connect these two values.
Substituting Φ_{pn} with ξ_{n}P_{fn} in Eq. (2) and function k(T_{pn}) in Eq. (8) gives
Global system thermal calculation
Combining Eq. (9) with the functions in Eqs. (11), (14), and (30) provides the following mathematical model for an arbitrary multipulley belt drive system:
The belt drive system in an engine generally contains N pulleys depending on the design of the engine and its parts. Equation (31) is a system of equations with a maximum of 2(N + 1) unknown variables (T_{p1},…, T_{pN}, ξ_{1},…,ξ_{N} and T_{b}). This study adopts LU decomposition^{47} to solve Eq. (31). In addition, the temperature distribution of a pulley in this system can be calculated layer by layer according to Eq. (24) when the value of T_{pn} is known.
Even though this model involves numerous ODEs and integrations, the analytical solution of each ODE is adopted. Moreover, the model is reduced to a system of equations in the final stage. Only a few seconds are required to obtain the predicted temperature distribution on a typical computer. Therefore, this model can predict the results efficiently.
Experiment
The purpose of this experiment is to validate the accuracy of the proposed thermal model by measuring the thermal distribution of the entire belt drive and comparing it with the results obtained by the thermal model under various operating conditions. A rigorous experimental arrangement is necessary before validating the thermal algorithm with the multipulley configuration at an elevated temperature.
Equipment and layout setup
First, a fivepulley belt drive system is designed (Fig. 9a) and installed on an engine simulator (Fig. 9b)^{48}. Table 1 summarizes the information of this belt drive layout. In this layout, pulley1 is used as the DR pulley, pulley4 is used as the DN pulley, pulley2 contains a hub load sensor under it to measure belt tension, pulley3 is used to measure belt speed, and pulley5 is a tensioner pulley used to reduce the vibration of the belt and keep belt tension constant. In addition, thermal measurement devices are installed close to the measurement locations. Five infrared sensors measure the temperatures at the selected locations, marked S1 to S5 (Fig. 9a). These sensors are the MicroEpsilon CSmiSF15C1^{49}. The measurement range of these sensors is from − 40 to 1030 °C. Its accuracy is within ± 1.5%. A thermal graphic camera is installed to target the DN pulley and measure its temperature distribution during the test. It is FLIR E60. The measurement range of this camera is from − 20 to 130 °C, and its accuracy is within ± 2%. Both types of measure devices meet the requirements of this study because the maximum temperature in this experiment is 121 °C. Further, an ancillary data acquisition system records the values measured by the sensors.
Load case
This experiment also considers the influence of operating conditions on temperature. A wide range of test operating conditions of the belt drive system is vital for validating the results of the thermal model in all aspects. Engine operation varies significantly from idle to the maximum revolution/torque conditions, and it is represented through three parameters: transmitted speed, transmitted torque, and ambient temperature. Table 2 summarizes the selected test points to fully consider the significant variation in engine operating conditions. Only one parameter is changed in each test, and hence, the total number of tests is 12. To investigate the influence of ambient temperature on the temperature distribution of the belt drive system, the ambient temperature in the warm start (WS) engine condition is twice as high as that in the cold start (CS) engine condition. The experiment controls the transmitted speed and torque of the belt drive system by altering the revolution and torque load of the DR pulley. In addition, a customized insulated chamber mounted on the engine dynamometer encloses the belt drive system, and engine ambient temperature is simulated by adjusting internal chamber temperature. Two thermocouples mounted at the two sides of the chamber monitor internal temperature, and a case fan inside the chamber ensures that chamber temperature is uniformly distributed. The temperature distributions of the belt system are measured under the designed working conditions for validating the thermal model using the above devices and equipment as a test bench.
Conduction of experiment
The operation of this experiment is to perform the temperature measurements under the different test operating conditions. There is 12 configurations base on Table 2. There are three tests performed in each configuration to reduce the error. In each test, the system requires about 10–20 min to reach the thermal equilibrium status regarding the various torque and speed loads. The temperature data are recorded three times after the readings from sensors are stable. The measured values are averaged before the comparison and validation for the thermal model.
At the same time, the thermal model uses the abovementioned information related to the belt layout and operating conditions, along with the pulley geometries and thermal properties (h_{p}, Λ_{p}, …) measured using standard procedures, to calculate the temperature distributions of the belt drive system under the designed operating conditions. The comparison between experimental and analytical results is presented in the following section.
Results and discussion
Assumption verification
The two assumptions used in the development of the thermal model are the fundamentals of this model. Hence, it is crucial to verify these two assumptions through experiments to ensure the feasibility of the method employed by this model.
The assumed uniform belt temperature in the thermal model can be verified by the temperatures measured at two different locations on the belt: sensor 2 and sensor 4 (Fig. 9a). Figure 10 shows the temperatures measured at these two locations under four different belt drive conditions. The differences in these measurements are less than 3 °C, indicating that it is feasible to assume that the belt has a uniform temperature during highspeed movement. Hence, these measured temperatures are considered as belt temperatures under different working conditions.
Figure 11 shows the temperature distribution of a running FRP pulley with the lowest rotation speed configuration (600 RPM). A perfectly rotational symmetric temperature distribution is measured. The temperature is the highest at the outer radius, where heat is generated, and it gradually decreases in the pulley radial direction but not in the tangential direction owing to heat dissipation. Thus, the thermal distribution on the pulley is in accordance with the assumption, and the governing equations based on this assumption are feasible for this work.
Experimental and analytical temperature verification
To validate the accuracy of the analytical results provided by the thermal model, this study compares the experimental and analytical temperatures at several selected critical locations under variations in transmitted torque and speed at the ambient conditions of CS and WS. This section describes the influence of the change in each parameter on the experimental temperature distribution and validates whether the model can predict the corresponding analytical temperature.
Temperature comparisons at different working loads
Increase in transmitted torque leads to high slip rate and heat generation, resulting in high temperatures at the belt–pulley engaged surface. The surface on the pulley side is at the outer radii of the pulleys. In this experiment, various torque loads are applied to the DR and DN pulleys. The critical positions are the belt surface and the outer radii of the DR and DN pulleys. The temperatures at these positions are plotted in detail under transmission torque loads of 12–17 N m at the ambient conditions of CS and WS, as shown in Fig. 12. The dashed column indicates the experimental temperatures, and the gray column indicates the analytical temperatures provided by the thermal model. This figure displays accurate temperature predictions by the model. The maximum difference between the experimental and analytical temperatures is 6 °C, and all differences are less than 10% of the measured temperatures.
Temperatures at different speeds
Another important aspect of evaluating the accuracy of the thermal model is to compare the analytical and experimental temperatures under different transmitted speeds. Unlike the effect of torque change, variation in speed alters not only heat generation but also the heat dissipation on exposed surfaces in the belt drive system. Figure 13a–c shows the analytical temperatures at the outer radius of the DR and DN pulleys and the belt surface under the ambient conditions CS and WS with a fixed transmission load of 12 N m and speed ranging from 600 to 1200 RPM. The temperatures on the pulley outer surfaces and belt surface increase with speed. This indicates that the increase in heat generation is more significant than that in the heat transfer coefficients. The small differences between the experimental and analytical results shown in Fig. 13 demonstrate that the thermal model can accurately capture this trend of temperature increase. The maximum temperature difference is only 5 °C, or 10% of the measured results.
Temperatures distributions of the pulley
The temperature distribution of a pulley is important because it is highly related to the heat dissipation of the pulley which influences the accuracy of temperature predictions. The output results from the analytical model are the temperature plot along the pulley’s radial direction. On the other hand, the experimental data is recorded from the five selected points at different radii on the pulley surfaces, as shown in Fig. 14. As for resultant comparisons, the predicted and experimental temperature distribution of the DR and DN pulleys under various speed conditions is displayed in Fig. 15. Based on these two comparisons, the model successfully predicts the decrease in temperature from the outer diameter to the inner diameter in all circumstances, and the predicted temperatures are similar to the measured values. With an accurate temperature plot for the pulley, the thermal model can calculate the heat dissipation and λpn for each pulley. Then, the system of equations can accurately predict the temperatures at pulley surfaces.
Thermal flow in the system
The thermal model provides useful information such as the thermal dissipation rate for each pulley, as shown in Table 3. The thermal dissipation rate is vital because a low thermal dissipation rate causes high heat accumulation and high temperature within a pulley, which finally leads to the thermal failure of the pulley. Hence, the thermal dissipation rate is a good indicator for a pulley in the thermal analysis. As shown in Table 3, pulleys 1 and 2 have high dissipation rates because they are fabricated from conventional steel, while the remaining three pulleys have low rates because of their FRP material, which has low thermal conductivity. Moreover, pulley 1 has a larger dissipation rate compared to pulley 2 because the diameter of pulley 1 is larger (Table 1). Furthermore, pulley 3 has a lower dissipation rate compared to pulley 2, but the values are close. This proves that pulley 3 consists of a fin section to increase heat dissipation, even though it is fabricated from FRP. In general, this table provides the impact of several parameters, including working conditions and pulley geometries and materials, on the thermal flow inside the belt drive system.
Calculation time comparison
The most important achievement of this work is that the model requires considerably less computational resources to accurately predict temperature. Table 4 indicates the times required by this developed method and other two existing methods from current literature. The multiple reference frame (MRF) and dynamic mesh models require hours and days to calculate temperature, while the proposed method only requires less than 5 s on a typical computer desktop with an Intel 4930 K CPU with 3.4 GHz and 32 GB memory. Therefore, this method is applicable for cases that consider large load, such as subsequent pulley material thermal fatigue analysis.
Concluding remarks
A novel analytical thermal model was developed to predict the temperature distributions of complex pulley geometries in a multipulley belt drive system in this work. This study established a model using classical ODEs for complex pulley structures to calculate the heat dissipation of pulleys. Combined with five types of heat generation within the belt drive system, this model can simulate the heat exchange among components and the thermal flows inside components. Moreover, the model was applicable for multipulley belt drive systems and for materials with a wide range of thermal conductivities. This work investigated a fivepulley FEAD system with steel and FRP pulleys as an example. Furthermore, the ODEs and integrations used in the model had analytical solutions; this significantly reduces computational cost and time. The time required for temperature prediction reduced from hours to seconds, thereby allowing for subsequent investigation such as thermal fatigue analysis or analyses with large load. Finally, the results provided by the model were validated by experiments performed on an engine dynamometer. The validation considered various conditions to cover a wide change of speeds and torque loads during engine operation. The predicted and experimental values generally showed good agreement in all conditions and at all points. The maximum difference between the values was less than 6 °C, which was less than 10% of the measured values. This was sufficient for ensuring the validity of the model. Future work will consider the influence of covering cases and the belt wear. In summary, this paper presented an accurate and efficient thermal model for predicting the temperature distribution of a multipulley belt drive system, which had considerable potential for research and industrial applications.
Abbreviations
 A _{1}, A _{2}, A _{3} :

Surface areas of pulley where it has a constant temperature, m^{2}
 A _{ b } :

Surface area of belt exposed to the environment, m^{2}
 A _{ fc } :

Connected area (green area in Fig. 7) between the fin and web, m^{2}
 A _{ exn } :

Contact area between the belt and the nth pulley, m^{2}
 C _{1}, C _{2} :

Constants for Kummer’s functions
 h _{ b } :

Heat transfer coefficient between belt and air, W/(m^{2} K)
 h _{ pn } :

Heat transfer coefficient between the nth pulley and air, W/(m^{2} K)
 h _{ exn } :

Thermal contact conductance coefficient at the nth engaged surface between the belt and the nth pulley in the system
 H _{ f } :

Height of the differential increase in the fin, m
 L _{ f } :

Length of the fin, m
 N :

Number of pulleys inside the belt drive system
 N _{ f } :

Number of fins of the pulley
 P _{ h } :

Heat flux generated per second in the belt in the power loss calculation, W
 P _{ fn } :

Heat flux generated per second by friction at the nth pulley–belt engaging surface in the power loss calculation, W
 R _{ i } :

Inner diameter of pulley, m
 R _{1} :

Inner diameter of first layer, m
 R _{2} :

Outer diameter of second layer, m
 R _{ o } :

Outer diameter of pulley, m
 r :

Differential increase in radius of pulley, m
 T _{ a } :

Ambient temperature of belt drive system, °C
 T _{ b } :

Uniform belt temperature, °C
 T _{ pn } :

Temperature at the outer radius of the nth pulley, °C
 W _{ f } :

Width of a fin, m
 W _{ i } :

Half width of the ith layer of the pulley, m
 x :

Differential increase along the fin, m
 D _{ b } :

Matrix storing the dimensions of the belt in the system
 D _{ pn } :

Matrix storing the dimensions of the nth pulley in the system
 η _{ f } :

Coefficient for fins to increase the heat transfer coefficient, \(h_{p}\), for their connected areas between a fin and the web
 λ _{ pn } :

Coefficient representing the heat dissipation ability of the nth pulley within the system
 Λ_{b} :

Thermal conductivity of belt, W/(m K)
 Λ_{p} :

Thermal conductivity of pulley, W/(m K)
 ξ _{ n } :

Percentage of frictional heat flux flowing into the nth pulley within the system
 τ _{ s } :

Transmitted torque of the belt drive system, N m
 Φ_{Afc} :

Heat dissipation flux per second through the fin connecting area without augmented dissipation on the fin, W/s
 Φ_{b} :

Heat dissipation flux per second from all exposed surfaces of the belt to the environment, W/s
 Φ_{exn} :

Heat exchange per second at the engaging surfaces between the belt and the nth pulley, W/s
 Φ_{fin} :

Heat flux dissipated through the fin outer surfaces per second, W/s
 Φ_{fe} :

Heat flux dissipated through the fin farend surface per second, W/s
 Φ_{fs} :

Heat flux dissipated through the fin side surfaces per second, W/s
 Φ_{pn} :

Frictional heat flux flow into the pulley per second or total heat dissipation from the surfaces of a pulley or the nth pulley in the system per second, W/s
 Φ_{pb} :

Heat dissipated per second at the surfaces on a pulley where it has a constant temperature, W/s
 Φ_{ps} :

Heat dissipated per second at the side surfaces of all three layers on a pulley, W/s
 ω _{ s } :

Transmitted speed of the belt drive system, RPM
References
 1.
Kim, H. S., Bae, H. S., Yu, J. & Kim, S. Y. Thermal conductivity of polymer composites with the geometrical characteristics of graphene nanoplatelets. Sci. Rep. 6, 1–9 (2016).
 2.
Yang, Q. et al. Study of the microclimate and bacterial distribution in the deadspace of N95 filtering face respirators. Sci. Rep. 8, 1–13 (2018).
 3.
Liu, X. & Behdinan, K. Analyticalnumerical model for temperature prediction of a serpentine belt drive system. Appl. Sci. 10, 2709 (2020).
 4.
Gerbert, G. Force and slip behaviour in Vbelt drives. Mechanical Engineering Series 67, (Finnish Academy of Technical Sciences, 1972).
 5.
Manin, L., Liang, X. & Lorenzon, C. Power losses prediction in polyv belt transmissions: application to front engine accessory drives. in International Gear Conference 2014: 26th–28th August 2014, Lyon 1162–1171 (Elsevier, 2014).
 6.
Silva, C. A. F. et al. Modeling of power losses in polyV belt transmissions: hysteresis phenomena (standard analysis). J. Adv. Mech. Des. Syst. Manuf. 11, JAMDSM0085 (2017).
 7.
Silva, C. A. F. et al. Modeling of power losses in polyV belt transmissions: hysteresis phenomena (enhanced analysis). Mech. Mach. Theory 121, 373–397 (2018).
 8.
Bertini, L., Carmignani, L. & Frendo, F. Analytical model for the power losses in rubber Vbelt continuously variable transmission (CVT). Mech. Mach. Theory 78, 289–306 (2014).
 9.
Chen, T. F., Lee, D. W. & Sung, C. K. An experimental study on transmission efficiency of a rubber Vbelt CVT. Mech. Mach. Theory 33, 351–363 (1998).
 10.
Balta, B., Sonmez, F. O. & Cengiz, A. Speed losses in Vribbed belt drives. Mech. Mach. Theory 86, 1–14 (2015).
 11.
Zhu, H., Hu, Y. & Pi, Y. Transverse hysteretic damping characteristics of a serpentine belt: Modeling and experimental investigation. J. Sound Vib. 333, 7019–7035 (2014).
 12.
Kim, H. & Marshek, K. M. Force distribution for a flat belt drive with a concentrated contact load. Mech. Mach. Theory 25, 667–677 (1990).
 13.
Qin, Z., Cui, D., Yan, S. & Chu, F. Hysteresis modeling of clamp band joint with macroslip. Mech. Syst. Signal Process. 66–67, 89–110 (2016).
 14.
Oliver, L. R., Song, G., Shen, Y., Breig, W. F. & Chandrashekhara, K. Accessory serpentine belt stress analysis using hyperelastic model. SAE Tech. Pap. Ser. 111, 1434–1445 (2002).
 15.
Julió, G. & Plante, J. S. An experimentallyvalidated model of rubberbelt CVT mechanics. Mech. Mach. Theory 46, 1037–1053 (2011).
 16.
Zhu, C., Liu, H., Tian, J., Xiao, Q. & Du, X. Experimental investigation on the efficiency of the pulleydrive CVT. Int. J. Automot. Technol. 11, 257–261 (2010).
 17.
Wurm, J., Fitl, M., Gumpesberger, M., Väisänen, E. & Hochenauer, C. Advanced heat transfer analysis of continuously variable transmissions (CVT). Appl. Therm. Eng. 114, 545–553 (2017).
 18.
Wurm, J., Fitl, M., Gumpesberger, M., Väisänen, E. & Hochenauer, C. Novel CFD approach for the thermal analysis of a continuous variable transmission (CVT). Appl. Therm. Eng. 103, 159–169 (2016).
 19.
Gerbert, G. Heat in Vbelt Drives (Lund Technical University, Lund, 1981).
 20.
Abe, S., Tokoro, M., Yaegashi, T. & Ogawa, K. Thermal analysis of timing belt. in SAE Technical Paper Series (1989).
 21.
Merghache, S. M. & Ghernaout, M. E. A. Experimental and numerical study of heat transfer through a synchronous belt transmission type AT10. Appl. Therm. Eng. 127, 705–717 (2017).
 22.
Safaei, B., MoradiDastjerdi, R., Qin, Z., Behdinan, K. & Chu, F. Determination of thermoelastic stress wave propagation in nanocomposite sandwich plates reinforced by clusters of carbon nanotubes. J. Sandw. Struct. Mater. https://doi.org/10.1177/1099636219848282 (2019).
 23.
Wu, Y., Wu, M., Zhang, Y. & Wang, L. Experimental study of heat and mass transfer of a rolling wheel. Heat Mass Transf. 50, 151–159 (2014).
 24.
Krane, R. J., Jischke, M. C. & Rasmussen, M. L. The thermal analysis of a belt type radiator by the method of matched asymptotic expansions. Int. J. Heat Mass Transf. 16, 1165–1174 (1973).
 25.
Singh, S., Negi, J. S., Bisht, S. & Sah, H. Thermal performance and frictional losses study of solid hollow circular disc with rectangular wings in circular tube. Heat Mass Transf. 55, 2975–2986 (2019).
 26.
Zhang, J. Z., Tan, X. M. & Zhu, X. D. Investigation on convective heat transfer over a rotating disk with discrete pins. Heat Mass Transf. 50, 85–94 (2014).
 27.
Kayhani, M. H., Shariati, M., Nourozi, M. & Karimi Demneh, M. Exact solution of conductive heat transfer in cylindrical composite laminate. Heat Mass Transf. 46, 83–94 (2009).
 28.
Merghache, S. M. & Ghernaout, M. E. A. Measuring the stiffness of a timing belt type binder magnetic AT10. in Recueil de mécanique 71–78 (2016).
 29.
McPhee, A. D. & Johnson, D. A. Experimental heat transfer and flow analysis of a vented brake rotor. Int. J. Therm. Sci. 47, 458–467 (2008).
 30.
Belhocine, A. & Bouchetara, M. Thermalmechanical coupled analysis of a brake disk rotor. Heat Mass Transf. 49, 1167–1179 (2013).
 31.
Talati, F. & Jalalifar, S. Analysis of heat conduction in a disk brake system. Heat Mass Transf. 45, 1047–1059 (2009).
 32.
Sellami, A., Kchaou, M., Elleuch, R. & Desplanques, Y. Thermal analysis of padondisc contact under tribological solicitations: a coupled numerical–experimental approach to identify surface temperatures and flow partition coefficient. Heat Mass Transf. 52, 1923–1934 (2016).
 33.
Hennecke, D. K. Thermal analysis of a highpressure compressor rotor of an aeroengine—Venting as a means for life improvement. Wärmeund Stoffübertragung 18, 227–230 (1984).
 34.
Song, G., Chandrashekhara, K., Breig, W. F., Klein, D. L. & Oliver, L. R. Analysis of Cordreinforced polyrib serpentine belt drive with thermal effect. J. Mech. Des. 127, 1198–1206 (2005).
 35.
Chen, G., Lee, J. H., Narravula, V. & Kitchin, T. Friction and noise of rubber belt in low temperature condition: The influence of interfacial ice film. Cold Reg. Sci. Technol. 71, 95–101 (2012).
 36.
Pyr’yev, Y. Analytical solution of thermal conduction in a twolayer cylinder modeling oscillator roller in an offset machine. Int. J. Therm. Sci. 136, 433–443 (2019).
 37.
Liu, J. & Ma, Y. S. 3D levelset topology optimization: A machining featurebased approach. Struct. Multidiscip. Optim. 52, 563–582 (2015).
 38.
Ning, J., Nguyen, V., Huang, Y., Hartwig, K. T. & Liang, S. Y. Inverse determination of Johnson–Cook model constants of ultrafinegrained titanium based on chip formation model and iterative gradient search. Int. J. Adv. Manuf. Technol. 99, 1131–1140 (2018).
 39.
Gerbert, B. G. Pressure distribution and belt deformation in Vbelt drives. J. Eng. Ind. (1975).
 40.
Pellé, J. & Harmand, S. Heat transfer study in a rotor–stator system airgap with an axial inflow. Appl. Therm. Eng. 29, 1532–1543 (2009).
 41.
Lallave, J. C., Rahman, M. M. & Kumar, A. Numerical analysis of heat transfer on a rotating disk surface under confined liquid jet impingement. Int. J. Heat Fluid Flow 28, 720–734 (2007).
 42.
Polikarpova, M. et al. Thermal effects of stator potting in an axialflux permanent magnet synchronous generator. Appl. Therm. Eng. 75, 421–429 (2015).
 43.
Lim, D. H. & Kim, S. C. Thermal performance of oil spray cooling system for inwheel motor in electric vehicles. Appl. Therm. Eng. 63, 577–587 (2014).
 44.
Gerbert.G. Traction Belt Mechanics. (Machine and Vehicle Design, Chalmers Univ. of Technology, 1999).
 45.
Shiming, Y. & Wenquan, T. Heat Transfer (Higher Education Press, Beijing, 1998).
 46.
Barton, D. E., Abramovitz, M. & Stegun, I. A. Handbook of mathematical functions with formulas, graphs and mathematical tables. J. R. Stat. Soc. Ser. A 128, 593 (1965).
 47.
SchwarzenbergCzerny, A. On matrix factorization and efficient least squares solution. Astron. Astrophys. Suppl. 110, 405 (1995).
 48.
Liu, X. & Behdinan, K. A novel analytical model for the thermal behavior of a fiberreinforced plastic pulley in a frontend accessory drive. Adv. Mech. Eng. 12, 168781402092449 (2020).
 49.
ThermoMETER Temperature Sensors Datasheet. https://www.instrumart.com/assets/CSmicrothermoMETERdatasheet.pdf
Funding
This work was supported by Mitacs Canada and Litens Automotive Group.
Author information
Affiliations
Contributions
X.L. investigated the problem, developed the methodology, performed experiments, validations, and wrote the paper; K.B. provided the supervision, and reviewed and edited the final draft.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Liu, X., Behdinan, K. Innovative analytical model for temperature prediction of frontend accessory drive. Sci Rep 11, 1476 (2021). https://doi.org/10.1038/s41598020799865
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598020799865
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.