Dynamika (10)

background image

7.1.1. Przedmiot dynamiki

Dynamika jest działem mechaniki, który zajmuje się badaniem zależności
między ruchem ciał materialnych i siłami wywołującymi ten ruch. Podstawą
dynamiki są prawa Newtona przytoczone w punkcie 1.2. Aby prawa te były
słuszne, w mechanice newtonowskiej ruch odnosimy do układów inercjalnych.
Z tych praw wynika, że dotyczą one punktu materialnego. W dynamice prawa
te będziemy stosować nie tylko do punktu materialnego, ale także

− po ich

odpowiednim przekształceniu

− do układu punktów materialnych, ciała sztywnego

i bryły sztywnej.
Badanie ruchu punktu materialnego o masie m i przyśpieszeniu a, na który
działa siła F, sprowadza się do analizy drugiego prawa Newtona:

F

a

=

m

. (7.1)

Powyższe równanie jest dynamicznym
równaniem ruchu punktu materialnego.
Jeżeli wektor wodzący
rozpatrywanego punktu materialnego
poprowadzony z

początku O

nieruchomego układu współrzędnych x,
y, z (rys. 7.1) oznaczymy przez r, to,
jak wiadomo z kinematyki,
przyśpieszenie a jest drugą pochodną
względem czasu wektora wodzącego.
Zatem równanie (7.1) przyjmie postać:

z

y

x

O

F

m

r

Rys. 7.1. Ruch punktu materialnego pod

działaniem siły

F

r =

2

2

t

d

d

m

. (7.2)

Jest to wektorowe równanie różniczkowe ruchu punktu materialnego. W
prostokątnym układzie współrzędnych, przedstawionym na rys. 7.1, równaniu temu
odpowiadają trzy skalarne dynamiczne równania ruchu punktu materialnego.

z

2

2

y

2

2

x

2

2

F

t

d

z

d

m

,

F

t

d

y

d

m

,

F

t

d

x

d

m

=

=

=

.

(7.3)

W równaniach tych x, y, z są współrzędnymi wektora wodzącego r, czyli
współrzędnymi punktu materialnego, a F

x

, F

y

, F

z

współrzędnymi siły F w

przyjętym układzie współrzędnych.


background image

Dynamiczne równania ruchu punktu materialnego (7.3) są w ogólnym

przypadku układem trzech równań różniczkowych i stanowią podstawę analizy
dynamiki punktu materialnego. Rozróżniamy tutaj dwie grupy zagadnień, które
omówimy w następnych punktach.

background image

7.1.2. Pierwsze podstawowe zagadnienie dynamiki

Pierwsze podstawowe zagadnienie dynamiki polega na wyznaczaniu siły
działającej na poruszający się znanym ruchem punkt materialny. Jest ono również
znane jako zagadnienie proste dynamiki. Jego rozwiązanie wynika bezpośrednio
z drugiego prawa Newtona i nie nastręcza większych trudności. Jeżeli znamy
równanie ruchu punktu materialnego w postaci:

( )

,

t

r

r

=

to w wyniku dwukrotnego różniczkowania względem czasu otrzymujemy
przyśpieszenie tego punktu:

t

d

d

2

r

a

=

i po podstawieniu tej zależności do równania (7.1) otrzymujemy siłę, a właściwie
wypadkową wszystkich sił działających na dany punkt:

2

2

t

d

d

m

r

F

=

. (7.4)

Przykład 7.1. Punkt materialny o masie m porusza się w płaszczyźnie xy
zgodnie z równaniami ruchu:

t

4sin

=

y

t,

cos2

3

x

π

π

=

, gdzie t jest czasem.

Wyznaczyć współrzędne siły działającej na ten punkt w funkcji współrzędnych
punktu x, y.

Rozwiązanie. Po zrzutowaniu wektorów występujących w równaniu (7.4) na
osie x i y otrzymujemy współrzędne siły działającej na nasz punkt materialny,
które wyrażają wzory:

.

t

d

y

d

m

F

t

d

x

d

m

F

2

2

y

2

2

x

=

=

,

(a)

Po dwukrotnym zróżniczkowaniu względem czasu równań ruchu otrzymujemy:

.

,

y

t

sin

4

dt

y

d

x

4

=

t

cos2

12

dt

x

d

2

2

2

2

2

2

2

2

π

=

π

π

=

π

π

π

=



background image

Po podstawieniu otrzymanych wyników do wzorów (a) otrzymujemy ostatecznie:

.

y

m

F

,

x

m

4

=

F

2

y

2

x

π

=

π

background image

7.1.3. Drugie podstawowe zagadnienie dynamiki


Drugie podstawowe zagadnienie dynamiki polega na wyznaczaniu ruchu punktu
materialnego poddanego działaniu znanej siły. Widzimy, że zagadnienie to jest
odwróceniem pierwszego zagadnienia dynamiki i stąd jest ono również znane pod
nazwą

zagadnienie odwrotne dynamiki.

Zagadnienie to jest znacznie trudniejsze niż pierwsze, ponieważ aby wyznaczyć
równanie ruchu punktu

( )

t

r

r

=

przy znanej sile F, należy scałkować równanie

różniczkowe (7.2) lub równoważny temu równaniu układ trzech skalarnych równań
różniczkowych (7.3). Z kursu matematyki wiadomo, że operacja taka nie jest
jednoznaczna i aby otrzymać rozwiązanie jednoznaczne, należy wyznaczyć stałe
całkowania. W tym celu musimy znać wartości funkcji i jej pochodnej (zwane
warunkami początkowymi) w pewnej chwili t

0

(w chwili początkowej):

( )

( )

0

0

0

0

t

d

t

d

,

t

v

r

r

r

=

=

. (7.5)

Znacznie większe trudności przy poszukiwaniu równania ruchu punktu
materialnego mogą wynikać z faktu, że w przypadku ogólnym siła F działająca na
punkt może być jednocześnie funkcją czasu t, położenia punktu r i prędkości v
punktu. Wtedy dynamiczne równanie ruchu punktu (7.2) należy zapisać w postaci:

(

v

r

F

r

,

,

t

t

d

d

m

2

2

=

)

. (7.6)

Rozwiązanie ogólne tego równania różniczkowego lub równoważnego mu
układu równań skalarnych w przyjętym układzie współrzędnych jest bardzo trudne
i tylko w nielicznych przypadkach udaje się otrzymać rozwiązanie ścisłe. Jeżeli nie
znamy rozwiązania równań różniczkowych, stosujemy metody przybliżone lub
numeryczne. W dalszym ciągu ograniczymy się do rozpatrzenia prostych
przykładów, w których siła F będzie stała oraz będzie funkcją tylko jednej
zmiennej

− czasu, położenia lub prędkości.


Przykład 7.2. Punkt materialny o masie m porusza się pod wpływem stałej siły
F = const. Wyznaczyć jego prędkość v = v(t) oraz równanie ruchu r = r(t); jeżeli
czas t = 0, to r(0) = r

0

i v(0) = v

0

.





background image

Rozwiązanie. Dynamiczne równanie ruchu punktu (7.2) możemy

przedstawić w postaci:

.

m

dt

d

lub

m

t

d

d

2

2

F

v

F

r

=

=

Po scałkowaniu otrzymamy:

1

t

m

dt

m

C

F

F

v

+

=

=

. (a)

Po podstawieniu w tym równaniu v = dr/dt oraz ponownym całkowaniu mamy:

.

t

t

m

2

dt

t

m

2

1

2

1

C

C

F

C

F

r

+

+

=

+

=

(b)

Stałe całkowania C

1

i C

2

wyznaczamy z podanych warunków początkowych przez

podstawienie do równań (a) i (b) r(0) = r

0

oraz v(0) = v

0

dla t = 0

C

1

= v

0

, C

2

= r

0

.

Ostatecznie prędkość punktu oraz równanie ruchu mają postać:

+

+

=

+

=

.

,

2

0

0

0

t

m

2

t

t

m

F

v

r

r

F

v

v

(c)

Z otrzymanych rezultatów wynika, że gdy siła F będzie równa zeru, to punkt
będzie się poruszał zgodnie z pierwszym prawem Newtona, czyli ruchem
jednostajnym po linii prostej.

Przykład 7.3. Punkt materialny o masie m = 1 kg porusza się po linii prostej
wzdłuż osi Ox (rys. 7.2) pod wpływem siły

( )

[ ]

F

t

=

10 1

N , gdzie t jest czasem

liczonym w sekundach. Po ilu sekundach punkt zatrzyma się i jaką drogę
przebędzie w tym czasie, jeżeli w chwili początkowej t = 0 jego prędkość
v

0

= 20 cm/s.







background image

Rozwiązanie

. Ponieważ punkt materialny porusza się wzdłuż osi Ox,

dynamiczne równanie jego ruchu możemy zapisać w postaci skalarnego równania
różniczkowego
.

( )

t

1

10

t

d

x

d

m

,

F

t

d

x

d

m

2

2

2

2

=

=

F

x

m

s

x

0

Rys. 7.2. Wyznaczenie drogi punktu

materialnego

lub

(

t

1

m

10

t

d

x

d

2

2

=

)

. (a)

Po scałkowaniu tego równania otrzymujemy prędkość punktu:

1

2

C

2

t

t

m

10

dt

dx

v

+

⎟⎟

⎜⎜

=

=

. (b)

Po podstawieniu do równania (b) warunku początkowego v = v

0

dla t = 0

wyznaczamy stałą całkowania C

1

= v

0

. Zatem prędkość punktu wyraża wzór:

⎥⎦

⎢⎣

⎟⎟

⎜⎜

+

=

⎟⎟

⎜⎜

+

=

=

s

m

2

t

t

10

2

,

0

2

t

t

m

10

v

dt

dx

v

2

2

0

.

(c)

Czas, po którym punkt się zatrzyma, obliczymy, podstawiając we wzorze (c) v = 0.
Stąd otrzymujemy równanie kwadratowe ze względu na czas t:

0

04

,

0

t

2

t

2

=

. (d)

Po obliczeniu pierwiastków tego równania i odrzuceniu pierwiastka ujemnego
otrzymujemy czas, po którym punkt się zatrzyma: t

1

= 2,02 s. Drogę przebytą przez

punkt materialny obliczymy, całkując równanie (b) w granicach od 0 do t

1

.

.

m

74

,

10

3

t

1

t

m

5

t

v

dt

2

t

t

m

10

v

s

1

1

1

0

t

0

2

0

1

=

⎛ −

+

=

⎟⎟

⎜⎜

+

=


Przykład 7.4. Punkt materialny o masie m jest przyciągany do środka O z siłą
o wartości P =

αm/x

4

(rys. 7.3), gdzie

α jest wartością stałą. Wyznaczyć prędkość

punktu w chwili, gdy jego odległość x = OM od punktu O będzie równa x

0

/2,

jeżeli w chwili początkowej (dla t = 0) x = x

0

, v = v

0

= 0.

background image

Rozwiązanie. Na rozpatrywany punkt
działa tylko siła P, wobec tego jego
równanie różniczkowe ma postać:

m

P

x

o

x

0

x

M

M

o

Rys. 7.3. Wyznaczenie prędkości

punktu materialnego

,

x

m

t

d

x

d

m

4

2

2

α

=


czyli

4

2

2

x

t

d

x

d

α

=

.

(a)

Po podstawieniu w powyższym równaniu:

v

dx

dv

dt

dx

dx

dv

dt

dv

t

d

x

d

2

2

=

=

=

otrzymamy:

,

x

dx

dv

v

4

α

=

a po rozdzieleniu zmiennych

4

x

dx

vdv

α

=

. (b)

Po scałkowaniu tego równania w granicach od 0 do v oraz od x

0

do x

0

/2

otrzymamy:

.

x

3

7

2

v

,

x

dx

vdv

3
0

2

x

2

1

x

4

v

0

0

0

α

=

α

=

Stąd prędkość punktu

3
0

x

3

14

v

α

=

. (c)

Czytelnikowi pozostawiamy wyznaczenie równania ruchu punktu.


background image

Przykład 7.5.

Ciało o masie m = 2 kg rzucone pionowo do góry z prędkością

początkową v

0

= 30 m/s pokonuje opór powietrza R, którego wartość przy

prędkości v [m/s] wynosi 0,4v [N]. Obliczyć, po ilu sekundach ciało osiągnie
najwyższe położenie H. Przyjąć przyśpieszenie ziemskie g =10 m/s

2

.

Rozwiązanie. Na ciało działają siły ciężkości
i oporu powietrza i obie są skierowane
w kierunku przeciwnym do kierunku ruchu (rys.
7.4). Zatem równanie różniczkowe ruchu ma
postać:

O

x

G

R

V

0

H

v=0

m

v

Rys. 7.4. Rzut pionowy z

uwzględnieniem oporu powietrza

,

v

4

,

0

mg

t

d

z

d

m

2

2

=

a po podstawieniu danych liczbowych możemy
napisać:

(

)

v

2

,

0

10

dt

dv

+

=

. (a)

Po rozdzieleniu zmiennych w równaniu (a) mamy:

dt

v

2

,

0

10

dv

=

+

. (b)

Po scałkowaniu tego równania w granicach od v

0

do 0 oraz od 0 do t,

uporządkowaniu i zastąpieniu różnicy logarytmów logarytmem ilorazu
otrzymujemy czas, po którym ciało osiągnie najwyższe położenie:

s.

2,35

=

ln1,6

5

10

0,2v

+

10

ln

5

t

,

dt

v

2

,

0

10

dv

2

,

0

2

0

1

0

t

0

0

v

0

=

=

=

+

,

background image

7.1.4. Zasada d’Alemberta


Po przeniesieniu obu wyrazów występujących w dynamicznym równaniu ruchu
punktu materialnego (7.1) na jedną stronę otrzymamy:

.

0

m

=

a

F

Po wprowadzeniu do tego równania zamiast

−ma fikcyjnej siły zwanej siłą

bezwładności lub siłą d’Alemberta,

P

b

m a

= −

, otrzymamy zasadę d’Alemberta

dla punktu materialnego:

0

b

=

+ P

F

, (7.7)

którą słownie wyrażamy następująco:
Suma

sił rzeczywistych i siły bezwładności działających na punkt materialny

jest w każdej chwili równa zeru.

Z zasady tej wynika, że poprzez formalne wprowadzenie siły bezwładności

zagadnienie dynamiczne można sprowadzić do zagadnienia statycznej równowagi
sił.
Przedstawioną wyżej zasadę d’Alemberta dotyczącą swobodnego punktu
materialnego zastosujemy do układu n punktów materialnych. W tym celu
rozpatrzmy układ n punktów materialnych o masach m

k

i przyśpieszeniach a

k

. Na

poszczególne punkty rozpatrywanego układu materialnego mogą działać siły
zewnętrzne i wewnętrzne.
Zgodnie z podziałem wprowadzonym w statyce (p. 3.1.2)
siłami wewnętrznymi ami rozpatrywanego układu materialnego, a siłami
zewnętrznymi siły pochodzące od innych punktów lub ciał nie należących do
naszego układu materialnego. Na rysunku 7.5 zaznaczono siły działające na dwa
punkty o masach m

k

i m

l

. Siły zewnętrzne zastąpiono tutaj siłami wypadkowymi P

k

i P

l

, a siły wzajemnego oddziaływania między tymi punktami oznaczono przez F

kl

i

F

lk

. Zgodnie z trzecim prawem Newtona siły te są równe co do wartości, ale mają

przeciwne zwroty:

.

F

F

kl

lk

= −

background image

x

z

y

r

k

m

k

O

-m

k

a

k

-m

l

a

l

r

l

F

kl

F

lk

m

l

P

k

P

l

Rys. 7.5. Siły zewnętrzne, wewnętrzne i bezwładności działające na punkty układu

materialnego


Siłę F

k

działającą na k-ty punkt możemy przedstawić w postaci sumy siły

zewnętrznej P

k

i wypadkowej wszystkich sił wewnętrznych P

wk

:

wk

k

k

P

P

F

+

=

, (7.8)

gdzie

P

wk

kl

=

=

l

l k

n

1

F

. (7.9)

Po oznaczeniu siły bezwładności działającej na rozważany punkt przez

P

a

bk

k

k

m

= −

zasadę d’Alemberta dla dowolnego punktu układu materialnego możemy
przedstawić w postaci równania:

(

)

n

,.

..

,

2

,

1

k

0

bk

wk

k

=

=

+

+

P

P

P

.

(7.10)

Suma

sił zewnętrznych, wewnętrznych oraz siły bezwładności działających na

dowolny punkt układu materialnego jest w każdej chwili równa zeru.

Jeżeli równanie (7.10) napiszemy dla każdego punktu materialnego i dodamy
stronami, to otrzymamy:

=

=

=

=

+

+

n

1

k

n

1

k

bk

n

1

k

wk

k

0

P

P

P

. (a)

background image

Występująca w tym równaniu suma wszystkich sił wewnętrznych dowolnego
układu materialnego zgodnie ze wzorem (3.3) jest równa zeru:

=

=

n

1

k

wk

0

P

. (7.11)

Zatem równanie (a) przyjmie postać:

=

=

=

+

n

1

k

n

1

k

bk

k

0

P

P

. (7.12)

Pomnóżmy teraz wektorowo każde z n równań (7.10) przez wektor wodzący r

k

i dodajmy wszystkie równania stronami. W wyniku tej operacji otrzymamy
równanie momentów:

0

n

1

k

bk

k

n

1

k

wk

k

n

1

k

k

k

=

×

+

×

+

×

=

=

=

P

r

P

r

P

r

.

(b)

Ponieważ siły wewnętrzne, zgodnie z trzecim prawem Newtona, działają parami

