# An Analytical Model of Multiarc Sprag Clutch considering Geometry and Internal Interaction during Engagement

Authors: Chuang Huang; Ming Liu; Yongqiang Zhao
Journal: Shock and Vibration
 (2017)

---

## Abstract

A new multiarc sprag clutch model considering geometry and internal interaction during engagement is proposed in this paper. To increase the accuracy of the model, an improved model of geometric deformation coordination is presented to describe the basic geometrical quantities after rotation of the sprag. Then, based on this model, a novel nonlinear iteration method focusing on the varied contact radius is proposed to compute the normal contact force and show a good agreement with the FEM model. In addition, the alternate friction model considering stationary and rate-dependency friction is formulated and applied in the contacts of both the inner race and outer race. The proposed model is verified by comparison with other published results and experimental results. In subsequent analysis, the new model fully reveals the dynamical behaviors of the multiarc clutch during preload engagement under torque excitation; therefore, it will be useful for the performance analysis and dynamic design of multiarc sprag clutch.


---

## Body

## 1. Introduction

The circular-arc sprag clutch is the type most commonly used today in applications such as the helicopter transmission and passing and cutting torque from the power source. To provide the greatest possible life and torque capacity, the sprag has attained a high degree of refinement with precision-machined compound curves such as multiarc, which provide optimum contact angles for driving and permitting higher transmitted torques for allowable Hertz and hoop stresses. Therefore, the performance of the sprag clutch is directly affected by geometric parameters, especially the contact compound curve of the sprag. Lynwander et al. [1] presented the design procedure for multiarc sprag clutch from the static view, including the deformations of the components and the geometric relationship. Kish [2] further proposed the iterative design method for single-arc sprag clutch containing additional centrifugal deformation of the outer race. However, the dynamic behavior of the sprag clutch cannot be explained from the perspective of statics, and failure still occurs during operation for the excessive dynamic pressure on the axial ends of the sprags or the wear caused by stick-slip during engagement. To address these problematic dynamic behaviors, it is essential to thoroughly understand the mechanism of the sprag clutch and to recognize the factors that affect its performance.Figure1 shows a manufacturers cut-out view of a complete sprag clutch, while Figure 2 gives a simplified schematic of its basic components. Such a sprag clutch is composed of five components: accurately formed sprags, a cylindrical inner race, a cylindrical outer race, energizing springs, and a sprag cage. Each sprag is essentially a strut placed between the races in such a way that it transmits power from one race to the other by an engaging action when either race is rotated in the driving direction. Rotation in the other direction disengages the sprags and the clutch is free or overruns. Either race may be the driven member or driving member. Note that the energizing spring plays a key role in the transition between the engagement and overrunning, as the contact normal force is very small at this time. In most circumstance, the specific compound curves such as multiarc are applied on the inner contact profile of the sprag, while the single circular of the outer contact profile remains unchanged.Figure 1
A cross-sectional view of Formsprag indexing clutch model.Figure 2
A schematic diagram of a multiarc sprag clutch and engaging conditions.Because the sprag clutch is a type of friction component with interactions in components due to the engaging action, its geometric deformation coordination, contact nonlinearity, and friction are the main attributes responsible for its dynamic performance during engagement. These attributes have been studied extensively over the last twenty years.Initially, total torsional stiffness and damping models are used to describe the dynamic characteristics of the single-arc sprag clutch, and then stepwise refined models have been proposed for single-arc sprag clutch concerning the geometry and internal interactions. The interactions include the normal contact forces and the tangential frictional forces between the sprags and races. Chassapis and Lowen [3] presented the torsional stiffness-damping model of a single-arc sprag clutch to simulate the dynamic behavior of a press-feed mechanism, in which energy dissipation was described by a linear damping model. Based on the work, Xu and Lowen [4] put forward the improved total torsional stiffness-damping model with a nonlinear torsional damping, which gave better agreement with experiment results. Further, Xu and Lowen [5] proposed a complete mathematical model, using an epicyclical train model to describe the geometry deformation coordination and motions of components and deriving the calculation method of the normal contact force, whereas it assumed pure rolling in contacts without slipping. Based on the same epicyclical train model, Liu et al. [6] presented a mathematical model of logarithmic spiral type sprag clutch for design and analysis, taking advantage of Hunt-Crossley contact impact theory to calculate the normal force but also neglected the slip in contacts. To investigate the dynamic behavior of stick-slip, Vernay et al. [7, 8] adapted a Coulomb friction model to describe the tangential interactions between the sprags and the inner race, and the structural deformations of the races were firstly considered in the geometry deformation coordination model; however, no slip was assumed in the contacts between the sprags and the outer race, and the edge contact was neglected. In addition, researchers constructed the model of single-arc sprag clutch on the software platform such as ADAMS [9, 10].From the above discussions, to establish accurate dynamic modelling of the sprag clutch, three key issues should be discussed, including geometry deformation coordination model, the calculation method of the normal contact force, and the friction model in contacts. Compared to the simplified epicyclic train model only including local contact deformations, the geometry deformation coordination model developed by Vernay further took structural deformations into account. However, the centroid position of the sprag was assumed and the edge contact was not considered, which may lead to the axial end crushing on the sprag due to excessive contact stress. To the best of the author’s knowledge, public literature mainly predicted about the dynamic performance of the single-arc sprag clutch. No studies have constructed a multiarc sprag model concerning the additional contact nonlinearity of the multiarc compound curve. Additionally, existing work on the friction model adopted the Coulomb friction model in the inner race contacts, which could reflect the behavior of stick-slip but did not consider the rate-dependency friction in boundary lubrication and the friction in the contacts of the outer race.The purpose of this paper is to develop a multiarc sprag clutch model considering geometry and internal interaction, fully reflecting the dynamic behavior of multiarc sprag clutch during engagement. First, to improve the accuracy of the model, an improved geometric deformation coordination model is presented, considering the actual centroid position of the sprag, the non-Hertzian local contact deformations by the slice method, and structural deformation of the races. Second, to consider the contact nonlinearity of the multiarc compound curve, a novel nonlinear iteration method focusing on the varied contact radius is proposed to compute the normal contact force. Third, the alternate friction considering stationary and rate-dependency friction is formulated and applied in the contacts of both the inner race and the outer race, increasing a degree of revolution for the whole model of the sprag clutch.This paper will be organized as follows. In the following section, the calculation method of the normal contact force and friction force in contacts are proposed, and a multiarc sprag clutch model considering geometry and internal interaction is established. In the validation section, the mathematical equations of the model are solved and the theoretical responses are compared to published results and experiment results. In the results and discussion section, the dynamic performance of multiarc sprag clutch is discussed during the preload engaging with torque excitation. In the last section, some conclusions are given.


