Мехатронні системи: шпаргалка-підсумок

Аркуші з формулами, ключові кроки методів для кожної теми та 12 задач для практики, з повними розв’язками наприкінці.

Нові задачіНазад до курсу

У вікні друку виберіть «Зберегти як PDF», щоб зберегти копію на пристрої.

1Формули та методи

Двигуни постійного струму

Коло якоря
V = E + IaRa
Зворотна ЕРС і момент
E = k·ω, T = k·Ia
Швидкість
ω = 2π·n/60, ω = (V − IaRa)/k
Швидкість неробочого ходу й стопорний момент
ω0 = V/k, Tstall = k·V/Ra
Потужності
Pin = V·Ia, Pdev = E·Ia = T·ω, Pout = Pdev − Prot
Пуск
Istart = V/(Ra + Rst)
Регулювання швидкості
SR = (n0 − nFL)/nFL
Максимальна розвинена потужність
Pmax = V²/(4Ra)
Привід із чопером
Va = D·Vdc

Метод

  1. Накресліть коло якоря: джерело V, опір Ra і зворотна ЕРС E послідовно.
  2. Скористайтеся V = E + Ia·Ra, щоб зв’язати напругу, струм і зворотну ЕРС.
  3. Використовуйте E = k·ω для швидкості й T = k·Ia для моменту, де ω = 2π·n/60.
  4. Для потужності та ККД почніть із V·Ia, відніміть втрати в міді Ia²·Ra, а потім втрати на обертання.

Зубчасті колеса й передачі

Зачеплення коліс
n1·N1 = n2·N2
Крутний момент
T2 = T1·(N2/N1)·η
Потужність
P = T·ω, ω = 2π·n/60
Багатоступенева передача
загальне передаточне відношення = добуток передаточних відношень ступенів

Збір даних і дискретизація

Критерій Найквіста
fs > 2·fmax, fN = fs/2
Частота аліасу
fa = |f − k·fs|, k = round(f/fs)
Розмір кроку
1 LSB = VFSR/2ᴺ
Похибка квантування (округлення)
±½ LSB
Код АЦП з відсіканням
D = ⌊Vin/LSB⌋
Кількість розрядів для роздільності r
2ᴺ ≥ span/r
Ідеальне SQNR, синус на повну шкалу
SQNR = 6,02·N + 1,76 dB
Швидкість передавання даних
канали × fs × N біт/с
АЦП з мультиплексуванням
fs = 1/(n·(tc + tset))
UART
Tb = 1/бод; 8N1 = 10 біт на байт; байт/с = бод/біт у кадрі
Дільник швидкості передавання (бод)
baud = fclk/(16·D)

Метод

  1. Випишіть, що дано: частоту дискретизації, частоти сигналу, кількість розрядів, вхідний діапазон, швидкість у бодах.
  2. У задачах на дискретизацію порівнюйте кожну частоту з fs/2; вища частота виглядає як fa = |f − k·fs| з найближчим кратним fs.
  3. У задачах на роздільність знайдіть 1 LSB = VFSR/2ᴺ або найменше N, для якого 2ᴺ ≥ діапазон/роздільність.
  4. У задачах на швидкість даних і лінії UART враховуйте все, що передається: канали, відліки за секунду, біти чи цілі байти на відлік, а також стартові, паритетні й стопові біти кожного кадру.

Міст Вітстона й тензорезистори

Напруги дільників
Va = VsR2/(R1 + R2), Vb = VsR4/(R3 + R4)
Виведення
Vab = Va − Vb
Баланс
R1/R2 = R3/R4
Реохордний міст
R/Rx = l/(100 − l)
Опір Тевенена
Rth = R1∥R2 + R3∥R4
Коефіцієнт тензочутливості
ΔR/R = GF·ε
Вихід мосту
чвертьміст Vs·GF·ε/4, напівміст Vs·GF·ε/2, повний міст Vs·GF·ε
Мала зміна, одне плече
Vab ≈ Vs·ΔR/(4R)
Напруження
σ = E·ε

