Microsoft Word - numero_30_art_57 G. Belingardi et alii, Frattura ed Integrità Strutturale, 30 (2014) 469-477; DOI: 10.3221/IGF-ESIS.30.57 469 Focussed on: Fracture and Structural Integrity related Issues Dynamic additional loads influencing the fatigue life of gears in an electric vehicle transmission G.Belingardi, V. Cuffaro, F. Curà Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Torino (Italy) giovanni.belingardi@polito.it, vincenzo.cuffaro@polito.it, francesca.cura@polito.it ABSTRACT. In recent years the implementation of the electric engine in the automotive industries has been increasingly marked. The speed of the electric motors is much higher than the combustion engine ones, bringing transmission gears to be subjected to high dynamic loads. For this reason the dynamic effects on fatigue life of these components have be taken into account in a more careful way respect to what is done with the usual gears. In the present work the overload effects due to both speed and meshing in a gear couple of an electric vehicle transmission have been analyzed. The electric vehicle is designed for urban people mobility and presents all the requirements to be certified as M1 vehicle (a weight less than 600 kg and a maximum speed more than 90 Km/h). To investigate the overload effects of teeth in contact, the reference gear design Standards (ISO 6336) introduce a specific multiplicative factor to the applied load called Internal Dynamic Factor (Kv). Aim of this work is to evaluate how dynamic overloads may influence the fatigue life of the above quoted gears in term of durability. To this goal, Kv values have been calculated by means of the analytical equations (ISO 6336 Methods B and C) and then they have been compared with the results coming from multibody simulations, involving full rigid and rigid-flexible models. KEYWORDS. Multibody simulations; Electric vehicle; Gear design. INTRODUCTION n recent years, greater emphasis has been placed on the design of high-speed, lightweight, precision systems. Both performance analysis and design of such systems can be greatly enhanced through transient dynamic simulations. The need for a better design has made necessary the insertion of many factors that have not been considered in the past. Systems such as engines, robotics, machine tools and space structures may operate at high speeds and in very high temperature environments. In particular, the implementation of the electric engine in the automotive industries has been increasingly marked. The speed of the electric motors is much higher than the combustion engine ones, bringing transmission gears to be subjected to high dynamic loads. For this reason the dynamic effects on fatigue life of these components have be taken into account in a more careful way respect to what is done with the usual gears. I G. Belingardi et alii, Frattura ed Integrità Strutturale, 30 (2014) 469-477; DOI: 10.3221/IGF-ESIS.30.57 470 So, the behaviour of new automotive transmissions may be simulated with specific multibody programs, being these systems suitable to be described as a collection of subsystems called bodies and the motion of them is kinematically constrained. Basic to any analysis of multibody mechanics is the understanding of the motion of subsystems (bodies or components). Referring to this environment, a gear is a complex body with a specific geometry, as an example the modulus, the number of teeth, the face width and the contact ratio. Once a gear is designed as a function of the power to be transmitted, it is necessary to analyse it from both kinematic and dynamic point of view [1]. Some papers have been published about the numerical analysis of the dynamic meshing gears. Most of them are based on lumped parameter models [2-5], where the teeth contact is represented by an equivalent mass, a stiffness and damping of the teeth. Another way to address the contact analysis is related to FEM simulations; this approach is helpful to determine the gear contact stiffness and the corresponding strain and stress distributions in the teeth. In particular the multibody kinematic approach uses a contact stiffness which depends on parameters that are not yet well understood [6]. In the present work the overload effects due to both speed and meshing in a gear couple of an electric vehicle transmission have been analyzed. The electric vehicle is designed for urban people mobility and presents all the requirements to be certified as M1 vehicle (a weight less than 600 kg and a maximum speed more than 90 Km/h). The first aim of this research is to evaluate how the dynamic overloads may influence the fatigue life of the above quoted gears in term of durability. To this goal, a specific multiplicative factor to the applied load, called Internal Dynamic Factor (Kv), involved in ISO 6336 Standards [7-10] has been taken into account; Kv values have been calculated by means of the related analytical equations (ISO 6336 Methods B and C [7]).Then, the dynamic analysis of the above quoted two gears has been performed by means of the multibody software RecurDyn (Functionbay), in order to obtain the real response of the components. The model has been developed by tuning step by step the contact parameters of engaging pair of gears and by optimizing the integration parameters. Forces values obtained by the simulations have been utilised to determine the dynamic factor Kv involved in ISO Standard 6336 [7] for the calculation of gears load capacity. Analytical results in terms of Kv values have been compared with those coming from multibody simulations, involving full rigid and rigid-flexible models. GEAR MODEL bject of this paper consists of two specific gears of a differential architecture for an electric vehicle. Its design has been developed in order to keep as low as possible the transmission weight; in this work, a single speed stage gearbox has been chosen. Transmission is constituted by an ordinary gear system and an epicyclical one named differential (see Fig. 1). The ordinary gear system consists of four helical gears, the differential of four bevel gears (two pinions and two suns). From the dynamic point of view, only the ordinary gears system have to be taken into account; the transmitted torque passes directly to the wheel axes through the pinion and the vehicle has been analysed in a rectilinear motion condition. Fig. 1 shows the complete transmission model that has been set up for the simulation. Figure 1: complete transmission model. O G. Belingardi et alii, Frattura ed Integrità Strutturale, 30 (2014) 469-477; DOI: 10.3221/IGF-ESIS.30.57 471 The main design data of the four helical gears are reported in Tab. 1 (engagement 1, Gears 1 and 2; engagement 2, Gears 3 and 4). All gears have been calculated according to ISO 6336 Standard [7-10]. In this work the engagement 2 has been considered (gears 3 and 4). Gear 1 Gear 2 Gear 3 Gear 4 Pressure angle α 20 20 20° 20° Helix angle β 30 (L.H.) 30 (R.H.) 30° (R.H.) 30° (L.H.) Normal module mn 1.5 1.5 2 mm 2 mm Number of teeth Z 22 68 34 75 Face width b 42 42 30 mm 30 mm Table 1: Main design dimensions of gear’s engagement 1 (input drive) and engagement 2 (final drive) MULTIBODY MODEL he simulation of meshing gears has been analysed (see Fig. 2). To obtain a correct simulation of this complex system, a lot of preliminary tests has been done in order to tune both kinematic and dynamic parameters. The final goal is to achieve a good matching between the multibody response and the corresponding theoretical one. In this paper two different type of contact have been realized; the first (Fig. 2 on the left) consists of two rigid bodies in contact, the second (Fig. 2 on the right) consists of two rigid-flexible bodies, simulating the gears web as a rigid body and the rim as flexible. Since the system is very complex, the computational time is very high. So, it was necessary to properly calibrate the simulation before its starting. The best compromise between parameter’s calibration and the computational time has been searched. Contact parameters have been changed from the default ones because the consequent increase in computation time is widely justified by the obtained improvement in the quality of the results. Figure 2: Multibody automotive transmission system, full-rigid on the left and rigid-flexible on the right. GEARS INTERNAL DYNAMIC FACTOR tandard design procedures compare the gears strength with a calculated stress value, in both bending and pitting cases. The most common approach used in calculation of these stress values, as widely described in ISO Standard 6336 part 1 [7], involves the use of influence factors, derived from results of research and field service. The influence factors may be distinguished between two categories, factors which are determined by gear geometry or which have been established by convention and factors which account for several influences and which are dealt with as independent of each other; in this last group is included the Internal Dynamic Factor Kv. As a matter of fact, even when the input torque and speed are constant, significant vibration of the gear masses and resultant dynamic tooth loads can exist. These loads result from the relative displacements between the mating gears as they vibrate in response to an excitation known as transmission error. T S G. Belingardi et alii, Frattura ed Integrità Strutturale, 30 (2014) 469-477; DOI: 10.3221/IGF-ESIS.30.57 472 The internal dynamic factor makes allowance for the effects of gear tooth accuracy grade as related to speed and load. High accuracy gearing requires less derating than low accuracy gearing.It is generally accepted that the internal dynamic load on the gear teeth is influenced by both design and manufacturing. “Perfect” gears are defined as having zero quasi-static transmission error at the nominal transmitted (design)mesh torque. They can only exist for a single load and, with proper modifications, have zero dynamic effectszero transmission error (perfect conjugate action), zero excitation, no fluctuation at tooth mesh frequencyand no fluctuation at rotational frequencies. With zero excitation from the gears, there is zero response at any speed [7]. The Internal Dynamic Factor Kv takes into account the effects due to the rotating masses; ISO 6336 part 1 [7] suggests three methods for calculating this factor. Method A (Kv-A) derives from the results of full scale load tests, precise measurements or comprehensive mathematical analysis of the transmission system and all gear and loading data shall be available, then this method, in this work, corresponds to the dynamic multibody analysis results. Method A generally results the most sophisticated. Method B (Kv-B) is suited for all types of transmission, spur and helical gearing with any basic rack profile and any gear accuracy grade and, in principle, for all operating conditions. Method C (Kv-C) supplies average values which can be used for industrial transmissions and gear systems with similar requirements, with restriction in the application field. In the present paper Methods B (Kv-B) and C (Kv-C) have been taken into account as the resolution of the corresponding equations described in detail in [7]. Kv-B coefficient has been calculated on the basis of the different operating ranges (subcritical, main resonance, intermediate and supercritical ranges) related to the resonance ratio N of the mating gears [7]. Firstly, the resonance running speed of the gear pair nE1 has been determined as a function of the reduced gear pair mass per unit face width, of the mesh stiffness and of the pinion number of teeth, as indicated in detail in [7]. Then the resonance ratio N, where N=n1/nE1, has been calculated, being n1 is the rotational speed of the pinion in rpm. Once the resonance ratio N has been obtained, the operating range has been determined as a function of the specific load Ft· KA / b, according to [7],where Ft is the tangential load, KA the application factor, b the gear face width and Ns is the lower limit of the main resonance range. The dynamic factor Kv-B has been computed for the different ranges, following the involved relationships indicated in [7]; non-dimensional parameters which take into account the effect of tooth deviations and profile modifications on the dynamic load are considered. In particular, it has been calculated for: subcritical range N ≤ NS (the majority of industrial gears operate in this range), main resonance range NS