Dyn bryly 5 6 5 7


5.6. Równania Eulera.

0x08 graphic
Moment sił działających na bryłę względem środka układu inercjalnego

XYZ jest:

0x01 graphic

Wygodniej jest prowadzić obliczenia względem układu sztywno związanego z ciałem - układu nieinercjalnego X'Y'Z'. Transformacja wektora 0x01 graphic
przy przejściu z układu inercjalnego do układu obracającego się z prędkością kątową ω:

0x01 graphic
gdzie 0x01 graphic

Wektor 0x01 graphic
jest więc rozpatrywany w układzie wirującym.

0x01 graphic

Zatem równanie wektorowe na moment siły możemy zapisać dla trzech składowych:

0x01 graphic

0x01 graphic

0x01 graphic

Gdy Mx= My= Mz= 0 , to ciało wykonuje swobodny ruch obrotowy (precesja swobodna).

Kula.

Ixx= Iyy= Izz wówczas równania Eulera są następujące:

0x08 graphic
0x01 graphic

0x01 graphic

0x01 graphic

0x08 graphic
Obręcz.

Cienka obręcz o promieniu R. 0x01 graphic

Ixx= Iyy

Izz=mR2

0x01 graphic

0x01 graphic
0x01 graphic

0x01 graphic

równania Eulera

dla osi X : 0x01 graphic

dla osi Y : 0x01 graphic
0x01 graphic

dla osi Z : 0x01 graphic

co możemy odpowiednio zapisać:

0x01 graphic

0x01 graphic

0x01 graphic
ωz = const; Lz = const

podstawiając 0x01 graphic
otrzymujemy: 0x01 graphic
(1)

oraz 0x01 graphic
(2)

różniczkując te równania po czasie dostajemy: 0x01 graphic

0x01 graphic

skoro z (1) wynika, że 0x01 graphic
to ostatecznie: 0x01 graphic

Jest to równanie analogiczne do równania ruchu harmonicznego:

ma = - kx czyli 0x01 graphic
gdzie 0x01 graphic

rozwiązaniem tego równania jest: x =Asin(ωt + ϕ)

W przypadku obręczy: rozwiązaniem jest: 0x01 graphic

oraz 0x01 graphic

0x08 graphic
Jeżeli więc 0x01 graphic
to wektor 0x01 graphic
jest stały w przestrzeni, natomiast składowe x i y wektora 0x01 graphic
są prostopadłe do osi symetrii obręczy i wektor 0x01 graphic
obraca się wokół osi Z ze stałą prędkością kątową 0x01 graphic
gdzie 0x01 graphic
jest tzw. częstością precesji obręczy.

5.7. Precesja - przykłady: bąk i spin elektronu.

0x08 graphic

Bąk symetryczny.

0x01 graphic

gdzie 0x01 graphic
jest momentem pędu w układzie własnym 0x01 graphic

0x01 graphic
czyli:

skoro M = mgRsinα to

mgRsinα = ω'Lsinα

stąd 0x01 graphic

0x08 graphic
oznacza to, że prędkość kątowa precesji jest stała !!

0x01 graphic
0x01 graphic

0x01 graphic
czyli

0x01 graphic

prędkość kątowa precesji:

0x01 graphic

0x01 graphic

0x08 graphic
Precesja spinu.

Elektron krążący po orbicie.

Moment pędu (orbitalnego):

0x01 graphic
czyli

L = rmV =r2mω

Moment magnetyczny elektronu o wartości: 0x01 graphic

ponieważ 0x01 graphic
więc 0x01 graphic

0x01 graphic
ostatecznie: 0x01 graphic

współczynnik giromagnetyczny 0x01 graphic
czyli 0x01 graphic

0x08 graphic
Jeżeli teraz przyłożymy zewnętrzne pole magnetyczne B (0, 0, Bz)

Działa wówczas moment siły:

0x01 graphic

czyli jest precesja spinu !!

0x01 graphic

stąd 0x01 graphic

a więc Mx= γLyBz; My= -γLxBz; Mz= 0 czyli rzut 0x01 graphic
na oś Z jest zachowany.

Częstość precesji spinu Ω = γ B

Lx(t) = A sinΩt; Ly(t) = A cosΩt; Lz = const

Rezonans spinowy ⇒ elektronowy paramagnetyczny

jądrowy ferromagnetyczny

NMR

Pole magnetyczne o o zmiennej z częstością ω składowej w kierunku X:

0x01 graphic

Równania ruchu obliczamy na podstawie działającego momentu siły:

0x01 graphic
skąd Mx = γLyBz;

My = -γLxBz + LzγB1sinωt czyli My = -ΩLx + γLzB1sint

Mz = -γLyB1sinωt

Założenie: Lz >> Ly to kąt precesji Θ jest mały

Bz >> Bx = B1sinωt

Szukamy rozwiązań Lx(t) i Ly(t)

Lx= Asinωt; Ly= C cosωt gdzie ω jest takie samo jak ω zmiennego pola w kierunku X !

Warunek rezonansu: A i C maja wartość maksymalną !

0x01 graphic
0x01 graphic

0x01 graphic
0x01 graphic

0x08 graphic
stąd : 0x01 graphic

0x01 graphic

Lz = const

gdzie Ω = γ B

Rezonans spinowy zachodzi gdy

ω = Ω

0x01 graphic

Skoro więc 0x01 graphic
to 0x01 graphic

czyli np. Ixxωx= const a zatem ωx= const czyli ogólnie

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic

0x01 graphic



Wyszukiwarka

Podobne podstrony:
Dyn bryly 5.1-5.3, Fizyka
Dyn bryly 5 4 5 5
MECHANIKA II DYN
A dane,inf,wiedza,uj dyn stat proc inf w zarz 2008 9
dynamika bryly sztywnej(1)
4 Dynamika bryly sztywnej id 37 Nieznany (2)
bryły, METODY OBLICZENIOWE
bryły obrotowe powtorzenie - lekcja otwarta w III g, Matematyka dla Szkoły Podstawowej, Gimnazjum
Bryły
Bryły klasa 6 test
dynamika bryły sztywnej
Blender 3D Bryły Podstawowe Podstawowe Operacje Tutorial
Bryły obrotowe cz II
Badanie ruchu bryły sztywnej na równi pochyłej
Dynamiki bryły sztywnej

więcej podobnych podstron