Метод

  1. Розглядайте міст як два дільники напруги на одному джерелі живлення.
  2. Для рівноваги покладіть R₁/R₂ = R₃/R₄ і розв’яжіть відносно невідомого резистора.
  3. Для незбалансованого мосту знайдіть Va і Vb з дільників і відніміть; коли вимірювач навантажує вихід, скористайтеся еквівалентом Тевенена.
  4. Для тензорезисторів переведіть деформацію в ΔR/R = GF·ε і користуйтеся Vs·GF·ε/4, /2 або ×1 для чвертьмосту, напівмосту чи повного мосту.

Усталена похибка та ПІД-регулювання

Тип системи
кількість полюсів G(s) у точці s = 0
Константи похибок
Kp = lim G(s), Kv = lim s·G(s), Ka = lim s²·G(s) (s → 0)
Східчастий сигнал A
ess = A/(1 + Kp)
Лінійно наростальний сигнал A·t
ess = A/Kv
Параболічний сигнал A·t²/2
ess = A/Ka
ПІД-регулятор
C(s) = Kp + Ki/s + Kd·s
Цільова система другого порядку
s² + 2ζωns + ωn²
Перерегулювання та часові показники
%OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd, ts ≈ 4/(ζωn)

Метод

  1. Підрахуйте полюси розімкненого контуру C(s)·G(s) у точці s = 0: це тип системи.
  2. Знайдіть константу похибки, що відповідає входу: Kp для східчастого, Kv для лінійно наростального, Ka для параболічного сигналу.
  3. Застосуйте ess = A/(1 + Kp), A/Kv або A/Ka, або розв’яжіть це відносно потрібного коефіцієнта підсилення.
  4. Для вимоги щодо демпфування чи швидкодії запишіть характеристичний поліном замкненої системи й порівняйте його з s² + 2ζω_n·s + ω_n².

2Задачі для практики

1.Двигуни постійного струмуСередня
T (N·m)n (rpm)0,3 N·mn0​n = ?

Двигун постійного струму (k = 0,15 N·m/A, Ra = 2 Ω) працює від 48 V і приводить навантаження з моментом 0,3 N·m. Знайдіть його швидкість n у rpm.

V = 48 Vk = 0,15 N·m/ARa = 2 ΩT = 0,3 N·m
n = rpm
2.Двигуни постійного струмуСкладна
+−48 VRs​ = ?Ra​ = 1,5 ΩM+−En = 1000 rpmT = 0,4 N·m

Двигун постійного струму (k = 0,2 V·s/rad, Ra = 1,5 Ω) на сталій напрузі 48 V приводить навантаження 0,4 N·m. Знайдіть послідовний опір Rs, який зменшує його швидкість до 1000 rpm.

V = 48 Vk = 0,2 V·s/radRa = 1,5 ΩT = 0,4 N·mn = 1000 rpm
Rs = Ω
3.Зубчасті колеса й передачіСередня
N1​ = 40N2​ = 50

Двигун потужністю 7,5 kW із частотою обертання 1450 rpm приводить шестірню з кількістю зубців 40, яка перебуває в зачепленні з колесом з кількістю зубців 50. ККД зубчастої пари — 96 %. Знайдіть вихідний крутний момент T2 на валу колеса.

P = 7,5 kWn1 = 1450 rpmN1 = 40N2 = 50η = 96%
T2 = N·m
4.Зубчасті колеса й передачіСкладна
d1​ = 120 mmd2​ = 315 mm1500 rpmC = 1000 mmT1​ = 800 NT2​P = ?

Відкритий пас з’єднує ведучий шків 120 мм при 1500 об/хв зі шківом 315 мм; центри розташовані на відстані 1000 мм. Коефіцієнт тертя 0,35, а натяг ведучої гілки не повинен перевищувати 800 Н. Знайдіть найбільшу потужність P, яку пас може передати без прослизання, у кВт (знехтуйте відцентровими ефектами).

d1 = 120 mmn1 = 1500 rpmd2 = 315 mmC = 1000 mmμ = 0,35T1 = 800 N
P = kW
5.Збір даних і дискретизаціяСередня
Vin​N = ?200 mm2ᴺ − 10≤ 0,1 mm

Датчик вимірює відстань у діапазоні 200 mm, а його вихід заповнює весь вхідний діапазон АЦП. Показ має розрізняти 0,1 mm. Знайдіть найменшу кількість розрядів N, потрібну АЦП.

