Mechanics & Industry
Volume 21, Number 5, 2020
|Number of page(s)||12|
|Published online||10 August 2020|
Contact pressure in misaligned spline couplings
Haute Ecole d'Ingénierie et de Gestion du Canton de Vaud, HEIG-VD, 1 route de Cheseaux, 1400 Yverdon-les-bains, Suisse
* e-mail: firstname.lastname@example.org
Accepted: 9 June 2020
Spline couplings are increasingly being studied because of their growing use, particularly in the field of electric cars. Due to clearances, geometrical obligations and elastic deformation, spline couplings operate mainly in slightly misaligned situations. From these kinds of situations, problems could appear such as vibrations or contact damages. The purpose of this article is to analyse the misalignment of splined couplings and the effects it induces. In order to identify behavioural trends, this work offers analytical approaches for estimating contact pressures in spline connections. In the chosen approach, the simplifying hypotheses necessary for equations are detailed as well as the various analytical modelling. Despites these assumptions, comparisons with pure numerical approaches illustrates the quality of analytical approach adopted in this paper. This method of modelling contact pressure allows the identification of the main mechanical factors to be identified including stiffness, thus deceasing problems due to misalignments.
Key words: Spline couplings / misalignment / contact pressures / modeling / stiffness
© AFM, EDP Sciences 2020
Spline coupling is a current component used to transmit a rotating torque from an input (a shaft) to an output (a gear hub). The torque transmission is done thanks to external teeth on the shaft engage with the internal teeth in the hub.
Use areas of splined joint include aerospace, manufacturing, wind power and automotive. In this last domain, current emergence of electrical cars will increase a lot its use and do necessary to get a better understanding of its specific behaviors.
From a mechanical point of view, the main reasons for using such component are its ability to transmit large torque (better mechanical resistance than keyed shaft), to operate at high speed, to accept reverse operating conditions, ease of assembly or disassembly and maintenance. Due to its ability to accommodate some slight misalignment, to accept some axial displacements, this simple and inexpensive joint do not need accurate positioning of the coupled shafts.
Despite its various characteristics, often spline couplings are over dimensioned and remain subject to damages or failures. The most common failure or damage mode observed in spline couplings are wear, fretting-wear , fretting-corrosion, fretting fatigue, teeth fatigue  and over loading of spline coupling shaft bearing . Some noise and vibration problems have also been encountered in industrial applications .
Optimal design and use of spline couplings are always tricky because it is quite difficult to assess contact behavior and load distribution characteristics , particularly when spline teeth have manufacturing errors. Numerical simulations are very hard and require specific codes or development integrating Boundary Element method. First developments with boundary method are twenty years old and have been done by Adey .
From this quick overview, it is clear that it is a great challenge to design a coupling joint exactly as is necessary in the early stages of design. To go deeper in detail, we could summarize to four main problems to address to this challenge.
In the first instance, the real operating conditions are never well known, particularly, when the misalignment and geometrical defects of teeth are present. As mentioned by Alfares  “we have to live with”. Curà and Mura  show experimentally the difficulty to correctly evaluate misalignment in true conditions. On the one hand, the importance of the defects is underestimated. On the other hand, the tempering effect of the contact elasticity and of the structure is not properly evaluated.
Secondly, the criteria for ranking the operating conditions of the coupling in relation to the four major damages are not efficient enough to be perfectly predictive. For wear, it is usual to use Archard relationship. But it is well known that this one needs to be adapted and adjusted experimentally as underlined by Fillot . Concerning fretting, and its various forms of damage, there is currently no unanimous criterion. For fretting fatigue, be it the Ruiz criterion or the SWT criterion, experimental calibration is necessary to have realistic predictions . For fretting wear this type of comment is confirmed by the large number of publications related to experimental results on this theme [1,2,11–13]. For teeth fatigue, except the fact multiaxial criteria needs to be used (Dang Van, Sines, Crossland, etc.), the difficulty lies in having access to the right level of stress.
Thirdly, the calculation of internal stresses in the teeth of this highly hyperstatic contact problem needs specific code. With commercial codes, for instance, it is almost impossible to simulate contact pressures and stresses in the contact areas between tooth with a good level of accuracy when the spline coupling is misaligned. Only a few authors, having developed specific calculation tools, have published on this subject, including Hong and al. [5,14] and Medina and Olver [15,16].
Finally, for dynamic analysis, the global stiffness of the spline coupling is misaligned condition is generally overestimated. Due to its self-alignment capability, the spline coupling in a misaligned situation generates high reaction moments and forces at the shaft supports. These forces, especially on insufficiently preloaded roller bearings, intensify their nonlinear behavior and the generation of vibration over-harmonics [17,18].
Each of the different aspects mentioned above is the subject of research. However, for the engineer who has to design a spline coupling, it is important for him to have a simpler global approach that can eventually constitute his know-how, once this approach has been correlated with the feedback from experience. The present work attempts to be part of this process. For that, we have chosen to use an analytical approach to fully understand the behavior sensitivities regarding the key parameters defining the tooth geometry of the spline coupling. This choice allows us to keep in mind the main mechanical aspects which govern the peaks of contact pressure between teeth. On the other hand, to get a full understanding of this subject interest, it is useful to highlight the fact that there are very few publications of analytical results [7,19–21]. The main and accurate results on misaligned spline coupling were obtained with dedicated simulation tools that are not commercially available [5,14], however, these results could be used as a reference for our approach. Having simple tools allowing evaluating contact pressure peaks, stiffness, slipping, support reactions would be an important help for the designer. The analytical approach makes it possible to provide answers of this type.
By adopting this approach, this implies that several simplifying hypotheses are necessary to treat the subject. These hypotheses will be presented below.
First, it was chosen to model misalignment from a geometric perspective. To this end, some definitions concerning the geometry of splined connections have been laid:
The material is elastic.
Misalignment angle is very low ≪ 1∘.
Spline coupling clearance is supposed very small.
Coulomb friction is taken into consideration.
The profile of the teeth with developing circle follow the norm NF ISO 4156-1.
The rotation speed on the splined shaft, written ω is constant.
In this paper, misalignment will be characterized by:
Angular misalignment, the acute γ angle between the axes of rotation of the female and male parts.
The parallel misalignment of δ p (offset), that is, the distance between the shaft axis and the female axis at the centre of the common part of the male and female teeth.
Thus, from the angle γ and the active length of the tooth L (Fig. 1), it is possible to express the maximum height δ max of the toothing end in the case of a misalignment. For any gamma angle misalignment and a parallel δ p offset on the y-axis, the maximum local misalignment will be:(1)
Simplified representation of misalignment.
In this section we construct the modelling of the geometric interference between the teeth of the male and female parts. These interferences are directly related to contact pressures between teeth.
Thus, as a first step, we have chosen to assume that there are internal loadings on all teeth. These efforts are unknown, but we can express the produced resulting torsor (reduced torsor at the centre of the splined bond). Forces and moments acting on the toothing in the presence of misalignment are considered. Shear loadings result from the sliding rubbing between the teeth. Each contact loading has been assumed to act on the primary diameter for the radial direction and are identified by the variable din for their axial positioning. As for the direction of these efforts in the plan (y, z) and their position, this is defined in Figures 2 and 3 for a spline coupling having n max teeth.
Relationships for moment and force resultant at the point O can be expressed with reference to the axis x, y, z. Therefore, we consider the angle of offset between the teeth θ = 2π/n max and μ the kinematic friction coefficient.
With the representation given in Figure 4, we obtain the following expression for the torsor:(2) (3) (4) (5) (6) (7)
In this contact problem, it can be noted that this is a highly hyperstatic problem which requires consideration of the deformation aspects of the teeth through stiffness. Stiffness modelling is given below.
By plotting this relative speed on a circular graph, we obtain the speed profile represented by Figure 5. This representation illustrates the friction acting on the shaft as friction forces will have the same orientation than the slippage.
Force application distances.
Geometrical variables and loads on the shaft.
Forces and torques on the splined shaft.
3D Representation of the slippage profile.
From splined gear literature  it seems more appropriate to look successively at how the contact pressures vary axially (according to x), then along the profile (according to the curvilinear abscissa s), as is presented in Figure 6.
Contact pressure application area.
To analyse the axial evolution, we have used a common modelling from literature [23–25]. Thus, the contact lineic force in the axial direction can be expressed as the torque variation over an infinitesimal interval along the x axis.
Therefore, on the assumption that all the teeth of the coupling are in perfect contact and that the forces can be linearly summed, we can investigate the equilibrium of a slice of thickness dx according to the diagram in Figure 7.
The K notation is used to define the torsional stiffness of teeth contact per unit width. This stiffness is multiple of the resulting stiffness for two antagonistic teeth of a shaft/hub torque in contact.
By defining the factors K i = K/G i I i and , the solution of this system can be simplified. The boundary conditions use for solving the differential equations are:
For x = 0 : T 2 = 0 et T 1 = T t .
For x = L : T 1 = 0 et T 2 = T t .
For 0 ≤ x ≤ L : T t = T 1 (x) + T 2 (x).
Due to the hypothesis of perfect teeth contact we can analyse the loading on each tooth.
With an axial lineic contact force P a (x) dx applied to each of the teeth of the fluted shaft, as represented in Figure 7, the torsional balance of these two slices allows the expression of axial lineic force:(16)
The full development can be found on the literature .
Representation of an infinitesimal portion of the splined connection.
To confirm the pertinence of this approach we have compared our results with results from the literature. One of the more accurate result has been produced by Hong . The contact pressure field on shaft tooth extracted from this work is presented in Figure 8.
From these results, with an integration on the curvilinear abscissa, we can get the lineic contact force along the axis x. Taking in consideration the characteristics of the Hong modelling in our own modelling and considering tooth stiffness taken at 1.04 × 109N/rad, lineic force from the two approaches are very close (Fig. 9).
Pressure field on the toothing with a Torque applied T = 2200 Nm.
Comparison of the two model.
As several authors have considered  a punch model to describe the contact pressure variation along the curvilinear axis of the teeth, we shall adopt the same approach. For each tooth in contact, deformation results from a translation and a rotation (Fig. 11), then it seems adequate to consider the contact pressure given by Sackfield  for a punch subject to normal load and a moment (Fig. 10).(18)
With the factor and G come from the material and are respectively the Kolosov constant (χ = 3 − 4v) and the shear modulus G.
Since this model gives an infinite contact pressure at the ends of the contact area, we introduce some corrections proposed by Barrot  with a ' > a and P in > P to get finite realistic value of pressure. The expression of the pression then becomes:(19)
In summary, a soon as the lineic contact force is determined, maximum pressure can be extrapolated from the contact expression of the corrected punch model.
Representation on our model.
Because of the multiple hyperstatic contacts, we described the teeth stiffness per unit of length . In this approach we regard the stiffness as independent of the position x on the teeth. Thus, we consider that this stiffness can be split in three main contributions with their mathematic expressions.
First, the embedment stiffness can be described by the following displacement due to a force F applied on the primary radius : (20)
Second, the displacement due to contact stiffness can be expressed with Hertz relationships :(21)
With the semi-width of the Hertzian contact.
The different characteristics used in these formulas are relatives to Figure 12.
With k sr the reduced cross-sectional coefficient in shear.
Thus, with the geometrical characteristics from Hong model, the stiffness is:
K t = 0.99 × 109N/rad.
This value obtained from elasticity consideration is in perfect agreement with the previous one, K = 1.04 × 109N/rad, obtained by curve fitting.
According to the ψ angle of the tooth (angular position in the plane (y, z)) and according to the face that is in contact (right or left), the contact area is on the front or back of the tooth. For modelling, it is necessary to establish a relation between the γ angular misalignment angle and the normal misalignment of each tooth n. For a coordinate point (x, θ) located on the fictitious primitive cylinder of the shaft, the misalignment according to the tangential direction is expressed as follows:(24)
We then consider the profiles of the spline sections in Figure 13. On it, the profile of the shaft perfectly aligned with hub profile is represented with the green line. The misaligned shaft appears in blue line, while the profile of the misaligned shaft subjected to torque is represented in red line. From this representation, we can express the condition of contact between teeth by the following relation:(25) (26)
Misalignment and torque situations.
Axial lineic forces.
Thus, for a left or right contact of an index tooth n whose one side touches its counterpart, the contact linear force developed according to the normal left or right contact direction, for a portion of tooth dx can be expressed as:(29) (30)
Using Figure 7, it is possible to express the torsional balance, which enables the relationship between the sum of the linear forces and the torque derivative to be obtained as follows:(31)
Considering the complexity of this model, we use a numerical solver, and did some manual iteration according to the number of unknown variable, to get a result with the different expressions and the Hong's model characteristics . It corresponds to introduce the following parameters on our model:
Thus, we find the following profiles of lineic forces on the 25 teeth (Fig. 14) where the lineic contact force is plotted along the vertical axis with respect to the horizontal x axis.
From these lineic forces results, we calculate the corresponding torsors with the relationships (2)–(3)–(4)–(5)–(6)–(7) and find respectively T = 2197 Nm, M F = 108 Nm and M T = − 561 Nm for the torques and N x = 12.7 kN, N y = − 50.5 kN and N z = − 38.3 kN for the forces. This means that to get an angular misalignment, the hub and the shaft must withstand extra moment and load. In the present case, with a torque of 2200Nm and a misalignment of 0.12∘, an extra moment of 560Nm is generated (≃25% of the torque). By knowing the nominal stiffness of the bearings supporting the shaft and the hub it is possible to estimate what will be the true misalignment.
In comparison, the results of Hong have the following form (Fig. 15).
A detail of the comparison of our model with Hong's results is presented in Figure 16 for the maximum force peaks.
Because of the importance of moments and reaction forces, and in order to confirm the benefits of our approach, a comparison was made with Hong's reference results using a 0.12° misalignment angle. We choose the largest angle of misalignment because it induces a broader dispersion on our/the results. After all computations, the comparison displays a good correlation, as illustrated in Figure 17.
Comparison between Hong peak results and our model.
Comparison between Hong results and our model.
For the design engineer, defining an optimal spline coupling in terms of dimensions remains a real challenge because there are no commercial codes allowing him to simulate the behavior of this joint in a misaligned situation. However, any spline coupling works more or less in a misaligned situation and it is in this type of situation that problems can arise such as the destruction of the component. In this work modelling the contact between the teeth of the shaft and the hub of the spline coupling, we have deliberately focused on an analytical approach. The main target was to provide answers that are easy for a design engineer to integrate in his design process. With this analytical approach, it has been shown that this type of approach makes it possible to have access, almost instantaneously, to the estimation of various key magnitudes of the operating conditions of the coupling. The lineic contact forces obtained from an approach mainly analytical are in very good agreement with one of the most refined numerical approaches for spline coupling without and with angular misalignment. From these lineic forces, an estimation of the contact peak pressures is done and can be used in wear model or fretting models. In the same way, one simply calculates from these lineic force distributions the reaction moments and the reaction forces.
The stiffness per unit length of the assembly as well as the tooth stiffness is also in agreement with the results of strictly numerical approaches. By integrating the stiffness along the spline joint, a global stiffness for torsional vibration analysis is available.
For a given level of misalignment, the calculation of the bearing reactions (forces and moments) is easily accessible and thus makes it possible to estimate what the bearing dimensions should be. By knowing the stiffness of these bearings, the theoretical misalignment value can be corrected.
- F. Curà, A. Mura, Evaluation of the fretting wear damage on crowned splined couplings, in 2nd International Conference on Structural Integrity, ICSI 2017, 4–7 September 2017, Funchal, Madeira, Portugal [Google Scholar]
- S. Monti, M. Guagliano, Failure analysis of an involute spline coupling of an overhead bridge crane, Eng. Fail. Anal. 104, 321–330 (2019) [Google Scholar]
- A. Mura, F. Curà, A. Di Molfetta, Investigation of bearings overloads due to misaligned splined shafts, Proc. Struct. Integr. 12, 52–57 (2018) [CrossRef] [Google Scholar]
- J. Wang, J. Zhang, Z. Yao, R. Sun, Y. Zhao, Nonlinear characteristics of a multi-degree-of-freedom spur gear system with bending-torsional coupling vibration, Mech. Syst. Sig. Process. 121, 810–827 (2019) [CrossRef] [Google Scholar]
- J. Hong, D. Talbot, A. Kahraman, Load distribution analysis of clearance-fit spline joints using finite elements, Mech. Mach. Theory 74, 42–57 (2014) [Google Scholar]
- R.A. Adey, J. Baynham, J.W. Taylor, Development of analysis tools for spline couplings, Proc. Inst. Mech. Eng. G 214, 347–357 (2000) [CrossRef] [Google Scholar]
- M.A. Alfares, A.H. Falah, A.H. Elkholy, Clearance distribution of misaligned gear coupling teeth considering crowning and geometry variations, Mech. Mach. Theory 41, 1258–1272 (2006) [Google Scholar]
- F. Curà, A. Mura, Experimental and theoretical investigation about reaction moment in misaligned spline coupling, Mech. Syst. Signal Process. 45, 504–512 (2014) [Google Scholar]
- N. Fillot, Etude mécanique de l'usure, Thesis INSA Lyon, 2004 [Google Scholar]
- J. Ding, W.S. Sum et al., Fretting fatigue predictions in a complex coupling. Int. J. Fatigue 29, 1229–1244 (2007) [Google Scholar]
- W. Qureshi, F. Curà, A. Mura, Characterization of fretting wear experiments on spline couplings by principal component analysis, Inst. Mech. Eng. Part J: J. Eng. Tribol. 231, 860–868 (2017) [CrossRef] [Google Scholar]
- F. Curà, W. Qureshi, A. Mura, A methodological approach for incremental fretting wear formulation, Tribol. Lett. 64, 20 (2016) [Google Scholar]
- V. Cuffaro, F. Curà, A. Mura, Surface Characterization of spline coupling teeth subjected to fretting wear, Proc. Eng. 74, 135–142 (2014) [CrossRef] [Google Scholar]
- J. Hong, D. Talbot, A. Kahraman, A semi-analytical load distribution model for side-fit involute splines, Mech. Mach. Theory 76, 39–55 (2014) [Google Scholar]
- S. Medina, A.V. Olver, An analysis of misaligned spline couplings, Proc. Inst. Mech. Eng. J 216, 269–278 (2002) [Google Scholar]
- S. Medina, A.V. Olver, Regimes of contact in spline coupling, J. Tribol. 124, 351–357 (2000) [Google Scholar]
- V. Hariharan, P. Srinivasan, Vibration analysis of parallel misaligned shaft with ball bearings system, S.J. Sci. Technol. 33, 61–68 (2011) [Google Scholar]
- M. Byrtus, V. Zeman, On modeling and vibration of gear drives influenced by nonlinear couplings, Mech. Mach. Theory 46, 375–397 (2011) [Google Scholar]
- A.H. Elkholy, M.A. Alfares, Misalignment loads in splined gear couplings, Int. J. Comput. Appl. Technol. 15, 128–137 (2002) [CrossRef] [Google Scholar]
- J.R. Mancuso, Moments and forces imposed on power transmission systems due to misalignment of a crowned tooth coupling, Master's thesis, Pennsylvania State University, 1971 [Google Scholar]
- Y. Guo, A. Lambert, R. Wallen, R. Errichello, J. Keller, Theoretical and experimental study on gear-coupling contact and loads considering misalignment, torque and friction influences, Mech. Mach. Theory 98, 242–262 (2016) [Google Scholar]
- P. Barrot, Analyse et modélisation du comportement des liaisons cannelées, Thesis INSA Toulouse, Novembre 2006 [Google Scholar]
- A. Barrot, M. Paredes, M. Sartor, Extended equations of load distribution in the axial direction in a spline coupling, Eng. Fail. Anal. 16, 200–211 (2009) [Google Scholar]
- A. Tjernberg, Load distribution in the axial direction in a spline coupling, Eng. Fail. Anal. 8, 357–570 (2001) [Google Scholar]
- J.-L. Ligier, M. Benoit, N. Pakpoor Gilani, Rapport interne COMATEC HEIG-VD, Mars 2018 [Google Scholar]
- C. Dupertuis, Désalignement des liaisons cannelées et effets induits, Travail de Bachelor, HEIG-VD, septembre 2019 [Google Scholar]
- A. Sackfield, C.E. Truman, D.A. Hills, The tilted punch under normal and shear load (with application to fretting tests), Int. J. Mech. Sci. 43, 1881–1892 (2001) [CrossRef] [Google Scholar]
- R. Fargere, Simulation du comportement dynamique des transmissions par engrenages sur paliers hydrodynamiques, Thesis INSA Lyon, 2012 [Google Scholar]
Cite this article as: C. Dupertuis, J.-L. Ligier, Contact pressure in misaligned spline couplings, Mechanics & Industry 21, 510 (2020)
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.