Zwangsbedingung: Unterschied zwischen den Versionen

Zwangsbedingung: Unterschied zwischen den Versionen

imported>Ajv39
K (Tippfehler entfernt)
 
imported>Jooonsen
K (Präzisierung vorherige Änderung + Orthographie)
 
Zeile 1: Zeile 1:
Als '''Zwangsbedingung''' wird in der [[Analytische Mechanik|analytischen Mechanik]] eine Einschränkung der Bewegungsfreiheit eines Ein- oder [[Mehrkörpersystem]]s bezeichnet, anders gesagt eine Bewegungsbeschränkung. Dadurch nimmt die Anzahl der [[Freiheitsgrad]]e eines Systems ab. Werden zu viele Zwangsbedingungen gestellt, so kann es passieren, dass keine physikalische Lösung existiert.
Als '''Zwangsbedingung''' wird in der [[Analytische Mechanik|analytischen Mechanik]] eine Einschränkung der Bewegungsfreiheit eines Ein- oder [[Mehrkörpersystem]]s bezeichnet, anders gesagt eine Bewegungsbeschränkung. Dadurch nimmt die Anzahl der [[Freiheitsgrad]]e eines Systems ab. Werden zu viele Zwangsbedingungen gestellt, so kann es passieren, dass keine mathematische Lösung existiert. Dann kann das Problem gegebenenfalls auch physikalisch nicht lösbar sein, sodass es beispielsweise zum Defekt eines Objektes kommt, oder aber die physikalische Lösbarkeit durch das Entstehen mechanischer Spannungen im Objekt gegeben ist.


Systeme mit Zwangsbedingungen können besonders gut beschreiben werden durch
Systeme mit Zwangsbedingungen können besonders gut beschrieben werden durch
* die [[Lagrange-Formalismus|Lagrangesche]] Formulierung der [[klassische Mechanik|klassischen Mechanik]]  
* die [[Lagrange-Formalismus|Lagrangesche]] Formulierung der [[klassische Mechanik|klassischen Mechanik]]  
* die [[Hamiltonsche Mechanik|Hamiltonsche]] Formulierung der klassischen Mechanik
* die [[Hamiltonsche Mechanik|Hamiltonsche]] Formulierung der klassischen Mechanik
Zeile 8: Zeile 8:


== Unterscheidung ==
== Unterscheidung ==
=== bezüglich Integrabilität ===
=== Bezüglich Integrabilität ===
Im Folgenden wird stets ein <math>N</math>-Teilchensystem in 3&nbsp;[[Raumdimension]]en betrachtet. Ohne Zwangsbedingungen bräuchte man für den [[Ortsvektor]] jedes Teilchens 3&nbsp;Raumkoordinaten, somit insgesamt <math>3N</math> Raumkoordinaten, um das gesamte System zu beschreiben. Diese Koordinaten werden fortlaufend durchnummeriert:
Im Folgenden wird stets ein <math>N</math>-Teilchensystem in 3&nbsp;[[Raumdimension]]en betrachtet. Ohne Zwangsbedingungen bräuchte man für den [[Ortsvektor]] jedes Teilchens 3&nbsp;Raumkoordinaten, somit insgesamt <math>3N</math> Raumkoordinaten, um das gesamte System zu beschreiben. Diese Koordinaten werden fortlaufend durchnummeriert:


Zeile 38: Zeile 38:


Denn notwendig für die Integrabilität ist, dass die Koeffizientenfunktionen <math>a_{li}</math> folgende [[Integrabilitätsbedingung]] erfüllen:
Denn notwendig für die Integrabilität ist, dass die Koeffizientenfunktionen <math>a_{li}</math> folgende [[Integrabilitätsbedingung]] erfüllen:
 
 
:<math>\begin{align}
:<math>\begin{align}
\frac{\partial a_{li}}{\partial x_{k}} & = \frac{\partial a_{lk}}{\partial x_{i}}\\
\frac{\partial a_{li}}{\partial x_{k}} & = \frac{\partial a_{lk}}{\partial x_{i}}\\
Zeile 56: Zeile 56:
Nicht-holonome oder auch anholonome Zwangsbedingungen können nicht als Gleichungen zwischen den Koordinaten formuliert werden. Die generalisierten Koordinaten, die in solchen anholonomen Zwangsbedingungen erscheinen, sind i.&nbsp;A. nicht unabhängig voneinander variierbar.
Nicht-holonome oder auch anholonome Zwangsbedingungen können nicht als Gleichungen zwischen den Koordinaten formuliert werden. Die generalisierten Koordinaten, die in solchen anholonomen Zwangsbedingungen erscheinen, sind i.&nbsp;A. nicht unabhängig voneinander variierbar.