діапазон = 200 mmроздільність = 0,1 mm
N =
6.Збір даних і дискретизаціяСкладна
T10 mV/°C×4∿01N = 10, 3,3 V? °C/LSB

Датчик температури дає 10 мВ/°C (0 В при 0 °C). Його сигнал підсилюється в 4 раз і подається на 10-розрядний АЦП із діапазоном від 0 до 3,3 В (1 LSB = VFSR/2ᴺ). Знайдіть зміну температури ΔTLSB, що відповідає одному LSB.

S = 10 mV/°CG = 4N = 10VFSR = 3,3 V
ΔTLSB = °C
7.Міст Вітстона й тензорезисториСередня
+−15 V100 Ω100 Ω + ΔR100 Ω100 Ωab+−Vab​ = ?

Міст має чотири плеча по 100 Ω і Vs = 15 V. Датчик R₂ збільшується на ΔR = 0,25 Ω. У наближенні малих змін знайдіть вихід Vab у mV.

Vs = 15 VR = 100 ΩΔR = 0,25 Ω
Vab = mV
8.Міст Вітстона й тензорезисториСкладна
+−2,5 V120 Ω120 Ω + ΔR120 Ω120 Ωab+−1,06 mV

Чвертьміст з одним активним тензорезистором (GF = 2,12) живиться від 2,5 V. Під навантаженням вихід показує 1,06 mV. Знайдіть деформацію ε у µε.

Vs = 2,5 VGF = 2,12Vo = 1,06 mV
ε = µε
9.Усталена похибка та ПІД-регулюванняСередня
r(t) = 1·t²/2+20(s + 1)s²(s + 1)C(s)−

Параболічний вхід r(t) = 1·t²/2 подано на контур з одиничним зворотним зв’язком і G(s) = 20(s + 1)/(s²(s + 1)). Знайдіть усталену похибку ess.

G(s) = 20(s + 1)/(s²(s + 1))r(t) = 1·t²/2
ess =
10.Усталена похибка та ПІД-регулюванняСкладна
r(t) = 2t+5 + 1/s2(s + 2)(s + 4)C(s)−

ПІ-регулятор C(s) = 5 + 1/s керує об’єктом G(s) = 2/((s + 2)(s + 4)) у контурі з одиничним зворотним зв’язком. Знайдіть усталену похибку ess для лінійно наростального входу r(t) = 2t.

C(s) = 5 + 1/sG(s) = 2/((s + 2)(s + 4))r(t) = 2t
ess =
11.Перехідна характеристикаСередня
t (s)y(t)0,5?

Система має передатну функцію G(s) = 6,5/(s² + 4s + 13) і отримує на вході одиничну сходинку. Знайдіть відсоток перерегулювання %OS перехідної характеристики.

G(s) = 6,5/(s² + 4s + 13)
%OS = %
12.Перехідна характеристикаСкладна
+−5 VR = 57 ΩL = 50 mHC = 10 µFvC​

Послідовне RLC-коло з R = 57 Ω, L = 50 mH і C = 10 µF підключають до сходинки 5 V у момент t = 0. Конденсатор незаряджений, і до сходинки струм не тече. Знайдіть момент tp першого піку напруги в ms.

R = 57 ΩL = 50 mHC = 10 µFV = 5 V
tp = ms

3Розчини

1.Двигуни постійного струму

Відповідь: n = 2801 rpm

Двигун постійного струму: якір — це опір Ra, увімкнений послідовно із зворотною ЕРС E, тож V = E + Ia·Ra. Зворотна ЕРС зростає зі швидкістю, E = k·ω, а момент — зі струмом, T = k·Ia, з однаковою константою двигуна k (V·s/rad = N·m/A) в одиницях SI. Швидкості переводьте за формулою ω = 2π·n/60. Розвинена потужність дорівнює E·Ia = T·ω; споживана потужність — V·Ia, а різницю становлять втрати в міді Ia²Ra та втрати на обертання. Швидкість спадає за прямою лінією зі збільшенням моменту, ω = V/k − Ra·T/k², від швидкості неробочого ходу V/k до нуля при стопорному моменті k·V/Ra. У стані спокою E = 0, тож пусковий струм V/Ra великий, якщо його не обмежує резистор.

Струм якоря

Ia = T/k = 0,3/0,15 = 2 A

Зворотна ЕРС і швидкість

