Next Article in Journal
In Vitro Toxicity of Bone Graft Materials to Human Mineralizing Cells
Next Article in Special Issue
Comparative Manufacturing of Hybrid Composites with Waste Graphite Fillers for UAVs
Previous Article in Journal
Effects of Organic Vehicle on the Rheological and Screen-Printing Characteristics of Silver Paste for LTCC Thick Film Electrodes
Previous Article in Special Issue
Analysis of Experimental Results Regarding the Selection of Spring Elements in the Front Suspension of a Four-Axle Truck
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Method for the Calculation of Characteristics of Disc Springs with Trapezoidal Cross-Sections and Rounded Edges

by
Dominik Sebastian Leininger
*,
Max Benedikt Geilen
,
Marcus Klein
and
Matthias Oechsner
Center for Engineering Materials (MPA-IfW), Technical University of Darmstadt, Grafenstraße 2, 64283 Darmstadt, Germany
*
Author to whom correspondence should be addressed.
Materials 2022, 15(5), 1954; https://doi.org/10.3390/ma15051954
Submission received: 1 February 2022 / Revised: 3 March 2022 / Accepted: 4 March 2022 / Published: 6 March 2022

Abstract

:
In the European standards specifying disc spring manufacturing, geometry, shape and characteristic, an edge rounding is prescribed. Common methods for the calculation of disc spring characteristics, even in these standards, are based on a rectangular cross-section. This discrepancy can lead to a considerable divergence of the computed characteristic from the characteristic determined by testing. In literature, this divergence has not yet been examined with regard to rounded edges. In this paper, a new method addressing this problem is introduced. For this purpose, the geometry of idealized disc springs is parameterized. Based on four edge radii and two angles of the inner and outer faces, equations to compute the initial cone angle and the lever arm are introduced. These equations are used to formulate an algorithm to adapt other computation methods to non-rectangular cross-sections and rounded edges. The method is applied to the formulas by AlmenLaszlo, CurtiOrlando, Zheng and those by Kobelev. FE simulations of disc springs with rounded edges and a non-rectangular cross-section were used to verify the new formulas. The results show that the introduced method can be applied to known characteristic computation methods and result in a model expansion taking cross-section variations into account. The adjusted characteristics show more accurate alignment to the FE simulation for the cross-section variations investigated. These findings not only close the geometric gap between the manufacturing guidelines and the computation on an analytical basis, they also define a new parameter space for designs of disc springs and a corresponding force computation method to optimize spring characteristics.

1. Introduction

Disc springs are axially loaded conical shaped washers with a wide range of applications [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34]. They are specified in the European standards DIN EN 16983 [35] and DIN EN 16984 [36]. The theoretical geometry of a spring is defined with a rectangular cross-section and sharp edges, whose characteristics can be calculated with formulas from Section 5 in DIN EN 16984 [36]. However, the manufacturing process defined by Table 7 in DIN EN 16983 [35] includes edge rounding.
As prescribed in Section 1 in DIN EN 16983 [35], proper function of a disc spring includes compliance with the spring force tolerances, among other requirements. Despite the development of different analytical approaches to compute force—deflection relations, the simplifications made in the geometry can cause inaccuracies in the calculation of disc spring characteristics. Since the manufacturing process can also lead to a non-rectangular, i.e., trapezoidally-skewed, cross-section, in addition to the radii, this non-conformance is also to be considered. Therefore, in this paper the idealized geometry of such disc springs is described and a new computational approach is introduced to take the geometric variations into account. The described method can be combined with different existing computation methods for disc springs with rectangular cross-section and sharp edges.
A finite element (FE) simulation was used in order to be able to evaluate the consequences of a wide range of different geometries. For setup and preprocessing of the FE models, the Spring_stack-module [37,38] was used.
The models used employ first order, reduced integration, axisymmertic finite elements sized approximately one percent of the spring thickness. Contact plates are implemented as analytical rigids. They are connected to the springs through a frictionsless, hard, surface to surface, finite sliding contact. The characteristics are computed by continuously increasing the displacement of the contact place close to edge I.
Literature on the characteristics of disc springs dates back to 1936 when Almen and Laszlo published their well-known formula [39]. The standard noted above also uses a modified version of this computational method. Almen and Laszlo used the tangential stresses to link the load of a disc spring with its deflection. Curti and Orlando published a similar approach incorporating radial stresses [40]. An energy-based approach to compute the characteristics was published by Zheng et al. [41]. Koblev [42] presents an analytical approach capable of approximating high working angle disc springs. Luxemburg and Givli [43] model non-rigid rotation.
Different approaches consider friction through an additional frictional moment [44], through a master curve concept [45] or through a coupling of axial and radial forces [42]. Some approaches allow the combination of different frictional coefficients [46,47,48] and rate-dependent frictional coefficients [47,48]. The latter [47,48] also considers contact stiffness.
Several models have been generated to model deviating geometries, including slotted disc springs [44,49,50,51,52], variable thickness disc springs [53,54,55] and disc springs with contact flats [56,57].
As an alternative to analytical approaches, disc springs have been modelled using FE simulations [2,3,4,5,8,10,13,30,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83]. FE simulations are computationally more expensive than analytical approaches. Furthermore, they often depend on software licenses. The main advantages, however, are the easy variation in geometry and the simulation of different material properties.

2. Idealized Geometry of Disc Springs

In the standards DIN EN 16983 [35] and DIN EN 16984 [36], the disc spring and its cross-section are characterized by the four edges I through IV as shown in Figure 1. The loading force F acts on edges I and III .
The dimensions of the spring are measured with the inner and outer diameter, D i and D e , the total height l 0 and the thickness t [35,36] as shown in Figure 2. The lever arm V measures the radial distance of the force application edges. Note that DIN EN 16983 [35] and DIN EN 16984 [36] the depiction of the lever arm V is different than in Figure 1. It is shown as the radial distance between the points I and IV , not between I and III , Figure 1a in [35,36]. This may be a consequence of the last revision correcting the depiction of D e , which shares the same extension line with the dimension of V. The vertical space below edge II is characterized by the cone height h 0 . It is approximated by
h 0 , DIN = l 0 t
using the small angle approximation [35,36]. The cone height is used to compute the deflection to flat position, s f , in which the upper and lower faces are horizontal in an idealized model.
Figure 3 shows an idealized cross-section of an exemplary real disc spring. The virtual position of the sharp edges are named I through IV , as, due to edge rounding, they are not material edges. The lever arm V in Figure 3 differs from the one in Figure 2. Furthermore, the deflection s f to flat position can not be directly measured as the cone height h 0 between II and the lower force application point, which moves on the circular arc of edge III during loading. The deflection s f to flat position is defined as
s f = l 0 t
which is generally not identical to the cone height h 0 , but matches the approximated cone height h 0 , DIN for rectangular cross-sections with sharp edges.
The non-rectangularity of the cross-section is idealized by a trapezoid with rounded corners described by the parameters shown in Figure 4. The four radii r I through r IV characterize the rounding of the corresponding edges with the center points I * through IV * . The trapezoidality caused by the inner and outer face of the disc spring is measured counterclockwise on the right cross-section of the spring by the angles β i and β e . Thus, positive values of these angles cause the inner and outer faces to align towards the axial direction of the disc spring. The upper and lower faces are assumed to be parallel.
Since the load is applied axially, the lever arm V (equal to the radial distance of the force application lines) can be measured as the radial distance of the points I * and III * . The initial slope angle of the disc spring is named φ .
Due to the fact that in the following investigations only geometric quantities of the cross-section are considered and thus the distance between the cross-section and the center-line of the disc spring can be neglected, the complexity is reduced by introducing the difference in radius R:
R = D e D i 2
The dimensions R and l 0 are measured as shown in Figure 4, so the points I through IV lie outside of the rectangle described by these quantities. Thus, these values are measurable for real disc springs using a caliper. While this is only true for β i , β e > φ , such large angles of the inner and outer faces are fairly uncommon in practice. In Figure 5 a cross-section with β i , β e > φ is shown. Note that the difference in radius R is measured from the arcs defining the edge rounding of the same edges as in Figure 4, II and IV , even though the points to measure are inside the cross-section and do not mark the radial extent of the cross-section in this case.

3. Calculation of Initial Slope Angle φ and Lever Arm V