Es handelt sich z.B. um [[Ungleichung]]en, wie Beschränkungen auf einen bestimmten Raumbereich:
Es handelt sich z.&nbsp;B. um [[Ungleichung]]en, wie Beschränkungen auf einen bestimmten Raumbereich:


:<math> \,f(q_1, q_2, \ldots, q_n, t) > 0</math>
:<math> \,f(q_1, q_2, \ldots, q_n, t) > 0</math>
Zeile 64: Zeile 64:
:<math>\sum_{i}a_{li} \cdot \mathrm{d}q_{i} + a_{lt} \cdot \mathrm{d}t = 0\ ,\quad l = 1, \ldots ,r</math>
:<math>\sum_{i}a_{li} \cdot \mathrm{d}q_{i} + a_{lt} \cdot \mathrm{d}t = 0\ ,\quad l = 1, \ldots ,r</math>


Nicht-integrabel heißt dabei, dass die Gleichung - anders als bei holonomen Zwangsbedingungen - nicht als vollständiges Differential einer Funktion darstellbar ist. Somit wird hier die Integrabilitätsbedingung von den Koeffizientenfunktionen nicht erfüllt:
Nicht-integrabel heißt dabei, dass die Gleichung anders als bei holonomen Zwangsbedingungen nicht als vollständiges Differential einer Funktion darstellbar ist. Somit wird hier die Integrabilitätsbedingung von den Koeffizientenfunktionen nicht erfüllt:


:<math>\frac{\partial a_{li}}{\partial q_k} \neq \frac{\partial a_{lk}}{\partial q_i} \quad\quad i,k \in \{ 1, \ldots ,n \}</math>
:<math>\frac{\partial a_{li}}{\partial q_k} \neq \frac{\partial a_{lk}}{\partial q_i} \quad\quad i,k \in \{ 1, \ldots ,n \}</math>


=== bezüglich Zeitabhängigkeit ===
=== Bezüglich Zeitabhängigkeit ===
Weiterhin werden Zwangsbedingungen bez. ihrer [[Zeit]]<nowiki/>abhängigkeit unterschieden in:
Weiterhin werden Zwangsbedingungen bez. ihrer [[Zeit]]<nowiki/>abhängigkeit unterschieden in:
* rheonom (fließend), wenn sie explizit von der Zeit abhängen.
* rheonom (fließend), wenn sie explizit von der Zeit abhängen.
Zeile 81: Zeile 81:
[[Datei:Pendel PT.svg|mini|Die Länge des Fadens bleibt konstant, das Pendel wird auf eine Kreisbahn gezwungen]]
[[Datei:Pendel PT.svg|mini|Die Länge des Fadens bleibt konstant, das Pendel wird auf eine Kreisbahn gezwungen]]
=== Das Pendel: holonom und skleronom ===
=== Das Pendel: holonom und skleronom ===
Der Stab eines ebenen [[Mathematisches Pendel|Pendel]]s (d.&nbsp;h. nur 2&nbsp;Raumdimensionen) soll stets die gleiche Länge <math>l</math> besitzen, muss also aufgrund des [[Satz des Pythagoras|Satzes von Pythagoras]] folgende Zwangsbedingung erfüllen (Anzahl der Zwangsbedingungen: <math>s = 1</math>):
Der Stab eines ebenen [[Mathematisches Pendel|Pendels]] (d.&nbsp;h. nur 2&nbsp;Raumdimensionen) soll stets die gleiche Länge <math>l</math> besitzen, muss also aufgrund des [[Satz des Pythagoras|Satzes von Pythagoras]] folgende Zwangsbedingung erfüllen (Anzahl der Zwangsbedingungen: <math>s = 1</math>):


:<math>f_1(x_1, x_2) = 0 \quad \mathrm{mit} \; x_1 = x, x_2 = y</math>
:<math>f_1(x_1, x_2) = 0 \quad \mathrm{mit} \; x_1 = x, x_2 = y</math>
:<math>\Leftrightarrow x^2 + y^2 - l^2 = 0</math>
:<math>\Leftrightarrow x^2 + y^2 - l^2 = 0</math>


