1. Introduction
1.1. Background
There is an industrial demand in the increased performance of mechanical power transmission devices. This need in high performance is driven by high load capacity, high endurance, low cost, long life, and high speed. For gears, this has lead to development of new designs, such as gears with asymmetric teeth. The geometry of these teeth is such that the drive side profile is not symmetric to the coast side profile. This type of geometry is beneficial for special applications where the loading of the gear is unidirectional. In such an instance, the loading on the gear tooth is not symmetric, thus calling for asymmetric teeth. Since one of the situations that demand high performance is the high rotational speeds, there is a need to understand the dynamic behavior of the gears with asymmetric teeth at such speeds. Such knowledge would shed light on detrimental characteristics like dynamic loads and vibrations. An efficient way in performing studies on the dynamic behavior of gears is using computer aided analysis on numerical models.
A number of studies on the design and stress analysis of asymmetric gears are available in literature. A large number of studies have been performed over the last two decades to assess whether asymmetric gears are an alternative to conventional gears in applications requiring high performance. In these studies, some standards (i.e., ISO 6336, DIN 3990), analytical methods (i.e., the Direct Gear Design method, the tooth contact analysis), and numerical methods (e.g., Finite element method) have been used to compare the performance of conventional and asymmetric gears under the same conditions (Cavdar et al., 2005; Kapelevich, 2000, Karpat, 2005; Karpat et al., 2005; Karpat & EkwaroOsire, 2008; Karpat et al., 2008; Karpat & EkwaroOsire, 2010). In the last ten years, the researches conducted in the area of gears with asymmetric teeth point to the potential impact of asymmetric gears on improving the reliability and performance requirements of gearboxes. The benefits of asymmetric gears which have been offered by researchers are: higher load capacity, reduced bending and contact stress, lower weight, lower dynamic loads, reduced wear depths on tooth flank, higher reliability, and higher efficiency. Each of the benefits can be obtained due to asymmetric teeth designed correctly by designers.
1.2. Dynamic analysis of involute spur gears with symmetric teeth
Gear dynamics has been a subject of intense interest to the gearing area during the last few decades. The dynamic response of a gear transmission system is becoming essential due to increased requirements for high speed, low vibration and heavy load in gear design. However, the numerous design parameters, manufacturing and assembly errors, tooth modifications, etc. make difficult to understand gear dynamic response. The dynamic load reducing in a gear pair may decrease noise, increase efficiency, improve pitting fatigue life, and prevent gear tooth failures. Thus far, many researchers have conducted theoretical and experimental studies on gear dynamics. Most of literature on mathematical models used to predict the gear dynamics have been reviewed by (Ozguven & Houser, 1988; Parey & Tandon, 2003). In these reviews, the theoretical studies use a numerical method which included the excitation terms due to errors and periodic variation of the mesh stiffness. This method was used by many researchers to calculate the dynamic contact load or the torsional response, depending on different gear parameters, i.e., tooth errors, addendum modification, mesh stiffness, lubrication, damping factor, gear contact factor, and friction coefficient.
In dynamic analysis of gears, the dynamic factor and static transmission are the two most important definitions. The dynamic factor is defined as the ratio of the maximum dynamic load to the maximum static load on the gear tooth. Dynamic loads of gears with low contact ratio (between 1 and 2) are affected by several parameters, namely: timevarying mesh stiffness, tooth profile error, contact ratio, friction, and sliding. Static transmission errors, which are defined as the difference between the position of an actual gear tooth and that of an idealized gear tooth, and dynamic loads, affect the gear vibrations, acoustic emissions, tooth fatigue, and surface failure. The static transmission errors change in a periodic manner, due to the variation of gear mesh stiffness during contact. This is the source of vibratory excitation in gear dynamics. The static transmission error has basic periodicities related to the shaft rotational frequencies and the gear mesh frequency. The mesh frequency and its first harmonics are the predominant contributors to the generation of noise. The Fast Fourier Transform (FFT) can be used to perform the frequency analysis of static transmission error.
1.3. Motivation and objectives
Involute spur gears with asymmetric teeth provide flexibility to designers for different application areas due to nonstandard design. If they are correctly designed, they can make important contributions to the improvement of designs in aerospace industry, automobile industry, and wind turbine industry. This often relates to improving the performance, increasing the load capacity, reduction of acoustic emission, and reduction of vibration. In the past, most of the analysis of gears with asymmetric teeth has been limited to cases under static loading.
Dynamic loads and vibration are a major concern for gears running at high speeds. Therefore, dynamic behavior should be analyzed to determine the feasibility of asymmetric gears in different applications. In order to utilize asymmetric gear designs more effectively, it is imperative to perform analyses of these gears under dynamic loading. This study offers designers preliminary results for understanding the response of asymmetric gears under dynamic loading. The effect of some design parameters, such as pressure angle or tooth height on dynamic loads, is shown. The asymmetric gears considered will have a larger pressure angle on the drive side compared to the coast side. In this study, to investigate the response of asymmetric gears under dynamic loading, the dynamic loads and static transmission errors were used. The first objective of this chapter is to use dynamic analysis to compare conventional spur gears with symmetric teeth and spur gears with asymmetric teeth. The second objective is to develop a MATLABbased virtual tool to analyze dynamic behavior of spur gears with asymmetric teeth. For this purpose a MATLAB based virtual tool called DYNAMIC is developed.
The first part of the study is focused on assymetric gear modelling. The second part focuses on the virtual tool parameters. In the third and the last part, the simulation results are given for different asymmetric gear parameters.
2. Dynamic model for involute spur gears with asymmetric teeth
There is an essential need to find the equations of motion for a gear tooth pair during a mesh to determine the variation of dynamic load with the contact position. A singledegreeoffreedom model of the gear system consists of a gear and a pinion shown in Error! Reference source not found.. The equations of motion can be expressed as follows:
where
The static tooth load is defined as:
The relative displacement, velocity, and acceleration can be writtehn as follows:
The effective gear masses are:
Including viscous damping, the equations of motion are reduced to:
The loaded static transmission errors can be obtained by dividing Eq. (8) by Eq. (7) to yield:
The equivalent stiffness of meshing tooth pairs, in Eq. (7) through (9), can be written as:
The friction experienced by the pinion and the gear can be expressed as:
The signs in the above expressions are positive (+) for the approach and negative (−) for the recess.
The coefficient of friction is expressed by formula:
where ν is the viscosity of lubricant (cSt). And
where
The value of the damping ratio, ξ, in Eq. (6), is commonly recommended in literature as one between 0.1 and 0.2. In this study, a constant value of 0.17 proposed in literature for the damping ratio, ξ, was adapted in the solution of equations.
The dynamic contact loads, which include tooth profile error, can then be written as:
where λ_{I} and λ_{II} are the tooth profile errors. In this study, the effects of profile errors on the dynamic response of gears are not considered. Thus, the tooth profile errors are assumed to be zero. The developed computer program has a capability of using any approach for the determination of errors.
It should be noted that the above equations are valid only when there is contact between two gears. When separation occurs between two gears, because of the relative errors between the teeth of gears, the dynamic load will be zero and equation of motion will be given by:
The meshing conditions are described as follows:
If 