, i wzdłuż jednej prostej, to suma momentów tych sił dla całego układu

materialnego względem dowolnego bieguna redukcji jest równa zeru:

F

F

kl

lk

= −

0

n

1

k

wk

k

=

×

=

P

r

(7.13)

i równanie (b) przyjmuje postać:

0

n

1

k

bk

k

n

1

k

k

k

=

×

+

×

=

=

P

r

P

r

. (7.14)

Orzymane równania (7.12) i (7.14) przedstawiają zasadę d’Alemberta dla
układów materialnych, którą można sformułować następująco:

Suma

sił zewnętrznych i sił bezwładności dla danego układu materialnego oraz

sumy momentów tych sił względem nieruchomego bieguna redukcji w każdej chwili
są równe zeru.

Przykład 7.6.

Punkt materialny M o ciężarze G = 10 N, zawieszony w

nieruchomym punkcie O na lince o długości OM = s = 30 cm, tworzy wahadło
stożkowe, tzn. zatacza okrąg w płaszczyźnie poziomej, przy czym linka tworzy z
pionem kąt

(rys. 7.6a). Wyznaczyć siłę F w lince oraz prędkość v punktu

M.

α = 60

o

background image

a

G

B

F

x

O

A

y

v

α

O

A

M

a)

b)

α

Rys. 7.6. Wyznaczenie siły w lince i prędkości punktu

Rozwiązanie. Na punkt materialny działa siła ciężkości G, siła w lince F oraz
siła bezwładności (odśrodkowa)

a

B

m

=

, gdzie a jest przyśpieszeniem

dośrodkowym (rys. 7.6b). Zgodnie z zasadą d’Alemberta (7.10) suma tych sił musi
być równa zeru:

0

=

+

+

B

F

G

.

Z rzutu tych sił na osie x i y otrzymujemy dwa równania równowagi:

⎪⎭

α

=

α

=

0.

=

G

cos

F

P

,

0

=

ma

+

sin

F

P

ky

kx

(a)

Z drugiego równania otrzymujemy siłę w lince:

N

20

cos60

10

cos

G

F

o

=

=

α

=

.

Po podstawieniu do pierwszego równania (a) wzoru na przyśpieszenie
dośrodkowe:

α

=

=

sin

s

v

AM

v

a

2

2

otrzymamy równanie:

0

=

sin

s

v

g

G

+

sin

F

2

α

α

.

Stąd prędkość punktu M

s

m

1

,

2

sin60

cos60

0,3

9,81

=

sin

cos

s

g

=

sin

s

g

G

F

v

o

o

/

=

α

α

α

=

.

background image

7.1.5. Praca. Praca w zachowawczym polu sił. Energia potencjalna


Pracą mechaniczną nazywamy energię dostarczoną z zewnątrz za pomocą
układu sił do rozpatrywanego układu
materialnego w czasie jego ruchu.
Celem ogólnego zdefiniowania
pracy rozpatrzymy

ruch punktu

materialnego po torze
krzywoliniowym pod wpływem siły
P. Punkt przyłożenia A siły P jest
opisany wektorem wodzącym r
(rys. 7.7).
Pracą elementarną siły P na
przesunięciu elementarnym ds,
równym przyrostowi promienia
wodzącego dr, nazywamy iloczyn
skalarny siły P i przemieszczenia dr:

x

z

y

O

P

A

A

1

A

2

dr

α

r

Rys. 7.7. Ilustracja do definicji pracy

r

P d

dL

=

(7.15)

lub korzystając z definicji iloczynu skalarnego

(

)

dr

cos

P

cos

dr

P

dL

α

=

α

=

.

(7.16)

Jednostką pracy w układzie SI jest dżul równy pracy 1 niutona na przesunięciu

1 metra:

J = N

⋅ m = kg ⋅ m

2

⋅ s

–2

,

a w układzie technicznym kilogram siły razy metr:

1 kG

⋅m = 9,81 J.

Mimo oznaczenia pracy elementarnej symbolem powszechnie używanym na
oznaczenie różniczki zupełnej należy pamiętać, że praca elementarna nie jest na
ogół różniczką zupełną żadnej funkcji.
Na podstawie wzorów (7.15) i (7.16) można sformułować poniższe wnioski.

a) Pracę wykonuje jedynie składowa siły styczna do toru, a praca składowej

normalnej jest równa zeru.
b) Wartość pracy może być zarówno dodatnia, jak i ujemna: dla D 

dodatnia, a dla

α> π/2 ujemna.

c) Jeżeli na punkt materialny działa układ sił P

k

, których suma jest równa

wypadkowej

, to praca tej siły na przesunięciu elementarnym d

r jest

równa sumie prac elementarnych poszczególnych sił na tym przesunięciu:

W

P

=

=

k

k

n

1

background image

r

P

r

P

r

P

r

W

d

d

d

d

dL

n

2

1

+

+

+

=

=

.

d) Praca elementarna siły

P na przesunięciu wypadkowym

jest

równa sumie prac elementarnych tej siły na przesunięciach składowych:

=

=

n

1

k

k

d

d

r

r

n

2

1

d

d

d

d

dL

r

P

r

P

r

P

r

W

+

+

+

=

=

.

Jeżeli wektory występujące po prawej stronie równania (7.15) przedstawimy za
pomocą współrzędnych:

,

dz

dy

dx

d

,

P

P

P

z

y

x

k

j

i

r

j

j

i

P

+

+

=

+

+

=

to pracę elementarną możemy przedstawić w postaci:

dz

P

dy

P

dx

P

dL

z

y

x

+

+

=

. (7.17)

Jeżeli punkt przyłożenia A siły

P przemieści się po krzywej od punktu A

1

do

A

2

, to na podstawie wzoru (7.17) praca wykonana przez siłę

P będzie całką

krzywoliniową:

(

)

+

+

=

=

2

1

2

1

A

A

A

z

y

x

12

dz

P

dy

P

dx

P

d

L

A

r

P

. (7.18)

Występująca w powyższym wzorze siła może w ogólnym przypadku być
funkcją czasu t, położenia w przestrzeni punktu A oraz prędkości tego punktu.
Współrzędne siły

P będą zatem funkcjami czasu, zmiennych x, y, z oraz ich

pochodnych względem czasu. Wtedy we wzorze (7.18) możemy podstawić:

dt

dt

dz

dz

,

dt

dt

dy

dy

,

dt

dt

dx

dx

=

=

=

i zamiast całki krzywoliniowej otrzymamy całkę oznaczoną w granicach
całkowania od t

1

do t

2

+

+

=

2

1

t

t

z

y

x

dt

dt

dz

P

dt

dy

P

dt

dx

P

L

. (7.19)

Ze

względu na zastosowania bardzo ważny jest przypadek, gdy siła

P jest

jedynie funkcją położenia (miejsca):

( )

r

P

P

=

,

a jej współrzędne są wziętymi ze znakiem minus pochodnymi cząstkowymi funkcji
U względem współrzędnych x, y, z:

background image

.

z

U

P

,

y

U

P

,

x

U

P

z

y

x

=

=

=

(7.20)

Wykażemy, że funkcja skalarna U(x, y, z) ma sens fizyczny energii. Praca
elementarna siły o współrzędnych (7.20)

⎟⎟

⎜⎜

+

+

=

⎟⎟

⎜⎜

+

+

=

=

dz

z

U

dy

y

U

dx

x

U

d

z

U

y

U

x

U

d

dL

r

k

j

i

r

P

.

Wyrażenie występujące w nawiasie po prawej stronie powyższego równania jest
różniczką zupełną funkcji U:

dz

z

U

dy

y

U

dx

x

U

dU

+

+

=

. (7.21)

Z matematyki wiadomo, że całka krzywoliniowa z różniczki zupełnej jest równa
różnicy wartości końcowej i początkowej zróżniczkowanej funkcji. Zatem pracę
wykonaną przez siłę

P na jej przemieszczeniu z punktu A

1

do A

2

wyraża wzór:

(

)

.

U

U

U

U

dU

L

2

1

1

2

A

A

12

2

1

=

=

=

(7.22)

Widzimy, że praca wykonana przez siłę opisaną wzorem (7.20) na
przemieszczeniu jej z położenia początkowego do końcowego jest równa ubytkowi
funkcji U. Funkcję tę nazywamy potencjałem albo energią potencjalną, siłę

P

spełniającą warunek (7.20) siłą potencjalną lub zachowawczą, a pole sił polem
potencjalnym
lub zachowawczym.

Potencjał w określonym punkcie przestrzeni jest równy pracy, którą wykonują

siły potencjalne przy przemieszczaniu punktu materialnego z danego punktu do
punktu, w którym potencjał jest równy zeru. Ponieważ punkt ten może być obrany
dowolnie, potencjał jest określony z dokładnością do dowolnej stałej C. Wnika to z
tego, że funkcja:

C

U

U

+

=

również spełnia zależności (7.20) i (7.22).
Ze wzoru (7.22) wynikają dwie ważne własności sił potencjalnych.
a) Praca siły potencjalnej nie zależy od toru jej punktu przyłożenia, lecz jedynie
od położenia tego punktu w chwilach początkowej i końcowej.

b) Praca wykonana przez siłę potencjalną jest równa ubytkowi energii

potencjalnej wynikającemu z przemieszczania się punktu przyłożenia siły. Wynika
stąd również, że praca po torze zamkniętym jest równa zeru.

background image

7.1.6. Przykłady sił potencjalnych

Siły sprężystości
Wykażemy obecnie, że siły odkształcenia sprężystego są siłami potencjalnymi.

W tym celu rozpatrzymy sprężynę śrubową, której koniec A jest unieruchomiony,
a koniec B może się przemieszczać wzdłuż osi Ox (rys. 7.8). Założymy, że w
chwili, gdy sprężyna nie jest napięta, koniec B pokrywa się z punktem O.

x

A

O

B

x

P

Rys. 7.8. Przykład siły sprężystej wykonującej pracę

Jeżeli wydłużymy sprężynę o wartość x, to zgodnie z prawem Hooke’a będzie
ona działać na punkt B siłą P proporcjonalną do wydłużenia:

i

P

x

k

=

,

gdzie współczynnik proporcjonalności k jest nazywany stałą sprężyny, a znak
minus oznacza, że siła P jest skierowana przeciwnie do kierunku odkształcenia
sprężyny.
Z

powyższego wzoru wynika, że współrzędna siły P jest funkcją tylko

współrzędnej x:

x

k

P

=

,

zatem potencjał U musi spełniać równanie:

x

k

P

dx

dU

x

U

=

=

=

.

Po scałkowaniu tego równania w granicach od O do x

1

otrzymujemy wzór na

potencjał siły sprężystej:

2

1

x

0

x

k

2

1

x

k

U

1

=

=

. (7.23)

Pracę siły sprężystej na skończonym przesunięciu, np. od 0 do x, można
obliczyć ze wzoru (7.22), przy czym dla x = 0 energia potencjalna U

1

= 0. Zatem

2

1

2

12

x

k

2

1

U

L

=

=

.

(7.24)


background image

Siły ciężkości
Jeżeli rozpatrzymy ograniczony obszar przestrzeni w pobliżu powierzchni

Ziemi o małych wymiarach w porównaniu z promieniem Ziemi, to można przyjąć,
że na każdy punkt materialny o masie m znajdujący się w tej przestrzeni działa
stała siła ciężkości:

G = mg,

gdzie g jest przyśpieszeniem ziemskim. Przy takim założeniu pole sił jest
jednorodnym polem sił ciężkości. Gdy w takim polu sił przyjmiemy układ
współrzędnych x, y, z o osi z skierowanej pionowo w górę, to zgodnie z rys. 7.9
współrzędne siły ciężkości G opisują zależności:

.

mg

G

,

0

G

G

z

y

x

=

=

=

(7.25)

Ze wzoru (7.20) wiadomo, że współrzędne sił potencjalnych są równe
pochodnym cząstkowym potencjału U względem współrzędnych wziętych ze
znakiem minus:

mg

z

U

G

,

0

y

U

G

,

0

x

U

G

z

y

x

=

=

=

=

=

=

. (7.26)

Z powyższych równań wynika, że potencjał U jest jedynie funkcją zmiennej z. Po
podstawieniu trzeciego równania (7.26) do wzoru (7.21) otrzymujemy różniczkę
potencjału pola sił ciężkości:

,

dz

mg

dU

=

a po scałkowaniu tego równania potencjał sił ciężkości

C

z

g

m

U

+

=

, (7.27)

gdzie C jest dowolną stałą.
Ze wzoru (7.27) wynika, że dla z = const potencjał U jest również stały. Zatem
w przypadku sił ciężkości wszystkie punkty każdej płaszczyzny poziomej mają
taką samą wartość potencjału. Powierzchnie, których punkty mają te same wartości
potencjału, nazywają się powierzchniami ekwipotencjalnymi.
Praca

siły ciężkości na dowolnym krzywoliniowym torze jest

− zgodnie ze

wzorem (7.22)

− równa różnicy potencjałów w położeniu początkowym i

końcowym:

(

)

h

g

m

z

z

g

m

U

U

L

2

1

2

1

12

=

=

=

,

(7.28)

gdzie h jest różnicą wysokości (rys. 7.9).

background image

x

y

z

O

A

1

A

2

G

h

A


Rys. 7.9. Praca siły ciężkości

x

z

y

r

P

A

M

O

m

Rys. 7.10. Siła wzajemnego przyciągania

Siły wzajemnego przyciągania
Wykażemy, że siła, z jaką nieruchomy punkt materialny o masie M działa na

dowolny punkt materialny o masie m, jest siłą potencjalną. Zgodnie z prawem
powszechnego ciążenia (1.2) punkt M działa na punkt m i odwrotnie z siłą P o
wartości

2

r

Mm

k

P

=

, (7.29)

gdzie k jest stałą grawitacji, a r jest odległością masy m od nieruchomej masy M.
Jeżeli masę M umieścimy w początku układu współrzędnych x, y, z, a masę m
w punkcie A o wektorze wodzącym r (rys. 7.10), to siłę P można opisać wzorem:

r

1

P

2

r

Mm

k

=

, (7.30)

gdzie 1

r

jest wektorem jednostkowym o kierunku wektora r.

Gdy współrzędne wektora wodzącego r oznaczymy przez x, y, z, to
współrzędne siły P będą następujące:

r

z

r

Mm

k

P

,

r

y

r

Mm

k

P

,

r

x

r

Mm

k

P

2

z

2

y

2

x

=

=

=

. (7.31)

Łatwo wykazać, że potencjałem omawianego pola sił jest funkcja

(

)

U x

k

Mm

r

C

k

Mm

x

y

z

C

, y, z

= −

+ = −

+

+

+

2

2

2

.

(7.32)

background image

przy czym C jest dowolną stałą. Aby siła P była potencjalna, jej współrzędne
(7.31) muszą spełniać wzory (7.20). Po zróżniczkowaniu funkcji (7.32) względem
x otrzymamy:

(

)

x

2

3

2

3

2

2

2

P

r

x

r

Mm

k

r

kMmx

z

y

x

x

2

2

1

kMm

x

U

=

=

=

+

+

⎛−

=

.

Postępując podobnie w odniesieniu do y i z, otrzymamy:

z

2

y

2

P

r

z

r

Mm

k

z

U

,

P

r

y

r

Mm

k

y

U

=

=

=

=

.

Pracę wykonaną przez siłę P na przemieszczenie masy m z położenia 1 do 2
zgodnie ze wzorem (7.22) i po uwzględnieniu równania (7.32) zapiszemy w
następującej postaci:

⎟⎟

⎜⎜

=

=

1

2

2

1

12

r

1

r

1

kMm

U

U

L

. (7.33)

background image

7.1.7. Moc i sprawność

Z technicznego punktu widzenia interesuje nas często nie tylko wartość pracy,
ale również czas, w jakim została ona wykonana. W tym celu wprowadzono
pojęcie mocy.

Mocą chwilową nazywamy stosunek pracy elementarnej dL do czasu dt.

t

d

L

d

N

=

. (7.34)

Po podstawieniu do tego wzoru pracy elementarnej zdefiniowanej wzorem
(7.15) otrzymujemy wzór na moc siły P.

v

P

r

P

=

=

t

d

d

N

.

(7.35)

Zatem

moc siły jest równa iloczynowi skalarnemu siły P i prędkości v jej punktu

przyłożenia.
Ze wzoru (7.34) widzimy, że między pracą elementarną dL i mocą N istnieje
prosty związek:

.

dt

N

L

d

=

Jeżeli siła P w chwili t

1

znajduje się w punkcie A

1

, a w chwili t

2

w punkcie A

2

(rys. 7.6), to praca L

12

wykonana przez tę siłę przy przemieszczeniu się po torze od

A

1

do A

2

będzie równa całce z mocy w granicach od t

1

do t

2

:

=

2

1

t

t

12

Ndt

L

. (7.36)

Gdy na układ materialny działa układ n sił, to moc tego układu jest równa sumie
mocy poszczególnych sił:

=

=

n

1

k

k

N

N

. (7.37)

Podstawową jednostką mocy w układzie SI jest wat (w skrócie W). Jest to moc
siły, która pracę jednego dżula wykonuje w ciągu jednej sekundy:

1 W = J

⋅ s

–1

.

W praktyce na określenie mocy silników i maszyn są używane większe
jednostki

− kilowaty (kW) i megawaty (MW):⋅

1

kW

=

1000

W,

background image

1 MW = 1000 kW = 1 000 000 W.

W technicznym układzie jednostek podstawową jednostką mocy jest kilogram
siły razy metr na sekundę:

1 kG

⋅ m ⋅ s

–1