Dabei bildet der Auslenkungswinkel <math>\varphi</math> des Pendels aus der Senkrechten, von dem die Koordinaten&nbsp;x und&nbsp;y des Kugelmittelpunktes abhängen (Anzahl der betrachteten Massepunkte: <math>N = 1</math>; Annahmen: x nach rechts, y nach unten, Ursprung im Aufhängungspunkt):
Dabei bildet der Auslenkungswinkel <math>\varphi</math> des Pendels aus der Senkrechten die generalisierte Koordinate. (Es gibt nur eine, da <math>n = 2N - s = 1</math>.) Die Koordinaten <math>x</math> und <math>y</math> des Kugelmittelpunktes hängen von <math>\varphi</math> ab (Annahmen: <math>x</math> nach rechts, <math>y</math> nach unten, Ursprung im Aufhängungspunkt):


::<math>\begin{align}
::<math>\begin{align}
x = l \cdot sin \, \varphi\\
x = l \cdot \sin\varphi\\
y = l \cdot cos \, \varphi
y = l \cdot \cos\varphi
\end{align}</math>
\end{align}</math>


die generalisierte Koordinate (Anzahl der generalisierten Koordinaten: <math>n = 2N - s = 1</math>).
Die generalisierte Koordinate erfüllt automatisch die Zwangsbedingung:
 
Die generalisierte Koordinate <math>\varphi</math> erfüllt automatisch die Zwangsbedingung:


:<math>\begin{align}
:<math>\begin{align}
(l \cdot sin \, \varphi)^2 + (l \cdot cos \, \varphi)^2 - l^2 = 0\\
(l \cdot \sin\varphi)^2 + (l \cdot \cos\varphi)^2 - l^2 = 0\\
\Leftrightarrow l^2 \cdot (sin^2 \, \varphi + cos^2 \, \varphi - 1) = 0
\Leftrightarrow l^2 \cdot (\sin^2\varphi + \cos^2\varphi - 1) = 0
\end{align}</math>
\end{align}</math>


da allgemein gilt:
da allgemein gilt:


::<math>sin^2 \, \varphi + cos^2 \, \varphi = 1</math>
::<math>\sin^2\varphi + \cos^2\varphi = 1</math>


Dies ist ein Beispiel für eine holonome Zwangsbedingung und, da sie nicht explizit von der Zeit abhängt (<math>f_1\neq f_1(t)</math>), für eine skleronome Zwangsbedingung.
Dies ist ein Beispiel für eine holonome Zwangsbedingung und, da sie nicht explizit von der Zeit abhängt (<math>f_1\neq f_1(t)</math>), für eine skleronome Zwangsbedingung.
Zeile 119: Zeile 117:


::<math>\begin{align}
::<math>\begin{align}
\frac{\mathrm{d}x}{\mathrm{d}t} & =  v \cdot cos \, \varphi\\
\frac{\mathrm{d}x}{\mathrm{d}t} & =  v \cdot \cos\varphi\\
\frac{\mathrm{d}y}{\mathrm{d}t} & = - v \cdot sin \, \varphi
\frac{\mathrm{d}y}{\mathrm{d}t} & = - v \cdot \sin\varphi
\end{align}</math>
\end{align}</math>


Zeile 127: Zeile 125:
Einsetzen der generalisierten Koordinate in die Zwangsbedingung in Form des vollständigen Differentials:
Einsetzen der generalisierten Koordinate in die Zwangsbedingung in Form des vollständigen Differentials:


:<math> \Rightarrow (l \cdot sin \, \varphi) \cdot (v \cdot cos \, \varphi) + (l \cdot cos \, \varphi) \cdot (- v \cdot sin \, \varphi)= 0</math>
:<math> \Rightarrow (l \cdot \sin\varphi) \cdot (v \cdot \cos\varphi) + (l \cdot \cos\varphi) \cdot (- v \cdot \sin\varphi)= 0</math>


das somit ebenfalls automatisch erfüllt ist.
das somit ebenfalls automatisch erfüllt ist.
Zeile 145: Zeile 143:
== Siehe auch ==
== Siehe auch ==
* [[Zwangskraft]]
* [[Zwangskraft]]
* [[Zwanglauf]]


