Rechner für Deckenschwingungen durch Waschmaschinen und Isolierpads
Simulieren Sie Schleuderschwingungen einer Waschmaschine auf einer Holzbalkendecke. Vergleichen Sie Isolierpads, einen Masseblock, Balkenaufdopplung und Querhölzer anhand eines gekoppelten Maschine-Decke-Modells.
Deckenschwingungen durch eine Waschmaschine sind häufig während des Hochlaufs am stärksten, nicht bei maximaler Schleuderdrehzahl. Die rotierende Beladungsunwucht erzeugt eine Kraft, die mit dem Quadrat der Drehzahl zunimmt. Zugleich besitzen Maschine, Gerätefüße und Holzbalkendecke Eigenfrequenzen. Beim Durchfahren jeder der beiden gekoppelten Resonanzen können die Deckenbalken stark schwingen, selbst wenn der Kraftübertragungsgrad bei Enddrehzahl gering ist.
Dieser Rechner richtet sich an Ingenieure, Studierende und Hauseigentümer, die Maßnahmen zur Verringerung der Schwingungen eines Waschmaschinensockels auf einer Holzbalkendecke vergleichen möchten. Er integriert ein gekoppeltes Maschine-Decke-Modell mit zwei Freiheitsgraden über einen linearen Drehzahlhochlauf und eine anschließende Haltephase bei Enddrehzahl. Die Ergebnisse unterscheiden zwischen der höchsten transienten Deckenbeschleunigung und der stationären Antwort bei der gewählten Enddrehzahl.
Aufstellung eingeben
Beginnen Sie mit der Masse der beladenen Maschine, der wirksamen Restunwucht und dem Drehzahlprofil. Die Unwucht ist auf das Gehäuse bezogen: Eine elastisch aufgehängte Trommel kann wesentlich weniger Kraft übertragen als eine starr rotierende Masse. Leiten Sie diesen Eingabewert möglichst aus Messungen ab. Geben Sie Balkenabmessungen, lichte Spannweite, Balkenachsabstand, verteilte Deckenmasse und mitwirkende Breite ein. Die durch eine halbe Sinuswelle beschriebene Deckeneigenform bildet gelenkig gelagerte Balken mit gemeinsamer Lastabtragung ab; eine Verbundwirkung der Beplankung oder nachgiebige Unterzüge werden nicht automatisch berücksichtigt.
Verwenden Sie für die Berechnung mit Isolierpads die gemeinsame dynamische Steifigkeit aller eingebauten Gerätefüße und Pads sowie deren Dämpfungsmaß. Die statische Einfederung wird als Plausibilitätskontrolle angezeigt; die dynamische Steifigkeit von Gummi kann jedoch von seiner statischen Steifigkeit abweichen. Vergleichen Sie den verfügbaren Federweg um die belastete Gleichgewichtslage mit der berechneten Relativbewegung.
Maßnahmen vergleichen
Die Tabelle vergleicht Originalfüße, Isolierpads, Pads mit zusätzlichem Masseblock, Pads mit idealer seitlicher Balkenaufdopplung über die volle Spannweite sowie Pads mit erhöhter Lastverteilung durch Querhölzer. Wählen Sie die jeweilige Konfiguration aus, um ihren transienten Zeitverlauf und ihre Animation zu betrachten. Zeilen mit Kennzeichnung für überschrittenen Pad-Federweg oder Abhebegefahr werden bei der Ermittlung des niedrigsten gültigen Spitzenwerts ausgeschlossen. Ohne Unwuchtanregung sind normierte Spitzenwertverhältnisse nicht verfügbar; alle Varianten sind dann gleichwertig. Ein Masseblock oberhalb der Pads bewegt sich mit der Maschine; ein Masseblock unterhalb der Pads erhöht die modale Deckenmasse. Die mitwirkende Breite mit Querhölzern ist eine ausdrücklich zu überprüfende Annahme. Sie ist keine zugesicherte Wirkung einer bestimmten Anzahl eingebauter Querhölzer.
Die stationären Antwortkurven über der Drehzahl helfen, Resonanzen und zu vermeidende Betriebsdrehzahlen zu erkennen. Zum Vergleich reichen sie unabhängig vom simulierten Profil bis mindestens 1.500 rpm. Das Hochlaufdiagramm zeigt den unbewerteten gleitenden Effektivwert der Beschleunigung über die jeweils letzte Sekunde; die Komfortbewertung verwendet den tatsächlichen momentanen Spitzenwert. Diese Größen werden bewusst getrennt betrachtet.
Überschlägige Bewertung einordnen
Der voreingestellte Grenzwert der Spitzenbeschleunigung beträgt 0,5 % der Erdbeschleunigung und entspricht einem AISC-Richtwert für ruhige Aufenthaltsbereiche. Er dient bei diesem Maschinenschwingungsproblem als Vergleichswert, nicht als Bewertung nach ISO 10137 oder als Garantie für die Akzeptanz durch die Nutzer. Dauer, Frequenzbewertung, Geräusche und benachbarte Räume können die wahrgenommene Störung beeinflussen.
Das Modell berücksichtigt weder Kippschwingungen, Wandern, lose Gerätefüße, Abheben, nichtlineare Pads noch die interne Trommelaufhängung oder höhere Deckeneigenformen. Montagearbeiten und bauliche Änderungen sind von einer qualifizierten Fachperson auszuführen. Für einen schwereren Sockel sind Tragfähigkeit der Decke, Auflager und Verbindungen gesondert nachzuweisen.
So wird gerechnet
Die Gleichungen werden hier aus den Newtonschen Gesetzen und der Balkenenergie nach Euler–Bernoulli hergeleitet. Die klassische Interpretation der Schwingungsisolierung folgt Den Hartog, Mechanical Vibrations, 4th edition (1956; Dover reprint 1985), sowie den frei zugänglichen Purdue-Unterlagen zur Schwingungsisolierung.
Dabei bezeichnet die Maschinenverschiebung und die Deckenverschiebung an den Gerätefüßen, jeweils bezogen auf die statische Gleichgewichtslage. Für die Lagersteifigkeit und den Dämpfungskoeffizienten gilt:
Die Unwucht beträgt . Im Hochlauf gelten und ; während der Haltephase gilt . Die vertikale Erregerkraft lautet und enthält damit auch den tangentialen Anteil während des Hochlaufs.
Für mitwirkende Balken, , und die längenbezogene Masse wird angesetzt:
Die Beschleunigung am Beobachtungspunkt wird mit skaliert. Der gesonderte statische Vergleich verwendet die exakte Balkensteifigkeit in Feldmitte anstelle der Näherung über die erste Eigenform. Der Canadian Guide to Mid-Rise Wood Construction des CWC, S. 130 der beigefügten Fassung, erläutert die Bedeutung der gesamten Deckenkonstruktion und der Durchbiegung unter Einzellast; dieses Werkzeug setzt das vollständige CSA/NBC-Verfahren für schwingungsbegrenzte Spannweiten nicht um.
Die fünf Konfigurationen werden mit dem impliziten Newmark-Verfahren der mittleren Beschleunigung (, ) integriert. Die zugeführte Arbeit wird aus der mittleren Kraft an den Intervallenden und dem Verschiebungsinkrement berechnet; die Dämpfungsarbeit verwendet die Geschwindigkeit in der Intervallmitte. Daraus ergibt sich eine unabhängige diskrete Energiebilanz. Der Solver verwendet ausdrücklich normierte SI-Koordinaten und stellt an seinen Ein- und Ausgängen die physikalischen Einheiten wieder her.
Die exakte komplexe dynamische Steifigkeit des Zweimassensystems liefert die stationären Antworten bei der Drehzahlvariation und die übertragene Lagerkraft. Die gekoppelten ungedämpften Eigenfrequenzen ergeben sich aus . Für den Vergleich bei starrer Decke gilt ; die Isolierwirkung beginnt oberhalb von .
Der Erregerkraftverlauf wird einmal auf dem gemeinsamen Integrationszeitraster berechnet und für alle fünf Konfigurationen wiederverwendet. Bei exakt verschwindender Unwucht wird der analytische Zeitverlauf ohne Bewegung angesetzt. Die Hüllkurve der Lagerkraft wird aus dem äquivalenten Maschinengleichgewicht bestimmt. Dadurch werden wiederholte Rechenoperationen vermieden, ohne Zeitschritt oder Gleichungen zu ändern.
Annahmen
- Eine vertikale Gehäusekoordinate und eine Biegeeigenform der Decke; starre Auflager, identische gelenkig gelagerte Balken, lineare kleine Bewegungen und eine anfangs ruhende statische Gleichgewichtslage.
- Die Masse der beladenen Maschine enthält die wirksame rotierende Masse. Unwucht, Steifigkeit, Dämpfung, Deckenmasse und mitwirkende Breite sind anpassbare Szenarioannahmen, keine Kennwerte eines bestimmten Geräts oder einer bestimmten Holzfestigkeitsklasse.
- Bei der seitlichen Balkenaufdopplung wird jedem mitwirkenden Balken ein identischer, über die volle Spannweite vollständig mitwirkender Beibalken hinzugefügt. Biegesteifigkeit und Masse der Balken verdoppeln sich; die Masse der Beplankung bleibt unverändert. Beibalken über Teilspannweiten und Verbindungsschlupf erfordern ein gesondertes Modell.
- Querhölzer verändern die angenommene mitwirkende Breite und erhöhen die Masse um die eingegebene äquivalente modale Masse. Ein zusätzlicher Multiplikator für die Biegesteifigkeit der Balken wird ihnen nicht zugeschrieben. Der optionale Masseblock ist starr mit der gewählten Seite der Pad-Lagerung verbunden.
- Jede Konfiguration verwendet das eingegebene Dämpfungsmaß mit für Lagerung und Decke. Eine Änderung der getragenen Masse bei konstantem Dämpfungsmaß verändert daher auch den äquivalenten viskosen Dämpfungskoeffizienten. Die tatsächlichen Pad-Eigenschaften sollten bei dieser Belastung gemessen werden.
- Aus der Durchbiegung unter 1 kN wird keine Bewertung der Akzeptanz von gehinduzierten Schwingungen nach CWC/Onysko abgeleitet. Das Einzellastergebnis dient zur Beurteilung der Nachgiebigkeit der Decke.
Anwendungsbereich
Eine überschlägige Komfortbewertung für eine vorgegebene wirksame Unwucht, keine Vorhersage für ein bestimmtes Gerät und kein Tragfähigkeitsnachweis. AISC Design Guide 11, 2nd edition (2016, with subsequent errata), behandelt menschliche Aktivitäten auf Stahltragwerken. Der voreingestellte Zielwert der Spitzenbeschleunigung von 0,005 g stammt aus der öffentlich zugänglichen AISC-Veröffentlichung Facts for Steel Buildings 5: Vibration (2018), Tabelle 2-1; seine Anwendung auf Waschmaschinenschwingungen ist ein ausdrücklich ingenieurmäßiger Vergleich. Wiederkehrende Maschinenschwingungen können einen anderen Grenzwert erfordern. Es erfolgen weder eine ISO-Frequenzbewertung noch eine Bewertung der Schwingungsdosis oder Exposition; prüfen Sie die Anforderungen der jeweils gültigen Norm.
Das Deckenmodell berücksichtigt weder Plattenwirkung, Verbundsteifigkeit der Beplankung, höhere Biege- und Torsionseigenformen, nachgiebige Auflager noch benachbarte Deckenfelder. Das Maschinenmodell berücksichtigt weder die interne Bottichdynamik, Kippschwingungen, horizontale Bewegungen, Gleiten noch Abheben. Die Pads bleiben linear und in Kontakt; passen Sie das Modell an, wenn Bewegungs- oder Kontaktbedingungen eine Warnung auslösen. Mitwirkende Breite, Dämpfung und die auf das Gehäuse bezogene Unwucht bestimmen häufig den größten Teil der Unsicherheit. Bestätigen Sie diese Größen durch Messungen vor Ort, bevor Sie sich auf den Vergleich von Nachrüstmaßnahmen stützen.
Der maximale Zeitschritt wird so reduziert, dass mindestens 32 Abtastwerte je Periode der schnellsten gekoppelten Eigenform oder je Schleuderumdrehung bei Enddrehzahl vorliegen. Die Voreinstellung wird mit halbiertem benutzerseitigem Zeitschritt überprüft. Die Newmark-Zeitdiskretisierung kann scharfe Resonanzen verschieben; verringern Sie den Zeitschritt unter „Erweitert“ bei geringer Dämpfung oder Haltephasen nahe einer Resonanz. Profile mit mehr als 30.000 Intervallen liegen außerhalb des numerischen Bereichs dieses Werkzeugs und werden gesperrt. Die Drehzahlvariation umfasst 1.501 berechnete Drehzahlen; sehr schmale Resonanzspitzen können dazwischen liegen. Die rückblickende Effektivwertkurve verwendet Trapezintegration und Interpolation des gespeicherten Zeitverlaufs; sie ist unbewertet und wird nicht mit dem Komfortgrenzwert für Spitzenwerte verglichen.
Die Animation überhöht die dynamischen Verschiebungen mit einem gemeinsamen automatischen Maßstab und lässt die statische Einfederung weg. Sie veranschaulicht die berechnete Schwingungsform, nicht die tatsächlich sichtbare Bewegung. Bei der statischen Einfederung wird angenommen, dass die eingegebene dynamische Steifigkeit auch der statischen Steifigkeit entspricht; prüfen Sie dies anhand der Last-Verformungs-Kennlinie der Pads.
Jeder Vergleich prüft den verfügbaren Federweg der vier Pad-Konfigurationen. Zur Kontaktprüfung wird konservativ der höchste Betrag der dynamischen Lagerkraft mit der statisch getragenen Gewichtskraft verglichen. Für die Originalfüße gilt der eingegebene Pad-Federweg nicht. Gekennzeichnete Zeilen bleiben als lineare Extrapolationen sichtbar, werden aber bei der Ermittlung des niedrigsten gültigen Spitzenwerts ausgeschlossen; auch die angezeigte Konfiguration erhält außerhalb dieser Grenzen keine Komfortbewertung. Bleibt keine zulässige Zeile übrig, wird keine bevorzugte Variante ausgegeben. Ohne Anregung sind alle Konfigurationen gleichwertig; normierte Spitzenwertverhältnisse sind dann nicht definiert und werden als nicht verfügbar angezeigt.
Quellen
- Den Hartog, J. P., Mechanical Vibrations, 4th ed. (1956), Dover reprint (1985), ISBN 9780486647852. Classical vibration isolation background; open derivation verified against Purdue's notes: https://www.purdue.edu/freeform/ervibrations/chapter-v-animations/vibration-isolation/
- Purdue University / Embry-Riddle, Vibration isolation: frequency-ratio crossover and damping trends. https://www.purdue.edu/freeform/ervibrations/chapter-v-animations/vibration-isolation/
- AISC, Facts for Steel Buildings No. 5: Vibration (2018), Table 2-1 and quiet-occupancy discussion; companion to Design Guide 11, 2nd ed. (2016). https://www.aisc.org/media/a24dtcdo/facts-for-steel-buildings-5-vibrations.pdf
- Canadian Wood Council, Canadian Guide to Mid-Rise Wood Construction, section 5.3.5.2, printed p. 128 (attached p. 130): floor construction and vibration serviceability. https://cwc.ca/wp-content/uploads/2025/02/CDN_GuidetoMid-Rise_Wood_Construction.pdf
- ANSYS Mechanical APDL Verification Manual, VM50: Thomson simply-supported-beam example, target 28.766 Hz and BEAM189 result 28.761 Hz. https://ansyshelp.ansys.com/public/Views/Secured/corp/v252/en/ans_vm/Hlp_V_VM50.html
Eingaben
| Eingabe | Symbol | Einheit | Vorgabe |
|---|---|---|---|
| Masse der beladenen Maschine | M_w | kg | 80 |
| Wirksame Unwuchtmasse Auf das Gehäuse bezogene Restunwucht; verwenden Sie nach Möglichkeit die gemessene übertragene Erregerkraft. | m_u | kg | 0.15 |
| Unwuchtexzentrizität | e_u | mm | 20 |
| Enddrehzahl | n_top | rpm | 1200 |
| Hochlaufdauer | t_r | s | 12 |
| Haltezeit bei Enddrehzahl | t_h | s | 4 |
| Angezeigte Konfiguration | remedy | Isolierpads | |
| Gesamtsteifigkeit der Originalfüße | k_feet | kN/m | 800 |
| Dämpfungsmaß der Originalfüße | z_feet | 0.06 | |
| Gesamtsteifigkeit der Pad-Lagerung Gemeinsame dynamische Steifigkeit sämtlicher Gerätefüße und Pads bei der vorhandenen Belastung. | k_pads | kN/m | 80 |
| Dämpfungsmaß der Pads | z_pads | 0.15 | |
| Verfügbarer dynamischer Pad-Federweg Verbleibender Federweg um die belastete Gleichgewichtslage, bevor die Pads auf Block gehen oder den linearen Bereich verlassen. | d_travel | mm | 8 |
| Zusätzlicher Masseblock | M_b | kg | 80 |
| Position des Masseblocks | base_side | Unter den Pads, auf der Decke | |
| Lichte Balkenspannweite | L_j | m | 4 |
| Balkenbreite | b_j | mm | 45 |
| Balkenhöhe | h_j | mm | 195 |
| Balkenachsabstand | s_j | mm | 400 |
| Mitwirkende Breite In der ersten Deckeneigenform mitwirkende Breite; anhand von Versuchen oder einer Tragwerksberechnung kalibrieren. | W_f | m | 1.2 |
| Biege-Elastizitätsmodul der Balken | E_j | GPa | 11 |
| Rohdichte der Balken | rho_j | kg/m^3 | 450 |
| Flächenmasse von Beplankung und ständiger Einrichtung Verteilte Masse ohne Balken, Maschine und optionalen Masseblock. | q_f | kg/m^2 | 35 |
| Dämpfungsmaß der Decke | z_floor | 0.03 | |
| Maschinenposition / Spannweite Von einem Auflager aus gemessen; die Eigenform ist eine halbe Sinuswelle. | r_m | 0.5 | |
| Beobachtungsposition / Spannweite | r_o | 0.5 | |
| Mitwirkende Breite mit Querhölzern Angenommene mitwirkende Breite nach Einbau der Querhölzer; diese verdoppeln die Biegesteifigkeit der Balken nicht. | W_block | m | 2.4 |
| Zusätzliche modale Masse der Querhölzer Auf die Verschiebung am Maschinenstandort bezogene Masse. | M_block | kg | 5 |
| Komfortgrenzwert (Spitzenwert) Voreinstellung: 0,5 % g, ein AISC-Richtwert für ruhige Aufenthaltsbereiche, kein normatives Abnahmekriterium für Maschinen. | a_limit | m/s^2 | 0.04903325 |
| Maximaler Integrationszeitschritt Wird automatisch so reduziert, dass die schnellste modellierte Schwingung mit mindestens 32 Werten je Periode abgetastet wird. | dt_max | s | 0.0015 |
| Wiedergabezeitpunkt | t_view | s | 8.46 |
Ergebnisse
Mit den obigen Vorgabewerten liefert der Rechner:
| Ergebnis | Wert |
|---|---|
| Überschlägige Komfortbewertung | Über dem gewählten Komfortgrenzwert |
| Spitzenbeschleunigung der Decke | 0,251 m/s² |
| Spitzenwert / gewählter Grenzwert | 5,12 |
| Eigenfrequenz der unbelasteten Decke | 12,82 Hz |
| Untere gekoppelte Eigenfrequenz | 4,73 Hz |
| Obere gekoppelte Eigenfrequenz | 13,64 Hz |
| Schleuderdrehzahl bei Beschleunigungsspitze | 846 1/min |
| Maximaler Federweg der Lagerung | 0,13 mm |
| Kraftübertragungsgrad bei Enddrehzahl | 0,112 |
| Stationäre Beschleunigung bei Enddrehzahl | 0,083 m/s² |
| Statische Einfederung der Lagerung | 9,81 mm |
| Deckendurchbiegung unter 1 kN in Feldmitte | 1,45 mm |
| Variante mit niedrigstem gültigem Spitzenwert | Pads + Masseblock |
Mit Quellen abgeglichen
Jeder Fall unten wurde mit diesem Rechner berechnet und stimmt mit dem Wert aus seiner Quelle überein (innerhalb der Toleranz, wo eine angegeben ist).
Unabhängige Balkeneigenfrequenz nach ANSYS VM50
Quelle: ANSYS VM50, Thomson Vibration Theory and Applications (1965), p. 18, Ex. 1.5-1: 80 in square 2 in beam, E=30×10⁶ psi, w=1.124 lbf/in, g=386 in/s². Density uses the source's g, not standard gravity. Published target 28.766 Hz; BEAM189 FEM gives 28.761 Hz. The reference checks the unloaded floor beam; short ramp is irrelevant to that quantity. https://ansyshelp.ansys.com/public/Views/Secured/corp/v252/en/ans_vm/Hlp_V_VM50.html
Eingaben: Lichte Balkenspannweite = 2,032 m; Balkenbreite = 50,8 mm; Balkenhöhe = 50,8 mm; Biege-Elastizitätsmodul der Balken = 206,8 GPa; Rohdichte der Balken = 7.780 kg/m³; Flächenmasse von Beplankung und ständiger Einrichtung = 0 kg/m²; Mitwirkende Breite = 0,4 m; Hochlaufdauer = 3 s; Haltezeit bei Enddrehzahl = 1 s; die übrigen wie oben.
| Ergebnis | Quelle | Dieser Rechner |
|---|---|---|
| Eigenfrequenz der unbelasteten Decke | 28,766 Hz | 28,766 Hz |
Veröffentlichte Isolationsgrenze bei √2
Quelle: Purdue University / Embry-Riddle, Vibration isolation: at frequency ratio r=√2, transmissibility is unity regardless of damping. The rigid-floor base-motion and force-transmissibility expressions are identical. The speed is selected to impose that published ratio; unity is the independent published target. https://www.purdue.edu/freeform/ervibrations/chapter-v-animations/vibration-isolation/
Eingaben: Enddrehzahl = 427,1 rpm; Hochlaufdauer = 3 s; Haltezeit bei Enddrehzahl = 1 s; die übrigen wie oben.
| Ergebnis | Quelle | Dieser Rechner |
|---|---|---|
| Kraftübertragungsgrad bei starrer Decke | 1 | 1 |
Related Resources
Related articles
-
Numerical Investigation of Nonlinear Vibrations in a 2-DOF Mechanical System Using MATLAB -
Advanced Dragon Flight Dynamics, MATLAB-Based 6-DOF Simulation and Performance Evaluation -
Nonlinear Multi-Machine Transient Stability Assessment Using Classical Swing Dynamics and Fault-Dependent Network Topologies in Matlab

Responses