Када се решава проблем кретања система више тела, користи се Лагранжев формализам који упрошћава праћење еволуције система. Тачније, користе се једначине:
− − -->
d
d
t
∂ ∂ -->
L
∂ ∂ -->
q
˙ ˙ -->
i
− − -->
∂ ∂ -->
L
∂ ∂ -->
q
i
=
0
{\displaystyle -{d \over dt}{\partial L \over \partial {\dot {q}}_{i}}-{\partial L \over \partial q_{i}}=0}
где је
L
=
T
− − -->
U
{\displaystyle L=T-U}
Лагранжијан или Лагранжева функција, док је Т — кинетичка енергија система, а U — потенцијална енергија система, док су
q
˙ ˙ -->
i
{\displaystyle {\dot {q}}_{i}}
— генералисане брзине, а
q
i
{\displaystyle q_{i}}
— генералисане координате .[ 1]
Извођење
Полазимо од Даламберовог принципа да је рад сила реакција подлоге при могућем или виртуелном померању тела једнак 0, тј. ако је:
m
ν ν -->
a
→ → -->
ν ν -->
− − -->
F
→ → -->
ν ν -->
=
R
→ → -->
ν ν -->
{\displaystyle m_{\nu }{\overrightarrow {a}}_{\nu }-{\overrightarrow {F}}_{\nu }={\overrightarrow {R}}_{\nu }}
, где је
m
ν ν -->
{\displaystyle m_{\nu }}
— маса ν-тог тела, аν — његово убрзање,
F
→ → -->
ν ν -->
{\displaystyle {\overrightarrow {F}}_{\nu }}
— резултанта дејствујућих сила на ν-то тело и
R
→ → -->
ν ν -->
{\displaystyle {\overrightarrow {R}}_{\nu }}
— реакција подлоге на ν-то тело, тада је:
∑ ∑ -->
ν ν -->
=
1
N
R
→ → -->
ν ν -->
⋅ ⋅ -->
δ δ -->
r
→ → -->
ν ν -->
=
0
{\displaystyle \sum _{\nu =1}^{N}{\overrightarrow {R}}_{\nu }\cdot \delta {\overrightarrow {r}}_{\nu }=0}
, односно
∑ ∑ -->
ν ν -->
=
1
N
(
m
ν ν -->
a
→ → -->
ν ν -->
− − -->
F
→ → -->
ν ν -->
)
⋅ ⋅ -->
δ δ -->
r
→ → -->
ν ν -->
=
0
{\displaystyle \sum _{\nu =1}^{N}{\Bigl (}m_{\nu }{\overrightarrow {a}}_{\nu }-{\overrightarrow {F}}_{\nu }{\Bigr )}\cdot \delta {\overrightarrow {r}}_{\nu }=0}
1.
δ δ -->
r
→ → -->
ν ν -->
=
∑ ∑ -->
i
=
1
n
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
i
δ δ -->
q
i
{\displaystyle \delta {\overrightarrow {r}}_{\nu }=\sum _{i=1}^{n}{\partial {\overrightarrow {r}}_{\nu } \over \partial q_{i}}\delta q_{i}}
;
v
→ → -->
ν ν -->
=
d
r
→ → -->
ν ν -->
d
t
=
∑ ∑ -->
i
=
1
n
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
i
q
˙ ˙ -->
i
{\displaystyle {\overrightarrow {v}}_{\nu }={d{\overrightarrow {r}}_{\nu } \over dt}=\sum _{i=1}^{n}{\partial {\overrightarrow {r}}_{\nu } \over \partial q_{i}}{\dot {q}}_{i}}
=>
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
i
=
∂ ∂ -->
v
→ → -->
ν ν -->
∂ ∂ -->
q
˙ ˙ -->
i
{\displaystyle {\partial {\overrightarrow {r}}_{\nu } \over \partial q_{i}}={\partial {\overrightarrow {v}}_{\nu } \over \partial {\dot {q}}_{i}}}
, сада израз 1. постаје
∑ ∑ -->
ν ν -->
=
1
N
(
m
ν ν -->
d
v
→ → -->
ν ν -->
d
t
− − -->
F
→ → -->
ν ν -->
)
⋅ ⋅ -->
δ δ -->
r
→ → -->
ν ν -->
=
0
{\displaystyle \sum _{\nu =1}^{N}{\Bigl (}m_{\nu }{d{\overrightarrow {v}}_{\nu } \over dt}-{\overrightarrow {F}}_{\nu }{\Bigr )}\cdot \delta {\overrightarrow {r}}_{\nu }=0}
уз
d
v
→ → -->
ν ν -->
d
t
=
∑ ∑ -->
i
=
1
n
(
∂ ∂ -->
v
→ → -->
ν ν -->
∂ ∂ -->
q
i
q
˙ ˙ -->
i
+
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
i
q
¨ ¨ -->
i
)
{\displaystyle {d{\overrightarrow {v}}_{\nu } \over dt}=\sum _{i=1}^{n}{\Bigl (}{\partial {\overrightarrow {v}}_{\nu } \over \partial q_{i}}{\dot {q}}_{i}+{\partial {\overrightarrow {r}}_{\nu } \over \partial q_{i}}{\ddot {q}}_{i}{\Bigr )}}
, а
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
i
q
¨ ¨ -->
i
=
d
(
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
i
q
˙ ˙ -->
i
)
d
t
− − -->
∂ ∂ -->
v
→ → -->
ν ν -->
∂ ∂ -->
q
i
q
˙ ˙ -->
i
{\displaystyle {\partial {\overrightarrow {r}}_{\nu } \over \partial q_{i}}{\ddot {q}}_{i}={d{\Bigl (}{\partial {\overrightarrow {r}}_{\nu } \over \partial q_{i}}{\dot {q}}_{i}{\Bigr )} \over dt}-{\partial {\overrightarrow {v}}_{\nu } \over \partial q_{i}}{\dot {q}}_{i}}
∑ ∑ -->
ν ν -->
=
1
N
m
ν ν -->
∑ ∑ -->
i
=
1
,
j
=
1
n
(
d
(
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
i
q
˙ ˙ -->
i
)
d
t
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
j
)
δ δ -->
q
j
− − -->
∑ ∑ -->
ν ν -->
=
1
N
F
→ → -->
ν ν -->
∑ ∑ -->
j
=
1
n
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
j
δ δ -->
q
j
=
0
{\displaystyle \sum _{\nu =1}^{N}m_{\nu }\sum _{i=1,j=1}^{n}{\Bigl (}{d{\Bigl (}{\partial {\overrightarrow {r}}_{\nu } \over \partial q_{i}}{\dot {q}}_{i}{\Bigr )} \over dt}{\partial {\overrightarrow {r}}_{\nu } \over \partial q_{j}}{\Bigr )}\delta q_{j}-\sum _{\nu =1}^{N}{\overrightarrow {F}}_{\nu }\sum _{j=1}^{n}{\partial {\overrightarrow {r}}_{\nu } \over \partial q_{j}}\delta q_{j}=0}
Како је кинетичка енергија,
T
=
∑ ∑ -->
ν ν -->
=
1
N
1
2
m
ν ν -->
v
ν ν -->
2
=
∑ ∑ -->
ν ν -->
=
1
N
m
ν ν -->
∑ ∑ -->
i
=
1
,
j
=
1
n
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
˙ ˙ -->
i
q
˙ ˙ -->
i
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
˙ ˙ -->
j
q
˙ ˙ -->
j
{\displaystyle T=\sum _{\nu =1}^{N}{\frac {1}{2}}m_{\nu }v_{\nu }^{2}=\sum _{\nu =1}^{N}m_{\nu }\sum _{i=1,j=1}^{n}{\partial {\overrightarrow {r}}_{\nu } \over \partial {\dot {q}}_{i}}{\dot {q}}_{i}{\partial {\overrightarrow {r}}_{\nu } \over \partial {\dot {q}}_{j}}{\dot {q}}_{j}}
=>
m
ν ν -->
∑ ∑ -->
j
=
1
n
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
j
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
i
q
˙ ˙ -->
j
=
∂ ∂ -->
T
∂ ∂ -->
q
˙ ˙ -->
i
{\displaystyle m_{\nu }\sum _{j=1}^{n}{\partial {\overrightarrow {r}}_{\nu } \over \partial q_{j}}{\partial {\overrightarrow {r}}_{\nu } \over \partial q_{i}}{\dot {q}}_{j}={\partial T \over \partial {\dot {q}}_{i}}}
=>
∑ ∑ -->
ν ν -->
=
1
N
m
ν ν -->
∑ ∑ -->
i
=
1
,
j
=
1
n
(
d
(
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
i
q
˙ ˙ -->
i
)
d
t
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
j
)
δ δ -->
q
j
=
∑ ∑ -->
j
=
1
n
(
d
(
∂ ∂ -->
T
∂ ∂ -->
q
˙ ˙ -->
j
)
d
t
− − -->
∂ ∂ -->
T
∂ ∂ -->
q
j
)
δ δ -->
q
j
{\displaystyle \sum _{\nu =1}^{N}m_{\nu }\sum _{i=1,j=1}^{n}{\Bigl (}{d{\Bigl (}{\partial {\overrightarrow {r}}_{\nu } \over \partial q_{i}}{\dot {q}}_{i}{\Bigr )} \over dt}{\partial {\overrightarrow {r}}_{\nu } \over \partial q_{j}}{\Bigr )}\delta q_{j}=\sum _{j=1}^{n}{\Bigl (}{d{\Bigl (}{\partial T \over \partial {\dot {q}}_{j}}{\Bigr )} \over dt}-{\partial T \over \partial q_{j}}{\Bigr )}\delta q_{j}}
Па ако су силе потенцијалне, тј. важи
F
→ → -->
ν ν -->
=
− − -->
∂ ∂ -->
U
∂ ∂ -->
r
→ → -->
ν ν -->
{\displaystyle {\overrightarrow {F}}_{\nu }=-{\partial U \over \partial {\overrightarrow {r}}_{\nu }}}
, то израз
F
→ → -->
ν ν -->
∑ ∑ -->
i
=
1
n
∂ ∂ -->
r
→ → -->
ν ν -->
∂ ∂ -->
q
i
=
− − -->
∂ ∂ -->
U
∂ ∂ -->
q
j
{\displaystyle {\overrightarrow {F}}_{\nu }\sum _{i=1}^{n}{\partial {\overrightarrow {r}}_{\nu } \over \partial q_{i}}=-{\partial U \over \partial q_{j}}}
и коначно једначина 1 постаје:
∑ ∑ -->
j
=
1
n
(
d
(
∂ ∂ -->
T
∂ ∂ -->
q
˙ ˙ -->
j
)
d
t
− − -->
∂ ∂ -->
T
∂ ∂ -->
q
j
+
∂ ∂ -->
U
∂ ∂ -->
q
j
)
δ δ -->
q
j
=
0
{\displaystyle \sum _{j=1}^{n}{\Bigl (}{d{\Bigl (}{\partial T \over \partial {\dot {q}}_{j}}{\Bigr )} \over dt}-{\partial T \over \partial q_{j}}+{\partial U \over \partial q_{j}}{\Bigr )}\delta q_{j}=0}
Како су могућа померања
δ δ -->
q
j
{\displaystyle \delta q_{j}}
произвољна, то следи:
− − -->
d
d
t
(
∂ ∂ -->
L
∂ ∂ -->
q
˙ ˙ -->
i
)
+
∂ ∂ -->
L
∂ ∂ -->
q
i
=
0
{\displaystyle -{\frac {d}{dt}}{\Bigl (}{\partial L \over \partial {\dot {q}}_{i}}{\Bigr )}+{\partial L \over \partial q_{i}}=0}
Примери
Приказ клатна са дејствујућим силама, везама и генералисаном координатом
Круто тело занемарљиве масе ограничава кретање тела масе ṁ занемарљивих димензија, тако да се кретање прати замо углом θ — отклона штапа од вертикале, па се добије:
d
r
→ → -->
=
d
r
e
→ → -->
r
+
r
d
θ θ -->
e
→ → -->
θ θ -->
{\displaystyle d{\overrightarrow {r}}=dr{\overrightarrow {e}}_{r}+rd\theta {\overrightarrow {e}}_{\theta }}
=>
d
r
→ → -->
=
r
d
θ θ -->
e
→ → -->
θ θ -->
;
d
r
=
0
{\displaystyle d{\overrightarrow {r}}=rd\theta {\overrightarrow {e}}_{\theta };dr=0}
v
=
r
ϕ ϕ -->
˙ ˙ -->
{\displaystyle v=r{\dot {\phi }}}
;
T
=
1
2
m
v
2
=
1
2
m
l
2
θ θ -->
˙ ˙ -->
2
{\displaystyle T={\frac {1}{2}}mv^{2}={\frac {1}{2}}ml^{2}{\dot {\theta }}^{2}}
;
U
=
m
g
l
(
1
− − -->
cos
-->
θ θ -->
)
{\displaystyle U=mgl(1-\cos \theta )}
L
=
T
− − -->
U
=
1
2
m
r
2
ϕ ϕ -->
˙ ˙ -->
2
− − -->
m
g
l
(
1
− − -->
cos
-->
θ θ -->
)
{\displaystyle L=T-U={\frac {1}{2}}mr^{2}{\dot {\phi }}^{2}-mgl(1-\cos \theta )}
, па из
− − -->
d
d
t
(
∂ ∂ -->
L
∂ ∂ -->
θ θ -->
˙ ˙ -->
)
+
∂ ∂ -->
L
∂ ∂ -->
θ θ -->
=
0
{\displaystyle -{\frac {d}{dt}}{\Bigl (}{\partial L \over \partial {\dot {\theta }}}{\Bigr )}+{\partial L \over \partial \theta }=0}
=>
l
2
θ θ -->
¨ ¨ -->
=
− − -->
l
g
sin
-->
θ θ -->
{\displaystyle l^{2}{\ddot {\theta }}=-lg\sin \theta }
, па за
θ θ -->
→ → -->
0
{\displaystyle \theta \rightarrow 0}
произилази решење
θ θ -->
=
A
cos
-->
(
ω ω -->
t
+
ϕ ϕ -->
0
)
{\displaystyle \theta =A\cos {(\omega t+\phi _{0})}}
ω ω -->
=
g
l
{\displaystyle \omega ={\sqrt {\frac {g}{l}}}}
Кулоново поље сила припада типу централних сила, код којих је момент импулса једнак 0.
r
→ → -->
× × -->
m
v
˙ ˙ -->
→ → -->
=
r
→ → -->
× × -->
F
(
r
)
r
→ → -->
r
{\displaystyle {\overrightarrow {r}}\times m{\overrightarrow {\dot {v}}}={\overrightarrow {r}}\times F(r){\frac {\overrightarrow {r}}{r}}}
, а по својству векторског производа
r
→ → -->
× × -->
r
→ → -->
=
0
{\displaystyle {\overrightarrow {r}}\times {\overrightarrow {r}}=0}
, па је
r
→ → -->
× × -->
m
v
→ → -->
=
m
r
2
ϕ ϕ -->
˙ ˙ -->
k
→ → -->
{\displaystyle {\overrightarrow {r}}\times m{\overrightarrow {v}}=mr^{2}{\dot {\phi }}{\overrightarrow {k}}}
константа кретања.
Исти резултат лако добијамо из Лангражевог формализма:
− − -->
d
d
t
∂ ∂ -->
L
∂ ∂ -->
ϕ ϕ -->
˙ ˙ -->
− − -->
∂ ∂ -->
L
∂ ∂ -->
ϕ ϕ -->
=
0
{\displaystyle -{d \over dt}{\partial L \over \partial {\dot {\phi }}}-{\partial L \over \partial \phi }=0}
L
=
T
− − -->
U
=
1
2
m
r
2
ϕ ϕ -->
˙ ˙ -->
2
− − -->
Z
e
2
4
π π -->
ϵ ϵ -->
0
r
{\displaystyle L=T-U={\frac {1}{2}}mr^{2}{\dot {\phi }}^{2}-{\frac {Ze^{2}}{4\pi \epsilon _{0}r}}}
=>
∂ ∂ -->
L
∂ ∂ -->
ϕ ϕ -->
˙ ˙ -->
=
m
r
2
ϕ ϕ -->
˙ ˙ -->
{\displaystyle {\partial L \over \partial {\dot {\phi }}}=mr^{2}{\dot {\phi }}}
,јер је
∂ ∂ -->
L
∂ ∂ -->
ϕ ϕ -->
=
0
{\displaystyle {\partial L \over \partial \phi }=0}
Ẕ — број протона у језгри атома или редни број атома, ṁ — маса електрона , е — наелектрисање електрона, ε0 — диелектричка пропустљивост вакуума.
Референце