Double tooth contact 
If 

Tooth separation 
If λ_{I }< 

Single tooth contact 
If λ_{II }< 

Single tooth contact 
3. Tooth stiffness
According to Equations (10) and (11), in order to calculate the equivalent stiffness of a meshing tooth pair, the tooth stiffness has to be known beforehand. In this study, a 2D finite element model was developed to calculate the deflections of both the asymmetric and the symmetric gear teeth. By using this model, nodal deflections are calculated for predetermined contact points. The load applied for each contact point is taken as a constant in order to determine tooth deflection under unit load. By putting the calculated nodal deflection values into Equations (2427), the tooth stiffness are calculated and then the approximate curves for the single tooth stiffness along the contact line are obtained with respect to the radius of the gears. This process was repeated for each gear previously designed for different gear parameters.
where
4. Computational procedure
The reduced equation of motion is solved numerically using a method that employs a linear iterative procedure. This involves dividing the mesh period into many equal intervals. In this study, the flowchart of this computational procedure developed in MATLAB, used for calculating the dynamic responses of spur gears, is shown in Fig. 2. The time interval, between the initial contact point and the highest point of single contact, is considered as a mesh period. In the numerical solution, each mesh period is divided into 200 intervals for good accuracy. Within each of the subintervals thus obtained, various parameters of equations of motion are taken as constants, and an analytical solution is obtained. The calculated values of the relative displacement and the relative velocity after one mesh period are compared with the initial values x_{r} and v_{r}. Unless the differences between them are smaller than a preset tolerance (0.000001), the iteration procedure is repeated by taking the previously calculated values of
5. DYNAMIC virtual tool
Physicsbased modeling and simulation is important in all engineering problems. The current mature stage of computer software and hardware makes it possible complex mechanical problems, such as gear design, to be solved numerically. Inhouse prepared codes to handle individual research projects, graduate, and/or PhD studies; commercial packages for engineers in industry are widely used to solve almost every engineering problem. Tailored with graphical user interfaces (GUIs) and easytouse design steps, anyoneeven a beginner can design a gear pair and obtain results, e.g Dynamic Load, Transmitted Torque, Static Transmission Error as a function of time, and Static Transmission Error Harmonics etc., just by pressing a command button. Lecturers have been increasingly using these packages to increase their teaching performance and student understanding.
Based on and triggered by these thoughts, a virtual tool DYNAMIC is prepared that can be used for educational and research purposes. The DYNAMIC is a general purpose gear analyzing tool (Fig. 3).
There are six blocks and a figure block on the front panel of the tool. Three blocks on the right side of the front panel, belong to the parameters which will be defined by the users (Fig. 2 a, b).
The two blocks above the figure are
The variation of dynamic load with respect to time can be seen in Fig. 6. The solutions for different variables can be plotted in one figure, for comparison. In Fig. 6 two different solutions for dynamic load are plotted for different revolution speed. Fig. 7 is an example for Transmitted Torque solution.
6. Results and discussions
The computer program developed has been used for the dynamic analysis of spur gears with symmetric and asymmetric teeth. In this study, seven different gear pairs are considered for the dynamic analysis of spur gears with asymmetric teeth. In order to simplify the analysis, all gear parameters are kept constant, apart from the pressure angle on the drive side and the tooth height. Since the effects of the tooth profile errors are not considered in this study, the analyzed gears are assumed to be “perfect gears” without tooth errors. The properties of these gear pairs are provided in table
In a previous work (Karpat, 2005), different approaches for minimizing the dynamic factors and the static transmission errors, in lowcontact ratio gears, were reviewed in details. In one of the approaches discussed, the usage of high gear contact ratio was included. It was observed that increasing the gear contact ratio reduced the dynamic load. In literature, minimum dynamic loads were obtained for contact ratios between 1.8 and 2.0. A way of increasing the contact ratio is by using higher addendum values. It should be noted that increasing the value of the addendum leads to a reduction in the bending stress at the tooth root. This occurs through the lowering of the location of the highest point of single tooth contact (HPSTC). The other gear characteristics impacted by high addendum are the thickness of tooth tip and undercut. In this study, for asymmetric gears, high addenda are analyzed, as a means of minimizing the dynamic factors and the static transmission errors (Gear Pair 4 and 5).









