Łańcuch kinematyczny

Z testwiki
Wersja z dnia 14:52, 3 maj 2024 autorstwa imported>Froscior (growthexperiments-addlink-summary-summary:2|1|0)
(różn.) ← poprzednia wersja | przejdź do aktualnej wersji (różn.) | następna wersja → (różn.)
Przejdź do nawigacji Przejdź do wyszukiwania

Łańcuch kinematyczny – część mechanizmu w postaci kilku połączonych ze sobą członów tworzących jedną lub wiele par kinematycznych[1], realizujący zdefiniowane przeniesienie ruchu.

Łańcuchy kinematyczne dzielą się na:

  • kinematyczne płaskie,
  • kinematyczne przestrzenne.

Własności ogólne

Podstawową cechą łańcucha kinematycznego jest jego ruchliwość. Ruchliwość określa ile stopni swobody posiada łańcuch, to znaczy ile różnych typów ruchu jest w stanie przenieść.

Ruchliwość w może być:

w=0 lub w<0 – łańcuch sztywny
w=1 – łańcuch normalny
w>1 – łańcuch swobodny

Ruchliwość łańcucha oblicza się ze wzoru strukturalnego (kryterium Czebszowa–Grüblera–Kutzbacha):

  • dla łańcucha przestrzennego:
w=6(n1)p12p23p34p45p5,
uwaga: od liczby ogniw n odejmuje się 1 w przypadku, gdy nie zalicza się do niej nieruchomej podstawy.
  • dla łańcucha płaskiego:
w=3np42p5,
gdzie:
n – liczba ogniw,
pn – ilość par kinematycznych n-tej klasy.

Interpretacja ruchliwości obliczonej ze wzoru strukturalnego wymaga pewnego doświadczenia, w szczególności gdy wskazuje on, iż łańcuch kinematyczny jest sztywny. Taka sytuacja jest oczywista w przypadku a. W pewnych przypadkach jednak, przy szczególnej geometrii, łańcuch teoretyczne sztywny może przenosić ruch (przypadek b). Strukturalnie jest on identyczny z a, lecz występują w nim trzy geometrycznie identyczne człony co umożliwia ruch łańcucha. Para kinematyczna, która powinna łańcuch usztywniać (wskazana czerwoną strzałką) jest węzłem biernym. Projektując mechanizm z węzłami biernymi, konstruktor musi zdawać sobie sprawę, że zużycie elementów mechanizmu, prowadzące do drobnych zmian w ich geometrii, może doprowadzić do usztywnienia mechanizmu. Łańcuch kinematyczny sztywny, aczkolwiek teoretycznie takim jest w praktyce jest stosowany jako konstrukcja służąca wyłącznie do przenoszenia obciążenia, a nie ruchu i nie jest przedmiotem teorii mechanizmów i maszyn, lecz wytrzymałości materiałów i teorii konstrukcji.

Łańcuch kinematyczny o ruchliwości równej jeden jest najczęstszym przypadkiem mechanizmu. Rodzaj ruchu członu czynnego determinuje wtedy ruch członu biernego i wszystkich członów pośredniczących. Rysunek c pokazuje typowy czworobok przegubowy o ruchliwości równej jeden.

W pewnych przypadkach wymagane jest by mechanizm miał większą ruchliwość. Typowym tego przykładem jest przekładnia obiegowa d, lub mechanizm kreślarski e. Istnieją przypadki, że w bardzo odpowiedzialnych mechanizmach powiększa się ruchliwość normalnie zabezpieczoną bezpiecznikiem, by uniknąć samousztywnienia się mechanizmu w wyniku zużycia lub odkształcenia elementów. Gdy sytuacja taka zaistnieje bezpiecznik uruchamia dodatkową parę. Mechanizm może wtedy stracić swoją funkcjonalność, lecz unika się wtedy trwałego zniszczenia mechanizmu lub środowiska w jakim pracuje.

Analiza ruchliwości

Badanie ruchliwości łańcucha kinematycznego może być przeprowadzone analitycznie. Ograniczając analizę do przypadku łańcuchów płaskich (dwuwymiarowych) założymy, że jego ogniwa są prętami prostymi, odkształcalnymi osiowo i połączonymi ze sobą w węzłach za pomocą idealnych przegubów. Do rozważań wprowadzimy dwa wektory