## 2. The Problem Formulation

Figure2 shows a simplified configuration of a complete sprag clutch in initial contact position and three sketch maps to illustrate the changes of engagement conditions. In initial contact position, uniform sprags are evenly distributed in annular space between the races and held in contact with the races by the energizing springs. The inner contact surface of the sprag is composed of three tangent circles with their own radius and center positions, and the switching of the circles in contact depends on the wedging position. The associated contact angles ω and υ are expressed with respect to a line connecting the inner and outer contact points A and B, and β and γ are described by the lines OCi and OCg with respect to the line OCe, respectively. Particularly, unlike the previous research, the centroid position of the sprag Cg has no assumption on it and only depends on the mass distribution of the sprag. Re is the radius of inner surface of the outer race; Ri is the radius of outer surface of the inner race; re is the radius of outer surface of the sprag, and Ce is the corresponding center position; the ri1, ri2, and ri3 are the different radius of inner surface of the sprag composed of three tangent circles, and Ci is a general term for the center positions C1, C2, and C3 changing with the wedging position.As shown in Figure2, there are three types of engagement conditions in the sprag clutch: Engaging into, preload engaging, and Engaging out. “Engaging into” sprags roll into a wedging position between the races when the driving torque on the outer race is sufficient to overcome the resisting torque on the inner race; otherwise, the “Engaging into” sprags become the “Engaging out” sprags until the sprags roll out of engagement. When “Engaging into” sprags stop the rotation with the driving torque, “Engaging into” sprags convert to “pre-load” engagement sprags, in which the races have no relative motion, and the driving torque is taken as the preload applied to the outer race.A mode describing the real engagement process of the sprag clutch would be extremely complex. In order to obtain an effective and practical mathematical model, the following assumptions are made.(1)
All sprags are wedged simultaneously.(2)
Sprags are circumferentially uniform and share the load evenly.(3)
The inner race and the outer race are concentric definitely, without considering the manufacturing precision and the installation error of the sprag.


### 2.1. The Geometry Deformation Coordination

#### 2.1.1. The Deformations in Components

In actual engagement of the sprag clutch, the sprags and the races deform with the rotation of the sprags, and substantial errors would occur if they are still considered as rigid components. To consider the deformations of components, the geometry deformation coordination model has been included in previous study, but the centroid position of the sprag was assumed to simplify the model and the local deformation did not consider edge contact. To improve the accuracy of the model, an improved model of geometry deformation coordination is proposed in this paper. The actual centroid position is included, and the edge contact is in consideration by the slice method, which has been effectively adapted in roller bearing [11].Figure3 shows the deformations of components after rotation of the sprag. Note that the sprags rotate and force the races to deform radially and compress the sprags from the orientation of contact points A and B, when the sprag clutch is in “Engaging into” condition. As a result of the wedging action, two types of deformations are generated: the local contact deformation and the structural deformation. The local deformations of the sprags interacting the outer and inner race are defined as δHE and δHI, respectively. The structural deformations for the inner and outer race are expressed as δBI and δBE, and the structural deformation for sprags is written as δg. Particularly, the additional centrifugal deformation of the outer race is defined asδce. The following discusses the calculation method used to obtain the deformations and basic geometric quantities of the sprag clutch after rotation.Figure 3
The deformations of components after rotation of the sprag.The local contact deformationsδHE and δHI are calculated by the slice method [12] shown in Figure 4, in which the sprag is axially divided into a certain number of slices with equal length, and the interaction among the deflections of slices is established by the influence coefficient method. The functional relationship between load and deflection on each slice can be obtained as follows:(1)Swi·qi=Δi,Swe·qe=Δe.Figure 4
The slice method applied on contact between the sprag and the races.In each of the above equations,Sw is the matrix of the weighted influence coefficients with specific expression:(2)Sw=n∑j,kwj,ks·wj,k⋯s·wj,n⋮⋱⋮s·wn,k⋯s·wn,nforj=1⋯n,k=1⋯n,where s are the elastic compliances associated with each slice in contact with either the inner or the outer race and wj,k is weighting function, meaning the influence of a unit pressure on j slice on the deflection of the k slice. Moreover, because the weighting function has the relationship wj,k=wk,j, therefore, Sw is a symmetric matrix.q is the vector of unknown forces on each slice element of the sprag.(3)q=qj⋮qnforj=1⋯n.Δ is the vector of deflections belonging to each slice element of the sprag, and the deflection is used to represent the interpenetration between rigid contact partners.(4)Δ=δj1/ex⋮δn1/exforj=1⋯n.The structural deformations for the races are calculated based on the Lame formula of thick-walled cylinder theory; the equations for structural deformation of the inner raceδBI and that of the outer race δBE can be obtained as follows:(5)δBE=Nben2πlERed2+Re2Red2-Re2+ν,δBI=Nbin2πlERi2+Rid2Ri2-Rid2-ν,where l is the length of the sprag, n is the number of the sprag, Red is the radius of the outer surface of the outer race, Rid is the radius of the inner surface of the inner race, E refer to the elastic modulus, ν represents Poisson’s ratio, and Nbi and Nbe represent the normal contact force on the inner race or outer race, respectively.With the structure deformation of the races, high radial loads are applied on both ends of the sprag and cause the deformation, which cannot be ignored. The compressive deflection for the sprag itself can be written as(6)δg=NnormalRe-RiblE,where Nnormal is the normal contact force on the inner race or outer race. Either one for the calculation is reasonable, since their values are very close due to the small contact angles ω and υ.Moreover, the calculation of the additional centrifugal deformation for the outer race is based on the hollow disk theory. For high speed application, the outer race can be seen as a rotating disc with a central aperture due to its large radial dimension. Therefore, the additional deformation of the outer race can be written as(7)δce=ReE3+υρ4πn302Red2+Re21-υ3+υ.


#### 2.1.2. Basic Geometric Quantities after Deformation

