Znalezienie funkcji przenoszenia sprężynowego układu tłumika drgań

9

Przeglądałem książkę Ogata Modern Control Engineering i pracowałem nad kilkoma ćwiczeniami, aby lepiej zrozumieć podstawowe zasady kontroli. Natknąłem się na następujący przykład, który staram się rozwiązać.

Muszę wymyślić funkcję przenoszenia, która modeluje ten przyrząd wibracyjny. Pytania są następujące:

W tym przykładzie będziesz analizować urządzenie testujące wibracje (ryc. 1). System ten składa się z tabeli masy M i cewki o masie m. Magnes stały sztywno przymocowany do podłoża zapewnia stabilne pole magnetyczne. Ruch cewki 𝑦 przez pole magnetyczne indukuje napięcie w cewce proporcjonalne do jego prędkości its, jak w równaniu. 1. 𝑒 = 𝛼𝑦̇ [równ. 1]

Przepływ prądu przez cewkę powoduje, że odczuwa ona siłę magnetyczną proporcjonalną do prądu jak w równaniu. 2. 𝐹 = 𝛽𝑖 [równ. 2]

wprowadź opis zdjęcia tutaj

Pytanie: Uzyskaj parametryczną funkcję przesyłania z wyjściem 𝑥 na wejście 𝑉.

Na niektóre pytania, na które trudno mi odpowiedzieć, ale dotyczą całego TF, są:

  • Jeśli K2 i B2 są ściskane na odległość Z, (gdy porusza się w górę z
    powodu interakcji cewki z polem magnetycznym), czy to oznacza, że ​​k1 i b1 są przedłużone o tę samą odległość Z?

  • Jeśli m(cewka) przesuwa się w górę o 2 cm, to czy M(stół) również przesuwa się w górę o 2 cm?


Co muszę zrobić:

  • Wymyśl dwa oddzielne schematy wolnych korpusów, jeden dla masy M stołu i jeden dla masy m cewki.
  • Naszkicuj jeden schemat obwodu, w tym tylny emf.
  • Przekształć w domenę s.
  • Rozwiązuj jednocześnie.

Co do tej pory zrobiłem:

  • Narysuj, aby oddzielić schematy dowolnych ciał i wyodrębnić równania.

  • Narysuj schemat połączeń i wyodrębnij równanie.

  • Konwertuj na domenę s.

Za pomocą funkcji MATLAB solveudało mi się uzyskać 2 różne funkcje przenoszenia 5. rzędu (po jednej dla każdej metody, którą proponuję poniżej), jednak nie jestem pewien, która z nich jest poprawna i dlaczego.


Ogólny system :

Jest to schematyczna reprezentacja tego, jak według mnie można modelować przyrząd do testowania wibracji, z wyłączeniem części elektrycznej.

wprowadź opis zdjęcia tutaj


Schemat swobodnego ciała 1 - Tabela - Konwencja w górę

Sprężyny k1i k2i amortyzatory b1i b2modelowane oddzielnie . Ponieważ nie można ich dodawać i traktować jako jednego, ich kompresja i rozszerzenie są osobne.

W górę siła pochodzi k2i b2które są dołączone do cewki. Te doświadczają ruchu w górę.

wprowadź opis zdjęcia tutaj

Równanie w domenie s:

Ms^2X + b1sX + k1X = b2s(X-Y) + k2(X-Y)


Schemat swobodnego ciała 2 - cewka - konwencja zwyżkowa

Cewka odczuwa siłę skierowaną do góry, jednak sprężyna i tłumik powstrzymują ją, działając w przeciwnym kierunku.

wprowadź opis zdjęcia tutaj

Równanie w domenie s:

Fem = Ms^2Y + b2s(X-Y) + k2(X-Y)


Dwie różne metody przedstawione powyżej dla FBD tabeli prowadzą do różnych równań w dziedzinie s i różnych funkcji przenoszenia.

Jaki jest prawidłowy schemat swobodnego ciała dla stołu i cewki?