As shown in Figure 6, the bounding box of the rounded off cross-section with dimensions R and l 0 can be enlarged to fit the cross-section without edge rounding. This outer bounding box touches the points I through IV and its offset is described by the four parameters a I through a IV which are added to l 0 and R, respectively.
The width of the lower face c l as well as the height of the outer and inner face, t ˜ i and t ˜ e , can be calculated to:
c l = II III ¯
t ˜ i = I II ¯ = t cos ( β i )
t ˜ e = III IV ¯ = t cos ( β e )
Thus, t ˜ e and t ˜ i are always greater t as long as the angles β i and β e are not equal to zero and their absolute values are in an ordinary range, e.g., | β i | , | β e | < 45 ° .
Therefore, the radial width of the outer bounding box can be calculated using the blue and the orange triangle in Figure 6:
a II + R + a IV = cos ( φ ) c l + sin ( φ β e ) t ˜ e
Analogously, the blue and green triangles describe the axial height of the cross-section’s outer bounding box:
a I + l 0 + a III = sin ( φ ) c l + cos ( φ β i ) t ˜ i
Eliminating c l , from Equations (7) and (8) follows:
a II + R + a IV sin ( φ β e ) t ˜ e cos ( φ ) = a I + l 0 + a III cos ( φ β i ) t ˜ i sin ( φ )
Thus, the initial slope angle φ can be calculated as the fixed point of the following equation:
φ = arctan a I + l 0 + a III cos ( φ β i ) t ˜ i a II + R + a IV sin ( φ β e ) t ˜ e
For known initial slope angle φ , the lever arm V can be calculated as follows. Similar to Equation (7), another formula including V can be used to calculate the radial distance of the points II and IV by using the orange and green triangles in Figure 6:
a II + R + a IV = sin ( φ β i ) t ˜ i + b I + V + b III + sin ( φ β e ) t ˜ e
V = R + a II + a IV b I b III sin ( φ β i ) t ˜ i sin ( φ β e ) t ˜ e
To determine the unknown lengths a I through a IV as well as the radial distance b I of the points I * , I and the radial distance b III of the points III * , III , the rounded edges are examined. Exemplarily, edge I is discussed below.
As shown in Figure 7, the distance a I measures the protrusion of I out of the inner bounding box which touches the rounded cross-section. It follows from the geometry of the blue triangle in Figure 7 that:
a I = sin ( φ ) r ˜ I
where r ˜ I can be calculated using the orange triangle in Figure 7.
r ˜ I = r ˜ I r ˜ I
r ˜ I = tan φ 2 r I
As depicted in Figure 7, r ˜ I is computed as the tangent of the bisecting angle of the sum of the right angle of a rectangular cross-section (dashed in Figure 7) combined with the deviation from that, β i , visualized by the green-stroked triangle in Figure 7:
r ˜ I = tan π 2 + β i 2 r I
Inserting r ˜ I from Equation (14) into Equation (13) by using Equations (16) and (15), yields
a I = sin ( φ ) tan π 2 + β i 2 tan φ 2 r I
To calculate b I (necessary for V), the already known r ˜ I from Equation (16) is used. As Figure 8 shows, b I can be described as the difference between b I and b I , which in turn can be computed using the blue and orange triangles in Figure 8:
b I = b I b I
b I = cos ( φ ) r ˜ I
b I = sin ( φ ) r I
This leads to the following equation to calculate b I :
b I = cos ( φ ) tan π 2 + β i 2 sin ( φ ) r I
= cos ( φ ) tan π 2 + β i 2 tan ( φ ) r I
The derivations of the parameters of the other edges, a II through a IV and b II through b IV , are outlined in Appendix A.

4. Calculation of Deformed Geometries of Loaded Disc Springs

In the previous sections, the disc spring was always assumed to experience neither deflection nor deformation. Under load, these assumptions obviously can’t hold true. In this section, a rotated cross-section of a deflected disc spring will be analyzed to overcome these restrictions. The cross-section is considered to be undistortable but able to rotate. This simplification is considered to be a good approximation and is used in several other computation approaches [39,40].
Since the beginning of Section 3, the cross-section is considered to be independent of the centerline of the disc spring. Since the loading does not affect this assumption, all measures are relative to the cross-section and the cross-section is assumed to be undistortable, the rotation center can be chosen freely, as the translation with respect to another rotation center, e.g., the neutral point in [39], is irrelevant. Here, the center of rotation is chosen as fixed at the center of edge rounding III * . This leads to a fixed position of III * during deflection and only the other point in question, the center of edge rounding I * , experiences displacement under load.
In Figure 9, the rotated cross-section is shown in blue after a rotation by ψ . This rotation results in a deflection s and an slightly lengthened lever arm V ψ . The total height l 0 of an unloaded (black) disc spring can be split up by introducing the axial distance Λ between the points I * and III * , as shown in Figure 9:
Λ = l 0 r I r III ,
The radial distance between I * and III * is measured with the already known lever arm V.
The angle θ , describing the angle between the horizontal and the connection line of the two points I * and III * , as well as their distance X can be computed as:
θ = arctan Λ V
X = Λ 2 + V 2
With these quantities, it follows from the definition of sine and cosine that:
V = cos ( θ ) X
Λ = sin ( θ ) X
They can be generalized for the loaded disc spring by subtracting the rotation angle ψ from the initial value θ :
V ψ = cos ( θ ψ ) X
= sin ( ψ ) Λ + cos ( ψ ) V
Λ ψ = sin ( θ ψ ) X
Furthermore, the following equation can be obtained from Figure 9 to describe the deflection s of the disc spring.
s = Λ Λ ψ
By using Equations (27) and (30), the following equation gives the axial deflection s as a function of the rotation angle ψ :
s = sin θ sin ( θ ψ ) X
The inverse of this,
ψ = θ arcsin sin ( θ ) s X ,
is used in Equation (28) to calculate the deflection-dependent lever arm V ψ as a function of the deflection s:
V ψ = cos arcsin sin ( θ ) s X X
= 1 sin ( θ ) s X 2 X
By setting ψ = φ the disc spring is in flat position and the deflection holds the condition s = s f . With Equation (32), s f can be calculated as follows:
s f = sin ( θ ) sin ( θ φ ) X
It can be shown that this is equivalent to Equation (2) for rectangular cross-sections with sharp edges.
It can be shown that s f is generally unequal to the cone height h 0 of point II above the lower force application point of the spring. This is due to the non-rectangularity and the rounded edges of the cross-section. The different definition of h 0 in the standard, which uses the small angle approximation and define the cone height h 0 , DIN as in Equation (1), Table 1 in [36], may cause an additional divergence of the two measures for non-rectangular cross-section and rounded edges.
The reduced height l ψ of a loaded disc spring can be computed analogous to Equation (23) as
l ψ = r I + Λ ψ + r III .
For a disc spring compressed to flat position ( s = s f , ψ = φ ), l ψ is equal to t.
Most analytical methods approximate the lever arm as the deflection-independent difference in radius R. It can be shown that for cross-sections satisfying condition (38), this leads to an slight overestimation of the lever arm and thus to a slight underestimation of the force.
R V ψ > 0 R 1 cos ( ψ ) > l 0 r I r III sin ( ψ ) + a II + a IV b I b III sin ( φ β i ) × t ˜ i sin ( φ β e ) × t ˜ e cos ( ψ )

5. Calculation of Spring Characteristics

In this section the basic approach to adjust the calculation of characteristics to non-rectangular cross-sections with rounded edges is presented. Subsequently, this is used to adapt the algorithm of AlmenLaszlo to those geometries as an example.

5.1. Concept

