Gelöste Aufgaben/Bike

Aus numpedia
Zur Navigation springen Zur Suche springen


Aufgabenstellung

Ein Hersteller von Scheibenbremsbelägen erhält wiederholt Reklamationen, weil Fahrräder beim Bremsen erhebliche Geräusche entwickeln oder die Gabel starke Schwingungen in der Laufradebene ausführt.

Vorderradgabel mit Scheibenbremse.

Gesucht ist ein mathematisches Modell, das die Entstehung von Geräuschen und Schwingungen an einer Vorderradgabel zu erklären hilft.

Wir simulieren dazu das dynamische Verhalten eines Vorderrades beim Bremsen mit Scheibenbremse.

Strukturieren

Wir suchen einen Mechanismus, bei dem sich eine Schwingung ohne äußere Erregung einstellt. In der Mechanik spricht man dann von einer selbsterregten Schwingung. Bei rein mechanischen Systemen wie unserem Vorderrad sind dabei oft Stick-Slip-Schwingungen die Ursache, bei denen z.B. die Bremse kurzzeitig auf der Bremsscheibe blockiert und sich dann wieder losreißt. Das charakteristische Merkmal, das diesen Mechanismus treibt, ist dabei eine fallende Kennlinie: bei fast allen Materialpaarungen ist die Haftkraft der Bremse größer als die Reibkraft.

Aufgabenstellung konkretisieren

Wir erwarten, den gesuchten Mechanismus nur anhand der Dynamik des Vorderrades beschreiben zu können. Das Vorderrad wird deshalb am Steuersatz vom Rahmen getrennt und dort – gedacht - in einer starren Führung mit konstanter horizontaler Geschwindigkeit v0 geführt. Um das Modell einfach zu halten, untersuchen wir Schwingungen nur in der Radebene. Das ist eine wesentliche Vereinfachung, weil wir bei Schwingungen unter Einwirkung der Bremse von Bewegungen des Rades um die Steuerrad-Achse rechnen müssen. Diese Vernachlässigen wir. Eine weitere Vereinfachung führen wir durch Setzen des Steuerwinkels von 90° ein. Damit vernachlässigen wir den Beitrag einer vertikalen Bewegung: beim Schwingen der Gabel in der Radebene bleibt der Steuersatz unverändert in der Ausgangshöhe. Diese Auf- und Ab-Bewegung berücksichtigen wir also nicht.

System strukturieren

Unser verbleibendes System besteht nun aus

  • Federgabel (Länge H) und Bremssattel (Abmessungen h und b),
  • Nabe (Außenradius r1)mit Bremsscheibe (Radius r2)
  • Felge (Innenradius R2)und Mantel (Außenradius R1) sowie den
  • Speichen.

Dieses System wollen wir als dynamisches System mit Bewegungen nur in der Radebene beschreiben. Wir machen folgende Annahmen für die Modellierung: Als starre Körper beschrieben wir die Teilkomponenten

  • Felge & Mantel sowie,
  • Nabe & Bremsscheibe.
  • Als rein elastischen Körper beschreiben wir die Speichen, die eine Verdrehung zwischen Nabe und Felge erlauben.

Als Kontinuum mit Masse und Elastizität beschrieben wir die Federgabel. Diese hat einen über ihre Länge veränderlichen Querschnitt.

Als Referenzpunkte wählen wir

  • A: den Aufstandspunkt des Rades auf dem Asphalt,
  • B: den Berührungspunkt zwischen Bremsscheibe und Bremsbelag,
  • C: die Achse des Rades, die Naben und Gabel verbindet sowie
  • N0,N1,N2: als Punkte entlang der x-Achse der Federgabel.

Um die Auslenkung der Federgabel zu erfassen, führen wir materielle Koordinaten x,z nach Vorgabe des Euler-Bernoulli-Balkens ein. Den Ursprung des Koordinatensystems legen wir passend in die Wurzel der Gabel bei N0.

Vorderradgabel mit Scheibenbremse.
Vorderradgabel mit Scheibenbremse.

Der Einfachheit halber sollen C und N2 zusammen fallen. Wir modellieren die Gabel mit zwei Finiten Elementen:

  1. von N0 nach N1 und
  2. von N1 nach N2.
Als Koordinaten der Drehbewegung führen wir
  • R1Ωt=v0 und relativ dazu
  • Ψ1(t) für Felge/Mantel und
  • Ψ2(t) für Nabe/Bremsscheibe

ein.