𝐝=[d1,d2,,dm] – liniowo niezależnych przemieszczeń węzłowych i
𝐬=[s1,s2,,sn] – wydłużeń/skróceń ogniw spowodowanych przemieszczeniami węzłowymi.

W celu uproszczenia zapisu wektory będziemy zapisywali albo wierszowo albo kolumnowo jeżeli nie będzie to prowadziło do nieporozumień.

Analiza ruchliwości węzłów łańcucha polega na badaniu rozwiązań układu równań

(a)𝐀𝐝=𝐬.

Macierz zgodności przemieszczeń węzłowych 𝐀, występująca w tym równaniu, może być utworzona na podstawie interpretacji jej kolumn. Rozważmy mianowicie węzeł łańcucha, w którym występują dwa przemieszczenia di i dj (rys. 1a-b). Związek (a) przyjmuje dla tego przypadku postać

𝐀𝐝=𝐬[skiskjslisljsmismj][didj]=[skslsm].
Rys. 1a-b – warunki zgodności przemieszczeń węzłowych

Elementy ski,sli,smi są wydłużeniami/skróceniami odpowiednio k-tego, l-tego i m-tego elementu ogniwa łańcucha, spowodowanymi przez i-te przemieszczenie węzłowe di=1 (rys. 1a). Analogicznie elementy skj,slj,smj są skutkami spowodowanymi przez przemieszczenie węzłowe dj=1 (rys. 1b). Elementy sk,sl,sm są sumarycznymi skutkami równoczesnego działania przemieszczeń di i dj.

Jeżeli wydłużeniom ogniw przypiszemy znak plus, a skróceniom – minus, to na podstawie rys. 1a-b otrzymamy

ski=+cos(k,i),sli=cos(l,i),smi=cos(m,i), – (rys. 1a),
skj=+cos(k,j),slj=cos(l,j),smj=+cos(m,j), – (rys. 1b).

W rzeczywistym łańcuchu kinematycznym jego ogniwa mają niezmienną długość, a to znaczy, że 𝐬=𝟎. Badaniu ruchliwości łańcucha kinematycznego należy więc poddać nie układ równań (a), lecz (b)

(b)𝐀𝐝=𝟎.

Z algebry liniowej wiadomo, że liczba liniowo niezależnych niezerowych rozwiązań takiego układu równań jest równa defektowi badanej macierzy. Pod tym terminem rozumiemy różnicę pomiędzy stopniem, a rzędem macierzy. W przypadku macierzy prostokątnych o rozmiarach m×n stopniem nazwiemy mniejszą z liczb m i n.

Przykład 1

Dla przykładu rozważymy układ pokazany na rys. 2a[2].

Rys. 2a – łańcuch kinematyczny o ruchliwości w=1. Rys. 2b-f – szkice ilustrujące budowanie kolumn macierzy 𝐀.

Na podstawie rys. 2b-f kolejno budujemy kolumny macierzy 𝐀:

𝐀𝐝=[sin(α)cos(α)000cos(β)sin(β)cos(β)sin(β)000sin(γ)cos(γ)000cos(β)sin(β)1][d1d2d3d4d5]=0.

Łatwo można sprawdzić, że przyjęcie d1=1 pozwala wyznaczyć jednoznacznie wartości pozostałych niewiadomych przemieszczeń węzłowych, które wynoszą:

d2=tg(α), d3=cos(β)tg(α)sin(β)cos(β)+tg(γ)sin(β),
d4=d3tg(γ), d5=cos(β)tg(α)sin(β).

Dla łańcucha o pionowych „słupkach” (α=γ=0) otrzymuje się

d1=d3=1,d2=d4=0,d5=cos(β).

Dla łańcucha o „słupkach” równoległych (γ=α) jest

d1=d3=1,d2=tg(α),d5=cos(β)tg(α)sin(β).

Dla łańcucha z „ryglem” poziomym (β=0)

d1=d3=1,d2=tg(α),d4=tg(γ),d5=1.

Przykład 2

W przypadku łańcuchów wielokrotnie przesuwnych sprawa nieco się komplikuje co zilustrujemy na przykładzie łańcucha „piętrowego” z rys. 3aSzablon:R. Wygenerowana kolumnowo macierz zgodności przemieszczeń węzłowych 𝐀 ma postać