Based on the above deformations, the essential geometric quantities of length in the sprag clutch have corresponding changes and can be expressed as(8)OA=Ri-δBI-δHI2,OB=Re+δBE+δHE2+δce,OCi=Ri+ri-δBI-δg2-δHI,OCe=Re-re+δBE+δg2+δHE.As shown in Figure3, the angles ∠OCGCi, β, and γ vary with the changes of OA, OB, OCi, and OCe, and the equations of the angles can be expressed as(9)∠OCGCi=θi-α=acos⁡a2+OCg2-OCi22×CiCg×OCG,β=acos⁡OCi2+OCe2-CiCe22×OCi×OCe,γ=acos⁡OCg2+OCe2-CgCe22×OCg×OCe,where θi is the initial value of the angle ∠OCGCi in initial contact position and α is the rotation angle of the sprag.


#### 2.1.3. Calculation of the Normal Contact Force

Based on the above equations of geometric deformation coordination, the normal contact force on the inner race would be obtained in each position and can be expressed as(10)Nbi=fα,Ri,Re,ri,re,OCi,δHE,δHI,δBE,δBI,δg,δce,α>0,0,else,where ri is a general term for ri1, ri2, and ri3, and OCi is a general term for OC1, OC2, and OC3. Their values depend on the position of the contact point A and change to the corresponding value when sprags roll across the junction point of tangent circles. Focusing on the varied contact radius on inner surface of the sprag, a novel nonlinear iteration method should be proposed to compute the normal contact force during whole engagement.Figure5 shows the mechanics of the novel nonlinear iteration method. First, a judgment is always carried out within each time step to determine whether the contact point has crossed the near junction point of tangent circles. If not, the calculation of normal contact force continues to proceed. Second, the normal contact force Nbi_junction is calculated when sprags rotate to the junction point, and the new rotation angle αnew is obtained using the new radius, corresponding center, and the normal contact force Nbi_junction. Third, the additional normal contact force Nbi_left is given corresponding to the left rotation angle α-αnew in the time step. At last, the total normal contact force is gained in the time step, passing through the junction point.Figure 5
The flow chart of the proposed method of normal contact force.In order to verify the effectiveness of the proposed method, the torsional stiffness of multiarc sprag clutch is calculated to compare with the other two methods, using the measured parameters of sprag clutch [13]. As shown in the appendix, the calculated torsional stiffness by the slice method has a better agreement with the finite element method (FEM) results than the model based on Hertzian theory, reflecting the discontinuity of its change rate at the junction points of circles.


### 2.2. Calculation of Friction Force

According to Vernay et al.’s research [7, 8], the motions between the sprags and races are pure rolling with no relative motions in most work time, and the friction state of interaction in sprag clutch belongs to the static friction phase with no lubrication. But when the required friction force keeping the pure rolling motion exceeds the maximum static friction, slip will occur in contacts and the relative motions exist between the sprags and races; because of the slip, lubrication should be considered in this phase, and the changeable sliding rate should take into account the friction force model to reflect the Stribeck effect. Therefore, as shown in Figure 6, the friction model should include static friction, rate-dependent friction, and the transformation between them [14, 15]. To overcome the discontinuity of the transition point, a narrow band of zero velocity [16] is introduced. Taking the friction in the inner race for example, the expression of the alternate friction force can be written as(11)Tbix,ri=Tbifx,ri,ifVelbi<DV,Tbifx,ri<TbisstickTbissgnTbifx,ri,ifVelbi<DV,Tbifx,ri>TbissticktosliptransitionNbiμbdVelbisgnVelbi,ifVelbi>DVslipFigure 6
Alternate friction model considering static and rate-dependent friction.with(12)μbdVelbi=μcb+μsb-μcbe-γbVelbi,where Velbi is the sliding velocity at the contact point of inner race and DV is a narrow band of zero velocity. x is the system state vector, Tbif(x,ri) is the static force to keep the pure rolling contact, Tbis is the maximum static friction force, and μbd is the velocity-depending friction coefficient under boundary lubrication, in which the exponential decaying behavior is consistent with the experimental results [17].Similar friction model is used to calculate the friction force on the outer raceTbe, and specific processes are followed to calculate the friction forces on the inner and outer races. First, the friction forces Tbi and Tbe are calculated to keep pure rolling at contact points, and Tbi is always larger than Tbe during engagement. Then, Tbif(x,ri) is compared to the maximum static friction force Tbis in each time step. It must have Tbif(x,ri)<Tbis for the engagement without slip; otherwise, slip begins. The direction of the slip depends on the relative velocity at contact point. Additionally, the similar process is carried out to calculate the friction force on the outer race Tbef(x,re).


### 2.3. Differential Equations of Motion

Based on Newton’s second law, the differential equations of the motion can be derived and describe the dynamic behavior of multiarc sprag clutch during engagement. To derive the equations of motion, first of all, it is necessary to analyze the forces acting on the sprag, the inner race, and the outer race.The forces acting on the sprag are shown in Figure7. Nbi and Nbe are the contact normal forces acting between the sprags and the races, Tbi and Tbe are the friction forces between the sprags and the races, Crl is the viscous damping between the races, Cbi and Cbe are the constant torques applied on the inner and outer races, respectively, Kr is the stiffness of the elastic blades holding the sprags and the races in contact, αmax is the assumed rotation angle of the sprag to ensure the initial contact torque formed by the energizing springs, and the Kr and αmax are calculated from the initial contact equilibrium equations.Figure 7
A schematic diagram of forces on the sprag.The application of Newton’s second law yields, for the sprags,(13)Igα¨+φ¨-beTbe-Krαmax-α-biTbi+hiNbi-heNbe=0,-mgOCgφ¨-cos⁡β+γTbi+sin⁡β+γNbi+Tbecos⁡γ-Nbesin⁡γ=0,with(14)bi=OCg·cos⁡β+γ-OA,be=OB-OCg∗cos⁡γ,hi=OCg·sin⁡β+γ,he=OCg·sin⁡γ,where Ig is the inertia of the sprag, Ibi and Ibe are the inertias of the inner race and outer races, respectively, φ and α are the two degrees of the sprag, φ is the azimuthal angle of the azimuth coordinate system oxGyG→ with respect to the absolute coordinate system oxoyo→, and α is the rotational angle of the body fitted coordinate system oxgyg→ with respect to the azimuth coordinate system oxGyG→.For the inner and outer race,(15)Ibiθ¨bi-nsTbiRi-δBI-δHI2+Crlθ˙bi-θ˙be-Cbi=0,Ibeθ¨be+Crlθ˙be-θ˙bi+nsTbeRe+δBE+δHE2+δce-Cbe=0,where θbi, θbe are the degrees of the inner and outer race, respectively.The system equations of motion are nonlinear ordinary second-order differential equations. They can be solved using a fourth-order Runge-Kutta algorithm with fixed time step. The time step for the numerical procedure is chosen to beΔt=1×10-5s.