.

Praktyczną jednostką mocy w tym układzie jest koń mechaniczny KM:

1 KM = 75 kG

⋅ m ⋅ s

–1

.

Między jednostkami mocy w układzie technicznym i w układzie SI istnieją
zależności:
1 kG

⋅ m ⋅ s

–1

= 9,81 W,

1 KM = 75

9,81 W = 0,736 kW,

1

W

=

0,102

kG

⋅ m ⋅ s

–1

,

1

kW

=

102

kG

⋅ m ⋅ s

–1

= 1,36 KM.

Do oceny stanu silnika czy maszyny wykorzystuje się pojęcie sprawności
mechanicznej.
Wiadomo, że część mocy dostarczonej do silnika (maszyny) jest
tracona na pokonanie oporów istniejących w samym silniku (maszynie), a tylko
część jest zamieniana na moc użyteczną.
Sprawnością mechaniczną nazywamy stosunek mocy użytecznej N

u

(lub pracy

L

u

) do mocy włożonej N

w

(lub pracy L

w

):

w

u

w

u

L

L

N

N

η

=

=

. (7.38)

Sprawność jest liczbą bezwymiarową spełniającą nierówność:

.

1

η

0

background image

7.1.8. Moc układu sił działających na bryłę sztywną


W poprzednim punkcie zdefiniowaliśmy moc siły P działającej na punkt
materialny. Obecnie obliczymy moc układu n sił zewnętrznych P

k

, gdzie

k = 1, 2, .... , n, przyłożonych odpowiednio w punktach A

1

, A

2

, .... , A

n

bryły

sztywnej, poruszającej się znanym ruchem względem nieruchomego układu
współrzędnych x, y, z (rys. 7.11). W dowolnym punkcie (biegunie redukcji)

umieścimy ruchomy układ współrzędnych

O

′ ′ ′

x , y , z poruszający się razem z bryłą.

Układ sił P

k

reprezentują wektor główny W i moment główny

umieszczone

w biegunie redukcji

, a ruch bryły jest określony za pomocą prędkości

bieguna

i prędkości kątowej

ω.

M

O

O

v

O

O

x

M

O

z

x

z

y

y

O

O

W

k

r

ω

v

O′

P

k

A

k

v

k

Rys. 7.11. Wyznaczenie mocy układu sił działających na bryłę sztywną

Zgodnie z definicją moc N

k

siły

P

k

k

k

k

N

v

P

=

.

Prędkość dowolnego punktu A

k

zgodnie ze wzorem (5.29) możemy zapisać

w następujący sposób:

k

O

k

r

ω

v

v

×

+

=

.

Po podstawieniu tego wzoru do wzoru na moc N

k

siły

P

k

oraz wykorzystaniu

własności iloczynu mieszanego (2.31) otrzymujemy:

(

)

(

)

(

)

k

k

O

k

k

O

k

k

O

k

k

N

P

r

ω

v

P

r

ω

P

v

P

r

ω

v

P

×

+

=

×

+

=

×

+

=

.

background image

Moc układu sił działających na bryłę sztywną otrzymamy po zsumowaniu

zgodnie ze wzorem (7.37)

− mocy poszczególnych sił:

(

)

[

]

k

n

1

k

k

n

1

k

k

O

n

1

k

k

k

O

k

n

1

k

k

N

N

P

r

ω

P

v

P

r

ω

v

P

=

=

=

=

×

+

=

×

+

=

=

.

Ostatecznie

ω

M

v

W

+

=

O

O

N

. (7.39)

Zgodnie z zależnościami (3.25) i (3.26) w powyższym wzorze

W jest wektorem

głównym, a

momentem głównym układu sił zewnętrznych zredukowanych

do bieguna redukcji

.

M

O

O

Wzór

(7.39)

można wyrazić słownie:

Moc układu sił zewnętrznych działających na bryłę sztywną jest równa sumie

iloczynu skalarnego wektora głównego i prędkości dowolnego bieguna

redukcji

oraz iloczynu skalarnego momentu głównego zredukowanego do tegoż bieguna
i prędkości kątowej.

background image

7.2.1. Pęd układu materialnego i bryły

Pędem punktu materialnego o masie m i prędkości v nazywamy iloczyn masy

punktu i jego prędkości:

p = mv. (7.40)

Z powyższej definicji
wynika, że pęd jest wektorem o
kierunku prędkości, a więc jest
wektorem stycznym do toru
punktu materialnego.

Dla układu n punktów

materialnych o masach m

k

i

prędkości v

k

(rys.

7.12) pęd

będzie równy sumie pędów
poszczególnych punktów
materialnych:

v

n

x

v

C

v

2

r

k

z

y

r

Ck

r

C

m

k

C

O

m

1

v

1

v

k

m

2

m

n

Rys. 7.12. Wyznaczenie pędu układu materialnego

=

=

n

1

k

k

k

m v

p

. (7.41)

Wzór (7.41) można przedstawić w postaci:

n

1

=

k

k

k

m

dt

d

=

r

p

. (a)

Widzimy,

że występująca pod znakiem pochodnej suma, zgodnie ze wzorem

(4.18), jest momentem statycznym S rozpatrywanego układu materialnego
względem początku nieruchomego układu współrzędnych x, y, z :

=

n

1

=

k

C

k

k

m

m

=

r

r

S

. (b)

Po podstawieniu wzoru (b) do wzoru (a) i wykonaniu różniczkowania wzór (7.41)
możemy zapisać w postaci:

dt

d

m

m

C

n

1

k

k

k

S

v

v

p

=

=

=

=

, (7.42)

gdzie m jest masą całkowitą układu materialnego.
Z otrzymanego wzoru wynika, że pęd układu materialnego jest równy
iloczynowi masy całkowitej m układu materialnego i prędkości v

C

środka masy C.

Ponadto wzór (7.42) pozwala na inne zdefiniowanie pędu.

background image

Pędem nazywamy pochodną względem czasu momentu statycznego układu
materialnego względem nieruchomego punktu
:

dt

d S

p

=

.

(7.43)

Ponieważ moment statyczny względem środka masy jest równy zeru (patrz p.
4.4), zatem pęd układu materialnego względem środka masy jest także równy zeru.
Pęd bryły sztywnej możemy obliczyć, dzieląc ją na elementy o masach

∆m

k

i

traktując ją jako układ punktów materialnych. Przybliżoną wartość pędu
otrzymamy po zsumowaniu pędów tych elementów, traktowanych jako punkty
materialne.
Z kolei wartość dokładną pędu otrzymamy po wyznaczeniu granicy sumy, gdy
liczba elementów dąży do nieskończoności

=

=

=

=

=

m

m

m

n

1

k

k

k

k

dm

dt

d

m

dt

d

dm

m

lim

r

r

v

v

p

.

Całka występująca w tym wzorze pod znakiem pochodnej jest momentem
statycznym bryły względem początku układu współrzędnych:

C

m

m

dm

r

r

=

.

Z uwzględnieniem powyższej zależności otrzymujemy wzór na pęd bryły:

(

)

C

C

C

m

dt

d

m

m

dt

d

v

r

r

p

=

=

=

. (7.44)

Widzimy

zatem,

że pęd bryły, podobnie jak pęd układu materialnego, jest

równy iloczynowi jej masy i prędkości środka masy.

background image

7.2.2. Zasada pędu i popędu. Zasada zachowania pędu

Rozpatrzymy obecnie układ składający się z n punktów materialnych o masach
m

k

i prędkości v

k

. Na poszczególne punkty rozpatrywanego układu materialnego

działają siły zewnętrzne i
wewnętrzne. Na rysunku 7.13
zaznaczono siły działające na dwa
punkty o masach m

k

i m

l

. Siły

zewnętrzne działające na te
punkty zastąpiono siłami
wypadkowymi P

k

i P

l

, siły

wzajemnego oddziaływania
między tymi punktami oznaczono
przez F

kl

i F

lk

.

v

k

v

2

x

z

y

r

C

r

k

m

k

C

O

v

c

r

l

F

kl

F

lk

m

l

P

k

P

l

Rys. 7.13. Siły zewnętrzne i wewnętrzne działające

na punkty układu materialnego

Wypadkowa

sił wewnętrznych

działających na punkt o masie m

k

=

=

n

k

l

1

l

kl

wk

F

P

, (7.45)

a wypadkowa wszystkich sił działających na ten punkt

F

k

= P

k

+ P

wk

. (7.46)

Zatem zgodnie z drugim prawem Newtona możemy dla dowolnego punktu
rozważanego układu materialnego napisać dynamiczne równanie ruchu w postaci:

(

n

,..

.,

2

,

1

k

t

d

d

m

wk

k

2

k

2

k

=

+

=

P

P

r

)

. (7.47)

Po założeniu, że masa m

k

jest wielkością stałą, lewą stronę tego równania możemy

przedstawić w postaci pochodnej względem czasu pędu m

k

v

k

punktu:

(

)

dt

m

d

dt

d

m

t

d

d

m

k

k

k

k

2

k

2

k

v

v

r

=

=

.

Równanie (7.47) można obecnie zapisać następująco:

(

)

(

.

n

,..

.,

2

,

1

k

dt

m

d

wk

k

k

k

=

+

=

P

P

v

)

(c)

Jeżeli dodamy stronami n powyższych równań, to otrzymamy:

background image

(

)

=

=

=

+

=

n

1

k

wk

n

1

k

k

n

1

k

k

k

dt

m

d

P

P

v

,

a jeżeli zastąpimy sumę pochodnych pędów pochodną ich sumy, to

=

=

=

+

=

n

1

k

kz

n

1

k

k

n

1

k

k

m

dt

d

P

P

v

k

. (d)

Lewa strona równania (d) jest pochodną względem czasu pędu układu
materialnego:

dt

d

m

dt

d

n

1

k

k

p

v

k

=

=

.

Pierwsza suma po prawej stronie równania (d) jest wektorem głównym sił
zewnętrznych:

=

=

n

1

k

k

P

W

,

a druga sumą wszystkich sił wewnętrznych działających w całym układzie
materialnym i zgodnie ze wzorem (3.3) jest równa zeru:

.

0

n

1

k

n

k

l

1

l

kl

n

1

k

wk

∑∑

=

=

=

=

=

F

P

Ostatecznie równanie (d) można zapisać w postaci:

W

p =

dt

d

. (7.48)

Równanie to przedstawia zasadę pędu układu punktów materialnych, którą
można wypowiedzieć następująco:

Pochodna

względem czasu pędu układu punktów materialnych jest równa

wektorowi głównemu sił zewnętrznych działających na ten układ.

W celu wyznaczenia zmiany pędu układu punktów materialnych w skończonym
przedziale czasu, np. od 0 do t, wywołanej przez siły zewnętrzne działające na ten
układ, scałkujmy równanie (7.48) w tym przedziale czasu. Otrzymamy wtedy:

( ) ( )

=

t

0

dt

0

t

W

p

p

. (7.49)

background image

Równanie to nazywamy zasadą pędu i popędu lub prawem zmienności pędu.

Przyrost

pędu układu materialnego w skończonym przedziale czasu jest równy

popędowi wektora głównego sił zewnętrznych działających na ten układ.
Całkę z prawej strony równania (7.49) nazywamy popędem wektora głównego
lub impulsem wektora głównego. Ta druga nazwa ma swoje uzasadnienie
zwłaszcza w przypadku sił krótkotrwałych, np. sił zderzeniowych. Łatwo
zauważyć, że gdy wektor główny układu sił zewnętrznych jest równy zeru:

W = 0,

popęd tego wektora jest również równy zeru, a z zasady pędu i popędu wynika, iż
pęd końcowy jest równy początkowemu:

( ) ( )

0

t

p

p

=

,

czyli pęd układu materialnego jest stały:

const

=

p

. (7.50)

Jest to zasada zachowania pędu:

Jeżeli wektor główny układu sił zewnętrznych działających na układ materialny
jest równy zeru, to pęd tego układu materialnego jest stały.

Gdy

pęd układu materialnego przedstawimy w postaci iloczynu masy m

i prędkości v

C

środka masy, to z zasady zachowania pędu:

const

m

C

=

v

wynika, że środek masy porusza się ruchem jednostajnym prostoliniowym.

Przykład 7.7. Klocek o masie m = 40 kg porusza się po równi pochyłej o kącie
nachylenia

pod działaniem siły będącej funkcją czasu P = P(t)

(rys. 7.14a). Miara tej siły zmienia się w czasie od 0 do P

α = 30

o

1

= 250 N zgodnie z

wykresem podanym na rys. 7.14b. Współczynnik tarcia między klockiem i równią

. Obliczyć prędkość v



= 0 1

,

1

, jaką osiągnie ciało w chwili t

1

= 3 s, jeżeli w chwili

t = 0 prędkość początkowa v

m s

0

10

=

/ .

background image

t

P

P

1

t

1

x

P(t)

N

T

α

G

a)

b)

0

Rys. 7.14. Wyznaczenie prędkości klocka

Rozwiązanie. Do rozwiązania zadania zastosujemy zasadę pędu i popędu (7.49).
W myśl tej zasady przyrost pędu klocka w czasie od t = 0 do t = t

1

będzie równy

popędowi wektora głównego sił zewnętrznych działających na niego:

( ) ( )

=

1

t

0

1

dt

0

t

W

p

p

.

Wektory z tego równania zrzutujemy na oś x równoległą do równi. Po
uwzględnieniu zależności (7.44) mamy:

=

1

t

0

x

0

1

dt

W

mv

mv

. (a)

Zgodnie z rysunkiem suma rzutów wszystkich sił działających na klocek na oś x

α

α

=

α

=

cos

mg

µ

sin

mg

(t)

P

T

sin

mg

(t)

P

W

x

, (b)

gdzie

α

=

=

cos

mg

µ

N

µ

T

. Po podstawieniu (b) do równania (a) mamy:

(

)

(

)

.

t

µcos

sin

mg

(t)dt

P

dt

µcos

sin

mg

(t)dt

P

mv

mv

1

t

0

t

0

t

0

0

1

1

1

1

α

+

α

=

=

α

+

α

=

(c)

Całka występująca w powyższym wzorze jest równa polu wykresu
przedstawionego na rys. 7.14b, czyli

1

1

t

0

t

P

2

1

(t)dt

P

1

=

.

Po podstawieniu tej równości do (c) otrzymujemy wzór na prędkość v

1

:

background image

(

)

1

1

1

0

1

t

µcos

sin

g

m

2

t

P

v

v

α

+

α

+

=

.

Po podstawieniu danych liczbowych otrzymujemy:

(

)

s

m

1

,

2

3

30

cos

1

,

0

sin30

81

,

9

40

2

3

250

10

v

o

o

1

/

=

+

+

=

.

background image

7.2.3. Twierdzenie o ruchu środka masy


Pęd p w wyprowadzonym w poprzednim punkcie równaniu (7.48),
wyrażającym zasadę pędu, możemy przedstawić za pomocą iloczynu całkowitej
masy m układu materialnego i prędkości v

C

środka jego masy C. Otrzymamy

wówczas:

(

)

W

v

v

p

=

=

=

dt

d

m

dt

m

d

dt

d

C

C

.

(e)

Występująca w tym równaniu pochodna prędkości środka masy względem czasu
jest przyśpieszeniem środka masy. Mamy więc:

W

a

=

C

m

. (7.51)

Po zapisaniu wektorów a

C

i W w układzie współrzędnych x, y, z:

+

+

=

+

=

k

j

i

W

k

j

i

a

y

y

x

Cz

Cy

Cx

C

W

W

W

,

a

a

+

a

(f)

wektorowe równanie (7.51) możemy przedstawić w postaci trzech równań
skalarnych:

z

Cz

y

Cy

x

Cx

W

ma

,

W

ma

,

W

ma

=

=

=

.

(7.52)

Wektorowe równania (7.51) i równoważne mu trzy równania skalarne (7.52) są
dynamicznymi równaniami ruchu środka masy. Pozwalają one na wyznaczenie
ruchu środka masy pod wpływem znanych sił zewnętrznych. Otrzymane równania
(7.51) lub (7.52) pozwalają na sformułowanie twierdzenia, znanego pod nazwą
twierdzenia o ruchu środka masy.

Środek masy układu materialnego porusza się tak jak punkt materialny o masie
równej całkowitej masie układu, na który działa siła równa wektorowi głównemu
sił zewnętrznych działających na ten układ.
Twierdzenie

o

ruchu

środka masy wynika również z pierwszej całki zasady

pędu, czyli z zasady pędu i popędu przedstawionej w postaci:

( )

( )

=

t

0

C

C

dt

0

m

t

m

W

v

v

. (7.53)

Twierdzenie to jest ważnym narzędziem badania ruchu środka masy, ale nie
pozwala na wyciągnięcie żadnych wniosków co do ruchu punktów należących do
układu względem środka masy.

background image

Z twierdzenia o ruchu środka masy wynika, że siły wewnętrzne nie mogą
zmienić ruchu środka masy ani jego położenia.
Twierdzenie to odnosi się nie tylko do układu punktów materialnych, ale
również do ciała sztywnego i bryły. Nałożywszy bowiem na układ punktów
materialnych warunek, aby odległość dowolnych punktów układu była niezmienna,
otrzymujemy model ciała sztywnego.

background image

7.2.4. Ruch układu o zmiennej masie


Do tej pory w rozważaniach dotyczących pędu układu materialnego

zakładaliśmy, że całkowita masa układu nie ulega zmianie w czasie ruchu. Obecnie
zajmiemy się ruchem układu materialnego, którego masa będzie się zmieniać z
upływem czasu poprzez odłączanie lub dołączanie elementów masy. Taka zmiana
masy układu będzie miała wpływ na jego ruch.
Typowym