Freikörperbild Vorderrrad.
Freikörperbild Gabel.
Koordinaten der Gabel für N1.

Die Bewegung der Gabel erfassen wir zunächst an den Punkten N0,N1,N2 mit den Koordinaten der Bewegung Wi(t) (horizontale Auslenkung) und Φi(t) (Kippen des Querschnitts) für i=0,1,2.

Schnittlasten.
Querschnitt.


Q~_=(W0Φ0W1Φ1W2Φ2,Ψ1,Ψ2)
δQ~_=(δW0δΦ0δW1δΦ1δW2δΦ2,δΨ1,δΨ2)


Y_=(Q_Q˙_)
Y˙_=f_(Y_)=(Q˙_M__1K__Q_+M__1P_)
M__Q¨_+K__Q_=P_(Q_,Q˙_)
δW=δWaδΠ=!0
δΠ=δΠG+δΠS
δΠS=KS(Ψ2Ψ1)(δΨ2δΨ1)
δΠG=δΠG1+δΠG2
δΠGi=0iMi(x)δwi(x)dx
Mi(x)=EI(x)wi(x)
wi(x)=Q_iTφ_
φ_=((ξ1)2(2ξ+1)iξ(ξ1)2ξ2(2ξ3)iξ2(ξ1))
Q~_i=(Wi1ΦiqWiΦi)
δΠGi=0iEIi(x)(Wi1φ1δWi1φ1+Φi1φ2δWi1φ1+Wiφ3δWi1φ1++Φiφ4δΦiφ4)
(.):=d(.)dx
δΠGi=δQiTK__iQi
ki,jk=0iEIi(x)φjφkdx
Ii(x)=π64(Di(x)4di(x)4)


Di(x)=Di1(ξ11)+Diξ1di(x)=di1(ξ11)+diξ1
ξ1=x1
ξ2=(xH)2
d(.)dx=d(.)dξi1i
K__=(k1,33+k2,11k1,34+k2,12k2,13k2,1400k1,34+k2,12k1,44+k2,22k2,23k2,2400k2,13k2,23k2,33k2,3400k2,14k2,24k2,34k2,44000000KSKS0000KSKS)
A__=(0K__000010R100)
b_=(P_0)
x_=(Q_AH)


K~__=(k1,11k1,12k1,13k1,140000k1,12k1,22k1,23k1,240000k1,13k1,23k1,33+k2,11k1,34+k2,12k2,13k2,1400k1,14k1,24k1,34+k2,12k1,44+k2,22k2,23k2,240000k2,13k2,23k2,33k2,340000k2,14k2,24k2,34k2,4400000000KSKS000000KSKS)


M~__=(m1,11m1,12km,13m1,140000m1,12m1,22km,23m1,240000m1,13m1,23km,33+m2,11m1,34+m2,12m2,13m2,1400m1,14m1,24km,34+m2,12m1,44+m2,22m2,23m2,240000m2,13m2,23m2,33+Mf+Mnm2,340000m2,14m2,24m2,34m2,4400000000Jf00000000Jn)


r_B,1=(Hhbsin(Φ1(t))v0t+W1(t)+bcos(Φ1(t)))
r_B,2=(Hr2cos(Ωt+Ψ2(t)+Θ2)v0t+W2(t)r2sin(Ωt+Ψ2(t)+Θ2))
r_A=(HR1cos(Ωt+Ψ1(t)+Θ1)v0t+W1(t)R1sin(Ωt+Ψ1(t)+Θ1))


v_B,rel=r˙_B,1r˙_B,2
v_A,rel=r˙_A=(0v0+R1(Ω+Ψ˙1(t))+W˙2(t))


(.)˙:=d(.)dt


δWa=δWAa+δWBa+δWdAlemberta
Θ1=ΩtΨ1(t)+πΘ2=ΩtΨ1(t)+ψ


δr_B,1=(δΦ1bδW1)
δr_B2=(r2δΨ2sin(ψ)δW2r2δΨ2cos(ψ))
δr_A=(0δW2+R1δΨ1)
δWAa=F_ATδr_A
δWBa=F_BTδr_B,1F_BTδr_B,2
F_B=(BxBz)
F_A=(AVAH)
h=r2cos(ψ)b=r2sin(ψ)


δWAa+δWBa=BzδW1bBxδΦ1+δW2(˙AHBz)+δΨ1R1AH+δΨ2(r2Bzcos(ψ)r2Bxsin(ψ))


tmp

Title

Text


1+1






Links

  • ...

Literature

  • ...