## 2.1. The Geometry Deformation Coordination

### 2.1.1. The Deformations in Components

In actual engagement of the sprag clutch, the sprags and the races deform with the rotation of the sprags, and substantial errors would occur if they are still considered as rigid components. To consider the deformations of components, the geometry deformation coordination model has been included in previous study, but the centroid position of the sprag was assumed to simplify the model and the local deformation did not consider edge contact. To improve the accuracy of the model, an improved model of geometry deformation coordination is proposed in this paper. The actual centroid position is included, and the edge contact is in consideration by the slice method, which has been effectively adapted in roller bearing [11].Figure3 shows the deformations of components after rotation of the sprag. Note that the sprags rotate and force the races to deform radially and compress the sprags from the orientation of contact points A and B, when the sprag clutch is in “Engaging into” condition. As a result of the wedging action, two types of deformations are generated: the local contact deformation and the structural deformation. The local deformations of the sprags interacting the outer and inner race are defined as δHE and δHI, respectively. The structural deformations for the inner and outer race are expressed as δBI and δBE, and the structural deformation for sprags is written as δg. Particularly, the additional centrifugal deformation of the outer race is defined asδce. The following discusses the calculation method used to obtain the deformations and basic geometric quantities of the sprag clutch after rotation.Figure 3
The deformations of components after rotation of the sprag.The local contact deformationsδHE and δHI are calculated by the slice method [12] shown in Figure 4, in which the sprag is axially divided into a certain number of slices with equal length, and the interaction among the deflections of slices is established by the influence coefficient method. The functional relationship between load and deflection on each slice can be obtained as follows:(1)Swi·qi=Δi,Swe·qe=Δe.Figure 4
The slice method applied on contact between the sprag and the races.In each of the above equations,Sw is the matrix of the weighted influence coefficients with specific expression:(2)Sw=n∑j,kwj,ks·wj,k⋯s·wj,n⋮⋱⋮s·wn,k⋯s·wn,nforj=1⋯n,k=1⋯n,where s are the elastic compliances associated with each slice in contact with either the inner or the outer race and wj,k is weighting function, meaning the influence of a unit pressure on j slice on the deflection of the k slice. Moreover, because the weighting function has the relationship wj,k=wk,j, therefore, Sw is a symmetric matrix.q is the vector of unknown forces on each slice element of the sprag.(3)q=qj⋮qnforj=1⋯n.Δ is the vector of deflections belonging to each slice element of the sprag, and the deflection is used to represent the interpenetration between rigid contact partners.(4)Δ=δj1/ex⋮δn1/exforj=1⋯n.The structural deformations for the races are calculated based on the Lame formula of thick-walled cylinder theory; the equations for structural deformation of the inner raceδBI and that of the outer race δBE can be obtained as follows:(5)δBE=Nben2πlERed2+Re2Red2-Re2+ν,δBI=Nbin2πlERi2+Rid2Ri2-Rid2-ν,where l is the length of the sprag, n is the number of the sprag, Red is the radius of the outer surface of the outer race, Rid is the radius of the inner surface of the inner race, E refer to the elastic modulus, ν represents Poisson’s ratio, and Nbi and Nbe represent the normal contact force on the inner race or outer race, respectively.With the structure deformation of the races, high radial loads are applied on both ends of the sprag and cause the deformation, which cannot be ignored. The compressive deflection for the sprag itself can be written as(6)δg=NnormalRe-RiblE,where Nnormal is the normal contact force on the inner race or outer race. Either one for the calculation is reasonable, since their values are very close due to the small contact angles ω and υ.Moreover, the calculation of the additional centrifugal deformation for the outer race is based on the hollow disk theory. For high speed application, the outer race can be seen as a rotating disc with a central aperture due to its large radial dimension. Therefore, the additional deformation of the outer race can be written as(7)δce=ReE3+υρ4πn302Red2+Re21-υ3+υ.


### 2.1.2. Basic Geometric Quantities after Deformation

Based on the above deformations, the essential geometric quantities of length in the sprag clutch have corresponding changes and can be expressed as(8)OA=Ri-δBI-δHI2,OB=Re+δBE+δHE2+δce,OCi=Ri+ri-δBI-δg2-δHI,OCe=Re-re+δBE+δg2+δHE.As shown in Figure3, the angles ∠OCGCi, β, and γ vary with the changes of OA, OB, OCi, and OCe, and the equations of the angles can be expressed as(9)∠OCGCi=θi-α=acos⁡a2+OCg2-OCi22×CiCg×OCG,β=acos⁡OCi2+OCe2-CiCe22×OCi×OCe,γ=acos⁡OCg2+OCe2-CgCe22×OCg×OCe,where θi is the initial value of the angle ∠OCGCi in initial contact position and α is the rotation angle of the sprag.


### 2.1.3. Calculation of the Normal Contact Force

Based on the above equations of geometric deformation coordination, the normal contact force on the inner race would be obtained in each position and can be expressed as(10)Nbi=fα,Ri,Re,ri,re,OCi,δHE,δHI,δBE,δBI,δg,δce,α>0,0,else,where ri is a general term for ri1, ri2, and ri3, and OCi is a general term for OC1, OC2, and OC3. Their values depend on the position of the contact point A and change to the corresponding value when sprags roll across the junction point of tangent circles. Focusing on the varied contact radius on inner surface of the sprag, a novel nonlinear iteration method should be proposed to compute the normal contact force during whole engagement.Figure5 shows the mechanics of the novel nonlinear iteration method. First, a judgment is always carried out within each time step to determine whether the contact point has crossed the near junction point of tangent circles. If not, the calculation of normal contact force continues to proceed. Second, the normal contact force Nbi_junction is calculated when sprags rotate to the junction point, and the new rotation angle αnew is obtained using the new radius, corresponding center, and the normal contact force Nbi_junction. Third, the additional normal contact force Nbi_left is given corresponding to the left rotation angle α-αnew in the time step. At last, the total normal contact force is gained in the time step, passing through the junction point.Figure 5
The flow chart of the proposed method of normal contact force.In order to verify the effectiveness of the proposed method, the torsional stiffness of multiarc sprag clutch is calculated to compare with the other two methods, using the measured parameters of sprag clutch [13]. As shown in the appendix, the calculated torsional stiffness by the slice method has a better agreement with the finite element method (FEM) results than the model based on Hertzian theory, reflecting the discontinuity of its change rate at the junction points of circles.


