Titel: Beitrag zur Berechnung zylindrischer Schraubenfedern.
Autor: Max Pilgram
Fundstelle: Band 337, Jahrgang 1922, S. 21
Download: XML
Beitrag zur Berechnung zylindrischer Schraubenfedern. Von Regierungsbaurat z. D. Max Pilgram.*) PILGRAM, Beitrag zur Berechnung zylindrischer Schraubenfedern. I. Allgemeine Voraussetzungen. Den üblichen Formeln der Hütte liegt die Annahme zugrunde, daß in den Querschnitten der Federn sich die Spannungen nach demselben Gesetz verteilen, wie im geraden Stabe gleichen Querschnitts. Ueber dieses Gesetz wird eine willkürliche Annahme gemacht. Eine andere zuweilen verwendete Methode benutzt das Gesetz für die Spannungsverteilung auf Grund der Theorie von de St. Vénant. In der vorliegenden Abhandlung wird nun vorausgesetzt, daß in allen Querschnitten der Feder sich die Spannungen nach demselben, vorläufig unbekannten Gesetz verteilen, und es wird festgestellt, welchen Einfluß der Wickelungshalbmesser auf dieses Gesetz und damit auf Formänderung und Beanspruchung hat. Auch wird über die Normalspannungen eine Theorie aufgestellt. Die Belastung ist achsial. Bei allen diesen Methoden wird eine kleine Steigung der Feder und Proportionalität zwischen Formänderungen und Spannungen sowie ein homogener Werkstoff vorausgesetzt. Zur Erzielung brauchbarer Formeln müssen die auftretenden Bessel'schen Funktionen durch ex bezw. e– x ersetzt werden. Der hierdurch entstehende Fehler beträgt weniger als 1,5 v. H., wenn das Verhältnis ag3 gewählt wird, s. Abb. 1.
[Textabbildung Bd. 337, S. 21]
Abb. 1.
Die Rechnung wird zunächst allgemein durchgeführt und auf den rechteckigen und runden Querschnitt angewandt. Es genügt, wenn nur der Fall berücksichtigt wird, wo hg1 ist, da beim stehenden Rechteck, hg<1, infolge der geringen radialen Querschnittausdehnung ohne weiteres nach de St. Vénant gerechnet werden kann.
[Textabbildung Bd. 337, S. 21]
Abb. 2.
II. Die Formänderung der Feder. Die Feder werde durch eine zentrale achsiale Last P auf die Steigung Null zusammengedrückt, und es werde dann aus der belasteten Feder ein Segment mit dem Zentriwinkel α herausgeschnitten. Es ist zu untersuchen, wie dieses Segment vor der Belastung ausgesehen haben muß. Es kann offenbar angenommen werden, daß eine kreisförmige Faser vom Querschnitt dr . dy und dem Halbmesser r im unbelasteten Zustande schraubenförmig verlief. Ferner war der Radius vor der Formänderung gleich r – u und der Zentriwinkel gleich α – Cα . α. Endlich wird sich der Abstand y, Abb. 1, um den Betrag v geändert haben. Die Größen u und v sind streng genommen Funktionen von r und y, während Cα konstant ist. Endlich hatten die Querschnitte AB, Abb. 1, vor der Formänderung eine gewölbte Form A1Bl, d.h. die kreis- bezw. schraubenförmigen Fasern haben sich ohne Längenänderung gegeneinander verschoben, wodurch im Sinne der Spannung Tr siehe Abb. 1 und 3, eine Gleitung entsteht: ϑ=[δφδrϑ0],       =rδδr(φr), und im Sinne von Ty sinngemäß die Gleitung δφδy.
[Textabbildung Bd. 337, S. 21]
Abb. 3.
Der Uebergang aus der Schraubenform in die Kreisform bedingt im Sinne der Spannung Tr eine weitere Gleitung f2rπ=cr, wenn f die Ganghöhe der Schraubenlinie. Man erhält demnach 1.      Tr=G(rδδr(φr)+cr) 2.      Ty=G(δφδy).
III. Die aus den Formänderungen abzuleitenden Spannungen.
[Textabbildung Bd. 337, S. 22]
Abb. 4.
Aus den vorstehend angenommenen Formänderungen ergeben sich Dehnungen und Gleitungen im Sinne der Spannungen T, R, Y, S, Tr, Ty, s. Abb. 3, (die Spannungen treten paarweise auf, im Interesse der Uebersichtlichkeit ist in Abb. 3 jede Spannung nur einmal eingetragen), die also nach Einführung eines Elastizitäts- bzw. Gleitmoduls E und G auf diese Formänderungen zurückzuführen sind. Die Spannungen werden in Richtung von r, y, a nach den Gleichgewichtsbedingungen untersucht, wobei 3 Gleichungen entwickelt werden. Hierbei ist in Abbild. 3 zu berücksichtigen, daß die Spannungen mit zunehmenden Koordinaten r und y ebenfalls einen Zuwachs erfahren, dagegen nicht mit wachsendem a, da die Spannungen vom Winkel a unabhängig sein sollen.
IV. Umformung der entwickelten Hauptgleichungen und Aufstellung der Grenzbedingungen. In tangentialer Richtung ergibt das Gleichgewicht 1r2δδr(r2Tr)+δTyδy=0 und mit Benutzung der Gl. 1 und 2 1r2δδrr3δδrφr+δ2φδy2=0 oder umgeformt, wenn r3δδrφr=δψδy, r3δδyφr=δψδrcr 3.      δ2ψδy2+δ2ψδr23rδψδr=2c; c=f2π, also nach Gleichung 1 und 2 4.      Ty=G1r2δψδr, Tr=Gr2δψδy Am Rande des Querschnitts muß, wenn y = f(r) die Gleichung der Randkurve, Abb. 4, auf Grund des Gleichgewichts in tangentialem Sinne – Ty dr + Tr f'(r) dr = 0 also      δψ\delrardr+δψdydy=0 Die Integration des vollständigen Differentials ergibt am Rande 5.      ψ = 0. Die Integrationskonstante kann Null gesetzt werden. Mit Hülfe des Gaußschen Integralsatzes wird bewiesen, daß die Momentensumme aller Spannungen Tr, Ty, wenn 0' als Bezugspunkt gewählt wird, verschwindet, ebenso die Summe aller wagerechten Spannungen Tr, während die Summe aller senkrechten Spannungen Ty gleich der Federbelastung P sein muß. Intergriert man also über den Querschnitt, so ist 6.      ∫∫ dr dy Ty P. Das Gleichgewicht der übrigen Spannungen R, S, Y, T führt auf zwei partielle Differentialgleichungen für u und v, deren allgemeine Lösung mit Berücksichtigung der Grenzbedingungen unmöglich ist. Es wird im nächsten Abschnitt ausführlich dargelegt, daß man mit großer Annäherung u = C (konstant) setzen und v vernachlässigen kann, womit zugleich die Spannungen Y, R, T verschwinden, und es bleibt 7.      S=E[c2r2+Cα+Cr]; c=f2π. c22r2 bedeutet die Verkürzung einer Schraubenlinie mit der Ganghöhe f und dem Halbmesser r auf einen Kreis gleichen Halbmessers: 1rα[rα(rα)2+(fα2π)2]=11(c)2r=∼c22r2. Die Konstanten in Gl. 7 ergeben sich aus: 8.       M = ∫∫ S dr dy . r = 0 9.              ∫∫ S dr dy         0, da sowohl das Moment aller Spannungen S in bezug auf die Federachse, als auch die Resultierende von S verschwindet. V. Die Berechnung der Spannungen R, Y, S, T mittelst Annäherungsverfahrens. Die Spannung S kann als Funktion nur von aufgefaßt werden. Die Integration ist zwar sehr umständlich, macht aber keine grundsätzlichen Schwierigkeiten und ergibt: 1. Für den rechteckigen Querschnitt, nenn ba=ζ: 10.      Cac2=1ζζ212ζlogζ2logζ(ζ+1)4(ζ1) 11.      Cαac2=1(ζ1)2ζlog2ζζ2(ζ1)2(ζ21)logζ C ist die Zunahme des Wickelungshalbmessers, Cα die des Zentriwinkels 1, also Cα . 2π die Zunahme des Zentriwinkels einer Windung. Im Grenzfalle, ζ = 1, ist l=a+b2=∼a 12.      C1c2=1. 13.      Cαl2c2=12. 2. Für den kreisförmigen Querschnitt: l=a+b2, δ=g2l s. Abb. 5
[Textabbildung Bd. 337, S. 22]
Abb. 5.
In Gleichung 7 ist einzusetzen: 14.      Clc2=11δ2=∼1+δ22 15. Cαl2c2=1δ2(111δ2)=∼12(1+3δ24) Im Grenzfalle, δ = 0, wird 16. Clc2=1, Cαl2c2=112, vergl. Gl. 12 u. 13. Die Querschnittsform hat also im Grenzfalle keinen Einfluß.
VI. Die Berechnungen der Schubspannungen im rechteckigen Querschnitt. Der Gleichung 3 genügt die Reihe 17. ψ=ΣAncosnπgyv(nπrg)+c(g24y2) wenn v ein Integral von 18.      v(x)3xv(x)v=0. Bedingung ist, daß ψ am Rande des Querschnittes verschwindet. Für y ± g/2 verschwindet ψ, nenn n eine beliebige ungerade Zahl. Da nun ψ auch für r = a und r = b verschwinden soll, Gl. 5, so müssen die Konstanten A nach der Theorie der Fourierschen Reihen passend bestimmt werden. Setzt man 19.      v(aπng)=v(bπng)=v so wird 20.      An.v0=8gsinnpi2g3n3π3c. Wenn I und K die beiden Integrale von I(x)+I(x)xI=0, so ist mit einem konstanten B 21.      v = 2xl' – x2l + B(2xK' – x2K). Setzt man x=nπag=α und x=nπbg=β, so läßt sich mit Hilfe der Gl. 19 B ermitteln. Dann erhält man aus Gl. 4 und 6 22.  TyG=8cgr2Σsinnπ21n2π2g21v0cosnπygv(nπrg). 23.  TrG=8cgr2Σsinnπ21n2π1g21v0sinnπygv(nπrg)+2cyr2. 24.      PG=16cgsin2nπ21n2π2g2w0v0. 25.      nenn w0=r=ar=bv(x)dxx2=[I(nπrg)+B.K.(nπrg)]r=ar=b Die Gleichungen sind praktisch nicht verwendbar, wenn man nicht vereinfacht: I(x)=exx, K(x)=exx, was mit einem Fehler von weniger als 1,5 v. H. möglich ist, wie nachgewiesen wird, sofern: a ≧ 3 g, Abb. 1. VII. Die Berechnung der Belastung P bei gegebener Durchbiegung f = 2πc. Setzt man γ=nπhg, α=nπag, β=nπbg so ist: PGc=16.g3g.π3(a+b)abΣ11n3 2(1+e2γ)γ(1e2γ)8aba+beγ(α+2)(β2)e2γ(2)(β+2); hier ist λ=hg1, γπ, γ ≧ π und man kann setzen 26. PGgc16π3g(a+b)abλpi8+8eλπaba+b(aπg+2)(bπg2)+λg2abπΣ31n42g2abπ2Σ31n5 wo Σ1n3=0,0147, Σ1n5=0,004.. Es wird zahlenmäßig bewiesen, daß angenähert, sofern hg5, gesetzt werden darf: 27.      PGgc=16π3g(a+b)ab(3,19λ2+4eλπ)(aπg+2)(bπg2) 28.      PGgc=32π5(g2ab)23(3,19λ2+4eλπ). VIII. Die Berechnung der größten achsialen Schubspannung Ty bei gegebener Durchbiegung f = 2πc. Die größte Spannung Ty ist den Koordinaten r = a, y = 0, Abb. 1, zugeordnet. Die Formeln werden in ähnlicher Weise vereinfacht, wie im vorigen Abschnitt, und man erhält 29.      TymaxgGc=4g2a2π2 [2aπg+1aπg+14abeλπa2π2g21(aπg+2)(bπg2)0,175] Die Formel ist ersetzbar durch 30.      TymaxgGc=4g2a2π2[1,8254abeλπ3g4a], wofür der Beweis zahlenmäßig für 1 ≦ λ ≦ 5 erbracht wird. IX. Die Berechnung der größten radialen Schubspannung Tr bei gegebener Durchbiegung f = 2πc. Dem größten Wert von Tr ist y=+g2 zugeordnet, dagegen macht es Schwierigkeiten, den zugehörigen Wert von r zu bestimmen. Es würde unzulässig sein, für r den mittleren Halbmesser l=a+b2 zu setzen. Dazu kommt, daß die Reihen schlecht konvergieren, wenn sich r dem Randwerte ra nähert. Es sei T'r der Wert von Tr, wenn y ± g/2. Dann stellt sich heraus, daß T' . r2 für einen gegebenen Querschnitt vom Wickelungshalbmesser nur wenig beeinflußt wird. Für ab= ergibt sich 31.      Tr2rGgc=18π2Σ1n21eγeγeγ. (enπg(rl)enπg(rl)) und diese Gleichung gilt allgemein für jeden Wert von ag. Da auch diese Gleichung nicht auf ein Maximum untersucht werden kann, wird an Hand umfangreichen Zahlenmaterials für 1 ≦ λ ≦ 5 ermittelt, daß 32.      Trmax=(11,45λπ)(0,08(λ21)ag+1g2ab) X. Unmittelbare Ableitung der Theorie de St. Vénant. Die bisher ermittelten Formeln gehen in der Tat in die Ergebnisse der Theorie von de St. Vénant über, wenn man zur Grenze übergeht und ag= bzw. lg= setzt. Die Theorie des Verfassers wird dadurch gestützt. XI. Die Berechnung von P und Tymax im kreisförmigen Querschnitt bei gegebener Durchbiegung f = 2πc. Eine geschlossene Lösung der Gleichung 3 mit Innehaltung der Grenzbedingung 5 ist nicht möglich, sie läßt sich aber sehr leicht durchführen, wenn man geringe Abweichungen von der Kreisform zuläßt. Die Genauigkeit dieser Abweichungen wird zahlenmäßig untersucht. Man kann mit Benutzung von Konstanten A, B, C setzen ψc=A+2Br2+Cr4+4\By2y2 und zwar verschwindet ψ am Rande, wenn y2=A+2Br2+Cr24B1 die Gleichung der Randkurve ist, die bei passender Wahl der Konstanten der des Kreises sehr nahe kommt. Die Durchführung der ziemlich umständlichen Integration liefert 33.      PGc=g2π4l2g24116l2+32g2=∼g3π32l3 n Uebereinstimmung mit der Hütte. Die größte achsiale Schubspannung ergibt sich für r = a 34.      TymaxgGc=lg22al2(1+3g232l2)=∼g22al während nach der Hütte TymaxgGc=g22l2. XII. Zusammenstellung der Formeln. Bezeichnungen: Es bedeutet P die achsiale Federbelastung, f = 2πc die achsiale Verkürzung, C die Zunahme des Halbmessers, Cα . 2π die Zunahme des Zentriwinkels einer Windung, a den inneren, b den äußeren, l=b+a2 den mittleren Halbmesser, h = b – a die radiale, g die achsiale Querschnittsausdehnung, λ=hg; S die Normalspannung, Trmax und Tymax die größte achsiale bzw. radiale Beanspruchung. Es kommen in Frage die Gleichungen 7, 10, 11, 14, 15, 27, 28, 29, 30, 32, 33, 34. XIII. Die Formeln des Verfassers gehen in die Formeln über, die sich aus der Theorie von de St. Venant ergeben würde, wenn man ag bzw. lg bzw. bg unendlich groß annimmt. In Fällen, wo ha, d.h. die radiale Ausdehnung des Querschnitts im Verhältnis zum Wickelungshalbmesser groß ist, ergeben sich zwischen den Formeln der Hütte von de St. Vénant und des Verfassers ganz erhebliche Abweichungen, wie an einem Beispiel bewiesen wird. Es wird darauf hingewiesen, daß alle Formeln nur zuverlässig sein können, wenn die Voraussetzungen zufolge Abschnitt I angenommen werden dürfen. Es ist nicht unmöglich, daß z.B. bei den Formeln der Hütte die Vernachlässigung der Steigung einerseits und die der genaueren Spannungsverteilung anderseits zu Fehlern führen, die einander aufheben. Für die praktische Berechnung insbesondere steil gewickelter und hochbeanspruchter Federn ist daher die Prüfung der entwickelten Formeln durch den Versuch erforderlich.