rrz0
źródło
2
Ładne pytanie, ale proszę zamieścić zdjęcie, w którym szczegóły są wyraźne, bez zmuszania nas do kliknięcia, aby je powiększyć. Np. Te znaki minus są ledwo dostrzegalne. Ponadto równanie w lewym dolnym rogu zostało częściowo przycięte. Na twoim arkuszu jest dużo wolnego miejsca, które można wykorzystać do powiększenia. W Internecie znajduje się wiele bezpłatnych programów do edycji obrazów (np. IrfanView lub FastSstone ImageViewer), dzięki czemu możesz także zrobić kilka zdjęć swojego arkusza (arkuszy) i wyciąć / przyciąć potrzebne części, aby opublikować ładne zdjęcia.
Lorenzo Donati - Codidact.org
@LorenzoDonati, dzięki za sugestię, natychmiast edytuje. Jeśli chodzi o równanie w lewym dolnym rogu, nie jest to interesujące, ponieważ moim zmartwieniem jest schemat swobodnego ciała. Jeśli jest to poprawne, równanie będzie poprawne. Spróbuję jednak odpowiednio edytować. Dziękuję za twój komentarz.
rrz0
Staraj się nie zakładać, co zrobiłeś źle. Opublikowanie zestawu ładnie narysowanych równań zgodnie z tokiem myślenia pokaże twoje wysiłki (a tym samym poprawi twoje pytanie - dając mu więcej szans na odpowiedź), a także może wskazać możliwe błędy. Wszelkie istotne informacje dotyczące twojego problemu mogą być przydatne dla potencjalnego użytkownika odpowiadającego.
Lorenzo Donati - Codidact.org
BTW, jeśli jesteś zadowolony ze składni LaTeX, edytor pytań może zrozumieć „notację dolara” formuł LaTeX (patrz pomoc online).
Lorenzo Donati - Codidact.org
1
Dzięki @LorenzoDonati, staram się przedstawić pytanie w bardziej uporządkowany i czytelny sposób.
rrz0

Odpowiedzi:

2

Wprowadzenie

M i m mają tylko jeden stopień swobody; oba mogą poruszać się tylko pionowo. Siła magnetyczna działa bezpośrednio na magnes m, a nie na masę M.

Aby trochę zdemistyfikować obraz, warto pomyśleć o magnesie ustawionym po drugiej stronie stołu. Obraz został narysowany w LTSPICE i nie ma żadnych strzałek. Więc najbliższe przybliżenie do strzałki to pin wyjściowy, a ponieważ mogą one wskazywać tylko poziomo w prawo, cały obraz jest obracany o w prawo. Z tego samego powodu strzałki „-y” i „-F” wskazują w prawo, podczas gdy wolałbym narysować strzałki „y” i „F” w lewo. Ponadto prawy powinien brzmieć .90ob1b2

Teraz jest jasne, że jest to szeregowe połączenie mas z elementami dynamicznymi między nimi, więc zaczynamy zapisywać równania ruchowe od prawej do lewej, zaczynając od równania elektrycznego dla m, które będzie zawierało najpierw V, y i F.
Następnie napiszemy równanie ruchowe dla m i dla M.
Ponieważ na M nie ma wpływu siła magnetyczna, to ostatnie równanie da nam y jako funkcję x, która zostanie użyta w pierwszym równaniu w odniesieniu do x V.

Elektryczny

Siła magnetyczna i ruch magnesu są sprzężone przez napięcie na cewce. A ponieważ przy założeniu, że L jest niezależne od y, mamy

e=αy˙,F=βi,Ve=Ri+Li˙
Ve=Vαy˙=Ri+Li˙=RβF+LβF˙

Teraz mamy pod względem (i ) i możemy zapisać równania ruchu, dodając wszystkie siły do ​​poruszających się obiektów i zmuszając je do zera (zgodnie z prawem).yFV

Magnes

F+my¨+b2(y˙x˙)+k2(yx)=0
Vαy˙=V(s)αsy=(R+Ls)i=(R+Ls)F/β
F=βR+Ls(V(s)αsy)
βV(s)R+LsαβR+Lssy+my¨+k2(yx)+b2(y˙x˙)=0
βV(s)R+LsαβR+Lssy+ms2y+k2(yx)+b2s(yx)=0
ms2y+(b2αβR+Ls)sy+k2yb2sxk2x=βV(s)R+Ls
xy
(ms2+b2sαβsR+Ls+k2)y(b2s+k2)x=βV(s)R+Ls

Ruchomy stół

Mx¨+k1x+b1x˙+k2(xy)+b2(x˙y˙)=0
Ms2x+k1x+b1sx+k2(xy)+b2s(xy)=0
b2syk2y+Ms2x+(b1+b2)sx+(k1+k2)x=0
xy
(b2s+k2)y+{Ms2+(b1+b2)s+k1+k2}x=0
y=Ms2+(b1+b2)s+k1+k2b2s+k2x

Ensemble

y=f(x)xyV

[(ms2+b2sαβsR+Ls+k2)Ms2+(b1+b2)s+k1+k2b2s+k2(b2s+k2)]x=βV(s)R+Ls

R+Ls

[{(R+Ls)(ms2+bs+k2)αβs}Ms2+(b1+b2)s+(k1+k2)b2s+k2(R+Ls)(b2s+k2)]x=βV(s)

b2s+k2

[{(R+Ls)(ms2+bs+k2)αβs}{Ms2+(b1+b2)s+(k1+k2)}(R+Ls)(b2s+k2)2]x=(b2s+k2)βV(s)

x(s)/V(s)

Joe Electro
źródło
Komentarze nie są przeznaczone do rozszerzonej dyskusji; ta rozmowa została przeniesiona do czatu .
Dave Tweed,