𝐀=[010000000001000010100000010001000001000100001010].
Rys. 3a – schemat łańcucha o ruchliwości w=2. Rys. 3b-c – dwa rodzaje przemieszczeń bazowych. Rys. d – nowy rodzaj przemieszczenia (przykładowa kombinacja liniowa).

Rozwiązując układ równań 𝐀𝐝=𝟎 za pomocą algorytmu Gaussa-Jordana z pełnym wyborem elementów wiodących stwierdzamy, że ma on rozwiązanie

d1=d3,d5=d7,d2=d4=d6=d8=0.

Oznacza to, że przemieszczenia d3 i d7 mogą spełniać rolę zmiennych swobodnych. Dzięki temu otrzymujemy w tym przypadku dwa rozwiązania podstawowe (bazowe)

d𝟏*=[10100000],d𝟐*=[00001010].

Dowolna kombinacja liniowa tych rozwiązań o postaci

d𝟏*=αd𝟏*+βd𝟐*,|α|+|β|0,

jest nowym rozwiązaniem równania 𝐀𝐝=𝟎.

Na rys. 3b-c pokazano mechaniczne interpretacje rozwiązań bazowych d𝟏* i d𝟐*, natomiast rys. 3d przedstawia przykładowe nowe rozwiązanie w postaci najprostszej kombinacji liniowej

d𝟑*=d𝟏*+d𝟐*=[10101010].

Przykład 3

Generowanie macierzy zgodności przemieszczeń węzłowych 𝐀 można przeprowadzić według takiego algorytmu postępowania, który bez żadnych zmian pozwala analizować ruchliwość zarówno łańcuchów płaskich, jak i przestrzennych. Podstawą tego algorytmu jest spostrzeżenie, że zmiana długości sij ogniwa i wywołana jednostkowym przemieszczeniem węzła j wyraża się prostym iloczynem skalarnym

(c)sij=±eidj=±cos(i,j),

w którym przez ei oznaczono wersor kierunku ogniwa, a przez dj wersor kierunku przemieszczenia węzłowego.

Znak we wzorze (c) przypiszemy według następującego kryterium:

(A) – gdy ogniwo i ma kierunek „do węzła” j, wtedy we wzorze (c) obowiązuje znak plus. Dla ogniwa skierowanego „od węzła” przypiszemy znak minus. To kryterium znakowania we wzorze (c) wymaga przyporządkowania każdemu ogniwu i konkretnego wersora kierunku ei.

W celu uproszczenia zapisu będziemy operować tylko tymi wektorami dwuwymiarowymi, które mają co najmniej jedną współrzędną niezerową (rys. 2a).

I tak dla łańcucha z przykładu 1 mamy

e1=[sin(α),cos(α)],e2=[cos(β),sin(β)],
e3=[sin(γ),cos(γ)],e4=[cos(β),sin(β)]

oraz

d1=[1,0],d2=[0,1],d3=[1,0],d4=[0,1],d5=[cos(β),sin(β)].

Na podstawie wzoru (c) otrzymujemy

𝐀=[+e1d1+e1d2000e2d1e2d2+e2d3+e2d4000e3d3e3d4000e4d3e4d4+e4d5].

Po obliczeniu iloczynów skalarnych możemy stwierdzić, że otrzymany wynik jest identyczny z uzyskanym w przykładzie 1.

Rys. 4a – przestrzenny łańcuch kinematyczny

Stosując to wektorowe podejście przeprowadzimy teraz analizę łańcucha przestrzennego z rys. 4a. Jedyna różnica będzie polegać tylko na tym, że wektory ei i dj będą teraz miały po trzy współrzędne.

Współrzędne wersorów dla poszczególnych ogniw można obliczyć na podstawie geometrii układu określonej wysokością OS=L i połową przekątnej kwadratowej podstawy OA=L.

e1=[1/21/22/2],e2=[1/21/22/2],
e3=[1/21/22/2],e4=[1/21/22/2],
d1=[100],d2=[010],d3=[001],d4=[100],
d5=[010],d6=[100],d7=[010].

Elementy macierzy 𝐀 obliczamy na podstawie wzoru (c) z przypisaniem znaków według kryterium (A):