E = 48 − 2·2 = 44 V
ω = E/k = 293,3 rad/s → n = 2801 rpm

Відповідь

n = 2801 rpm
2.Двигуни постійного струму

Відповідь: Rs = 12,03 Ω

Двигун постійного струму: якір — це опір Ra, увімкнений послідовно із зворотною ЕРС E, тож V = E + Ia·Ra. Зворотна ЕРС зростає зі швидкістю, E = k·ω, а момент — зі струмом, T = k·Ia, з однаковою константою двигуна k (V·s/rad = N·m/A) в одиницях SI. Швидкості переводьте за формулою ω = 2π·n/60. Розвинена потужність дорівнює E·Ia = T·ω; споживана потужність — V·Ia, а різницю становлять втрати в міді Ia²Ra та втрати на обертання. Швидкість спадає за прямою лінією зі збільшенням моменту, ω = V/k − Ra·T/k², від швидкості неробочого ходу V/k до нуля при стопорному моменті k·V/Ra. У стані спокою E = 0, тож пусковий струм V/Ra великий, якщо його не обмежує резистор.

Струм і зворотна ЕРС

Ia = 0,4/0,2 = 2 A
E = 0,2·104,7 = 20,94 V

Послідовний опір

Ra + Rs = (48 − 20,94)/2 = 13,53 Ω
Rs = 12,03 Ω

Відповідь

Rs = 12,03 Ω

Резистор розсіює 48,11 W; чопер керує швидкістю з набагато меншими втратами.

3.Зубчасті колеса й передачі

Відповідь: T2 = 59,27 N·m

Зубчасті передачі: у колесах, що зачеплені, швидкість у точці контакту зубців однакова, тож n₁·N₁ = n₂·N₂: колесо з меншою кількістю зубців обертається швидше. Колеса на одному валу обертаються з однаковою частотою. Потужність передається далі (за вирахуванням втрат), тож зі зменшенням частоти обертання крутний момент зростає: T₂ = T₁·(N₂/N₁)·η, а P = T·ω, де ω = 2π·n/60.

Момент двигуна

ω₁ = 2π·1450/60 = 151,8 rad/s
T1 = P/ω₁ = 7500/151,8 = 49,39 N·m

Вихідний крутний момент

T2 = 49,39·50/40·0,96 = 59,27 N·m

Відповідь

T2 = 59,27 N·m
4.Зубчасті колеса й передачі

Відповідь: P = 4,851 kW

Пасові передачі: без прослизання обидва шківи мають однакову колову швидкість обода v = π·d·n/60, тож n₁·d₁ = n₂·d₂, а передаточне відношення дорівнює d₂/d₁. Прослизання s робить ведений шків трохи повільнішим: n₂ = n₁·(d₁/d₂)·(1 − s). Передавана потужність P = (T₁ − T₂)·v, де T₁ — натяг ведучої гілки, а T₂ — веденої. Тертя обмежує відношення натягів: T₁/T₂ = eμθ, де θ — кут обхвату в радіанах на меншому шківі; для відкритого паса θ = π − 2·asin((d₂ − d₁)/(2C)) при міжосьовій відстані C.

Кут обхвату

θ = 2,946 rad

Натяг веденої гілки

eμθ = e1,031 = 2,804, T₂ = 800/2,804 = 285,3 N

Потужність

P = (800 − 285,3)·9,425 = 4,851 kW

Відповідь

P = 4,851 kW
5.Збір даних і дискретизація

Відповідь: N = 11

Роздільність АЦП: N-розрядний перетворювач ділить свій діапазон повної шкали VFSR на 2ᴺ кодів, від 0 до 2ᴺ − 1. Один крок дорівнює 1 LSB = VFSR/2ᴺ: найменша зміна, яку може показати перетворювач. АЦП з відсіканням дає код ⌊Vin/LSB⌋; АЦП із округленням ніколи не відхиляється від істинного значення більш ніж на ½ LSB, це його похибка квантування. Кожен додатковий розряд зменшує LSB удвічі й додає близько 6 дБ до ідеального відношення сигнал/шум квантування SQNR = 6,02·N + 1,76 дБ для синуса на повну шкалу. Щоб досягти роздільності r на діапазоні, оберіть найменше N, для якого 2ᴺ ≥ діапазон/r.