## 2.1.1. The Deformations in Components

In actual engagement of the sprag clutch, the sprags and the races deform with the rotation of the sprags, and substantial errors would occur if they are still considered as rigid components. To consider the deformations of components, the geometry deformation coordination model has been included in previous study, but the centroid position of the sprag was assumed to simplify the model and the local deformation did not consider edge contact. To improve the accuracy of the model, an improved model of geometry deformation coordination is proposed in this paper. The actual centroid position is included, and the edge contact is in consideration by the slice method, which has been effectively adapted in roller bearing [11].Figure3 shows the deformations of components after rotation of the sprag. Note that the sprags rotate and force the races to deform radially and compress the sprags from the orientation of contact points A and B, when the sprag clutch is in “Engaging into” condition. As a result of the wedging action, two types of deformations are generated: the local contact deformation and the structural deformation. The local deformations of the sprags interacting the outer and inner race are defined as δHE and δHI, respectively. The structural deformations for the inner and outer race are expressed as δBI and δBE, and the structural deformation for sprags is written as δg. Particularly, the additional centrifugal deformation of the outer race is defined asδce. The following discusses the calculation method used to obtain the deformations and basic geometric quantities of the sprag clutch after rotation.Figure 3
The deformations of components after rotation of the sprag.The local contact deformationsδHE and δHI are calculated by the slice method [12] shown in Figure 4, in which the sprag is axially divided into a certain number of slices with equal length, and the interaction among the deflections of slices is established by the influence coefficient method. The functional relationship between load and deflection on each slice can be obtained as follows:(1)Swi·qi=Δi,Swe·qe=Δe.Figure 4
The slice method applied on contact between the sprag and the races.In each of the above equations,Sw is the matrix of the weighted influence coefficients with specific expression:(2)Sw=n∑j,kwj,ks·wj,k⋯s·wj,n⋮⋱⋮s·wn,k⋯s·wn,nforj=1⋯n,k=1⋯n,where s are the elastic compliances associated with each slice in contact with either the inner or the outer race and wj,k is weighting function, meaning the influence of a unit pressure on j slice on the deflection of the k slice. Moreover, because the weighting function has the relationship wj,k=wk,j, therefore, Sw is a symmetric matrix.q is the vector of unknown forces on each slice element of the sprag.(3)q=qj⋮qnforj=1⋯n.Δ is the vector of deflections belonging to each slice element of the sprag, and the deflection is used to represent the interpenetration between rigid contact partners.(4)Δ=δj1/ex⋮δn1/exforj=1⋯n.The structural deformations for the races are calculated based on the Lame formula of thick-walled cylinder theory; the equations for structural deformation of the inner raceδBI and that of the outer race δBE can be obtained as follows:(5)δBE=Nben2πlERed2+Re2Red2-Re2+ν,δBI=Nbin2πlERi2+Rid2Ri2-Rid2-ν,where l is the length of the sprag, n is the number of the sprag, Red is the radius of the outer surface of the outer race, Rid is the radius of the inner surface of the inner race, E refer to the elastic modulus, ν represents Poisson’s ratio, and Nbi and Nbe represent the normal contact force on the inner race or outer race, respectively.With the structure deformation of the races, high radial loads are applied on both ends of the sprag and cause the deformation, which cannot be ignored. The compressive deflection for the sprag itself can be written as(6)δg=NnormalRe-RiblE,where Nnormal is the normal contact force on the inner race or outer race. Either one for the calculation is reasonable, since their values are very close due to the small contact angles ω and υ.Moreover, the calculation of the additional centrifugal deformation for the outer race is based on the hollow disk theory. For high speed application, the outer race can be seen as a rotating disc with a central aperture due to its large radial dimension. Therefore, the additional deformation of the outer race can be written as(7)δce=ReE3+υρ4πn302Red2+Re21-υ3+υ.


## 2.1.2. Basic Geometric Quantities after Deformation

Based on the above deformations, the essential geometric quantities of length in the sprag clutch have corresponding changes and can be expressed as(8)OA=Ri-δBI-δHI2,OB=Re+δBE+δHE2+δce,OCi=Ri+ri-δBI-δg2-δHI,OCe=Re-re+δBE+δg2+δHE.As shown in Figure3, the angles ∠OCGCi, β, and γ vary with the changes of OA, OB, OCi, and OCe, and the equations of the angles can be expressed as(9)∠OCGCi=θi-α=acos⁡a2+OCg2-OCi22×CiCg×OCG,β=acos⁡OCi2+OCe2-CiCe22×OCi×OCe,γ=acos⁡OCg2+OCe2-CgCe22×OCg×OCe,where θi is the initial value of the angle ∠OCGCi in initial contact position and α is the rotation angle of the sprag.


## 2.1.3. Calculation of the Normal Contact Force

Based on the above equations of geometric deformation coordination, the normal contact force on the inner race would be obtained in each position and can be expressed as(10)Nbi=fα,Ri,Re,ri,re,OCi,δHE,δHI,δBE,δBI,δg,δce,α>0,0,else,where ri is a general term for ri1, ri2, and ri3, and OCi is a general term for OC1, OC2, and OC3. Their values depend on the position of the contact point A and change to the corresponding value when sprags roll across the junction point of tangent circles. Focusing on the varied contact radius on inner surface of the sprag, a novel nonlinear iteration method should be proposed to compute the normal contact force during whole engagement.Figure5 shows the mechanics of the novel nonlinear iteration method. First, a judgment is always carried out within each time step to determine whether the contact point has crossed the near junction point of tangent circles. If not, the calculation of normal contact force continues to proceed. Second, the normal contact force Nbi_junction is calculated when sprags rotate to the junction point, and the new rotation angle αnew is obtained using the new radius, corresponding center, and the normal contact force Nbi_junction. Third, the additional normal contact force Nbi_left is given corresponding to the left rotation angle α-αnew in the time step. At last, the total normal contact force is gained in the time step, passing through the junction point.Figure 5
The flow chart of the proposed method of normal contact force.In order to verify the effectiveness of the proposed method, the torsional stiffness of multiarc sprag clutch is calculated to compare with the other two methods, using the measured parameters of sprag clutch [13]. As shown in the appendix, the calculated torsional stiffness by the slice method has a better agreement with the finite element method (FEM) results than the model based on Hertzian theory, reflecting the discontinuity of its change rate at the junction points of circles.


