Bertrandov izrék [bertránov ~] v klasični mehaniki pravi, da le za dva tipa potencialov obstajajo stabilni sklenjeni tiri (orbite), za obratno kvadratno centralno silo, kot sta gravitacijski ali elektrostatski potencial:
![{\displaystyle V({\vec {\mathbf {r} }})={\frac {-k}{r}}\!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/916a9e58a85496531eb42043491c6e254bc3d771)
in za preprost potencial radialnega harmoničnega oscilatorja:
![{\displaystyle V({\vec {\mathbf {r} }})={\frac {kr^{2}}{2}}\!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/4dfb9c443f2b234d59ce2e0762cdfbf82bce9599)
Izrek se imenuje po Josephu Louisu Françoisu Bertrandu, ki ga je leta 1873 objavil.[1][2]
Vse privlačne centralne sile lahko povzročajo krožne tire, ki so seveda sklenjene. Edina zahteva je, da je centralna sila točno enaka centripetalni sili, ki določa ustrezno kotno hitrost za dani krožni polmer. Necentralne sile – tiste, ki so odvisne od kotnih spremenljivk, in tudi od polmera – se tukaj ne upoštevajo, ker v splošnem ne povzročajo krožnih tirov.
Enačba gibanja na polmeru
za delec z maso
, ki se giblje v centralnem potencialu
, je dana z Euler-Lagrangeevimi enačbami:
![{\displaystyle m{\frac {\mathrm {d} ^{2}r}{\mathrm {d} t^{2}}}-mr\omega ^{2}=m{\frac {\mathrm {d} ^{2}r}{\mathrm {d} t^{2}}}-{\frac {\Gamma ^{2}}{mr^{3}}}=-{\frac {\mathrm {d} V}{\mathrm {d} r}}\!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/f01b01eb1d40a22a1a330cb34bfb20565c5d5c57)
Pri tem se
in vrtilna količina
ohranjata. Prvi člen na levi strani je za krožne tire enak 0, sila
, ki deluje navzven, je enaka centripetalni sili
, kot se pričakuje.
Definicija vrtilne količine omogoča spremembo odvisne spremenljivke iz
v
:
![{\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}={\frac {\Gamma }{mr^{2}}}{\frac {\mathrm {d} }{\mathrm {d} \theta }}\!\,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/75822c90bf6b5bd3e78d1b32444350156556b7fb)
kar da novo gibalno enačbo neodvisno od časa:
![{\displaystyle {\frac {\Gamma }{r^{2}}}{\frac {\mathrm {d} }{\mathrm {d} \theta }}\left({\frac {L}{mr^{2}}}{\frac {\mathrm {d} r}{\mathrm {d} \theta }}\right)-{\frac {\Gamma ^{2}}{mr^{3}}}=-{\frac {\mathrm {d} V}{\mathrm {d} r}}\!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/bd739303c0e6fb20d0b9d88f629af50afeab2ef9)
Ta enačba postane kvazilinearna pri zamenjavi spremenljivk
in množenju obeh strani z
(glej tudi Binetova enačba):
![{\displaystyle {\frac {\mathrm {d} ^{2}u}{\mathrm {d} \theta ^{2}}}+u=-{\frac {m}{\Gamma ^{2}}}{\frac {\mathrm {d} }{\mathrm {d} u}}V\left({\frac {1}{u}}\right)\!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/86f0a8904c72cc5b518ea8ee75ed3e0f22c702ca)
Kot je omenjeno zgoraj, lahko centralne sile za dano začetno hitrost povzročajo krožne tire. Vendar pri določeni radialni hitrosti ti tiri niso nujno stabilni ali sklenjeni. Stabilne in strogo sklenjene tire lahko povzročajo le obratne kvadratne sile in potencial radialnega harmoničnega oscilatorja (pokazana sta potreben in zadosten pogoj).
V enačbo za
se zaradi zgoščenega zapisa uvede funkcijo
:
![{\displaystyle {\frac {\mathrm {d} ^{2}u}{\mathrm {d} \theta ^{2}}}+u=J(u)\equiv -{\frac {m}{\Gamma ^{2}}}{\frac {\mathrm {d} }{\mathrm {d} u}}V(1/u)=-{\frac {m}{\Gamma ^{2}u^{2}}}f(1/u)\!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/d103e6788ffee9e94385b59194ecb44688a65337)
kjer
predstavlja radialno silo. Po kriteriju za popolno krožno gibanje pri polmeru
mora biti prvi člen na levi strani enak 0:
![{\displaystyle u_{0}=J(u_{0})=-{\frac {m}{\Gamma ^{2}u_{0}^{2}}}f(1/u_{0})\!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/f7335802e4262fc1fea3620e88f38846c6d608b7)
kjer je
.
V naslednjem koraku se obravnava enačbo za
pri majhnih motnjah
iz popolnoma krožnih tirov. Na desni strani se lahko funkcijo
razvije v standardno Taylorjevo vrsto:
![{\displaystyle J(u)\approx u_{0}+\eta J^{\prime }(u_{0})+{\frac {1}{2}}\eta ^{2}J^{\prime \prime }(u_{0})+{\frac {1}{6}}\eta ^{3}J^{\prime \prime \prime }(u_{0})+\ldots \!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/039f949f39f1b407c57635c102a9d419972d9c4f)
Če se vstavi ta razvoj v enačbo za
in odšteje konstantne člene, se dobi:
![{\displaystyle {\frac {\mathrm {d} ^{2}\eta }{\mathrm {d} \theta ^{2}}}+\eta =\eta J^{\prime }(u_{0})+{\frac {1}{2}}\eta ^{2}J^{\prime \prime }(u_{0})+{\frac {1}{6}}\eta ^{3}J^{\prime \prime \prime }(u_{0})\ldots \!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/0408248fa8c933bf71c4c23a4469cda1c2304fd3)
kar se lahko zapiše kot:
![{\displaystyle {\frac {\mathrm {d} ^{2}\eta }{\mathrm {d} \theta ^{2}}}+\beta ^{2}\eta ={\frac {1}{2}}\eta ^{2}J^{\prime \prime }(u_{0})+{\frac {1}{6}}\eta ^{3}J^{\prime \prime \prime }(u_{0})\ldots \!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/31bfb1428ef8de422afad5ddab971bf147612d98)
kjer je
konstanta.
mora biti nenegativna, drugače se bo polmer tira spreminjal eksponentno od vrednosti začetnega polmera. (Rešitev
odgovarja popolnoma krožnemu tiru.) Če se lahko zanemari desno stran (npr. pri zelo malih motnjah), so rešitve:
![{\displaystyle \eta (\theta )=h_{1}\cos \beta \theta \!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/d180d947c5029ccb7180293aa7c2b92091efff7f)
kjer je
integracijska konstanta. Da so tiri sklenjeni, mora biti
racionalno število. Mora biti tudi enako racionalno število za vse polmere, saj se
ne more vseskozi spreminjati; racionalna števila so med seboj popolnoma nepovezana. Ker morajo enačbe iz definicije:
![{\displaystyle J^{\prime }(u_{0})\equiv -2+{\frac {u_{0}}{f(1/u_{0})}}{\frac {\mathrm {d} f}{\mathrm {d} u}}=1-\beta ^{2}\!\,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/e6dd4a87bea7678a343dfebe9009f00a793e64bb)
veljati za poljubno vrednost
, se lahko zapiše:
![{\displaystyle {\frac {\mathrm {d} f}{\mathrm {d} r}}=\left(\beta ^{2}-3\right){\frac {f}{r}}\!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/b05ebcadbdd7c3de650838b5eddf871ae13dc8bc)
od koder sledi, da mora za silo veljati potenčni zakon:
![{\displaystyle f(r)=-{\frac {k}{r^{3-\beta ^{2}}}}\!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/4e3df6fc68fa559caf93c0331087cbdf31a1095f)
Zaradi tega mora
imeti splošno obliko:
![{\displaystyle J(u)={\frac {mk}{\Gamma ^{2}}}u^{1-\beta ^{2}}\!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/49015fbd037087e1e0f6e474bbfe5c79b047ea99)
Za bolj splošne razlike od krožnosti, (če se npr. v razvoju
v Taylorjevo vrsto ne da zanemariti člene višjih redov), se lahko
razvije v Fourierovo vrsto, na primer:
![{\displaystyle \eta (\theta )=h_{0}+h_{1}\cos \beta \theta +h_{2}\cos 2\beta \theta +h_{3}\cos 3\beta \theta +\ldots \!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/3055bca01e4705abcc96cbf2941563f7e22ad93e)
Če se to rešitev vstavi v obe strani enačbe za
in izenači koeficiente, ki pripadajo isti frekvenci, se dobi sistem enačb:
![{\displaystyle h_{0}=h_{1}^{2}{\frac {J^{\prime \prime }(u_{0})}{4\beta ^{2}}}\!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/a2101a77571e83c48f09ef6f9b5e46d825552a1e)
![{\displaystyle h_{2}=-h_{1}^{2}{\frac {J^{\prime \prime }(u_{0})}{12\beta ^{2}}}\!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/c2145a9db917835775e08ced6f251d0f481dfbef)
![{\displaystyle h_{3}=-{\frac {1}{8\beta ^{3}}}\left[h_{1}h_{2}{\frac {J^{\prime \prime }(u_{0})}{2}}+h_{1}^{3}{\frac {J^{\prime \prime \prime }(u_{0})}{24}}\right]\!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/40c6cb2f342501219c8d18893db28fd48375bd9e)
in, kar je najpomembnejše:
![{\displaystyle \left(2h_{1}h_{0}+h_{1}h_{2}\right){\frac {J^{\prime \prime }(u_{0})}{2}}+h_{1}^{3}{\frac {J^{\prime \prime \prime }(u_{0})}{8}}=0\!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/1d60ce287804c7b79c9e1b117ba8526b2381b9bb)
Zadnja enačba skupaj z enačbo za
, izražena z
, vodi do glavnega rezultata Bertrandovega izreka:
![{\displaystyle \beta ^{2}\left(1-\beta ^{2}\right)\left(4-\beta ^{2}\right)=0\!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/aeb9231e256bf96d704873bcfc787cbbfba3a7dc)
Tako so edini potenciali, ki lahko povzročajo stabilne, sklenjene, nekrožne tire, obratni kvadratni zakon sile (
) in potencial radialnega harmoničnega oscilatorja (
). Rešitev
odgovarja popolnoma krožnim tirom, kar je omenjeno zgoraj; negativni rešitvi
pa nimata fizikalnega pomena.
Za obratni kvadratni zakon sile, kot sta gravitacijski ali elekstrostatični potencial, se lahko zapiše potencial kot:
![{\displaystyle V({\vec {\mathbf {r} }})={\frac {-k}{r}}=-ku\!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/e410977a291a0ebff89e59e79f96f105c0a05f17)
Tir
se lahko izpelje iz splošne enačbe:
![{\displaystyle {\frac {\mathrm {d} ^{2}u}{\mathrm {d} \theta ^{2}}}+u=-{\frac {m}{L^{2}}}{\frac {\mathrm {d} }{\mathrm {d} u}}V(1/u)={\frac {km}{L^{2}}}\!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/5c73dc5eab245479ec56f1e99dff140501d975d6)
katere rešitev je konstanta
s preprosto sinusoido:
![{\displaystyle u\equiv {\frac {1}{r}}={\frac {km}{L^{2}}}\left[1+e\cos \left(\theta -\theta _{0}\right)\right]\!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/aafb8110e4ed864e1fad36c13157b1fab7f18de6)
Tu sta
(izsrednost) in
(fazni premik) integracijski konstanti.
To je splošna enačba za stožnico z goriščem v izhodišču;
odgovarja krožnici,
elipsi,
paraboli,
pa hiperboli. Izsrednost
je povezana s skupno energijo
(glej na primer Laplace-Runge-Lenzov vektor):
![{\displaystyle e={\sqrt {1+{\frac {2EL^{2}}{k^{2}m}}}}\!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/8330c155de4f892c6a3dc7d8533e41e388c2a04a)
Primerjava teh enačb kaže, da
odgovarja elipsi,
paraboli,
pa hiperboli. Posebni primer, ko je:
![{\displaystyle E=-{\frac {k^{2}m}{2L^{2}}}\!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/e62f680b64f2a6eaad55732f9b75f7e9c832bf4a)
odgovarja popolnoma krožnim tirom.
Računanje tira v potencialu radialnega harmoničnega oscilatorja je lažje z vektorskimi komponentami
. Potencialno energijo se lahko zapiše kot:
![{\displaystyle V({\vec {\mathbf {r} }})={\frac {1}{2}}kr^{2}={\frac {1}{2}}k\left(x^{2}+y^{2}+z^{2}\right)\!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/9c775da7ad7f2b3b349aaf287946cb0bcfde2fea)
Enačba gibanja za telo z maso
je dana s tremi neodvisnimi Euler-Lagrangeevimi enačbami:
![{\displaystyle {\frac {\mathrm {d} ^{2}x}{\mathrm {d} t^{2}}}+\omega _{0}^{2}x=0\!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/7262889a204faada16f13c9bb9a26922b0dcd30a)
![{\displaystyle {\frac {\mathrm {d} ^{2}y}{\mathrm {d} t^{2}}}+\omega _{0}^{2}y=0\!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/787aa8cae2c6480173dc8055a48ab43510a569e3)
![{\displaystyle {\frac {\mathrm {d} ^{2}z}{\mathrm {d} t^{2}}}+\omega _{0}^{2}z=0\!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/4c0a3b3b41c536515c0bf6b309f0e8ec6f91e12a)
kjer mora biti konstanta
pozitivna, oziroma
, da so tiri omejeni in sklenjeni. Drugače bi telo odneslo v neskončnost. Rešitve tega preprostega harmoničnega oscilatorja so si podobne:
![{\displaystyle x=A_{x}\cos \left(\omega _{0}t+\phi _{x}\right)\!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/cc90f1a2d01ea6c68510c9c61102db0c44824492)
![{\displaystyle y=A_{y}\cos \left(\omega _{0}t+\phi _{y}\right)\!\,,}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/a0ab242e34c145adc04ae966dad0cbbf1eac0f0f)
![{\displaystyle z=A_{z}\cos \left(\omega _{0}t+\phi _{z}\right)\!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/8607bf23b224ff8615ab39318c34a0d883f0bf64)
Tu pozitivne konstante
,
in
predstavljajo amplitude nihanj, koti
,
in
pa njihove faze. Tir
je sklenjen, ker se ponovi ravno po periodi:
![{\displaystyle T\equiv {\frac {2\pi }{\omega _{0}}}\!\,.}](https://cdn.statically.io/img/wikimedia.org/api/rest_v1/media/math/render/svg/3074b835eac4213cc6c8a4c71fae83104d8da0f6)
Sistem je tudi stabilen, ker majhna odstopanja amplitud in faz povzročajo odgovarjajoče majhne spremembe na celotnem tiru.