In the original formulas, the necessary force to create a specific deflection s is computed based on a set of basic geometry parameters of a rectangular cross-section with sharp edges described in Section 2. The algorithm adjustment described below is based on the use of a more precise lever arm and two geometric modifications to increase the accuracy of the original formulas when applied to non-rectangular cross-sections and rounded edges: the first shape modification adjusts the diameters and the total height to values as they would be of a rectangular cross-section with sharp edges. This leads to an increase of the area of the cross-section, which must be compensated by the second modification. These two modifications result in the description of a cross-section compatible with the original formulas with properties very similar to those of a spring with a non-rectangular cross-section and rounded edges under investigation.
As dicussed in Section 3 and Section 4, due to the rounded edges and the non-rectangularity of the cross-section, the lever arm of such springs can differ considerably from the one with a rectangular cross-section and sharp edges. This is accounted for by using the deflection-dependent lever arm V ψ described in Equation (34) instead of using e.g., the deflected independent difference in radius R, as Almen and Laszlo do to approximate the lever arm.
As the original computation methods are designed for rectangular cross-sections with sharp edges, the shape modifications noted above are used to imitate this kind of cross-section, while retaining the basic characteristics such as the initial slope angle. Figure 10 shows a comparison of the different cross-sections.
It can be observed that the diameters D e and D i as well as the total height l 0 need to be slightly modified to the adjusted diameters D e , D i and the adjusted total height l 0 to maintain the basic properties of the disc spring such as the initial slope angle. The outer diameter D e is increased by a II on both sides of the disc spring to imitate sharp edges. Furthermore, a negative angle β e of the outer face of the disc spring (as shown in Figure 10) causes the outer diameter to be larger, which is compensated on both sides of the disc spring (factor 2) by the last term in the equation below:
D e = D e + 2 a II 2 × cos ( φ ) tan ( β e ) t 2 = D e + 2 a II cos ( φ ) tan ( β e ) t
Analogously, the following equation is obtained for the adjusted inner diameter D i :
D i = D i 2 a IV + cos ( φ ) tan ( β i ) t
The total height l 0 must be adjusted by a I and a III to compensate the edge rounding and the influence of the angles β i and β e .
l 0 = l 0 + a I + a III sin ( φ ) tan ( β i ) + tan ( β e ) t 2
Thus, the above adjustments create a cross-section with the edge radii, outer face angle and inner face angle of zero, but retaining the slope angle and the length of the cross-section’s center-line. With this, the cross-section area is slightly larger which results in an increased resistance to deformation. This is compensated by the second adjustment of the diameter measures which result in D e and D i :
D e = D e 2 k e
D i = D i + 2 k i
As shown in Figure 10, the product of the cosine of the slope angles and the correction lengths k i and k e would be the accurate adjustment. Due to the fact that Almen and Laszlo [39] use the small angle approximation and the integration over the area of the cross-section uses the diameters as limits, this is neglected in the second shape adjustment. Since the limits of the second integral uses the thickness t, the second adjustment only affects the diameters and it has no influence on the adjusted total height ( l 0 = l 0 ).
When applying this approach to other characteristic calculation methods, the choice of whether to apply these simplifications or not, has to be made individually for those other methods, based on the use of the small angle approximation and the limits in the calculation of resistance to deformation. In three of the four methods used in Section 6, this simplification was applicable due to small angle approximation in the definition of the integration limits in the original formulas. In the adjustment of the Kobelev-formula, the cosine noted above must be considered (see Appendix B).
The correction lengths k i and k e used in the second shape adjustment are defined by comparing the outer and inner part of a rectangular cross-section with a cross-section with rounded edges. To simplify this, two averaged edge radii, r e and r i , are introduced, which are defined so that the area of the cross-section remains constant. Equation (44) shows the balance of area differences to a cross-section with sharp edges to define r i :
1 π 4 2 r i 2 = 1 π 4 r I 2 + r II 2
This can be solved for r i . The average outer edge radius r e can be computed analogously:
r i = r I 2 + r II 2 2
r e = r III 2 + r IV 2 2
With those average edge radii, r e and r i , the correction lengths k i and k e for the inner and outer diameter are defined such that the cross-sectional area remains unchanged.
t ( r i k i ) = 2 π 4 r i 2 + ( t 2 r i ) r i
This can be solved for the inner correction length k i . In this consideration the influence of the angles β i and β e is neglected. The outer correction length k e can be computed analogously:
k i = 1 2 t ( 4 π ) r i 2
k e = 1 2 t ( 4 π ) r e 2
The correction lengths k i and k e from Equations (48) and (49) indicate by which quantity the diameters need to be adjusted to keep the area of the cross-section the same, see Equations (42) and (43).

5.2. Applying the Method to Almen’s Formula

To apply the method described above and to compute the characteristic of a disc spring with non-rectangular cross-section and rounded edges, the adjusting method is applied to the formulas published in 1937 by Almen and Laszlo [39].
Almen and Laszlo based their calculation on the equality of an internal moment M i and an external moment M e [39]:
M e = M i
The internal moment M i results from the stress caused by the rotation of the cross-section [39]. Its computation is not modified here except for the use of the adjusted diameters D e , D i and the adjusted total height l 0 derived above.
The external moment on a sector d ϑ is calculated as the product of the part of the force F A , o applied on this sector and the lever arm of the force application [39]. Almen and Laszlo simplify the lever arm by measuring the distance between radii of the inner and outer edge of the neutral surface. According to [84], the use of the measured diameters instead of those of the neutral surface in the formulas results in a better match with characteristics determined experimentally. This measurement of the diameters is also done in the standard DIN EN 16984 [36] and by Zheng et al. [41], and thus is used in the current work for AlmenLaszlo-, CurtiOrlando and Zheng-based computation. It is not used for the Kobelev-based computation [85] as the results couldn’t be improved using different diameter definitions.
On a rectangular cross-section, the external moment is computed by
M e = F A , o d ϑ 2 π D e 2 sin ( φ ) t 2 D i 2 + sin ( φ ) t 2 = F A , o d ϑ 2 π R sin ( φ ) t
using the difference in radius defined in Equation (3).
Almen and Laszlo use this formula with the small angle approximation in a simplified form [39]:
M e = F A , o d ϑ 2 π R
In order to address the non-rectangular cross-section and the rounded edges, the lever arm is not simplified by the difference in radius but described by the deflection angle-dependent lever arm V ψ defined in Equation (34). The external moment is calculated with Equation (53), so the force F A differs from the one calculated with the original Almen formula, F A , o :
M e = F A d ϑ 2 π V ψ
Taking the above replacement of the diameters in the underlying Almen-method and the more accurate lever arm V ψ into account, the force F A at a given axial deflection s is computed according to the adjusted Almen formula, Equation (7) in [39] by:
F A = 4 E s t 1 μ 2 D e 2 M h 0 , DIN s h 0 , DIN s 2 + t 2
1 M = δ + 1 δ 1 2 ln ( δ ) π δ δ 1 2
where the adjusted diameter ratio δ and the adjusted cone height using the small angle approximation h 0 , DIN is defined analogously to [39] p. 311, and Equation (1), respectively:
δ = D e D i
h 0 , DIN = l 0 t
The Equations (52) and (53) lead to the following formula to replace the lever arm V ψ in the Almen-equation to compute the force F A , o :
F A = F A , o R V ψ
Again, the above replacement of the diameters in the underlying Almen-method can be used in Equation (58), so the computation of the force F A can simply be written as:
F A = F A , o ( D e , D i , l 0 ) R V ψ
with the adjusted difference in radius defined analogously to Equation (3):
R = D e D i 2
Equation (54) as well as Equation (59) derived above show an exemplary application of the new method introduced in Section 5.1 to adopt other algorithms to a non-rectangular cross-section with rounded edges. To show the effect of the adjustments, the method is also applied to the algorithm by Curti and Orlando [40], as well as the one by Zhang [41], by replacing the force F A , o computed using the original Almen formula with the force F C , o or F Z , o computed using Curti’s or Zheng’s original formulas, respectively. The algorithm by Kobelev [85] is also adapted by using the force F K , o instead of F A , o in the above equations and replacing the adjusted difference in radius R by the lever arm H r , K used by Kobelev. The final formulas are written out in full in Appendix B.

6. Results and Discussion