## 2.2. Calculation of Friction Force

According to Vernay et al.’s research [7, 8], the motions between the sprags and races are pure rolling with no relative motions in most work time, and the friction state of interaction in sprag clutch belongs to the static friction phase with no lubrication. But when the required friction force keeping the pure rolling motion exceeds the maximum static friction, slip will occur in contacts and the relative motions exist between the sprags and races; because of the slip, lubrication should be considered in this phase, and the changeable sliding rate should take into account the friction force model to reflect the Stribeck effect. Therefore, as shown in Figure 6, the friction model should include static friction, rate-dependent friction, and the transformation between them [14, 15]. To overcome the discontinuity of the transition point, a narrow band of zero velocity [16] is introduced. Taking the friction in the inner race for example, the expression of the alternate friction force can be written as(11)Tbix,ri=Tbifx,ri,ifVelbi<DV,Tbifx,ri<TbisstickTbissgnTbifx,ri,ifVelbi<DV,Tbifx,ri>TbissticktosliptransitionNbiμbdVelbisgnVelbi,ifVelbi>DVslipFigure 6
Alternate friction model considering static and rate-dependent friction.with(12)μbdVelbi=μcb+μsb-μcbe-γbVelbi,where Velbi is the sliding velocity at the contact point of inner race and DV is a narrow band of zero velocity. x is the system state vector, Tbif(x,ri) is the static force to keep the pure rolling contact, Tbis is the maximum static friction force, and μbd is the velocity-depending friction coefficient under boundary lubrication, in which the exponential decaying behavior is consistent with the experimental results [17].Similar friction model is used to calculate the friction force on the outer raceTbe, and specific processes are followed to calculate the friction forces on the inner and outer races. First, the friction forces Tbi and Tbe are calculated to keep pure rolling at contact points, and Tbi is always larger than Tbe during engagement. Then, Tbif(x,ri) is compared to the maximum static friction force Tbis in each time step. It must have Tbif(x,ri)<Tbis for the engagement without slip; otherwise, slip begins. The direction of the slip depends on the relative velocity at contact point. Additionally, the similar process is carried out to calculate the friction force on the outer race Tbef(x,re).


## 2.3. Differential Equations of Motion

Based on Newton’s second law, the differential equations of the motion can be derived and describe the dynamic behavior of multiarc sprag clutch during engagement. To derive the equations of motion, first of all, it is necessary to analyze the forces acting on the sprag, the inner race, and the outer race.The forces acting on the sprag are shown in Figure7. Nbi and Nbe are the contact normal forces acting between the sprags and the races, Tbi and Tbe are the friction forces between the sprags and the races, Crl is the viscous damping between the races, Cbi and Cbe are the constant torques applied on the inner and outer races, respectively, Kr is the stiffness of the elastic blades holding the sprags and the races in contact, αmax is the assumed rotation angle of the sprag to ensure the initial contact torque formed by the energizing springs, and the Kr and αmax are calculated from the initial contact equilibrium equations.Figure 7
A schematic diagram of forces on the sprag.The application of Newton’s second law yields, for the sprags,(13)Igα¨+φ¨-beTbe-Krαmax-α-biTbi+hiNbi-heNbe=0,-mgOCgφ¨-cos⁡β+γTbi+sin⁡β+γNbi+Tbecos⁡γ-Nbesin⁡γ=0,with(14)bi=OCg·cos⁡β+γ-OA,be=OB-OCg∗cos⁡γ,hi=OCg·sin⁡β+γ,he=OCg·sin⁡γ,where Ig is the inertia of the sprag, Ibi and Ibe are the inertias of the inner race and outer races, respectively, φ and α are the two degrees of the sprag, φ is the azimuthal angle of the azimuth coordinate system oxGyG→ with respect to the absolute coordinate system oxoyo→, and α is the rotational angle of the body fitted coordinate system oxgyg→ with respect to the azimuth coordinate system oxGyG→.For the inner and outer race,(15)Ibiθ¨bi-nsTbiRi-δBI-δHI2+Crlθ˙bi-θ˙be-Cbi=0,Ibeθ¨be+Crlθ˙be-θ˙bi+nsTbeRe+δBE+δHE2+δce-Cbe=0,where θbi, θbe are the degrees of the inner and outer race, respectively.The system equations of motion are nonlinear ordinary second-order differential equations. They can be solved using a fourth-order Runge-Kutta algorithm with fixed time step. The time step for the numerical procedure is chosen to beΔt=1×10-5s.


## 3. Validation

To verify the effectiveness of the proposed model, the parameters and operating condition of Vernay’s model are adapted, and a comparative analysis of the sprag clutch is carried out. By setting the same radius and center of the three circles, the multiarc sprag clutch became the single-arc sprag clutch shown in Table1. The operating condition is that firstly the inner race accelerates and moves to the free direction and then reverses direction and engages suddenly when the inner race meets the elastic obstacle on the torsion bench [8].Table 1
The parameters of the sprag clutch applied in the simulation.
Parameter
Value
nbg
12
R
i
d (mm)
7.21
R
i (mm)
11.11
R
e (mm)
19.44
R
e
d (mm)
26.23
r
i
1 (mm)
4.02
r
i
2 (mm)
4.02
r
i
3 (mm)
4.02
m
g (kg)
1.8×10-3
I
1 (kg⋅m2)
1.97
×
1
0
-
4
I
b
i (kg⋅m2)
1.41
×
1
0
-
3
I
b
e (kg⋅m2)
1.41×10-3
I
g (kg⋅m2)
1.18
×
1
0
-
8
K
r (Nm/rad)
0.01
C
r
l (Nms/rad)
0.1
C
f (Nm)
1.66
×
1
0
-
2The comparison of the inner race speed process against the Vernay’s model and the experimental data are shown in Figure8. In general, the shape and amplitude of the curve obtained from the proposed model approximate the experiment results better than the curves from Vernay’s model. This good agreement shows the effectiveness of the new model. Note that the damping in the simulation is considered constant equal to an average value obtained from Vernay’s experiments, because the nonlinearity of stiffness is more important than the nonlinear damping on the dynamic behavior of the freewheel. This is the reason why the new model and the Vernay model agree with each other well and still have similar error profiles with respect to the experimental results. In addition, the friction forces are calculated using the alternate friction model to obtain the torque transmitted, focusing on the problem that whether slip exists in contacts during engagement. However, slip in initial engagement lasts only half thousandth of a second when the sprag clutch initially engages and has little effect on the torque transmitted [7, 8].Figure 8
Comparison of results of the new model with results of Vernay’s model and experiment.