2 mm  2 mm  2 mm  2 mm  2 mm 

20  20  32  32  32 

20°  20°  20°  20°  20° 

20°  24°  32°  24°  32° 

2  2  2  2  2 

1 kg  1 kg  1 kg  1 kg  1 kg 

2 kg  2 kg  2 kg  2 kg  2 kg 

Steel  Steel  Steel  Steel  Steel 

100 cSt  100 cSt  100 cSt  100 cSt  100 cSt 

0,17  0,17  0,17  0,17  0,17 

20 mm  20 mm  20 mm  20 mm  20 mm 

1 m_{n}  1 m_{n}  1 m_{n}  1.32 m_{n}  1.17 mn 

1.64  1.49  1.31  1.90  1.52 
For the sample gear pair whose dimensions and properties are given in Table 1, variations of dynamic loads are determined for various pinion speeds between 1000 rpm and 20 000 rpm. As an example, the dynamic load variation of gear pair 1 for 1000 rpm, 3000 rpm, 10 000 rpm and 18 000 rpm is shown in Figure 8.
Fig. 9 shows the relationship between the dynamic factors and the rotational speed. When comparing the maximum dynamic factors in the corresponding gear pairs in Fig. 9. (e.g., Gear Pair 1 versus Gear Pair 3), it is generally stated that the dynamic factor for spur gears with asymmetric teeth increases with increasing pressure angles on the drive side. Furthermore, it is obvious that the sample Gear Pair 4, which is the gear pair with the highest gear contact ratio 1.90, has a lower dynamic load, at all speeds; this indicates that the impact of gear contact ratio on dynamic loads. The highest dynamic factor is observed at the resonant rotational speed (about 12 000). Beyond this speed, the asymmetric teeth have consistently higher dynamic factors than symmetric teeth. One of reasons for that may be the effect of contact ratio on dynamic loads. As the pressure angle on drive side increases, the contact ratio decreases. However, the dynamic factor in gear systems decreases with increasing the contact ratio. This result may be due to the narrow single contact zone. Because of the narrow single contact zone, this zone is passed speedily as gear rotate and system can not respond. Other reason may be seen by analyzing the variation of mesh stiffness with respect to time. As can be seen from this figure, in the single contact zone, the asymmetric gear (Gear Pair 4) has higher mesh stiffness than the symmetric gear (Gear Pair 1). The high mesh stiffness is one of the reasons for the high dynamic factor observed in Fig.9.
Fig. 10 shows the impact of increasing the pressure angle, on the drive side, on the static transmission error. Generally, changing the pressure angle will impact the tooth mesh characteristics, such as the tooth contact zone and contact ratio. Error! Reference source not found.1 indicates that the single tooth contact zone increases with increased pressure angle. Thus, compared to gears with symmetric teeth, gears with asymmetric teeth have a larger single tooth contact zone.
Furthermore, the static transmission error, at the center of the single tooth contact zone, decreases with increasing of pressure angle. The frequency spectra of the static transmission errors are depicted in Error! Reference source not found. 11. In these figures, the sum of first five harmonics slightly increases with increasing pressure angle.
Fig. 12 (d) and (e) shows the static transmission error for increased values of addendum for asymmetric teeth. Increasing the addendum, the amplitude of the static transmission errors is decreased for a comparable pressure angle. Additionally, the single tooth contact zone decreased for a comparable pressure angle. In Fig. 12 (d), it is noted that the asymmetric tooth with
7. Conclusions
Virtual tools have become very effective in teaching engineering problems. The user need not to know graphical user interface details, programming tips, etc. Instead, such tools have the capability of handling a variety of different gear problems.
A MATLABbased virtual tool, DYNAMIC, is introduced to analyze dynamic behavior of spur gears with asymmetric tooth design. The DYNAMIC is used to compare conventional spur gears with symmetric teeth and spur gears with asymmetric teeth in this study. The results for dynamic load, dynamic factor, transmitted torque, static transmission error and static transmission error harmonics are obtained for various tooth parameters to show the powerful aspects of asymmetric teeth.
References
 1.
Ozguven HN Houser DR 1988 The mathematical models used in gear dynamics A review.  2.
Parey A. Tandon N. 2003 Spur gear dynamic models including defects: A review.  3.
Cavdar K. Karpat F. Babalik F. C. 2005 Computer aided analysis of bending strength of involute spur gears with asymmetric profile.  4.
Kapelevich A. 2000 Geometry and design of involute spur gears with asymmetric teeth.  5.
Karpat F. 2005 Analysis of Involute Spur Gears with Asymmetric Teeth. Ph.D Thesis [BursaTurkey: Uludag University.  6.
Karpat F. EkwaroOsire S. 2010 Influence of Tip Relief Modification on the Wear of Spur Gears with Asymmetric Teeth,  7.
Karpat F. EkwaroOsire S. Cavdar K. Babalik F. C. 2008 Dynamic Analysis of Involute Spur Gears with Asymmetric Teeth,  8.
Karpat F. EkwaroOsire S. 2008 Influence of Tip Relief Modification on the Wear of Spur Gears with Asymmetric Teeth,  9.
Karpat F. EkwaroOsire S. Khandaker M. P. H. 2008 Probabilistic Analysis of MEMS Asymmetric Gear Tooth,