Потрібна кількість кроків

200/0,1 = 2000

Розряди

log₂(2000) = 10,97
210 = 1024 < 2000 ≤ 211 = 2048
N = 11

Відповідь

N = 11
6.Збір даних і дискретизація

Відповідь: ΔTLSB = 0,08057 °C

Роздільність АЦП: N-розрядний перетворювач ділить свій діапазон повної шкали VFSR на 2ᴺ кодів, від 0 до 2ᴺ − 1. Один крок дорівнює 1 LSB = VFSR/2ᴺ: найменша зміна, яку може показати перетворювач. АЦП з відсіканням дає код ⌊Vin/LSB⌋; АЦП із округленням ніколи не відхиляється від істинного значення більш ніж на ½ LSB, це його похибка квантування. Кожен додатковий розряд зменшує LSB удвічі й додає близько 6 дБ до ідеального відношення сигнал/шум квантування SQNR = 6,02·N + 1,76 дБ для синуса на повну шкалу. Щоб досягти роздільності r на діапазоні, оберіть найменше N, для якого 2ᴺ ≥ діапазон/r.

Розмір кроку

1 LSB = 3300/1024 = 3,223 mV

Чутливість на АЦП

10·4 = 40 mV/°C

Температурна роздільність

ΔTLSB = 3,223/40 = 0,08057 °C

Відповідь

ΔTLSB = 0,08057 °C

АЦП досягає повної шкали при 82,5 °C.

7.Міст Вітстона й тензорезистори

Відповідь: Vab = 9,375 mV

Міст Вітстона: два дільники напруги живляться від одного джерела Vs. Лівий дільник R1–R2 дає Va = Vs·R2/(R1 + R2), а правий дільник R3–R4 дає Vb = Vs·R4/(R3 + R4); вихід — Vab = Va − Vb. Міст збалансований (Vab = 0, у вимірювачі немає струму), коли R1/R2 = R3/R4, тобто R1·R4 = R2·R3. З боку виходу міст — це джерело Тевенена з Rth = R1∥R2 + R3∥R4.

Відносна зміна

ΔR/R = 0,25/100 = 0,0025

Лінеаризований вихід

Vab ≈ 15·0,0025/4 = 9,375 mV

Відповідь

Vab = 9,375 mV

Точне значення становить 9,363 mV; наближення добре, поки ΔR/R мале.

8.Міст Вітстона й тензорезистори

Відповідь: ε = 800 µε

Тензорезистори: тензорезистор змінює опір під дією деформації: ΔR/R = GF·ε, де коефіцієнт тензочутливості GF для металевої фольги близько 2. У мосту, всі плечі якого дорівнюють R, один активний тензорезистор (чвертьміст) дає Vo ≈ Vs·GF·ε/4. Два активні тензорезистори з протилежними деформаціями в сусідніх плечах (напівміст) дають Vs·GF·ε/2, а чотири (повний міст) — Vs·GF·ε. Напруження випливає із закону Гука, σ = E·ε. Тензометричні датчики сили вказують чутливість у mV/V: вихід при повному навантаженні на один вольт живлення.

Розв’яжіть рівняння чвертьмосту

ε = 4·0,00106/(2,5·2,12) = 0,0008

У мікродеформаціях

ε = 800 µε

Відповідь

ε = 800 µε
9.Усталена похибка та ПІД-регулювання

Відповідь: ess = 0,05

Усталена похибка: у контурі з одиничним зворотним зв’язком і передавальною функцією розімкненої системи G(s) похибка дорівнює E(s) = R(s)/(1 + G(s)), а за теоремою про кінцеве значення ess = lims→0 s·E(s). Тип системи — це кількість полюсів G(s) у точці s = 0. Константи похибок: Kp = lim G(s), Kv = lim s·G(s) і Ka = lim s²·G(s), усі при s → 0. Східчастий сигнал величини A дає ess = A/(1 + Kp), лінійно наростальний A·t дає A/Kv, а параболічний A·t²/2 дає A/Ka. Контур типу 0 відпрацьовує східчастий сигнал зі скінченною похибкою, контур типу 1 відпрацьовує його точно, а лінійно наростальний — зі скінченною похибкою.

Константа прискорення

Ka = 20·1/1 = 20