## 4. Results and Discussion

Based on the proposed model of multiarc sprag clutch, the dynamic behavior focusing on edge contact and stick-slip are investigated during preload engaging. To investigate the dynamic behavior, it is necessary to define the structural parameters and operating conditions for the sprag clutch. First, based on the parameters of single-circle sprag clutch listed in Table1, A multiarc sprag clutch is formed by changing the suitable value of ri2 and ri3 to 4.08 and 4.15, respectively. Then, the centroid position coordinates are automatically obtained through the CAD model of the clutch using the macro programming method. Subsequently, the operating conditions are defined as follows: with the preload Ce1 on the inner race, a sine torque excitation of Ce2=Amp×sin⁡(Ω×t) is applied to the outer race, and a high equivalent stiffness spring is assumed to connect the outer race to a fixed abutment [7]. Two different combinations for Ce1 and Ce2 are selected for preload engaging: one is the high torsional vibration excitation Ce2 with respect to preload Ce1, in which the amplitude of Ce2 is assigned to 6 Nm, 8 Nm, and 10 Nm, respectively, and Ce1 is set to 3.34 Nm; the other is low torsional vibration excitation Ce2 with respect to preload Ce1, in which the amplitude of Ce2 is assigned to 6 Nm, 8 Nm, and 10 Nm respectively, and Ce1 is set to 16.8 Nm.


### 4.1. Preload Engaging Suffering Torsional Vibration Excitation

#### 4.1.1. High Torsional Vibration Excitation with respect to Preload

Figure9 shows abnormal disengagement phenomenon of multiarc sprag clutch, including the time varying of the rotation angle of the inner race, the contact stress distribution under different excitation torques, and the friction coefficient on the inner contact points.Figure 9
Abnormal disengagement phenomenon of multiarc sprag clutch.
(a)
Slip to the disengagement direction
(b)
Excitation torque 6 Nm
(c)
Variation of friction coefficient
(d)
Excitation torque 10 NmThe rotation angle of the inner race has a clear offset to the overrunning direction whenCe2 is 8 Nm and 10 Nm, but no offset when Ce2 is 6 Nm. With the offset, the lower limit of the contact stress moves toward to zero, and the friction coefficient changes its value, transforming between the static and sliding friction coefficient. The reason for the offset is that slip appears in contacts, which moves the angle of the inner race to the overrunning direction until out of the engagement. As shown in Figure 9(d), the contact stress is close to zero when the slip appears, indicating the sprags are about to depart from the races. At this time, the spring force exerted by the spring blade plays a key role to end the slip, holding the sprag to reengaging into the races, since the contact forces are small at this moment.Note that the negative value of the friction coefficient means that the direction of the friction torque for the sprags is opposite to the engagement direction at the moment. As shown in Figure9(c), the slip appears as the value of friction coefficient is more than the maximum static friction and the friction force changes its direction when the values of Ce2 are 8 Nm and 10 Nm. There is no slip when the value Ce2 is 6 Nm, because the corresponding friction coefficient is no more than the maximum static friction.


#### 4.1.2. Low Torsional Vibration Excitation with respect to Preload

According to Vernay’s study, the maximum friction coefficient for steel-steel contact is dynamic and varies in a range of 0.04 to 0.12 during mixed lubrication. For this operating condition, slip does not occur when the maximum friction coefficient is near its upper limit but appears when the maximum friction coefficient is toward to its lower limit. As shown in Figure10, slip occurs until the dynamic maximum friction coefficient is reduced to 0.082. In addition, the loads on the races also influence the slip in contacts. When excitation torque Ce2 is 8 Nm and 10 Nm, the curves reveal the offset of the inner race moving to the direction of engagement, rather than the direction of overrunning, but no offset when the value Ce2 is 6 Nm.Figure 10
Slip phenomenon of multiarc sprag clutch under low torsional vibration excitation with respect to preload.
(a)
Slip of angle of inner race to engagement direction
(b)
Excitation torque 6 Nm
(c)
Variation of friction coefficient
(d)
Excitation torque 10 NmNote that the dynamic contact stress fluctuates around the reference static contact stress determined by the preload, while it would be of the same order of magnitude with allowable contact stress under rated preload torque. If slip occurs with high contact stress, serious wear can be caused on surface of the sprag. Therefore, the friction coefficient cannot be lower than a certain value in this operating condition, which puts forward requirements for surface material, surface technology, and lubrication design for the sprag clutch.


## 4.1. Preload Engaging Suffering Torsional Vibration Excitation

### 4.1.1. High Torsional Vibration Excitation with respect to Preload

Figure9 shows abnormal disengagement phenomenon of multiarc sprag clutch, including the time varying of the rotation angle of the inner race, the contact stress distribution under different excitation torques, and the friction coefficient on the inner contact points.Figure 9
Abnormal disengagement phenomenon of multiarc sprag clutch.
(a)
Slip to the disengagement direction
(b)
Excitation torque 6 Nm
(c)
Variation of friction coefficient
(d)
Excitation torque 10 NmThe rotation angle of the inner race has a clear offset to the overrunning direction whenCe2 is 8 Nm and 10 Nm, but no offset when Ce2 is 6 Nm. With the offset, the lower limit of the contact stress moves toward to zero, and the friction coefficient changes its value, transforming between the static and sliding friction coefficient. The reason for the offset is that slip appears in contacts, which moves the angle of the inner race to the overrunning direction until out of the engagement. As shown in Figure 9(d), the contact stress is close to zero when the slip appears, indicating the sprags are about to depart from the races. At this time, the spring force exerted by the spring blade plays a key role to end the slip, holding the sprag to reengaging into the races, since the contact forces are small at this moment.Note that the negative value of the friction coefficient means that the direction of the friction torque for the sprags is opposite to the engagement direction at the moment. As shown in Figure9(c), the slip appears as the value of friction coefficient is more than the maximum static friction and the friction force changes its direction when the values of Ce2 are 8 Nm and 10 Nm. There is no slip when the value Ce2 is 6 Nm, because the corresponding friction coefficient is no more than the maximum static friction.


