Simulation > Creo Flow Analysis > Preprocessing > Physik definieren > Dynamik > Physik > Bewegungen eines Starrkörpers
Bewegungen eines Starrkörpers
In Simulationen sind die Flächen eines Volumenkörperobjekts normalerweise Wandberandungen in einer Strömungsdomäne. Wenn ein Volumenkörperobjekt oder eine Fläche dynamischen und mechanischen Kräften und einer thermischen Wirkung ausgesetzt ist, kann das Ungleichgewicht der Nettokräfte dazu führen, dass sich der Körper bewegt und verformt. Ein Volumenkörperobjekt wird in Fluss-Simulationen häufig als Starrkörper betrachtet. Daher wird für ein Volumenkörperobjekt, das Kraftungleichgewichten ausgesetzt ist, davon ausgegangen, dass es linear (Translation), radial (Rotation) oder sowohl linear als auch radial ohne Verformung bewegt werden kann. Bei einer CFA-Rechendomäne kann die Berandungsbewegung jedoch zur Domänenänderung führen, und daher kann sich das Volumennetz verformen, wie im Modul Flow beschrieben.
Für einen Starrkörper werden die Gleichungen, die für seine Bewegungen gelten, direkt aus der Erhaltung des linearen Impulses und des Drehimpulses abgeleitet:
Linearer Impuls (Translation)
Gleichung 2.426
Drehimpuls (Rotation)
Gleichung 2.427
In Gleichung 2.426 ist die Masse des sich bewegenden Objekts; ⃗ ist die lineare/translatorische Geschwindigkeit; und ⃗ sind die Gesamt-/Nettokräfte, die auf den Körper unter der Translation ausgeübt werden. In Gleichung 2.427 ist das Trägheitsmoment, ⃗ die Winkelgeschwindigkeit und ⃗ das Gesamt-/Netto-Drehmoment, das auf den rotierenden Körper wirkt.
Gleichung 2.426 und Gleichung 2.427 gelten für die allgemeinen Bewegungen eines Volumenkörpers, die sechs Freiheitsgrade (6-DOF) mit jeweils drei Graden für Translation (3-DOF) und Rotation (3-DOF) haben. Creo Flow Analysis berücksichtigt nur 1-FG-Translation und -Rotation, die in diesem Abschnitt erläutert werden.
1-FG-Translation
Unter der Annahme, dass sich ein Volumenkörper linear in eine beliebig angegebene Richtung bewegt (bleibt unverändert), definiert durch einen Einheitsvektor , wird die translatorische Bewegung des Körpers auf einen Freiheitsgrad (1-DOF) reduziert. Anschließend wird Gleichung 2.426 zur linearen Impulserhaltung zu einer skalaren Gleichung entlang der Bewegungsrichtung, da die Bewegungsgeschwindigkeit und -kraft in Bezug auf ausgedrückt werden:
Gleichung 2.428
Gleichung 2.429
Gleichung 2.430
Dabei ist der Betrag des Positionsvektors an einem Punkt von Interesse auf dem Volumenkörper entlang der Bewegungsrichtung . In einem kartesischen Koordinatensystem gilt
Gleichung 2.431
Wenn die Masse des Volumenkörpers konstant bleibt und der Kraftterm so erweitert wird, dass explizit alle auf den Körper angewendeten Kräfte eingeschlossen werden, ergibt sich als skalare lineare Impulsgleichung:
Gleichung 2.432
Die Kräfte auf der rechten Seite geben Folgendes an:
Strömungskraft – Besteht aus Druck- und Schubkräften. Sie werden durch die relative Bewegung zwischen dem Flüssigkeitsfluss und den Flächen des Volumenkörpers verursacht, die in Kontakt mit dem Fluss stehen. Die Druck- und Schubkräfte werden aus den Flow-Lösungen (Ausgabemengen) ermittelt:
Gleichung 2.433
Dämpfungskraft – Eine Verzögerungskraft, die durch die Reibungsdämpfungswirkung verursacht wird. Sie wird durch die Bewegung des Volumenkörperobjekts und den benutzerdefinierten Dämpfungskoeffizienten bestimmt:
Gleichung 2.434
Federkraft – Ist abhängig von der Federverschiebung , der Federkonstante und der Federvorspannkraft :
Gleichung 2.435
Dabei ist die Federverschiebung definiert als:
Gleichung 2.436
Dabei ist der Betrag des Positionsvektors an der vorherigen Position .
Reibungskraft – Das Kontaktreibungsmodell wird übernommen, um die Wirkung der Reibung in einem dynamischen System zu berücksichtigen. Die Reibungskraft wird modelliert als:
Gleichung 2.437
Dabei ist die Normalkomponente der Kontaktkraft, die auf die Körperfläche von Interesse ausgeübt wird. Für den Reibungskoeffizienten werden der Haftreibungskoeffizient und der Gleitreibungskoeffizient weiter für die stationären bzw. sich bewegenden Körper eingeführt:
Gleichung 2.438
Zusätzliche Kraft – Für zusätzliche benutzerdefinierte Kräfte hinzugefügt.
1-FG-Rotation
Wenn eine beliebige rotierende Achse durch einen Punkt (Mittelpunkt der Achse) und den Richtungseinheitsvektor definiert ist, wird die Rotation des Volumenkörpers um die Achse ebenfalls auf eine 1-FG-Rotation reduziert. Entsprechend wird Gleichung 2.427 zur Drehimpulserhaltung ebenfalls zu einer skalaren Gleichung entlang der tangentialen Richtung , definiert als:
Gleichung 2.439
Dabei ist der Vektor, der vom Mittelpunkt der Achse zu einem beliebigen Punkt auf dem Volumenkörper zeigt:
Gleichung 2.440
Die Winkelgeschwindigkeit und das Drehmoment am Punkt werden umformuliert als:
Gleichung 2.441
Gleichung 2.442
Gleichung 2.443
Dabei ist der Rotationswinkel des Punkts relativ zur Start- oder Referenzposition.
Wenn das Trägheitsmoment konstant bleibt und der Kraftterm so erweitert wird, dass explizit alle auf den Körper angewendeten Drehmomente eingeschlossen werden, ergibt sich als skalare Drehimpulsgleichung:
Gleichung 2.444
Die Drehmomentterme auf der rechten Seite sind wie folgt definiert:
Hydrodynamisches Drehmoment – Kombination des Drehmoments aufgrund von Druck- und Schubkräften:
Gleichung 2.445
Dämpfungsmoment – Hängt von der Rotationsgeschwindigkeit und dem benutzerdefinierten Dämpfungskoeffizienten ab:
Gleichung 2.446
Federdrehmoment – Drehmoment induziert durch Torsion, die vom Verschiebungswinkel , vom benutzerdefinierten Vorspannungsmoment und von der Torsionskonstante abhängt.
Gleichung 2.447
Dabei ist der Referenzwinkel. Dieser ist in der Regel die Position der Berandung oder des Volumens während des Modell-Setups, kann jedoch einer anderen Position entsprechen. Bei einer Nullwinkelverschiebung ist der Referenzwinkel beispielsweise nicht gleich der Anfangswinkelposition.
Reibungsdrehmoment – Drehmoment, das durch die Reibungskraft entsteht, die bei der Bewegung zweier Objekte in Kontakt auftritt. In Experimenten wird dieses durch die Differenz zwischen dem angewendeten Drehmoment und dem beobachteten oder Netto-Drehmoment bestimmt. Es ist abhängig vom Reibungskoeffizienten und dem Kontaktmoment aus der Anwendung der Normalkraft auf die Kontaktfläche:
Gleichung 2.448
Dabei ist ein benutzerdefinierter Parameter, der in Gleichung 2.438 definiert ist.
Zusätzliche Drehmomente – Für zusätzliche benutzerdefinierte Drehmomente hinzugefügt.
Ablenkungsmodell
In vielen Situationen wird ein Volumenkörper nur in einem begrenzten Raum (begrenzter Abstand oder Winkel) verschoben oder rotiert oder sowohl verschoben als auch rotiert, d.h. er hat eine maximale, minimale oder sowohl eine maximale als auch eine minimale Position. Beispiel: Wenn, wie in der folgenden Abbildung gezeigt, ein einfaches Schwerkraftpendel von der ursprünglichen Position mit dem Winkel freigegeben wird, verursacht die auf seine Masse wirkende Rückstellkraft, dass es um die Gleichgewichtsposition herum schwingt. Der maximale Winkel auf jeder der Seiten der Gleichgewichtsposition hängt von seiner Freigabeposition ab. Wenn keine Reibung (reibungsfreier Drehpunkt und Vakuum) vorhanden ist, bleibt der maximale Winkel unverändert, und das Pendel schwingt dauerhaft zwischen denselben Extrempositionen hin und her. Wenn sich ein Pendel dagegen beispielsweise in der Atmosphäre befindet, führt der Luftwiderstand (Dämpfung) dazu, dass der maximale Schwingwinkel mit der Zeit kleiner wird und das Pendel schließlich an der Gleichgewichtsposition anhält.
Abbildung
1. Reibungsloser Neutralpunkt
2. Masseloser Stab
3. Massives Pendelgewicht
4. Gleichgewichtsposition
5. Leitkurve des Pendelgewichts
6. Amplitude
Außerdem ändert in einem Schwingzyklus (Periode) das Pendel, wenn es die höchste Position erreicht, die Richtung mit dem Gesamtverlust seiner kinetischen Energie. Bei dem einfachen Schwerkraftpendel wird die kinetische Energie vollständig in die potentielle Energie übertragen, während dann, wenn ein Widerstand des Mediums berücksichtigt wird, ein Teil der kinetischen Energie verloren geht, um die viskose Dämpfung zu überwinden. Die Nettokraft oder die potentielle Energie treibt das Pendel jedoch dazu an, sich in die entgegengesetzte Richtung hin zur Gleichgewichtsposition zu bewegen, wo die kinetische Energie (Geschwindigkeit) maximal und das Potenzial am niedrigsten ist. In diesem Fall gibt eine No-Bounce-Bedingung für die Gleichung 2.444, den 1-FG-Drehimpuls, an.
Zusätzlich zur No-Bounce-Bedingung kann ein beweglicher Körper an der Grenzposition möglicherweise keinerlei kinetische Energie verlieren und zurückprallen (perfekte Ablenkung) oder nur einen Teil seiner kinetischen Energie verlieren (partielle Ablenkung). Daher werden die folgenden drei Ablenkungsbedingungen angewendet, wenn die Dynamikgleichungen für Translation (1 FG) und Rotation (1 FG), Gleichung 2.432 und Gleichung 2.444 gelöst werden, um die Bewegungen eines Volumenkörpers oder einer Wandberandung für die Strömungsdomäne zu ermitteln:
No Bounce – Standardmodell in Creo Flow Analysis. Dadurch wird bestimmt, dass ein Volumenkörper (oder eine Berandung) bei Erreichen seiner Bewegungsgrenze die Richtung ändert und dabei seine kinetische Energie vollständig verloren geht. Wenn und Ablenkung und Inzidenz und und die Translations- und Rotationsgeschwindigkeit (nur Betrag) darstellen, wird dieses Ablenkungsmodell wie folgt ausgedrückt:
Translation
Gleichung 2.449
Rotation
Gleichung 2.450
Partielle Ablenkung – Modell, das vorgibt, dass ein Volumenkörper (oder eine Berandung) bei Erreichen seiner Bewegungsgrenze die Richtung ändert und dabei seine kinetische Energie teilweise verloren geht, wie durch einen benutzerspezifizierten Faktor angegeben:
Translation
Gleichung 2.451
Rotation
Gleichung 2.452
Perfekte Ablenkung – Modell, das vorgibt, dass ein Volumenkörper (oder eine Berandung) bei Erreichen seiner Bewegungsgrenze die Richtung ändert und dabei seine kinetischen Energie in keiner Weise verloren geht, :
Translation
Gleichung 2.453
Rotation
Gleichung 2.454