The characteristics computed according to the equations derived above are compared to those of AlmenLaszlo, CurtiOrlando, Zheng and Kobelev as well as results of FE simulations by applying them to multiple cross-sectional geometries. For all springs, Young’s modulus E is defined as 206 G Pa and the Poisson’s ratio μ is set to 0.3, as prescribed in [35].
In Figure 11, the characteristic of disc spring EN 16983 —C 50 with a rectangular cross-section and sharp edges is shown. The deviation between the characteristic computed by the original formulas of Almen (Equation (A14), Figure 11, ---) and the adjusted ones based on Almen’s formulas (Equation (A7), Figure 11, ) are caused by the correction of the lever arm, which is simplified by Almen as the difference in radius R and thus constant during deflection in the original equations. In the new approach the deflection-dependent lever arm V ψ is often (see condition (38)) smaller than the difference in radius R. Due to the use of the slightly smaller deflection-dependent lever arm V ψ instead of the difference in radius R, the calculated forces are slightly larger, compared to the results of the original formulas. For large deflections, the lever arm elongates and gets closer to the difference in radius resulting in a decreasing relative difference between the adjusted and the original characteristics for rectangular cross-sections. The same applies to the original characteristics by Curti (Equation (A15), Figure 11, ---), Zheng (Equation (A16), Figure 11, ---) and Kobelv (Equation (A24), Figure 11, ---) and the adjusted characteristics based on those (Equations (A8)–(A10), Figure 11, , and ).
The increasing force of the characteristics obtained by FE simulation (Figure 11, ) starting at about s h 0 , DIN = 0.9 is caused by the setup of the simulation where the spring is compressed between two flat surfaces. Since this effect is neglected in the analytical approaches their characteristics don’t reflect this phenomenon.
It is immediately apparent that the characteristics based on CurtiOrlando and Kobelev better match the FE simulation. The characteristics based on AlmenLaszlo and Zheng overestimate the force, while Zheng still shows a good fit for small relative deflections below 0.15.
In Figure 12 and Figure 13, the same disc spring as in Figure 11 is shown but with rounded edges ( r I , r II , r III , r IV = 0.5 mm, Figure 12) and non-rectangular cross-section ( β i , β e = 5 ° , Figure 13), respectively. In Figure 14, both variations of the cross-section are combined. To highlight the influence of the geometric variations, the FE simulation of the rectangular cross-section (Figure 11, ) is added as a black dashed line (Figure 12, ---) to the plots. Since the original formulas are based on a rectangular cross-sections, the dashed lines are equivalent to the ones of same color in Figure 11.
Figure 12 shows a substantial increase in the force with rounded edges with a 0.5 mm radius, i.e., two fifth of the sheet thickness each. The characteristics based on CurtiOrlando and Kobelev show a good fit to the FE simulation, while both show a slight underestimation around a relative deflection of 0.35. The characteristics based on AlmenLaszlo and Zheng et al. show the same divergence from the simulated characteristics as for a rectangular cross-section. This can be explained by the fact that the new approach is based on the original formulas and doesn’t account for the deviation of the characteristics computed with the original formulas from those computed by FE simulations.
In Figure 13, the angles of the inner and outer faces, β i and β e , are both set to 5°. This causes a small decrease in force, both in the FE simulation and in the computed characteristics. Both show smaller forces compared to the results from the rectangular cross-section in Figure 11. The difference between the two FE simulated characteristics is slightly larger than the difference between the analytically computed (colored, solid) ones from the two figures. This indicates that the effect from angled inner and outer faces is not completely covered by the analytical description derived in the current work.
A combination of the two cross-sectional variations from Figure 12 and Figure 13 is applied in one single disc spring geometry. The characteristics of the cross-section with both, radii of 0.5 mm and angles of the inner and outer face of 5°, are shown in Figure 14. It shows again the slight overestimation in force along with β i and β e at 5°, but the general matches between the characteristics based on CurtiOrlando and Kobelev are still good, especially for relative deflections below 0.5.
To compare the results quantitatively, the maximum relative error ( | F i , o F FE | / F FE ) was used, where F FE is the force computed by FE simulation and F i is the force computed by the new method based on formulas by different authors abbreviated as i. Relative deflections above 0.8 were not included in the analysis since the FE simulation shows a progressively increasing force for large deflections, as discussed above. The results for the CurtiOrlando and Kobelev based characteristics are shown in Table 1. At less than 8, the maximum relative errors of the characteristics of disc springs with geometric variations are in the same order of magnitude as those of springs with geometries with a rectangular cross-section and sharp edges.
In Appendix C, characteristics and maximum relative errors for more geometric variations are shown. The correction made by the adjustments is not large enough to reflect the changes in the simulated characteristic caused by negative angles β e and β i . For negative angles larger in magnitude (−10°), this divergence increases over-proportionally and leads to a better fit of the FE simulated characteristics with the results of Zheng, which is likely due to the compensation of errors. Large angles (10°) cause the characteristics based on Curti and Kobelev to overestimate the force over a wide range of relative deflections (here, 0.2 to 0.8).
The new analytical method is a computationally inexpensive alternative to FE simulations of complex cross-sectional disc spring geometries. The new approach shows good characteristic adjustments while depending on an accurate characteristic computation method for disc springs with rectangular cross-section with sharp edges. The only weaknesses not attributable to the original formula were found with large angles of the inner and outer face. While those cases cannot be generally discounted, they are unlikely to occur. For large angles, the geometry quantities should be interpreted as outlined at the end of Section 2.
The best results were obtained by using the formulas of CurtiOrlando or Kobelev as base formulas, with the formula of Kobelev showing marginally better results. The method to adapt an existing formula to special cross-sections outlined in the current work may be adapted to other existing or future methods.

7. Conclusions

This paper presents the following findings:
  • A set of formulas were introduced to describe the geometric relations of the cross-section of a disc spring with non-rectangular cross-section and rounded edges. They enable the formulation of a more accurate way to compute the lever arm, taking its elongation during deflection into account.
  • A method adapting characteristic calculation formulas and algorithms to non-rectangular cross-sections with rounded edges was outlined. This new adaption method can be applied to different traditional characteristic calculation algorithms designed for disc springs with a rectangular cross-section. The method was exemplarily applied to the algorithms of AlmenLaszlo, CurtiOrlando, Zheng and Kobelev to show its applicability.
  • The adjusted formulas were used to calculate characteristics of different disc springs. The results were compared to those of the original formulas and those obtained by FE simulations. Provided that the original characteristics approximate those of rectangular cross-sectional disc springs well, the adjusted characteristics based on them fit well for most investigated geometries. Based on the results of exemplarily computed characteristics, we recommend the use of the proposed method, utilizing formulas introduced by CurtiOrlando or Kobelev. The relative errors made with the new method are in the same order of magnitude as those of geometries with rectangular cross-section.

Author Contributions

Conceptualization, D.S.L. and M.B.G.; methodology, D.S.L. and M.B.G.; software, D.S.L. and M.B.G.; investigation, D.S.L.; resources, M.O.; writing—original draft preparation, D.S.L.; writing—review and editing, D.S.L., M.B.G., M.K. and M.O.; visualization, D.S.L.; supervision, M.O.; project administration, M.K.; funding acquisition, M.O. All authors have read and agreed to the published version of the manuscript.

Funding

This article was created as part of the research project AVIF A 309 “Bewertung des Einflusses realer Bauteilgeometrien auf die Beanspruchbarkeit von Tellerfedern anhand numerischer Simulation” (assessment of the influence of real geometries on the load capacity of disc springs by numerical simulation). This project is funded by Stiftung Stahlanwendungsforschung, which is part of Stifterverband für die Deutsche Wissenschaft e. V. (Donors’ Association for the Promotion of Science and Humanities in Germany). The Association’s mission is the promotion of research into the manufacturing and utilization of steel in Germany. The research proposal was audited by apanel of experts from the Research Association of the Working Group of the Iron- and Metal-processing Industries (AVIF), which is composed of specialists from the steelworking industry and academia. The project is accompanied by a working group from the Association of the German Spring Industry (VDFI). We acknowledge support by the German Research Foundation and the Open Access Publishing Fund of Technical University Darmstadt.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Derivation of a II through a IV and b II through b IV

Analogously to a I and b I , a II and b II can be obtained from Figure A1.
Figure A1. Geometry of edge II to calculate a II and b II ( β i > 0 ).
Figure A1. Geometry of edge II to calculate a II and b II ( β i > 0 ).
Materials 15 01954 g0a1
The following equations emerge from the relations on edge II :
a II = sin ( φ β i ) tan π 2 β i 2 tan φ β i 2 r II
b II = cos ( φ β i ) tan π 2 β i 2 tan ( φ β i ) r II
As can be seen by comparing Equations (A1) and (A2) with Equations (17) and (21), respectively, the above equations can be obtained by replacing φ by φ β i and π 2 + β i by π 2 β i as well as substituting the radius r I by r II in the equations.
Essentially, the equations of the remaining two edges III and IV can be obtained by simply rotating Figure 7, Figure 8 and Figure A1 by 180° and replacing β i by β e , r I by r III and r II by r IV .
In this way the following equations can be obtained to compute a III , b III :
a III = sin ( φ ) tan π 2 + β e 2 tan φ 2 r III
b III = cos ( φ ) tan π 2 + β e 2 tan ( φ ) r III
as well as a IV and b IV :
a IV = sin ( φ β e ) tan π 2 β e 2 tan φ β e 2 r IV
b IV = cos ( φ β e ) tan π 2 β e 2 tan ( φ β e ) r IV

