^23 = ~C2Sl(l2yy(c\ ~ 1) - I2yyc\ + I2xxc\(fif - 1) + 2I2zzĄ(4 ~ 1) + hxx4C2sl) (29)
(30)
^31 — —ClC2S2(I2zzC2 — I2xxSi)
d32 = -C2S\{l2yy{c\ - 1) - hyyĄ + I2xxc\{Ą - 1) + 2I2zzc\{Ą - 1) + ifccscfesi) (31) ^33 — mą d" ^2yyC2 d" d2xxC]C2 ^2xxC\C2 l2yyC\C2 + l2zzC\C2 hzzCify (32)
W układzie jedyną siłą zachowawczą jest oddziaływanie grawitacyjne. Wektor przyspieszenia grawitacyjnego ma kierunek równoległy a zwrot przeciwny do osi z. Z tego powodu energia potencjalna elementu równa się iloczynowi jego masy i z-towej składowej pozycji środka masy. Energia potencjalna opisana jest wzorem:
V = -mili + m2(-c2l2 + li) + ms(c2(l2 d- d3 — -dmax) + li) (33)
Należy następnie obliczyć pochodne cząskowe po zmiennych stanu:
SV <50 2
SV
<50!
= 0
(34)
-m2s2l2 — m3s2(l2 + d3 —
:)
Sd3
Co potrzebne jest do zapisania równania dynamiki (następny podrozdział).
(35)
(36)
Równanie końcowe dynamiki zapisano w skróconej formie z wykorzystaniem symboli Chri-stofella.
SV
(37)
Gdzie dkj są elementami macierzy K podanej w rozdziale XXX. Natomiast cijk dane jest wzorem:
_ 1 / <5dkj Sdki Sdij.
^-2
(38)
Poszczególne elementy dane są wzorami:
Cm — 0
(39)
11