### 4.1.2. Low Torsional Vibration Excitation with respect to Preload

According to Vernay’s study, the maximum friction coefficient for steel-steel contact is dynamic and varies in a range of 0.04 to 0.12 during mixed lubrication. For this operating condition, slip does not occur when the maximum friction coefficient is near its upper limit but appears when the maximum friction coefficient is toward to its lower limit. As shown in Figure10, slip occurs until the dynamic maximum friction coefficient is reduced to 0.082. In addition, the loads on the races also influence the slip in contacts. When excitation torque Ce2 is 8 Nm and 10 Nm, the curves reveal the offset of the inner race moving to the direction of engagement, rather than the direction of overrunning, but no offset when the value Ce2 is 6 Nm.Figure 10
Slip phenomenon of multiarc sprag clutch under low torsional vibration excitation with respect to preload.
(a)
Slip of angle of inner race to engagement direction
(b)
Excitation torque 6 Nm
(c)
Variation of friction coefficient
(d)
Excitation torque 10 NmNote that the dynamic contact stress fluctuates around the reference static contact stress determined by the preload, while it would be of the same order of magnitude with allowable contact stress under rated preload torque. If slip occurs with high contact stress, serious wear can be caused on surface of the sprag. Therefore, the friction coefficient cannot be lower than a certain value in this operating condition, which puts forward requirements for surface material, surface technology, and lubrication design for the sprag clutch.


## 4.1.1. High Torsional Vibration Excitation with respect to Preload

Figure9 shows abnormal disengagement phenomenon of multiarc sprag clutch, including the time varying of the rotation angle of the inner race, the contact stress distribution under different excitation torques, and the friction coefficient on the inner contact points.Figure 9
Abnormal disengagement phenomenon of multiarc sprag clutch.
(a)
Slip to the disengagement direction
(b)
Excitation torque 6 Nm
(c)
Variation of friction coefficient
(d)
Excitation torque 10 NmThe rotation angle of the inner race has a clear offset to the overrunning direction whenCe2 is 8 Nm and 10 Nm, but no offset when Ce2 is 6 Nm. With the offset, the lower limit of the contact stress moves toward to zero, and the friction coefficient changes its value, transforming between the static and sliding friction coefficient. The reason for the offset is that slip appears in contacts, which moves the angle of the inner race to the overrunning direction until out of the engagement. As shown in Figure 9(d), the contact stress is close to zero when the slip appears, indicating the sprags are about to depart from the races. At this time, the spring force exerted by the spring blade plays a key role to end the slip, holding the sprag to reengaging into the races, since the contact forces are small at this moment.Note that the negative value of the friction coefficient means that the direction of the friction torque for the sprags is opposite to the engagement direction at the moment. As shown in Figure9(c), the slip appears as the value of friction coefficient is more than the maximum static friction and the friction force changes its direction when the values of Ce2 are 8 Nm and 10 Nm. There is no slip when the value Ce2 is 6 Nm, because the corresponding friction coefficient is no more than the maximum static friction.


## 4.1.2. Low Torsional Vibration Excitation with respect to Preload

According to Vernay’s study, the maximum friction coefficient for steel-steel contact is dynamic and varies in a range of 0.04 to 0.12 during mixed lubrication. For this operating condition, slip does not occur when the maximum friction coefficient is near its upper limit but appears when the maximum friction coefficient is toward to its lower limit. As shown in Figure10, slip occurs until the dynamic maximum friction coefficient is reduced to 0.082. In addition, the loads on the races also influence the slip in contacts. When excitation torque Ce2 is 8 Nm and 10 Nm, the curves reveal the offset of the inner race moving to the direction of engagement, rather than the direction of overrunning, but no offset when the value Ce2 is 6 Nm.Figure 10
Slip phenomenon of multiarc sprag clutch under low torsional vibration excitation with respect to preload.
(a)
Slip of angle of inner race to engagement direction
(b)
Excitation torque 6 Nm
(c)
Variation of friction coefficient
(d)
Excitation torque 10 NmNote that the dynamic contact stress fluctuates around the reference static contact stress determined by the preload, while it would be of the same order of magnitude with allowable contact stress under rated preload torque. If slip occurs with high contact stress, serious wear can be caused on surface of the sprag. Therefore, the friction coefficient cannot be lower than a certain value in this operating condition, which puts forward requirements for surface material, surface technology, and lubrication design for the sprag clutch.


## 5. Conclusion

To understand the dynamic behavior of multiarc sprag clutch, a new multiarc sprag clutch is proposed to consider the geometry and interaction in this paper. The following conclusions can be drawn from the work presented in this paper.(1)
To improve the accuracy of the model, an improved model of geometric deformation coordination considering centroid position and edge contact have been presented. Based on the improved model, the calculated torsional stiffness has better agreement with FEM results than the previous model, and the dynamic behavior of the sprag clutch reveals the contact stress concentration throughout the engagement.(2)
To reflect the contact nonlinearity of the multiarc compound curve, a novel nonlinear iteration method of normal contact force is proposed. Using the method, the calculated torsional stiffness of multiarc sprag clutch has good agreement with FEM results, and the change rate of it presents the discontinuity at the junction points of circles.(3)
The alternate friction model considering stationary and rate-dependency friction is formulated and applied in the contacts of both the inner race and outer race, increasing a degree of revolution for the whole model of the sprag clutch. Due to the slip in contacts, the sprags may disengage with the races in the condition of high torsional vibration excitation with respect to preload, and serious wear would be caused in the condition of low torsional vibration excitation with respect to preload, when the maximum friction coefficient is toward its lower limit.(4)
The proposed model of multiarc sprag clutch is effective to reflect the dynamic behavior during engagement, especially the edge contact and stick-slip. Using this model, further research can be performed to determine the key influencing factors and help to modify the design and control parameters of multiarc sprag clutch.



---

*Source: 2393578-2017-03-21.xml*