Appendix B. Full Equations

Below the adapted formulas based on Almen and Laszlo, Curti and Orlando, Zheng et al. as well as Kobelev are written out.
F A = F A , o ( D e , D i , l 0 ) R V ψ
F C = F C , o ( D e , D i , l 0 ) R V ψ
F Z = F Z , o ( D e , D i , l 0 ) R V ψ
F K = F K , o ( D e , D i , l 0 ) H r , K V ψ
where H r , K = r i , K cos ( ψ K ) cos ( α K ) ( Δ K 1 )
where the argument ( D e , D i , l 0 ) denotates that the original forces need to be computed using D e , D i and l 0 instead of D e , D i and l 0 , in quantities defined using them or directly in the formulas below, e.g., h 0 , DIN is replaced by h 0 , DIN = l 0 t . Note the different calculation of D e and D i discussed in Section 5.1 for the Kobelev-based adjustment, in which the following definition should be used instead of Equations (42) and (43):
D e = D e 2 k e cos ( φ )
D i = D i + 2 k i cos ( φ )
The original forces of AlmenLaszlo [39] F A , o , CurtiOrlando [40] F C , o and Zheng [41] F Z , o are calculated as follows:
F A , o = 4 E 1 μ 2 t 3 s D e 2 1 M A · 1 + h 0 , DIN t s t h 0 , DIN t s 2 t
F C , o = 4 E 1 μ 2 t 3 s D e 2 1 M C · 1 + h 0 , DIN t s t h 0 , DIN t s 2 t
F Z , o = 4 E 1 μ 2 π α 1 , Z · 1 + α 0 , Z · h 0 , DIN t s t h 0 , DIN t s 2 t
where 1 M A = δ + 1 δ 1 2 ln δ π δ δ 1 2
1 M C = 1 μ 2 2 π 1 μ δ 2 ( δ 1 ) 3 · 1 + δ 2 + μ 1 + μ + δ μ + 1 1 1 δ μ
α 0 , Z = 1 μ 2 6 δ 1 2 ln δ · 1 4 δ 2 1 δ 2 δ 2 1 ln δ 2
α 1 , Z = 1 6 δ δ 1 2 ln δ
To compute the force F K , o with the disc spring equations by Kobelev [85], the following algorithm was used: Since the geometry definition by Kobelev is based on the radii r i , K , r e , K at the middle of the inner and outer face, these quantities must first be computed. This is done iteratively with the following three equations until convergence is achieved:
α K = arctan l 0 t r e , K r i , K
r i , K = D i 2 + sin α K × t 2
r e , K = D e 2 sin α K × t 2
According to Equations (4.33) through (4.35) in [85], those quantities can then be used to compute the force:
F K , o = π E r i , K 2 F e , K μ K + F f , K μ K 3 cos ψ K
F e , K = 2 ( 1 Δ K ) + ( 1 + Δ K ) ln ( Δ K ) ln ( Δ K ) · cos ( ψ K ) cos ( α K ) cos 2 ( α K ) / sin ( ψ K )
F f , K = ln Δ K 6 Δ K 1 · ( sin α K sin ψ K ) cos ψ K
where
μ K = t r i , K
ψ K = arcsin H K x i , K x e , K
H K = l 0 t s
x i , K = c i , K r i , K cos α K
x e , K = c i , K r e , K cos α K
c i , K = Δ K 1 ln Δ K × r i , K
Δ K = r e , K r i , K ,
according to Equations (4.4), (4.13), (4.14), (4.30), (4.31) and Equation on p. 100 in [85].

Appendix C. Force vs. Deflection Curves for Additional Geometric Variations

Below, more examples of geometric variations are shown. The maximum relative errors are shown in Table A1. The base geometry with a rectangular cross-section and sharp edges can be found in Figure 11.
Figure A2. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0 m m and β i , β e = 5 ° .
Figure A2. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0 m m and β i , β e = 5 ° .
Materials 15 01954 g0a2
Figure A3. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0 m m and β i , β e = 10 ° .
Figure A3. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0 m m and β i , β e = 10 ° .
Materials 15 01954 g0a3
Figure A4. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0 m m and β i , β e = 10 ° .
Figure A4. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0 m m and β i , β e = 10 ° .
Materials 15 01954 g0a4
Figure A5. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0.5 m m and β i , β e = 5 ° .
Figure A5. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0.5 m m and β i , β e = 5 ° .
Materials 15 01954 g0a5
Figure A6. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0.5 m m and β i , β e = 10 ° .
Figure A6. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0.5 m m and β i , β e = 10 ° .
Materials 15 01954 g0a6
Figure A7. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0.5 m m and β i , β e = 10 ° .
Figure A7. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0.5 m m and β i , β e = 10 ° .
Materials 15 01954 g0a7
Figure A8. Characteristic comparison of a disc spring with r I = 0.3 m m , r II = 0.8 m m , r III = 0.5 m m , r IV = 0.3 m m and β i , β e = 5 ° .
Figure A8. Characteristic comparison of a disc spring with r I = 0.3 m m , r II = 0.8 m m , r III = 0.5 m m , r IV = 0.3 m m and β i , β e = 5 ° .
Materials 15 01954 g0a8
Table A1. Maximum relative errors of the adjusted characteristics for additional geometric variations.
Table A1. Maximum relative errors of the adjusted characteristics for additional geometric variations.
CurtiOrlando BasedKobelev Based
β i , β e r I , , r IV = 0 mm r I , , r IV = 0.5 mm r I , , r IV = 0 mm r I , , r IV = 0.5 mm
10°9.64%4.06%10.18%4.68%
7.36%1.73%17.88%2.35%
5.09%4.22%15.58%3.69%
−5°3.98%6.37%13.37%5.87%
−10°5.41%8.47%14.82%8.00%
r I = 0.3 mm , r II = 0.8 mm , r III = 0.5 mm , r IV = 0.3 mm
1.97%1.43%