przykładem ruchu układu o zmiennej masie są rakiety, z których

w czasie pracy silnika następuje wypływ gazów spalinowych, a tym samym
zmniejsza się masa rakiety. Innym przykładem mogą być urządzenia do transportu
ciągłego ze zmieniającą się w czasie ilością przenoszonego materiału.
W

dalszych

rozważaniach ze zrozumiałych względów ograniczymy się

jedynie do wyprowadzenia równania ruchu ciała o zmiennej masie. Do ułożenia
równania ruchu wykorzystamy zasadę pędu (7.48) zapisaną w postaci:

(

)

W

v

=

dt

m

d

C

.

(g)

Przyjmijmy,

ze

środek układu materialnego o masie m porusza się

względem układu odniesienia z prędkością v

C

i w pewnej chwili masa układu

zaczyna się zmieniać w sposób ciągły. Zakładając, że w czasie dt od układu
odrywa się (lub przyłącza do niego) masa elementarna dm z prędkością
bezwzględną v

b

, określimy elementarną zmianę pędu. W chwili początkowej t pęd

układu wynosi

C

m v

,

(h)

a w chwili t + dt

(

)(

)

b

C

dm

d

dm

m

v

v

v

+

.

(i)

Elementarną zmianę pędu otrzymamy przez odjęcie zależności (i) od (h).

(

)

(

)(

)

[

]

(

)

.

d

dm

dm

md

dm

d

dm

dm

md

m

m

dm

d

dm

m

m

m

d

C

b

b

C

C

C

b

C

C

C

v

v

v

v

v

v

v

v

v

v

v

v

v

v

v

=

=

+

+

=

=

+

=

Po pominięciu iloczynu różniczek dmdv jako małej wartości drugiego rzędu
elementarna zmiana pędu

(

)

w

C

dm

md

m

d

v

v

v

=

,

(j)

background image

gdzie

v

w

= v

b

v

C

i jest prędkością masy dm względem masy m, czyli prędkością względną. Po
uwzględnieniu wyrażenia (h) w równaniu (e) otrzymamy równanie ruchu układu
o zmiennej masie nazywane równaniem Mieszczerskiego:

W

v

v

+

=

dt

dm

dt

d

m

w

C

lub w postaci

W

R

v

+

=

dt

d

m

C

,

(7.54)

gdzie

dt

dm

w

v

R

=

(7.55)

i jest reakcją cząstki elementarnej.
Jeżeli występująca we wzorze (7.55) pochodna

> 0, czyli masa

układu wzrasta z upływem czasu, to wektor R ma zwrot prędkości względnej v

dt

dm /

w

i

jest siłą hamującą. Gdy masa układu materialnego będzie maleć z upływem czasu,
czyli dm/dt < 0, to wektor R będzie miał zwrot przeciwny do prędkości względnej
v

w

,

a więc będzie siłą napędową.
Jeżeli równanie (7.54) zastosujemy do badania ruchu rakiety i założymy,
że wektor prędkości względnej v

w

wypływających z rakiety gazów jest styczny do

trajektorii lotu, to wektor R będzie siłą ciągu rakiety (rys. 7.15).

R

v

C

v

w

W

Rys. 7.15. Ruch układu o zmiennej masie

Przykład 7.8. Rakieta o masie początkowej m

0

porusza się w przestrzeni

międzyplanetarnej z prędkością początkową v

C0

. Po włączeniu silnika prędkość

względna v

w

wypływających z rakiety produktów spalania paliwa jest stała, a jej

wektor jest

styczny do trajektorii lotu. Wyznaczyć prędkość rakiety po

zmniejszeniu się jej masy do m oraz równanie jej ruchu s = s(t).

background image

Rozwiązanie. Ponieważ rakieta porusza się w przestrzeni
międzyplanetarnej, siły zewnętrzne na nią działające można pominąć, zatem W =
0, a dynamiczne równanie ruchu rakiety na podstawie (7.54) po uwzględnieniu
(7.55) można zapisać w postaci:

dt

dm

dt

d

m

w

C

v

v =

lub

m

dm

d

w

C

=

v

v

, lub

m

dm

d

w

C

v

v

=

.

(a)

Po scałkowaniu tego równania w granicach wyznaczonych przez warunki
początkowe, czyli dla t = 0 v

C

(0) = v

C0

i m(0) = m

0

, otrzymujemy:

=

m

m

w

v

v

C

o

co

m

dm

d

v

v

,

a po obliczeniu całek

0

0

C

C

m

m

ln

w

v

v

v

+

=

.

(b)

Ponieważ wektory prędkości v

C

i v

w

działają wzdłuż jednej prostej i mają

zwroty przeciwne (rys. 7.13), wektorowy wzór (b) można zapisać jednym wzorem
skalarnym:

v

v

v

C

C

w

=

0

ln

m

m

0

.

(c)

Powyższy wzór został po raz pierwszy wyprowadzony przez rosyjskiego
uczonego polskiego pochodzenia K. Ciołkowskiego.

Wektorowy wzór (b) lub równoważny mu (c) przedstawia prawo zmiany

prędkości rakiety. Ze wzorów tych wynika, że prędkość rakiety zależy od stosunku
masy końcowej rakiety m do jej masy początkowej m

0

.

Teraz wyznaczymy równanie drogi rakiety w funkcji czasu. Podstawiwszy do

wzoru (c)

dt

ds

v

C

=

,

otrzymujemy równanie różniczkowe o postaci:

dt

m

m

ln

v

dt

v

ds

0

w

0

C

=

.

Po scałkowaniu tego równania w granicach od s

0

do s i od 0 do t otrzymujemy

równanie ruchu rakiety:

background image

dt

m

m

ln

v

t

v

s

s

t

0

C

w

0

C

0

+

=

.

(d)

Aby

obliczyć występującą w tym równaniu całkę, należy znać funkcję

zmiany masy w czasie. Załóżmy, że w czasie pracy silnika rakiety jej masa maleje
wykładniczo według wzoru:

t

0

e

m

m

α

=

,

gdzie D jest stałym współczynnikiem. W tym przypadku

2

t

0

t

0

t

t

0

0

t

2

1

tdt

lne

dt

m

m

ln

α

=

=

=

α

.

Po podstawieniu otrzymanego wyniku do wzoru (d) otrzymujemy równanie ruchu
rakiety w funkcji czasu:

2

w

0

C

0

t

v

2

1

t

v

s

s

α

+

+

=

.

(e)

background image

7.3.1. Definicja krętu i kręt układu materialnego


Krętem k

O

punktu materialnego o masie m względem punktu O nazywamy

moment pędu

v

p m

=

tego punktu materialnego względem punktu O:

v

r

p

r

k

m

O

×

=

×

=

. (7.56)

Z

powyższej definicji wynika, że kręt

− zdefiniowany podobnie jak moment siły

względem punktu

− jest wektorem

prostopadłym do płaszczyzny
wyznaczonej przez punkt O i wektor
prędkości v (rys. 7.16).
Kręt punktu będzie równy zeru,
poza przypadkami trywialnymi (r = 0 i
v = 0), gdy wektory r i v będą
współliniowe.
Jeżeli będziemy mieli układ n
punktów materialnych o masach m

k

opisanych wektorami wodzącymi r

k

i

poruszających się z prędkością v

k

(rys. 7.17), to kręt tego układu materialnego

względem nieruchomego punktu O będzie równy sumie krętów (sumie momentów
pędów) nieruchomego punktu O będzie równy sumie krętów (sumie momentów
pędów)

k

o

m

mv

O

r

Rys. 7.16. Kręt (moment pędu) punktu

materialnego

=

=

×

=

×

=

n

1

k

k

k

k

n

1

k

k

k

O

m v

r

p

r

k

. (7.57)

background image

7.3.2. Redukcja krętu do środka masy


Wzór (7.57) opisuje kręt układu materialnego obliczony względem dowolnego
nieruchomego punktu O. Zadajmy sobie pytanie, jaki będzie kręt tego samego
układu materialnego względem środka masy C. W tym celu przyjmijmy w środku
masy C początek ruchomego układu współrzędnych o osiach

równoległych do odpowiednich osi nieruchomego układu współrzędnych x, y, z
(rys. 7.17). W tej sytuacji układ

′ ′ ′

x , y , z

′ ′ ′

x , y , z będzie się poruszał ruchem postępowym

względem układu nieruchomego x, y, z z prędkością środka masy v

C

.

v

1

v

2

r

C

r

Ck

m

k

z

x

z

y

y

x

r

k

C

O

m

1

v

k

m

2

m

n

v

n

v

C

v

C

v

Ck

Rys. 7.17. Rozkład prędkości układu punktów materialnych

Przy

takim

założeniu prędkość bezwzględna v

k

każdego punktu materialnego

względem układu nieruchomego x, y, z będzie sumą prędkości unoszenia równej
prędkości środka masy v

C

i prędkości względnej v

Ck

wzgędem układu ruchomego

, nazywanej dalej prędkością względem środka masy:

′ ′ ′

x , y , z

Ck

C

k

v

v

v

+

=

. (a)

Kręt rozpatrywanego układu punktów materialnych względem środka masy wyrazi
wzór:

=

×

=

n

1

k

k

Ck

C

m

k

v

r

k

, (7.58)

gdzie r

Ck

jest promieniem wodzącym punkt materialny o masie m

k

w układzie

. Z rysunku 7.17 wynika, że promień wodzący r

′ ′ ′

x , y , z

k

jest równy sumie

promienia wodzącego środka masy r

C

i promienia r

Ck

:

Ck

C

k

r

r

r

+

=

.

background image

Po wyznaczeniu z tej zależności

C

k

Ck

r

r

r

=

i podstawieniu do wzoru (7.58) otrzymamy:

(

)

=

=

=

×

×

=

×

=

n

1

k

n

1

k

k

C

n

1

k

k

k

C

C

m

m

m

k

k

k

k

k

v

r

v

r

v

r

r

k

.

(b)

Pierwsza suma po prawej stronie tego wzoru, zgodnie ze wzorem (7.57), jest
krętem k

O

względem nieruchomego punktu O, druga zaś jest pędem omawianego

układu materialnego. Na podstawie wzoru (7.42) możemy zapisać:

C

n

1

k

k

m

m

v

v

p

k

=

=

=

,

gdzie m jest masą całego układu. Zatem równanie (b) przyjmie postać:

C

C

O

C

mv

r

k

k

×

=

lub

C

C

C

O

mv

r

k

k

×

+

=

. (7.59)

Kręt k

O

układu punktów materialnych względem dowolnego nieruchomego

punktu O jest równy krętowi k

C

tego układu względem środka masy

powiększonemu o kręt

masy całkowitej skupionej w środku masy.

C

C

mv

r

×

Wzór (7.58) przedstawia kręt układu materialnego względem środka masy
obliczony dla ruchu bezwzględnego, ponieważ występująca w tym wzorze
prędkość v

k

jest prędkością względem nieruchomego układu odniesienia.

Zastanówmy się, czemu będzie równy kręt tego układu materialnego względem
środka masy wyznaczony dla ruchu względnego. W tym celu podstawmy do wzoru
(7.58) zależność (a).

(

)

.

m

m

m

m

m

m

m

m

n

1

k

Ck

k

Ck

n

1

k

k

Ck

C

n

1

k

Ck

k

Ck

C

n

1

k

k

Ck

n

1

k

n

1

k

Ck

k

Ck

C

k

Ck

n

1

k

Ck

C

k

Ck

n

1

k

k

k

Ck

C

=

=

=

=

=

=

=

=

×

+

×

=

×

+

×

⎟⎟

⎜⎜

=

=

×

+

×

=

+

×

=

×

=

v

r

r

v

v

r

v

r

v

r

v

r

v

v

r

v

r

k

Ale suma

=

=

n

1

k

k

Ck

0

m

r

,

background image

ponieważ moment statyczny układu względem środka masy jest równy zeru.
Ostatecznie mamy:

=

=

×

=

×

=

n

1

k

Ck

k

Ck

n

1

k

k

k

Ck

C

m

m

v

r

v

r

k

.

(7.60)

Z otrzymanej zależności wynika stwierdzenie:
Kręt układu punktów materialnych względem środka masy wyznaczony dla
ruchu bezwzględnego jest równy krętowi względem środka masy wyznaczonemu
dla ruchu względnego.

background image

7.3.3. Kręt bryły

Wyznaczmy kręt bryły o masie m poruszającej się ruchem dowolnym, a więc

bryły swobodnej. Podobnie jak w kinematyce bryły (p. 5.3.2) przyjmiemy dwa
układy współrzędnych

− jeden nieruchomy o początku w nieruchomym punkcie

O i osiach x, y, z, a drugi ruchomy, sztywno związany z bryłą o osiach

(rys. 7.18) i początku nie w dowolnym punkcie

′ ′ ′

x , y , z

O , lecz w środku masy C. W

bryle wydzielmy myślowo element masy dm o wektorze wodzącym

r

r

r

+

=

C

, (c)

gdzie

.

z

y

x

,

z

y

x

C

C

C

C

k

j

i

r

k

j

i

r

+

+

=

+

+

=

Znając prędkość v

C

środka masy C i prędkość kątową

ω, możemy obliczyć

prędkość v dowolnego punktu bryły (wzór 5.32). Zatem prędkość elementarnej
masy dm

r

ω

v

v

×

+

=

C

. (d)

Zgodnie z definicją kręt elementu
masy dm względem nieruchomego
punktu O

d

dm

O

k

r v

= ×

×

=

m

O

dm

v

r

k

.

Kręt bryły będzie równy całce z
powyższej zależności rozciągniętej
na całą masę m bryły:

.

Po podstawieniu do tego wzoru
zależności (c) i (d) otrzymamy:

(

) (

)

(

)

(

)

.

dm

dm

dm

dm

dm

m

C

m

C

m

m

C

C

C

C

×

ω

×

+

×

+

×

ω

×

+

×

=

×

ω

+

×

+

=

r

r

v

r

r

r

v

r

r

v

r

r

m

O

+

k

v

x

z

x

z

y

y

r

C

r

r

dm

C

O

Rys. 7.18. Opis położenia dowolnego elementu

bryły sztywnej

background image

Występujące pod całkami wielkości r

C

, v

C

i

ω nie podlegają całkowaniu i mogą

być wyciągnięte przed znaki całek:

(

)

.

dm

dm

dm

dm

m

m

C

m

C

m

C

C

O

×

×

+

×

⎟⎟

⎜⎜

×

×

+

×

=

r

ω

r

r

v

r

ω

r

v

r

k

Dwie środkowe całki są momentami statycznymi bryły względem środka masy,
a więc są równe zeru:

,

0

dm

m

=

r

a pierwsza całka jest masą całkowitą bryły:

.

=

m

dm

m

Ostatecznie kręt bryły możemy zapisać w postaci:

(

)

C

C

m

O

m

dm

v

r

r

ω

r

k

×

+

×

×

=

. (7.61)

Całka występująca w tym wzorze jest krętem

bryły w jej ruchu względem

środka masy C z prędkością kątową

ω.

k

C

(

)

×

×

=

m

C

dm

r

ω

r

k

. (7.62)

Zatem wzór (7.61) możemy zapisać w postaci:

C

C

C

O

m v

r

k

k

×

+

=

. (7.63)

Kręt k

O

bryły względem dowolnego nieruchomego punktu O jest równy krętowi

k

C

bryły względem środka masy C (w jej ruchu względem środka masy z

prędkością kątową

ω) powiększonemu o kręt r

v

C

m

C

×

masy m bryły

poruszającej się z prędkością v

C

środka masy.

Obecnie obliczymy współrzędne wektora k

C

w ruchomym układzie

współrzędnych

o początku w środku masy C (rys. 7.18). W tym układzie

współrzędnych wektory występujące we wzorze (7.62) mają następujące
współrzędne:

′ ′ ′

x , y , z

,

z

y

x

,

k

k

k

z

C

y

C

x

C

C

k

j

i

r

k

j

i

k

+

+

=

+

+

=

background image

=

ω

.

z

y

x

k

j

i

ω

+

ω

+

ω

Po rozpisaniu podwójnego iloczynu wektorowego ze wzoru (7.62), zgodnie ze
wzorem (2.34) otrzymamy:

(

)

( )

(

)

.

dm

dm

r

dm

dm

m

2

m

m

C

⋅′

=

⋅′

⋅′

=

m

r

ω

r

ω

ω

r

r

r

r

ω

k

Pierwsza całka występująca po prawej stronie powyższego równania jest
biegunowym momentem bezwładności względem środka masy C:

( )

=

m

2

C

dm

r

I

,

a więc

(

)

dm

I

m

C

C

⋅′

=

r

ω

r

ω

k

. (7.64)

Współrzędne krętu k

C

otrzymamy po zrzutowaniu tego wektora na osie

:

′ ′ ′

x , y , z

ω

=

=

C

x

C

x

C

I

k

i

k

(

)

,

dm

x

m

⋅′ ω

r

ω

=

=

C

y

C

y

C

I

k

j

k

(

)

,

dm

y

m

⋅′ ω

r

ω

=

=

C

z

C

z

C

I

k

k

k

(

)

.

dm

z

m

⋅′ ω

r

Po podstawieniu do tych wzorów iloczynu skalarnego:

=

⋅′ω

r

z

y

x

z

y

x

ω

+

ω

+

ω

oraz wyłączeniu przed całki współrzędnych prędkości kątowej otrzymujemy:

( )

( )

( )

.

dm

z

dm

z

y

dm

x

z

I

k

,

dm

z

y

dm

y

dm

y

x

I

k

,

dm

x

z

dm

y

x

dm

x

I

k

2

z

m

y

m

x

C

z

z

C

m

z

m

2

y

m

x

C

y

y

C

m

z

m

y

m

2

x

C

x

x

C

ω

ω

ω

ω

=

ω

ω

ω

ω

=

ω

ω

ω

ω

=

Całki występujące w powyższych wzorach są zdefiniowanymi w p. 6.1.2
momentami bezwładności bryły względem odpowiednich płaszczyzn i momentami
dewiacyjnymi. Po wykorzystaniu zależności (6.7) i (6.9) między momentami

background image

bezwładności względem bieguna, płaszczyzn i osi oraz odpowiednim
uporządkowaniu wyrazów współrzędne krętu k

C

bryły opisują wzory:

ω

+

ω

ω

=

ω

ω

+

ω

=

ω

ω

ω

=

.

I

D

D

k

,

D

I

D

k

,

D

D

I

k

z

z

z

y

y

x

z

x

z

C

z

y

z

y

y

y

x

x

y

C

x

z

z

y

x

y

x

x

x

C

(7.65)

Z

powyższych wzorów wynika, że do obliczenia krętu k

C

bryły swobodnej

względem środka masy C musimy znać wszystkie osiowe momenty bezwładności
i wszyskie momenty dewiacyjne, czyli tensor bezwładności. Wzory (7.65) znacznie
się upraszczają, gdy osie

′ ′ ′

x , y , z są głównymi centralnymi osiami bezwładności.

W tym przypadku, jak wiadomo z p. 6.5, wszystkie momenty dewiacyjne są równe
zeru i kręt

k

j

i

k

ω

+

ω

+

ω

=

z

z

y

y

x

x

C

I

I

I

. (7.66)

Jeżeli założymy, że osią obrotu bryły jest np. oś

z′ , to prędkość kątowa

ω pokryje

się z osią obrotu:

ω =

k

k

ω

=

ω

z

.

Wówczas kręt wyznaczony ze wzorów (7.65) ma postać:

k

j

i

k

ω

+

ω

ω

=

z

z

y

x

z

C

I

D

D

,

(7.67)

a na podstawie wzoru (7.66)

k

k

ω

=

z

z

C

I

. (7.68)

Z porównania wzorów (7.67) i (7.68) wynika, że jeżeli oś obrotu jest główną
centralną osią bezwładności, to wektor krętu leży na tej osi; gdy tak nie jest,
kierunek wektora krętu nie pokrywa się z osią obrotu.

Przykład 7.9. Korba OA o masie m

m

1

=

obraca się z prędkością kątową

ω

0

wokół osi z przechodzącej przez punkt O i prostopadłej do płaszczyzny rys. 7.19.

Na końcu A korby jest osadzona cienka
jednorodna tarcza o masie

i

promieniu r, która toczy się bez poślizgu
po nieruchomym kole o promieniu R.
Wyznaczyć kręt układu względem osi z.
Korbę OA uważać za pręt jednorodny.

m

2

2

=

ω

2

O

A

r

ω

0

R

C

v

A

Rys. 7.19. Wyznaczenie krętu układu

m



background image

Rozwiązanie

. Kręt układu względem osi z składa się z krętu

korby OA

poruszającej się ruchem obrotowym wokół osi z oraz krętu

tarczy poruszającej

się ruchem postępowym środka ciężkości A tarczy z prędkością

oraz ruchem

obrotowym z prędkością

względem osi

k

1z

k

z

2

v

A

ω

2

z równoległej do osi z i

przechodzącej przez środek tarczy:

z

2

z

1

z

k

k

k

+

=

. (a)

Kręt korby OA względem osi z

0

z

z

1

I

k

ω

=

. (b)

Kręt tarczy względem tej samej osi na podstawie wzoru (7.63) możemy wyrazić
zależnością:

(

)

A

2

2

z

z

2

v

m

r

R

I

k

+

+

ω

=

. (c)

We wzorach (b) i (c)

I i

są odpowiednio momentami bezwładności korby

względem osi z przechodzącej przez punkt O i tarczy względem osi
przechodzącej przez jej środek A. Zgodnie ze wzorami (f) i (a) z przykładu 6.2:

I

z

z

z

(

)

(

)

2

2

2

z

2

2

1

z

r

m

r

m

2

1

I

,

r

R

m

3

1

r

R

m

3

1

I

=

=

+

=

+

=

. (d)

Prędkość środka tarczy

(

)

0

A

r

R

v

ω

+

=

. (e)

Ponieważ punkt C (rys. 7.19) styku tarczy z nieruchomym kołem jest chwilowym
środkiem obrotu tarczy, mamy również:

,

r

v

2

A

ω

=

stąd

(

)

0

A

2

r

r

R

r

v

ω

+

=

=

ω

.

(f)

Po uwzględnieniu w związkach (b) i (c) wzorów (d), (e) i (f) oraz po ich
podstawieniu do równania (a) otrzymujemy kręt układu względem osi z.

(

)

(

)

(

)(

)

(

)(

)

.

r

10

R

7

r

R

m

3

1

r

R

r

R

m

2

r

r

R

r

m

r

R

m

3

1

k

0

0

0

2

0

2

z

ω

+

+

=

=

ω

+

+

+

ω

+

+

ω

+

=

background image

7.3.4. Zasada krętu i pokrętu. Zasada zachowania krętu


Załóżmy, że mamy układ materialny składający się z n punktów materialnych
o masach m

k

poruszających się z prędkością v

k

(rys. 7.17). Na każdy punkt niech

działa siła zewnętrzna P

k

oraz siły wewnętrzne F

kl

. Zgodnie z drugim prawem

Newtona możemy dla dowolnego punktu rozważanego układu materialnego
napisać dynamiczne równanie ruchu:

wk

P

P

r

+

=

k

2

k

2

k

dt

d

m

lub

(

)

n

,

2

,

1

k

t

d

d

m

k

k

k

,

.

.

.

P

P

v

wk

=

+

=

W powyższym równaniu zgodnie ze wzorem (7.45) P

wk

jest wypadkową sił

wewnętrznych działających na punkt o masie m

k

. Pomnóżmy wektorowo każde z n

równań obustronnie przez wektor wodzący r

k

i dodajmy wszystkie równania

stronami. Otrzymamy:

(

)

=

=

=

×

+

×

=

+

×

=

×

n

1

k

wk

k

k

n

1

=

k

k

n

1

k

n

1

k

wk

k

k

k

k

k

t

d

d

m

P

r

P

r

P

P

r

v

r

. (e)

Druga suma po prawej stronie tego równania jest sumą momentów sił
wewnętrznych względem punktu O i jak wykazano w p. 7.1.4 (wzór 7.13), jest
równa zeru. Z kolei suma momentów sił zewnętrznych względem punktu O jest
równa momentowi głównemu (3.26):

k

n

1

=

k

k

o

P

r

M

×

=

.

Sumę występującą po lewej stronie równania (e) możemy przekształcić:

(

)

(

)

.

dt

d

m

dt

d

m

dt

d

dt

d

m

dt

d

m

dt

d

m

O

n

1

k

n

1

k

k

k

k

k

k

k

k

k

n

1

k

k

n

1

k

n

1

=

k

k

k

k

k

k

k

k

k

k

v

r

v

r

v

r

v

r

v

v

v

r

=

×

=

×

=

=

×

=

×

+

×

=

×

=

=

=

=

Wynika z tego, że lewa strona równania (e) jest pochodną krętu całego układu
materialnego względem nieruchomego punktu O. Ostatecznie otrzymujemy:

background image

O

O

dt

d

M

k =

. (7.69)

Otrzymana zależność różniczkowa jest zasadą krętu.

Pochodna

względem czasu krętu układu punktów materialnych względem

dowolnego nieruchomego punktu jest równa momentowi głównemu wszystkich sił
zewnętrznych względem tego samego punktu.

Po obustronnym scałkowaniu równania (7.69) w granicach od 0 do t

otrzymamy:

( )

( )

=

t

0

O

O

O

dt

0

t

M

k

k

. (7.70)

Całka występująca w tym równaniu nosi nazwę pokrętu momentu głównego, a
samo równanie jest zasadą krętu i pokrętu.
Przyrost

krętu układu materialnego względem dowolnego nieruchomego punktu

jest równy pokrętowi momentu głównego sił zewnętrznych względem tego samego
punktu.

Równania (7.69) i (7.70) są słuszne nie tylko dla układu punktów materialnych,
ale i dla bryły.
Często się zdarza, że moment główny układu sił zewnętrznych względem
obranego nieruchomego bieguna redukcji O jest stale równy zeru bądź jest
pomijalnie mały,

. Wtedy całka po prawej stronie równania (7.70) jest

równa zeru i zasada krętu i pokrętu przechodzi w zasadę zachowania krętu:

0

O

M

( )

( )

( )

( )

const

0

t

czyli

0

0

t

O

O

O

O

=

=

=

k

k

,

k

k

lub

const

to

0

li

jeże

O

O

=

=

k

,

M

. (7.71)

Otrzymaną zasadę zachowania krętu można wyrazić słownie:

Jeżeli moment główny sił zewnętrznych względem nieruchomego punktu
redukcji O jest równy zeru, to kręt układu materialnego (bryły) względem tego
punktu jest wielkością stałą.


background image

7.3.5. Redukcja zasady krętu i pokrętu do środka masy


Zastanówmy

się, jaką postać przyjmie zasada krętu i pokrętu (7.70), jeżeli za

biegun redukcji przyjmiemy nie dowolny punkt O, lecz środek masy układu
materialnego C. W celu udzielenia odpowiedzi na postawione pytanie podstawmy
do równania (7.69) wzór (7.59):

C

C

C

O

mv

r

k

k

×

+

=

oraz twierdzenie o momencie głównym (3.29):

W

r

M

M

×

+

=

C

C

O

i dokonajmy różniczkowania:

(

)

W

r

M

v

r

k

×

+

=

×

+

C

C

C

C

C

dt

m

d

dt

d

,

(

)

W

r

M

v

r

v

r

k

×

+

=

×

+

×

+

C

C

C

C

C

C

C

dt

m

d

m

dt

d

dt

d

.

(f)

Drugi wyraz po lewej stronie powyższego równania jest równy zeru, ponieważ jest
to iloczyn wektorowy wektorów równoległych:

0

m

m

dt

d

C

C

C

C

=

×

=

×

v

v

v

r

,

a pochodna występująca w trzecim wyrazie jest pochodną względem czasu pędu
układu materialnego, równą wektorowi głównemu układu sił zewnętrznych (7.48):

(

)

W

p

v

=

=

dt

d

dt

m

d

C

.

Po uwzględnieniu powyższych zależności w równaniu (f) i uproszczeniu
otrzymamy zasadę krętu przy redukcji do środka masy:

C

C

dt

d

M

k

=

. (7.72)

Z kolei po scałkowaniu tego równania od zera do t otrzymamy zasadę krętu
i pokrętu zredukowaną do środka masy układu:

( )

( )

=

t

0

C

C

C

dt

0

t

M

k

k

. (7.73)

background image

Widzimy,

że formalna postać otrzymanych równań (7.72) i (7.73) jest taka

sama jak równań (7.69) i (7.70), ale równania (7.72) i (7.73) nie opisują ruchu
środka masy C. Do opisu ruchu środka masy C należałoby zastosować zasadę pędu
(7.48).
Jeżeli założymy teraz, że moment sił zewnętrznych względem środka masy C
układu materialnego będzie stale równy zeru,

M

C

≡ 0 , to zasada krętu i pokrętu

(7.73) zredukowana do środka masy przejdzie w zasadę zachowania krętu
względem środka masy, co można zapisać w następujący sposób:

const

to

0

li

jeże

C

C

=

=

k

,

M

(7.74)

lub ująć słownie:

Jeżeli moment główny sił zewnętrznych względem środka masy układu
materialnego jest równy zeru, to kręt tego układu materialnego względem środka
masy jest wielkością stałą.

Przykład 7.10. Punkt materialny A o masie m

1

zaczął się poruszać wzdłuż

cięciwy BC (rys. 7.20a) poziomej jednorodnej tarczy kołowej o promieniu R i
masie m według równania:

sinkt

b

x

=

,

gdzie x oznacza współrzędną odmierzoną jak na rys. 7.20, k pewną stałą, a

2b BC

. Tarcza może się obracać bez tarcia wokół osi pionowej z przechodzącej

przez środek tarczy O. Wyznaczyć prędkość kątową

ω tarczy w funkcji czasu t,

jeżeli odległość cięciwy od środka tarczy wynosi b, a tarcza w chwili początkowej

była nieruchoma.

t

= 0

O

A

R

O

A

R

v

w

r

ω

b

A

0

α

x

x

b

v

u

A

0

a)

b)

α

B

C

Rys. 7.20. Wyznaczenie prędkości kątowej tarczy

Rozwiązanie. Na układ działają siły zewnętrzne ciężkości tarczy i punktu
materialnego oraz reakcje w łożyskach osi obrotu tarczy. Siły ciężkości są
równoległe do osi obrotu, więc ich momenty względem osi obrotu są zawsze

background image

równe zeru. Nie dają momentu względem tej osi również reakcje w łożyskach.
Zatem zgodnie z zasadą zachowania krętu (7.71) kręt układu względem osi nie
ulega zmianie. Ponieważ w chwili początkowej t

= 0 , gdy punkt A był jeszcze

nieruchomy, kręt układu był równy zeru, zatem w dowolnej chwili t kręt tego
układu również będzie równy zeru. Po rozpoczęciu ruchu punktu A tarcza zacznie
się poruszać ruchem obrotowym z prędkością kątową

do ruchu punktu (rys. 7.20b). Prędkość punktu tarczy, w którym w chwili t
znajduje się punkt A, czyli prędkość unoszenia punktu A

kt

sin

1

b

x

b

r

v

2

2

2

u

+

ω

=

+

ω

=

ω

=

.

Prędkość punktu A względem tarczy (prędkość względna)

coskt

bk

dt

dx

v

w

=

=

.

Z kolei prędkość bezwzględna punktu A jest równa sumie wektorowej prędkości
unoszenia i prędkości względnej:

w

u

A

v

v

v

+

=

.

Rzut wektora prędkości bezwzględnej punktu A na kierunek prostopadły do
promienia OA r

= jest równy

u

w

v

cos

v

α

.

Kręt układu w chwili t względem osi obrotu z składa się z krętu

punktu A

i krętu

tarczy względem tej osi. Kręt punktu A

k

1z

k

z

2

(

)

(

)

(

) (

)

(

)

[

]

,

kt

sin

1

b

coskt

k

b

m

x

b

kt

sin

1

b

coskt

k

b

m

r

kt

sin

1

b

b

v

m

r

v

cos

r

v

m

v

cos

v

r

m

k

2

2

2

1

2

2

2

2

1

2

w

1

u

w

1

u

w

1

z

1

+

ω

=

=

+

+

ω

=

+

ω

=

=

α

=

α

=

a kręt tarczy względem osi obrotu

ω

=

ω

=

2

z

z

2

mR

2

1

I

k

.

Ponieważ kręt całkowity układu jest w każdej chwili równy zeru, otrzymujemy:

(

)

[

]

0

mR

2

1

kt

sin

1

b

coskt

k

b

m

2

2

2

2

1

=

ω

+

ω

.

background image

Z powyższego równania znajdujemy prędkość kątową tarczy:

(

)

2

2

2

1

2

1

mR

2

1

kt

sin

1

b

m

coskt

k

b

m

+

+

=

ω

.

background image

7.4.1. Energia kinetyczna układu punktów materialnych


Energią kinetyczną punktu materialnego o masie m, poruszającego się z
prędkością v, nazywamy połowę iloczynu masy punktu i kwadratu jego prędkości:

2

mv

E

2

=

.

Dla układu n punktów materialnych o masach m

k

poruszających się

z prędkością v

k

energia kinetyczna będzie równa sumie energii kinetycznych

poszczególnych punktów materialnych:

=

=

=

=

n

1

k

2
k

k

n

1

k

2
k

k

v

m

2

1

2

v

m

E

. (7.75)

Podobnie jak w przypadku krętu układu punktów materialnych (7.3.2),
prędkość bezwzględną v

k

każdego punktu materialnego rozłożymy na prędkość

unoszenia v

C

, wywołaną ruchem postępowym ruchomego układu współrzędnych

o początku w środku masy C względem układu nieruchomego x, y, z,

i prędkość względną v

′ ′ ′

x , y , z

Ck

względem układu ruchomego (rys. 7.17):

Ck

C

k

v

v

v

+

=

.

Po podstawieniu tej zależności do wzoru (7.75) oraz przedstawieniu kwadratu
prędkości w postaci iloczynu skalarnego

k

k

2
k

v

v

v

=

otrzymamy:

(

) (

)

(

)

=

+

+

=

=

+

+

=

=

=

=

=

n

1

k

2
Ck

Ck

C

2
C

k

Ck

C

n

1

k

Ck

C

k

n

1

k

k

k

k

v

2

v

m

2

1

m

2

1

m

2

1

E

v

v

v

v

v

v

v

v

=

=

=

+

+

=

n

1

k

2
Ck

k

1

Ck

k

C

n

1

k

k

2
C

v

m

2

1

m

m

v

2

1

n

k

v

v

. (a)

Drugi wyraz po prawej stronie powyższego równania jest równy zeru, ponieważ
występująca w nim suma jest pędem układu punktów materialnych w jego ruchu
względem ruchomego układu współrzędnych

′ ′ ′

x , y , z . Wiadomo jednakże, że

pęd jest równy iloczynowi masy całkowitej i prędkości środka masy (7.44), która w
stosunku do ruchomego układu odniesienia

′ ′ ′

x , y , z jest równa zeru. Zatem

background image

0

m

n

1

k

Ck

k

=

=

v

.

Ostatni wyraz jest energią kinetyczną układu punktów materialnych w jego ruchu
względem ruchomego układu odniesienia

′ ′ ′

x , y , z :

=

=

n

1

k

2
Ck

k

c

v

m

2

1

E

. (7.76)

Po oznaczeniu masy całkowitej rozpatrywanego układu materialnego przez

=

=

n

1

k

k

m

m

równanie (a) przyjmuje postać:

2
C

C

mv

2

1

E

E

+

=

. (7.77)

Zależność (7.77) nosi nazwę twierdzenia Koeniga.

Energia kinetyczna układu punktów materialnych jest równa energii tegoż
układu w jego ruchu względem środka masy oraz energii kinetycznej masy
całkowitej poruszającej się z prędkością środka masy.

background image

7.4.2. Energia kinetyczna bryły


W celu wyznaczenia energii kinetycznej bryły o masie m poruszającej się
ruchem ogólnym postąpimy podobnie jak przy wyznaczaniu krętu bryły (p. 7.3.3).
W bryle myślowo wydzielimy element masy dm (rys. 7.18) poruszający się z
prędkością zgodną ze wzorem (5.32):

r

ω

v

v

×

+

=

C

. (b)

Energia kinetyczna tego elementu

dm

2

1

dE

v

v

=

,

a energia bryły jest równa całce względem całej masy z tego wyrażenia:

=

m

dm

2

1

E

v

v

. (c)

Po podstawieniu do wzoru (c) prędkości w postaci (b) otrzymamy:

(

) (

)

=

×

+

×

+

=

m

C

C

dm

2

1

E

r

ω

v

r

ω

v

(

)

(

) (

)

dm

2

1

dm

dm

v

2

1

m

m

C

m

2
C

×

×

+

×

+

=

r

ω

r

ω

r

ω

v

.

(d)

Po przekształceniu wyrażeń podcałkowych w drugiej i trzeciej całce do postaci:

(

) (

)

(

) (

)

(

)

[

]

r

ω

r

ω

r

ω

r

ω

r

ω

v

r

ω

v

×

×

=

×

×

×

=

×

,

C

C

oraz wyłączeniu przed całki v

C

i

ω, jako wielkości niezależnych od zmiennych

całkowania

, wzór (d) możemy zapisać:

′ ′ ′

x , y , z

(

)

(

)

×

×

+

×

+

=

m

m

C

m

2
C

dm

2

1

dm

dm

v

2

1

E

r

ω

r

ω

r

ω

v

.

(e)

Pierwsza całka jest masą bryły, druga momentem statycznym względem środka
masy, a trzecia krętem bryły w ruchu względem środka masy (7.62), czyli

(

)

×

ω

×

=

=

=

m

C

m

m

dm

oraz

0

dm

,

dm

m

r

r

k

r

.

Po uwzględnieniu powyższych zależności we wzorze (e) otrzymujemy:

background image

2
C

C

mv

2

1

2

1

E

+

=

k

ω

. (7.78)

Pierwszy wyraz w powyższym wzorze jest energią kinetyczną bryły w jej
chwilowym ruchu obrotowym względem środka masy:

.

C

C

2

1

E

k

ω

=

(7.79)

Zatem energię kinetyczną bryły możemy przedstawić w postaci identycznej ze
wzorem (7.77):

E

E

C

=

+

1
2

mv

C

2

. (7.80)

Jest to twierdzenie Koeniga dla bryły.

Energia kinetyczna bryły w ruchu ogólnym jest sumą energii kinetycznej bryły w
jej chwilowym ruchu obrotowym względem środka masy i energii kinetycznej masy
całkowitej poruszającej się z prędkością środka masy.

Aby

obliczyć energię E

C

we wzorze (7.79), przedstawimy iloczyn skalarny za

pomocą współrzędnych wektorów

ω i k

C

danych w układzie ruchomym

:

′ ′ ′

x , y , z

C

C

2

1

E

k

ω

=

=

(

)

z

C

z

y

C

y

x

C

x

k

k

k

2

1

ω

+

ω

+

ω

.

Po podstawieniu w tym wzorze współrzędnych krętu danych wzorami (7.65)
i uporządkowaniu wyrazów energię kinetyczną bryły w jej ruchu względem środka
masy możemy przedstawić w postaci:

(

)

ω

+

ω

+

ω

=

2
z

z

2
y

y

2
x

x

C

I

I

I

2

1

E

(

)

x

z

x

z

z

y

z

y

y

x

y

x

D

D

D

ω

ω

+

ω

ω

+

ω

ω

(7.81)

Zatem, podobnie jak w przypadku krętu k

C

, do obliczenia energii kinetycznej

bryły w jej ruchu względem środka masy musimy znać wszystkie osiowe i
dewiacyjne momenty bezwładności.
Gdy

osie

′ ′ ′

x , y , z są głównymi centralnymi osiami bezwładności, momenty

dewiacyjne znikają, a wzór (7.81) upraszcza się do postaci:

(

)

2
z

z

2
y

y

2
x

x

C

I

I

I

2

1

E

ω

+

ω

+

ω

=

.

(7.82)

background image

Jeżeli ruch bryły jest ruchem obrotowym wokół stałej osi obrotu, np. l, z
prędkością kątową

ω, to energia ruchu obrotowego

2

l

I

2

1

E

ω

=

, (7.83)

gdzie I

l

jest momentem bezwładności

względem osi obrotu l.

Przykład 7.11. Kołowrót o masie
m

1

= 5m i promieniach r oraz R = 1,5r

toczy się bez poślizgu małym obwodem
po poziomej listwie (rys. 7.17). Środek
masy C tego kołowrotu znajduje się na
osi symetrii obrotowej i ma stałą
prędkość v

C

. Na duży obwód nawinięto

linkę, na której końcu zawieszono
ciężarek o masie m

2

= m. Promień

bezwładności kołowrotu względem osi symetrii prostopadłej do płaszczyzny
rysunku jest równy

. Obliczyć energię kinetyczną tego układu.

i

C

ω

v

2

v

A

v

C

v

A

A

C

R

r

S

v

C

m

2

Rys. 7.21. Wyznaczenie energii kinetycznej

kołowrotu


Rozwiązanie. Energia kinetyczna układu jest równa sumie energii kinetycznej
kołowrotu E

1

poruszającego się ruchem płaskim i energii kinetycznej ciężarka E

2

poruszającego się ruchem postępowym:

2

1

E

E

E

+

=

.

Wzór na energię kinetyczną kołowrotu, zgodnie z równaniem (7.80) wynikającym
z twierdzenia Koeniga, po uwzględnieniu zależności (7.83) ma postać:

C

1

2

C

1

v

m

2

1

I

2

1

E

+

ω

=

, (a)

gdzie moment bezwładności kołowrotu względem osi symetrii obrotowej

2
C

2
C

1

C

mi

5

i

m

I

=

=

. (b)

Energia kinetyczna ciężarka

2
2

2
2

2

2

mv

2

1

v

m

2

1

E

=

=

. (c)

background image

Ponieważ kołowrót toczy się bez poślizgu, chwilowy środek obrotu znajduje się w
punkcie S styku kołowrotu z listwą. Korzystając z własności chwilowego środka
obrotu, możemy napisać:

(

)

C

C

A

C

v

2

5

v

r

r

R

r

R

v

,

r

v

=

+

=

+

ω

=

=

ω

.

(d)

Zgodnie z rysunkiem prędkość ciężarka v

2

jest równa sumie geometrycznej

prędkości v

C

i v

A

. Stąd kwadrat prędkości v

2

2
C

2
C

2
A

2
2

v

4

29

v

v

v

=

+

=

. (e)

Po dodaniu wzoru (c) do (a) i uwzględnieniu zależności (b), (d) i (e) otrzymujemy
całkowitą energię kinetyczną układu:

2
C

2

C

2

2
C

2

C

2
C

mv

8

49

r

i

2

5

mv

8

29

mv

2

5

r

v

mi

2

5

E

+

=

+

+

=

.

background image

7.4.3. Zasada pracy i energii kinetycznej


Dla

każdego z n punktów materialnych układu omówionego w p. 7.2.2 i

przedstawionego na rys. 7.12 napiszemy, tak jak poprzednio, dynamiczne równanie
ruchu (7.47):

wk

k

2

k

2

k

dt

d

m

P

P

r

+

=

albo

(

)

n

,.

..

,

2

,

1

k

t

d

d

m

wk

k

k

k

=

+

=

P

P

v

.

Pomnóżmy skalarnie każde z tych równań przez prędkość v

k

i dodajmy je

stronami:

(

)

=

=

=

=

+

=

+

=

n

1

k

k

wk

n

1

k

k

k

k

n

1

k

wk

k

n

1

k

k

k

k

t

d

d

m

v

P

v

P

v

P

P

v

v

. (e)

Zgodnie z definicją podaną w p. 7.1.7 pierwsza suma w równaniu (e) jest mocą
układu sił zewnętrznych:

=

=

n

1

k

k

k

z

N

v

P

,

(7.84)

a druga podwójna suma mocą wszystkich sił wewnętrznych:

=

=

n

1

k

k

wk

w

N

v

P

. (7.85)

Wykażemy, że lewa strona równania (e) jest pochodną względem czasu energii
całkowitej układu punktów materialnych:

(

)

.

dt

dE

m

2

1

dt

d

dt

m

d

2

1

dt

d

m

dt

d

m

2

1

t

d

d

m

n

1

k

k

k

k

n

1

k

k

k

k

n

1

k

k

k

k

k

k

k

n

1

k

k

k

k

=

=

=

=

+

=

=

=

=

=

v

v

v

v

v

v

v

v

v

v

Ostatecznie równanie (e) przyjmuje postać:

w

z

N

N

dt

dE

+

=

.

(7.86)

background image

Zatem pochodna względem czasu energii kinetycznej układu materialnego jest
równa sumie mocy wszystkich sił zewnętrznych i wewnętrznych. Po scałkowaniu
obustronnie równania (7.86) od 0 do t otrzymamy:

( ) ( )

+

=

t

0

w

t

0

z

dt

N

dt

N

0

E

t

E

. (f)

Całki występujące w powyższym równaniu, zgodnie ze wzorem (7.36),
przedstawiają odpowiednio pracę sił zewnętrznych i wewnętrznych:

=

=

t

0

w

w

t

0

z

z

dt

N

L

,

dt

N

L

. (g)

Po wprowadzeniu oznaczeń (g) do równania (f) otrzymujemy zasadę pracy
i energii kinetycznej dla układu punktów materialnych:

( ) ( )

w

z

L

L

0

E

t

E

+

=

lub po wprowadzeniu oznaczeń E(t) = E

2

, E(0) = E

1

w

z

1

2

L

L

E

E

+

=

. (7.87)

Przyrost energii kinetycznej układu punktów materialnych w skończonym
przedziale czasu jest równy pracy wykonanej w tym samym czasie przez wszystkie
siły zewnętrzne i wewnętrzne.

Bez przeprowadzania dowodu metodą analityczną można zauważyć, że praca sił
wewnętrznych jest ściśle związana ze zmianą odległości między punktami układu
materialnego. Gdy odległości między punktami układu materialnego nie ulegają
zmianie, praca sił wewnętrznych będzie równa zeru. Zatem dla bryły sztywnej lub
ciała sztywnego praca sił wewnętrznych jest równa zeru, L

w

= 0. W tej sytuacji

zasadę pracy i energii kinetycznej dla bryły sztywnej można zapisać w postaci:

z

1

2

L

E

E

=

. (7.88)

Przyrost energii kinetycznej bryły sztywnej w skończonym przedziale czasu jest
równy pracy wykonanej w tym samym czasie przez wszystkie siły zewnętrzne
działające na tę bryłę.

Przykład 7.12. Do bębna kołowrotu o promieniu r i masie m

1

jest przyłożony

stały moment obrotowy M. Do końca wiotkiej liny nawiniętej na bęben
przymocowano ciężar o masie m

2

, który przesuwa się po równi pochyłej o kącie

nachylenia Dα(rys. 7.22). Współczynnik tarcia między masą m

2

a równią wynosi

µ.

Jaką prędkość kątową

ω osiągnie bęben po obróceniu się o ϕ radianów, jeżeli w

background image

chwili początkowej układ był w spoczynku? Masę liny pominąć, a bęben uważać
za jednorodny walec.

r

ϕ

v

2

ϕ,ω

M

N

T

α

G

2

O

r

Rys. 7.22. Wyznaczenie prędkości kątowej bębna

Rozwiązanie. Do rozwiązania zadania zastosujemy zasadę pracy i energii
kinetycznej (7.88):

L

E

E

1

2

=

.

Z uwagi na to, że układ w chwili początkowej znajdował się w spoczynku, jego
energia kinetyczna była równa zeru, E

1

= 0. Otrzymujemy więc:

L

E

2

=

. (a)

Energia kinetyczna układu składa się z energii kinetycznej ruchu postępowego
masy m

2

oraz ruchu obrotowego bębna:

2

O

2
2

2

2

I

2

1

v

m

2

1

E

ω

+

=

.

Ponieważ moment bezwładności bębna I

O

względem osi obrotu i prędkość v

2

równe:

r

v

,

r

m

2

1

I

2

2

1

O

ω

=

=

,

mamy:

(

)

2

2

2

1

2

2

1

2

2

2

2

r

m

2

m

4

1

r

m

4

1

r

m

2

1

E

ω

+

=

ω

+

ω

=

.

(b)

Pracę L wykonują: moment obrotowy M, składowa siły ciężkości G

2

równoległa

do równi oraz siła tarcia T. Jeżeli zauważymy, że przy obrocie bębna o kąt

ϕ ciężar

o masie m

2

przesunie się w górę równi o r

ϕ, możemy napisać:

background image

(

)

L M

m g

T r

=

+

ϕ

α

2

sin

ϕ .

Po podstawieniu do tego wzoru

α

=

=

cos

g

m

µ

N

µ

T

2

wykonana praca

(

ϕ

⎥⎦

⎢⎣

α

+

α

=

r

cos

µ

sin

g

m

r

M

L

2

)

.

(c)

Po podstawieniu zależności (b) i (c) do wzoru (a) otrzymujemy równanie:

(

)

(

)

ϕ

⎥⎦

⎢⎣

α

+

α

=

ω

+

r

cos

µ

sin

g

m

r

M

r

m

2

m

4

1

2

2

2

2

1

,

skąd

(

)

ϕ

+

α

α

=

ω

2

1

2

m

2

m

cos

µ

+

sin

r

g

m

M

r

2

.

background image

7.4.4. Zasada zachowania energii

Obecnie rozpatrzymy ruch układu materialnego, na który działają siły
potencjalne, zarówno zewnętrzne jak i wewnętrzne. W punkcie 7.1.5
udowodniono, że jeżeli na punkt materialny działa siła potencjalna, to praca
wykonana przez tę siłę jest równa ubytkowi energii potencjalnej. Przyjmiemy bez
dowodu, że zależność ta jest słuszna nie tylko dla każdego punktu, ale i dla całego
układu materialnego. Zatem pracę sił zewnętrznych i wewnętrznych możemy
zapisać w postaci:

=

=

,

U

U

L

,

U

U

L

2

w

1

w

w

2

z

1

z

z

(h)

gdzie U

z1

i U

z2

oznaczają energię potencjalną sił zewnętrznych w położeniu

początkowym i końcowym, a U

w1

i U

w2

energię potencjalną sił wewnętrznych

w położeniu początkowym i końcowym. Po podstawieniu wzorów (h) do równania
zasady pracy i energii kinetycznej (7.87) otrzymamy:

E

2

– E

1

= U

z1

– U

z2

+ U

w1

– U

w2

lub

E

2

+ U

z2

+ U

w2

= E

1

+ U

z1

+ U

w1

.

(i)

Z równania (i) wynika, że suma energii kinetycznej i energii potencjalnej sił
zewnętrznych i wewnętrznych jest w każdym położeniu układu wielkością stałą.
Po wprowadzeniu do równania (i) oznaczeń:

U

2

= U

z2

+ U

w2

i U

1

= U

z1

+ U

w1

otrzymamy:

E

2

+ U

2

= E

1

+ U

1

albo ogólnie

E + U = const. (7.89)

Jest to zasada zachowania energii mechanicznej.

Gdy na układ materialny działają siły potencjalne, wtedy suma energii
kinetycznej i potencjalnej tego układu jest wielkością stałą.

Zasada zachowania energii mechanicznej jest słuszna również w przypadku,
gdy działające siły można rozłożyć na siły potencjalne i siły, które nie są
potencjalne, ale nie wykonują pracy, np. reakcje gładkich powierzchni.
Układy materialne, do których odnosi się zasada zachowania energii
mechanicznej, nazywamy układami zachowawczymi, a siły siłami
zachowawczymi. Układy, których nie dotyczy ta zasada, nazywamy układami
rozpraszającymi lub dyssy-patywnymi, np. układy z tarciem.

background image

Zasada

zachowania

energii

mechanicznej jest trzecią zasadą zachowania

w dynamice, po zasadzie zachowania pędu i zasadzie zachowania krętu. Należy
pamiętać, że zasady zachowania są słuszne tylko wówczas, gdy są spełnione
odpowiednie założenia poczynione przy ich wyprowadzaniu.

Przykład 7.13. Cienki jednorodny pręt OA o długości L i masie m może się
obracać bez tarcia wokół osi poziomej prostopadłej do osi pręta przechodzącej

przez jego koniec O (rys. 7.23). Jaką
prędkość należy nadać końcowi A w
chwili, gdy pręt jest w spoczynku w
położeniu równowagi stałej, aby wykonał
on ćwierć obrotu?

L/2

L

O

ω

A

C

mg

v

A

U = 0

Rys. 7.23. Wyznaczenie prędkości

początkowej końca pręta

Rozwiązanie. Na pręt działa siła
ciężkości, która jest siłą potencjalną.
Zatem do rozwiązania zadania możemy
zastosować zasadę zachowania energii
mechanicznej (7.89):

2

2

1

1

U

E

U

E

+

=

+

.

(a)

Jeżeli poziom zerowej energii potencjalnej przyjmiemy na wysokości środka
ciężkości C, jak na rysunku, to

U

1

0

= . Po wykonaniu ćwierć obrotu pręt zajmie

położenie poziome i zatrzyma się. Jego energia kinetyczna będzie równa zeru,

. Równanie (a) będzie miało więc postać:

E

2

0

=

2

1

U

E

=

. (b)

W chwili początkowej energia kinetyczna

2

O

1

I

2

1

E

ω

=

.

Moment bezwładności pręta jednorodnego względem jego końca (patrz przy-
kład 6.2)

3

mL

I

2

O

=

.

Z kolei prędkość kątowa pręta

L

v

A

=

ω

.

background image

Energia kinetyczna pręta ma więc postać:

6

mv

L

v

3

mL

2

1

E

2
A

2

A

2

1

=

=

. (c)

Energia potencjalna pręta w położeniu końcowym

2

L

mg

U

2

=

. (d)

Po podstawieniu wzorów (c) i (d) do równości (b) otrzymujemy równanie:

2

mgL

6

mv

2
A

=

.

Stąd prędkość początkowa końca A pręta

L

g

3

v

A

=

.

Czytelnikowi pozostawiamy wyznaczenie prędkości, jaką należy nadać
końcowi A pręta, aby wykonał on pełen obrót.

background image

7.5.1. Ruch bryły swobodnej


Swobodna

bryła sztywna ma w przestrzeni sześć stopni swobody i do określenia

jej ruchu potrzeba sześciu równań ruchu. Ruch bryły możemy rozbić na ruch
środka masy, wywołany przez działanie wektora głównego sił zewnętrznych, i
obrót bryły względem środka masy, wywołany przez moment główny sił
zewnętrznych zredukowany do środka masy.
Do

ułożenia równań ruchu bryły wykorzystamy wyprowadzone poprzednio

zasady pędu i krętu. W punkcie 7.2.3 wykazano, że pochodna pędu względem
czasu równa wektorowi głównemu sił zewnętrznych opisuje ruch środka masy, a w
punkcie 7.3.5, że pochodna krętu zredukowanego do środka masy względem czasu
równa momentowi głównemu sił zewnętrznych opisuje obrót bryły względem
środka masy. Mamy więc dwa równania wektorowe opisujące ruch bryły
swobodnej:

C

C

dt

d

,

dt

d

M

k

W

p

=

=

.

(7.90)

Te dwa równania wektorowe są równoważne sześciu równaniom skalarnym.
Otrzymamy je po zrzutowaniu wektorów występujących w powyższych

równaniach na

osie prostokątnego

układu współrzędnych. Podobnie jak
przy obliczaniu krętu bryły
przyjmiemy dwa układy
współrzędnych: jeden nieruchomy x,
y, z o początku w dowolnym punkcie
O i drugi ruchomy

′ ′ ′

x , y , z sztywno

związany z bryłą o początku w
środku masy C (rys. 7.24). Ponadto
dla uproszczenia obliczeń założymy,
że osie

z

,

y

,

x

układu ruchomego są

głównymi centralnymi osiami
bezwładności. Przy takim założeniu
zgodnie ze wzorem (7.66) kręt bryły

k

j

i

k

ω

+

ω

+

ω

=

z

z

y

y

x

x

C

I

I

I

, (a)

gdzie

są głównymi centralnymi momentami bezwładności,

a

współrzędnymi wektora prędkości kątowej

ω w układzie

ruchomym.

z

y

x

I

,

I

,

I

ω ω ω

x

,

,

y

z

x

z

x

z

y

y

r

C

C

O

M

C

W

Rys. 7.24. Ruch swobodny bryły sztywnej

background image

W pierwszej kolejności obliczymy pochodną krętu k

C

względem czasu z

wykorzystaniem podanych w kinematyce bryły wzorów na pochodne względem
czasu wersorów układu ruchomego (5.31).

k

ω

k

j

ω

j

i

ω

i

×

=

×

=

×

=

t

d

d

,

t

d

d

,

t

d

d

.

+

ω

+

ω

+

ω

=

=

ω

+

ω

+

ω

+

+

ω

+

ω

+

ω

=

k

j

i

k

j

i

k

j

i

k

dt

d

I

dt

d

I

dt

d

I

dt

d

I

dt

d

I

dt

d

I

dt

d

I

dt

d

I

dt

d

I

dt

d

z

z

y

y

x

x

z

z

y

y

x

x

z

z

y

y

x

x

C

+

(

)

k

j

i

ω

ω

+

ω

+

ω

×

z

z

y

y

x

x

I

I

I

.

Wyrażenie w nawiasie w powyższym wzorze jest krętem bryły względem
środka masy. Zatem pochodna krętu k

C

względem czasu

C

z

z

y

y

x

x

C

dt

d

I

dt

d

I

dt

d

I

dt

d

k

ω

k

j

i

k

×

+

ω

+

ω

+

ω

=

. (7.91)

Po obliczeniu iloczynu wektorowego występującego w tym wzorze oraz
odpowiednim pogrupowaniu wyrazów otrzymamy ostatecznie:

(

)

(

)

(

)

.

I

I

dt

d

I

I

I

dt

d

I

I

I

dt

d

I

dt

d

y

x

x

y

z

z

z

x

z

x

y

y

z

y

y

z

x

x

C

k

j

i

k

⎥⎦

⎢⎣

ω

ω

+

ω

+

+

ω

ω

+

ω

+

+

⎥⎦

⎢⎣

ω

ω

+

ω

=

(7.92)

Po zapisaniu występującego w równaniach (7.90) wektora głównego W i
momentu głównego M

O

w ruchomym układzie współrzędnych:

k

j

i

M

k

j

i

W

+

+

=

+

+

=

z

C

y

C

x

C

C

z

y

x

M

M

M

,

W

W

W

oraz podstawieniu do drugiego równania (7.90) wzoru (7.92) i porównaniu
wyrażeń przy wersorach otrzymamy sześć skalarnych równań ruchu bryły:

background image

(

)

(

)

(

)

⎪⎪

=

ω

ω

+

ε

=

ω

ω

+

ε

=

ω

ω

+

ε

=

=

=

,

M

I

I

I

,

M

I

I

I

,

M

I

I

I

,

W

ma

,

W

ma

,

W

ma

z

C

y

x

x

y

z

z

y

C

z

x

z

x

y

y

x

C

z

y

y

z

x

x

z

z

C

x

x

C

y

y

C

(7.93)

w których zamiast pochodnych względem czasu współrzędnych prędkości kątowej
ω wprowadzono odpowiednie współrzędne przyśpieszenia kątowego ε:

dt

d

,

dt

d

,

dt

d

z

z

y

y

x

x

ω

=

ε

ω

=

ε

ω

=

ε

,

a

są współrzędnymi przyśpieszenia a

a

a

a

Cx

Cy

Cz

,

,

C

środka masy C.

Powyższe równania różniczkowe wraz z warunkami początkowymi
jednoznacznie opisują ruch bryły pod wpływem przyłożonego do niej układu sił.
Przy wyprowadzaniu równań ruchu bryły (7.93) za biegun redukcji przyjęto
środek masy C bryły. Początek ruchomego układu współrzędnych można przyjąć
poza środkiem masy, pod warunkiem że punkt ten jest nieruchomy. Jeżeli w
poruszającej się bryle istnieje nieruchomy punkt, np. O, to obierając go za biegun
redukcji, otrzymamy równania o postaci (7.93), ale wtedy zamiast współrzędnych

momentu głównego M

z

C

y

C

x

C

M

,

M

M

C

zredukowanego do środka masy C

należy podstawić współrzędne

momentu M

z

O

y

O

x

O

M

,

M

M

O

zredukowanego do

tego nieruchomego punktu. Występujące w tych równaniach momenty
bezwładności muszą być głównymi momentami bezwładności.

background image

7.5.2. Obrót bryły wokół stałej osi obrotu

Obrót dowolny bryły wokół głównej osi bezwładności

Ważnym zagadnieniem w dynamice maszyn jest ruch obrotowy bryły wokół

stałej osi obrotu. Z tym zagadnieniem mamy do czynienia we wszystkich
maszynach wirnikowych. Aby taki ruch można było zrealizować, bryła (wirnik)
musi być ograniczona więzami. Są nimi najczęściej łożyska, w których w czasie
ruchu bryły powstają odpowiednie reakcje.

Z kinematyki wiadomo, że obracająca się bryła wokół stałej osi obrotu ma jeden

stopień swobody. Taki ruch bryły można jednoznacznie opisać jednym równaniem
ruchu w postaci kąta obrotu w funkcji czasu

ϕ = ϕ(t).

x

x

y

y

z = z

O

C

ω

ε

r

c

l

M

o

W

ϕ

ϕ


Rys. 7.25. Ruch obrotowy bryły sztywnej wokół głównej osi bezwładności

Dla wyprowadzenia dynamicznego równania ruchu bryły założymy, że bryła
obraca się ruchem dowolnym, czyli że prędkość kątowa bryły nie jest stała, wokół
osi

będącej główną osią bezwładności (rys. 7.25). Ponadto przyjmujemy, że

początki nieruchomego układu współrzędnych x, y, z i ruchomego znajdują się w
nieruchomym punkcie O znajdującym się na osi obrotu l. Poza tym dla
uproszczenia wzorów założymy, że środek masy C bryły leży na osi

.

z z

= ′

x

Ponieważ dla takiego ruchu prędkość kątowa

ω i przyśpieszenie kątowe ε leżą

na osi obrotu, zatem

background image

0

ε

ε

i

0

y

x

y

x

=

=

=

ω

=

ω

, (b)

a wektory

ω i ε można zapisać wzorami:

,

dt

d

z

k

k

k

k

ω

z

ϕ

=

ω

=

ω

=

ω

=

.

dt

d

dt

d

ε

ε

ε

2

2

z

k

k

k

k

k

ε

z

ϕ

=

ω

=

=

=

=

Przyśpieszenie a

C

środka masy C obliczymy ze wzoru (5.37) podanego w p.

5.3.4 dotyczącym kinematyki ruchu obrotowego:

(

)

.

C

C

C

r

ω

ω

r

ε

a

×

×

+

×

=

Po podstawieniu do tego wzoru zależności

i

r

=

C

C

r

, wynikającej wprost z

rys. 7.25, otrzymamy:

(

)

i

j

i

k

k

i

k

a

ω

=

×

ω

×

ω

+

×

=

C

2

C

C

C

C

r

r

ε

r

r

ε

,

czyli współrzędne przyśpieszenia środka masy wynoszą:

0

a

,

r

ε

a

,

r

a

z

C

C

y

C

C

2

x

C

=

=

ω

=

. (c)

Wyprowadzone w poprzednim punkcie równania (7.93) po podstawieniu
zależności (b) oraz wzorów (c) redukują się do postaci (7.94):

.

M

ε

I

,

W

r

ε

m

,

W

r

m

z

O

z

z

y

C

x

C

2

=

=

=

ω

(7.94)

Stąd

M

M

W

Ox

Oy

z

=

=

=

0

0

,

oraz

0 .

(d)

Z

zależności (d) oraz z równań (7.94) wynika, że w przypadku obrotu bryły

wokół głównej osi bezwładności układ sił zewnętrznych redukuje się do momentu
głównego M

O

leżącego na osi obrotu l i wektora głównego W leżącego w

płaszczyźnie

′ ′

x y i prostopadłego do tej osi.

Trzecie

z

równań (7.94) jest dynamicznym równaniem ruchu obrotowego bryły

i przy znanych warunkach początkowych pozwala na wyznaczenie równania jej
ruchu

ϕ = ϕ(t). Z dwóch pierwszych równań możemy wyznaczyć siły wywołane

tym, że środek masy leży poza osią obrotu, czyli oś obrotu nie jest główną
centralną osią bezwładności, albo

− używając terminologii z dynamiki maszyn −

background image

bryła jest niewyważona statycznie. Równania te pozwalają na wyznaczenie reakcji
więzów (reakcji łożysk).
Jeżeli oś obrotu l będzie główną centralną osią bezwładności, czyli środek masy
C będzie leżał na osi obrotu (r

C

= 0), co będzie oznaczało idealne wyważenie

bryły, równania (7.94) redukują się do jednego równania:

z

O

z

z

M

ε

I

=

, (7.95)

a po uwzględnieniu (d) widzimy, że wszystkie współrzędne wektora głównego
oraz dwie współrzędne momentu głównego są równe zeru:

0

M

M

oraz

0

W

W

W

y

O

x

O

z

y

x

=

=

=

=

=

.

(e)

Z dynamicznego równania ruchu obrotowego bryły (7.95) wynika, że jeżeli
suma momentów wszystkich sił zewnętrznych (sił czynnych i reakcji łożysk osi
obrotu) względem osi obrotu będzie równa zeru, M

Oz

= 0 , to również

, zatem prędkość kątowa będzie stała,

ω = const, czyli bryła będzie

się poruszać ruchem jednostajnie obrotowym. Z takim przypadkiem będziemy
mieli do czynienia, gdy bryła będzie się obracać wokół pionowej osi obrotu
osadzonej w idealnie gładkich łożyskach. Siłami zewnętrznymi są wówczas siły
ciężkości i reakcje gładkich łożysk, których momenty względem osi obrotu są
równe zeru.

0

dt

d

ε

=

ω

=

/

Przykład 7.14. Jednorodna tarcza o masie m i
promieniu r obraca się wokół nieruchomej osi
przechodzącej przez środek O tej tarczy (rys.
7.26) pod wpływem przyłożonego
momentu o stałej wartości, M = const. Na tarczę
działa moment oporu

proporcjonalny do

prędkości kątowej

ω (

M

O

,

k

M

O

ω

=

gdzie k jest

znanym współczynnikiem). Wyznaczyć prędkość
kątową tarczy w funkcji czasu,

( )

t

ω

=

ω

, oraz jej

wartość maksymalną,

.

max

ω

=

ω

M

O

M

O

r

ω

Rys. 7.26. Wyznaczenie prędkości

kątowej tarczy


Rozwiązanie. Po podstawieniu do dynamicznego równania ruchu obrotowego
bryły (7.95), zgodnie z treścią zadania,

dt

d

ε

ε

,

I

I

z

O

z

ω

=

=

=

oraz

O

z

O

M

M

M

=

otrzymamy równanie ruchu obrotowego tarczy w postaci:

background image

ω

=

ω

=

ω

k

M

dt

d

I

lub

M

M

dt

d

I

O

O

O

.

Moment bezwładności tarczy względem osi obrotu

. Zatem

2

r

m

I

2

O

/

=

ω

=

ω

k

M

dt

d

r

m

2

1

2

.

Po rozdzieleniu zmiennych powyższe równanie różniczkowe możemy zapisać
w postaci:

dt

k

M

d

2

r

m

2

=

ω

ω

albo

dt

k

M

d

k

k

2

r

m

2

=

ω

ω

.

Scałkujemy to równanie w granicach od 0 do

ω oraz od 0 do t:

=

ω

ω

ω

t

0

0

2

dt

k

M

d

k

k

2

r

m

.

Po wykonaniu całkowania i zastąpieniu różnicy logarytmów logarytmem ilorazu
otrzymamy:

t

M

k

M

ln

k

2

r

m

2

=

ω

lub

2

r

m

t

k

2

M

k

M

ln

=

ω

.

Stąd prędkość kątowa



=

ω

2

mr

kt

2

e

1

k

M

.

Z otrzymanego wzoru widzimy, że z upływem czasu t do nieskończoności drugi
wyraz w nawiasie będzie dążył do zera, czyli prędkość kątowa

ω będzie dążyć do

wartości maksymalnej równej:

k

M

max

=

ω

.

background image

Obrót jednostajny bryły wokół osi dowolnej. Reakcje dynamiczne

Obecnie rozpatrzymy ruch bryły obracającej się ze stałą prędkością kątową

ω

wokół dowolnej osi podpartej w łożyskach, jak na rys. 7.27. Wskutek działania sił
czynnych na rozpatrywaną bryłę w łożyskach powstaną reakcje statyczne, które
można wyznaczyć z poznanych w statyce warunków równowagi. Zagadnienia tego
nie będziemy tutaj rozpatrywać, zajmiemy się natomiast siłami i momentami
wywołanymi przez zadany ruch. Innymi słowy, rozpatrzymy ruch bezwładny bryły
poruszającej się ze stałą prędkością kątową bez udziału sił zewnętrznych.
Na osi obrotu w punkcie O przyjmiemy początek nieruchomego układu
współrzędnych x, y, z oraz początek układu ruchomego

′ ′ ′

x , y , z sztywno

związanego z bryłą. Założymy przy tym, że osie układu ruchomego są głównymi
osiami bezwładności, a środek masy nie leży na osi obrotu, czyli bryła jest
niewyważona zarówno dynamicznie, jak i statycznie.
Ponieważ prędkość kątowa

ω jest stała i jej rzuty ω ω ω

x

,

,

y

z

na osie

ruchomego układu współrzędnych również są stałe, więc współrzędne
przyśpieszenia kątowego są równe zeru:

0

ε

ε

ε

z

y

x

=

=

=

.

(f)

x

z

x

y

z

y

O

C

ω

r

c

Rys. 7.27. Ruch obrotowy bryły sztywnej wokół osi dowolnej

Zatem przyśpieszenie a

C

środka masy bryły wyrazi wzór:

(

)

(

)

2

C

C

C

C

ω

r

r

ω

ω

r

ω

ω

a

=

×

×

=

. (g)

Jeżeli wektor wodzący r

C

środka masy zapiszemy za pomocą współrzędnych

w układzie ruchomym:

,

z

y

x

C

C

C

C

k

j

i

r

+

+

=

background image

to po zrzutowaniu wektora (g) na osie

′ ′ ′

x , y , z i odpowiednim pogrupowaniu

wyrazów otrzymamy wzory na współrzędne przyśpieszenia a

C

środka masy:

ω

ω

ω

ω

ω

ω

=

ω

ω

ω

ω

ω

ω

=

ω

ω

ω

ω

ω

ω

=

).

y

z

(

z

x

a

),

x

y

(

y

z

a

),

z

x

(

x

y

a

C

z

C

y

y

C

x

C

z

x

z

C

C

y

C

x

x

C

z

C

y

z

y

C

C

x

C

z

z

C

y

C

x

y

x

C

)

(

)

(

)

(

(h)

Po podstawieniu zależności (f) oraz wzorów (h) do równań (7.93) i zmianie
bieguna redukcji z C na O otrzymamy sześć równań opisujących omawiany ruch
bryły:

⎪⎪

=

ω

ω

=

ω

ω

=

ω

ω

=

ω

ω

ω

ω

ω

ω

=

ω

ω

ω

ω

ω

ω

=

ω

ω

ω

ω

ω

ω

.

M

I

I

,

M

I

I

,

M

I

I

,

W

y

z

z

x

m

,

W

x

y

y

z

m

,

W

z

x

x

y

m

z

O

y

x

x

y

y

O

z

x

z

x

x

O

z

y

y

z

z

C

z

C

y

y

C

x

C

z

x

y

C

y

C

x

x

C

z

C

y

z

x

C

x

C

z

z

C

y

C

x

y

)

(

)

(

)

(

)]

(

)

(

[

)]

(

)

(

[

)]

(

)

(

[

(7.96)

Po

uwzględnieniu we wzorze (7.91) zależności (f) oraz przyjęciu za biegun

redukcji zamiast punktu C nieruchomego punktu O pochodna krętu k

O

względem

czasu

O

O

dt

d

k

ω

k

×

=

, (i)

a po uwzględnieniu zasady krętu możemy napisać:

O

O

k

ω

M

×

=

. (j)

Po pomnożeniu skalarnie obu stron powyższego wzoru przez prędkość kątową

otrzymamy:

(

)

(

)

(

)

0

O

O

O

O

=

×

=

×

=

×

=

ω

ω

k

k

ω

ω

ω

k

ω

ω

M

. (k)

Warunek ten można przedstawić w postaci:

=

ω

M

O

0

M

M

M

z

z

O

y

y

O

x

x

O

=

ω

+

ω

+

ω

.

(l)

Z

powyższego równania wynika, że moment główny M

O

wywołany przez siły

bezwładności jest w czasie obrotu bryły zawsze prostopadły do prędkości kątowej
ω, czyli do osi obrotu. Gdy tak nie jest, obrót jednostajny bryły nie jest możliwy.

background image

Ponadto z warunku (l) wynika, iż tylko dwa z trzech ostatnich równań (7.96) są
niezależne, czyli z równań (7.96) możemy w układzie

′ ′ ′

x , y , z wyznaczyć pięć

składowych reakcji spowodowanych omawianym ruchem bryły. Ponieważ układ

wiruje razem z bryłą wokół osi obrotu z prędkością kątową

ω, z tą samą

prędkością wirują reakcje w łożyskach względem układu nieruchomego x, y, z.
Reakcje te nazywamy reakcjami dynamicznymi.

′ ′ ′

x , y , z

Gdy

środek masy bryły będzie się znajdował na osi obrotu, czyli bryła będzie

wyważona statycznie, wtedy

0

z

y

=

x

C

C

C

=

=

i lewe strony trzech pierwszych

równań (7.96) będą równe zeru, a tym samym znikną siły wywołane przez
niewyważenie statyczne

0

W

W

W

z

y

x

=

=

=

. W tym przypadku z trzech

ostatnich równań (7.96) wynika, że reakcje dynamiczne w łożyskach będą
spowodowane przez moment M

O

związany z działaniem sił bezwładności.

Ponieważ na podstawie warunku (l) moment ten jest prostopadły do osi obrotu,
zatem reakcje dynamiczne w łożyskach będą tworzyć parę sił wirującą z
prędkością równą prędkości kątowej

ω. Mówimy wtedy, że bryła jest niewyważona

dynamicznie

.

Jeżeli oś obrotu będzie główną centralną osią bezwładności, np. oś

pokryje

się z osią z, to pozostałe osie

z

x i y układu ruchomego będą do niej prostopadłe,

czyli

. Wynika z tego, że trzy pozostałe równania (7.96) znikają, a

tym samym znikają reakcje dynamiczne w łożyskach. Na podstawie powyższych
rozważań możemy sformułować następujący wniosek:

0

y

x

=

ω

=

ω

Jeżeli oś obrotu bryły jest główną centralną osią bezwładności, czyli bryła jest
wyważona statycznie i dynamicznie, to reakcje dynamiczne są równe zeru.

Z przeprowadzonych w tym punkcie rozważań wynika, że ruch wirującej bryły
wywołuje okresowo zmienne siły działające na łożyska, które przenosząc się na
korpus maszyny, a dalej na fundament wywołują drgania. Drgania te powodują
przyśpieszone zużycie elementów maszyny, a także niekorzystnie wpływają na
otoczenie. Aby temu zapobiec, wirujące części maszyn projektuje się tak, aby oś
obrotu była główną centralną osią bezwładności. Jednak np. ze względu na błędy
wykonawcze spełnienie tego warunku nie zawsze jest możliwe. Dlatego wirujące
części maszyn są sprawdzane po wykonaniu i ewentualnie wyważane przez
odpowiednią korektę masy.

Przykład 7.15. Cienka jednorodna płyta prostokątna o masie m i bokach h oraz
b obraca się wokół przekątnej ze stałą prędkością kątową

ω. Obliczyć reakcje

dynamiczne łożysk A i B, jeżeli odległość między nimi wynosi L (rys. 7.28).

background image

C

z

ω

R

B

R

B

x

L

h

b

z

x

α

A

B

Rys. 7.28. Wyznaczenie reakcji dynamicznych łożysk

Rozwiązanie. Ponieważ środek ciężkości C płyty leży na osi obrotu, która nie
jest główną centralną osią bezwładnóści, reakcje w łożyskach A i B będą
spowodowane niewyważeniem dynamicznym. W środku ciężkości przyjmiemy
ruchomy układ współrzędnych sztywno związany z płytą w ten sposób, że osie

są osiami symetrii płyty, a oś

x i z

y jest prostopadła do płaszczyzny rysunku.

W tym układzie współrzędnych prędkość kątowa

ω ma współrzędne:

α

ω

=

ω

=

ω

α

ω

=

ω

cos

,

0

,

sin

z

y

x

.

Po podstawieniu tych wzorów do trzech ostatnich równań (7.96) i po zastąpieniu
punktu O punktem C otrzymujemy:

0

M

,

0

M

z

C

x

C

=

=

oraz

(

)

(

)

α

α

ω

=

ω

ω

=

cos

sin

I

I

I

I

M

2

z

x

z

x

z

x

y

C

.

(a)

Momenty bezwładności prostokątnej płyty względem osi symetrii otrzymamy ze
wzorów (d) wyprowadzonych w przykładzie 6.3:

12

mb

I

,

12

mh

I

2

z

2

x

=

=

. (b)

Z rysunku wynika, że

2

2

2

2

b

h

h

=

cos

,

b

h

b

=

sin

+

α

+

α

. (c)

background image

Po podstawieniu oznaczeń (b) i (c) do wzoru (a) otrzymujemy:

(

)

2

2

2

2

2

C

y

C

b

h

bh

b

h

12

m

M

M

ω

+

=

=

. (d)

Z zależności (d) wynika, że wektor momentu

leży na osi

, czyli jest

prostopadły do płaszczyzny płyty i wiruje razem z nią. Moment ten jest wywołany
przez parę sił (reakcji) R

C

M

y

A

i R

B

prostopadłych do osi obrotu. Wartoci momentu i

reakcji są równe:

(

)

(

)

2

2

2

2

2

C

B

A

B

A

C

L

b

h

bh

b

h

12

m

L

M

R

R

L

R

L

R

M

ω

+

=

=

=

=

=

,

. (e)

W czasie obrotu reakcje R

A

i R

B

wirują razem z płytą. Ponadto są one

proporcjonalne do kwadratu prędkości kątowej i w przypadku zbyt szybko
obracającej się bryły mogą osiągać duże wartości.

background image

7.5.3. Ruch płaski bryły


W kinematyce ruchu bryły sztywnej ruchem płaskim nazwaliśmy ruch, w czasie
którego wszystkie punkty bryły zakreślają tory równoległe do pewnej płaszczyzny
nazywanej płaszczyzną ruchu lub płaszczyzną kierującą.

ϕ

y

z

x

z

y

x

r

C

O

M

C

W

ω

C

Rys. 7.29. Ruch płaski bryły sztywnej

Na rysunku 7.29 przedstawiono przekrój bryły płaszczyzną ruchu przechodzącą

przez środek masy C. W dowolnym punkcie O przyjęto nieruchomy układ
współrzędnych x, y, z tak, że osie x, y leżą w płaszczyźnie ruchu, a oś z jest do niej
prostopadła. Ruchomy układ współrzędnych

z

,

y

,

x

o początku w środku masy C

przyjęto w ten sam sposób, czyli osie

y

,

x ′

poruszają się w płaszczyźnie ruchu, a

jest do niej prostopadła. Wynika z tego, że osie

z

z i z

′ są do siebie

równoległe.
W dalszych rozważaniach dynamiki ruchu płaskiego bryły przyjmiemy
następujące założenia:

a) oś jest główną centralną osią bezwładności,

z

b) ruch bryły odbywa się pod wpływem sił działających w płaszczyźnie ruchu.

Bryła poruszająca się ruchem płaskim ma trzy stopnie swobody, a więc do jego
opisu wystarczy podać trzy równania ruchu

− dwóch współrzędnych środka masy

x

C

i y

C

oraz kąta obrotu

ϕ układu ruchomego względem nieruchomego.

Kinematyczne równania ruchu płaskiego (5.51) i (5.52) możemy zapisać w postaci:

( )

( )

( )

t

oraz

t

y

y

,

t

x

x

C

C

C

C

ϕ

=

ϕ

=

=

. (7.97)

background image

Zatem do opisu dynamiki ruchu płaskiego bryły niezbędne są trzy dynamiczne
równania ruchu. Do ich wyznaczenia wykorzystamy równania (7.93) opisujące
ruch bryły swobodnej.
Z

założenia b) na podstawie własności płaskiego układu sił (3.8) wynika, że

wektor główny W będzie leżał w płaszczyźnie sił, a moment główny M

C

będzie

prostopadły do tej płaszczyzny. Możemy w tej sytuacji zapisać:

0

M

oraz

0

x

C

=

=

=

y

C

z

M

W

. (m)

Ponadto w ruchu płaskim bryły (p. 5.3.8) prędkość kątowa

ω jest prostopadła do

płaszczyzny ruchu, czyli

0

y

x

=

ω

=

ω

. (n)

Po

uwzględnieniu zależności (m) i (n) równania (7.93) redukują się do trzech

dynamicznych równań ruchu płaskiego bryły.

z

C

z

z

y

y

C

x

x

C

M

ε

I

,

W

ma

,

W

ma

=

=

=

. (7.98)

Po

wyrażeniu przyśpieszenia a

C

środka masy oraz wektora głównego W

w nieruchomym układzie współrzędnych x, y oraz uwzględnieniu, że

(wzór 5.63), równania (7.98) można zapisać następująco:

ε

ε

=

z

z

C

z

y

Cy

x

Cx

M

ε

I

,

W

ma

,

W

ma

=

=

=

. (7.99)

Ponieważ współrzędne przyśpieszenia środka masy C w nieruchomym układzie
współrzędnych są równe drugim pochodnym względem czasu współrzędnych x

C

i

y

C

, powyższym równaniom można nadać postać

równań różniczkowych po uwzględnieniu drugiego wzoru (5.64):

z

C

2

2

z

y

C

2

x

2

C

2

M

t

d

d

I

,

W

t

d

y

d

m

,

W

t

d

x

d

m

=

ϕ

=

=

.

(7.100)


Przykład 7.16. Na poziomym szorstkim stole znajduje się szpula, której środek
masy C leży na osi symetrii obrotu. Szpula ma masę m oraz dwa promienie R i r.
Rysunek 7.30 przedstawia szpulę w rzucie na płaszczyznę prostopadłą do osi
symetrii. Moment bezwładności względem tej osi wynosi I

C

. Z obwodu

o promieniu r odwija się nić, do której końca przyłożono stałą siłę poziomą P.
Wyznaczyć maksymalną wartość siły P = P

max

, pod wpływem której szpula będzie

się toczyć bez poślizgu, jeżeli współczynnik tarcia statycznego między szpulą a
stołem jest równy

µ, a współczynnik tarcia tocznego f.

Dla tego przypadku wyznaczyć przyśpieszenie osi szpuli a

C

.

background image

Rozwiązanie. Na szpulę
działają dwie siły obciążające:
siła ciężaru szpuli G oraz siła P
powodująca ruch szpuli. Reakcję
stołu rozłożono na siłę tarcia T
skierowaną w kierunku
przeciwnym do kierunku ruchu
oraz reakcję normalną N
przesuniętą w kierunku toczenia
szpuli o wartość współczynnika
tarcia tocznego f (rys. 3.11b).
Rozważany ruch szpuli jest
ruchem płaskim, zatem na
podstawie wzoru (7.99)
dynamiczne równania ruchu
szpuli będą następujące :

y

x

O

f

P

G

N

T

C

R

r

a

C

ε

Rys. 7.30. Ruch szpuli z uwzględnieniem oporu

toczenia

=

=

=

.

f

N

r

P

R

T

ε

I

,

G

N

0

,

T

P

a

m

C

C

(a)

Jeżeli szpula toczy się bez poślizgu, to między przyśpieszeniem środka szpuli i
przyśpieszeniem kątowym musi być spełniona następująca C zależność
kinematyczna:

ε

R

a

C

=

. (b)

Z drugiego z równań (a) wynika, że reakcja normalna jest równa ciężarowi szpuli:

g

m

G

N

=

=

, (c)

gdzie g jest przyśpieszeniem ziemskim.
Maksymalną wartość siły P otrzymamy, założywszy, że siła tarcia T jest
graniczną siłą tarcia o wartości (wzór 3.5):

g

m

µ

N

µ

T

=

=

. (d)

Jeżeli do pierwszego i trzeciego równania (a) podstawimy wzory (c) i (d), a w
trzecim uwzględnimy zależność (b), otrzymamy dwa równania:

background image

⎪⎭

=

=

.

g

m

f

r

P

R

g

m

µ

R

a

I

,

g

m

µ

P

a

m

C

C

C

(e)

W równaniach tych mamy dwie niewiadome:

max

C

P

P

i

a

=

. W celu

wyeliminowania przyśpieszenia

podzielimy równania stronami i otrzymamy:

a

C

g

m

f

r

P

R

g

m

µ

g

m

µ

P

I

R

m

C

=

.

Stąd

(

)

r

R

m

I

g

m

I

µ

R

m

f

R

m

µ

P

P

C

C

2

max

+

+

=

=

.

(f)

Po podstawieniu tego wzoru do pierwszego równania (e) wyznaczymy
przyśpieszenie osi szpuli.

(

)

[

]

r

R

m

I

R

g

m

f

r

R

a

C

C

+

µ

=

. (g)

Z otrzymanego wzoru wynika, że oś szpuli porusza się ze stałym przyśpiesze-
niem, czyli ruchem jednostajnie przyśpieszonym. Czytelnikowi pozostawiamy
wyznaczenie równania ruchu

( )

?

t

x

x

C

C

=

=

dla warunków początkowych, np.

dla

0

v

i

0

x

,

0

t

C

C

=

=

=

.


Wyszukiwarka

Podobne podstrony:
dynamika 10 zadań z rozwiązaniami
Uklad Dynamiczny 10 Wikipedia
dynamika 10 zadań z rozwiązaniami
pamięci dynamiczne (10 str)
10 Dynamika cw3id 10826
Kopia (10) WYZNACZANIE MODUŁU SZTYWNOŚCI METODĄ DYNAMICZNĄ
Ćwiczenie nr 10 Bloki Dynamiczne
10 Dynamiczne przydzielanie pamieci
Mechanika - zestaw 10, Dynamika II
10 Dynamika cw3
Wyklad 10 Dynamika MS
Mechanika - Dynamika, cwiczeniadynamika10, Ćwiczenia 10
Zadania do ćwiczeń nr 10 – Dynamika punktu
Mechanika - Dynamika, dynamikawyklad10, Zasady ruchu dla punktu materialnego Wykład 10
10 Dynamiczna alokacja pamiecii Nieznany (2)
10 Dynamika konstrukcji
Dynamika Foliacji 10 Kulczycki p2

więcej podobnych podstron