Похибка за параболічного сигналу

ess = 1/20 = 0,05

Відповідь

ess = 0,05
10.Усталена похибка та ПІД-регулювання

Відповідь: ess = 8

П-, ПІ- і ПІД-регулювання: регулятор C(s) = Kp + Ki/s + Kd·s діє на похибку. Пропорційний коефіцієнт підсилення робить контур швидшим і зменшує похибку об’єкта типу 0, але ніколи не усуває її. Інтегральна складова додає полюс у точці s = 0: тип системи зростає на одиницю, і східчастий сигнал відпрацьовується без похибки. Диференціальна складова додає демпфування. Щоб підібрати коефіцієнти, запишіть характеристичний поліном замкненої системи й порівняйте його з s² + 2ζωns + ωn². Тоді %OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd при ωd = ωn√(1 − ζ²), а ts ≈ 4/(ζωn).

Швидкісна константа

s·C(s) = 5s + 1 → 1 при s → 0.

Kv = 1·2/8 = 0,25

Похибка за лінійно наростального сигналу

ess = 2/0,25 = 8

Відповідь

ess = 8

ПІ-регулятор повністю усуває похибку за східчастого сигналу; лінійно наростальний сигнал відпрацьовується з незмінним відставанням.

11.Перехідна характеристика

Відповідь: %OS = 12,31 %

Перехідна характеристика: система першого порядку G(s) = K/(τs + 1) відповідає на сходинку розміром U як y(t) = K·U·(1 − e−t/τ): вона досягає 63,2 % кінцевого значення K·U за одну сталу часу τ і 98 % за 4τ. Система другого порядку G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) має полюси s = −ζωₙ ± jωd, де ωd = ωₙ√(1 − ζ²). За 0 < ζ < 1 перехідна характеристика має перерегулювання: відсоток перерегулювання %OS = 100·e−ζπ/√(1−ζ²), час піку tp = π/ωd, а 2-відсотковий час встановлення ts ≈ 4/(ζωₙ). Порівняйте знаменник із s² + 2ζωₙs + ωₙ², щоб знайти ωₙ і ζ. Для системи маса–пружина–демпфер ωₙ = √(k/m) і ζ = c/(2√(km)); для послідовного RLC-кола ωₙ = 1/√(LC) і ζ = (R/2)·√(C/L).

Порівняйте зі стандартною формою

Знаменник має збігатися з s² + 2ζωₙs + ωₙ².

ωₙ² = 13 → ωₙ = 3,606 rad/s
2ζωₙ = 4 → ζ = 4/(2·3,606) = 0,5547

Результат

%OS = 100·e−ζπ/√(1−ζ²) = 100·e−2,094 = 12,31 %

Відповідь

%OS = 12,31 %
12.Перехідна характеристика

Відповідь: tp = 2,427 ms

Перехідна характеристика: система першого порядку G(s) = K/(τs + 1) відповідає на сходинку розміром U як y(t) = K·U·(1 − e−t/τ): вона досягає 63,2 % кінцевого значення K·U за одну сталу часу τ і 98 % за 4τ. Система другого порядку G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) має полюси s = −ζωₙ ± jωd, де ωd = ωₙ√(1 − ζ²). За 0 < ζ < 1 перехідна характеристика має перерегулювання: відсоток перерегулювання %OS = 100·e−ζπ/√(1−ζ²), час піку tp = π/ωd, а 2-відсотковий час встановлення ts ≈ 4/(ζωₙ). Порівняйте знаменник із s² + 2ζωₙs + ωₙ², щоб знайти ωₙ і ζ. Для системи маса–пружина–демпфер ωₙ = √(k/m) і ζ = c/(2√(km)); для послідовного RLC-кола ωₙ = 1/√(LC) і ζ = (R/2)·√(C/L).

Власна частота й коефіцієнт демпфування

ωₙ = 1/√(LC) = 1/√(0,05·1 × 10⁻⁵) = 1414 rad/s
ζ = (R/2)·√(C/L) = (57/2)·√(1 × 10⁻⁵/0,05) = 0,4031

Результат

ωd = 1294 rad/s
tp = π/ωd = 2,427 ms

Відповідь

tp = 2,427 ms

ExamProblems · examproblems.com/courses/mechatronic-systems