References

  1. Köhler, W. Die Tellerfeder, ein Hilfsmittel zur Klärung von Grenzflächenvorgängen. Mater. Corros. 1955, 6, 389–397. [Google Scholar] [CrossRef]
  2. Lee, C.Y.; Chae, Y.S.; Gwon, J.D.; Nam, U.H.; Kim, T.H. Finite element analysis and optimal design of automobile clutch diaphragm spring. Trans. Korean Soc. Mech. Eng. A 2000, 24, 1616–1623. [Google Scholar]
  3. Chen, Y.S.; Chiu, C.C.; Cheng, Y.D. Dynamic analysis of disc spring effects on the contact pressure of the collet—Spindle interface in a high-speed spindle system. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2009, 223, 1191–1201. [Google Scholar] [CrossRef]
  4. Bouzid, A.H.; Nechache, A. The modelling of bolted flange joints used with disc springs and tube spacers to reduce relaxation. Int. J. Press. Vessel. Pip. 2010, 87, 730–736. [Google Scholar] [CrossRef]
  5. Sfarni, S.; Bellenger, E.; Fortin, J.; Malley, M. Numerical and experimental study of automotive riveted clutch discs with contact pressure analysis for the prediction of facing wear. Finite Elem. Anal. Des. 2011, 47, 129–141. [Google Scholar] [CrossRef]
  6. Ha, S.H.; Seong, M.S.; Choi, S.B. Design and vibration control of military vehicle suspension system using magnetorheological damper and disc spring. Smart Mater. Struct. 2013, 22, 065006. [Google Scholar] [CrossRef]
  7. Castagnetti, D. A Belleville-spring-based electromagnetic energy harvester. Smart Mater. Struct. 2015, 24, 094009. [Google Scholar] [CrossRef]
  8. Aktir, Y.; Brunel, J.F.; Dufrenoy, P.; Mahé, H. Three-dimensional finite element model of an automotive clutch for analysis of axial vibrations. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2016, 230, 1324–1337. [Google Scholar] [CrossRef]
  9. Tsang, K.S.; Ion, W.; English, M. Investigation into the Numerical Model Behaviour of Belleville Washers in Cold Roll Forming. Key Eng. Mater. 2016, 716, 719–727. [Google Scholar] [CrossRef] [Green Version]
  10. Sgambitterra, E.; Maletta, C.; Furgiuele, F. Modeling and simulation of the thermo-mechanical response of NiTi-based Belleville springs. J. Intell. Mater. Syst. Struct. 2016, 27, 81–91. [Google Scholar] [CrossRef]
  11. Ramhormozian, S.; Clifton, G.C.; MacRae, G.A.; Davet, G.P. Stiffness-based approach for Belleville springs use in friction sliding structural connections. J. Constr. Steel Res. 2017, 138, 340–356. [Google Scholar] [CrossRef]
  12. Xu, L.; Fan, X.; Li, Z. Experimental behavior and analysis of self-centering steel brace with pre-pressed disc springs. J. Constr. Steel Res. 2017, 139, 363–373. [Google Scholar] [CrossRef]
  13. Krishnasamy, K.; Masse, F.; Simon, O. Fatigue analysis of Diaphragm spring in double dry clutch including manufacturing process. Procedia Eng. 2018, 213, 606–612. [Google Scholar] [CrossRef]
  14. Xu, L.; Fan, X.; Li, Z. Hysteretic analysis model for pre-pressed spring self-centering energy dissipation braces. J. Struct. Eng. 2018, 144, 04018073. [Google Scholar] [CrossRef]
  15. Xu, L.H.; Xie, X.S.; Li, Z.X. Development and experimental study of a self-centering variable damping energy dissipation brace. Eng. Struct. 2018, 160, 270–280. [Google Scholar] [CrossRef]
  16. Chang, Y.; Li, W.; Luo, H.; He, X.; Iannuzzo, F.; Blaabjerg, F.; Lin, W. A 3D thermal network model for monitoring imbalanced thermal distribution of press-pack IGBT modules in MMC-HVDC applications. Energies 2019, 12, 1319. [Google Scholar] [CrossRef] [Green Version]
  17. Fukasawa, T.; Okamura, S.; Somaki, T.; Miyagawa, T.; Uchita, M.; Yamamoto, T.; Watakabe, T.; Fujita, S. Research and Development of Three-Dimensional Isolation System for Sodium-Cooled Fast Reactor: Part 4—Proposal of Optimal Combination Method for Disc Spring Units and Newly Friction Model for Sliding Elements. In Pressure Vessels and Piping Conference; American Society of Mechanical Engineers: New York, NY, USA, 2019; p. V008T08A038. [Google Scholar]
  18. Shou, M.; Liao, C.; Zhang, H.; Xie, L. A design methodology based on full dynamic model for magnetorheological energy absorber equipped with disc springs. Smart Mater. Struct. 2019, 28, 065020. [Google Scholar] [CrossRef]
  19. Sapountzakis, I.E.; Tranakidis, P.G.; Antoniadis, I.A. Implementation of the KDamper concept using disc springs. J. Low Freq. Noise Vib. Act. Control 2019, 38, 168–186. [Google Scholar] [CrossRef]
  20. Sano, T.; Shirai, K.; Suzui, Y.; Utsumi, Y. Loading tests of a brace-type multi-unit friction damper using coned disc springs and numerical assessment of its seismic response control effects. Bull. Earthq. Eng. 2019, 17, 5365–5391. [Google Scholar] [CrossRef]
  21. Sun, Y.; Gong, D.; Zhou, J.; Sun, W.; Xia, Z. Low frequency vibration control of railway vehicles based on a high static low dynamic stiffness dynamic vibration absorber. Sci. China Technol. Sci. 2019, 62, 60–69. [Google Scholar] [CrossRef]
  22. Huang, X.; Su, Z.; Wang, S.; Wei, X.; Wang, Y.; Hua, H. High-frequency disturbance force suppression mechanism of a flywheel equipped with a flexible dynamic vibration absorber. J. Vib. Control 2020, 26, 2113–2124. [Google Scholar] [CrossRef]
  23. Xiao, S.J.; Xu, L.H.; Li, Z.X. Design and experimental verification of disc spring devices in self–centering reinforced concrete shear walls. Struct. Control Health Monit. 2020, 27, e2549. [Google Scholar] [CrossRef]
  24. Xiao, S.J.; Xu, L.H.; Li, Z.X. Development and experimental verification of self-centering shear walls with disc spring devices. Eng. Struct. 2020, 213, 110622. [Google Scholar] [CrossRef]
  25. Zhou, J.; Zhang, C.; Wang, Z.; Mao, K. Study on Dynamic Characteristics of the Disc Spring System in Vibration Screen. Shock Vib. 2020, 2020, 3518037. [Google Scholar] [CrossRef]
  26. Xu, L.H.; Xiao, S.J. A performance-based section design method of a self-centering concrete shear wall with disc spring devices. Eng. Mech. 2020, 37, 70–77. [Google Scholar]
  27. Menne, H.; Franck, C. Benchmarking the Suitability of a Bi-Stable Disc Spring as Novel Ultra-Fast Switchgear Actuation Principle. In 2020 CIGRE e-session. Papers and Proceedings; International Council on Large Electric Systems: Paris, France, 2020; p. A3302. [Google Scholar]
  28. Pal, U.; Palit, P.; Gokarn, P.; Kanrar, S.; Kumar, A. Improvement of wagon tippler disc spring life by design modification. Eng. Fail. Anal. 2021, 120, 105030. [Google Scholar] [CrossRef]
  29. Fang, C.; Wang, W.; Shen, D. Development and experimental study of disc spring–based self-centering devices for seismic resilience. J. Struct. Eng. 2021, 147, 04021094. [Google Scholar] [CrossRef]
  30. Mrázek, M.; Sedláček, F.; Skovajsa, M. Design of Composite Disc Spring for Automotive Suspension with using Numerical Simulation. Manuf. Technol. 2021, 21, 1213–2489. [Google Scholar] [CrossRef]
  31. Holtmann, C.; Köhler, C. (Eds.) Eddy Current and Friction Brake Fusion by a Disc-Spring to a Hybrid-Brake; International Council on Large Electric Systems (Cigré): Paris, France, 2021. [Google Scholar]
  32. Chen, R.; Li, X.; Xu, J.; Yang, Z.; Yang, H. Properties Analysis of Disk Spring with Effects of Asymmetric Variable Friction. Int. J. Appl. Mech. 2021, 13, 2150076. [Google Scholar] [CrossRef]
  33. Kumar, S. Design and Analysis of Weld Bead Cutting Tool Using FEA Technique with Experimental Verification of Cutting Force. In Advances in Mechanical Processing and Design; Springer: Singapore, 2021; pp. 581–592. [Google Scholar]
  34. Hong, G.; Zhongmin, D.; Yanlin, Z.; Hongbo, Y.; Xinjie, Z.; Lingzi, M.; Qi, L. Time-delayed feedback control of nonlinear dynamics in a giant magnetostrictive actuator. Nonlinear Dyn. 2022, 1–24. [Google Scholar] [CrossRef]
  35. DIN EN 16983:2017-09; Tellerfedern—Qualitätsanforderungen—Maße; Deutsche Fassung EN 16983:2016. DIN Deutsches Institut für Normung: Berlin, Germany, 2017.
  36. DIN EN 16984:2017-09; Tellerfedern—Berechnung; Deutsche Fassung EN 16984:2016. DIN Deutsches Institut für Normung: Berlin, Germany, 2017.
  37. Geilen, M.B.; Klein, M.; Oechsner, M. Spring_stack—Ein Modul zur numerischen Simulation von Tellerfedern und Tellerfedersäulen. In Ilmenauer Federntag 2019—Neueste Erkenntnisse zu Funktion, Berechnung, Prüfung und Gestaltung von Federn und Werkstoffen; ISLE: Ilmenau, Germany, 2019; pp. 77–86. [Google Scholar]
  38. Geilen, M.B.; Klein, M.; Oechsner, M. A Novel Algorithm for the Determination of Walker Damage in Loaded Disc Springs. Materials 2020, 13, 1661. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Almen, J.O.; Laszlo, A. The Uniform-Section Disc Spring. Trans. Am. Soc. Mech. Eng. 1936, 58, 305–314. [Google Scholar]
  40. Curti, G.; Orlando, M. Ein neues Berechnungsverfahren für Tellerfedern. Draht 1997, 30, 17–22. [Google Scholar]
  41. Zheng, E.; Jia, F.; Zhou, X. Energy-based method for nonlinear characteristics analysis of Belleville springs. Thin-Walled Struct. 2014, 79, 52–61. [Google Scholar] [CrossRef]
  42. Kobelev, V. Exact shell solutions for conical springs. Mech. Based Des. Struct. Mach. 2016, 44, 317–339. [Google Scholar] [CrossRef]
  43. Luxenburg, Y.; Givli, S. The Static Response of Axisymmetric Conical Shells Exhibiting Bistable Behavior. J. Appl. Mech. 2021, 88. [Google Scholar] [CrossRef]
  44. Curti, G.; Montanini, R. On the influence of friction in the calculation of conical disk springs. J. Mech. Des. 1999, 121, 622–627. [Google Scholar] [CrossRef]
  45. Ozaki, S.; Tsuda, K.; Tominaga, J. Analyses of static and dynamic behavior of coned disc springs: Effects of friction boundaries. Thin-Walled Struct. 2012, 59, 132–143. [Google Scholar] [CrossRef]
  46. Mastricola, N.P.; Dreyer, J.T.; Singh, R. Analytical and experimental characterization of nonlinear coned disk springs with focus on edge friction contribution to force-deflection hysteresis. Mech. Syst. Signal Process. 2017, 91, 215–232. [Google Scholar] [CrossRef]
  47. Chen, R.; Li, X.; Yang, Z.; Xu, J.; Yang, H. Nonlinear behavior of disk spring with complex contact state. Sci. Prog. 2021, 104, 00368504211052360. [Google Scholar] [CrossRef]
  48. Li, X.; Chen, R.; Yang, Z.; Yang, H.; Xu, J. Static behavior analysis of disc spring considering variable static friction coefficient. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2021, 235, 5583–5593. [Google Scholar] [CrossRef]
  49. Schremmer, G. Die geschlitzte Tellerfeder. Konstruktion 1972, 24, 226–229. [Google Scholar]
  50. Curti, G.; Orlando, M. Geschlitzte Tellerfedern. Draht 1981, 32, 608–615. [Google Scholar]
  51. Fawazi, N.; Lee, J.Y.; Oh, J.E. A load-displacement prediction for a bended slotted disc using the energy method. Proc. Inst. Mech. Engeneers, Part C J. Mech. Eng. Sci. 2012, 226, 2126–2137. [Google Scholar] [CrossRef]
  52. Wang, J.; Jian, L.; Leng, W.; She, C.; Sun, Z. An Analytical Model for the Stiffness of Slotted Disk Springs. Int. J. Appl. Mech. 2020, 12, 2050120. [Google Scholar] [CrossRef]
  53. La Rosa, G.; Messina, M.; Risitano, A. Stiffness of Variable Thickness Belleville Springs. J. Mech. Des. 2001, 123, 294. [Google Scholar] [CrossRef]
  54. Chaturvedi, R.; Trikha, M.; Simha, K.R. Theoretical and numerical analysis of stepped disk spring. Thin-Walled Struct. 2019, 136, 162–174. [Google Scholar] [CrossRef]
  55. Du, X.; Liao, C.; Gan, B.; Zhang, Y.; Xie, L.; Zhang, H. Analytical modeling and experimental verification for linearly gradient thickness disk springs. Thin-Walled Struct. 2021, 167, 108153. [Google Scholar] [CrossRef]
  56. Muhr, K.H.; Niepage, P. Zur Berechnung von Tellerfedern mit rechteckigem Querschnitt und Auflageflächen. Konstruktion 1966, 18, 24–27. [Google Scholar]
  57. Ferrari, G. A New Calculation Method for Belleville Disc Springs with Contact Flats and Reduced Thickness. Int. J. Manuf. Mater. Mech. Eng. 2013, 3, 63–73. [Google Scholar] [CrossRef]
  58. Curti, G.; Appendino, D. Vergleich von Berechnungsverfahren für Tellerfedern. Draht 1982, 33, 38–40. [Google Scholar]
  59. Schiffner, K.; Dietrich, M. Numerische Simulation der Entstehung von Eigenspannungen 1. Art an Tellerfedern. ZAMM 1987, 67, 235–237. [Google Scholar]
  60. Wagner, W.; Wetzel, M. Berechnung von Tellerfedern mit Hilfe der Methode der Finiten Elemente. Konstruktion 1987, 39, 147–150. [Google Scholar]
  61. Bagavathiperumal, P.; Chandrasekaran, K.; Manivasagam, S. Elastic load-displacement predictions for coned disc springs subjected to axial loading using the finite element method. J. Strain Anal. 1991, 26, 147–152. [Google Scholar] [CrossRef]
  62. Curti, G.; Raffa, F.A. Material nonlinearity effects in the stress analysis of conical disk springs. J. Mech. Des. 1992, 114, 238–244. [Google Scholar] [CrossRef]
  63. Yahata, N.; Watanabe, M.; Ishii, N. Analysis of coned disc spring by finite element method. Jpn. Soc. Mech. Eng. 1993, 59, 260–265. [Google Scholar] [CrossRef] [Green Version]
  64. Curti, G.; Montanini, R. Theoretical, Numerical and Experimental Analysis of Conical Disk Spring. XXV AIAS Gall. Lecce Sept. 1996, 1, 573. [Google Scholar]
  65. Lee, B.S. Load Characteristics and Sensitivity Analysis for an Automotive Clutch Diaphragm Spring. Trans. Korean Soc. Automot. Eng. 2006, 14, 54–59. [Google Scholar]
  66. Shen, W.; Fang, W. Design of a Friction Clutch Using Dual Belleville Structures. J. Mech. Des. 2007, 129, 986. [Google Scholar] [CrossRef]
  67. Berger, C.; Veleva, D. AVIF A 244: Experimentelle Untersuchung und numerische Simulation des Relaxationsverhaltens von Tellerfedern. In Abschlussbericht; TU Darmstadt: Darmstadt, Germany, 2010. [Google Scholar]
  68. Veleva, D.; Kaiser, B.; Berger, C. Experimentelle Untersuchung und numerische Simulation des Relaxationsverhaltens von Tellerfedern. In Ilmenauer Federntag 2010; ISLE: Illmenau, Germany, 2010; pp. 71–78. [Google Scholar]
  69. Wehmann, C.; Neidnicht, M.; Nützel, F.; Roith, B.; Rieg, F. Finite Elemente Analysen zur Berechnung von Maschinenelementen mit nichtlinearem Verhalten. In 8. Gemeinsames Kolloquium Konstruktionstechnik; Docupoint: Magdeburg, Germany, 2010; pp. 47–52. [Google Scholar]
  70. Karakaya, Ş. Investigation of hybrid and different cross-section composite disc springs using finite element method. Trans. Can. Soc. Mech. Eng. 2012, 36, 399–412. [Google Scholar] [CrossRef]
  71. Veleva, D. Experimentelle Untersuchung und numerische Simulation des Relaxationsverhaltens von Tellerfedern. Ph.D. Thesis, TU Darmstadt, Darmsatdt, Germany, 2012. [Google Scholar]
  72. Maletta, C.; Filice, L.; Furgiuele, F. NiTi Belleville washers: Design, manufacturing and testing. J. Intell. Mater. Syst. Struct. 2013, 24, 695–703. [Google Scholar] [CrossRef]
  73. Jia, F.; Zhang, F.C. A Study on the Mechanical Properties of Disc Spring Vibration Isolator with Viscous Dampers. Adv. Mater. Res. 2014, 904, 454–459. [Google Scholar] [CrossRef]
  74. Fang, C.; Zhou, X.; Osofero, A.I.; Shu, Z.; Corradi, M. Superelastic SMA Belleville washers for seismic resisting applications: Experimental study and modelling strategy. Smart Mater. Struct. 2016, 25, 105013. [Google Scholar] [CrossRef] [Green Version]
  75. Lonkwic, P.; Różyło, P.; Usydus, I. Numerical Analysis and Experimental Investigation of Disk Spring Configurations with Regard to Load Capacity of Safety Progressive Gears. Appl. Comput. Sci. 2016, 12, 5–16. [Google Scholar]
  76. Baak, N.; Garlich, M.; Schmiedt, A.; Bambach, M.; Walther, F. Characterization of residual stresses in austenitic disc springs induced by martensite formation during incremental forming using micromagnetic methods. Mater. Test. 2017, 59, 309–314. [Google Scholar] [CrossRef] [Green Version]
  77. Yubing, G.; Defeng, Z. Mechanism for the Forced Strengthening on the Diaphragm Spring’s Load-Deflection Characteristic. Int. J. Eng. Technol. 2017, 9, 287–292. [Google Scholar] [CrossRef] [Green Version]
  78. Fang, C.; Yam, M.C.H.; Chan, T.M.; Wang, W.; Yang, X.; Lin, X. A study of hybrid self-centring connections equipped with shape memory alloy washers and bolts. Eng. Struct. 2018, 164, 155–168. [Google Scholar] [CrossRef]
  79. Zala, R.G.; Gawande, K.R. Disc Spring Washer: Design and Axisymmetric analysis using ANSYS. Trends Mach. Des. 2018, 4, 21–25. [Google Scholar]
  80. Zhu, D.; Ding, F.; Liu, H.; Zhao, S.; Liu, G. Mechanical property analysis of disc spring. J. Braz. Soc. Mech. Sci. Eng. 2018, 40, 230. [Google Scholar] [CrossRef]
  81. Koplin, C.; Schröder, C.; Beckert, W.; Stockmann, J.; Kailer, A. Demonstration der Zuverlässigkeit von keramischen Schraubendruck- und Tellerfedern für korrosive Umgebungen und hohe Temperaturen. In Ilmenauer Federntag 2019; ISLE: Illmenau, Germany, 2019; pp. 11–21. [Google Scholar]
  82. Khader, I.; Koplin, C.; Schröder, C.; Stockmann, J.; Beckert, W.; Kunz, W.; Kailer, A. Characterization of a silicon nitride ceramic material for ceramic springs. J. Eur. Ceram. Soc. 2020, 40, 3541–3554. [Google Scholar] [CrossRef]
  83. Zhou, J.; Zhang, C.; Wang, Z.; Mao, K.; Wang, X. A Study on the Influence of Different Constraint Modes and Number of Disc Springs on the Dynamics of Disc Spring System. Shock Vib. 2021, 2021, 8866159. [Google Scholar] [CrossRef]
  84. Niepage, P. Vergleich verschiedener Verfahren zur Berechnung von Tellerfedern—Teil I. Draht 1983, 34, 105–108; 255–258. [Google Scholar]
  85. Kobelev, V. Durability of Springs, 2nd ed.; Springer: Cham, Switzerland, 2021. [Google Scholar]