(d)𝐀=[1/21/22/200001/21/22/2001/21/21/21/22/21/21/2001/21/22/20000].

Rozważymy teraz kilka różnych przypadków:

Przypadek 1

Załóżmy, że d4=d5=d6=d7=0. Oznacza to, że węzły B i C są unieruchomione i układ równań 𝐀𝐝=𝟎 ma w tym przypadku postać

𝐀=[1/21/22/21/21/22/21/21/22/21/21/22/2][d1d2d3]=𝟎.

Łatwo można sprawdzić, że istnieje taki minor 𝐌(3×3) macierzy 𝐀, że Det(𝐌)0. Wynika stąd wniosek, że jedynymi rozwiązaniami równania 𝐀𝐝=𝟎 są w tym przypadku d1=d2=d3=0. Oznacza to, że unieruchomienie węzłów B i C powoduje unieruchomienie całego łańcucha (brak ruchliwości).

Przypadek 2

Uwzględnimy teraz występowanie w węzłach B i C przemieszczeń tylko w kierunku osi 0x2 tzn. przyjmiemy, że d4=d6=0. Wynika stąd taka postać układu równań

(e)𝐀=[1/21/22/2001/21/22/201/21/21/22/21/201/21/22/200][d1d2d3d5d7]=𝟎.

Przyjmijmy, że zostało wywołane przemieszczenie d5=1. Przyjęcie to umożliwia jednoznaczne rozwiązanie układu równań (e)

𝐝=[d1d2d3d5d7]=[01/22/411],

co oznacza, że łańcuch ma jeden stopień swobody ruchu.

Przypadek 3

Przy tym samym założeniu co w przypadku 2, tzn. że d4=d6=0 wymusimy przemieszczenie d3=1. I tym razem istnieje jednoznaczne rozwiązanie

𝐝=[d1d2d3d5d7]=[0212222].

Obydwa rozwiązania z przypadków 2 i 3 różnią się tylko stałym mnożnikiem 22 co tylko potwierdza podobieństwo obu stanów przemieszczenia w układzie o jednym stopniu swobody ruchu.

Algorytm alternatywny

Przytoczymy teraz inny, statyczny sposób generowania macierzy 𝐀. Skorzystamy w tym celu z równań równowagi węzłów wyciętych z łańcucha, które można zapisać w postaci

(f)𝐁𝐒=𝐃,

przy czym

𝐒=[S1,S2,,Sn] – jest wektorem sił osiowych w ogniwach, a
𝐃=[D1,D2,,Dm] – wektorem sił węzłowych.
Rys. 5a-b – Statyczny wariant analizy ruchliwości łańcucha kinemtycznego

Dla przykładu rozważymy równanie równowagi sił działających na wycięty węzeł W (rys. 5a) w kierunku określonym przez jego wersor di=1. Równanie to otrzymamy sumując rzuty na ten kierunek wszystkich sił Sα,α(k,l,m).

Di=+(diek)Sk(diel)Sl(diem)Sm=sikSksilSlsimSm.

Podobnie mamy dla kierunku dj (rys. 5b)

Dj=+(djek)Sk(djel)Sl(djem)Sm=sjkSksjlSlsjmSm.

Dla łańcucha o n ogniwach i m liniowo niezależnych przemieszczeniach węzłowych macierz 𝐁 możemy zapisać w postaci

𝐁=[d1e1d1e2d1end2e1d2e2d2endme1dme2dmen],

która jednak nie uwzględnia znakowania określonego przez kryterium (A).

Analogicznie możemy napisać na podstawie wzoru (c)

𝐀=[e1d1e1d2e1dme2d1e2d2e2dmend1end2endm].

Wobec tego, że elementy aij oraz bij macierzy 𝐀 i 𝐁 spełniają związki

aij=eidj=djei=bji

otrzymujemy istotny związek

𝐁T=𝐀

pozwalający generować elementy macierzy 𝐀 na dwa równorzędne sposoby.

Zobacz też

Przypisy

Szablon:Przypisy

  1. A. Gawęcki, Podstawy mechaniki konstrukcji prętowych, Wyd. Politechniki Poznańskiej, Poznań 1985.
  2. B. Olszowski, M. Radwańska, Mechanika budowli, t. 1, s. 147, Kraków 2010, Wyd. Politechniki Krakowskiej.