[[Kategorie:Klassische Mechanik]]
[[Kategorie:Klassische Mechanik]]


[[en:Constraint (classical mechanics)]]
[[en:Constraint (classical mechanics)]]

Aktuelle Version vom 22. Juni 2021, 08:30 Uhr

Als Zwangsbedingung wird in der analytischen Mechanik eine Einschränkung der Bewegungsfreiheit eines Ein- oder Mehrkörpersystems bezeichnet, anders gesagt eine Bewegungsbeschränkung. Dadurch nimmt die Anzahl der Freiheitsgrade eines Systems ab. Werden zu viele Zwangsbedingungen gestellt, so kann es passieren, dass keine mathematische Lösung existiert. Dann kann das Problem gegebenenfalls auch physikalisch nicht lösbar sein, sodass es beispielsweise zum Defekt eines Objektes kommt, oder aber die physikalische Lösbarkeit durch das Entstehen mechanischer Spannungen im Objekt gegeben ist.

Systeme mit Zwangsbedingungen können besonders gut beschrieben werden durch

Unterscheidung

Bezüglich Integrabilität

Im Folgenden wird stets ein $ N $-Teilchensystem in 3 Raumdimensionen betrachtet. Ohne Zwangsbedingungen bräuchte man für den Ortsvektor jedes Teilchens 3 Raumkoordinaten, somit insgesamt $ 3N $ Raumkoordinaten, um das gesamte System zu beschreiben. Diese Koordinaten werden fortlaufend durchnummeriert:

$ {\begin{array}{l}{\vec {r}}_{1}=(x_{1},x_{2},x_{3})\\{\vec {r}}_{2}=(x_{4},x_{5},x_{6})\\\vdots \\{\vec {r}}_{N}=(x_{3N-2},x_{3N-1},x_{3N})\end{array}} $

Holonome Zwangsbedingungen

Holonome Zwangsbedingungen können als Gleichungen zwischen den Koordinaten $ x_{i} $ des Systems formuliert werden ($ s $: Anzahl der holonomen Zwangsbedingungen):

$ f_{l}(x_{1},x_{2},\ldots ,x_{3N},t)=0\quad \mathrm {mit} \;l=1,\ldots ,s $

Die $ 3N $ Koordinaten lassen sich mit $ s $ unabhängigen holonomen Zwangsbedingungen auf $ n=3N-s $ unabhängige generalisierte Koordinaten $ q_{i} $ reduzieren, die automatisch die Zwangsbedingungen erfüllen müssen:

$ x_{i}=x_{i}(q_{1},q_{2},\ldots ,q_{n},t)\quad \mathrm {mit} \;i=1,\ldots ,3N\;\mathrm {und} \;n=3N-s $
$ \Rightarrow f_{l}(q_{1},q_{2},\ldots ,q_{n},t)=0\quad \mathrm {mit} \;l=1,\ldots ,s\;\mathrm {und} \;n=3N-s $

Holonome Zwangsbedingungen sind mit dem vollständigen Differential einer Funktion darstellbar:

$ {\begin{aligned}{\frac {\partial f_{l}}{\partial x_{1}}}\mathrm {d} x_{1}+{\frac {\partial f_{l}}{\partial x_{2}}}\mathrm {d} x_{2}+\ldots +{\frac {\partial f_{l}}{\partial x_{3N}}}\mathrm {d} x_{3N}+{\frac {\partial f_{l}}{\partial t}}\mathrm {d} t&=0\\\Leftrightarrow a_{l1}\cdot \mathrm {d} x_{1}+a_{l2}\cdot \mathrm {d} x_{2}+\ldots +a_{l,3N}\cdot \mathrm {d} x_{3N}+a_{lt}\cdot \mathrm {d} t&=0\end{aligned}} $

und somit integrierbar.

Denn notwendig für die Integrabilität ist, dass die Koeffizientenfunktionen $ a_{li} $ folgende Integrabilitätsbedingung erfüllen:

$ {\begin{aligned}{\frac {\partial a_{li}}{\partial x_{k}}}&={\frac {\partial a_{lk}}{\partial x_{i}}}\\\Leftrightarrow {\frac {\partial ^{2}f_{l}}{\partial x_{i}\cdot \partial x_{k}}}&={\frac {\partial ^{2}f_{l}}{\partial x_{k}\cdot \partial x_{i}}}\quad \quad i,k\in \{1,\ldots ,3N\},\end{aligned}} $