Figure 1. Half-section of a disc spring with rectangular cross-section with sharp edges.
Figure 1. Half-section of a disc spring with rectangular cross-section with sharp edges.
Materials 15 01954 g001
Figure 2. Dimensioned right half of a disc spring section.
Figure 2. Dimensioned right half of a disc spring section.
Materials 15 01954 g002
Figure 3. Dimensioned right half of a disc spring section with non-rectangular cross-section and rounded edges.
Figure 3. Dimensioned right half of a disc spring section with non-rectangular cross-section and rounded edges.
Materials 15 01954 g003
Figure 4. Definition of an idealized cross-section of a real disc spring ( β i , β e < 0 ).
Figure 4. Definition of an idealized cross-section of a real disc spring ( β i , β e < 0 ).
Materials 15 01954 g004
Figure 5. Definition of the geometric quantities for β i , β e > φ .
Figure 5. Definition of the geometric quantities for β i , β e > φ .
Materials 15 01954 g005
Figure 6. Cross-section with bounding boxes to calculate the initial slope angle φ and the lever arm V ( β i , β e < 0 ).
Figure 6. Cross-section with bounding boxes to calculate the initial slope angle φ and the lever arm V ( β i , β e < 0 ).
Materials 15 01954 g006
Figure 7. Geometry of edge I to calculate a I ( β i > 0 ). The dashed lines show the geometry with β i = 0 .
Figure 7. Geometry of edge I to calculate a I ( β i > 0 ). The dashed lines show the geometry with β i = 0 .
Materials 15 01954 g007
Figure 8. Geometry of edge I to calculate b I ( β i > 0 ). The dashed lines show the geometry with β i = 0 .
Figure 8. Geometry of edge I to calculate b I ( β i > 0 ). The dashed lines show the geometry with β i = 0 .
Materials 15 01954 g008
Figure 9. Cross-section of a unloaded (black) and loaded (blue) disc spring.
Figure 9. Cross-section of a unloaded (black) and loaded (blue) disc spring.
Materials 15 01954 g009
Figure 10. Comparison of an idealized cross-section (black, β i , β e < 0 ), a rectangular cross-section based on the basic geometry parameters (green, same R and l 0 as black but a rectangular cross-section and sharp edges, the cross-section area and the slope angle differs) and a rectangular cross-section (orange) to have similar properties as the idealized cross-section (black). The dashed blue lines show the cross-section after the second shape adjustment with the corrections k i and k e to keep the cross-section area approximately the same. Its slope angle is the same as the orange and the black one’s.
Figure 10. Comparison of an idealized cross-section (black, β i , β e < 0 ), a rectangular cross-section based on the basic geometry parameters (green, same R and l 0 as black but a rectangular cross-section and sharp edges, the cross-section area and the slope angle differs) and a rectangular cross-section (orange) to have similar properties as the idealized cross-section (black). The dashed blue lines show the cross-section after the second shape adjustment with the corrections k i and k e to keep the cross-section area approximately the same. Its slope angle is the same as the orange and the black one’s.
Materials 15 01954 g010
Figure 11. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0 mm and β i , β e = 0 ° .
Figure 11. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0 mm and β i , β e = 0 ° .
Materials 15 01954 g011
Figure 12. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0.5 mm and β i , β e = 0 ° .
Figure 12. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0.5 mm and β i , β e = 0 ° .
Materials 15 01954 g012
Figure 13. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0 mm and β i , β e = 5 ° .
Figure 13. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0 mm and β i , β e = 5 ° .
Materials 15 01954 g013
Figure 14. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0.5 mm and β i , β e = 5 ° .
Figure 14. Characteristic comparison of a disc spring with r I , r II , r III , r IV = 0.5 mm and β i , β e = 5 ° .
Materials 15 01954 g014
Table 1. Maximum relative errors of the adjusted characteristics.
Table 1. Maximum relative errors of the adjusted characteristics.
CurtiOrlando BasedKobelev Based
β i , β e r I , , r IV = 0 mm r I , , r IV = 0.5 mm r I , , r IV = 0 mm r I , , r IV = 0.5 mm
5.09%4.22%5.58%3.69%
7.36%1.73%7.88%2.35%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Leininger, D.S.; Geilen, M.B.; Klein, M.; Oechsner, M. A New Method for the Calculation of Characteristics of Disc Springs with Trapezoidal Cross-Sections and Rounded Edges. Materials 2022, 15, 1954. https://doi.org/10.3390/ma15051954

AMA Style

Leininger DS, Geilen MB, Klein M, Oechsner M. A New Method for the Calculation of Characteristics of Disc Springs with Trapezoidal Cross-Sections and Rounded Edges. Materials. 2022; 15(5):1954. https://doi.org/10.3390/ma15051954

Chicago/Turabian Style

Leininger, Dominik Sebastian, Max Benedikt Geilen, Marcus Klein, and Matthias Oechsner. 2022. "A New Method for the Calculation of Characteristics of Disc Springs with Trapezoidal Cross-Sections and Rounded Edges" Materials 15, no. 5: 1954. https://doi.org/10.3390/ma15051954

APA Style

Leininger, D. S., Geilen, M. B., Klein, M., & Oechsner, M. (2022). A New Method for the Calculation of Characteristics of Disc Springs with Trapezoidal Cross-Sections and Rounded Edges. Materials, 15(5), 1954. https://doi.org/10.3390/ma15051954

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop