Curs Mecanica
Curs Mecanica
Curs Mecanica
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
INTRODUCERE ÎN MECANICA CLASICĂ:<br />
TEORIE S¸I APLICAT¸ II<br />
Stan CHIRIT¸ Ă
Cuprins<br />
1 Cinematica 1<br />
1.1 Cinematica punctului material . . . . . . . . . . . . . . . . . . . 1<br />
1.1.1 Spat¸iu ¸si timp . . . . . . . . . . . . . . . . . . . . . . . . 1<br />
1.1.2 Mi¸scarea unui punct . . . . . . . . . . . . . . . . . . . . . 4<br />
1.1.3 Viteză ¸si accelerat¸ie . . . . . . . . . . . . . . . . . . . . . 6<br />
1.1.4 Mi¸scări plane . . . . . . . . . . . . . . . . . . . . . . . . . 12<br />
1.1.5 Viteza areolară . . . . . . . . . . . . . . . . . . . . . . . . 17<br />
1.1.6 Mi¸scări centrale . . . . . . . . . . . . . . . . . . . . . . . . 19<br />
1.1.7 Mi¸scări uniform variate ¸si periodice . . . . . . . . . . . . . 21<br />
1.1.8 Mi¸scări circulare ¸si uniforme . . . . . . . . . . . . . . . . 22<br />
1.1.9 Mi¸scări armonice . . . . . . . . . . . . . . . . . . . . . . . 25<br />
1.1.10 Mi¸scări elicoidale . . . . . . . . . . . . . . . . . . . . . . . 29<br />
1.2 Cinematica sistemelor materiale ¸si corpurilor rigide . . . . . . . . 30<br />
1.2.1 Legături ¸si sisteme olonome . . . . . . . . . . . . . . . . . 30<br />
1.2.2 Cinematica sistemelor rigide . . . . . . . . . . . . . . . . . 38<br />
1.2.3 Mi¸scări particulare ale rigidului . . . . . . . . . . . . . . . 41<br />
1.2.4 Unghiurile lui Euler . . . . . . . . . . . . . . . . . . . . . 44<br />
1.2.5 Starea de mi¸scare . . . . . . . . . . . . . . . . . . . . . . . 49<br />
1.2.6 Formula lui Poisson . . . . . . . . . . . . . . . . . . . . . 51<br />
1.2.7 Teorema lui Mozzi . . . . . . . . . . . . . . . . . . . . . . 53<br />
1.2.8 Aplicat¸ii . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58<br />
1.2.9 Cinematica mi¸scărilor relative . . . . . . . . . . . . . . . . 59<br />
1.2.10 Mi¸scări de transport speciale . . . . . . . . . . . . . . . . 64<br />
1.2.11 Mi¸scări relative pentru corpurile rigide . . . . . . . . . . . 65<br />
1.2.12 Aplicat¸ii ¸si exemple . . . . . . . . . . . . . . . . . . . . . 66<br />
1.2.13 Mi¸scări rigide plane . . . . . . . . . . . . . . . . . . . . . 69<br />
1.2.14 Traiectorii în coordonate polare . . . . . . . . . . . . . . . 74<br />
1.2.15 Accelerat¸ia unei mi¸scări rigide plane . . . . . . . . . . . . 79<br />
1.2.16 Mi¸scarea unui corp rigid cu un punct fix . . . . . . . . . . 80<br />
Bibliografie 82<br />
Index 86<br />
iii
iv CUPRINS
Capitolul 1<br />
Cinematica<br />
1.1 Cinematica punctului material<br />
1.1.1 Spat¸iu ¸si timp<br />
Cinematica studiază mi¸scarea corpurilor dintr-un punct de vedere pur descriptiv.<br />
Astfel, mi¸scarea este reprezentată ¸si studiată folosind mijloace matematice<br />
adecvate pornind de la legile fizice, care pun în legătură mi¸scarea cu cauzele<br />
(fort¸ele) care o determină.<br />
Fiecare fenomen de mi¸scare are loc într-un mediu spat¸io–temporal. Prin<br />
urmare, prima întrebare care urmează a fi discutată este descrierea conceptelor<br />
de spat¸iu ¸si timp. Este cunoscut faptul că acestea sunt not¸iuni primare, adică,<br />
ele nu sunt deduse din alte cantităt¸i, dar nu acesta este motivul pentru care nu<br />
este posibil să se obt¸ină o reprezentare matematică precisă a lor. Presupunem<br />
că spatiul ¸si timpul sunt continue, în sensul că este semnificativ de spus că un<br />
eveniment are loc într-un un anumit punct din spatiu ¸si la un anumit moment de<br />
timp ¸si că există standarde universale de lungime ¸si timp; cu alte cuvinte, observatori<br />
din locuri diferite la momente diferite de timp pot compara măsuratorile<br />
lor.<br />
Presupunem în continuare că există o scală universală pentru timp, ceea ce<br />
înseamnă că doi observatori care ¸si-au sincronizat ceasurile lor, vor fi întotdeauna<br />
de acord cu privire la timpul de producere a oricărui eveniment, în plus, noi<br />
presupunem că geometria spatiului este euclidiană ¸si faptul că, în principiu, nu<br />
există nicio limită a preciziei cu care putem masura pozit¸iile ¸si momentele.<br />
Astfel, în acest cadrul al mecanicii clasice, spat¸iu înconjurător este descris<br />
matematic ca un spat¸iu afin euclidian tridimensional ( 1 ). Aceasta înseamna un<br />
spat¸iu metric particular E, ale cărui elemente P, Q, . . . sunt numite puncte ¸si<br />
pentru care distant¸a are unele proprietăt¸i particulare ( 2 ).<br />
1 Această alegere, care, în contextul actual poate părea a fi evidentă, este de o mare<br />
important¸ă pentru dezvoltarea teoriei, a¸sa cum este strict legată de pricipiile mecanicii clasice.<br />
De fapt, diferite reprezentări ale conceptului de spat¸iu pot duce la descrieri diferite.<br />
2<br />
1
2 CAPITOLUL 1. CINEMATICA<br />
Asociem E cu spat¸iul vectorial tridimensional V , ale cărui elemente u, v, . . .<br />
sunt numite vectori. Fiecare vector u poate poate fi individualizat ca diferenta<br />
a două puncte ale spat¸iului E, adică<br />
u = P − Q.<br />
Pe V , considerăm not¸iunile obi¸snuite de produs scalar ¸si produs vectorial,<br />
care vor fi notate · ¸si respectiv × .<br />
În cadrul mecanicii clasice, not¸iunea de timp este definită ca un concept absolut,<br />
adică, derularea sa este independentă de obiectele ¸si entităt¸ile exterioare.<br />
Acest fapt ne permite să dăm o reprezentare relativ simplă a acestei not¸iuni.<br />
De fapt, folosind omogenitatea timpului (adică, faptul că momente privilegiate<br />
de timp nu există), este posibil să-l reprezintăm prin intermediul unui spat¸iu<br />
afin euclidian unudimensional R, ale cărui elemente sunt momente. Să notăm<br />
că această abordare a conceptului de timp nu este compatibilă cu principiile<br />
mecanicii relativiste, deoarece, în acest caz, durata unui fenomen depinde de<br />
cadrul de referint¸ă.<br />
În cele din urmă, să introducem o scală pentru a masura distant¸ele ¸si intervalele<br />
de timp folosind din nou proprietatea de omogenitate a spat¸iului ¸si timpului,<br />
sau, chiar mai bine spus, structura lor ca spat¸ii afine euclidiene. De fapt,<br />
este posibil, pe de o parte, să introducem scala pentru măsurarea distant¸elor<br />
prin intermediul unui e¸santion considerat a fi o unitate de lungime, ¸si, pe de<br />
alta parte, să folosim fenomene periodice pentru a reproduce unitatea de masurare<br />
a timpului ¸si, în consecint¸ă, pentru a defini ceasul. Comunitatea ¸stiintifica<br />
folose¸ste ca etalon de masură pentru lungime, care este, de lungimea un bar<br />
fabricat dintr-un aliaj de platina si iridiu păstrat la Bureau International des<br />
Poids et Mesures de Sèvres care ar trebui sa corespundă la 10 −7 din distanta de<br />
la Ecuator la Polul Nord masurată de-a lungul meridianul care trece prin Paris.<br />
Cel de-al doilea etalon a fost ales pentru o unitate de masură a timpului ¸si a<br />
fost init¸ial definit ca 24 −1 × 60 −2 dintr-o zi solare.<br />
Ar trebui să fie clar că reprezentarea conceptelor de spat¸iu ¸si timp prin<br />
intermediul unor modele matematice este o fază foarte delicată în construct¸ia<br />
Definit¸ie 1.1.1 Un spat¸iu metric E este numit spat¸iu euclidian afin tridimensional dacă<br />
metrica sa d : E × E → R + este astfel încât mult¸imea H a izometriilor, definită astfel<br />
H = {α : E → E , inversabilă; d(P, Q) = d[α(P ), α(Q)], pentru orice P, Q ∈ E} ,<br />
are următoarele proprietăt¸i:<br />
1. H este grup în raport cu legea de compunere;<br />
2. H cont¸ine un subgrup V , numit grupul translat¸iilor, care este abelian în raport cu legea<br />
de compunere;<br />
3. există operat¸ia de înmult¸ire cu scalari λ : R×V → V care face din V un spat¸iu vectorial<br />
tridimensional, cu operat¸ia de adunare (+) ca lege de compozit¸ie;<br />
4. există produsul scalar pe V, notat prin “punct” (·), astfel ca, pentru orice P, Q ∈ E ¸si<br />
pentru orice u ∈ V astfel ca u(P ) =Q, avem<br />
(d(P, Q)) 2 = u · u.
1.1. CINEMATICA PUNCTULUI MATERIAL 3<br />
principiilor mecanicii clasice. De fapt, diferite reprezentări ar putea conduce<br />
fie la complicat¸ii formale ale teoriei sau unele contradict¸ii logice care duc la<br />
dezvoltari complet diferite ale acesteia. De asemenea, ar trebui sa fie clar faptul<br />
că astfel de modele matematice sunt doar reprezentări ideale ale lumii fizice ¸si<br />
acestea ar trebui să fie considerate a fi în bună corespondent¸ă cu realitatea numai<br />
pentru studiul unor fenomene ¸si într-o aproximare adecvată.<br />
În particular, ele<br />
sunt adecvate pentru cazul în care vitezele sunt “mici”, fat¸ă de viteza luminii ¸si<br />
distant¸ele sunt “mari” cu privire la distante atomice.<br />
Pentru a specifica pozit¸iile ¸si momentele, fiecare observator poate alege o<br />
origine pe scala temporală, o origine în spat¸iu ¸si un set de trei axe de coordonate<br />
carteziene. Ne referim la toate acestea impreună spunând că s-a ales un cadru<br />
de referint¸ă.<br />
Pozit¸ia ¸si timpul fiecărui eveniment pot fi specificate fat¸ă de acest sistem<br />
cartezian de coordonate ¸si timp. Deoarece spat¸iul euclidian E este tridimensional,<br />
este posibil să fixăm un punct O ¸si să coniderăm trei direct¸ii mutual<br />
ortogonale x1, x2, x3 pornind din O. Asociem aceste direct¸ii cu trei vectori unitari<br />
i1, i2, i3 care formează un triplet drept, adică, direct¸iile lor coincid cu cele<br />
ale degetul mare, arătătorului ¸si degetului mijlociu de la mâna dreapta. Acest<br />
triplet centrat în O define¸ste un sistem de referint¸ă ( 3 )<br />
Prin urmare, este necesar să introducem pentru început conceptul de sistem<br />
material B. De fapt, acesta este definit ca o mult¸ime constituită dintr-un număr<br />
finit (sau infinit) de elementeX1, X2, X3, . . . , numite puncte materiale, înzestrat<br />
cu o familie P de aplicat¸ii injective ¸si netede ˜ P : B → E. O aplicat¸ie ˜ P este<br />
numită localizare a corpului B ¸si determină configurat¸ia specifică lui B în spat¸iul<br />
E.<br />
Definit¸ie 1.1.2 Un sistem material B este numit corp rigid dacă, pentru orice<br />
pereche de localizări ˜ P1, ˜ P2 ∈ P , avem<br />
<br />
d ˜P1 (X1) , ˜ <br />
P1 (X2) = d ˜P2 (X1) , ˜ <br />
P2 (X2)<br />
for all X1, X2 ∈ B. (1.1)<br />
Punctul ˜ P (X) poate fi acum identificat cu vectorul x = ( ˜ P (X) − O).<br />
Cu alte cuvinte, un sistem de material este un corp rigid, dacă toate posibilele<br />
configurat¸ii păstrează distant¸ele dintre punctele de material cu trecerea<br />
timpului.<br />
În mod natural, în cadrul mecanicii clasice, este presupus că astfel<br />
de sisteme rigide material există. Prin utilizarea acestei ipoteze fundamentale,<br />
este posibil să considerăm un sistem referint¸ă cartezian invariant cu trecerea<br />
timpului, ca ¸si cum ar fi fixat într-un corp rigid.) în spat¸iu (Figure 1.1). Astfel,<br />
fiecare vector x = P − O poate fi reprezentat în următoarea formă:<br />
x = x1i1 + x2i2 + x3i3,<br />
3 Cu scopul de a da definit¸ia corectă a “sistemului de referinta”, să ne reamintim faptul că<br />
notiunea de mi¸scare este un concept relativ care implică prezent¸a altor obiecte sau corpuri<br />
capabile a fi observate, astfel încât este posibil să ne referim la ea ca miscare fat¸ă de aceste<br />
obiecte.
4 CAPITOLUL 1. CINEMATICA<br />
x 1<br />
i 3<br />
O<br />
i 1<br />
x 2<br />
x 3<br />
i 2<br />
P<br />
x 3<br />
Figura 1.1:<br />
x 1<br />
unde x1, x2, x3 sunt numite coordonate ale lui x în raport cu cele trei axe.<br />
Ele sunt obt¸inute ca proiect¸ii ale lui x pe axele i1, i2, i3, adică, x1 = x · i1,<br />
x2 = x · i2, x3 = x · i3. Un sistem de referint¸ă există în afara not¸iunii de timp<br />
¸si, în acest moment, are un sens pur matematic. Cu alte cuvinte, un reper<br />
cartezian ortogonal nu poate păstra propriile caracteristici cu trecerea timpului.<br />
Definit¸ie 1.1.3 Numim cadru de referintă, un set de trei axe de coordonate<br />
(sistemul de referinta), fixat într-un corp rigid împreună cu un sistem de măsurare<br />
a timpului (ceasul).<br />
Este clar că, în scopul de a introduce not¸iunea de cadru de referint¸ă, este<br />
necesar să presupunem existent¸a unui corp rigid. Mai mult, sistemul de referint¸ă<br />
fix într-un corp rigid va fi uneori confundat cu tripletul ortogonal (O, x1, x2, x3).<br />
1.1.2 Mi¸scarea unui punct<br />
În prima parte a cinematicii, studiem mi¸scarea unui punct material, aceasta<br />
înseamnă, mi¸scarea unui sistem de material format dintr-un unic punct P . Un<br />
astfel de sistem de material reprezintă foarte frecvent un bun model pentru<br />
studiul corpurilor ale căror dimensiuni sunt suficient de mici pentru a permite<br />
aceasta reprezentare; de asemenea, poate reprezenta, un anumit punct al sistemului<br />
material.<br />
Mi¸scarea unui punct P este definită complet de aplicat¸ia<br />
x 2<br />
ˆP : I → E, (1.2)
1.1. CINEMATICA PUNCTULUI MATERIAL 5<br />
unde I ⊂ R este intervalul de timp în care este definită mi¸scarea, ¸si E este<br />
spat¸iul euclidian care poate fi asociat cu un cadru de referint¸ă sau cu un triplet<br />
rigid. Mi¸scarea va fi notată cu simbolul<br />
P = ˆ P (t),<br />
sau prin funct¸ia vectorială x(t) definită de<br />
x (t) = ˆ P (t) − O, (1.3)<br />
sau prin intermediul funct¸iilor componente ale ecuat¸iilor vectoriale (1.3)<br />
x1 = ˆx1 (t) , x2 = ˆx2 (t) , x3 = ˆx3 (t) . (1.4)<br />
Ulterior, funct¸iile care definesc miscarea sunt presupuse a fi cel put¸in de<br />
clasă C 2 . Imaginea intervalului I în E define¸ste traiectoria punctului P relativ<br />
la mi¸scarea P (t). Putem descrie această traiectorie intrinsec, independent de<br />
variabila temporală. Fie un punct fix O1 ¸si o direct¸ie pozitivă pe traiectorie ¸si<br />
să indicăm prin s abscisa curbilinie a lui P , care reprezintă distant¸a cu semn<br />
de la P la O1 măsurată de-a lungul traiectoriei (a se vedea Figura 1.2), Atunci,<br />
traiectoria este dată de funct¸ia<br />
P = ˆ P (s). (1.5)<br />
Mai mult, în mi¸scare, abscisa curbilinie s este o funct¸ie de timp t, exprimată<br />
printr-o funct¸ie s = ˆs(t), pe care o numim ecuat¸ie orară a mi¸scării lui P .<br />
Prin urmare, mi¸scarea punctului P poate fi descrisă de către sistemul<br />
P = ˆ P (s), s = ˆs(t), (1.6)<br />
unde prima funct¸ie define¸ste traiectoria, în timp ce legea temporală asociată<br />
punctului P oferă pozit¸ia instantanee a punctului de-a lungul traiectoriei.<br />
Traiectoria punctului P cu privire la cadrul de referint¸ă ales poate fi descris<br />
prin funct¸iile de ˆx1(s), ˆx2(s), ˆx3(s) definite ca proiect¸ii ale relat¸iei vectoriale<br />
(1.5), adică<br />
x1 = ˆx1 (s) , x2 = ˆx2 (s) , x3 = ˆx3 (s) . (1.7)<br />
Exercit¸iu 1.1.1 Mi¸scarea unui punct este dată de x1 = R cos √ t, x2 = R sin √ t,<br />
x3 = R √ 3t, t ∈ [0, π 2 ], R > 0. Să se determine traiectoria ¸si ecuat¸ia orară a<br />
acestei mi¸scări.<br />
<br />
Solut¸ie. Deoarece ds = |dx| = (dx1) 2 + (dx2) 2 + (dx3) 2 , deducem că,<br />
pentru mi¸scarea noastră,<br />
t t<br />
R<br />
s = ds(z) = √ dz = 2R<br />
z √ t,<br />
0<br />
¸si deci ecuat¸ia orară este s = 2R √ t. Substituind √ t = s<br />
0<br />
2R<br />
în ecuat¸iile de mi¸scare,<br />
obt¸inem următoarea expresie a traiectoriei: x1 = R cos s<br />
2R , x2 = R sin s<br />
2R ,<br />
x3 = s√3 2 , s ∈ [0, 2πR].
6 CAPITOLUL 1. CINEMATICA<br />
x 1<br />
O<br />
x 3<br />
O 1<br />
s<br />
Figura 1.2:<br />
P<br />
Exercit¸iu 1.1.2 Presupunem că traiectoria unui punct P este cercul de intersect¸ie<br />
dintre sfera x 2 1 + x 2 2 + x 2 3 = R 2 , R > 0, cu planul x3 = R cos θ0, θ0 fixat<br />
în (0, π). Să se determine mi¸scarea punctului P având ecuat¸ia orară s = t,<br />
t ∈ [0, √ 2πR sin θ0].<br />
Solut¸ie. Traiectoria este descrisă de<br />
x1 = R sin θ0 cos ϕ, x2 = R sin θ0 sin ϕ, x3 = R cos θ0, ϕ ∈ [0, 2π].<br />
Deoarece s = (R sin θ0) ϕ ¸si ecuat¸ia orară este s = t, rezultă că ϕ =<br />
t ∈ [0, √ 2πR sin θ0]. Dacă alegem a = R sin θ0, atunci mi¸scarea este dată de<br />
x1 = a cos t<br />
a , x2 = a sin t<br />
a , x3 = R cos θ0, t ∈ [0, √ 2πa].<br />
1.1.3 Viteză ¸si accelerat¸ie<br />
x 2<br />
t<br />
R sin θ0 ,<br />
Considerăm un punct material a cărui mi¸scare este descrisă de sistemul (1.6).<br />
Presupunând traiectoria P = ˆ P (s) fixată, mi¸scarea este definită simplu de<br />
ecuat¸ia orară s = ˆs(t).<br />
Definit¸ie 1.1.4 Numim viteză a punctului P de-a lungul traiectoriei ˆ P (s)<br />
derivata lui ˆs în raport cu timpul, ¸si o notăm cu ( 4 ) prin<br />
v(t) def<br />
= dˆs(t)<br />
. (1.8)<br />
dt<br />
4 Pentru a evita confuziile, vom nota întotdeauna derivata în raport cu timpul printr-un<br />
punct plasat deasupra funct¸iei considerate, adică, dˆs/dt = ˙s, ¸si nu vom face distinct¸ie dintre<br />
punctul P ¸si funct¸ia corespunzătoare ˆ P .
1.1. CINEMATICA PUNCTULUI MATERIAL 7<br />
Dacă ˙s > 0, atunci mi¸scarea este numită directă, în timp ce pentru ˙s < 0,<br />
mi¸scarea este numită retrogradă; momentele la care ˙s = 0 sunt numite momente<br />
de stat¸ionare. În final, dacă ˆs este o funct¸ie liniară în timp, adică,<br />
atunci mi¸scarea este numită uniformă.<br />
ˆs(t) = vt + s0,<br />
(1.9)<br />
În general, atunci când vorbim despre viteza unui punct, întotdeauna înt¸elegem<br />
o cantitate vectorială. Într-adevar, presupunând că am fixat cadrul de referint¸ă<br />
¸si că mi¸scarea este descrisă în raport cu acest cadru de ecuat¸ia de mi¸scare<br />
P = ˆ P (t), definim vectorul viteză astfel:<br />
Definit¸ie 1.1.5 Vectorul<br />
v(t) def<br />
= d ˆ P (t) d<br />
<br />
= ˆP (t) − O<br />
dt dt<br />
(1.10)<br />
este numit viteza punctului P fat¸ă de cadrul de referint¸ă considerat, reprezentat<br />
de sistemul de referint¸ă (O, x1, x2, x3).<br />
Mai mult, dacă i1, i2, i3 sunt vectorii unitari ortogonali ai sistemului (O, x1, x2, x3),<br />
atunci, folosind relat¸ia ˆ P − O = ˆx1i1 + ˆx2i2 + ˆx3i3 = x, obt¸inem<br />
v(t) = ˙x1(t)i1 + ˙x2(t)i2 + ˙x3(t)i3 = d<br />
x(t), (1.11)<br />
dt<br />
unde ˙x1(t), ˙x2(t), ˙x3(t) sunt derivatele funt¸iilor ˆx1, ˆx2, ˆx3 în raport cu timpul ¸si<br />
reprezintă componentele vectorului v de-a lungul axelor x1, x2, x3.<br />
Din definit¸ia de mai sus, folosind expresiile (1.6) pentru mi¸scarea punctului<br />
P , adică<br />
P = ˆ P (ˆs(t)), (1.12)<br />
rezultă că<br />
v(t) = d ˆ P<br />
ds (s)dˆs . (1.13)<br />
dt<br />
Acum, considerăm următoarea cantitate:<br />
Rata de cre¸stere ˆ P (s+h)− ˆ P (s)<br />
h<br />
d ˆ P<br />
(s) = lim<br />
ds h→0<br />
ˆP (s + h) − ˆ P (s)<br />
. (1.14)<br />
h<br />
define¸ste un vector a cărui direct¸ie este de-a lungul<br />
coardei care une¸ste ˆ P (s) cu ˆ P (s + h) ¸si direct¸ia coincide cu cea de cre¸stere a<br />
arcelor. Când h se apropie de zero, această rat¸ie tinde spre un vector a cărui<br />
direct¸ie este este paralelă tangenta la curbă în punctul P (s) (Figura 1.3).<br />
Considerăm acum mărimea dată de expresia (1.14)<br />
<br />
<br />
d<br />
<br />
<br />
ˆ P<br />
ds (s)<br />
<br />
<br />
<br />
<br />
=<br />
<br />
<br />
<br />
<br />
lim<br />
ˆP (s + h) −<br />
h→0<br />
ˆ <br />
P (s)<br />
<br />
<br />
<br />
h .
8 CAPITOLUL 1. CINEMATICA<br />
x 1<br />
O<br />
x 3<br />
O 1<br />
Rat¸ia | ˆ P (s+h)− ˆ P (s)|<br />
|h|<br />
s<br />
Figura 1.3:<br />
t<br />
P(s)<br />
h<br />
P(s + h)<br />
x 2<br />
nu este nimic altceva decât raportul dintre lungimea coardei<br />
ce une¸ste punctele ˆ P (s) ¸si ˆ P (s+h), ¸si arcul corespunzător. Este cunoscut faptul<br />
că acest raport tinde spre unitate atunci când lungimea arcului se apropie de<br />
zero.<br />
Astfel, putem concluziona că<br />
d ˆ P<br />
(s) = t(s), (1.15)<br />
ds<br />
unde t este versorul tangentei la traiectoria punctului ˆ P (s) ¸si a cărui direct¸ie<br />
coincide cu direct¸ia de cre¸stere a arcului. Prin urmare, din (1.13), obt¸inem<br />
v(t) = ˙st. (1.16)<br />
Rezultă din argumentele de mai sus că viteza este întotdeauna îndreptată dea<br />
lungul tangentei la traiectoria punctului considerat. În particular, observăm<br />
că mărimea vitezei v, pe care o notăm cu v sau cu |v|, este definită de<br />
<br />
v = | ˙s| = ˙x 2 1 + ˙x2 2 + ˙x2 3 .<br />
Dacă viteza punctului P este constantă pentru un interval de timp, mi¸scarea<br />
este numită rectilinie ¸si uniformă în acest interval. După cum se ¸stie, expresia<br />
(1.10) poate fi scrisă în forma echivalentă<br />
ˆP (t + ∆t) − ˆ P (t)<br />
= v(t) + ε(∆t), (1.17)<br />
∆t<br />
unde ε este un vector, astfel încât lim∆t→0 ε(∆t) = 0. În consecint¸ă, din (1.17),<br />
obt¸inem<br />
∆P def<br />
= ˆ P (t + ∆t) − ˆ P (t) = v(t)∆t + ε(∆t)∆t.
1.1. CINEMATICA PUNCTULUI MATERIAL 9<br />
Vectorul ∆P reprezintă deplasarea punctului P în intervalul de timp ∆t. Dacă<br />
∆t este “suficient de mic”, v∆t dă o bună aproximare a deplasării, adică<br />
∆P v∆t.<br />
Luând în considerare această observat¸ie, introducem vectorul dP def<br />
= vdt pe<br />
care îl numim deplasare elementară. În general, acesta nu corespunde cu deplasarea<br />
reală, dar el reprezintă o bună aproximare pentru ea, sub presupunerea<br />
că valorile lui dt sunt “mici”.<br />
Definit¸ie 1.1.6 Numim vectorul<br />
a(t) def<br />
= dv<br />
(t) (1.18)<br />
dt<br />
accelerat¸ie a punctului P fat¸ă de cadru de referint¸ă ales.<br />
sau<br />
Folosind formula vitezei, este posibil să arătăm că<br />
a(t) = d2Pˆ (t),<br />
dt2 a(t) = ¨x1(t)i1 + ¨x2(t)i2 + ¨x3(t)i3 = d2x (t),<br />
dt2 unde ¨x1 = d2 ˆx1<br />
dt2 , ¨x2 = d2 ˆx2<br />
dt2 , ¨x3 = d2 ˆx3<br />
dt2 pot fi considerate ca fiind accelerat¸iile<br />
proiect¸iilor punctului P de-a lungul axelor x1, x2, x3.<br />
Folosind formula (1.16) a vitezei, obt¸inem următoarea expresie importantă<br />
pentru accelerat¸ia unui punct:<br />
a = d<br />
( ˙st) = ¨st + ˙sdt(s(t)) = ¨st + ˙s<br />
dt dt<br />
Este necesar acum să considerăm limita<br />
dt<br />
= lim<br />
ds h→0<br />
2 dt<br />
. (1.19)<br />
ds<br />
t(s + h) − t(s)<br />
. (1.20)<br />
h<br />
Considerăm planul π definit de triunghiul (P (s), A, B) (a se vedea Figura<br />
1.4). Dacă curba este într-un plan, atunci planul π va coincide cu planul care<br />
cont¸ine curba. Dacă curba nu este într-un plan, atunci, când h se apropie de<br />
zero, π tinde spre un plan care trece prin P (s) ¸si cont¸ine vectorul t(s). Numim<br />
acest plan plan osculator. Deoarece rat¸ia t(s+h)−t(s)<br />
h are aceea¸si direct¸ie cu<br />
B − A, limita acestei rat¸ii, adică dt/ds, trebuie să apart¸ină planului osculator.<br />
Mai mult, deoarece mărimea vectorului t(s) este constantă, concluzionăm că<br />
dt/ds trebuie să fie ortogonal pe t, ¸si, în consecint¸ă, ortogonal curbei ( 5 ); în<br />
final, el trebuie să fie orientat spre centrul curbei. Mai mult,<br />
<br />
<br />
<br />
dt <br />
<br />
|t(s + h) − t(s)| |t(s + h) − t(s)| ∆α<br />
ds<br />
= lim<br />
= lim<br />
, (1.21)<br />
h→0 |h|<br />
h→0 ∆α |h|<br />
5 dt<br />
dt<br />
Reamintim faptul că, dacă t· t = 1, atunci 2 · t = 0, deci este ortogonal pe t.<br />
ds ds
10 CAPITOLUL 1. CINEMATICA<br />
P(s)<br />
t(s)<br />
A<br />
∆α<br />
B<br />
t(s + h)<br />
∆α<br />
Figura 1.4:<br />
P(s + h)<br />
t(s + h)<br />
unde ∆α este unghiul dintre vectorii t(s) ¸si t(s+h), măsurat în radiani. Alegem<br />
unde 1<br />
ρ<br />
∆α 1<br />
lim =<br />
h→0 |h| ρ ,<br />
este numită curbura, iar ρ este numită raza de curbură. În plus,<br />
|t(s + h) − t(s)|<br />
lim<br />
= 1,<br />
h→0 ∆α<br />
deoarece considerăm limita dintre lungimea arcului ¸si coarda de sprijin a acestuia.<br />
Cercul de rază ρ, situat în planul osculator, tangent la traiectoria lui P (s)<br />
¸si ales din două cercuri posibile tangente la traiectoria lui P (s) ca cel situat pe<br />
partea concavă a traiectoriei, este numit de obicei cerc osculator. Mai mult,<br />
prin n(s) notăm vectorul unitar normal principal al traiectoriei lui P (s), adică,<br />
vectorul unitar care este ortogonal la curbă, se afla în planul osculator ¸si este<br />
îndreptat spre centru curbei. Atunci, din relat¸iile (1.20) ¸si (1.21), obt¸inem ( 6 )<br />
dt<br />
ds<br />
1<br />
= n. (1.22)<br />
ρ<br />
Acum, putem oferi o reprezentare importantă pentru vectorul accelerat¸ie. Astfel,<br />
din relat¸iile (1.19) ¸si (1.22), obt¸inem<br />
a = ¨st + ˙s2<br />
n. (1.23)<br />
ρ<br />
Este convenabil să considerăm de asemenea – împreuna cu cei doi vectori unitari<br />
t ¸si n – vectorul unitar b, numit binormală, care este ortogonal pe t ¸si n astfel<br />
6 Aceasta ecuat¸ie este numită prima formulă a lui Frènet. Pentru detalii, a se vedea<br />
sect¸iunea A.8 din Appendix A.
1.1. CINEMATICA PUNCTULUI MATERIAL 11<br />
încât t, n, ¸si b formează un reper drept. Acest sistem de vectori este numit<br />
triplet intrisec de vectori pentru traiectoria punctului P.<br />
Folosind expresia (1.23), se poate observa că, spre deosebire de viteza, în<br />
afara de componenta at = ¨st orientată de-a lungul tangentei ¸si numită accelerat¸ie<br />
tangent¸ială, accelerat¸ia are ¸si altă componentă an = ˙s2<br />
ρ n orientată de-a lungul<br />
normalei, numită accelerat¸ie normală sau centripetă. Dacă accelerat¸ia tangent¸ială<br />
se anulează, atunci este necesar ca<br />
¨s = 0,<br />
deci ˙s = constant, ¸si mi¸scarea este uniformă. Dacă accelerat¸ia centripetă se<br />
anulează în intervaul de timp, adică ˙s2<br />
1<br />
ρ n = 0, ¸si ˙s(t) = 0, atunci curbura ρ = 0,<br />
¸si deci mi¸scarea este rectilinie. În sfâr¸sit, dacă a = at + an = 0 pentru un<br />
interval de timp, atunci mi¸scarea este rectilinie ¸si uniformă.<br />
Definit¸ie 1.1.7 Mi¸scarea unui punct este numită accelerată la un moment dat<br />
dacă mărimea vitezei la acel moment este o funct¸ie de t crescătoare; mi¸scarea<br />
este numită încetinită dacă mărimea vitezei la acel moment este o funct¸ie de t<br />
descrescătoare;.<br />
Deoarece<br />
d<br />
dt ˙s2 = 2¨s ˙s, (1.24)<br />
rezultă că mi¸scarea poate fi accelerată sau încetinită în funct¸ie de cum ˙s ¸si ¨s,<br />
ambele diferite de zero, au sau nu au acela¸si semn.<br />
d<br />
În primul caz, avem dt<br />
˙s 2 > 0, ¸si deci ˙s 2 este o funct¸ie crescătoare în timp, în timp ce, în cel de-al<br />
doilea caz, d<br />
dt ˙s2 < 0 ¸si atunci ˙s 2 este o funct¸ie de t descrescătoare.<br />
Exercit¸iu 1.1.3 Punctul P se mi¸scă pe curba x1 = 2e 2t , x2 = 3 sin 2t, x3 =<br />
2 cos 2t, t ∈ R. Să se determine vectorul viteză ¸si vectorul accelerat¸ie la momentul<br />
t. Calculat¸i mărimile vitezei ¸si accelerat¸iei la momentul t = 0.<br />
Solut¸ie. Vectorul deplasare este x = 2e 2t i1 + 3 sin 2ti2 + 2 cos 2ti3 ¸si deci<br />
deducem că v(t) = ˙x(t) = 4e 2t i1 + 6 cos 2ti2 − 4 sin 2ti3 ¸si a(t) = ¨x(t) = 8e 2t i1 −<br />
12 sin 2ti2 − 8 cos 2ti3. Pentru t = 0, avem v(0) = 4i1 + 6i2 ¸si a(0) = 8i1 − 8i3<br />
¸si prin urmare v = √ 16 + 36 = 2 √ 13 ¸si a = √ 64 + 64 = 8 √ 2.<br />
Exercit¸iu 1.1.4 Un punct P porne¸ste din pozit¸ia P0(−3, 2, 1) la timpul t = 0<br />
cu viteza init¸ialăv0 = −i1 + 2i2 + 3i3 ¸si se deplasează cu accelerat¸ia a = e −t i1 +<br />
4 cos 2ti2 + 8 sin 2ti3. Să se găsească vectorul viteză a punctului ¸si ecuat¸iile de<br />
mi¸scare.<br />
Solut¸ie. Din relat¸ia ˙v(t) = a(t), prin integrare în raport cu timpul t,<br />
deducem că v(t) = −e −t i1 + 2 sin 2ti2 − 4 cos 2ti3 + c1, unde c1 este un vector<br />
constant arbitrar. Deoarece v(0) = v0, obt¸inem −i1−4i3+c1 = −i1+2i2+3i3 ¸si<br />
deci avem c1=2i2+7i3 ¸si viteza este v(t) = −e −t i1+2(sin 2t+1)i2+(−4 cos 2t+<br />
7)i3.
12 CAPITOLUL 1. CINEMATICA<br />
Din relat¸ia ˙x(t) = v(t), obt¸inem x(t) = e −t i1 +(− cos 2t+2t)i2 +(−2 sin 2t+<br />
7t)i3 + c, unde c este o constantă oarecare. Deoarece x(0) = −3i1 + 2i2 + i3,<br />
rezultă că i1 − i2 + c = −3i1 + 2i2 + i3 ¸si deci c = −4i1 + 3i2 + i3. Astfel,<br />
mi¸scarea punctului este descrisă de x(t) = (e −t − 4)i1 + (− cos 2t + 2t + 3)i2 +<br />
(−2 sin 2t + 7t + 1)i3.<br />
Exercit¸iu 1.1.5 Mi¸scarea unui punct P este x(t) = ti1 + 1<br />
2t2i2 + 1<br />
6t3i3. Să<br />
se determine accelerat¸ia tangent¸ială ¸si accelerat¸ia normală a punctului la un<br />
moment t.<br />
Solut¸ie. Deoarece ˙x = i1 + ti2 + 1<br />
2 t2 i3 ¸si ¨x = i2 + ti3, rezultă că ds = | ˙x| dt<br />
¸si deci ˙s = 1<br />
2 (t2 + 2), ¨s = t. Mai mult, versorul tangentei la curbă este t =<br />
2i1 + 2ti2 + t2 <br />
i3 ¸si<br />
1<br />
t 2 +2<br />
¸si deci<br />
dt<br />
ds<br />
dt dt<br />
=<br />
dt ds =<br />
=<br />
2<br />
(t 2 + 2) 2 [−2ti1 − (t 2 − 2)i2 + 2ti3] dt<br />
ds =<br />
4<br />
(t2 + 2) 3 [−2ti1 − (t 2 − 2)i2 + 2ti3] = 1<br />
ρ n,<br />
n = 1<br />
t2 + 2 [−2ti1 − (t 2 − 2)i2 + 2ti3], 1<br />
ρ =<br />
4<br />
(t2 .<br />
+ 2) 2<br />
Deci, putem concluziona că accelerat¸ia tangent¸ială este at = tt ¸si accelerat¸ia<br />
centripetă este an = n.<br />
1.1.4 Mi¸scări plane<br />
Considerăm punctul P care se mi¸scă într-un plan. Este posibil să descriem<br />
mi¸scarea lui P prin intermediul unui sistem de ecuat¸ii de următoarea formă:<br />
x1 = ˆx1(t),<br />
x2 = ˆx2(t),<br />
unde x1, x2 sunt coordonatele carteziene ale lui P relativ la un sistem de<br />
referint¸ă din planul de mi¸scare. Atunci, viteza ¸si accelerat¸ia pot fi determinate<br />
aplicând formulele din sect¸iunea de mai sus.<br />
Un interesant capitol particular în studiul mi¸scării plane este cel al sistemului<br />
de coordonate polare (ρ, θ), unde ρ = |P − O| ¸si θ = ∠ (OP, Ox1) sunt numite<br />
distant¸ă polară ¸si respectiv unghi polar. Mi¸scarea punctului P va fi descrisă în<br />
coordinate polare de sitemul (Figura 1.5)<br />
ρ = ˆρ(t), θ = ˆ θ(t).<br />
Eliminând variabila t din acest ultim sistem, obt¸inem ecuat¸ia polară ρ = ˆρ(θ)<br />
a traiectoriei punctului P . Deci, dacă O este originea sistemului de referint¸ă ¸si<br />
r = , atunci, prin derivarea identităt¸ii<br />
P −O<br />
|P −O|<br />
(P − O) = ρr, (1.25)
1.1. CINEMATICA PUNCTULUI MATERIAL 13<br />
h<br />
O<br />
x 2<br />
r<br />
θ<br />
Figura 1.5:<br />
P<br />
obt¸inem<br />
d(P − O)<br />
v = = ˙ρr+ρ<br />
dt<br />
dr<br />
. (1.26)<br />
dt<br />
Vectorul unitar r depinde de timp prin variabila θ, adică,<br />
¸si prin urmare, din (1.26), obt¸inem<br />
x 1<br />
r(t) = r(θ(t)), (1.27)<br />
v = ˙ρr+ρ ˙ θ dr<br />
. (1.28)<br />
dθ<br />
Dacă considerăm reperul cartezian (O, x1, x2) cu originea în O, ¸si axa x1 coincizând<br />
cu axa polară, avem<br />
Dacă h = − sin θi1 + cos θi2, rezultă din (1.29) că<br />
r = cos θi1 + sin θi2. (1.29)<br />
h = dr<br />
, (1.30)<br />
dθ<br />
¸si, în consecint¸ă, |h| = 1 ¸si h · r = 0. Prin urmare, h este un vector unitar<br />
ortogonal pe r inclus în planul (x1, x2). Pe de altă parte, este u¸sor de observat<br />
că<br />
dh<br />
= −r.<br />
dθ<br />
Întorcându-ne la (1.28), obt¸inem<br />
(1.31)<br />
v = ˙ρr+ρ ˙ θh. (1.32)<br />
Observat¸ie 1.1.1 Viteza punctului P , exprimată în coordonate polare, poate<br />
fi reprezentată ca suma a doi termeni: primul termen, vρ = ˙ρr, este numit<br />
viteză radială, iar cel de-al doilea termen, vθ = ρ ˙ θh, este numit vectorul viteză
14 CAPITOLUL 1. CINEMATICA<br />
unghiulară (transversală). Deoarece vρ ¸si vθ sunt ortogonali, mărimea vitezei<br />
este dată de formula<br />
<br />
v = ˙ρ 2 + ρ2 ˙ θ2 .<br />
Prin derivare directă a relat¸iei (1.32), obt¸inem următoarea expresie a accelerat¸iei<br />
în coordonate polare:<br />
a = dv<br />
dt = ¨ρr+ ˙ρ ˙ θ dr<br />
dθ + ˙ρ ˙ θh + ρ ¨ θh + ρ ˙ θ<br />
Folosind relat¸iile (1.30) ¸si (1.31) în (1.33), deducem că<br />
2 dh<br />
. (1.33)<br />
dθ<br />
a = (¨ρ − ρ ˙ θ 2 )r + (ρ ¨ θ + 2 ˙ρ ˙ θ)h. (1.34)<br />
Observat¸ie 1.1.2 Accelerat¸ia punctului P , exprimată în coordonate polare,<br />
poate fi reprezentată ca suma a doi termeni: primul termen, aρ = (¨ρ−ρ ˙ θ2 )r, este<br />
numit accelerat¸ie radială, ¸si cel de-al doilea termen, aθ = (ρ¨ θ+2 ˙ρ ˙ θ)h, este numit<br />
accelerat¸ie unghiulară (sau transversală). Deoarece ρ¨ θ + 2 ˙ρ ˙ θ = 1<br />
ρ (ρ2 ¨ θ + 2ρ ˙ρ ˙ θ),<br />
putem de asemenea exprima aθ ca aθ = 1 d<br />
ρ dt (ρ2 ˙ θ)h.<br />
Exercit¸iu 1.1.6 Mi¸scarea unui punct este descrisă de x1 = e t cos t, x2 =<br />
e t sin t, t ∈ R. Determinat¸i vectorii accelerat¸ie radială¸si transvesală ai punctului.<br />
Solut¸ie. Trebuie să introducem sistemul de coordonate polare (ρ, θ), astfel<br />
ca x1 = ρ cos θ, x2 = ρ sin θ. Luând în considerare ecuat¸ia de mi¸scare, deducem<br />
că<br />
<br />
ρ(t) = x2 1 + x22 = et , θ(t) = arctan x2<br />
= t.<br />
Mai mult, avem<br />
¸si<br />
x1<br />
r = cos ti1 + sin ti2, h = − sin ti1 + cos ti2,<br />
aρ = (¨ρ − ρ ˙ θ 2 )r = 0, aθ = 1 d<br />
ρ dt<br />
<br />
ρ 2 <br />
θ˙<br />
h = 2e t h.<br />
Exercit¸iu 1.1.7 Determinat¸i traiectoria punctului P care se mi¸scă într-un<br />
plan cu mărimea vitezei constante, ¸si astfel încât mărimea vitezei radiale fat¸ă<br />
de punctul O este de asemenea constantă.<br />
Solut¸ie. Introducem coordonatele polare (ρ, θ) în planul considerat. Atunci,<br />
viteza este dată de formula v = ˙ρr + ρ ˙ θh. Luând în considerare ipotezele<br />
problemei, obt¸inem<br />
˙ρ 2 + ρ 2 ˙ θ 2 = c1, ˙ρ = c2,<br />
unde constantele c1 ¸sic2 îndeplinesc în mod evident c1 > c 2 2. Astfel, avem<br />
ρ = c2t + ρ0, unde ρ0 = ρ(0) ¸si prin urmare, rezultă din ˙ρ 2 + ρ 2 ˙ θ 2 = c1 că
1.1. CINEMATICA PUNCTULUI MATERIAL 15<br />
P<br />
O<br />
Figura 1.6:<br />
x 2<br />
θ<br />
P 0<br />
ρ ˙ θ = k, unde k = ± c1 − c 2 2 este constant. Atunci, avem ˙ θ = k<br />
c2t+ρ0<br />
consecint¸ă, dacă θ(0) = 0, deducem că<br />
θ = k<br />
c2<br />
Rezultă din ultima expresie că<br />
log(c2τ + ρ0)| t 0 = k<br />
ρ = ρ0 exp( c2<br />
k θ),<br />
¸si traiectoria este spirala logaritmică (Figura 1.6).<br />
c2<br />
x 1<br />
log ρ<br />
.<br />
ρ0<br />
¸si în<br />
2 θ p<br />
Exercit¸iu 1.1.8 Traiectoria unei mi¸scări este parabola ρ cos 2 = 2 , p > 0.<br />
Un punct P se mi¸scă pe această parabolă asfel încât v = kρ, unde k este o<br />
constantă pozitivă. La momentul t = 0 punctul este în vârful parabolei ¸si se<br />
mi¸scă în sensul în care θ cre¸ste. Determinat¸i ecuat¸iile de mi¸scare ¸si vectorii<br />
accelerat¸ie radială ¸si transversală.<br />
Solut¸ie. Avem următoarele condit¸ii init¸iale:<br />
ρ(0) = p<br />
, θ(0) = 0,<br />
2<br />
¸si, în plus, θ(t) ˙ 2 2 > 0. Din relat¸ia v = kρ, deducem că ˙ρ + ρ θ˙ 2 2 2 = k ρ . În<br />
continuare vom determina ˙ρ ¸si θ. ˙ Pentru aceasta, derivăm ecuat¸ia parabolei<br />
pentru a obt¸ine ˙ρ cos θ<br />
2 − ρ ˙ θ sin θ<br />
2 = 0. Astfel, din aceste două relat¸ii de mai sus,<br />
deducem că<br />
˙ρ = ±kρ sin θ<br />
2 , θ ˙<br />
θ<br />
= ±k cos<br />
2 ,
16 CAPITOLUL 1. CINEMATICA<br />
din care, luând în considerare pozitivitatea lui k ¸si a lui ˙ θ(t), obt¸inem<br />
˙ρ = kρ sin θ<br />
2 , ˙ θ = k cos θ<br />
2 .<br />
Prin integrare, din ecuat¸ia diferent¸ială ˙ θ = k cos θ<br />
2 , obt¸inem<br />
<br />
<br />
ln <br />
tan <br />
θ π <br />
+ + c =<br />
4 4<br />
k<br />
t, c = constant,<br />
2<br />
¸si deci, din condit¸iile init¸iale θ(0) = 0, obi¸ntem c = 0 ¸si prin urmare<br />
Deoarece<br />
găsim<br />
tan θ<br />
4<br />
cos θ<br />
2<br />
cos θ<br />
2 =<br />
kt<br />
e 2 − 1<br />
=<br />
e kt<br />
2 + 1<br />
θ 1 − tan2 4<br />
=<br />
2 θ 1 + tan<br />
1<br />
cosh kt<br />
2<br />
Dacă substituim ˙ρ = kρ sin θ<br />
2<br />
dρ<br />
ρ<br />
kt<br />
kt<br />
e 4 − − e 4<br />
=<br />
4<br />
e kt<br />
kt<br />
4 + e− 4<br />
, sin θ<br />
2<br />
= tanh kt<br />
4 .<br />
θ 2 tan 4 =<br />
1 + tan<br />
2 θ<br />
4<br />
, sin θ<br />
2 =<br />
<br />
θ kt<br />
1 − cos2 = tanh<br />
2 2 .<br />
în această relat¸ie, obt¸inem<br />
kt<br />
p<br />
= k tanh dt, ρ(0) =<br />
2 2 .<br />
Astfel, obt¸inem următoarele ecuat¸ii de mi¸scare:<br />
Din aceste relat¸ii obt¸inem<br />
ρ = p<br />
2 cosh2 kt<br />
, θ = 2 arccos<br />
2<br />
1<br />
cosh kt<br />
2<br />
aρ = (¨ρ − ρ ˙ θ 2 )r = k2p [cosh (kt) − 2] r<br />
4<br />
aθ = 1 d<br />
<br />
ρ<br />
ρ dt<br />
2 <br />
θ˙<br />
h = 3k2 <br />
p kt<br />
sinh h.<br />
4 2<br />
Exercit¸iu 1.1.9 Punctul P se află într-o mi¸scare plană în care componenta<br />
radială a vitezei este direct proport¸ională cu timpul t ¸si componenta transversală<br />
este constantă. La momentul t = 0 punctul ocupă pozit¸ia P0(1, 0) fat¸ă de un<br />
sistem de referint¸ă. Să se determine traiectoria unui punct ¸si vectorii accelerat¸ie<br />
radială ¸si transversală.<br />
Solut¸ie. Alegem sistemul de coodonate polare (ρ, θ) cu polul în originea<br />
sistemului ¸si axa polară să coincidă cu axa x1. Din ipoteze avem că<br />
˙ρ = 2c 2 1t, ρ ˙ θ = c2,<br />
.<br />
,
1.1. CINEMATICA PUNCTULUI MATERIAL 17<br />
O<br />
Figura 1.7:<br />
O 1<br />
P(t)<br />
P(t + ∆t)<br />
P(t*)<br />
unde c1 ¸si c2 sunt constante pozitive prescrise. Notăm că avem următoarele<br />
condit¸ii init¸iale: ρ(0) = 1, θ(0) = 0. Atunci, prin integrare, obt¸inem ρ = c 2 1t 2 +c,<br />
c = constant, ¸si deci, din condit¸iile init¸iale ρ(0) = 1, avem ρ = c 2 1t 2 + 1.<br />
Apoi, avem ˙ θ = c2<br />
ρ<br />
= c2<br />
c 2 1 t2 +1<br />
¸si deci θ = c2<br />
c1 arctan (c1t) + c ∗ , c ∗ = constant.<br />
Din condit¸iile init¸iale θ(0) = 0, obt¸inem c∗ = 0 ¸si deci θ(t) = c2<br />
arctan (c1t).<br />
c1<br />
Eliminând parametrul t din relat¸iile ρ = c2 1t2 +1, θ(t) = c2<br />
c1 arctan (c1t), deducem<br />
ecuat¸ia traiectoriei ρ = 1 + tan2 <br />
.<br />
c1<br />
c2 θ<br />
Accelerat¸iile radială ¸si transversală sunt<br />
1.1.5 Viteza areolară<br />
aρ = 2c41t 2 + 2c2 1 − c2 2<br />
c2 1t2 r, aθ =<br />
+ 1<br />
2c21c2t c2 1t2 + 1 h.<br />
Pentru o mi¸scare plană, introducem not¸iunea de viteză areolară. Dacă un punct<br />
O1 este fixat pe traiectorie, notăm cu A(t) aria măturată de raza vectoare<br />
(P − O), aceasta este aria regiunii delimitate de vectorii (O1 − O), (P (t) − O)<br />
¸si arcul O1P (t) al traiectoriei, unde O este originea sistemului de coordonate<br />
(ρ, θ) (Figura 1.7).<br />
Definit¸ie 1.1.8 Numim viteză areolară ˙<br />
A a punctului P fact¸ă de polul O derivata<br />
funct¸iei A(t) în raport cu timpul.<br />
Rezultă din definit¸ia vitezei areolare că<br />
˙A<br />
A(t + ∆t) − A(t)<br />
= lim<br />
= lim<br />
∆t→0 ∆t<br />
∆t→0<br />
∆A<br />
. (1.35)<br />
∆t<br />
Este u¸sor de demonstrat că aria măturată între momentele t ¸si t + ∆t este dată<br />
de formula<br />
∆A = 1<br />
2 ρ2 (t ∗ )∆θ, (1.36)
18 CAPITOLUL 1. CINEMATICA<br />
unde t ∗ ∈ [t, t + ∆t] ¸si ∆θ = θ(t + ∆t) − θ(t). Cu alte cuvinte, există un moment<br />
t ∗ astfel ca aria ∆A este egală cu aria sectorului circular cu unghiul la centru<br />
∆θ ¸si raza ρ(t ∗ ). Astfel, din (1.35) ¸si (1.36), obt¸inem<br />
˙A(t) = 1<br />
2 ρ2 (t) ˙ θ(t). (1.37)<br />
În coordonate carteziene, deoarece x1 = ρ cos θ, x2 = ρ sin θ, avem<br />
˙x1 = ˙ρ cos θ − ρ ˙ θ sin θ, ˙x2 = ˙ρ sin θ + ρ ˙ θ cos θ,<br />
x1 ˙x2 − x2 ˙x1 = ρ 2 ˙ θ cos 2 θ + ρ ˙ρ sin θ cos θ − ρ ˙ρ sin θ cos θ + ρ 2 ˙ θ sin 2 θ<br />
= ρ 2 ˙ θ,<br />
¸si deci ecuat¸ia vitezei areolare poate fi scrisă ca<br />
˙A = 1<br />
2 (x1 ˙x2 − x2 ˙x1). (1.38)<br />
Exercit¸iu 1.1.10 Mi¸scarea unui punct P pe suprafat¸a plană (O, x1, x2) este<br />
dată de ecuat¸iile carteziane<br />
x1 = C exp(−pt), x2 = C exp(pt), C > 0, p > 0.<br />
Să se determine traiectoria, viteza areolară fat¸ă de O, ¸si componentele radială<br />
¸si transversală a vectorului accelerat¸ie.<br />
Solut¸ie. Eliminând timpul t din ecuat¸iile de mi¸scare, obt¸inem<br />
x1 · x2 = C 2 .<br />
Prin urmare, traiectoria este o ramură a unei hiperbolei (Figura 1.8) situată în<br />
primul cadran al sistemului de coordonate. Mai mult, viteza areolară este dată<br />
de formula<br />
˙A = 1<br />
2 (x1 ˙x2 − x2 ˙x1)<br />
= 1 2 2 2<br />
C p exp(−pt) exp(pt) + C p exp(−pt) exp(pt) = C p.<br />
2<br />
Deoarece viteza areolară este constantă, componenta transversală a accelerat¸iei<br />
este aθ = 0, în timp ce componenta radială dă accelerat¸ia totală ¸si deci<br />
<br />
aρ = a = ¨x 2 1 + ¨x2 2 = Cp2exp(−2pt) + exp(2pt) = p 2 ρ,<br />
unde ρ este distant¸a dintre P ¸si O, care este dată de formula<br />
<br />
ρ = x2 1 + x22 = Cexp(−2pt) + exp(2pt).
1.1. CINEMATICA PUNCTULUI MATERIAL 19<br />
O<br />
x 2<br />
Figura 1.8:<br />
1.1.6 Mi¸scări centrale<br />
Considerăm o mi¸scare care este nu este neapărat plană.<br />
Definit¸ie 1.1.9 Mi¸scarea unui punct P este numită centrală dacă accelerat¸ia<br />
sa este întotdeauna direct¸ionată de-a lungul vectorului P − O, unde O este un<br />
punct fixat numit centrul mi¸scării.<br />
Teoremă 1.1.1 Orice mi¸scare centrală cu centrul O este plană ¸si viteza areolară<br />
fat¸ă de O este constanta, ¸si vice versa.<br />
Demonstrat¸ie. Din definit¸ia mi¸scării centrale, obt¸inem<br />
Din ultima egalitate, rezultă că<br />
¸si deci<br />
x 1<br />
a(t) × (P (t) − O) = 0 pentru orice t.<br />
d<br />
d(P − O)<br />
[v × (P − O)] − v × =<br />
dt dt<br />
d<br />
[v × (P − O)] = 0,<br />
dt<br />
v × (P − O) = k, (1.39)<br />
unde k este un vector constant. Presupunem că k = 0, ¸si apoi, din (1.39),<br />
obt¸inem<br />
0 = v × (P − O) · (P − O) = k · (P − O).<br />
Prin urmare, punctul P trebuie că rămână în planul ortogonal la k ¸si care trece<br />
prin punctul O. Dacă k = 0, atunci<br />
v × (P − O) = 0 pentru orice t,
20 CAPITOLUL 1. CINEMATICA<br />
deci v ¸si a sunt întotdeauna paralelei, ¸si de asemenea ambii sunt paraleli cu<br />
(P − O). Ultima implică că an = ˙s2<br />
1<br />
ρ = 0 pentru orice t, ¸si deoarece ρ = 0,<br />
mi¸sarea este rectilinie.<br />
Prin urmare, mi¸scarea centrala este una plană. Astfel, putem să o reprezentăm<br />
în coordonare polare cu polul O. Mai mult, vectorul accelerat¸ie este radial,<br />
deoarece are aceea¸si direct¸ie cu (P − O), ¸si va implica că aθ = 1<br />
<br />
d<br />
ρ dt ρ2 ˙ <br />
θ = 0.<br />
Prin urmare,<br />
ρ 2 θ ˙ = c (1.40)<br />
implică că viteza areolară a mi¸scării lui P fat¸ă de O este constantă ¸si valoarea<br />
sa este dată de formula<br />
˙A = c<br />
,<br />
2<br />
(1.41)<br />
unde c este numită constanta ariilor.<br />
Să demonstrăm acum că, dacă viteza areolară fat¸ă de polul O pentru o<br />
mi¸scare plană este constantă, atunci mi¸scarea este centrală. Într-adevăr, deoarece<br />
aθ = 0, faptul că viteza areolară este constantă implică că accelerat¸ia a = aρ<br />
este mereu îndreptată spre O.<br />
Teoremă 1.1.2 Pentru o mi¸scare centrală având constanta ariilor c, accelerat¸ia<br />
a poate fi determinată, cunoscând doar traiectoria punctului (ρ = ˆρ (θ)), prin<br />
intermediul formulei lui Binet<br />
a = − c2<br />
ρ2 2 d<br />
dθ2 <br />
1<br />
+<br />
ρ<br />
1<br />
<br />
r. (1.42)<br />
ρ<br />
Demonstrat¸ie. Fat¸ă de un sistem de coordonate polare având originea<br />
în O, ecuat¸ia traiectoriei este ρ = ˆρ(θ). Dacă mi¸scarea este centrală, viteza<br />
areolară este constantă ¸si prin urmare avem, ρ 2 ˙ θ = c; unde acceleralt¸ia este<br />
Pe de altă parte, avem<br />
¸si deci<br />
a = aρr = (¨ρ − ρ ˙ θ 2 )r. (1.43)<br />
˙ρ = dρ<br />
dθ ˙ θ = c<br />
ρ2 dρ d<br />
= −c<br />
dθ dθ<br />
¨ρ = −c d2<br />
dθ 2<br />
<br />
1<br />
˙θ = −<br />
ρ<br />
c2<br />
ρ2 d 2<br />
dθ 2<br />
<br />
1<br />
, (1.44)<br />
ρ<br />
<br />
1<br />
. (1.45)<br />
ρ<br />
Mai mult, folosind relat¸iile (1.40) ¸si (1.45), din relat¸iile (1.43) deducem formula<br />
lui Binet<br />
aρ = − c2<br />
ρ2 2 d<br />
dθ2 <br />
1<br />
+<br />
ρ<br />
1<br />
<br />
,<br />
ρ<br />
care ne permite să determinăm accelerat¸ia folosind ecuat¸ia traiectoriei ρ = ˆρ(θ)<br />
¸si presupunând cunoscută constanta ariilor c.
1.1. CINEMATICA PUNCTULUI MATERIAL 21<br />
Exercit¸iu 1.1.11 Punctul P descrie o curbă plană astfel încât accelerat¸ia sa<br />
trece mereu printr-un punct fix O. Demonstrat¸i că<br />
a = v dv<br />
dρ ,<br />
unde v = |v| ¸si a este componenta radială a accelerat¸iei.<br />
Solut¸ie. Mi¸scarea unui punct P este centrală. Folosind sistemul polar de<br />
coordonate cu polul în O, avem ρ2 ˙ θ = c, unde c este constanta ariilor. Astfel,<br />
avem aθ = 1<br />
<br />
d ρ2 ˙ <br />
θ h = 0 ¸si a = aρ = (¨ρ − ρ ˙ θ2 )r. Pe de altă parte, avem<br />
ρ dt<br />
v 2 = ˙ρ 2 + ρ 2 ˙ θ 2 ,<br />
¸si deci, prin derivare directă în raport cu t, obt¸inem<br />
Astfel, obt¸inem<br />
2v dv<br />
dt = 2 ˙ρ¨ρ + 2ρ ˙ρ ˙ θ 2 + 2ρ 2 <br />
θ˙ θ ¨ = 2 ˙ρ ¨ρ − ρ ˙ θ 2<br />
+ 2ρ ˙ <br />
θ ρ¨ θ + 2 ˙ρ ˙ <br />
θ .<br />
¸si deci relat¸ia cerută.<br />
v dv<br />
dt<br />
= dρ<br />
dt a,<br />
1.1.7 Mi¸scări uniform variate ¸si periodice<br />
Numim uniformă orice mi¸scare a cărui viteză este constanta în timp; o astfel de<br />
definit¸ie nu depinde de traiectoria punctului. Prin urmare, notând cu v0 = ˙s(t)<br />
această valoare constantă, ecuat¸ia orară devine<br />
s(t) = v0t + s0, (1.46)<br />
unde cei doi parametri s0 ¸si v0 reprezintă abscisa curbilinie init¸ială ¸si, respectiv,<br />
viteza punctului P .<br />
Să considerăm o mi¸scare care nu este în mod necesar rectilinie.<br />
Definit¸ie 1.1.10 Mi¸scarea unui punct P se nume¸ste uniform variată dacă mărimea<br />
accelerat¸iei tangent¸iale este constantă, adică, există o constantă a0 astfel ca<br />
¨s(t) = a0.<br />
Prin urmare, prin integrarea ultimei relat¸ii de două ori în raport cu timpul<br />
t, obt¸inem următoarea ecuat¸ie orară pentru mi¸scarea uniform variată:<br />
s(t) = 1<br />
2 a0t 2 + v0t + s0, (1.47)<br />
unde s0 ¸si v0 reprezintă abcisa curbilinie ¸si, respectiv viteza la momentul t = 0.<br />
Este evident din (1.47) că ecuat¸ia orară pentru mi¸scarea uniform variată este<br />
reprezentată grafic ca o parabolă care este concavă pentru a0 < 0, ¸si convexă
22 CAPITOLUL 1. CINEMATICA<br />
pentru a0 > 0. Astfel, independent de concavitatea sau convexitatea parabolei,<br />
mi¸scarea în direct¸ia de cre¸stere a arcului parabolei este numită directă ¸si cea în<br />
direct¸ia de descre¸stere a arcului parabolei este numită retrogradă.<br />
Înainte de a considera mi¸scarea circulară ¸si uniformă, explicăm ce întelegem<br />
prin mi¸scare periodică a punctului P care se mi¸scă pe o traiectorie asociată.<br />
Definit¸ie 1.1.11 Spunem că mi¸scarea unui punct t P este periodică cu perioda<br />
T dacă ecuat¸ia orară ˆs(t) define¸ste o funct¸ie periodică de t cu perioada T , adică<br />
ˆs(t + T ) = ˆs(t). (1.48)<br />
Observat¸ie 1.1.3 Dacă mi¸scarea este periodică, atunci viteza ¸si accelerat¸ia<br />
scalară ( 7 ) sunt periodice în t.<br />
1.1.8 Mi¸scări circulare ¸si uniforme<br />
Mi¸scarea circulară este o mi¸scarea plană particulară definită astfel:<br />
Definit¸ie 1.1.12 Mi¸scarea unui punct P este numită circulară dacă traiectoria<br />
sa este un cerc sau un arc de cerc. În plus, dacă viteza este constantă, atunci<br />
este numită circular ¸si uniformă.<br />
Considerăm o mi¸scarea circulară relativ la un cerc de rază R (Figura 1.9).<br />
Dacă notăm cu s abscisa curbilie astfel încât 0 ≤ s ≤ 2πR, ¸si dacă presupunem<br />
că s = ˆs(t) este ecuat¸ia orară corespunzătoare, atunci vectorii viteza<br />
¸si accelerat¸ie sunt date de formulele (1.16) ¸si (1.23 ), adică<br />
unde R este raza cercului. Deoarece n = −<br />
ecuat¸ie din (1.49) poate fi rescrisă ca<br />
v = ˙st, a = ¨st + ˙s2<br />
n, (1.49)<br />
R<br />
P −O<br />
|P −O|<br />
P −O = − R , ce-a de a doua<br />
a = ¨st − ˙s2<br />
(P − O). (1.50)<br />
R2 Teoremă 1.1.3 Mi¸scarea circulară uniformă reprezintă un exemplu important<br />
de mi¸scare periodică. Dacă ˙s = v0, atunci perioada unei astfel de mi¸scări este<br />
T = 2πR<br />
. (1.51)<br />
Mai mult, accelerat¸ia este centripetă ¸si este dată de formula a = −ω 2 (P − O),<br />
unde ω = v0/R.<br />
v0<br />
7 Prin accelerat¸ie scalară, înt¸elegem componenta accelerat¸iei directe de-a lungul tangentei<br />
la curbă.
1.1. CINEMATICA PUNCTULUI MATERIAL 23<br />
O<br />
x 2<br />
n<br />
Figura 1.9:<br />
r<br />
Demonstrat¸ie. Avem s(t) = Rθ(t) + s0 ¸si prin urmare obt¸inem<br />
P<br />
s<br />
O 1<br />
x 1<br />
v = ˙st = R ˙ θt. (1.52)<br />
Deoarece ˙s este constantă, rezultă că ˙ θ este constantă ¸si astfel, alegând ˙ θ = ω,<br />
obt¸inem<br />
ˆθ(t) = ωt + θ0, (1.53)<br />
unde θ0 este valoarea unghiului θ la momentul t = 0. Din (1.52) deducem că<br />
v0 = R ˙ θ = Rω,<br />
¸si deci ω = v0/R. Urmează din (1.53) că funct¸ia ˆ θ satisface relat¸ia<br />
ˆθ(t + 2π<br />
ω ) = ˆ θ(t) + 2π,<br />
¸si prin urmare mi¸scarea este periodică cu perioada 2π/ω (a se vedea Figura 1.9).<br />
Prin urmare, deoarece ω = v0/R, rezultă că perioada mi¸scării circulare este<br />
T = 2π<br />
ω<br />
2πR<br />
= . (1.54)<br />
v0<br />
Inversa acestei perioade este numită frecvent¸ă ν = 1 ω<br />
T = 2π .<br />
Deoarece t = 1<br />
Rk × (P − O), unde k este vector unitar, ortogonal cercului ¸si<br />
direct¸ionat astfel ca t, k, (P −O) să formeze un triplet drept, din (1.52) obt¸inem<br />
v = ˙ θk × (P − O) = ω × (P − O), (1.55)
24 CAPITOLUL 1. CINEMATICA<br />
unde ω = ˙ θk este numită viteză unghiulară.<br />
Expresia (1.55) pentru viteza lui P în termenii vectorului ω poate fi de<br />
asemenea scrisă folosind matricea antisimetrică W = (Whk) legată de vectorul<br />
ω = (ω1, ω2, ω3) ¸si definită astfel<br />
⎛<br />
W = ⎝ 0 −ω3<br />
⎞<br />
ω2<br />
ω3 0 −ω1 ⎠ .<br />
−ω2 ω1 0<br />
Este u¸sor de verificat ( 8 ) că , dacă x = (x1, x2, x3), atunci ω×(P −O) = Wx,<br />
¸si deci<br />
v = Wx.<br />
Mai mult, deoarece<br />
ˆs(t) = R ˆ θ(t) = R(ωt + θ0),<br />
rezultă că funct¸ia ˆs este de asemenea periodică cu perioada T = 2πR , ¸si deci<br />
v0<br />
mi¸scarea este periodică cu aceea¸si perioada T .<br />
În final, deoarece mi¸scarea este uniformă (¨s = 0), accelerat¸ia este centripetă<br />
Ecuat¸iile carteziene a mi¸scării circulare sunt<br />
a = v2 0<br />
R n = −ω2 (P − O). (1.56)<br />
x1 = R cos θ, x2 = R sin θ, θ = ˆ θ(t).<br />
Dacă mi¸scarea este uniformă, atunci ˆ θ(t) = ωt + θ0, ¸si astfel<br />
x1 = R cos(ωt + θ0), x2 = R sin(ωt + θ0).<br />
Exercit¸iu 1.1.12 Mi¸scarea unui punct P este descrisă de x(t) = 3 cos ωti1 +<br />
3 sin ωti2, unde ω este o constantă prescrisă. Să se demonstreze că mi¸scarea<br />
este centrală. Calculat¸i x · v ¸si x × v.<br />
Solut¸ie. Avem v = −3ω sin ωti1 + 3ω cos ωti2 and a = −3ω 2 cos ωti1 −<br />
3ω 2 sin ωti2. Apoi, avem a = −ω 2 x ¸si deci mi¸scarea este centrală. Obt¸inem<br />
x·v = −9ω sin ωt cos ωt+9ω sin ωt cos ωt = 0 ¸si x×v = (3 cos ωti1 +3 sin ωti2)×<br />
(−3ω sin ωti1 + 3ω cos ωti2) = 9ωi3.<br />
Exercit¸iu 1.1.13 Un punct P se mi¸scă pe un cerc a cărui rază este R cu<br />
accelerat¸ia tangent¸ială constantă at. Punctul P porne¸ste din P0 la momentul<br />
t = 0. Determinat¸i intervalul de timp în care accelerat¸ia centripetă an devine<br />
egală cu accelerat¸ia tangent¸ială at.<br />
Solut¸ie. Avem at = ¨s ¸si an = ˙s2<br />
ρ<br />
= ˙s2<br />
R . Astfel, deducem că ˙s = att + c ¸si,<br />
din condit¸ia init¸ială ˙s(0) = 0, obt¸inem ˙s = att. Avem an = at unde ˙s2<br />
deci (att) 2 <br />
R<br />
= Rat. Prin urmare, intervalul cerut este t = at .<br />
8 A se vedea Teorema D.2 din Appendix D.<br />
R = at ¸si
1.1. CINEMATICA PUNCTULUI MATERIAL 25<br />
Exercit¸iu 1.1.14 Un punct se mi¸scă pe cercul de raza R după următoarea<br />
ecuat¸ie orară s = v0t − c<br />
2t2 , unde v0 ¸si c sunt constante. Să se determine<br />
mărimea accelerat¸iei.<br />
Solut¸ie. Avem ˙s = v0 − ct ¸si ¨s = −c, ¸si astfel obt¸inem<br />
Deci, avem<br />
a = ¨st + ˙s2<br />
ρ n = −ct + (v0 − ct) 2<br />
n.<br />
R<br />
a =<br />
<br />
c2 + (v0 − ct) 4<br />
R2 .<br />
Exercit¸iu 1.1.15 Un punct se mi¸scă pe un cerc de rază R cu accelerat¸ia unghiulară<br />
constantă α. La momentul t = 0, punctul porne¸ste din repaus. Să se<br />
demonstreze că la momentul t viteza unghiulară este ω = αt ¸si că a parcurs<br />
lungimea de arc s = 1<br />
2 Rαt2 .<br />
Solut¸ie. Deoarece ¨ θ = α, rezultă că θ = 1<br />
2αt2 + θ1t + θ0, unde θ0, θ1 sunt<br />
constante. Luând în considerare condit¸iile init¸iale, deducem că θ1 = 0 ¸si deci<br />
θ − θ0 = 1<br />
2αt2 . Apoi, avem ω = αt ¸si s = R[θ(t) − θ0] = 1<br />
2Rαt2 .<br />
1.1.9 Mi¸scări armonice<br />
Începem cu studiul unei miscari circulare uniforme a unui punct P pe un cerc<br />
de centru O ¸si rază R. Notăm cu P ∗ proiect¸ia lui P pe un diametru fixat<br />
AB. Atunci, în timp ce P descrie cercul, P ∗ se mi¸scă pe diametrul AB după<br />
următoarea lege (Figura 1.10):<br />
x = R cos( ˙ θt + θ0),<br />
unde x este componenta lui P −O de-a lungul diametrului AB, ¸si θ este unghiul<br />
P OB. În final, θ0 este valoarea lui θ la t = 0. Deoarece mi¸scarea este uniformă,<br />
˙θ = ω este constant ¸si avem<br />
x = R cos(ωt + θ0), (1.57)<br />
¨x = −Rω 2 cos(ωt + θ0). (1.58)<br />
Definit¸ie 1.1.13 O mi¸scare rectilinie este numită oscilat¸ie armonică dacă ecuat¸ia<br />
orară este dată de<br />
s(t) = C cos(ωt + γ), (1.59)<br />
unde constantele C, ω ¸si γ sunt numite amplitudine, pulsat¸ie (sau frecvent¸ă<br />
unghiulară) ¸si fază.<br />
Rezultă din (1.57) că mi¸scarea lui P ∗ de-a lungul diametrului AB este armonică.<br />
În plus, din (1.57) ¸si (1.58), obt¸inem proprietatea importantă descrisă<br />
de ecuat¸ia<br />
¨s = −ω 2 s, (1.60)
26 CAPITOLUL 1. CINEMATICA<br />
A O x P* B x<br />
Figura 1.10:<br />
θ<br />
P<br />
adică, într-o mi¸scare armonică, accelerat¸ia scalară ¨s este proport¸ională cu distant¸a<br />
parcursă s, are semn opus ¸si coeficient¸ul său de proport¸ionalitate este egal cu<br />
pătratul frecvent¸ei unghiulare ω.<br />
Notăm că expresia (1.60) nu cont¸ine nici amplitudinea, nici faza init¸ială a<br />
mi¸scării armonice. Într-adevăr, avem următorul rezultat:<br />
Teoremă 1.1.4 Orice mi¸scare armonică de frecvent¸ă unghiulară ω (cu amplitudine<br />
a ¸si fază arbitrară) satisface ecuat¸ia diferent¸ială (1.60), ¸si vice versa.<br />
Demontrat¸ie. Dacă A este amplitudinea ¸si γ este faza init¸ială a unei<br />
mi¸scări armonice date<br />
s(t) = A cos(ωt + γ),<br />
atunci, în baza relat¸iilor (1.57) ¸si (1.58), satisface ecuat¸ia (1.60).<br />
Vice versa, dată ecuat¸ia diferent¸ială (1.60), concluzionăm că ecuat¸ia caracteristică<br />
este<br />
λ 2 + ω 2 = 0,<br />
a cărui solut¸ii sunt λ1 = −iω, λ2 = iω; astfel, solut¸ia generală este dată de<br />
formula<br />
s(t) = C1 cos ωt + C2 sin ωt. (1.61)<br />
Alegând două constante A ¸si γ astfel ca<br />
din (1.61) obt¸inem<br />
C1 = A cos γ, C2 = −A sin γ,<br />
s(t) = A cos ωt cos γ − A sin ωt sin γ = A cos(ωt + γ).<br />
Astfel, ecuat¸ia (1.60) este caracteristică mi¸scărilor armonice cu frecvent¸a unghiulară<br />
ω.
1.1. CINEMATICA PUNCTULUI MATERIAL 27<br />
Observat¸ie 1.1.4 O mi¸scare armonică cu frecvent¸a unghiulară ω are aceea¸si<br />
perioadă cu mi¸scarea circulară uniformă, adică, T = 2π/ω; în timp ce amplitudinea<br />
oscilat¸iei coincide cu rasa cercului, ¸si faza coincide cu valoarea unghiului<br />
θ la t = 0.<br />
Exercit¸iu 1.1.16 Un punct P are o mi¸scare oscilatorie armonică descrisă de<br />
ecuat¸ia<br />
<br />
2π<br />
x = A sin<br />
T t<br />
<br />
.<br />
Pentru x = x1 viteza punctului este v1, în timp ce pentru x = x2 viteza este<br />
v2. Să se determine amplitudinea A ¸si perioada T a mi¸scării oscilatorii a<br />
punctului P .<br />
Solut¸ie. Viteza punctului P este v = A 2π<br />
T cos 2π<br />
T t . Atunci, din ipoteză,<br />
avem<br />
<br />
2π<br />
x1 = A sin<br />
T t1<br />
<br />
, v1 = A 2π<br />
T cos<br />
<br />
2π<br />
T t1<br />
<br />
,<br />
Eliminând t1 ¸si t2, obt¸inem<br />
din care deducem<br />
<br />
2π<br />
x2 = A sin<br />
T t2<br />
<br />
, v2 = A 2π<br />
T cos<br />
<br />
2π<br />
T t2<br />
<br />
.<br />
x 2 1 +<br />
A =<br />
2 T<br />
4π2 v2 1 = A 2 , x 2 2 T<br />
2 +<br />
4π2 v2 2 = A 2 ,<br />
<br />
x2 1v2 2 − x22 v2 1<br />
v2 2 − v2 , T = 2π<br />
1<br />
<br />
x2 1 − x22 v2 2 − v2 .<br />
1<br />
Exercit¸iu 1.1.17 Legea de mi¸scare a unui lift este x = H (1 − cos ϕ), unde H<br />
2<br />
2k<br />
H<br />
este cea mai mare înălt¸ime la care ajunge liftul ¸si ϕ = t, k = constant. Să<br />
se determine viteza ¸si accelerat¸ia liftului. Determinat¸i timpul necesar liftului pe<br />
tru a ajunge la înălt¸imea H.<br />
Solut¸ie. Prin derivări succesive, obt¸inem<br />
v =<br />
kH<br />
2<br />
sin ϕ, a = k cos ϕ.<br />
În plus, pentru x = H, deducem H<br />
2 (1 − cos ϕ) = H ¸si deci ϕ = π. Astfel,<br />
<br />
2k<br />
H<br />
relat¸ia π = H t ne oferă timpul t = π 2k .
28 CAPITOLUL 1. CINEMATICA<br />
x 1<br />
O<br />
θ<br />
x 3<br />
P<br />
P<br />
P *<br />
Figura 1.11:<br />
x 2
1.1. CINEMATICA PUNCTULUI MATERIAL 29<br />
1.1.10 Mi¸scări elicoidale<br />
Considerăm un cilindru circular de rază R. Numim elice circulară o formă<br />
descrisă de o curbă care intersectează mereu generatoarea cilindrului sub un alt<br />
unghi (Figura 1.11).<br />
Definit¸ie 1.1.14 Mi¸scarea unui punct P pe o suprafat¸ă cilindrică este numită<br />
elicoidală dacă punctul P se mi¸scă pe cilindru după o elice.<br />
Alegem un sistem cartezian ortogonal (O, x1, x2, x3) astfel ca axa x3− să<br />
coincidă cu axa cilindrului. Apoi, notăm cu θ unghiul dintre proiect¸ia (P ∗ − O)<br />
a lui (P − O) pe planul x1Ox2 ¸si axa x1.<br />
Este posibil să reprezentăm mi¸scarea punctului folosind următoarea expresie<br />
a vectorului (P − O) :<br />
P − O = (P − P ∗ ) + (P ∗ − O). (1.62)<br />
Deoarece mi¸scarea punctului P ∗ este una circulară, alegând convenabil sistemul<br />
de referint¸ă (O, x1, x2, x3), obt¸inem<br />
P − O = R cos θi1 + R sin θi2 + hθi3,<br />
unde h este un parametru ales astfel ca |P − P ′ | = 2πh. Observă că |P − P ′ |<br />
reprezintă distant¸a dintre două puncte consecutive ale elicei, situate pe aceea¸si<br />
generatoare ¸si numită pasul elicei. Din ultima relat¸ie obt¸inem<br />
x1 = R cos θ, x2 = R sin θ, x3 = hθ.<br />
Acest sistem reprezintă (elicea) drumul lui P , în timp ce ecuat¸ia orară este dată<br />
în termenii lui θ, de funct¸ia θ = ˆ θ(t).<br />
Dacă mi¸scarea este uniformă, atunci avem ˙ θ = constant, ¸si mi¸scarea va fi<br />
numită elicoidală ¸si uniformă.<br />
Folosind formula (1.62), este posibil să obt¸inem următoarea descompunere<br />
a vitezei:<br />
v =<br />
d(P − O)<br />
dt<br />
= d(P − P ∗ )<br />
dt<br />
Deoarece mi¸scarea lui P ∗ este circulară, avem<br />
Astfel, din ω(t) = ˙ θ(t)i3, obt¸inem<br />
v = h ˙ θi3 + ˙ θi3 × (P ∗ − O).<br />
v = hω + ω × (P ∗ − O).<br />
+ d(P ∗ − O)<br />
.<br />
dt<br />
Ultima relat¸ie demonstrează că viteza are două componente, prima corespunde<br />
unei mi¸scări rectilinii de-a lungul axei x3− ¸si ce-a de a doua corespunde unei
30 CAPITOLUL 1. CINEMATICA<br />
mi¸scări circulare. Să definim acum vectorul tangent t ¸si normala principală n.<br />
Deoarece s = (R 2 + h 2 ) 1/2 θ, avem dθ<br />
ds = (R2 + h 2 ) −1/2 . Mai mult, obt¸inem<br />
t = dP dθ dθ<br />
=<br />
dθ ds<br />
n = ρ dt dt<br />
= ρ<br />
ds dθ<br />
ds (−R sin θi1 + R cos θi2 + hi3),<br />
<br />
dθ dθ<br />
= ρ<br />
ds ds<br />
2<br />
R(− cos θi1 − sin θi2).<br />
Deoarece n este un vector unitar, deducem din ultima relat¸ie că ρ = <br />
dθ −2 1<br />
ds R =<br />
R 2 +h 2<br />
R , ¸si în consecint¸ă avem n = − P ∗ −O<br />
R . Astfel, normala principală la curbă<br />
coincide cu normala la suprafat¸ă ¸si astfel elicile sunt geodezice ( 9 ) ale cilindrului.<br />
Trebuie punctat faptul că o descriere generală a mi¸scării folosind coordonatele<br />
curbilinii este prezentată în Appendix A.<br />
Exercit¸iu 1.1.18 Mi¸scarea unui punct este dată de x = a cos e −t i1+a sin e −t i2+<br />
be −t i3, unde a, b sunt constante pozitive. Să se determine componentele tangent¸ială<br />
¸si normală ale accelarat¸iei punctului.<br />
Solut¸ie. Avem o mi¸scare elicoidală. Prin derivare directă, avem dx =<br />
−e −t (−a sin e −t i1 + a cos e −t i2 + bi3) dt, ¸si deci ds = e −t√ a 2 + b 2 dt. Mai mult,<br />
avem<br />
¸si<br />
¸si<br />
t = dx<br />
ds =<br />
1 −t √ a sin e i1 − a cos e<br />
a2 + b2 −t <br />
i2 − bi3 ,<br />
dt<br />
ds<br />
Apoi, deducem că<br />
dt dt a<br />
= = −<br />
dt ds a2 + b2 −t<br />
cos e i1 + sin e −t <br />
i2 .<br />
v = ˙st = e −t a 2 + b 2 t,<br />
at = ¨st = −e −t a2 + b2 2 dt<br />
t, an = ˙s<br />
ds = −ae−2th, unde h = cos e −t i1 + sin e −t i2.<br />
1.2 Cinematica sistemelor materiale ¸si corpurilor<br />
rigide<br />
1.2.1 Legături ¸si sisteme olonome<br />
Să considerăm un sistem material B de N puncte materiale, care sunt notate cu<br />
P1, P2, . . ., PN. Dacă punctele sunt libere să ocupe pozit¸ii arbitrare din spat¸iu,<br />
atunci sistemul material este numit liber. Configurat¸ia sistemului material liber<br />
a N puncte dată într-un sistem de referint¸ă (O, x1, x2, x3) este cunoscută atunci<br />
9 Reamintim că geodezica la o suprafat¸ă este acea curbă de pe suprafat¸ă a cărui normală<br />
este direct¸ionată de-a lungul normalei la suprafat¸ă (a se vedea Appendix A, definit¸ia A.13).
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE31<br />
când sunt cunoscut¸i vectorii de pozit¸ie ai fiecărui punct relativ la un punct fixat<br />
în (O, x1, x2, x3). Astfel, N cantităt¸i vectoriale sau, echivalent, 3N cantităt¸i<br />
scalare sunt cerute pentru a specifica configurat¸ia sistemului material liber întrun<br />
sistem de referint¸ă fixat.<br />
Dacă, spre deosebire, mi¸scarea unui sistem material este afectată de prezent¸a<br />
corpurilor care vin în contact cu câteva dintre punctele lui B, legături pot fi<br />
impuse asupra pozit¸iilor pe care punctele materiale le pot ocupa sau asupra<br />
manierei în care aceste pozit¸ii se pot schimba. În acest caz, clasa P, reprezentând<br />
toate pozit¸iile posibile ale sistemului material B, nu este destul de largă pentru<br />
a permite corpului B să aibă o configurat¸i arbitrară în E. Se spune astfel<br />
că sistemul material B este supus la legături. Dacă, pornind de la cunoa¸sterea<br />
câtorva componente ale deplasărilor sistemului material, putem afirma ceva despre<br />
deplasările rămase, putem spune că această legătură este activă.<br />
Definit¸ie 1.2.1 Numim legătură orice mecanism care impune restrict¸ii privind<br />
pozit¸ia ¸si viteza ale celor N puncte care formează sistemul material. Aceste<br />
restrict¸ii pot fi exprimate analitic prin intermediul unei relat¸ii între coordonatele<br />
¸si vitezele punctelor sistemului de material în forma<br />
ψ (x1, x2, . . . , xN, ˙x1, ˙x2, . . . , ˙xN, t) ≥ 0. (1.63)<br />
În relat¸iile de mai sus, (xi, ˙xi) reprezintă pozit¸ia ¸si viteza punctului Pi, i =<br />
1, 2, ..., N. Mai mult, noi vom presupune ulterioe că funct¸ia ψ este suficient de<br />
regulată.<br />
Ca un prim exemplu de sistem material care este supus legăturilor, putem<br />
considera un sistem rigid, care este un sistem material P1, P2, ..., PN pentru<br />
care distant¸ele dintre punctele rămân invariabile în raport cu timpul, adică<br />
d(Ph(t), Pk(t)) = dhk, h, k = 1, 2, ..., N,<br />
unde d reprezintă distant¸a dintre puncte ¸si dhk sunt independente de timp.<br />
Un alt exemplu de sistem constrâns poate fi găsit în studiul mi¸scării unui<br />
punct P fort¸at să se mi¸ste pe un cerc. De fapt, dacă O este centrul cercului de<br />
rază R, avem<br />
(P − O) 2 = R 2 . (1.64)<br />
Definit¸ie 1.2.2 Spunem că o legătură este bilaterală când restrict¸iile sistemului<br />
materil pot fi reprezentate de o relat¸ie de tipul (1.63) dar cu egalitate.<br />
Spunem că avem o legătură unilaterală când relat¸ia ce o descrie este o inegalitate.<br />
Exemplele de mai sus reprezentând un sistem rigid ¸si un punct ce se mi¸scă<br />
pe cerc descriu legături bilaterale. Un exemplu de legătură unilaterlă este acela<br />
a unui punct fort¸at să rămână într-un plan sau cel a unui punct constrâns să<br />
rămână în interiorul unei sfere.
32 CAPITOLUL 1. CINEMATICA<br />
O<br />
x 2<br />
R(t)<br />
Figura 1.12:<br />
Definit¸ie 1.2.3 Spunem că o legătură este scleronomă sau independentă de<br />
timp dacă relat¸ia care descrie legătură nu cont¸ine timpul în mod explicit. O<br />
legătură este reonomă sau dependentă de timp dacă relat¸ia care descrie legătura<br />
depinde explicit de timp.<br />
Un exemplu de legătură reonomă este descrisă de Figura 1.12 de un punct<br />
constrâns să rămână pe un cerc de rază R(t), variabilă în timp, adică este<br />
legătura reprezentată de<br />
P<br />
x 1<br />
(P − O) 2 = x 2 1 + x 2 2 = R 2 (t).<br />
Definit¸ie 1.2.4 O legătură este numită olonomă sau geometrică sau de pozit¸ie<br />
dacă ea restrict¸ionează doar pozit¸iile sistemului ¸si deci aceste legături sunt independente<br />
de vitezele punctelor, adică legăturăa are următoarea formă<br />
ψ (x1, x2, . . . , xN, t) ≥ 0. (1.65)<br />
În general, o legătură este numită neolonomă sau cinematică sau de mi¸scare<br />
dacă relat¸iile care descriu legătura sunt dependente de vitezele punctelor ¸si deci<br />
au forma (1.63).<br />
Exemplele pe care le-am considerat mai sus sunt toate referitoare la legături<br />
olonome. Un exemplu de legătură neolonomă este reprezentat de legăturăa<br />
care determină rostogolirea unei sfere pe un plan fară să alunece. Vom discuta<br />
ulterior modul în care această legătură poate fi definită de o ecuat¸ie de forma<br />
(1.63).
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE33<br />
În domeniul mecanicii, legăturile neolonome nu sunt foarte frecvent întâlnite.<br />
Aceste legături sunt de obicei exprimate prin relat¸ii care sunt liniare în raport cu<br />
vitezele punctelor care formează sistemul. Pentru cazul unui sistem constituit<br />
dintr-un num˘r finit N de puncte supuse unor legături bilaterale, aceste relat¸ii<br />
au următoarea forma:<br />
N<br />
as(x1, . . . , xN , t) · ˙xs + α (x1, . . . , xN , t) = 0. (1.66)<br />
s=1<br />
De fapt, pentru a fi un sistem neolonom, forma diferent¸ială (1.66) trebuie să nu<br />
fie integrabilă. Aceasta înseamna că nu există nicio funct¸ie ˆ F care depinde de<br />
coordonatele punctelor sistemului material, astfel ca<br />
d<br />
dt ˆ F (x1, . . . , xN , t) =<br />
N<br />
as(x1, . . . , xN , t) · ˙xs + α (x1, . . . , xN , t) = 0. (1.67)<br />
s=1<br />
Într-adevăr, dacă o astfel de funct¸ie ar exista, din (1.67), putem obt¸ine<br />
ˆF (x1, . . . , xN , t) = constant, (1.68)<br />
care este o relat¸ie de tipul (1.65), ¸si care caracterizează o legătură olonomă.<br />
Prin urmare, pentru a avea o legătură neolonomă, este esent¸ial ca relat¸ia (1.66)<br />
să nu fie o formă diferentială integrabilă . În caz contrar, legătura olonomă s-ar<br />
putea exprima ca relat¸ii reductibile la ecuat¸ia (1.68).<br />
Definit¸ie 1.2.5 Un sistem material este numit olonom dacă posibilele sale<br />
legături sunt toate olonome ¸si posibilele sale configurat¸ii pot fi identificate în<br />
mod unic de un numă r finit n de parametri independent¸i, q1, q2,. . . , qn, numit¸i<br />
coordonate generalizate sau coordonate lagrangiane. Numărul n este numit<br />
numărul gradelor de libertate al sistemului, sau se spune că sistemul are n<br />
grade de libertate.<br />
Un punct material liber (adică, un punct a cărui mi¸scarea nu este supusă<br />
la nicio legătură, ¸si în consecint¸ă la relat¸ii de tipul (1.63)) reprezintă un sistem<br />
olonom cu trei grade de libertate. Un punct material fort¸at să se mi¸ste pe o<br />
suprafat¸ă, adică este supus unei legături definite de o relat¸ie de următorul tip:<br />
ϕ (x, y, z, t) = 0,<br />
c reprezintă un sistem olonom cu două grade de libertate. Într-adevăr, pentru a<br />
determina pozit¸ia unui punct pe o suprafat¸ă, ne trebuie doi parametri. În final,<br />
un punct constrâns să se mi¸ste pe o curbă reprezintă un sistem olonom cu doar<br />
un grad de libertate, deoarece un parametru este suficient pentru a identifica<br />
pozit¸ia punctului pe o curbă dată.<br />
Este posibil să prezentăm conceptul de grad de libertate ¸si numărul lor<br />
pornind cu un sistem constituit dintr-un număr finit N de puncte materiale
34 CAPITOLUL 1. CINEMATICA<br />
cu constrangeri olonome bilaterale. Dacă sistemul este supus la r < 3N legături<br />
bilaterale, definite de r ecuat¸ii independente de următorul tip:<br />
ψh(x1, x2, . . . , xN, t) = 0, h = 1, 2, . . . , r, (1.69)<br />
atunci există numai n = 3N − r parametri independent¸i, deoarecem folosind<br />
sistemul (1.69), putem, pentru exemplu, exprima r coordonatele în termenii<br />
a celor n = 3N − r rămase. Vorbind mai general, putem găsi n parametri<br />
independent¸i q1, q2,. . . , qn care determină pozit¸ia oricărui punct al sistemul,<br />
adică avem<br />
Ps = Ps(q1, . . . , qn, t), s = 1, 2, . . . , N. (1.70)<br />
Prin urmare, numărul n de grade de libertate a sistemului poate fi obt¸inut<br />
scăzând numărul r al ecuat¸iilor legăturilor din 3N, care este numărul gradelor<br />
de libertate ale unui sistem constituit din puncte libere.<br />
Pentru un sistem material dat, este posibil să asociem n coordonate la-<br />
grangiane q1, q2, . . . , qn într-un număr infinit de moduri.<br />
Într-adevăr, orice<br />
transformare χ : R n → R n care este injectivă ¸si suficient de regulată poate<br />
determina un nou n–uplu de parametri lagrangiani.<br />
Presupunem că, pe lângă legăturile bilaterale, un sistem olonom este de<br />
asemenea supus la legături unilaterale de următoarea formă<br />
ψ ′ (x1, x2, . . . , xN , t) ≥ 0. (1.71)<br />
Este clar că, dacă luăm în calcul doar legături bilaterale, folosind argumentat¸ii<br />
similare cu cele folosite mai sus, putem defini de asemenea, în acest caz, n coordonate<br />
lagrangiane q1, q2, . . . , qn, ¸si în consecint¸ă obt¸inem ecuat¸iile (1.70).<br />
Deoarece coordonatele punctelor sistemului trebuie să satisfacă inegalităt¸ile<br />
(1.71), ace¸sti parametri lagrangiani vor satisface de asemenea inegalităt¸i de<br />
următoarea formă:<br />
ϕ (q1, q2, . . . , qn, t) ≥ 0. (1.72)<br />
Este posibil să explicăm de ce aceste inegalităt¸i nu pot reduce numărul de<br />
grade de libertate ¸si prin ramâne egal cu cea a sistemului care este supus numai<br />
la legături olonome ¸si bilaterale. Spre exemplu, numărul de parametri<br />
independent¸i pentru un punct de constrâns să se mi¸ste într-o camera numărul<br />
gradelor de libertate rămâne egal cu trei, chiar dacă ace¸sti parametri sunt legat¸i<br />
reciproc prin intermediul relatiilor de tipul (1.72), deoarece punctele nu pot<br />
parăsi sala.<br />
Definit¸ie 1.2.6 Un sistem material este numit neolonom dacă este supus la cel<br />
put¸in o legătură neolonomă.<br />
De¸si legăturile neolonome impun unele restrict¸ii cu privire la vitezele punctelor<br />
din sistem, nu le interzice să aibă orice pozitie decât dacă este supusă unei<br />
legături olonome, astfel legăturile neolonome nu reduc numărul de parametri<br />
lagrangiani ai sistemului. Prin urmare, în studiul sistemelor neolonome, trebuie<br />
mai întâi să considerăm sistemul material care este supus doar la legături
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE35<br />
olonome pentru a determina gradele de libertate ¸si a parametrilor lagrangiani;<br />
apoi, trebuie să introducem noi legături neolonome cum ar fi noi ecuat¸ii sau inegalităt¸i.<br />
Dacă q1, q2, . . . , qn sunt parametri lagrangiani a unui sistem cu n grade<br />
de libertate, atunci legătura neolonomă de tipul (1.66) poate fi reprezentată ca<br />
n<br />
Ai(q; t) ˙qi + α(q; t) = 0. (1.73)<br />
i=1<br />
Exercit¸iu 1.2.1 Pentru oricare din următoarele cazuri, determinat¸i dacă legătură<br />
este olonomă sau neolonomă: a) un punct material care se mi¸scă pe un cerc;<br />
b) un punct material greu care se mi¸scă pe un plan înclinat; c) o placă rigidă<br />
alunecând pe un plan fixat x1Ox2; d) un punct material P de coordonate (x1, x2, x3)<br />
este fort¸at să se mi¸ste în a¸sa fel încât componentele vitezei satisfac următoarea<br />
relat¸ie ˙x1 = f(x2, x3) ( ˙x2 + ˙x3), cu ∂f ∂f<br />
= ; e) o lama subt¸ire rigidă fixată<br />
∂x2 ∂x3<br />
pe o placă rigidă care alunecă pe un plan fixat x1Ox2.<br />
Solut¸ie. a) Punctul se mi¸scă pe o curbă ¸si deci este o legătură olonomă.<br />
b) Punctul se mi¸scă pe o suprafat¸ă ¸si deci legătura este olonomă.<br />
c) Placa rigidă se mi¸scă pe un planul înclinat fixat ¸si legătura este olonomă.<br />
d) Dacă această legătură ar fi olonomă, atunci poate fi scrisă în următoarea<br />
formă F (x1, x2, x3) = 0. Din această ipoteză, deducem că<br />
∂F<br />
∂x1<br />
¸si aceasta coincide cu relat¸ia de legătură<br />
dx1 + ∂F<br />
dx2 +<br />
∂x2<br />
∂F<br />
dx3 = 0,<br />
∂x3<br />
dx1 − f(x2, x3)dx2 − f(x2, x3)dx3 = 0,<br />
dacă ¸si numai dacă există λ(x1, x2, x3) astfel ca<br />
∂F<br />
= λ,<br />
∂x1<br />
∂F<br />
= −λf,<br />
∂x2<br />
∂F<br />
= −λf.<br />
∂x3<br />
Dacă vom egala derivatele mixte de ordinul doi ale funct¸iei F , obt¸inem<br />
∂λ<br />
∂x2<br />
¸si deci deducem că<br />
= ∂λ<br />
= −<br />
∂x3<br />
∂λ<br />
f, −<br />
∂x1<br />
∂λ<br />
f − λ<br />
∂x3<br />
∂f<br />
= −<br />
∂x3<br />
∂λ<br />
f − λ<br />
∂x2<br />
∂f<br />
,<br />
∂x2<br />
∂f<br />
∂x2<br />
o relat¸ie care contrazice ipoteza că ∂f<br />
∂x2<br />
= ∂f<br />
,<br />
∂x3<br />
∂f<br />
= . Prin urmare, avem o legătură<br />
∂x3<br />
neolonomă.<br />
e) Deoarece avem o mi¸scare de alunecare, fără rotat¸ie ¸si fără pivotare, rezultă<br />
că viteza unui punct P al lamei subt¸iri este tangentă la lamă. Dacă vom nota<br />
cu x1 ¸si x2 coordonatele punctului ¸si cu θ unghiul format de lamă cu axa Ox1,
36 CAPITOLUL 1. CINEMATICA<br />
atunci x1, x2 ¸si θ constituie coordonatele generalizate pentru lama subt¸ire ¸si, în<br />
plus, avem<br />
dx2<br />
= tan θ. (1.74)<br />
dx1<br />
Legătura de mai sus este neolonomă, deoarece nu este integrabilă. De fapt, dacă<br />
presupunem că există o relat¸ie de tipul F (x1, x2, θ) = 0, deducem că<br />
∂F<br />
∂x1<br />
dx1 + ∂F<br />
dx2 +<br />
∂x2<br />
∂F<br />
dθ = 0, (1.75)<br />
∂θ<br />
¸si deci, luând în considerare legătura (1.74) ¸si deoarece dθ este arbitrar, obt¸inem<br />
∂F<br />
∂θ<br />
= 0,<br />
∂F<br />
∂x1<br />
+ ∂F<br />
tan θ = 0. (1.76)<br />
∂x2<br />
Prima relat¸ie din (1.76) implică că F este independent de θ, în timp ce a doua<br />
relat¸ie din (1.76) conduce la<br />
∂F<br />
= 0,<br />
∂x1<br />
∂F<br />
= 0,<br />
∂x2<br />
deoarece tan θ este arbitrar. Astfel, putem concluziona că F este independent<br />
de x1, x2 ¸si θ ¸si nu poate fi o legătură. Aceasta constradict¸ie provide din faptul<br />
că am presupuns că (1.74) este o legătură olonomă. A¸sadar, legătura (1.74) este<br />
neolonomă.<br />
Exercit¸iu 1.2.2 Determinat¸i numărul gradelor de libertate pentru următoarele<br />
cazuri: a) un punct care se mi¸scă pe o curbă din spat¸iu: b) trei punct care se<br />
mi¸scă liber într-un plan: c) patru puncte care se mi¸scă liber în spat¸iu: d) două<br />
puncte care se mi¸scă în spat¸iu, unite printr-o bară rigidă<br />
Solut¸ie. a) Curba din spat¸iu poate fi dată de reprezentarea naturală x1 =<br />
x1(s), x2 = x2(s), x3 = x3(s). Prin urmare, pozit¸ia punctului pe curbă poate fi<br />
descrisă de parametrul specificats. Astfel, sistemul are un grad de libertate.<br />
b) Fiecare punct cere două coordontate pentru a îi specifica pozit¸ia sa în<br />
plan. Astfel, sunt necesare 3 · 2 = 6 coordonate pentru a specifica pozit¸ia<br />
tuturor celor trei puncte ¸si astfel sistemul are 6 grade de libertate.<br />
c) Pentru a specifica pozit¸ia unui punct material din sistemul considerat,<br />
avem nevoie de trei coordonate. Astfel, sistemul necesită 4 · 3 = 12 coordonate<br />
¸si deci are 12 grade de libertate.<br />
d) Coordonatele (x1, x2, x3) ¸si (y1, y2, y3) ale celor două punct sunt în a¸sa fel<br />
încât distant¸a dintre ele rămâne constantă, adică (x1 − y1) 2 + (x2 − y2) 2 + (x3 −<br />
y3) 2 = constant. Prin urmare, una din cele ¸sase coordonate poate fi exprimată<br />
în termenii celorlalte cinci ¸si deci sistemul are cinci grade de libertate.<br />
Exercit¸iu 1.2.3 Câte grade de libertate are un corp rigid când: a) se mi¸scă<br />
liber în spat¸iul tridimensional; b) are un punct fix ¸si se mi¸scă în jurul lui; c)<br />
are două puncte fixe ¸si se mi¸scă în jurul axei ce trece prin aceste două puncte<br />
distincte; d) se mi¸scă în jurul unei axe fixe; e) se mi¸scă în a¸sa fel încât trei<br />
puncte necoliniare rămân într-un plan fix.
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE37<br />
Solut¸ie. Dacă trei puncte necoliniare Pi, i = 1, 2, 3, ale rigidului sunt fixate,<br />
atunci întreg rigidul este fixat. Astfel, mi¸scarea corpului rigid va fi cunoscută<br />
când cunoa¸stem cum se mi¸scă trei puncte necoliniare ae corpului rigid.<br />
a) Într-un sistem fix de coordonate (O, x1, x2, x3), fie (x1, x2, x3), (y1, y2, y3),<br />
(z1, z2, z3) coordonatele a trei puncte. Deoarece d(P1, P2) = constant, d(P2, P3) =<br />
constant, d(P3, P1) = constant, rezultă că<br />
(x1 − y1) 2 + (x2 − y2) 2 + (x3 − y3) 2 = constant,<br />
(y1 − z1) 2 + (y2 − z2) 2 + (y3 − z3) 2 = constant, (1.77)<br />
(z1 − x1) 2 + (z2 − x2) 2 + (z3 − x3) 2 = constant,<br />
¸si deci putem exprima trei coordonate în termenii a celorlalte ¸sase. Prin urmare,<br />
avem nevoie de ¸sase coordonate independente pentru a descrie mi¸scarea corpului<br />
rigid ¸si deci această mi¸scarea are ¸sase grade libertate.<br />
b) Presupunem punctul fix ca fiind P1 de la punctul anterior a) ¸si deci avem<br />
x1 = a1, x2 = a2, x3 = a3. (1.78)<br />
Astfel, din sistemul de ecuat¸ii descris de relat¸iile (1.77) ¸si (1.78), putem exprima<br />
¸sase coordonate în termenii a altor trei. Prin urmare, mi¸scarea poate fi descrisă<br />
prin trei parametri independent¸i ¸si are trei grade de libertate.<br />
c) Presupunem că punctele fixe sunt P1 ¸si P2 de la punctul a), adică<br />
x1 = a1, x2 = a2, x3 = a3, (1.79)<br />
y1 = b1, y2 = b2, y3 = b3.<br />
Atunci, din sistemul de ecuat¸ii descris de (1.77) ¸si (1.79) putem exprima opt<br />
coordonate în termenii unei singure coordonate. Prin urmare, această mi¸scare<br />
poate fi descrisă doar de un parametru ¸si deci are un grad de libertate. <br />
0 d) Fie u versorul director al unei drepte fixe (d) ¸si fie P0 x1, x0 2, x0 <br />
3 un<br />
punct fixat pe (d). Atunci dreapta fixă (d) are ecuat¸ia vectorială P − P0 = λu,<br />
λ ∈ R. Luăm P1 ¸si P2 pe dreapta (d) ¸si P3 /∈ (d). Atunci, avem P1 − P0 = λ1u,<br />
P2 − P0 = λ2u, cu λ1, λ2 ∈ R ¸si deci<br />
x1 − x 0 1 = λ1u1, x2 − x 0 2 = λ1u2, x3 − x 0 3 = λ1u3,<br />
y1 − x 0 1 = λ2u1, y2 − x 0 2 = λ2u2, y3 − x 0 3 = λ2u3. (1.80)<br />
Prin urmare, din relat¸iile (1.77) ¸si (1.80), puteam exprima nouă coordonate în<br />
termenii parametrilor λ1 ¸si λ2. Bazându-ne pe aceasta, putem concluziona că<br />
această mi¸scare poate fi descrisă prin doi parametri independet¸i ¸si deci are două<br />
grade de libertate.<br />
e) Fie (π) : ax1 + bx2 + cx3 + d = 0 un plan fix. Presupunând că Pi ∈ (π),<br />
i = 1, 2, 3, avem<br />
ax1 + bx2 + cx3 + d = 0,<br />
ay1 + by2 + cy3 + d = 0, (1.81)<br />
az1 + bz2 + cz3 + d = 0.<br />
Din relat¸iile (1.77) ¸si (1.81), putem exprima ¸sase coordonate în termenii altor<br />
trei coordonate ¸si deci putem concluziona că mi¸scarea are trei grade de libertate.
38 CAPITOLUL 1. CINEMATICA<br />
x 1<br />
O<br />
x 3<br />
y 1<br />
O ′<br />
y 2<br />
Figura 1.13:<br />
1.2.2 Cinematica sistemelor rigide<br />
P<br />
Multe sisteme materiale sunt constituite de către un corp rigid sau dintr-un<br />
număr de corpuri rigide conectate între ele. Am observat deja că un corp<br />
rigid este un sistem material supus unor legăruti care conservă distant¸ele între<br />
punctele corpului, adică, dacă P ¸si Q sunt două puncte arbitrare ale corpului,<br />
avem<br />
(P (t) − Q(t)) 2 = constant.<br />
Este important să ret¸inem că un corp rigid este definit ca un model matematic<br />
pentru a descrie multe alte corpuri solide într-un mod suficient de exact.<br />
Astfel de corpuri nu există în natura, deoarece ultimele particule componente<br />
ale oricărui corp (atomi) sunt întotdeauna supuse unor mi¸sci relative. Această<br />
mi¸scare este microscopică ¸si poate fi neglijată atunci când descriem mi¸scarea<br />
macroscopică a corpului. Pe de altă parte, măsurători precise ale acestor corpuri<br />
pot pune în evident¸ă prezent¸a unor mici deformări. Prin urmare, considerăm<br />
corpului a fi rigid numai în cazul în care astfel de deformări nu influent¸ează<br />
mi¸scarea sa.<br />
Considerăm un corp rigid liber, adică, un rigid supus numai la legăturile<br />
de rigiditate. Pentru a studia mi¸scarea sa, vom introduce un sistem ortogonal<br />
de referint¸ă drept (O, x1, x2, x3), pe care îl numim fix în spatiu, fat¸ă de un<br />
observator la care referim mi¸scarea, ¸si un sistem de referint¸ă (ortogonal drept)<br />
(O ′ , y1, y2, y3) fixat în corp (a se vedea Figura 1.13).<br />
Propozit¸ie 1.2.1 Pozit¸ia fiecărui punct al corpului rigid poate fi identificată<br />
dacă se cunoa¸ste configurat¸ia tripletului fixat în corp fat¸ă de cel fixat în spat¸iu.<br />
y 3<br />
x 2
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE39<br />
Demonstrat¸ie. Fie (c1, c2, c3) un sistem de coordonate carteziene cu originea<br />
în punctul O ′ fat¸ă de un sistem de referint¸ă fixat în spact¸iu, fie i1, i2, i3 vectorii<br />
unitari ai axelor x1, x2, x3 ¸si fie j1, j2, j3 vectorii directori ai axelor y1, y2, y3.<br />
Atunci, cosinusurile αhk ale unghiurilor axelor y1, y2, y3 cu axele triedrului fixat<br />
sunt date de matricea<br />
αhk = ih · jk, jk =<br />
3<br />
αhkih. (1.82)<br />
Astfel, este posibil să deducem formula care definet¸e transformarea ce face<br />
legătura dintre coordonatele punctului P calculate în cele două sisteme de<br />
referint¸ă. Astfel, scriem<br />
h=1<br />
P − O = (P − O ′ ) + (O ′ − O). (1.83)<br />
Dacă (x1, x2, x3) and (y1, y2, y3) sunt coordonatele lui P fat¸ă de sistemele<br />
de refet¸ă cu originea în O ¸si respectiv O ′ , atunci putem scrie mai departe (1.83)<br />
astfel<br />
x1i1 + x2i2 + x3i3 = y1j1 + y2j2 + y3j3 + c1i1 + c2i2 + c3i3<br />
sau în forma echivalentă<br />
3<br />
xhih =<br />
h=1<br />
3<br />
ykjk +<br />
k=1<br />
Înmult¸ind ultima ecuat¸ie cu il, obt¸inem<br />
xl = cl +<br />
3<br />
chih.<br />
h=1<br />
(1.84)<br />
3<br />
αlkyk. (1.85)<br />
k=1<br />
Astfel, din (1.85), rezultă că pozit¸ia fiecărui punct P a corpului rigid este determinată<br />
odată cu coordonatele punctului O ′ ¸si matricei (αhk), ale cărui componente<br />
sunt cosinusurile unghiurilor axelor y1, y2, y3,. Folosind notat¸ia<br />
x = x1i1 + x2i2 + x3i3, y = y1j1 + y2j2 + y3j3,<br />
c = c1i1 + c2i2 + c3i3, A = (αhk) ,<br />
este posibil să exprimăm ecuat¸iile (1.84), (1.85) în urmăoarea formă compactă<br />
x = c +<br />
3<br />
ykjk, (1.86)<br />
k=1<br />
x = c + Ay. (1.87)<br />
Rezultă din definit¸ia lui αhk, din a doua ecuat¸ie a relat¸iilor (1.82), că<br />
3<br />
k=1<br />
αikαjk = δij, unde δij =<br />
1 pentru i = j<br />
0 pentru i = j<br />
. (1.88)
40 CAPITOLUL 1. CINEMATICA<br />
În formă matriceală, (1.88) poate fi reprezentată astfel<br />
AA T = 1, unde 1 este matricea unitate.<br />
Prin urmare, matricea A = (αhk) este ortogonală, ¸si astfel reprezintă rotat¸ia,<br />
numită rotat¸ia tripletului (O ′ , y1, y2, y3) fat¸ă de sistemul de referint¸ă centrat în<br />
O ′ ¸si având axele paralele cu axele reperului x1, x2, x3. Deci, pentru a determina<br />
pozit¸ia unui rigid, este suficient să spunem configurat¸ia tripletului fix din el.<br />
Pentru aceasta, este necesar să definim în mod precis cele trei coordonate ale<br />
punctului O ′ ¸si cele nouă cosinusuri αhk, care, totu¸si, sunt legate între ele prin<br />
¸sase relat¸ii (1.88).<br />
Observat¸ie 1.2.1 Un rigid are ¸sase grade de libertate. Totu¸si, pentru a determina<br />
configurat¸ia tripletei solidare cu rigidul, avem nevoie de nouă parametri<br />
independent¸i. Ace¸sti parametri pot fi coordonatele originii O ′ ¸si trei unghiuri<br />
independente, astfel sunt unghiurile lui Euler, care, după cum vom demonstra,<br />
pot fi folosit¸i pentru a defini componentele matricei de rotat¸ie.<br />
În final, mi¸scarea sistemului rigid este cunoscută dacă mi¸scarea punctului<br />
O ′ ¸si legea de schimbare a cosinusurilor αhk sunt determinate; adică<br />
xh(t) = ch(t) +<br />
3<br />
αhk(t)yk, (1.89)<br />
unde yk sunt coordonatele punctului P relativ la sistemul de referint¸ă (O ′ , y1, y2, y3)<br />
¸si care nu depind de timp. Astfel, mi¸scarea unui rigid poate fi considerată ca<br />
suma a două mi¸scări independente, o translat¸ie a unui punct al corpului plus o<br />
rotat¸ie în jurul acestui punct.<br />
Exercit¸iu 1.2.4 O lamă dreptunghiulară ABCD cu dimensiunile AB = 10 ¸si<br />
BC = 20 se mi¸scă astfel încât rămâne mereu paralelă cu un plan fixat x1Ox2.<br />
Mi¸scarea punctului O ′ (c1, c2, c3), de intersect¸ie a diagonalelor lamei, fat¸ă de un<br />
sistem de referint¸ă fix (O, x1, x2, x3) este descrisă de c1(t) = t 2 + 1, c2(t) =<br />
t 2 − 1, c3(t) = 2t. Reperul solidar cu rigidul (O ′ , y1, y2, y3) are o mi¸scare<br />
descrisă de j1 = cos θi1 + sin θi2, j2 = − sin θi1 + cos θi2, j3 = i3, ¸si θ = πt. Să<br />
se determine coordonatele x1, x2, x3 ale vârfurilor lamei la momentul t = 1.<br />
Solut¸ie. Fat¸ă de un sistem de referint¸ă (O ′ , y1, y2, y3), punctele A, B, C, D<br />
pot fi date (de exemplu) de A(−5, −10, 0), B(5, −10, 0), C(5, 10, 0), D(−5, 10, 0).<br />
La momentul t = 1, pozit¸ia sistemului de referint¸ă (O ′ , y1, y2, y3) este dată de<br />
O ′ − O = 2i1 + 2i3 ¸si j1 = −i1, j2 = −i2, j3 = i3, ¸si astfel avem<br />
Deci, deducem că<br />
k=1<br />
x1i1 + x2i2 + x3i3 = 2i1 + 2i3 + y1j1 + y2j2.<br />
x1 = 2 − y1, x2 = −y2, x3 = 2.<br />
Prin urmare, înlocuid y1 = −5, y2 = −10 în relat¸iile de mai sus, deducem<br />
coordonatele punctului A fat¸ă de un sistem de referint¸ă (O, x1, x2, x3) ca fiind<br />
Ax(7, 10, 2). Similar, deducem că Bx(−3, 10, 2), Cx(−3, −10, 2) ¸si Dx(7, −10, 2).
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE41<br />
1.2.3 Mi¸scări particulare ale rigidului<br />
Înainte de a studia mi¸scarea unui corp rigid în forma sa generală, vom considera<br />
câteva forme particulare ale mi¸scării. În plus, atunci când vorbim despre<br />
mi¸scare, ne vom referi întotdeauna la un interval de timp alocat pe care îl vom<br />
nota cu I ⊂ R.<br />
Mi¸scarea de translat¸ie<br />
Definit¸ie 1.2.7 Spunem că un rigid execută o mi¸scare de translat¸ie dacă orice<br />
reper solidar cu rigidul are o mi¸scare de translat¸ie fat¸ă de un reper fix în spat¸iu,<br />
adică matricea cosinusurilor directoare αhk este constantă în tot timpul mi¸scării.<br />
Deoarece<br />
xh(t) = ch(t) +<br />
¸si αhk ¸si yk sunt constante, deducem<br />
3<br />
k=1<br />
αhkyk,<br />
˙xh = ˙ch, ¨xh = ¨ch.<br />
Observat¸ie 1.2.2 În timpul unei mi¸scări de translat¸ie, toate punctele rigidului<br />
au aceea¸si viteză ¸si aceea¸si accelerat¸ie. Viteza comună a tuturor punctelor<br />
corpului poartă numele de câmpul vitezei de translat¸ie. Reciproca este de asemenea<br />
adevărată: dacă la fiecare moment toate punctele din rigid au aceea¸si viteză<br />
atunci rigidul execută o mi¸scare de translat¸ie. Astfel, o mi¸scare de translat¸ie<br />
poate fi definită prin formula<br />
vP (t) = u(t), t ∈ I, (1.90)<br />
unde vP reprezintă viteza unui punct arbitrar P , ¸si u este un vector care nu<br />
depinde de P , pe care îl vom alege să fie egal cu viteza v(O ′ ) a originii O ′ .<br />
Cu ajutorul relat¸iei (1.90), putem deduce următoare formulă pentru deplasarea<br />
relativă elementară a punctului P :<br />
dP = udt = dO ′ .<br />
Exemplu 1.2.1 Un paralelogram articulat ABCD (figure 1.14) este format din<br />
trei bare AB, BC ¸si CD. Punctele A ¸si D sunt fixe iar barele AB ¸si CD<br />
se pot roti în jurul lor, în timp ce barele sunt conenctate în punctele B ¸si C<br />
prin legături articulate. Astfel mi¸scarea barei BC este de translat¸ie, pentru că<br />
sistemul de referint¸ă (B, y1, y2) fixat de bara BC nu-¸si schimbă orientarea în<br />
raport cu sistemul fix de referint¸ă (A, x1, x2).<br />
Definit¸ie 1.2.8 O mi¸scare de translat¸ie se nume¸ste translat¸ie rectilinie (uniformă)<br />
dacă mi¸scarea unui punct arbitrar din rigid este rectilinie (uniformă).
42 CAPITOLUL 1. CINEMATICA<br />
x 2<br />
y 2<br />
A D<br />
B C<br />
Mi¸scarea de rotat¸ie<br />
Figura 1.14:<br />
Cosiderăm un corp rigid care cont¸ine o axă fixă care face obiectul următoarelor<br />
constrângeri: două puncte O ¸si O1 rămân fixe, ¸si deci, în particular, prin proprietăt¸ile<br />
corpurilor rigide, întregul segment cuprins între punctele O ¸si O1<br />
rămân fixe, de asemenea.<br />
Observat¸ie 1.2.3 Un corp rigid care cont¸ine o axă fixă formează un sistem<br />
cu numai un grad de libertate. Este posibil să alegem unghiul ϕ format de un<br />
plan fix al corpului care cont¸ine axa O − O1cu un alt plan fix care cont¸ine de<br />
asemenea axa fixă O − O1 (figura 1.15) ca unic parametru.<br />
Definit¸ie 1.2.9 Mi¸scarea unui corp rigid se nume¸ste de rotat¸ie dacă toate<br />
punctele care se află pe o dreaptă fixă din corp rămân fixe. Această dreaptă<br />
poartă numele de axă de rotat¸ie.<br />
Să alegem, ca un triplet fix în spat¸iu, un sistem de axe ortogonale cu originea<br />
în O ¸si axa x3 având aceea¸si direct¸ie ca (O1 −O). Ca de obicei, se alege tripletul<br />
ata¸sat corpului cu originea în O cu axa y3 să coincidedă cu x3, care deci va<br />
avea aceea¸si direct¸ie cu (O1 − O). Notând cu ϕ unghiul x1y1, putem exprima<br />
cosinusurile directoare αhk ca funct¸ii de acest unghi. Astfel avem<br />
(αhk) =<br />
⎛<br />
⎝<br />
x 1<br />
y 1<br />
cos ϕ − sin ϕ 0<br />
sin ϕ cos ϕ 0<br />
0 0 1<br />
⎞<br />
⎠ . (1.91)<br />
Prin urmare, ecuat¸iile mi¸scării de rotat¸ie ale unui corp rigid au următoarea<br />
formă:<br />
x1(t) = cos ϕ(t) y1 − sin ϕ(t) y2,<br />
x2(t) = sin ϕ(t) y1 + cos ϕ(t) y2, (1.92)<br />
x3(t) = y3.
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE43<br />
x 1<br />
O 1<br />
O<br />
ϕ<br />
x 3 ≡ y 3<br />
y 1<br />
Figura 1.15:<br />
y 2<br />
Folosind sistemul (1.92), putem concluziona că mi¸scare unui punct arbitrar din<br />
rigid este circulară. Ridicând la puterea a doua ¸si sumând primele două ecuat¸ii<br />
ale sistemului (1.92), obt¸inem<br />
x 2<br />
x 2 1(t) + x 2 2(t) = y 2 1 + y 2 2 = constant,<br />
x3 = constant. (1.93)<br />
Sistemul (1.93) denotă ecuat¸ia unui cerc. Prin urmare, viteza unui punct<br />
arbitrar dintr-un corp rigid care execută o mi¸scare de rotat¸ie este dată de formula:<br />
vP = ˙ϕi3 × (P − O), (1.94)<br />
unde O este un punct fix de pe axa de rotat¸ie. Vectorul ω = ˙ϕi3 poartă numele<br />
de viteza unghiulară a corpului rigid.<br />
Folosind (1.94), putem imediat să determinăm următoarea formulă pentru<br />
deplasarea elementara a lui P :<br />
Mi¸scarea de roto–translatie ¸<br />
dP = dϕi3 × (P − O).<br />
Definit¸ie 1.2.10 Mi¸scarea unui corp rigid în care o dreaptă fixă a corpului se<br />
mi¸scă de-alungul unei drepte din spat¸iu se nume¸ste de roto–translat¸ie.<br />
Să alegem din nou tripletul (O ′ , y1, y2, y3) cu axa y3 având aceea¸si direct¸ie<br />
cu drepta fixă a corpului ¸si cu originea O ′ într-un punct de pe această dreaptă,
44 CAPITOLUL 1. CINEMATICA<br />
de coordonate O ′ = (0, 0, c3). Este clar că, în acest caz, αhk sunt date de relat¸ia<br />
(1.91), în timp ce ecuat¸iile de mi¸scare au următoarea formă:<br />
x1(t) = cos ϕ(t) y1 − sin ϕ(t) y2,<br />
x2(t) = sin ϕ(t) y1 + cos ϕ(t) y2, (1.95)<br />
x3(t) = c3(t) + y3.<br />
Astfel, deducem că proiect¸ia mi¸scării punctului P pe planul (x1, x2) este un<br />
cerc. Din (1.95), avem<br />
vP = ˙c3(t)i3 + ˙x1i1 + ˙x2i2. (1.96)<br />
În cele din urmă, luând în considerarea expresia vitezei de rotat¸ie, obt¸inem<br />
vP = ˙c3(t)i3 + ˙ϕi3 × (P − O ′ ) = v(O ′ ) + ω × (P − O ′ ). (1.97)<br />
Din (1.97), deducem formula pentru deplasarea elementară relativă:<br />
dP = dO ′ + dϕi3 × (P − O ′ ). (1.98)<br />
Mi¸scarea de roto–translat¸ie se numne¸ste elicoidală dacă viteza v(O ′ ) din<br />
expresia (1.97) este proport¸ională cu ω.<br />
Exercit¸iu 1.2.5 Un rigid laminat în formă dreptunghiulară ABCD se mi¸scă<br />
în plan în pozit¸ia A ′ B ′ C ′ D ′ , adică vârfurile A, B, C, D se deplasează în<br />
vârfurile A ′ , B ′ , C ′ , D ′ , respectiv. Demonstrat¸i că mi¸scarea se poate scrie ca o<br />
sumă a unor mi¸scări de translat¸ie ¸si de rotat¸ie în jurul unui punct corespunzător<br />
al rigidului.<br />
Solut¸ie. Fie E un punct în dreptunghiul ABCD care corespunde punctului<br />
În primul rând se execută translat¸ia din punctul<br />
E ′ din dreptunghiul A ′ B ′ C ′ D ′ .<br />
E în punctul E ′ , astfel că dreptunghiul ABCD devine A1B1C1D1. Mai departe,<br />
folosind pe E ′ ca punct de rotat¸ie, executăm rotat¸ia de unghi θ a dreptunghiului<br />
A1B1C1D1, unde θ este unghiul dintre dreptele suport ale laturilor AB ¸si<br />
respectiv A ′ B ′ . Astfel mi¸scarea este compusă dintr-o translat¸ie ¸si o rotat¸ie.<br />
1.2.4 Unghiurile lui Euler<br />
Presupunem că, pe lângă sistemul de referint¸ă (O ′ , y1, y2, y3) ata¸sat rigidului,<br />
este dat un nou sistem de referint¸ă (O ′ , z1, z2, z3), cu originea în acela¸si punct O ′ ,<br />
dar care are axele paralele cu axele x1, x2, x3 (Figura 1.16). Configurat¸ia corpu-<br />
lui rigid se va defini din nou în funct¸ie de coordonatele lui O ′ ¸si de cosinusurile<br />
directoare ale axelor y1, y2, y3 în funct¸ie de tripletul z1, z2, z3. În mod normal,<br />
ace¸sti cosinu¸si directori coincid cu αhk. Astfel, matricea de rotat¸ie A = (αhk)<br />
descrie complet orientarea relativă a celor două sisteme. Matricea de rotat¸ie<br />
A cont¸ine trei unghiuri independente. Sunt multe posibilităt¸i de a alege aceste
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE45<br />
x 1<br />
x 3<br />
O<br />
z 1<br />
ψ<br />
n<br />
y 3<br />
ϕ<br />
θ<br />
O ′<br />
y 1<br />
Figura 1.16:<br />
z 3<br />
unghiuri. O posibilitate foarte populară este reprezentată de schema de rotat¸ie<br />
a lui Euler ( 10 ).<br />
Fie drepata n, numită linia nodurilor, obt¸inută ca intersect¸ia dintre planul<br />
(O ′ , y1, y2) cu planul (O ′ , z1, z2). Alegem orientarea acestei drepte astfel încât<br />
(O ′ , z3), (O ′ , y3) ¸si n să formeze un triplet compatibil cu regula mâinii drepte.<br />
Definit¸ie 1.2.11 Numim unghiul de nutat¸ie θ unghiul format de (O ′ , y3) ¸si<br />
(O ′ , z3); de asemenea numim unghiul de precesie ψ unghiul format de (O ′ , z1)<br />
cu n; ¸si în final, numin unghiul de rotat¸ie proprie ϕ unghiul format de n cu<br />
(O ′ , y1). Cele trei unghiuri θ, ψ ¸si ϕ poartă numele de unghiurile lui Euler, ¸si<br />
direct¸iile lor pozitive se găsesc cu regula mâinii drepte aplicată pentru n, (O ′ , z3)<br />
¸si (O ′ , y3).<br />
Observat¸ie 1.2.4 După ce cele două sisteme (O ′ , y1, y2, y3) ¸si (O ′ , z1, z2, z3)<br />
au fost date, unghiurile lui Euler se determină în mod unic. Reciproc, dacă<br />
se dă tripletul (O ′ , z1, z2, z3) ¸si unghiurile lui Euler θ, ψ, ϕ, atunci tripletul<br />
(O ′ , y1, y2, y3) se determină în mod unic. Prin urmare, coeficient¸ii αhk pot fi<br />
determinat¸i din moment ce ψ determină linia nodurilor, θ determină planul care<br />
cont¸ine pe n, ¸si ϕ determină planul definit de (O ′ , y3) ¸si (O ′ , y1).<br />
Unghiurile lui Euler sunt generate prin următoarea serie de rotat¸iiare, care<br />
duc tripletul (O ′ , z1, z2, z3) în tripletul (O ′ , y1, y2, y3) :<br />
I. Prima rotat¸ie este de unghi ψ în sens invers acelor de ceasornic în jurul<br />
axei z3 ¸si transformă sistemul (z1, z2, z3) în sistemul (z ′ 1, z ′ 2, z ′ 3) = (z ′ 1, z ′ 2, z3).<br />
10 Euler (1776).<br />
x 2<br />
P<br />
y 2<br />
z 2
46 CAPITOLUL 1. CINEMATICA<br />
Din moment ce rotat¸ia are loc în planul z1O ′ z2, matricea transformării este<br />
⎛<br />
cos ψ sin ψ 0<br />
⎞<br />
Aψ = ⎝ − sin ψ cos ψ 0 ⎠ (1.99)<br />
0 0 1<br />
¸si<br />
¸si<br />
z ′ = Aψz, (1.100)<br />
i1 = cos ψi ′ 1 − sin ψi ′ 2,<br />
i2 = sin ψi ′ 1 + cos ψi ′ 2, (1.101)<br />
i3 = i ′ 3.<br />
Viteza unghiulară, ωψ care corespunde acestei rotat¸ii infinitezimale în jurul axei<br />
care-l cont¸ine pe i3 este dată de<br />
ωψ = ˙ ψi3. (1.102)<br />
II. A doua rotat¸ie este de unghi θ în sens invers acelor de ceasornic în jurul<br />
axei z ′ 1 ¸si transformă (z ′ 1, z ′ 2, z ′ 3) în (z ′′<br />
1 , z ′′<br />
2 , z ′′<br />
3 ) = (z ′ 1, z ′′<br />
2 , y3). Pentru că rotat¸ia<br />
are loc în planul z ′ 2O ′ z ′ 3, în jurul liniei nodurilor, matricea transformării este<br />
¸si<br />
¸si<br />
⎛<br />
Aθ = ⎝<br />
1 0 0<br />
0 cos θ sin θ<br />
0 − sin θ cos θ<br />
⎞<br />
⎠ (1.103)<br />
z ′′ = Aθz ′ , (1.104)<br />
i ′ 1 = i ′′<br />
1,<br />
i ′ 2 = cos θi ′′<br />
2 − sin θi ′′<br />
3, (1.105)<br />
i ′ 3 = sin θi ′′<br />
2 + cos θi ′′<br />
3.<br />
Dcaă notăm prin n versorul liniei nodurilor, adică n = i ′ 1, atunci vectorul<br />
vitezei unghiulare, ωθ corespunzător acestei rotat¸ii infinitezimale este dat de<br />
ωθ = ˙ θn = ˙ θi ′ 1. (1.106)<br />
III. A treia rotat¸ie este de unghi ϕ în sens invers acelor de ceasornic în jurul<br />
axei z ′′<br />
3 ¸si transformă sistemul (z ′′<br />
1 , z ′′<br />
2 , z ′′<br />
3 ) în sistemul (y1, y2, y3). Pentru că<br />
rotat¸ia are loc în planul z ′′<br />
1 O ′ z ′′<br />
2 , matricea transformării este<br />
⎛<br />
Aϕ = ⎝<br />
cos ϕ sin ϕ 0<br />
− sin ϕ cos ϕ 0<br />
0 0 1<br />
⎞<br />
⎠ (1.107)
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE47<br />
iar<br />
¸si<br />
y = Aϕz ′′ , (1.108)<br />
i ′′<br />
1 = cos ϕj1 − sin ϕj2,<br />
i ′′<br />
2 = sin ϕj1 + cos ϕj2, (1.109)<br />
i ′′<br />
3 = j3.<br />
Pentru că rotat¸ia sistemului are loc în jurul axei j3, rezultă că viteza unghiulară<br />
este dată de<br />
ωϕ = ˙ϕj3. (1.110)<br />
Prin compunerea celor trei rotat¸ii prezentate mai sus, deducem că transformarea<br />
completă din sistemul zi în sistemul yi este dată de<br />
¸si matricea de rotat¸ie A este<br />
y = Aϕz ′′ = AϕAθz ′ = AϕAθAψz, (1.111)<br />
A = AϕAθAψ. (1.112)<br />
T¸ inând cont de relat¸iile (1.99), (1.103), (1.107) ¸si (1.112), deducem că această<br />
matrice are următoarele componente<br />
α11 = cos ϕ cos ψ − cos θ sin ψ sin ϕ,<br />
α21 = − sin ϕ cos ψ − cos θ sin ψ cos ϕ,<br />
α31 = sin θ sin ψ,<br />
α12 = cos ϕ sin ψ + cos θ cos ψ sin ϕ,<br />
α22 = − sin ϕ sin ψ + cos θ cos ψ cos ϕ,<br />
α32 = − sin θ cos ψ,<br />
α13 = sin ϕ sin θ, α23 = cos ϕ sin θ, α33 = cos θ. (1.113)<br />
Pentru ceea ce urmează, este cel mai convenabil să exprimăm cele trei viteze<br />
unghiulare în funct¸ie de sistemul ata¸sat corpului, adică în funct¸ie de versorii<br />
bazei j1, j2, j3. Astfel, prin intermediul relat¸iilor (1.101), (1.102), (1.105),<br />
(1.106), (1.109) ¸si (1.110), avem<br />
ωψ = ˙ ψ (sin θ sin ϕj1 + sin θ cos ϕj2 + cos θj3) , (1.114)<br />
ωθ = ˙ θ (cos ϕj1 − sin ϕj2) , (1.115)<br />
ωϕ = ˙ϕj3. (1.116)<br />
Exercit¸iu 1.2.6 În schema de rotat¸ie a lui Euler, găsit¸i relat¸iile dintre vectorii<br />
unitari i1, i2, i3 ¸si j1, j2, j3.
48 CAPITOLUL 1. CINEMATICA<br />
Solut¸ie. T¸ inând cont de rotat¸iile efectuate prin schema lui Euler, date de<br />
relat¸iile (1.101), (1.105) ¸si (1.109), deducem că<br />
i1 = (cos ϕ cos ψ − cos θ sin ψ sin ϕ) j1 − (sin ϕ cos ψ + cos θ sin ψ cos ϕ) j2 +<br />
+ sin θ sin ψj3,<br />
i2 = (cos ϕ sin ψ + cos θ cos ψ sin ϕ) j1 + (− sin ϕ sin ψ + cos θ cos ψ cos ϕ) j2 −<br />
− sin θ cos ψj3,<br />
Din aceste relat¸ii, vom deduce că<br />
i3 = sin ϕ sin θj1 + cos ϕ sin θj2 + cos θj3. (1.117)<br />
j1 = (cos ϕ cos ψ − cos θ sin ψ sin ϕ) i1 + (cos ϕ sin ψ + cos θ cos ψ sin ϕ) i2 +<br />
+ sin ϕ sin θi3,<br />
j2 = − (sin ϕ cos ψ + cos θ sin ψ cos ϕ) i1 + (− sin ϕ sin ψ + cos θ cos ψ cos ϕ) i2 +<br />
+ cos ϕ sin θi3,<br />
j3 = sin θ sin ψi1 − sin θ cos ψi2 + cos θi3. (1.118)<br />
Exercit¸iu 1.2.7 În schema de rotat¸ie a lui Euler, exprimat¸i vectorii viteză<br />
unghiulară ωψ, ωθ, ωϕ în sistemul fix de coordonate.<br />
Solut¸ie. Considerând defiit¸iile unghiului de presesie ψ ¸si a vectorilui unitate<br />
n, avem (vezi Figura 1.16)<br />
n = cos ψi1 + sin ψi2.<br />
Astfel, relat¸iile (1.102), (1.106), (1.110) ¸si (1.118), ne conduc la<br />
ωψ = ˙ ψi3, ωθ = ˙ θ (cos ψi1 + sin ψi2) ,<br />
ωϕ = ˙ϕ (sin θ sin ψi1 − sin θ cos ψi2 + cos θi3.) (1.119)<br />
Exercit¸iu 1.2.8 Transformarea din sistemul fix de coordonate (O ′ , z1, z2, z3)<br />
în sistemul ata¸sat corpului (O ′ , y1, y2, y3) este descrisă prin următoarea matrice<br />
de rotat¸ie<br />
A = 1<br />
⎛<br />
⎝<br />
8<br />
2 √ 6 − √ 2 2 √ 6 + √ 2 2 √ 3<br />
− √ 6 − 2 √ 2 √ 6 − 2 √ 2 6<br />
2 √ 6 −2 √ 6 4<br />
Utilizând schema de rotat¸ie a lui Euler de mai sus, găsit¸i unghiurile lui Euler<br />
care descriu orientarea relativă a corpului în sistemele de mai sus.<br />
⎞<br />
⎠ .
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE49<br />
Solut¸ie. Transformarea între două baze corespondente este dată de<br />
i1 =<br />
√<br />
2<br />
<br />
2<br />
8<br />
√ √ <br />
3 − 1 j1 − 3 + 2 j2 + 2 √ i2 =<br />
<br />
3j3 ,<br />
√<br />
2<br />
<br />
2<br />
8<br />
√ √ <br />
3 + 1 j1 + 3 − 2 j2 − 2 √ i3 =<br />
<br />
3j3 ,<br />
1<br />
√ <br />
3j1 + 3j2 + 2j3 .<br />
4<br />
Astfel, vectorul director unitar al liniei nodurilor este<br />
Apoi, avem<br />
n = i3 × j3 1<br />
√ <br />
= 3j1 − j2 .<br />
|i3 × j3| 2<br />
cos θ = i3 · j3 = 1<br />
2 , cos ψ = i1 · n =<br />
¸si deci unghiurile lui Euler sunt θ = π<br />
3<br />
1.2.5 Starea de mi¸scare<br />
√<br />
2<br />
2 , cos ϕ = j1<br />
√<br />
3<br />
· n =<br />
2 ,<br />
π π<br />
, ψ = 4 , ϕ = 6 .<br />
După cum am observat deja, mi¸scarea se raportează în permanent¸ă la un anumit<br />
interval de timp. În particular, este de asemenea important să cunoa¸stem<br />
comportarea corpului la un timp t din intervalul I.<br />
Definit¸ie 1.2.12 Numim stare de mi¸scare sau stare cinetică a rigidului la timpul<br />
t, mult¸imea vitezelor tuturor punctelor singulare ale corpului la acel moment.<br />
Definit¸ie 1.2.13 Numim stare de mi¸scare de translat¸ie sau stare cinetică de<br />
translat¸ie la momentul t următoarea distribut¸ie a vitezelor pentru un corp rigid<br />
vP (t) = vO ′(t), (1.120)<br />
adică vitezele tuturor punctelor corpului la momentul t sunt egale cu viteza punctului<br />
particular O ′ .<br />
Observat¸ie 1.2.5 Dacă mi¸scare corpului rigid este una în care starea de mi¸scare<br />
la fiecare moment este de translat¸ie, atunci mi¸scarea este de asemenea de translat¸ie<br />
¸si reciproc.<br />
Definit¸ie 1.2.14 Numim stare de mi¸scare de rotat¸ie sau stare cinetică de rotat¸iela<br />
momentul t următoarea distribut¸ie a vitezelor pentru un corp rigid:<br />
vP (t) = ω(t) × (P − O ′ ), (1.121)<br />
adică distribut¸ia vitezelor la momentul t este aceea¸si ca la mi¸scare de rotat¸ie.<br />
Vectorul ω este numit viteză unghiulară instantanee.
50 CAPITOLUL 1. CINEMATICA<br />
Observat¸ie 1.2.6 Dacă mi¸scarea unui rigid este de rotat¸ie, atunci la fiecare<br />
moment corpul se află într-o stare cinetică de rotat¸ie. Reciproca acestei remarci<br />
nu este adevărată în general.<br />
Într-adevăr, mai târziu, vom arăta că o mi¸scare rigidă plană sau un corp<br />
rigid cu un punct fix la fiecare moment se află într-o stare cinetică de rotat¸ie,<br />
chiar dacă mi¸scarea nu este în general de rotat¸ie.<br />
Definit¸ie 1.2.15 Numim stare de mi¸scare de roto–translat¸ie sau stare cinetică<br />
de roto–translat¸ie sau stare elicoidală la momentul t următoarea distribut¸ie a<br />
vitezelor pentru un corp rigid:<br />
vP (t) = vO ′(t) + ω(t) × (P − O′ ), (1.122)<br />
unde O ′ este un punctcare are viteza vO ′ paralelă cu ω.<br />
Prin urmare, distribut¸ia vitezelor este aceea¸si ca în cazul mi¸scării de roto–<br />
translat¸ie.<br />
Dacă un corp rigid se rote¸ste în jurul unui punct fix O, atunci, la fiecare<br />
moment, este o linie L de puncte din corp sau dintr-o prelungire a sa care sunt<br />
instantaneu în repaus. Linia L poartă numele de axă instantanee de rotat¸ie a<br />
corpului. Se poate arăta că linia L trece prin O ¸si este paralelă cu ω. Rezultă că,<br />
în orice moment, mi¸scarea unui corp rigid cu un punct fix poate fi considerată<br />
o rotat¸ie în jurul unei linii care trece printr-un punct fix.<br />
Dacă rigidul nu are nici un punct fix, atunci în general nu există nici o linie<br />
de puncte care instantaneu să se afle în repaus, dar există o linie L de puncte<br />
care se mi¸scă instantaneu de-alungul liniei L, adică o linie L de-alungul cărei nu<br />
există mi¸scare perpendiculară pe L. Axa care trece prin O ′ ¸si care este paralelă<br />
cu ω poartă numele de axă instantanee de roto–translat¸ie. Punctele acestei linii<br />
au vitezele paralele cu viteza unghiulară ω, presupusă a fi nenulă. Rezultă că,<br />
în orice moment, mi¸scarea unui corp rigid, când nici un punct nu este fix, poate<br />
fi considerată instantaneu o mi¸scare elicoidală: translat¸ia de-alungul unei axe ¸si<br />
rotat¸ia în jurul ei.<br />
Exercit¸iu 1.2.9 Un con se rostogole¸ste fără să alunece pe o suprafat¸ă dură<br />
perfect orizontală. Găsit¸i axa instantanee de roto-translat¸ie a conului.<br />
Solut¸ie. Avem un corp rigid cu un punct fix. Alegem ca originile sistemului<br />
ata¸sat corpului ¸si a sistemului fix să coincidă, aceasta înseamnă O ≡ O ′ . Alegem<br />
sistemul fix astfel încât planul x1Ox2 să coincidă cu planul orizontal ¸si ca Ox3<br />
să corespundă cu direct¸ia verticală. Axa y1 ata¸ată corpului este aleasă astfel<br />
încât să fie axă de simetrie pentru con ¸si planul y2Oy3 să fie ortogonal cu Oy1.<br />
Apoi, folosind schema de rotat¸ie a lui Euler, obt¸inem ϕ = constant ¸si, mai mult,<br />
ψ = θ sin ϕ pentru că rostogolirea este fără alunecare. Astfel, viteza unghiulară<br />
pentru această mi¸scare este dată de ω = ˙ θ (sin ϕi3 + j1). Din moment ce avem<br />
j1 = cos ϕn − sin ϕi3, unde n este versorul liniei nodurilor, rezultă că ω =<br />
˙θ cos ϕn ¸si deci axa instantanee de rotat¸ie are ecuat¸ia P − O = λω, λ ∈ R, adică<br />
P −O = λ ˙ θ cos ϕn. Prin urmare, axa instantanee de rotat¸ie este linia de contact<br />
dintre con ¸si planu orizontal.
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE51<br />
1.2.6 Formula lui Poisson<br />
Să determină starea generală de mi¸scare pentru un corp rigid. În acest scop,<br />
vom diferent¸ia ecuat¸ia (1.83) în raport cu timpul pentru a obt¸ine<br />
d<br />
vP (t) = vO ′(t) +<br />
dt (P − O′ ). (1.123)<br />
Observând că P − O ′ = y1j1 + y2j2 + y3j3, unde yi (i = 1, 2, 3) nu se schimbă în<br />
timp, putem scrie (1.123) în formă următoare:<br />
vP (t) = vO ′(t) +<br />
3<br />
h=1<br />
yh<br />
djh<br />
. (1.124)<br />
dt<br />
Teoremă 1.2.1 (Formula lui Poisson) Fie j1, j2, j3 un triplet ortonormat<br />
de vectori care variază în timp. Atunci, există un vector ω depinzând de timp<br />
astfel încât<br />
djh<br />
dt (t) = ω(t) × jh(t), h = 1, 2, 3. (1.125)<br />
Demonstrat¸ie. Pentru că jh are modulul constant, adică<br />
jh · jh = 1 pentru orice h = 1, 2, 3,<br />
printr-o diferent¸iere directă în raport cu timpul, obt¸inem<br />
djh<br />
dt · jh = 0 pentru orice h = 1, 2, 3.<br />
Pentru că djh/dt ¸si jh sunt ortogonali, rezultă din egalitatea din urmă că există<br />
o familie de vectori ωh (astfel încât componentele lor de-alungul jh pot fi alese<br />
în mod arbitrar) care satisface (vezi A.6, (A.41))<br />
djh<br />
dt = ωh × jh, pentru orice h = 1, 2, 3. (1.126)<br />
Folosindu-ne de posibilitatea de a alege în mod arbitrar pe ωh, vom arăta că<br />
există un ω astfel încât ωh = ω pentru h = 1, 2, 3. În acest scop, pentru că<br />
jh · jk = 0 pentru h = k, deducem că<br />
<br />
<br />
djh<br />
djk<br />
· jk = −jh · . (1.127)<br />
dt<br />
dt<br />
Astfel, folosind (1.126) ¸si pentru h = 1, k = 2, din (1.127) obt¸inem<br />
ω1 × j1 · j2 = −j1 · ω2 × j2.<br />
Folosind proprietăt¸ile produsului mixt, din ultima egalitate rezultă că<br />
ω1 · j3 = ω2 · j3.
52 CAPITOLUL 1. CINEMATICA<br />
Mai mult, pentru că de-alungul lui j1 componenta ω1 este arbitrară, putem alege<br />
ca ea să fie egală cu componenta ω2 de-alungul lui j1, astfel ca (ω1·j1) = (ω2·j1).<br />
În mod analog, pentru că ω2 poate fi ales arbitrar de-alungul lui j2, vom alege<br />
ω2 astfel ca (ω2 · j2) = (ω1 · j2).<br />
În acest mod, ω1 ¸si ω2 au acelea¸si componente de-alungul j1, j2, j3. Să alegem<br />
acum ω3 de-alungul j3 să fie egală cu ω1 = ω2 de-alungul j3. Atunci, pentru<br />
h = 1, k = 3, din (1.127), avem ω1 × j1 · j3 = −j1 · ω3 × j3, ¸si în consecint¸ă,<br />
ω1 · j2 = ω3 · j2. În mod analog, pentru h = 2, k = 3 (1.127), putem arăta că<br />
ω2 · j1 = ω3 · j1. Astfel, putem concluziona că<br />
ω1 = ω2 = ω3,<br />
¸si notăm acest vector comun cu ω. Prin urmare, din (1.126), obt¸inem formula<br />
lui Poisson.<br />
În cele din urmă, din (1.124) ¸si (1.125), obt¸inem expresia generală a câmpului<br />
viteză:<br />
vP (t) = vO ′(t) + ω(t) × (P − O′ ). (1.128)<br />
Observat¸ie 1.2.7 Cea mai generală stare de mi¸scare a corpului rigid poate<br />
fi reprezentată în forma (1.128) ¸si poartă numele de formula fundamentală a<br />
cinematicii sistemelor rigide. Vectorul ω este unic ¸si nu depinde pe punctul P .<br />
De fapt, dacă doi vectori ω ¸si ω ′ există astfel încât (1.128) este satisfăcută,<br />
atunci ar trebui să avem<br />
(ω − ω ′ ) × (P − O ′ ) = 0<br />
pentru toate punctele P ceea ce implică faptul că ω = ω ′ .<br />
Vectorul ω nici nu depinde de punctul O ′ . Într-adevăr, dacă alegem alt<br />
punct O ′′ , rezultă din (1.128) că<br />
Utilizând încă o dată (1.128), avem<br />
vO ′′ = vO ′ + ω(t) × (O′′ − O ′ ).<br />
vP = vO ′′ − ω × (O′′ − O ′ ) + ω × (P − O ′ ) = vO ′′ + ω × (P − O′′ ).<br />
Putem de asemenea să concluzionăm, din (1.128), că fiecare stare de mi¸scare<br />
a rigidului este compusă dintr-o stare de mi¸scare de translat¸ie ¸si dintr-o stare<br />
de mi¸scare de rotat¸ie. Putem vorbi de starea de roto–translat¸ie pentru că viteza<br />
vO ′ nu este în general paralelă cu ω.<br />
Putem să obt¸inem formula (1.128) pentru viteza punctelor unui sistem rigid<br />
folosind relat¸ia (1.87)<br />
x(t) = c(t) + A(t)y. (1.129)<br />
Într-adevăr, prin diferent¸ierea relat¸iei (1.129) în raport cu timpul, obt¸inem<br />
˙x(t) = ˙c(t) + ˙A(t)y. (1.130)
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE53<br />
Pentru că A este o matrice ortogonală, avem că A −1 = A T . Astfel, rezolvând<br />
ecuat¸ia (1.129) pentru y ¸si înlocuind rezultatul în (1.130), avem<br />
˙x(t) = ˙c(t) + ˙A(t)A T (t) [x(t) − c(t)] , (1.131)<br />
unde ˙AA T este o matrice antisimetrică. Pentru a vedea acest lucru, vom<br />
diferent¸ia egalitatea AA T = 1, ¸si astfel avem<br />
Pentru că A T · = ( ˙A T ), din (1.132) deducem<br />
˙AA T + A A T · = 0. (1.132)<br />
˙AA T = −A ˙A T <br />
= − ˙AA T T .<br />
În plus, pentru că matricea ˙AA T ω1i1 + ω2i2 + ω3i3 astfel încât<br />
este antisimetrică, există un vector ω =<br />
˙AA T ⎛<br />
= ⎝ 0 −ω3<br />
ω3 0<br />
ω2<br />
−ω1<br />
⎞<br />
⎠ . (1.133)<br />
−ω2 ω1 0<br />
În plus, considerând x−c = P −O ′ , ˙x = vP , ˙c = vO ′, din (1.131) concluzionăm<br />
vP = vO ′ + ω × (P − O′ ). (1.134)<br />
Observat¸ie 1.2.8 Din (1.128), rezultă că cea mai generală expresie pentru deplasarea<br />
elementară arbitrară a unui rigid este dată de<br />
1.2.7 Teorema lui Mozzi<br />
dP = dO ′ + ωdt × (P − O ′ ).<br />
Din formula fundamentală a sistemelor rigide (1.128), rezultă că starea de<br />
mi¸scare a unui corp rigid poate fi totodeauna scrisă ca o sumă de stări de mi¸scare<br />
de translat¸ie ¸si de rotat¸ie. Cu alte cuvinte, nu rezultă din această formulă că<br />
starea de mi¸scare a unui rigid, cea mai generală este cea de roto–translat¸ie.<br />
Această concluzie se bazează pe următoarea teoremă.<br />
Teoremă 1.2.2 (Mozzi) La fiecare moment, cea mai generală stare de mi¸scare<br />
pentru un sistem rigid, este cea de roto–translat¸ie sau elicoidală. În particular,<br />
poate fi de translat¸ie sau rotat¸ie.<br />
Demonstrat¸ie. Dacă avem ω = 0 în formula (1.128) atunci viteza vO ′<br />
poate fi scrisă întotdeauna ca sumă a două componente: prima este aleasă să<br />
fie paralelă cu ω ¸si este notată prin v <br />
O ′(t), ¸si cea cealaltă este ortogonală cu ω<br />
¸si o notăm cu v⊥ O ′(t). Reamintim că, dând doi vectori ortogonali v⊥ O ′(t) ¸si ω, va
54 CAPITOLUL 1. CINEMATICA<br />
exista întotdeauna un vector notat prin O ′ − O ′′ , care este ortogonal cu v ⊥ O ′(t),<br />
astfel încât<br />
v ⊥ O ′(t) = ω × (O′ − O ′′ ). (1.135)<br />
De aceea, din (1.128), obt¸inem<br />
vP (t) = v <br />
O ′ + v⊥ O ′ + ω × (P − O′ ).<br />
Prin intermediul relat¸iei (1.135), din ultima egalitate rezultă că<br />
vP (t) = v <br />
O ′ + ω × (P − O′′ ). (1.136)<br />
Dacă P = O ′′ , putem conluziona din (1.136) că v ′′ O ′ este paralel cu ω ¸si egal<br />
cu vO ′′. Prin urmare, starea de mi¸scare exprimată de către (1.136) este una de<br />
roto–translat¸ie.<br />
Dacă ω(t) = 0, atunci rezultă că starea de mi¸scare la momentul t este de<br />
translat¸ie. În cele din urmă, dacă v<br />
O ′(t) = 0, atunci va rezulta că atarea de<br />
mi¸scare la momentul t este de rotat¸ie.<br />
Observat¸ie 1.2.9 Dacă, în timpul mi¸scării unui corp rigid, un anumit punct<br />
O ′ rămâne fix, t¸inând cont de faptul că vO ′ = 0, din (1.128) concluzionăm<br />
vP (t) = ω(t) × (P − O ′ ). (1.137)<br />
Astfel, starea relativă de mi¸scare este de rotat¸ie la fiecare moment t. Cu toate<br />
acestea, este clar că mi¸scarea nu este, în general, de rotat¸ie, pentru că vectorul<br />
ω(t) î¸si poate schimba direct¸ia în timp.<br />
Studiul mi¸scării unui rigid se poate exprima întotdeauna fat¸ă de sistemul<br />
(O ′ , z1, z2, z3) cu originea în O ′ ¸si cu axele z1, z2, z3 paralele cu axele din tripletul<br />
x1, x2, x3 fixat în spat¸iu. În raport cu acest sistem, mi¸scarea poate fi văzută<br />
ca o mi¸scare a unui copr rigid cu un punct fix. Din (1.137), pentru starea de<br />
mi¸scare de rotat¸ie, obt¸inem<br />
v ′ P = ω × (P − O ′ ),<br />
unde v ′ P este viteza punctului P în raport cu (O′ , z1, z2, z3). În plus, vectorul ω<br />
este acela¸si ca în (1.128), din moment ce versorii tripletului (O ′ , z1, z2, z3) coincid<br />
cu cei din (O, x1, x2, x3). De aceea, vectorul ω din formula (1.128) define¸ste<br />
starea de mi¸scare de rotat¸ie a corpului în raport cu sistemul de coordonate<br />
(O ′ , z1, z2, z3).<br />
Observat¸ie 1.2.10 Din moment ce avem identitatea<br />
ω × (ω × v0) = (ω · v0) ω − ω 2 v0,<br />
deducem că, pentru viteze unghiulare nenule,<br />
v0 =<br />
ω · v0<br />
ω 2<br />
ω − 1<br />
ω 2 ω × (ω × v0) . (1.138)
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE55<br />
Substituind relat¸ia (1.138) în (1.122), obt¸inem următoarea formă pentru câmpul<br />
viteză:<br />
ω · v0<br />
v =<br />
ω2 <br />
ω + ω× P − O ′ − 1<br />
<br />
(ω × v0) ,<br />
ω2 (1.139)<br />
care este o suma a două componente: una este paralelă cu viteza unghiulară<br />
iar cealaltă se află într-un plan ortogonal pe ω. Componenta ω·v0<br />
ω2 ω poartă<br />
numele de viteza de lunecare, în timp ce cealaltă componentă, perpendiculară<br />
pe ω, reprezintă o viteză de rotat¸ie. Din moment ce viteza punctelor de pe<br />
axa instantanee de roto–translat¸ie L este paralelă cu viteza unghiulară ω(t), din<br />
relat¸ia (1.139) deducem că ecuat¸ia vectorială a dreptei L este<br />
P − O ′ = 1<br />
ω 2 (ω × v0) + λω(t), λ ∈ R. (1.140)<br />
Luând produsul interior a relat¸iei (1.139) cu ω, obt¸inem următoarea proprietate<br />
v(t) · ω(t) = v0(t) · ω(t), (1.141)<br />
adică produsul v(t) · ω(t) este independent de punctul P pe axa instantanee de<br />
roto–translat¸ie.<br />
Axa instantanee de roto–translat¸ie a unui copr rigid de obicei se schimbă<br />
în timp ¸si nu este compusă dintr-o mult¸ime de puncte fixe din rigid. Locurile<br />
geometrice ale axei instantanee de roto–translat¸ie în raport cu sistemele de<br />
coordonate fix ¸si mobil sunt două suprafet¸e plane, ce poartă numele de axoid<br />
fix ¸si mobil, respectiv.<br />
Definit¸ie 1.2.16 Mi¸scarea unui corp rigid se nume¸ste mi¸scare rigidă plană<br />
dacă vitezele puntelor corpului sunt întotdeauna paralele cu un plan fix π.<br />
Teoremă 1.2.3 Pentru cazul mi¸scării rigide plane, starea de mi¸scare este totdeauna<br />
de rotat¸ie ¸si de translat¸ie.<br />
Demonstrat¸ie. Alegem ca sistemul de referint¸ă (O, x1, x2, x3) să fie fix în<br />
spat¸iu, astfel încât axa x3 să fie ortogonală pe π. Sistemul de axe ata¸sat corpului<br />
(O ′ , y1, y2, y3) poate fi ales totdeauna astfel încât axa y3 să aibe accea¸si direct¸ie<br />
precum x3. Prin urmare, versorul j3 este constant. Dacă ω = 0, din formula lui<br />
Poisson avem că<br />
dj3<br />
dt = ω × j3 = 0.<br />
Astfel, ω trebuie să fie paralel cu j3, în timp ce viteza vO ′ din (1.128) este<br />
ortogonală pe j3. Prin urmare, luând în considerare că în aceste ipoteze există<br />
un punct O ′′ astfel încât<br />
din (1.128) obt¸inem<br />
vO ′ = ω × (O′ − O ′′ ),<br />
vP (t) = ω(t) × (O ′ − O ′′ ) + ω(t) × (P − O ′ ) =<br />
= ω(t) × (P − O ′′ ), (1.142)
56 CAPITOLUL 1. CINEMATICA<br />
care este o stare de mi¸scare de rotat¸ie.<br />
Dacă, din contra, ω = 0, atunci, din formula fundamentală a cinematicii<br />
sistemelor rigide (1.128), imediat rezultă că aceasta corespunde unei starări de<br />
mi¸scare de translat¸ie.<br />
Observat¸ie 1.2.11 Amintim că formula fundamentală a cinematicii sistemelor<br />
rigide (1.128) poate fi de asemenea scrisă sub forma (1.124), adică<br />
dj1 dj2 dj3<br />
vP (t) = vO ′(t) + y1 + y2 + y3 . (1.143)<br />
dt dt dt<br />
Dincolo de definit¸ia stării de mi¸scare ¸si, în consecint¸ă, dincolo de distribut¸ia<br />
vitezelor unui corp rigid, este folositor să discutăm despre mult¸imea accelerat¸iilor<br />
punctelor unui corp rigid. Pentru aceasta, vom diferent¸ia ecuat¸ia (1.128) în raport<br />
cu timpul, pentru a obt¸ine<br />
aP (t) = d2O ′ dω<br />
+<br />
dt2 dt × (P − O′ ) + ω × d(P − O′ )<br />
. (1.144)<br />
dt<br />
Folosind încă o dată ecuat¸ia (1.128) scrisă în forma:<br />
obt¸inem<br />
d(P − O ′ )<br />
dt<br />
= ω × (P − O ′ ), (1.145)<br />
aP (t) = d2O ′ dω<br />
+<br />
dt2 dt × (P − O′ ) + ω × [ω × (P − O ′ )] , (1.146)<br />
unde primul termen reprezintă accelerat¸ia punctului O ′ , al doilea apare datorită<br />
variat¸iei vectorului ω, iar al treilea termen poate fi exprimat sub următoarea<br />
formă:<br />
ω × [ω × (P − O ′ )] = ω × [ω × ((P − P ∗ ) + (P ∗ − O ′ ))] =<br />
= ω × [ω × (P − P ∗ )] , (1.147)<br />
unde P ∗ este proiect¸ia lui P pe axa care trece prin O ′ ¸si este paralel ω. Folosind<br />
proprietatea dublului produs vectorial , din (1.147) obt¸inem<br />
ω × [ω × (P − O ′ )] = [ω · (P − P ∗ )]ω − ω 2 (P − P ∗ ) =<br />
= −ω 2 (P − P ∗ ). (1.148)<br />
Din pricina acestui rat¸ionament, ultimul termen în (1.146) reprezintă accelerat¸ia<br />
punctului P într-o mi¸scare uniformă de rotat¸ie în jurul axei care trece prin O ′<br />
¸si este paralelă cu ω. În acest caz, obt¸inem<br />
d2O ′<br />
= 0,<br />
dt2 dω<br />
dt<br />
= 0.
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE57<br />
Observat¸ie 1.2.12 Distribut¸ia accelerat¸iilor pentru un corp rigid poate fi exprimată<br />
nu numai prin ecuat¸ia (1.146), ci ¸si printr-o relat¸ie de tipul (1.143),<br />
adică<br />
aP (t) = aO ′(t) + y1<br />
d2j1 + y2<br />
dt2 d2j2 + y3<br />
dt2 d2j3 . (1.149)<br />
dt2 Definit¸ie 1.2.17 Punctul Q ce apart¸ine rigidului în care accelerat¸ia se anulează<br />
la momentul t poartă numele de pol al accelerat¸iei.<br />
Teoremă 1.2.4 Dacă ω × ˙ω = 0, atunci există un singur pol Q al accelerat¸iei<br />
dat de<br />
Q−O ′ =<br />
1<br />
(ω × ˙ω) 2<br />
2<br />
ω × ˙ω+ω ω · aO ′<br />
<br />
ω + ( ˙ω · aO ′) ˙ω + (ω · aO ′) ˙ω × ω .<br />
(1.150)<br />
Demonstrat¸ie. Să notăm prin Q polul accelerat¸iei, ceea ce înseamnă că<br />
aQ(t) = 0, ¸si să considerăm ρ = Q−O ′ . Atunci, relat¸ia (1.146) oferă următoarea<br />
ecuat¸ie pentru determinarea lui ρ = Q − O ′ :<br />
˙ω × ρ + ω × (ω × ρ) = −aO ′. (1.151)<br />
Demonstrăm faptul că această ecuat¸ie determină un unic ρ.<br />
observăm că, având în ipoteză faptul că ω × ˙ω = 0, avem<br />
În acest scop,<br />
ω × ˙ω · (ω × ˙ω) = (ω × ˙ω) 2 > 0, (1.152)<br />
¸si deci tripletul {ω, ˙ω, ω × ˙ω} constituie o bază în V . Astfel, putem descompune<br />
ρ în baza {ω, ˙ω, ω × ˙ω} după cum urmează:<br />
T¸ inând cont de (1.152), din (1.153), obt¸inem<br />
λ1 =<br />
λ2 =<br />
λ3 =<br />
În plus, prin intermediul relat¸iei (1.151), obt¸inem<br />
ρ = λ1ω + λ2 ˙ω + λ3ω × ˙ω. (1.153)<br />
1<br />
2<br />
(ω × ˙ω)<br />
˙ω × (ω × ˙ω) · ρ,<br />
1<br />
2 (ω × ˙ω) × ω · ρ,<br />
(ω × ˙ω)<br />
1<br />
2 ω × ˙ω · ρ.<br />
(ω × ˙ω)<br />
(1.154)<br />
˙ω × (ω × ˙ω) · ρ = ω × ˙ω · aO ′ + ω2 (ω · aO ′) ,<br />
(ω × ˙ω) × ω · ρ = ˙ω · aO ′,<br />
ω × ˙ω · ρ = − ˙ω · aO ′. (1.155)<br />
Dacă înlocuim relat¸ia (1.155) în (1.154) ¸si rezultatul în (1.153), atunci obt¸inem<br />
relat¸ia (1.150) ¸si astfel demonstrat¸ia este completă.
58 CAPITOLUL 1. CINEMATICA<br />
x 1<br />
O 1<br />
O<br />
x 3<br />
θ<br />
A<br />
A*<br />
G<br />
Figura 1.17:<br />
B*<br />
B<br />
Observat¸ie 1.2.13 Dacă înlocuim (P −O ′ ) = (P −Q)+(Q−O ′ ) = (P −Q)+ρ<br />
în (1.146), obt¸inem<br />
aP (t) = ˙ω × (P − Q) + ω × [ω × (P − Q)] . (1.156)<br />
Astfel, pentru ω × ˙ω = 0, relat¸ia (1.156) dovede¸ste că mi¸scarea unui corp rigid,<br />
la un moment dat, este echivalentă cu o mi¸scare de rotat¸ie instantanee în jurul<br />
polului accelerat¸ie.<br />
1.2.8 Aplicat¸ii<br />
Să considerăm o bară AB de lungime 2l, constrânsă a se mi¸sca astfel încât<br />
capetele sale A ¸si B să rămână în două plane paralele la distant¸a 2h (h < l) unul<br />
de celălalt (Figura 1.17). Alegem sistemul de referint¸ă fix astfel încât originea<br />
O să fie echidistantă fat¸ă de cele două plane, mai mult, axa x3 trebuie să fie<br />
ortogonală pe cele două plane. Din alegerea făcută ¸si din constrângerea impusă,<br />
rezultă că, în timpul mi¸scării , mijlocul G al barei AB rămâne în planul (x1, x2)<br />
tot timpul. De aceea, acest sistem are trei grade de libertate, din moment ce<br />
este posibil să determină pozit¸ia barei AB în funct¸ie de coordonate lui G în<br />
planul (x1, x2) ¸si unghiul θ pe care proiect¸ia lui AB pe planul (x1, x2) îl face,<br />
de exemplu, cu axa x1. Folosindu-ne de formula (1.128) ¸si alegând punctul G<br />
ca origine a sistemului ata¸sat barei, este u¸sor să exprimăm viteza unui punct<br />
arbitrar P al barei cu ajutorul<br />
x 2<br />
vP (t) = vG + ω × (P − G). (1.157)<br />
În plus, deoarece vitezele punctelor barei AB sunt paralele cu cele două plane<br />
fixe, mi¸scarea va fi plană. Mi¸scarea barei AB, relativă la sistemul de referint¸ă<br />
cu originea în G ¸si cu axele paralele cu cele ale sistemului fix, este o rotat¸ie în
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE59<br />
jurul axei care trece prin G, paralelă cu x3, ¸si viteza unghiulară a acestei mi¸scări<br />
este ω = ˙ θi3. Notând coordonatele lui G prin (x1G, x2G), din (1.157) ont¸inem<br />
vP = ˙x1Gi1 + ˙x2Gi2 + ˙ θi3 × [(x1 − x1G)i1 + (x2 − x2G)i2 + x3i3]. (1.158)<br />
Pentru că vG este ortogonal pe ω, întotdeauna va exista un punct C astfel încât<br />
vG = ˙ θi3 × (G − C) ¸si, în consecint¸ă,<br />
vP = ˙ θi3 × (P − C). (1.159)<br />
Aceasta înseamnă că starea de mi¸scare este de rotat¸ie în jurul axei care trece<br />
prin C, paralelă cu i3. Pentru a determinaa coordonatele (x1C, x2C, x3C) ale lui<br />
C, să compară expresiile (1.158) ¸si (1.159). Astfel obt¸inem<br />
[ ˙x1G − ˙ θ(x2 − x2G)]i1 + [ ˙x2G + ˙ θ(x1 − x1G)]i2 = − ˙ θ(x2 − x2C)i1 + ˙ θ(x1 − x1C)i2,<br />
astfel egalitatea de mai sus implică<br />
˙x1G + ˙ θx2G = ˙ θx2C,<br />
˙x2G − ˙ θx1G = − ˙ θx1C.<br />
De aceea, toate punctele C care satisfac ecuat¸ia (1.159) apart¸in dreptei paralele<br />
cu x3 care trece prin punctul din planul (x1, x2) ¸si care are coordonatele<br />
x1C = x1G − ˙x2G<br />
˙θ ,<br />
1.2.9 Cinematica mi¸scărilor relative<br />
x2C = x2G + ˙x1G<br />
. (1.160)<br />
˙θ<br />
Considerând doi observatori distinct¸i reprezentat¸i prin două sisteme Cartesiene<br />
compatibile cu regula mâinii drepte, (O, x1, x2, x3) ¸si (O ′ , y1, y2, y3), care se<br />
mi¸scă unul fat¸ă de altul ¸si sunt înzestrate cu acela¸si sistem de măsurare al<br />
timpului (acela¸si ceas) (Figura 1.18).<br />
O astfel de reprezentare este posibilă pentru că, în cadrul mecanicii clasice,<br />
se presupune că:<br />
1 distant¸a dintre două puncte fixe nu depinde de sistemul de referint¸ă ales.<br />
Acest lucru garantează existent¸a a două triplete ortogonale care se pot<br />
mi¸sca unul fat¸ă celălalt;<br />
2 timpul absolut (adică, timpul este independent de observatori) există. Prin<br />
urmare, dacă t ¸si t ′ sunt momente de timp relative la acela¸si eveniment,<br />
măsurat de două sisteme de referint¸ă diferite, există totdeauna posibilitatea<br />
să setăm ceasurile astfel încât<br />
t = t ′ .
60 CAPITOLUL 1. CINEMATICA<br />
x 1<br />
O<br />
x 3<br />
y 3<br />
O ′<br />
y 1<br />
Figura 1.18:<br />
Pentru u¸surint¸a notării, For the convenience of notation only, vom numi<br />
ulterior sistemele de referint¸ă (O, x1, x2, x3) ¸si (O ′ , y1, y2, y3) fix ¸si respectiv<br />
mobil. Mai mult, considerăm sistemul material constituit dintr-un unic punct<br />
P . Mi¸scarea punctului P fat¸ă de sistemul de referint¸ă fix, pe care-l vom numi<br />
absolut, este determinată de sistemul<br />
sau de ecuat¸ia vectorială echivalentă<br />
x1 = ˆx1(t),<br />
y 2<br />
x 2<br />
x2 = ˆx2(t), (1.161)<br />
x3 = ˆx3(t),<br />
P (t) − O = ˆx1(t)i1 + ˆx2(t)i2 + ˆx3(t)i3, (1.162)<br />
unde i1, i2, i3 sunt versorii axelor x1, x2, x3. În mod analog, mi¸scarea punctului<br />
P în raport cu sistemul de referint¸ă mobil , pe care-l vom numi relativ, este<br />
determinată de<br />
P (t) − O ′ = ˆy1(t)j1 + ˆy2(t)j2 + ˆy3(t)j3, (1.163)<br />
unde ˆy1, ˆy2, ˆy3 ¸si j1, j2, j3 sunt coordonatele punctului P ¸si respectiv versorii<br />
sistemului de coordonate (O ′ , y1, y2, y3).<br />
În cele din urmă, mi¸scarea punctului O ′ , în raport cu sistemul de referint¸ă<br />
(O, x1, x2, x3), este determinată de<br />
O ′ (t) − O = c1(t)i1 + c2(t)i2 + c3(t)i3. (1.164)
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE61<br />
Definit¸ie 1.2.18 Numim viteză absolută va ¸si viteză relativă vr a punctului<br />
P viteza lui P în raport cu sistemul de referint¸ă (O, x1, x2, x3) ¸si respectiv<br />
(O ′ , y1, y2, y3), . De aceea, din (1.162) ¸si (1.163), avem că<br />
va = ˙x1i1 + ˙x2i2 + ˙x3i3, (1.165)<br />
vr = ˙y1j1 + ˙y2j2 + ˙y3j3. (1.166)<br />
Numim viteză de transport vτ viteza punctului în tripletului mobil care coincide<br />
cu P .<br />
Observat¸ie 1.2.14 Pentru că mi¸scarea tripletului mobil este mi¸scarea unui<br />
corp rigid în raport cu (O, x1, x2, x3), viteza de transport a punctului este viteza<br />
pe care P ar putea-o avea în cazul ar fi considerat un punct fix în sistemul de<br />
referint¸ă mobil, ¸si este dată de formula (1.128), sau de relat¸ia echivalentă<br />
dj1 dj2 dj3<br />
vτ = vO ′ + y1 + y2 + y3 . (1.167)<br />
dt dt dt<br />
Teoremă 1.2.5 (Compunerea vitezelor) Viteza absolută a unui punct la<br />
fiecare moment este dată de suma dintre viteza relativă ¸si viteza de transport,<br />
adică<br />
va = vr + vτ . (1.168)<br />
Demonstrat¸ie. Considerăm identitatea<br />
P (t) − O = (P (t) − O ′ (t)) + (O ′ (t) − O).<br />
Printr-o diferent¸iere directă în raport cu t, ¸si t¸inând cont de relat¸iile (1.162),<br />
(1.163) ¸si (1.164), deducem că<br />
dP (t)<br />
dt<br />
= dO′ (t)<br />
dt + ˙y1j1 + ˙y2j2 + ˙y3j3 +<br />
+y1<br />
dj1<br />
dt<br />
dj2 dj3<br />
+ y2 + y3 . (1.169)<br />
dt dt<br />
Prin urmare, utilizând expresiile date de (1.165), (1.166), ¸si (1.167) pentru va,<br />
vr ¸si respectiv vτ , din (1.169) obt¸inem (1.168).<br />
Definit¸ie 1.2.19 Numim accelerat¸ie absolută aa ¸si accelerat¸ie relativă ar a<br />
punctului P accelerat¸iile lui P în sistemele de referint¸ă (O, x1, x2, x3) ¸si respectiv<br />
(O ′ , y1, y2, y3), , adică<br />
aa = ¨x1i1 + ¨x2i2 + ¨x3i3, (1.170)<br />
ar = ¨y1j1 + ¨y2j2 + ¨y3j3. (1.171)<br />
Definit¸ie 1.2.20 Numim accelerat¸ie de transport aτ a lui P accelerat¸ia punctului<br />
în tripletul mobil care coincide cu punctul P . Numim accelerat¸ie complementară<br />
sau accelerat¸ie Coriolis cantitatea definită de formula<br />
unde ω este acela¸si vector ca în formula lui Poisson.<br />
ac = 2ω × vr, (1.172)
62 CAPITOLUL 1. CINEMATICA<br />
Observat¸ie 1.2.15 Deoarece accelerat¸ia de transport a punctului P este aceea<br />
pe care punctul P ar avea-o dacă ar fi considerat fix fat¸ă de sistemul mobil de<br />
referint¸ă, utilizând formula accelerat¸iei punctelor unui rigid, deducem că aτ este<br />
dată de (1.146) sau (1.149), adică<br />
aτ = aO ′ + y1<br />
d2j1 + y2<br />
dt2 d2j2 + y3<br />
dt2 d2j3 . (1.173)<br />
dt2 Teoremă 1.2.6 (Compunerea accelerat¸iilor) Accelerat¸ia absolută a unui punct<br />
la fiecare moment este reprezentată de suma accelerat¸iilor relative, de transport<br />
¸si accelerat¸ia Coriolis, adică<br />
aa(t) = ar(t) + aτ (t) + ac(t). (1.174)<br />
Demonstrat¸ie. Diferent¸iind expresia (1.169) în raport cu timpul, obt¸inem<br />
aa(t) = aO ′ + ¨y1j1<br />
<br />
<br />
dj1 dj2 dj3<br />
+ ¨y2j2 + ¨y3j3 + 2 ˙y1 + ˙y2 + ˙y3<br />
dt dt dt<br />
+y1<br />
d2j1 + y2<br />
dt2 d2j2 + y3<br />
dt2 d2j3 .<br />
dt2 Din ultima ecuat¸ie, folosind expresiile (1.171), (1.173) ¸si formula lui Poison,<br />
obt¸inem<br />
aa(t) = ar + aτ + 2( ˙y1 ω × j1 + ˙y2ω × j2 + ˙y3ω × j3). (1.175)<br />
Prin urmare, folosind expresia (1.172) pentru accelerat¸ia Coriolis, din (1.175)<br />
obt¸inem relat¸ia (1.174) ¸si demonstrat¸a este completă.<br />
Observat¸ie 1.2.16 Într-un sistem neinert¸ial, viteza absolută este dată de expresia<br />
¸si viteza absolută are forma<br />
va = vO ′ + ω × (P − O′ ) + vr<br />
aa = aO ′ + ˙ω × (P − O′ ) + 2ω × vr + ω × [ω × (P − O ′ )] + ar.<br />
Termenul ω×[ω × (P − O ′ )] este cunoscut ca accelerat¸ia centripetă a punctului.<br />
Exercit¸iu 1.2.10 Două puncte materiale P1 ¸si P2 au urnătorii vectori de pozit¸ie:<br />
x1 = 2ti1 − t 2 i2 + (3t 2 − 4t)i3 ¸si respectiv x2 = (5t 2 − 12t + 4)i1 + t 3 i2 − 3ti3, .<br />
Determinat¸i viteza ¸si accelerat¸ia relativă al celui de al doilea punct în raport cu<br />
primul la momentul t = 2.<br />
Solut¸ie. Vitezele celor două puncte materiale sunt<br />
v1 = ˙x1 = 2i1 − 2ti2 + (6t − 4)i3,<br />
v2 = ˙x2 = (10t − 12)i1 + 3t 2 i2 − 3i3,
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE63<br />
¸si astfel, la momentul t = 2, avem v1 = 2i1 − 4i2 + 8i3 ¸si v2 = 8i1 + 12i2 − 3i3.<br />
Astfel, viteza relativă a punctului P2 în raport cu P1 este v2 − v1 = 6i1 + 16i2 −<br />
11i3.<br />
Accelerat¸iile celor două puncte materiale sunt<br />
a1 = ¨x1 = −2i2 + 6i3, a2 = ¨x2 = 10i1 + 6ti2,<br />
astfel că, la momentul t = 2, avem a1 = −2i2 + 6i3, a2 = 10i1 + 12i2. Deci,<br />
accelerat¸ia relativă cerută este a2 − a1 = 10i1 + 14i2 − 6i3.<br />
Exercit¸iu 1.2.11 Un unghi invariant y1Oy2 se rote¸ste în planul său în jurul<br />
punctului O cu viteza unghiulară ω. Un punct material P se mi¸scă în planul<br />
acestui unghi ¸si are coordonatele y1 ¸si y2 în raport cu sistemul de referint¸ă plan<br />
(O, y1, y2). 10 ) S˘ se determine componentele vitezei absolute a punctului P<br />
în raport cu axele Oy1 ¸si Oy2; 2 0 ) Dacă y1Oy2 = π<br />
2<br />
absolute sunt C1<br />
y1<br />
¸si componentele vitezei<br />
¸si C2<br />
y2 , unde C1 ¸si C2 sunt constante arbitrare, să se arate că<br />
y2 1 + y2 2 este o funct¸ie liniară în raport cu timpul; 30 ) Dacă y1Oy2 = π<br />
2 ¸si ω este<br />
o constantă, să se determine traiectoria punctului P atunci când accelerat¸iile<br />
absolută ¸si relativă sunt egale.<br />
Solut¸ie. 1 0 ) Fie j1 ¸si j2 versorii axelor Oy1 ¸si respectiv Oy2. Dacă introducem<br />
notat¸ia α = y1Oy2 ¸si considerăm vectorul unitar u ortogonal versorului<br />
j1 atunci avem<br />
j2 = cos αj1 + sin αu<br />
¸si astfel deducem<br />
d<br />
dt j1 = ωu,<br />
d<br />
u = −ωj1,<br />
dt<br />
Deoarece P − O = y1j1 + y2j2, rezultă că<br />
¸si deci<br />
va =<br />
va = ˙y1j1 + ˙y2j2 + y1<br />
d<br />
dt j2 = ω(cos αu − sin αj1).<br />
d<br />
dt j1<br />
d<br />
+ y2<br />
dt j2<br />
<br />
˙y1 − ωy1 cot α − ωy2<br />
<br />
j1 + ˙y2 + ωy2 cot α +<br />
sin α<br />
ωy1<br />
<br />
j2.<br />
sin α<br />
2 0 ) Utilizând ipotezele problemei ¸si rezultatele de mai sus, deducem ecuat¸iile<br />
care implică<br />
Prin integrare, obt¸inem<br />
˙y1 − ωy2 = C1<br />
, ˙y2 + ωy1 = C2<br />
,<br />
y1<br />
y1 ˙y1 + y2 ˙y2 = C1 + C2.<br />
y 2 1 + y 2 2 = 2(C1 + C2)t + C, C = constant.<br />
y2
64 CAPITOLUL 1. CINEMATICA<br />
3 0 ) Pe baza ipootezelor, găsim<br />
ar = ¨y1j1 + ¨y2j2, at = −ω 2 (y1j1 + y2j2), ac = 2ω(− ˙y2j1 + ˙y1j2).<br />
Dacă aa = ar, atunci at + ac = 0 ¸si deci putem deduce<br />
2 ˙y1 − ωy2 = 0, 2 ˙y2 + ωy1 = 0,<br />
din care rezultă că punctul se mi¸scă pe cercul y 2 1 +y 2 2 = c 2 , c = constant. Astfel,<br />
traiectoria este cercul de ecuat¸ie y 2 1 + y 2 2 = c 2 .<br />
1.2.10 Mi¸scări de transport speciale<br />
Mi¸scarea unui sistem mobil de coordonate fat¸ă de sistemul fix poartă numele<br />
de mi¸scare de transport. În această sect¸iune, vom considera câteva mi¸scări de<br />
transport speciale. Începem cu cea de translat¸ie, adică presupunem că sistemul<br />
mobil execută o mi¸scare de translat¸ie fat¸ă de sistemul fix. Considerând aceste<br />
ipoteze, toate punctele ata¸sate sistemului mobil au aceea¸si viteză ¸si accelerat¸ie.<br />
Prin urmare, putem alege ca viteză ¸si accelerat¸ie de transport a oricărui punct<br />
P , viteza ¸si accelerat¸ia originii O ′ a sistemului mobil. Teorema compunerii<br />
vitezelor are următoarea formă:<br />
va(t) = vr(t) + vO ′(t). (1.176)<br />
În ceea ce prive¸ste teorema de compunere a accelerat¸iilor, observă că ac =<br />
2ω × vr = 0, pentru că mi¸scarea de transport este de translat¸ie ¸si astfel ω = 0.<br />
Ecuat¸ia (1.174) devine<br />
aa(t) = ar(t) + aO ′(t). (1.177)<br />
În particular, dacă mi¸scare de transport este de translat¸ie, rectilinie ¸si uniformă,<br />
atunci relat¸ia (1.177) se reduce la<br />
aa(t) = ar(t), (1.178)<br />
adică accelerat¸ia lui P este aceea¸si în raport cu cei doi observatori.<br />
Definit¸ie 1.2.21 Două sisteme de referint¸ă se numesc echivalente dacă mi¸scare<br />
fiecărui sistem în raport cu celălalt este de translat¸ie, rectilinie ¸si uniformă.<br />
Observat¸ie 1.2.17 Cu toate că accelerat¸ia are un caracter relativ, este invariantă<br />
fat¸ă de clasa sistemelor de referint¸ă echivalente. Nu acela¸si lucru se<br />
întâmplă cu viteza, care întotdeauna are un caracter relativ.<br />
În cele din urmă, considerăm mi¸scarea uniformă de transport de rotat¸ie în<br />
jurul unei axe care trece prin O. Astfel, avem<br />
vτ = ω × (P − O),<br />
aτ = −ω 2 (P − P ∗ ),<br />
unde P ∗ este proiect¸ia lui P pe axa de rotat¸ie. Prin urmare, găsim<br />
va = vr + ω × (P − O),<br />
aa = ar − ω 2 (P − P ∗ ) + 2ω × vr.
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE65<br />
x 1<br />
y 3<br />
O<br />
x 3<br />
z 3<br />
O ′′<br />
z 1<br />
O ′<br />
Figura 1.19:<br />
1.2.11 Mi¸scări relative pentru corpurile rigide<br />
y 1<br />
P<br />
Considerăm un corp rigid B care se mi¸scă în raport cu două sisteme de coordonate<br />
(O, x1, x2, x3) ¸si (O ′ , y1, y2, y3). Numim mi¸scare absolută a rigidului B<br />
mi¸scarea fat¸ă de sistemul fix de referint¸ă , ¸si mi¸scare relativă, mi¸scarea în raport<br />
cu sistemul mobil de referint¸ă. Mai mult, putem asocia corpului rigid B un<br />
nou sistem fix de referint¸ă ortogonal în B pe care-l vom nota cu (O ′′ , z1, z2, z3)<br />
(Figure 1.19).<br />
Definit¸ie 1.2.22 Numim viteză unghiulară absolută ωa ¸si viteză unghiulară<br />
relativă ωr vitezele unghiularea ale lui B în timpul mi¸scării sale în raport cu<br />
(O, x1, x2, x3) ¸si respectiv (O ′ , y1, y2, y3).<br />
Observat¸ie 1.2.18 Pentru fiecare punct P ∈ B, avem<br />
y 2<br />
z 2<br />
v (a)<br />
P (t) = v(a)<br />
O ′′(t) + ωa(t) × (P − O ′′ ), (1.179)<br />
v (r)<br />
P (t) = v(r)<br />
O ′′(t) + ωr(t) × (P − O ′′ ), (1.180)<br />
unde v (a)<br />
P ¸si v (a)<br />
O ′′, v(r)<br />
P ¸si v (r)<br />
O ′′ sunt vitezele lui P ¸si O′′ fat¸ă de sistemul de<br />
referint¸ă fix în spat¸iu (O, x1, x2, x3) ¸si respectiv fat¸ă de sistemul mobil de referint¸ă<br />
(O ′ , y1, y2, y3).<br />
x 2
66 CAPITOLUL 1. CINEMATICA<br />
y 2<br />
O<br />
x 2<br />
ωt<br />
x 0<br />
Figura 1.20:<br />
P<br />
r<br />
Definit¸ie 1.2.23 Numim viteză unghiulară de transport ωτ viteza unghiulară<br />
a lui B considerat a fi ata¸sat sistemului mobil de referint¸ă (O ′ , y1, y2, y3).<br />
Observat¸ie 1.2.19 Dacă P este un punct al lui B, atunci<br />
unde v (τ)<br />
P<br />
y 1<br />
x 1<br />
v (τ)<br />
P (t) = v(τ)<br />
O ′′(t) + ωτ (t) × (P − O ′′ ), (1.181)<br />
¸si v(τ)<br />
O ′′ sunt vitezele de transport ale lui P ¸si respectiv O′′ .<br />
Teoremă 1.2.7 Viteza unghiulară absolută ωa a sistemului rigid la fiecare moment<br />
este reprezentată de suma vitezei unghiulară relativă ωr cu viteza de transport<br />
ωτ ; ceea ce înseamnă<br />
ωa = ωr + ωτ . (1.182)<br />
Demonstrat¸ie. Pe baza definit¸iilor lui v (a)<br />
P ¸si v(a)<br />
O ′′, v(r)<br />
P<br />
utilizând Teorema Compunerii Vitezelor, găsim<br />
v (a)<br />
P = v (r)<br />
P + v(τ)<br />
P ,<br />
v (a)<br />
O ′′ = v(r)<br />
O ′′ + v(τ)<br />
O ′′.<br />
¸si v(r)<br />
O ′′, v(τ)<br />
P ¸si v(τ)<br />
O ′′,<br />
Astfel, scăzând (1.180) ¸si (1.181) din (1.179), pentru fiecare P , deducem că<br />
(ωa − ωr − ωτ ) × (P − O ′′ ) = 0. (1.183)<br />
Prin urmare, pentru că îl putem alege pe P în mod arbitrar, expresia (1.182)<br />
din Teorema Compunerii Vitezelor Unghiulare rezultă din (1.183).<br />
1.2.12 Aplicat¸ii ¸si exemple<br />
Un punct P se mi¸scă uniform de-alungul dreptei r. În raport cu (O, x1, x2, x3),<br />
linia r trece prin originea O ¸si se rote¸ste cu viteza unghiulară ω în jurul axei x3<br />
(Figura 1.20).<br />
Prin urmare, cunoa¸stem mi¸scare lui P în raport cu dreapta r ¸si mi¸scare<br />
dreptei r în raport cu sistemul de referint¸ă (O, x1, x2, x3) pe care-l considerăm
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE67<br />
fix în spat¸iu. Utilizând teorema copunerii vitezelor ¸si accelerat¸iilor, putem determina<br />
mi¸scare lui P în raport cu sistemul fix de referint¸ă. De fapt, dacă v<br />
este modului vitezei lui P de-alungul dreptei r, atunci avem<br />
vr = vj1, vτ = ωi3 × (P − O), (1.184)<br />
unde j1 este versorul dreptei r care are aceea¸si direct¸ie cu vr. Prin urmare, dacă<br />
momentul init¸ial al măsurării timpului este ales acela când unghiul P Ox1 = ωt,<br />
atunci avem<br />
j1 = cos ωti1 + sin ωti2,<br />
y1 = vt + x0,<br />
unde x0 semnifică pozit¸ia init¸ială a lui P pe r, când r coincide cu axa x1. Astfel,<br />
prin Teorema Compunerii Vitezelor, putem concluziona că<br />
va = vj1 + ωi3 × (vt + x0)j1<br />
= v(cos ωti1 + sin ωti2) + ω(vt + x0) cos ωti2 − ω(vt + x0) sin ωti1.<br />
Din ultima egalitate deducem<br />
˙x1 = v cos ωt − ω(vt + x0) sin ωt,<br />
˙x2 = v sin ωt + ω(vt + x0) cos ωt. (1.185)<br />
Folosind condit¸iile init¸iale din nou, este u¸sor să verificăm că (1.185) implică<br />
¸si, în consecint¸ă,<br />
x1 = (vt + x0) cos ωt,<br />
x2 = (vt + x0) sin ωt, (1.186)<br />
x1<br />
x2<br />
= cot ωt.<br />
În cele din urmă, dacă introducem sistemul polar de referint¸ă (ρ, θ) astfel<br />
, atunci, din (1.186), obt¸inem<br />
încât ρ = y1 ¸si θ = ωt = arccot x1<br />
x2<br />
ρ =<br />
<br />
x 2 1 + x2 2 = vt + x0 = v<br />
ω θ + x0. (1.187)<br />
Din (1.187) rezultă că traiectoria este spirala lui Arhimede.<br />
În ce prive¸ste accelerat¸ia absolută, observăm că<br />
ar = 0, a τ = −ω 2 (P − O), ac = 2ωi3 × vj1.<br />
Astfel, în funct¸ie de componentele radială ¸si transversală, găsim<br />
aa = −ω 2 y1r + 2ωvh,<br />
unde am introdus notat¸iile r = j1, h = i3 × j1.
68 CAPITOLUL 1. CINEMATICA<br />
x ′<br />
1<br />
n<br />
y 3<br />
ψ<br />
y 1<br />
O<br />
ϕ<br />
θ<br />
′<br />
′<br />
x3<br />
Figura 1.21:<br />
x ′<br />
2<br />
Ca o aplicat¸ie la Teorema Compunerii Vitezelor Unghiulare, vom deduce<br />
relat¸ia dintre viteza unghiulară ω a rigidului ¸si unghiurile lui Euler. Considerăm<br />
un corp rigid liber ¸si expresia corespunzătoare a vitezei punctului P<br />
y 2<br />
n ′<br />
′′ n<br />
vP (t) = vO ′(t) + ω × (P − O′ ),<br />
unde O ′ este un punct al rigidului pe care îl alegem drept origine pentru sistemele<br />
de referint¸ă (O ′ , y1, y2, y3) ata¸sat corpului ¸si a unui nou sistem de referint¸ă<br />
(O ′ , x ′ 1, x ′ 2, x ′ 3) ale cărui axe sunt paralele sau invariabile fat¸ă de (O, x1, x2, x3)<br />
fixat în spat¸iu. Să notăm prin θ, ϕ, ψ unghiurile lui Euler (ca în Figura 1.21).<br />
Vom nota de asemenea prin n ′ o nouă axă, astfel încât (O ′ , n, n ′ , y3) este un<br />
triplet ortogonal care respectă regula mâinii drepte ¸si prin n ′′ , o altă axă, astfel<br />
încât (O ′ , n, n ′′ , x ′ 3) este de asemenea un triplet ortogonal care respectă regula<br />
mâinii drepte. Dacă vom considera mi¸scarea corpului relativ la cele două triplete<br />
(O ′ , x ′ 1, x ′ 2, x ′ 3) ¸si (O ′ , n, n ′ , y3), din Teorema Compunerii Vitezelor Unghiulare<br />
deducem<br />
ω = ωr + ωτ ,<br />
unde ω reprezintă vectorul rotat¸ie instantanee a mi¸scării corpului în raport<br />
cu (O ′ , x ′ 1, x ′ 2, x ′ 3). Mai mult, avem ωr = ˙ϕj3, pentru că, fat¸ă de tripletul<br />
(O ′ , n, n ′ , y3), corpul se mi¸scă în jurul axei y3. În cele din urmă, ωτ este viteza<br />
unghiulară a mi¸scării tripletului (O ′ , n, n ′ , y3) în raport cu (O ′ , x ′ 1, x ′ 2, x ′ 3). De<br />
aceea, considerând mi¸scarea lui (O ′ , n, n ′ , y3) ca fiind una rigidă, o putem studia<br />
în raport cu cei doi observatori asociat¸i celor două sisteme de referint¸ă<br />
(O ′ , n, n ′′ , x ′ 3) ¸si (O ′ , x ′ 1, x ′ 2, x ′ 3). Astfel, obt¸inem<br />
ωτ = ω ′ r + ω ′ τ ,<br />
unde ω ′ r = ˙ θn (n este versorul liniei nodurilor) este viteza unghiulară a mi¸scării<br />
lui (O ′ , n, n ′ , y3) relativ la (O ′ , n, n ′′ , x ′ 3), ¸si ω ′ τ = ˙ ψi3 este viteza unghiulară a
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE69<br />
lui (O ′ , n, n ′′ , x ′ 3) relativ la (O ′ , x ′ 1, x ′ 2, x ′ 3). Prin urmare,<br />
ω = ˙ϕj3 + ˙ θn+ ˙ ψi3. (1.188)<br />
Dacă ω = ω1j1 + ω2j2 + ω3j3 ¸si folosim relat¸iile (1.105), (1.109) ¸si (1.117),<br />
atunci din relat¸ia (1.188) găsim<br />
ω1 = ˙ θ cos ϕ + ˙ ψ sin θ sin ϕ,<br />
ω2 = − ˙ θ sin ϕ + ˙ ψ sin θ cos ϕ, (1.189)<br />
ω3 = ˙ϕ + ˙ ψ cos θ.<br />
Exercit¸iu 1.2.12 În schema de rotat¸ie a lui Euler, găsit¸i relat¸iile care exprimă<br />
componentele vitezei unghiulare ω a rigidului, în funct¸ie de tripletul i1, i2, i3,<br />
în termenii unghiurilor lui Euler.<br />
Solut¸ie. Dacă ω = ω 0 1i1 + ω 0 2i2 + ω 0 3i3, atunci, prin intermediul relat¸iilor<br />
(1.118) ¸si (1.188) ¸si considerând n = cos ψi1 + sin ψi2, obt¸inem următoarele<br />
rela¸tii:<br />
ω 0 1 = ˙ θ cos ψ + ˙ϕ sin θ sin ψ,<br />
ω 0 2 = ˙ θ sin ψ − ˙ϕ sin θ cos ψ,<br />
ω 0 3 = ˙ ψ + ˙ϕ cos θ.<br />
1.2.13 Mi¸scări rigide plane<br />
Ne reamintim că am definit mi¸scarea rigidă plană ca fiind mi¸scarea rigidă în<br />
care vitezele punctelor trebuie să rămână permanent paralele cu un plan fix π.<br />
Mai mult, pentru că axa y3 a sistemului de referint¸ă ata¸sat corpului poate fi<br />
aleasă astfel încât să fie permanent paralelă cu axa x3 a sistemului de referint¸ă<br />
fix, este de asemenea posibil să alegem originea O ′ astfel ca, în timpul mi¸scării,<br />
planele (x1, x2) ¸si (y1, y2) să coincidă întotdeauna cu planul π (Figura 1.22).<br />
Deoarece configurat¸ia corpului este determinată atunci când pozit¸ia tripletului<br />
ata¸sat corpului este determinată, putem spune:<br />
Observat¸ie 1.2.20<br />
În timpul unei mi¸scări rigide plane, corpul posedă trei grade<br />
de libertate. Doi parametri care descriu această mi¸scarea sunt coordonatele lui<br />
O ′ în planul (x1, x2) iar al treilea este unghiul format de axa y1 cu axa x1. Ace¸sti<br />
parametri determină în mod unic configurat¸ia corpului. Astfel, mi¸scare corpului<br />
este determinată de mi¸scare sistemului (O ′ , y1, y2), sau, dintr-un punct de<br />
vedere mai general, de mi¸scarea figurii plane rigide dată de intersect¸ia corpului<br />
cu panul (O, x1, x2).<br />
Mai mult, după cum am demonstrat în Theorem 1.2.3, la fiecare moment,<br />
starea de mi¸scare este de rotat¸ie sau de translat¸ie. Deoarece cazul de mai sus<br />
poate fi considerat un caz degenerat al cazului stării de mi¸scare de rotat¸ie cu axa<br />
instantanee de rotat¸ie la infinit, vom presupune în consecint¸ă că starea cinetică<br />
este întotdeauna de rotat¸ie.
70 CAPITOLUL 1. CINEMATICA<br />
x 1<br />
π<br />
x 3<br />
O<br />
θ<br />
O ′<br />
y 3<br />
Figura 1.22:<br />
y 1<br />
Definit¸ie 1.2.24 Numim centru instantaneu al rotat¸ie C punctul de intersect¸ie<br />
al axei instantanee de rotat¸ie cu planul π.<br />
Observat¸ie 1.2.21 Viteza fiecărui punct P al planului (O ′ , y1, y2) sau cea a<br />
figurii plane rigide ce se află în acest plan este dată de<br />
y2<br />
x 2<br />
vP (t) = ω(t) × (P (t) − C(t)), (1.190)<br />
unde C(t) este centru instrantaneu al rotat¸ie la momentul t.<br />
După ce facem produsul interior al expresiei (1.190) cu (P − C), obt¸inem<br />
vP · (P − C) = ω × (P − C) · (P − C) = 0, (1.191)<br />
¸si astfel, din (1.191) rezultă că viteza fiecărui punct P al figurii plane este ortogonală<br />
pe segmentul cu capetele în P ¸si C. De aceea, putem formula următoarea<br />
propozit¸ie:<br />
Propozit¸ie 1.2.2 Viteza fiecărui punct al figurii plane este determinată de<br />
îndată ce am determinat pozit¸ia centrului instantaneu ¸si viteza unui punct al<br />
figurii plane.<br />
Vom considera două puncte P0 ¸si P ale figurii plane. Presupunem că viteza<br />
vP0 a punctului P0 ¸si pozit¸ia centrului instantaneu C sunt date (Figura 1.23).<br />
Atunci, direct¸ia vitezei unui punct generic P este determinată de vectorul (P −<br />
C), pe care trebuie să fie ortogonală ¸si sensul ei este acela¸si cu cel al lui vP0 . În<br />
cele din urmă, putem afla modulul lui vP by (1.190), t¸inând cont de faptul că<br />
vP = ωr, vP0 = ωr0, (1.192)
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE71<br />
x 2<br />
O<br />
v P0<br />
P 0<br />
P<br />
v P<br />
r 0<br />
Figura 1.23:<br />
r<br />
unde r ¸si r0 sunt distant¸ele dintre P ¸si C, ¸si respectiv dintre P0 ¸si C. Astfel,<br />
folosind (1.192), obt¸inem valoarea modulului ca fiind vP = vP0 r/r0.<br />
Este folositor să observăm că, din (1.190), punctul figurii plane care coincide<br />
cu centrul instantaneu are viteza zero. Mai mult, dacă centrul C se mi¸scă către<br />
infinit, (P0 − C) ¸si (P − C) tind să devină paralele cu vP0 ¸si vP , ¸si rat¸ia<br />
|P0 − C|<br />
|P − C|<br />
C<br />
x 1<br />
r0<br />
=<br />
r ≤ |P0 − P | + |P − C|<br />
= 1 +<br />
|P − C|<br />
|P0 − P |<br />
|P − C|<br />
tinde spre 1, deoarece, în timp ce centrul C se îndreaptă spre infinit, |P0 − P | / |P − C|<br />
tinde către zero. De aceea, pentru că direct¸iile lui vP0 ¸si vP coincid, putem concluziona<br />
că toate punctele au aceea¸si viteză ¸si deci acest lucru este în concordant¸ă<br />
cu observat¸ia anterioară care spunea că starea de mi¸scare este de translat¸ie când<br />
C este la infinit.<br />
Deoarece viteza fiecărui punct al figurii plane este ortogonală cu raza care<br />
une¸ste punctul cu centrul instantaneu, este posibil să determinăm pozit¸ia centrului<br />
C folosind următoarea metodă(vezi Figura 1.23):<br />
Observat¸ie 1.2.22 Dacă, la un moment dat, direct¸iile traiectoriilor sau vitezele<br />
a două puncte ale figurii plane sunt cunoscute, atunci centrul instantaneu al<br />
mi¸scării relative este dat de intersect¸ia dreptelor care sunt ortogonale pe direct¸iile<br />
traiectoriilor sau pe vitezele celor două puncte considerate.<br />
Exemplu 1.2.2 Considerăm o bară AB ale cărei capete sunt constrânse să<br />
se mi¸ste de-alungul a două axe mutual ortogonale x1 ¸si x2. Atunci, centrul<br />
instantaneu C al acestei mi¸scări rigide plane la fiecare moment este punctul de<br />
intersect¸ie al dreptelor normale traiectoriilor (adică al dreptelor ortogonale pe<br />
axele x1 ¸si x2) în punctele A ¸si B (vezi Figura 1.24).<br />
Exemplu 1.2.3 Un alt exemplu interesant vine din studiul ansamblului arbore<br />
cotit - biela - mecanism cu piston, prezentat într-o formă schematică în figura
72 CAPITOLUL 1. CINEMATICA<br />
A<br />
O<br />
x 2<br />
B<br />
C<br />
Figura 1.24:<br />
1.25. Barele OA ¸si AB sunt conectate în punctul A printr-o legătură articulată,<br />
în timp ce arborele cotit OA se rote¸ste în jurul punctului fix O, iar tija de<br />
legătură AB este conectat la piston în punctul B. Deoarece sistemul se mi¸scă<br />
cu viteze paralele cu un acela¸si plan, putem concluziona că centrul instantaneu<br />
se află pe tija AB, este localizat pe vectorii care sunt normali traiectoriilor<br />
punctelor A ¸si B, ceea ce ne dă că centrul instantaneu se află la intersect¸ia<br />
dreptei OA cu dreapta care trece prin B ¸si care este perpendiculară pe OB.<br />
<br />
Mai mult, dacă θ este unghiul format de OA cu OB, atunci vA = ˙ <br />
<br />
θ<br />
|A − O|.<br />
Dacă C este centrul instantaneu, atunci avem vA/vB = |C − A| / |C − B| ¸si<br />
prin urmare<br />
|A − O| |C − B|<br />
<br />
<br />
vB = <br />
|C − A|<br />
˙ <br />
<br />
θ<br />
.<br />
În cele din urmă, dacă D este punctul situat pe dreapta AB ¸si care se află ¸si pe<br />
dreapta care trece prin O ¸si este perpendiculară pe OB, atunci<br />
<br />
<br />
vB = |D − O| ˙ <br />
<br />
θ<br />
= v ′ D, (1.193)<br />
unde v ′ D<br />
este viteza punctului care este solidar cu manivela ¸si coincide cu D.<br />
Exercit¸iu 1.2.13 Un cilindru circular drept de rază R = 10 se mi¸scă astfel<br />
încât baza sa rămâne întotdeauna în planul fix x1Ox2. Mi¸scarea sistemului de<br />
referint¸ă (O ′ , y1, y2, y3) ata¸sat cilindrului este descrisă prin ecuat¸iile c1(t) =<br />
t 3 + 2, c2(t) = 1 − t 3 , c3(t) = 0 ¸si θ(t) = π(t 2 − t). Determinat¸i: a) coordonatele<br />
x1, x2, x3 ale punctelor P1(y1 = 10, y2 = 0, y3 = 0) ¸si P2(y1 = 0, y2 = 5, y3 = 0)<br />
la momentul t = 2; b) viteza punctului P3(y1 = 1, y2 = 3, y3 = 0).<br />
Solut¸ie. Cilindru execută o mi¸scare plană. Astfel, avem<br />
j1 = cos θi1 + sin θi2, j2 = − sin θi1 + cos θi2, j3 = i3,<br />
x 1
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE73<br />
D<br />
O<br />
x 2<br />
A<br />
Figura 1.25:<br />
¸si, pentru orice punct P (x1, x2, x3), găsim<br />
¸si astfel deducem că<br />
C<br />
B<br />
x1i1 + x2i2 + x3i3 = c1i1 + c2i2 + c3i3 + y1j1 + y2j2 + y3j3<br />
x1 = c1 + y1 cos θ − y2 sin θ, x2 = c2 + y1 sin θ + y2 cos θ, x3 = y3.<br />
Substituind funct¸iile c1, c2, c3 ¸si θ în formula de mai sus, obt¸inem<br />
x1 = t 3 + 2 + y1 cos π(t 2 − t) − y2 sin π(t 2 − t) ,<br />
x2 = 1 − t 3 + y1 sin π(t 2 − t) + y2 cos π(t 2 − t) ,<br />
x3 = y3.<br />
a) La momentul t = 2, avem<br />
x1 = 10 + y1, x2 = −7 + y2, x3 = y3<br />
¸si prin urmare coordonatele punctului P1(y1 = 10, y2 = 0, y3 = 0) sunt x1 = 20,<br />
x2 = −7, x3 = 0, în timp ce coordonatele punctului P2(y1 = 0, y2 = 5, y3 = 0)<br />
sunt x1 = 10, x2 = −2, x3 = 0.<br />
b) Printr-o diferent¸iere directă, obt¸inem următorul câmp al vitezelor:<br />
˙x1 = 3t 2 − π(2t − 1)y1 sin π(t 2 − t) − π(2t − 1)y2 cos π(t 2 − t) ,<br />
˙x2 = −3t 2 + π(2t − 1)y1 cos π(t 2 − t) − π(2t − 1)y2 sin π(t 2 − t) ,<br />
˙x3 = 0.<br />
Înlocuind y1 = 1, y2 = 3, y3 = 0 în relat¸ia de mai sus, obt¸inem viteza punctului<br />
P3 :<br />
v = 3t 2 − π(2t − 1) sin π(t 2 − t) − 3π(2t − 1) cos π(t 2 − t) i1 +<br />
+ −3t 2 + π(2t − 1) cos π(t 2 − t) − 3π(2t − 1) sin π(t 2 − t) i2.<br />
x 1
74 CAPITOLUL 1. CINEMATICA<br />
y 2<br />
O<br />
x 2<br />
O ′<br />
H<br />
y 1<br />
Figura 1.26:<br />
γ<br />
γ<br />
Exercit¸iu 1.2.14 Găsit¸i centrul instantaneu al rigidului a cărui mi¸scare este<br />
paralelă cu un plan fix π.<br />
Solut¸ie. Alegem reperul fix astfel încât planul x1Ox2 să coincidă cu planul<br />
dat π. Fie O ′ un punct fix al corpului rigid. Alegem sistemul de referint¸ă ata¸sat<br />
rigidului astfel încât planul y1Oy2 să fie paralel cu planul π. Astfel, viteza unui<br />
punct generic P al corpului rigid este<br />
vP (t) = vO ′ + ω × (P − O′ ).<br />
Dacă P coincide cu centrul instantaneu C, atunci vC(t) = 0 ¸si prin urmare<br />
avem<br />
ω × (C − O ′ ) = −vO ′.<br />
Astfel, obt¸inem<br />
−ω × vO ′ = ω × [ω × (C − O′ )] ,<br />
¸si prin urmare, folosind formula de dezvoltare a dublului produs vectorial, găsim<br />
[ω · (C − O ′ )] ω − ω 2 (C − O ′ ) = −ω × vO ′.<br />
Pentru că ω este perpendicular pe planul π, în timp ce (C − O ′ ) este paralel cu<br />
π, atunci ω · (C − O ′ ) = 0. Centrul C ne este dat de relat¸ia de mai sus ca fiind<br />
′<br />
x 1<br />
C − O ′ = 1<br />
ω × vO ′. (1.194)<br />
ω2 1.2.14 Traiectorii în coordonate polare<br />
Să considerăm o mi¸scare în care curba γ ′ a planului în mi¸scare (O ′ , y1, y2) rulează<br />
pe curba γ a planului fix (O, x1, x2) (Figura 1.26). Presupunem că cele două<br />
curbe sunt destul de regulate ¸si că admit o tangentă comună în punctul de<br />
intersect¸ie H. Vom spune că γ ′ se mi¸scă pe γ.
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE75<br />
Definit¸ie 1.2.25 Viteza acelui punct al curbei γ ′ care coincide la orice moment<br />
cu punctul de contact H poartă numele de viteza de alunecare. Dacă viteza de<br />
alunecare devine zero, spunem că γ ′ se rostogole¸ste fără alunecare pe γ.<br />
Teoremă 1.2.8 Dacă, în timpul mi¸scării, o curbă γ ′ se rostogole¸ste pe γ,<br />
atunci punctul ce apart¸ine curbei γ ′ care coincide cu punctul de contact H are<br />
viteza orientată spre dreapta tangentă la curba γ în punctul H.<br />
Demonstrat¸ie. Punctul de contact H al celor două curbe γ ¸si γ ′ se mi¸scă<br />
în timp de-alungul curbei γ sau de-alungul curbei γ ′ , acest lucru depinzând<br />
de observator, dacă este conenctat cu sistemul de referint¸dă (O, x1, x2) sau<br />
(O ′ , y1, y2). Din Teorema Compunerii Vitezelor, obt¸inem<br />
va (H) = vr (H) + vτ (H) , (1.195)<br />
unde va (H) ¸si vr (H) sunt vectorii viteză absolută ¸si respectiv viteză relativă<br />
ai punctului H, în timp ce vτ (H), este viteza de transport a lui H ¸si reprezintă<br />
viteza unui punct ce apart¸ine figurii în mi¸scare care coincide cu punctul de<br />
contact H, numită viteză de alunecare. Mai mult, pentru că va ¸si vr sunt<br />
orintat¸i după tangenta la curbele γ si ¸ γ ′ în punctul H, vτ ar trebui să aibe aceea¸si<br />
direct¸ie, ¸si, în consecint¸ă, ar trebui să fie direct¸ionate de-alungul tangentei la γ<br />
în H.<br />
Exemplu 1.2.4 În multe probleme, este posibil să determinăm centru instantaneu<br />
utilizând concluzia teoremei precedente. Drept exemplu, considerăm o bară<br />
AB sust¸inută de o axă în punctul A ¸si de un cerc fix de rază R ¸si cu centrul O,<br />
presum arată Figura 1.27.<br />
Întrucât, în timpul mi¸scării, bara AB alunecă pe cerc, vectorul viteză al<br />
punctului H de pe bara AB care coincide cu punctul de contact, are aceea¸si<br />
direct¸ie cu AB. În acest fel, centrul instantaneu C este punctul de intersect¸ie a<br />
dreptei OH cu normala la axa de sust¸inere în punctul A.<br />
În timpul mi¸scării, centrul instrantaneu î¸si schimbă pozit¸ia în raport cu ambele<br />
sisteme de referint¸ă, descriind două parabole distincte.<br />
Definit¸ie 1.2.26 Locul geometric al pozit¸iilor ocupate de centrul instantaneu,<br />
în timpul mi¸scării plane a unui rigid, fat¸ă de sistemul fix de referint¸ă (O, x1, x2)<br />
poartă numele de bază. Locul geometric al pozit¸iilor ocupate de centrul instantaneu,<br />
în timpul mi¸scării plane a unui rigid, fat¸ă de sistemul de referint¸ă mobil<br />
(O ′ , y1, y2) este o curbă numită ruletă. Baza ¸si ruleta se numesc traiectorii<br />
polare.<br />
În general, baza ¸si ruleta sunt două curbe, aflate în planele O, x1, x2 ¸si respectiv<br />
O ′ , y1, y2. Un exemplu interesant în acest sens provine din studiul mic¸scării<br />
reprezentate în Figura 1.24. Centrul instantaneu C rămâne la o distant¸ă constantă<br />
fat¸ă de originea O în raport cu (O, x1, x2), deoarece |C − O| = |A − B|.<br />
De aceea, baza va fi cercul cu centrul în O ¸si rază egală cu |A − B| (Figura 1.28).
76 CAPITOLUL 1. CINEMATICA<br />
B<br />
B<br />
O<br />
H<br />
O<br />
x 2<br />
A<br />
Figura 1.27:<br />
′ x2 A<br />
Figura 1.28:<br />
C<br />
C<br />
Pentru a determina ruleta, care este locul geometric al pozit¸iilor ocupate de centrul<br />
instantaneu fat¸ă de sistemul de referint¸ă (B, y1, y2) (F igura1.28) ata¸sat de<br />
bara AB, putem observa că AB este ipotenuza triunghiului drept ABC care<br />
se mi¸scă astfel încât latura AB să rămână fixă iar vârful C să varieze. Astfel,<br />
ruleta este cercul cu diametrul AB.<br />
Teoremă 1.2.9 În timpul mi¸scării plane a rigidului, ruleta se rostogole¸ste fără<br />
să alunece peste curba bază.<br />
Demonstrat¸ie. La fiecare moment, baza ¸si ruleta au punctul C în comun.<br />
În plus, din (1.190), viteza vτ (C) a punctului figurii mobile care coincide cu<br />
C se anulează. De aceea, dacă va(C) ¸si vr(C) sunt vectorii viteză absolută ¸si<br />
′ x1 x 1
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE77<br />
viteza relativă a centrului instantaneu, atunci avem<br />
va(C) = vr(C). (1.196)<br />
Deoarece vitezele va ¸si vr ar trebui să fie tangente curbelor bază ¸si ruletă, în C,<br />
deducem din (1.196) că aceste două curbe sunt tangente. De aceea, în timpul<br />
mi¸scării, ruleta se rostogole¸ste peste curba bază. În cele din urmă, pentru că<br />
vτ (C) = 0, viteza de alunecare se anulează. Astfel, ruleta se rostogole¸ste fără<br />
să alunece peste curba bază.<br />
Observat¸ie 1.2.23 Baza ¸si ruleta sunt singurele curbe fixe în sistemele de<br />
referint¸ă (O, x1, x2) ¸si respectiv (O ′ , y1, y2), care se rostogolesc fără să alunece<br />
una pe alta. De fapt, dacă există alte două curbe care au acelea¸si proprietăt¸i,<br />
viteza punctului care coincide cu punctul de contact din figura în mi¸scare fără<br />
alunecare, se anulează. Totu¸si, din (1.190), doar centrul instantaneu are această<br />
proprietate.<br />
Exemplu 1.2.5 Considerăm un disc care se rostogole¸ste fără alunecare peste o<br />
axa (Figura 1.29). Această mi¸scare poartă numele de cicloida, pentru că fiecare<br />
punct descrie o cicloidă. Din Observat¸ia 1.2.23, pentru această mi¸scare, axa ¸si<br />
discul vor fi baza ¸si respectiv ruleta. Astfel, la fiecare moment, centrul instantaneu<br />
C va coincide cu punctul de contact dintre disc ¸si axa pe care discul se<br />
mi¸scă. Să determinăm numărul gradelor de libertate ale acestui sistem. Notăm<br />
prin x distant¸a dintre C ¸si O, ¸si prin θ unghiul pe care diametrul AA ′ al discului<br />
îl formează cu vectorul (C − O ′ ). Evident, x ¸si θ determină pozit¸ia discului, dar<br />
este u¸sor să observăm că impunerea condit¸iei de rostogolire fără alunecare duce<br />
la o relat¸ie între x ¸si θ. Viteza punctului O ′ este dată de<br />
vO ′ = ˙xi1. (1.197)<br />
În plus, ¸stiind că discul se află într-o stare de mi¸scare de rotat¸ie în jurul lui C,<br />
putem obt¸ine următoarea expresie pentru vO ′<br />
vO ′ = ω × (O′ − C), (1.198)<br />
unde ω este acela¸si vector care apare în formula fundamentală<br />
vP = vO ′ + ω × (P − O′ ).<br />
Prin urmare, ω reprezintă viteza unghiulară a mi¸scării discului fat¸ă de sistemul<br />
de referint¸ă cu originea O ′ ¸si axele paralele sau fixe fat¸ă de (x1, x2). În acest<br />
caz, avem<br />
ω = ˙ θi3. (1.199)<br />
De aceea, comparând (1.197), (1.198) ¸si (1.199), obt¸inem<br />
˙x = R ˙ θ, (1.200)
78 CAPITOLUL 1. CINEMATICA<br />
O<br />
x 2<br />
x<br />
A<br />
θ<br />
C<br />
P<br />
′ O<br />
′<br />
A<br />
Figura 1.29:<br />
unde R este raza discului. Această relat¸ie reprezintă constrângerea legată de<br />
rulare fără alunecare ¸si exprimă o constrângere nonholonomică. De fapt, relat¸ia<br />
(1.200) poate fi integrată u¸sor ¸si astfel obt¸inem ecuat¸ia<br />
x 1<br />
x = Rθ + x0, (1.201)<br />
unde x0 este o constanta convenabilă. Din (1.201) vedem că relat¸ia dintre θ ¸si<br />
x este de tip holonomic. Prin urmare, discul care se rostogole¸ste fără alunecare<br />
are doar un singur grad de libertate. Dacă în schimb considerăm o sferă care<br />
se rostogole¸ste fără alunecare pe un plan, este posibil să demonstrăm că ecuat¸ia<br />
diferent¸ială care rezultă din definit¸ia constrângerii de rostogolire fără alunecare<br />
nu este integrabilă. Astfel, nu este posibil să efectuăm acelea¸si operat¸ii ca în<br />
cazul discului pentru a trece de la (1.200) la (1.201). Prin urmare, acest sistem<br />
nu este holonomic.<br />
Exercit¸iu 1.2.15 Un cilindru de rază a (> 0) se mi¸scă pe un plan orizontal.<br />
Găsit¸i baza ¸si ruleta acestei mi¸scări.<br />
Solut¸ie. Presupunând că cilindrul se mi¸scă pe planul x1Ox2, de-alungul axei<br />
pozitive x2 cu viteza v0 (viteza centrului O ′ a sect¸iunii transversale a cilindrului<br />
situat în planul x2Ox3), rostogolindu-se în jurul lui O ′ , viteza unghiulară ω.<br />
Alegem sistemul de referint¸ă ata¸sat corpului cu originea în O ′ ¸si x3 = a ca plan<br />
y1O ′ y2. Făcând aceatsă alegere, mi¸scarea este în planul x2Ox3. Mai mult, avem<br />
ω = −ωi1, vO ′ = v0i2, ¸si din această cauză, obt¸inem ω × vO ′ = −ωv0i3. Astfel,<br />
relat¸ia (1.194), care ne dă centrul instantaneu, deducem<br />
C − O ′ = − v0<br />
ω i3.<br />
Observăm că această relat¸ie implică faptul că |C − O ′ | = v0<br />
¸si prin urmare<br />
ω<br />
ruleta este cercul de rază v0<br />
ω . Presupunând că nu există alunecare, v0 = aω ¸si<br />
astfel ruleta este circumferint¸a cilindrului.<br />
În raport cu reperul de referint¸ă fix, relat¸ia C − O ′ = − v0<br />
ω i3 determină<br />
x2 − x 0 2 = 0, x3 − x 0 3 = − v0<br />
ω ,<br />
unde O ′ − O = x0 2i2 + x0 3i3, x0 3 = a. De aceea, centrul instantaneu descrie<br />
dreapta de ecuat¸ie x3 = x0 3 − v0<br />
ω care se află în planul x2Ox3, adică baza este<br />
o dreaptă paralecă cu axa x2.
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE79<br />
x 2<br />
O<br />
dw<br />
dt<br />
− a M<br />
× (A − M)<br />
− ω 2 (A − M)<br />
A<br />
a M<br />
M<br />
dω<br />
dt h − ω2 r<br />
P<br />
Figura 1.30:<br />
dw<br />
dt<br />
× (P − M)<br />
− ω 2 (P − M)<br />
1.2.15 Accelerat¸ia unei mi¸scări rigide plane<br />
Pentru a determina distribut¸ia accelerat¸iilor punctelor unei figuri plane care se<br />
mi¸scă în plan, pornim de la relat¸ia (1.146), pe care o scriem în raport cu un<br />
punct M a figurii plane (Figura 1.30), sub următoarea formă:<br />
aP = aM + dω<br />
dt × (P − M) − ω2 (P − M). (1.202)<br />
Teoremă 1.2.10 Pentru o mi¸scare rigidă plană, distribut¸ia accelerat¸iilor este<br />
aceea¸si ca pentru mi¸scarea în jurul unui punct A, numit pol al accelerat¸iei, cu<br />
viteza unghiulară ω ¸si accelerat¸ia unghiulară dω<br />
dt .<br />
Demonstrat¸ie. Deoarece ω are o direct¸ie constantă, dω/dt are aceea¸si<br />
direct¸ie ca ω. Astfel, dω<br />
dt × (P − M) apart¸ine aceluia¸si plan cu figura ¸si este<br />
P −M<br />
ortogonal pe (P − M), a cărui direct¸ie o vom nota cu h. Stabilim r = |P −M| ¸si<br />
observăm că diferent¸a<br />
dω<br />
dt × (P − M) − ω2 <br />
dω<br />
(P − M) = |P − M|<br />
dt h − ω2 <br />
r (1.203)<br />
este un vector care are modulul proport¸ional cu |−M| ¸si are direct¸ia astfel încât<br />
unghiul pe care-l formează cu (P − M) să nu depindă de P , deoarece dω/dt<br />
¸si ω 2 nu depind de P . Vectorii h ¸si r sunt totdeauna mutul ortogonali. De<br />
aceea, există un punct A al figurii plane astfel încât expresia (1.203) calculată<br />
în raport cu acel punct să fie egală cu −aM . Astfel, accelerat¸ia aA a acelui<br />
punct se anulează. Dacă stabilim O ′ = A în relat¸ia (1.146), obt¸inem<br />
x 1<br />
aP = dω<br />
dt × (P − A) − ω2 (P − A). (1.204)<br />
Distribut¸ia accelerat¸iilor este aceea¸si ca în cazul mi¸scării de rotat¸ie în jurul lui<br />
A cu viteza unghiulară ω ¸si cu accelerat¸ia unghiulară dω/dt.
80 CAPITOLUL 1. CINEMATICA<br />
y 3<br />
x 1<br />
y 1<br />
O<br />
x 3<br />
Figura 1.31:<br />
Observat¸ie 1.2.24 Merită să ment¸ionăm în acest context că plul accelerat¸iei<br />
în general nu coincide cu centrul instantaneu C, deoarece accelerat¸ia lui C nu<br />
este în mod obligatoriu zero.<br />
1.2.16 Mi¸scarea unui corp rigid cu un punct fix<br />
Conuri Poinsot<br />
Este cunoscut că starea de mi¸scare a unui corp rigid cu un punct fix O este de<br />
rotat¸ie, cu axa instantanee de rotat¸ie trecând prin O. Pe parcursul mi¸scării,<br />
această axă î¸si schimbă orientarea fat¸ă de ambele sisteme de referint¸ă, cel fix<br />
(O, x1, x2, x3) ¸si fat¸ă de cel mobil. Deoarece axa trece întotdeauna prin O,<br />
descrie pe parcurs două conuri, unul fix în (O, x1, x2, x3), ¸si unul de asemenea<br />
fix, dar în (O, y1, y2, y3), pe care le vom numi conurile Poinsot (Figura 1.31).<br />
Prin utilizarea unui rat¸ionament similar cu cel din studiul mi¸scărilor rigide<br />
plane, putem demonstra că, pe parcursul mi¸scării, cele două conuri se rostogolesc<br />
fără alunecare unul pe celălalt, ¸si generatoarea pe care o au în comun, la fiecare<br />
moment, este axa instantanee de rotat¸ie.<br />
x 2<br />
y 2
1.2. CINEMATICA SISTEMELOR MATERIALE S¸I CORPURILOR RIGIDE81<br />
x 1<br />
O<br />
x 3<br />
w 2<br />
w 1<br />
Mi¸scarea de precesie<br />
Figura 1.32:<br />
Presupunând că starea de mi¸scare a corpului rigid cu un punct fix este reprezentată<br />
de suma a două stări de mi¸scare de rotat¸ie, dacă O este un punct fix, atunci<br />
x 2<br />
vP = ω1 × (P − O) + ω2 × (P − O), (1.205)<br />
unde ω1 este un vector a cărui direct¸ie este de-alungul unei dreptei f fixe în corp<br />
¸si care trece prin O, pe care o vom numi axa figurii. Vectorul ω2 are direct¸ia<br />
de-alungul unei drepte p care trece prin O ¸si atribuită în raport cu reperul<br />
(O, x1, x2, x3) fix în spat¸iu, care poartă numele de axă de precesie. Dacă ω1 ¸si<br />
ω2 au modulul constant, spunem că mi¸scarea este o precesie regulată.<br />
Pentru punctele P ∗ ale axei din figură, relat¸ia (1.205) ia forma (Figura 1.32)<br />
vP ∗ = ω2 × (P ∗ − O)<br />
¸si, în consecint¸ă, aceste puncte se rotesc în jurul axei de precesie. Prin urmare,<br />
în timpul mi¸scării, corpul se rote¸ste în jurul axei din figură cu viteza unghiulară<br />
ω1 iar aceste axe se rotesc ¸si ele în jurul axei de precesie.<br />
După cum se poate vedea în relat¸ia (1.205), starea de mi¸scare este de rotat¸ie<br />
în jurul axei care trece prin O ¸si având direct¸ia lui ω1+ ω2. Deoarece ω1 ¸si ω2<br />
au modulele constante¸si unghiul făcut de ω1 ¸si ω2 este de asemenea constant,<br />
unghiurile α ¸si β formate de vectorul ω1+ω2 cu vectorii ω1 ¸si respectiv ω2,, sunt<br />
constante (Figura 1.33). Astfel, în raport cu sistemul de referint¸ă (O, x1, x2, x3)<br />
fix în spat¸iu, axa instantanee de rotat¸ie se rote¸ste în jurul direct¸ie lui ω2 în<br />
acela¸si mod ca ¸si vectorul ω1+ ω2, în timp ce, fat¸ă de reperul ata¸sat corpului,<br />
axa instantanee de rotat¸ie se rote¸ste în jurul directiei lui ω1. De aceea, conurile<br />
lui Poinsot sunt două conuri circulare care au, corespunzător, axa de precesie ¸si<br />
axele figurii drept axe de simetrie.
82 CAPITOLUL 1. CINEMATICA<br />
x 1<br />
O<br />
x 3<br />
w 2<br />
β<br />
α<br />
w 1 + w 2<br />
w 1<br />
Figura 1.33:<br />
Mi¸scarea titirezului ¸si mi¸scarea Pământului sunt exemple de mi¸scări de precesie.<br />
În raport cu un observator cu originea în centrul Pământului ¸si axele<br />
orientate spre stelele fixe, Pământul se rote¸ste în jurul axei sale proprii polare,<br />
care este axa figurii. La rândul său, această axă se rote¸ste în jurul unei axe<br />
ortogonale la planul eclipticii cu viteza unghiulară ω2 = 2π/T , unde perioada<br />
T este aproximativ egală cu 26 000 ani. Această rotat¸ie lentă a axei polare<br />
cauzează a¸sa-numitul fenomen de ”precesie a echinoct¸iului;” dat de, trecerea<br />
Soarelui de-a lungul liniei nodurilor (determinată ca intersect¸ia dintre planul de<br />
ecliptică cu planul Ecuatorului), care este verificată în fiecare an, cu un avans<br />
anumit. Aceasta duce la deplasarea relativă a stelelor fixe ¸si, prin urmare, este<br />
necesar să se actualizeze calendarul în conformitate cu un program special.<br />
x 2
Bibliografie<br />
[1] Abraham, R. H. and Marsden, J. E., Foundations of Mechanics: A Mathematical<br />
Exposition of Classical Mechanics, 2nd ed. New York: Benjamin,<br />
1978.<br />
[2] Abramowitz, M. and Stegun, I., Handbook of Mathematical Functions. New<br />
York, Dover, 1965.<br />
[3] Ames, J. S. and Murnaghan, F. D., Theoretical Mechanics: An Introduction<br />
to Mathematical Physics. New York: Dover, 1958.<br />
[4] Arnold, V. I., Mathematical Methods of Classical Mechanics, 2nd ed. New<br />
York: Springer–Verlag, 1989.<br />
[5] Bowen, R. M. and Wang, C. C., Introduction to Vectors and Tensors, Vol.<br />
1. New York: Plenum, 1976.<br />
[6] Brouwer, D. and Clemence, G. M., Methods of Celestial Mechanics. New<br />
York: Academic Press, 1961.<br />
[7] Cercignani, C., Theory and Application of the Boltzmann Equation. Edinburgh:<br />
Scot. Academic Press, 1975.<br />
[8] Chow, T. L., Classical Mechanics. New York: Wiley, 1995.<br />
[9] Corben, H. C. and Stehle, P., Classical Mechanics, 2nd ed. New York:<br />
Wiley, 1960.<br />
[10] Fabrizio, M., Introduzione alla Meccanica Razionale e ai suoi Metodi<br />
Matematici. Bologna: Zanichelli, 1994.<br />
[11] Fabrizio, M., Elementi di Meccanica Classica. Bologna: Zanichelli, 2002.<br />
[12] Fowles, G. R. and Cassiday, G. L., Analytical Mechanics, 5th ed. Orlando:<br />
Saunders, 1993.<br />
[13] French, A. P., Newtonian Mechanics. New York: W. W. Norton, 1971.<br />
[14] Gantmacher, F. R., Lectures in Analytical Mechanics. Moscow: Mir Publishers,<br />
1970.<br />
83
84 BIBLIOGRAFIE<br />
[15] Graffi, D., Elementi di Meccanica Razionale. Bologna: Pátron, 1970.<br />
[16] Greenwood, T. D., Classical Dynamics. New York: Dover, 1997.<br />
[17] Griffits, J. B., The Theory of Classical Dynamics. Cambridge: Cambridge<br />
University Press, 1985.<br />
[18] Gurtin, M. E., An Introduction to Continuous Mechanics. New York: Academic<br />
Press, 1981.<br />
[19] Hirsch, M. W. and Smale, S., Differential Equations, Dynamical Systems<br />
and Linear Algebra. New York: Academic Press, 1974.<br />
[20] Hunter, S. C., Mechanics of Continuous Media, 2nd ed. Ellis Horwood,<br />
1983.<br />
[21] Kellog, O. D., Foundations of Potential Theory. Berlin: Springer–Verlag,<br />
1929.<br />
[22] Kelvin, W. T. and Tait, P. G., Principles of Mechanics and Dynamics, 2<br />
vols. New York: Dover, 1962.<br />
[23] Kibble, T. W. B. and Berkshire, F. H., Classical Mechanics. Harlow: Longman,<br />
1996.<br />
[24] Kilmister, C. W. and Reeve, J. E., Rational Mechanics.London: Longmans,<br />
1966.<br />
[25] Kittel, C., Knight, W. D. and Ruderman, M. A., Mechanics, 2nd ed. New<br />
York: McGraw–Hill, 1973.<br />
[26] Kleppner, D. and Kolenkow, R. J., An Introduction to Mechanics. New<br />
York: McGraw–Hill, 1973.<br />
[27] Knops, R. J. and Wilkes, E. W., Theory of Elastic Stability. In Encyclopedia<br />
of Physics, vol. VIa/3, (C. A. Truesdell, ed.), Berlin: Springer–Verlag, 125–<br />
302, 1973.<br />
[28] Knudsen, J. M. and Hjorth, P. G., Elements of Newtonian Mechanics. New<br />
York: Springer–Verlag, 1995.<br />
[29] Lagrange, J. L., Mécanique Analitique, 4th ed., 2 vols. Paris: Gauthier–<br />
Villars et fils, 1888–89.<br />
[30] Lamb, H., Dynamics, 2nd ed. London: Cambridge University Press, 1961.<br />
[31] Lanczos, C., The Variational Principles of Mechanics, 4th ed. New York:<br />
Dover, 1986.<br />
[32] Landau, L. D. and Lifshitz, E. M., Mechanics 3rd ed. Oxford: Pergamon<br />
Press, 1976.
BIBLIOGRAFIE 85<br />
[33] Levi Civita, T. and Amaldi, U., Lezioni di Meccanica Razionale. Bologna:<br />
Zanichelli, 1927.<br />
[34] Liboff, R. L., Introduction to the Theory of Kinetic Equations. New York:<br />
John Wiley and Sons, Inc., 1969.<br />
[35] Mach, E., The Science of Mechanics: A Critical and Historical Exposition<br />
of its Principles.Chicago, IL: Open Court, 1893.<br />
[36] Macmillan, W. D., Dynamics of Rigid Bodies. New York: Dover, 1960.<br />
[37] Marion, J. B. and Thornton, S. T., Classical Dynamics of Particles and<br />
Systems, 4th ed. Philadelphia: Saunders, 1995.<br />
[38] Marsden, J. E. and Ratiu, T. S., Introduction to Mechanics and Symmetry:<br />
A Basic Exposition of Classical Mechanical Systems. New York: Springer–<br />
Verlag, 1994.<br />
[39] Mesarovic, M. D. and Takahara, Y., General Systems Theory: Mathematical<br />
Foundations. New York: Academic Press, 1975.<br />
[40] Moulton, F. R., An Introduction to Celestial Mechanics, 2nd rev. ed. New<br />
York: Dover, 1970.<br />
[41] Natanson, I. P., Theory of Functions of a Real Variable.New York: Ungar,<br />
1955.<br />
[42] Nusse, H. E. and Yorke, J. A., Dynamics: Numerical Explorations. New<br />
York: Springer–Verlag, 1994.<br />
[43] Osgood, W. F., Mechanics. New York: Macmillan, 1937.<br />
[44] Pars, L. A., A Treatise on Analytical Dynamics. New York: Wiley, 1965.<br />
[45] Percival, I. and Richards, D., Introduction to Dynamics. Cambridge: Cambridge<br />
University Press, 1987.<br />
[46] Pollard, H., Mathematical Introduction to Celestial Mechanics. Englewood<br />
Cliffs, NJ: Prentice Hall, 1966.<br />
[47] Riley, W. F. and Sturges, L. D., Engineering Mechanics: Statics, 2nd ed.<br />
New York: Wiley, 1966.<br />
[48] Rosenauer, N. and Willis, A. H., Kinematics of Mechanisms. New York:<br />
Dover, 1967.<br />
[49] Routh, E. J., A Treatise on the Dynamics of a System of Rigid Bodies:<br />
Elementary Part, 7th ed., rev. and enl. New York: Dover, 1960.<br />
[50] Siegel, C. L. and Moser, J. K., Lectures on Celestial Mechanics, 2nd ed.<br />
Berlin: Springer–Verlag, 1995.
86 BIBLIOGRAFIE<br />
[51] Slater, J. C. and Frank, N. H., Mechanics. New York: McGraw–Hill, 1947.<br />
[52] Sommerfeld, A., Mechanics. New York: Academic Press, 1952.<br />
[53] Stiefel, E. L. and Scheifele, G., Linear and Regular Celestial Mechanics:<br />
Perturbed Two–Body Motion, Numerical Methods, Canonical Theory. New<br />
York: Springer–Verlag, 1971.<br />
[54] Symon, K. R., Mechanics, 3rd ed. Reading, MA: Addison–Wesley, 1971.<br />
[55] Synge, J. L. and griffith, B. A., Principles of Mechanics, 3rd ed. New York:<br />
McGraw–Hill, 1959.<br />
[56] Thornton, M., Classical Dynamics of Particles and Systems. Philadelphia:<br />
Harcourt Brace College Publishers, 1995.<br />
[57] Timoshenko, S. and Young, D. H., Engineering Mechanics, 4th ed, 2 vols.<br />
New York: McGraw–Hills, 1956.<br />
[58] Tisserand, F. F., Traité de Mécanique Celeste, 4 vols. Paris: Gauthier–<br />
Villars et fils, 1889–96.<br />
[59] Verhulst, F., Nonlinear Differential Equations and Dynamical Systems.<br />
Berlin: Springer-Verlag, 1990.<br />
[60] Whittaker, E. T., Analytical Dynamics, 4th ed. Cambridge: Cambridge<br />
University Press, 1959.<br />
[61] Williams, D., Elements of Mechanics. Oxford: Oxford University Press,<br />
1997.
Glosar<br />
A<br />
absolute acceleration, 61<br />
– angular velocity, 65<br />
– concept, 2<br />
– force, ??<br />
– temperature, ??<br />
– velocity, 61<br />
accelerated motion, 11<br />
acceleration of a point, 9<br />
accoustic tensor, ??<br />
active force, ??<br />
– loads, ??<br />
Almansi strain tensor, ??<br />
altitude of shooting, ??<br />
amplitude, 25<br />
angle of static friction, ??<br />
angular acceleration, 14<br />
– frequency, 25, ??<br />
– momentum, ??, ??<br />
– transport velocity, 66<br />
– velocity vector, 13<br />
anisotropic elastic material, ??<br />
aperiodic motion with critical damping,<br />
??<br />
Appel’s equations, ??<br />
areal velocity, 17<br />
arm of power, ??<br />
– of the resistance, ??<br />
asymptotic stable equilibrium state, ??<br />
attractor, ??<br />
axial momentum, ??, ??<br />
axis of oscillation, ??<br />
– of precession, 81<br />
– of rotation, 42<br />
87<br />
B<br />
balance law of angular momentum, ??<br />
– law of linear momentum, ??<br />
basis, 75<br />
bending of a beam, ??<br />
Bernoulli theorem, ??<br />
Betti’s reciprocal theorem, ??<br />
bifurcation, ??<br />
bilateral constraint, 31<br />
Binet formula, 20<br />
binormal unit vector, 10<br />
body cone, ??, ??<br />
Boltzmann equation, ??<br />
C<br />
canonical equations, ??<br />
– invariants, ??<br />
– systems, ??<br />
– transformation, ??<br />
Cardan’s suspension, ??<br />
cardinal equations, ??, ??, ??<br />
Cauchy hypothesis, ??<br />
– stress, ??, ??<br />
causal dynamic system, ??<br />
central force, ??<br />
– motion, 19<br />
centre of mass, ??<br />
– of motion, 19<br />
centrifugal force, ??, ??<br />
centripetal acceleration, 11<br />
change in observer, ??<br />
chaos, ??<br />
characteristic polynomial, ??<br />
circular motion, 22<br />
– and uniform motion, 22
88 GLOSAR<br />
Clausius–Duhem inequality, ??<br />
closed cycle, ??, ??<br />
coefficient of dynamical friction, ??<br />
– of friction strength, ??<br />
– of static friction, ??<br />
– of static sliding friction, ??<br />
complementary acceleration, 61<br />
complete canonical transformation, ??<br />
compound pendulum, ??<br />
compression, ??<br />
conditional stability, ??<br />
cone of static friction, ??, ??<br />
conjugate variables, ??<br />
conservation of energy, ??, ??<br />
conservative force, ??<br />
– Hamiltonian system, ??<br />
– system of constitutive forces, ??<br />
constitutive equation, ??, ??<br />
– force, ??, ??, ??<br />
constrained material point, ??<br />
constraint, 31<br />
– of mobility, 32<br />
– of position, 32<br />
– reaction, ??<br />
continuation, ??<br />
continuity equation in Lagrangian form,<br />
??<br />
– equation in spatial form, ??<br />
continuous body, ??<br />
– material system, ??<br />
coordinates, 4<br />
Coriolis acceleration, 61<br />
– force, ??, ??<br />
couple, ??<br />
curvature, 10<br />
curvilinear abscissa, 5<br />
– coordinate, 30, ??, ??<br />
cyclic coordinate, ??<br />
cycloidal motion, 77<br />
cylindrical coordinates, ??<br />
D<br />
D’Alembert principle, ??<br />
damped aperiodic motion, ??<br />
– oscillatory motion, ??<br />
damping parameter, ??<br />
deformation, ??<br />
degree of freedom, 33<br />
density of mass, ??<br />
– of probability, ??<br />
dependent of time constraint, 32<br />
deviation moment, ??<br />
– of a heavy point, ??<br />
direct motion, 7, 22<br />
Dirichlet theorem, ??<br />
discrete material system, ??<br />
displacement, ??<br />
dissipation principle, ??<br />
double pendulum, ??<br />
duration of shooting, ??<br />
dynamic state, ??<br />
dynamical process, ??<br />
E<br />
effective force, ??<br />
eigenvalue, ??<br />
eigenvector, ??<br />
elastic solid, ??<br />
elasticity tensor, ??<br />
elementary displacement, 9<br />
– work, ??<br />
ellipsoid of inertia, ??<br />
energy equation, ??<br />
entropy, ??<br />
equilibrium configuration, ??<br />
– position, ??<br />
equivalent reference frames, 64<br />
– representations, ??<br />
– states, ??<br />
Eulerian angles, 45<br />
– coordinates, ??<br />
– dynamical process, ??<br />
Euler’s equations, ??, ??, ??<br />
extension, ??<br />
extra stress, ??<br />
F<br />
first integral, ??, ??, ??<br />
– Piola–Kirchhoff stress, ??
GLOSAR 89<br />
– principle of thermodynamics, ??<br />
fixed axoide, 55<br />
flexion point, ??<br />
flow region, ??<br />
forced oscillation, ??<br />
Foucault’s pendulum, ??<br />
frequency, 23<br />
friction force, ??<br />
– strength, ??<br />
function of distribution, ??<br />
– of transition of states, ??, ??, ??<br />
functional of response, ??, ??<br />
fundamental equation, ??<br />
– formula of kinematics of rigid systems,<br />
52<br />
– harmonic, ??<br />
G<br />
Galileian invariance principle, ??<br />
– observer, ??<br />
general integral, ??, ??<br />
generalized coordinates, 33, ??<br />
– integral of energy, ??<br />
– moments, ??<br />
– potential, ??<br />
– shear modulus, ??<br />
generating function, ??<br />
geodesic curve, ??<br />
geometric constraint, 32<br />
gravitational acceleration, ??<br />
Green strain tensor, ??<br />
group of translations, 2<br />
gyroscope, ??<br />
gyroscopic axis, ??<br />
– phenomena, ??<br />
– structure, ??<br />
– system, ??<br />
– system of forces, ??<br />
gyrostat, ??<br />
H<br />
Hamilton–Jacobi equation, ??<br />
– method, ??<br />
– reduced equation, ??<br />
– time dependent equation, ??<br />
Hamiltonian, ??, ??<br />
– action, ??<br />
– system, ??<br />
Hamilton’s equations, ??<br />
– principal function, ??<br />
– principle, ??<br />
harmonic oscillator, ??<br />
– oscillatory motion, 25<br />
heat equation, ??<br />
helical motion, 29<br />
– state of motion, 50<br />
holonomic constraint, 32<br />
– material system, 33<br />
homogeneous deformation, ??<br />
– elastic material, ??<br />
Hopf’s bifurcation, ??<br />
Huygens theorem, ??<br />
hydrostatic pressure, ??<br />
hyperelastic material, ??<br />
hysteresis cycle, ??<br />
I<br />
ideal fluid, ??<br />
ignorable coordinate, ??<br />
impulse of a force, ??, ??<br />
impulsive motion, ??<br />
incompressible material, ??<br />
independent of time constraint, 32<br />
inertia force, ??<br />
– matrix, ??<br />
inertial observer, ??<br />
inertia’s motion, ??<br />
infinitesimal rigid displacement, ??<br />
– theory, ??<br />
inflexion point, ??<br />
inhomogeneous elastic material, ??<br />
initial phase, 25<br />
Input–Output system, ??<br />
instantaneous angular velocity, 49<br />
– centre of rotation, 70<br />
integrable system, ??<br />
internal variable, ??, ??<br />
intrinsic vector triplet, 11, ??<br />
irreversible virtual displacement, ??
90 GLOSAR<br />
isochoric deformation, ??<br />
– dynamical process, ??<br />
isolated material point, ??<br />
isotropic elastic material, ??<br />
J<br />
Jacobi integral, ??<br />
Jacobi’s identity, ??<br />
jump, ??<br />
K<br />
Kepler’s laws, ??<br />
kinematical viscosity, ??<br />
kinetic constraint, 32, ??, ??<br />
– energy, ??, ??, ??<br />
– moments, ??<br />
– roto–translational state of motion,<br />
50<br />
– state, 49<br />
König theorem, ??<br />
L<br />
Lagrange generalized forces, ??<br />
Lagrange’s equations, ??<br />
– multipliers, ??<br />
Lagrangian, ??, ??<br />
– coordinates, 33, ??<br />
Lamé moduli, ??<br />
left Cauchy strain tensor, ??<br />
– Cauchy–Green strain tensor, ??<br />
limit cycle, ??<br />
linear momentum, ??, ??<br />
– theory, ??<br />
Liouville’s operator, ??<br />
– theorem, ??<br />
Lipschitz function, ??<br />
Lissajous curve, ??<br />
list of principal invariants, ??<br />
local integrability, ??<br />
localization of a body, 3<br />
longitudinal wave, ??<br />
Lyapunov function, ??<br />
– theorem, ??<br />
Lyapunov’s first method, ??<br />
– second method, ??<br />
M<br />
Mach number, ??<br />
material coordinates, ??<br />
– point, 3, ??<br />
– state, ??<br />
– system, 3<br />
maximum point, ??<br />
Maxwell–Boltzmann distribution, ??<br />
mean kinetic energy, ??<br />
mechanical system, ??<br />
meter, 2<br />
minimum point, ??<br />
mobile axoide, 55<br />
modulus of compression, ??<br />
moment, 2<br />
– of arrest, 7<br />
– of inertia, ??<br />
– of momentum, ??, ??<br />
motion of precession, 82<br />
Mozzi’s theorem, 53<br />
N<br />
Navier–Stokes equations, ??<br />
Newtonian attraction force, ??<br />
– fluid, ??<br />
Newton’s first law, ??<br />
– second law, ??<br />
– third law, ??<br />
nonholonomic constraint, 32<br />
– material system, 34<br />
nonlinear oscillations, ??<br />
normal acceleration, 11<br />
– distribution, ??<br />
– force, ??<br />
number of degrees of freedom, 33<br />
nutation, 45, ??<br />
O<br />
osculating circle, 10<br />
– plane, 9
GLOSAR 91<br />
P<br />
part of a body, ??<br />
percussion, ??<br />
perfectly smooth constraint, ??<br />
perimeter of support, ??<br />
period of motion, 22<br />
periodic motion, 22<br />
phase diagram, ??<br />
– plane, ??<br />
physical pendulum, ??<br />
pinnacle, ??<br />
plane of material symmetry, ??<br />
– polar coordinates, 12<br />
– motion, 12<br />
– rigid motion, 55<br />
Poincaré’s integral, ??<br />
Poinsot cones, 80<br />
Poinsot’s motion, ??<br />
point of application, ??<br />
points, 1<br />
Poiseuille motion, ??<br />
Poisson bracket, ??<br />
Poisson’s distribution, ??<br />
– formulae, 51<br />
– ratio, ??, ??<br />
polar angle, 12<br />
– distance, 12<br />
– trajectories, 75<br />
pole of acceleration, 57, 79<br />
positional force, ??<br />
positive definite elasticity tensor, ??<br />
potential, ??, ??<br />
– energy, ??<br />
– flow, ??<br />
power, ??, ??<br />
precession, 45, ??<br />
prescribed external forces, ??<br />
pressure, ??<br />
principal axes of inertia, ??<br />
– direction, ??<br />
– frequency, ??<br />
– moments of inertia, ??<br />
– normal unit vector, 10<br />
– oscillation, ??<br />
– stress, ??<br />
principle of conservation mass, ??<br />
– of determinism, ??<br />
– of the gyroscopic effect, ??<br />
– of the constraint reactions, ??<br />
– of virtual work, ??<br />
problem of torsion, ??<br />
products of inertia, ??<br />
proper rotation, 45<br />
pulsation, 25, ??<br />
pure shear, ??<br />
– tension, ??<br />
R<br />
radial acceleration, 14<br />
– velocity vector, 13<br />
radius of curvature, 10<br />
range of shooting, ??<br />
rectilinear and uniform motion, 8<br />
– translation motion, 41<br />
reduced length, ??<br />
– mass, ??<br />
reference frame, 3, 4<br />
– system, 3<br />
regular causal stationary dynamic system,<br />
??<br />
relative acceleration, 61<br />
– angular velocity, 65<br />
– equilibrium position, ??<br />
– velocity, 61<br />
resistance, ??<br />
– force, ??<br />
resonance, ??<br />
rest position, ??<br />
retarded motion, 11<br />
retrograde motion, 7, 22<br />
reversible virtual displacement, ??<br />
Reynolds number, ??<br />
rheonomic constraint, 32<br />
right Cauchy strain tensor, ??<br />
– Cauchy–Green strain tensor, ??<br />
rigid body, 3<br />
– body with a fixed axis, ??<br />
– body with a fixed point, 36, 80<br />
rigid material system, ??<br />
rolling without sliding, 75
92 GLOSAR<br />
rotational kinetic state, 49<br />
– motion, 42<br />
– state of motion, 49<br />
roto–translational motion, 43<br />
– state of motion, 50<br />
roulette, 75<br />
Routh’s equations, ??<br />
S<br />
Saint–Venant problem, ??<br />
Saint–Venant’s compatibility conditions,<br />
??<br />
scalar multiplication, 2<br />
scleronomic constraint, 32<br />
second, 2<br />
– cosmic speed, ??<br />
– principle of thermodynamics, ??<br />
secular equation, ??, ??<br />
shear modulus, ??<br />
– strain, ??<br />
shock, ??<br />
simple linear connected domain, ??<br />
– material, ??<br />
– pendulum, ??<br />
– shear, ??<br />
sinusoidal progressive wave, ??<br />
sliding velocity, 55, 75<br />
small oscillations, ??, ??<br />
sound speed, ??<br />
space of configurations, ??<br />
– of states, ??, ??, ??, ??<br />
– of minimal states, ??<br />
spatial cone, ??<br />
– coordinates, ??<br />
– gradient of deformation, ??<br />
speed, 6<br />
spherical polar coordinates, ??<br />
spontaneous motion, ??<br />
stable equilibrium position, ??<br />
– equilibrium state, ??<br />
state, ??, ??, ??<br />
– of motion, 49<br />
statically determined, ??<br />
– undetermined, ??<br />
stationary dynamic system, ??<br />
statistical mechanics, ??<br />
steady flow, ??<br />
Stokes flow, ??<br />
strain–displacement relation, ??<br />
strain energy, ??<br />
– energy density, ??<br />
stress power, ??<br />
strongly elliptic elasticity tensor, ??<br />
stuck point, ??<br />
successive harmonics, ??<br />
superpotential, ??<br />
surface traction, ??<br />
suspension axis, ??<br />
Symbolic Equation of Dynamics, ??<br />
– Equation of Statics, ??<br />
symmetric elasticity tensor, ??<br />
symmetry transformation, ??<br />
synchronous varied motions, ??<br />
T<br />
tangent unit vector, 8<br />
tangential acceleration, 11<br />
terminal velocity, ??<br />
theorem of axial moment, ??<br />
thermoelastic material, ??<br />
thermomechanic process, ??<br />
three–dimensional affine Euclidean space,<br />
2<br />
time law, 5<br />
torsional rigidity, ??<br />
trajectory, 5<br />
translation velocity field, 41<br />
translational kinetic state, 49<br />
– motion, 41<br />
– state of motion, 49<br />
transport acceleration, 61<br />
– force, ??<br />
– velocity, 61<br />
transverse acceleration, 14<br />
– velocity vector, 13<br />
– wave, ??<br />
two–body problem, ??
GLOSAR 93<br />
U<br />
uniform extension, ??<br />
– motion, 7<br />
– translation motion, 41<br />
uniformly varied motion, 21<br />
unilateral constraint, 31<br />
unit of force, ??<br />
– of mass, ??<br />
unstable equilibrium position, ??<br />
– equilibrium state, ??<br />
V<br />
Van der Pol equation, ??<br />
variation of a functional, ??<br />
vector, 2<br />
velocity of a point, 7<br />
vertical straight line, ??<br />
virtual displacement, ??<br />
– velocity, ??, ??<br />
– work, ??<br />
viscosity, ??<br />
viscous fluids, ??, ??<br />
Vlasov’s equation, ??<br />
Volterra–Lotka differential system, ??<br />
W<br />
warping function, ??<br />
Weierstrass method, ??<br />
weight force, ??<br />
work and kinetic energy theorem, ??<br />
– done by a system of forces, ??<br />
Y<br />
Young’s modulus, ??