was bei holonomen Bedingungen automatisch gegeben ist ($ f_{l} $ zweimal stetig differenzierbar, siehe Satz von Schwarz).

Das vollständige Differential läuft darauf hinaus, dass jede holonome Zwangsbedingung als eine Gleichung der Geschwindigkeiten darstellbar ist:

$ {\begin{alignedat}{5}&a_{l1}\cdot \mathrm {d} x_{1}&&+a_{l2}\cdot \mathrm {d} x_{2}+\ldots &&&+a_{l,3N}\cdot \mathrm {d} x_{3N}&&&&+a_{lt}\cdot \mathrm {d} t&&&&&=0\\\Leftrightarrow &a_{l1}\cdot {\frac {\mathrm {d} x_{1}}{\mathrm {d} t}}&&+a_{l2}\cdot {\frac {\mathrm {d} x_{2}}{\mathrm {d} t}}+\ldots &&&+a_{l,3N}\cdot {\frac {\mathrm {d} x_{3N}}{\mathrm {d} t}}&&&&+a_{lt}&&&&&=0\end{alignedat}} $

Anholonome Zwangsbedingungen

Nicht-holonome oder auch anholonome Zwangsbedingungen können nicht als Gleichungen zwischen den Koordinaten formuliert werden. Die generalisierten Koordinaten, die in solchen anholonomen Zwangsbedingungen erscheinen, sind i. A. nicht unabhängig voneinander variierbar.

Es handelt sich z. B. um Ungleichungen, wie Beschränkungen auf einen bestimmten Raumbereich:

$ \,f(q_{1},q_{2},\ldots ,q_{n},t)>0 $

oder um differentielle, nicht-integrable Zwangsbedingungen, wie Gleichungen zwischen den Geschwindigkeiten (Bsp. für $ r $ anholonome Zwangsbedingungen):

$ \sum _{i}a_{li}\cdot \mathrm {d} q_{i}+a_{lt}\cdot \mathrm {d} t=0\ ,\quad l=1,\ldots ,r $

Nicht-integrabel heißt dabei, dass die Gleichung – anders als bei holonomen Zwangsbedingungen – nicht als vollständiges Differential einer Funktion darstellbar ist. Somit wird hier die Integrabilitätsbedingung von den Koeffizientenfunktionen nicht erfüllt:

$ {\frac {\partial a_{li}}{\partial q_{k}}}\neq {\frac {\partial a_{lk}}{\partial q_{i}}}\quad \quad i,k\in \{1,\ldots ,n\} $

Bezüglich Zeitabhängigkeit

Weiterhin werden Zwangsbedingungen bez. ihrer Zeitabhängigkeit unterschieden in:

  • rheonom (fließend), wenn sie explizit von der Zeit abhängen.
  • skleronom (starr), wenn sie nicht explizit von der Zeit abhängen.

Skleronome Zwangsbedingungen führen bei Anwendung des Lagrange'schen Formalismus in der Regel zu der Feststellung, dass das Potential nicht implizit von der Zeit abhängt. Ist das Potential nun auch nicht explizit zeitabhängig, so sind die Kräfte konservativ und die Energie ist erhalten. In diesem Fall ist die Hamiltonfunktion – die Legendre-Transformierte der Lagrange-Funktion – gleich der Gesamtenergie.

Dagegen lassen holonom-rheonome Zwangsbedingungen nicht direkt den Schluss auf eine Nicht-Erhaltung der Energie zu.

Beispiele

Die Länge des Fadens bleibt konstant, das Pendel wird auf eine Kreisbahn gezwungen

Das Pendel: holonom und skleronom

Der Stab eines ebenen Pendels (d. h. nur 2 Raumdimensionen) soll stets die gleiche Länge $ l $ besitzen, muss also aufgrund des Satzes von Pythagoras folgende Zwangsbedingung erfüllen (Anzahl der Zwangsbedingungen: $ s=1 $):

$ f_{1}(x_{1},x_{2})=0\quad \mathrm {mit} \;x_{1}=x,x_{2}=y $
$ \Leftrightarrow x^{2}+y^{2}-l^{2}=0 $

Dabei bildet der Auslenkungswinkel $ \varphi $ des Pendels aus der Senkrechten die generalisierte Koordinate. (Es gibt nur eine, da $ n=2N-s=1 $.) Die Koordinaten $ x $ und $ y $ des Kugelmittelpunktes hängen von $ \varphi $ ab (Annahmen: $ x $ nach rechts, $ y $ nach unten, Ursprung im Aufhängungspunkt):

$ {\begin{aligned}x=l\cdot \sin \varphi \\y=l\cdot \cos \varphi \end{aligned}} $

Die generalisierte Koordinate erfüllt automatisch die Zwangsbedingung:

$ {\begin{aligned}(l\cdot \sin \varphi )^{2}+(l\cdot \cos \varphi )^{2}-l^{2}=0\\\Leftrightarrow l^{2}\cdot (\sin ^{2}\varphi +\cos ^{2}\varphi -1)=0\end{aligned}} $

da allgemein gilt:

$ \sin ^{2}\varphi +\cos ^{2}\varphi =1 $

Dies ist ein Beispiel für eine holonome Zwangsbedingung und, da sie nicht explizit von der Zeit abhängt ($ f_{1}\neq f_{1}(t) $), für eine skleronome Zwangsbedingung.

Vollständiges Differential der Zwangsbedingung:

$ {\begin{alignedat}{3}&{\frac {\partial (x^{2}+y^{2}-l^{2})}{\partial x}}\mathrm {d} x&&+{\frac {\partial (x^{2}+y^{2}-l^{2})}{\partial y}}\mathrm {d} y+{\frac {\partial (x^{2}+y^{2}-l^{2})}{\partial t}}\mathrm {d} t&&&=0\\\Leftrightarrow &2x\cdot {\frac {\mathrm {d} x}{\mathrm {d} t}}&&+2y\cdot {\frac {\mathrm {d} y}{\mathrm {d} t}}&&&=0\\\Leftrightarrow &x\cdot {\frac {\mathrm {d} x}{\mathrm {d} t}}&&+y\cdot {\frac {\mathrm {d} y}{\mathrm {d} t}}&&&=0\end{alignedat}} $

Die Geschwindigkeits-Komponenten des Pendels lassen sich in der generalisierten Koordinate wie folgt ausdrücken (aufgrund der Zwangsbedingung kann sich die Kugel nur senkrecht zum Stab bewegen; Annahme hier: Bewegung nach rechts oben):

$ {\begin{aligned}{\frac {\mathrm {d} x}{\mathrm {d} t}}&=v\cdot \cos \varphi \\{\frac {\mathrm {d} y}{\mathrm {d} t}}&=-v\cdot \sin \varphi \end{aligned}} $

mit dem Betrag $ v=l\cdot {\frac {\mathrm {d} \varphi }{\mathrm {d} t}} $ der gesamten Geschwindigkeit.

Einsetzen der generalisierten Koordinate in die Zwangsbedingung in Form des vollständigen Differentials:

$ \Rightarrow (l\cdot \sin \varphi )\cdot (v\cdot \cos \varphi )+(l\cdot \cos \varphi )\cdot (-v\cdot \sin \varphi )=0 $

das somit ebenfalls automatisch erfüllt ist.

Teilchen in Kugel: anholonom und skleronom

Ein Teilchen sei in einer Kugel eingesperrt. Das bedeutet mathematisch, dass die Entfernung des Teilchens vom Mittelpunkt der Kugel (Koordinatenursprung) stets kleiner sein muss als der Radius R der Kugel:

$ {\sqrt {x^{2}+y^{2}+z^{2}}}<R $
$ \Leftrightarrow x^{2}+y^{2}+z^{2}<R^{2} $

Da diese Zwangsbedingung aus einer Ungleichung besteht, ist sie nichtholonom, und darüber hinaus, da sie nicht explizit von der Zeit abhängt, auch skleronom.

Literatur

  • H. Goldstein: Klassische Mechanik. Wiley-VCH. ISBN 978-3527405893
  • T. Fließbach: Mechanik – Lehrbuch zur Theoretischen Physik I. Spektrum Akademischer Verlag Heidelberg 2009. ISBN 978-3-8274-2148-7

Siehe auch

en:Constraint (classical mechanics)