Przygotowanie do egzaminu FE Mechanical: zestaw na powtórkę

Wzory, najważniejsze kroki metod dla każdego tematu i 42 zadań do ćwiczenia, z pełnymi rozwiązaniami na końcu.

Nowe zadaniaWróć do przedmiotu

W oknie drukowania wybierz „Zapisz jako PDF”, aby zachować kopię na swoim urządzeniu.

1Wzory i metody

Reakcje podpór belek

Równowaga
ΣFx = 0, ΣFy = 0, ΣM = 0
Obciążenie równomierne
wypadkowa w·L w środku obciążenia
Moment
M = F·d

Metoda

  1. Narysuj schemat sił i zastąp każdą podporę jej reakcjami.
  2. Zastąp każde obciążenie ciągłe jedną siłą przyłożoną w jego środku ciężkości.
  3. Napisz równanie momentów względem jednej podpory, aby wyznaczyć reakcję w drugiej.
  4. Z ΣFy = 0 wyznacz pozostałą reakcję i sprawdź drugim równaniem momentów.

Naprężenie i odkształcenie osiowe

Naprężenie normalne
σ = F/A
Odkształcenie
ε = ΔL/L
Prawo Hooke’a
σ = E·ε
Wydłużenie
ΔL = FL/(EA)
Odkształcenie termiczne
εT = α·ΔT

Skręcanie wałów

Naprężenie styczne
τ = T·r/J
Wał pełny
J = πd⁴/32
Wał drążony
J = π(D⁴ − d⁴)/32
Kąt skręcenia
φ = TL/(GJ)
Moc
P = T·ω = 2πnT/60

Metoda

  1. Wyznacz biegunowy moment bezwładności: J = πd⁴/32 dla wału pełnego.
  2. Naprężenie styczne na powierzchni: τ = T·r/J, gdzie r to promień zewnętrzny.
  3. Kąt skręcenia: φ = TL/(GJ) w radianach.
  4. Jeśli podana jest moc, najpierw wyznacz moment z P = T·ω, ω = 2πn/60.

Siła tnąca i moment zginający

Siła tnąca i moment
V = ΣF po jednej stronie, M = ΣF·d po jednej stronie
Naprężenie od zginania
σ = M·y/I
Prostokąt
I = bh³/12
Wskaźnik wytrzymałości przekroju
W = I/ymax
Swobodnie podparta, obciążenie równomierne
Mmax = wL²/8

Metoda

  1. Wyznacz reakcje podpór z równowagi całej belki.
  2. Przetnij belkę w potrzebnym przekroju i zostaw jedną stronę.
  3. Siła tnąca: suma sił pionowych po tej stronie. Moment zginający: suma ich momentów względem przekroju.
  4. Aby znaleźć moment maksymalny, wyznacz miejsce, w którym siła tnąca jest równa zeru.

Koło Mohra i transformacja naprężeń

Środek
C = (σx + σy)/2
Promień
R = √(((σx − σy)/2)² + τxy²)
Naprężenia główne
σ1,2 = C ± R
Kierunek główny
tan 2θp = 2τxy/(σx − σy)
Płaszczyzna ukośna
σθ = C + (σx − σy)/2·cos 2θ + τxy·sin 2θ

Metoda

  1. Oblicz środek koła Mohra: C = (σx + σy)/2.
  2. Oblicz promień: R = √(((σx − σy)/2)² + τxy²).
  3. Naprężenia główne to C ± R, a największe naprężenie styczne w płaszczyźnie to R.
  4. Kąty na kole są dwa razy większe niż kąty na elemencie: tan 2θp = 2τxy/(σx − σy).

Ugięcie belki

Siła skupiona w środku rozpiętości
δ = PL³/(48EI)
Obciążenie równomierne, belka swobodnie podparta
δ = 5wL⁴/(384EI)
Wspornik, siła na końcu
δ = PL³/(3EI)
Wspornik, obciążenie równomierne
δ = wL⁴/(8EI)
Superpozycja
δ = δ1 + δ2

Wyboczenie słupów

Siła krytyczna Eulera
Pcr = π²EI/(KL)²
Współczynnik długości wyboczeniowej
K = 1 przegub–przegub, 0,5 utwierdzenie–utwierdzenie, 0,7 utwierdzenie–przegub, 2 utwierdzenie–wolny koniec
Naprężenie krytyczne
σcr = Pcr/A
Smukłość
λ = KL/r, r = √(I/A)

Rzut ukośny

Składowe
vx = v0cos θ, vy = v0sin θ
Położenie
x = vxt, y = h + vyt − ½gt²
Maksymalna wysokość
H = h + vy²/(2g)
Zasięg (teren płaski)
R = v0²·sin 2θ/g

Metoda

  1. Rozłóż prędkość początkową na składowe: vx = v₀·cos θ i vy = v₀·sin θ.
  2. Ruch poziomy ma stałą prędkość; ruch pionowy ma przyspieszenie −g.
  3. Wyznacz czas lotu z ruchu pionowego.
  4. Użyj tego czasu w ruchu poziomym, aby wyznaczyć zasięg.

Praca i energia

Energia kinetyczna
Ek = ½mv²
Energia potencjalna
Ep = mgh
Bilans energii
Ek1 + Ep1 + W = Ek2 + Ep2 + losses
Praca tarcia
W = μ·N·d
Moc
P = W/t

Drgania mechaniczne

Częstotliwość drgań własnych
ωn = √(k/m), fn = ωn/(2π)
Z ugięcia statycznego
ωn = √(g/δst)
Tłumienie krytyczne
cc = 2√(km) = 2mωn
Współczynnik tłumienia
ζ = c/cc
Częstotliwość drgań tłumionych
ωd = ωn√(1 − ζ²)
Logarytmiczny dekrement tłumienia
δ = (1/n)·ln(x0/xn), ζ = δ/√(4π² + δ²)
Sprężyny połączone równolegle
k = k1 + k2
Sprężyny połączone szeregowo
1/k = 1/k1 + 1/k2

Silniki cieplne i chłodziarki

Sprawność
η = W/QH = 1 − QC/QH
Sprawność Carnota
ηC = 1 − TC/TH (kelvin)
COP chłodziarki
COPR = QC/W
COP pompy ciepła
COPHP = QH/W = COPR + 1
Bilans energii
QH = W + QC

Obiegi silników cieplnych gazowych

Izentropa, stosunek objętości
T2/T1 = rk−1
Izentropa, stosunek ciśnień
T2/T1 = rp(k−1)/k
Sprawność obiegu Otto
η = 1 − 1/rk−1
Sprawność obiegu Diesla
η = 1 − (1/rk−1)·(rck − 1)/(k(rc − 1))
Sprawność obiegu Braytona
η = 1 − 1/rp(k−1)/k
Doprowadzone ciepło
Otto: cv(T3 − T2) Diesel, Brayton: cp(T3 − T2)
Praca netto obiegu Braytona
w = cp[(T3 − T4) − (T2 − T1)]
Sprawności izentropowe
wc = wc,s/ηc, wt = ηt·wt,s
Model zimnego powietrza
k = 1,4, cp = 1,005, cv = 0,718 kJ/(kg·K)

Ciągłość i Bernoulli

Ciągłość
A1v1 = A2v2
Bernoulli
p1 + ½ρv1² + ρgz1 = p2 + ½ρv2² + ρgz2
Natężenie przepływu
Q = A·v
Pole koła
A = πd²/4

Tarcie w rurach i straty wysokości

Liczba Reynoldsa
Re = V·D/ν = ρ·V·D/μ
Charakter przepływu
laminarny poniżej Re ≈ 2000, turbulentny powyżej Re ≈ 4000
Prędkość średnia
V = Q/A, A = π·D²/4
Darcy–Weisbach
hf = f·(L/D)·V²/(2g)
Straty miejscowe
hm = ΣK·V²/(2g)
Typowe wartości K
wlot 0,5, wylot 1,0, kolano 90° ≈ 0,3
Spadek ciśnienia
Δp = ρ·g·hf
Dwa otwarte zbiorniki
Δz = (f·L/D + ΣK)·V²/(2g)
Swobodny wypływ
H = (1 + f·L/D + ΣK)·V²/(2g)

Pompy i moc pomp

Moc hydrauliczna
Ph = ρ·g·Q·H = Q·Δp
Moc na wale
Pshaft = Ph/ηpump
Moc elektryczna
Pel = Pshaft/ηmotor
Wysokość podnoszenia pompy
H = Δz + Δp/(ρ·g) + hL
Straty w instalacji
hL = (f·L/D + ΣK)·V²/(2g)

Przewodzenie ciepła

Prawo Fouriera
Q̇ = kAΔT/L
Opór cieplny przewodzenia
R = L/(kA)
Opór cieplny konwekcji
R = 1/(hA)
Opory szeregowo
Q̇ = ΔTtotal/ΣR
Współczynnik U
U = 1/(ΣR·A)

Koła i przekładnie zębate

Zazębione koła
n1·N1 = n2·N2
Moment skręcający
T2 = T1·(N2/N1)·η
Moc
P = T·ω, ω = 2π·n/60
Przekładnia wielostopniowa
przełożenie całkowite = iloczyn przełożeń poszczególnych stopni

Uchyb ustalony i regulacja PID

Typ układu
liczba biegunów G(s) w s = 0
Stałe uchybowe
Kp = lim G(s), Kv = lim s·G(s), Ka = lim s²·G(s) (s → 0)
Skok A
ess = A/(1 + Kp)
Narastanie liniowe A·t
ess = A/Kv
Parabola A·t²/2
ess = A/Ka
Regulator PID
C(s) = Kp + Ki/s + Kd·s
Cel drugiego rzędu
s² + 2ζωns + ωn²
Przeregulowanie i czasy
%OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd, ts ≈ 4/(ζωn)

Metoda

  1. Policz bieguny układu otwartego C(s)·G(s) w s = 0: to typ układu.
  2. Wyznacz stałą uchybową pasującą do wejścia: Kp dla skoku, Kv dla narastania liniowego, Ka dla paraboli.
  3. Zastosuj ess = A/(1 + Kp), A/Kv lub A/Ka albo rozwiąż to względem potrzebnego wzmocnienia.
  4. Przy wymaganiu dotyczącym tłumienia lub prędkości zapisz wielomian charakterystyczny układu zamkniętego i porównaj go z s² + 2ζω_n·s + ω_n².

2Zadania do ćwiczenia

1.Reakcje podpór belekŚredni
w = 10 kN/mP = 50 kNABMA​ = ?8 m

Wspornik jest utwierdzony w ścianie w A. Przenosi obciążenie równomiernie rozłożone w na całej długości i siłę skupioną P w odległości 8 m od ściany. Wyznacz wartość momentu utwierdzenia MA w ścianie.

L = 8 mP = 50 kN w 8 mw = 10 kN/m od 0 do 8 m
MA = kN·m
2.Reakcje podpór belekTrudny
w0​ = 20 kN/mP = 70 kNABMA​ = ?1,5 m2,5 m

Wspornik jest utwierdzony w ścianie w A. Przenosi obciążenie liniowo zmienne o wartości 20 kN/m przy ścianie, malejące do zera na swobodnym końcu B, oraz siłę skupioną P w odległości 1,5 m od ściany. Wyznacz wartość momentu utwierdzenia MA w ścianie.

L = 4 mw0 = 20 kN/mP = 70 kN w 1,5 m
MA = kN·m
3.Naprężenie i odkształcenie osioweŚredni
P = ?12500 mm200 mmØ20Ø50δ = 0,2 mm

Pręt składa się z odcinka (materiał: stal, średnica 20 mm, długość 500 mm, E = 200 GPa) połączonego końcami z odcinkiem (materiał: mosiądz, średnica 50 mm, długość 200 mm, E = 100 GPa). Jest zamocowany na jednym końcu i rozciągany na drugim. Wyznacz obciążenie osiowe P, które wydłuża go łącznie o 0,2 mm.

odcinek 1: stal, Ø20 mm, 500 mm, E = 200 GPaodcinek 2: mosiądz, Ø50 mm, 200 mm, E = 100 GPaδ = 0,2 mm
P = kN
4.Naprężenie i odkształcenie osioweTrudny
Ø 25 mmAluminium2 mØ 12 mmMosiądz1 mP = 15 kN

Pręt schodkowy jest utwierdzony w ścianie i obciążony na swobodnym końcu. Wyznacz całkowite wydłużenie δ.

P = 15 kNodcinek 1: Ø25 mm, 2 m, aluminium (E = 70 GPa)odcinek 2: Ø12 mm, 1 m, mosiądz (E = 100 GPa)
δ = mm
5.Skręcanie wałówŚredni
T = 750 N·md = 40 mmL = 0,5 mφ = ?

Pełny wał (mosiądz) jest utwierdzony w ścianie i skręcany momentem T na swobodnym końcu. Wyznacz kąt skręcenia φ swobodnego końca w stopniach.

T = 750 N·md = 40 mmL = 0,5 mG = 39 GPa
φ = °
6.Skręcanie wałówTrudny
do​ = ?50 kW, 1500 rpm

Wał drążony ma przenosić moc 50 kW przy 1500 rpm. Jego średnica wewnętrzna ma wynosić 0,6 średnicy zewnętrznej, a naprężenie styczne nie może przekroczyć 75 MPa. Wyznacz najmniejszą średnicę zewnętrzną do.

P = 50 kWN = 1500 rpmdi/do = 0,6τallow = 75 MPa
do = mm
7.Siła tnąca i moment zginającyŚredni
40 kN60 kNABCDEFG2 m2 m2 m2 m2 m2 m

Punkty od A do G dzielą belkę na sześć równych części. W którym punkcie moment zginający ma największą wartość bezwzględną, dodatni lub ujemny? Kliknij go na rysunku.

40 kN @ 8 m60 kN @ 10 m
  1. 1A (x = 0 m)
  2. 2B (x = 2 m)
  3. 3C (x = 4 m)
  4. 4D (x = 6 m)
  5. 5E (x = 8 m)
  6. 6F (x = 10 m)
  7. 7G (x = 12 m)
8.Siła tnąca i moment zginającyTrudny
w = 4 kN/mP = 80 kNAB1,5 m3,5 m

Belka swobodnie podparta przenosi obciążenie równomiernie rozłożone w na całej rozpiętości oraz siłę skupioną P. Wyznacz maksymalny moment zginający Mmax w belce.

L = 5 mP = 80 kN w 1,5 mw = 4 kN/m od 0 do 5 m
Mmax = kN·m
9.Koło Mohra i transformacja naprężeńŚredni
σx​ = 110σy​ = 100τxy​ = 30MPa

W pewnym punkcie obciążonej części występują σₓ = 110 MPa, σᵧ = 100 MPa i τₓᵧ = 30 MPa (rozciąganie dodatnie). Wyznacz maksymalne naprężenie styczne w płaszczyźnie τmax.

σx = 110 MPaσy = 100 MPaτxy = 30 MPa
τmax = MPa
10.Koło Mohra i transformacja naprężeńTrudny
σx​τxy​powierzchnia wału

Pełny wał o średnicy 20 mm przenosi osiową siłę ściskającą 80 kN wraz z momentem skręcającym 250 N·m. Na powierzchni wału, wyznacz największe naprężenie główne σ1.

d = 20 mmF = −80 kNT = 250 N·m
σ1 = MPa
11.Ugięcie belkiŚredni
w = 6 kN/mδL = 8 m

Belka swobodnie podparta przenosi obciążenie równomiernie rozłożone w na całej rozpiętości. Wyznacz maksymalne ugięcie δmax w środku rozpiętości.

w = 6 kN/mL = 8 mE = 200 GPa (stal)I = 84 × 10⁶ mm⁴
δmax = mm
12.Ugięcie belkiTrudny
P = 20 kNAB2,5 m2,5 m

Belka jest utwierdzona w A i swobodnie podparta w B. Przenosi siłę skupioną P w środku rozpiętości. EI jest stałe. Wyznacz reakcję RB na podporze przesuwnej.

L = 5 mP = 20 kN w 2,5 m
RB = kN
13.Wyboczenie słupówŚredni
PL = 2 m40 × 50 mmprzekrój poprzeczny

Słup o długości L z materiału (aluminium) jest przegubowo podparty na obu końcach. Korzystając ze wzoru Eulera, wyznacz jego siłę krytyczną Pcr.

L = 2 mE = 70 GPa (aluminium)b × h = 40 × 50 mm
Pcr = kN
14.Wyboczenie słupówTrudny
PL = 3 md = 60 mmprzekrój poprzeczny

Słup o długości L z materiału (aluminium) jest utwierdzony u podstawy i wolny u góry. Wyznacz jego smukłość i naprężenie krytyczne σcr ze wzoru Eulera.

L = 3 mE = 70 GPa (aluminium)d = 60 mm
σcr = MPa
15.Rzut ukośnyŚredni
v0​ = 38 m/s50°x = ?

Piłka zostaje wyrzucona pod kątem 50° do poziomu z płaskiego terenu. Pomiń opór powietrza. Wyznacz odległość poziomą x pokonaną, zanim piłka uderzy w ziemię.

v0 = 38 m/sθ = 50°g = 9,81 m/s²
x = m
16.Rzut ukośnyTrudny
h = 40 mv0​ = 14 m/s40°x = ?

Piłka zostaje wyrzucona pod kątem 40° do poziomu ze szczytu klifu o wysokości 40 m. Pomiń opór powietrza. Wyznacz odległość poziomą x pokonaną, zanim piłka uderzy w ziemię.

v0 = 14 m/sθ = 40°h = 40 mg = 9,81 m/s²
x = m
17.Praca i energiaŚredni
35°d = 17 mμ = 0,15

Klocek o masie 2,5 kg zsuwa się na odcinku 17 m po zboczu nachylonym pod kątem 35° do poziomu, przy współczynniku tarcia kinetycznego 0,15. Wyznacz pracę Wf wykonaną przeciw tarciu.

m = 2,5 kgθ = 35°μ = 0,15d = 17 mg = 9,81 m/s²
Wf = J
18.Praca i energiaTrudny
k = 3000 N/m30°μ = 0,150,5 kg

Sprężynę o sztywności 3000 N/m u podstawy zbocza nachylonego pod kątem 30° ściśnięto o 5 cm, opierając o nią klocek o masie 0,5 kg, i zwolniono. Współczynnik tarcia kinetycznego między klockiem a zboczem wynosi 0,15. Jak daleko d w górę zbocza od punktu zwolnienia przesunie się klocek, zanim się zatrzyma? (Pomiń krótki odcinek, na którym działa sprężyna, tzn. uwzględnij tarcie i siłę ciężkości na całej drodze d).

m = 0,5 kgk = 3000 N/mx = 5 cmθ = 30°μ = 0,15g = 9,81 m/s²
d = m
19.Drgania mechaniczneŚredni
tx25 mmLiczba cykli: 55,2 mm

W próbie drgań swobodnych tłumionego układu masa–sprężyna po maksimum przemieszczenia 25 mm po 5 cyklach następuje maksimum 5,2 mm. Wyznacz współczynnik tłumienia ζ.

x0 = 25 mmx5 = 5,2 mmn = 5
ζ =
20.Drgania mechaniczneTrudny
k1​ = 5k2​ = 20k3​ = 20m = 50 kg

Masa 50 kg wisi na sprężynie k₃ = 20 kN/m. Górny koniec k₃ jest przymocowany do sztywnej belki, która wisi na dwóch sprężynach obok siebie, k₁ = 5 kN/m i k₂ = 20 kN/m; belka pozostaje pozioma. Pomiń tłumienie oraz masy sprężyn i belki. Wyznacz częstotliwość drgań własnych fn w kierunku pionowym w Hz.

m = 50 kgk1 = 5 kN/mk2 = 20 kN/mk3 = 20 kN/m
fn = Hz
21.Silniki cieplne i chłodziarkiŚredni
TH​ = 850 °CTL​ = 30 °CEQH​QL​η = 40%

Silnik cieplny pobiera ciepło ze źródła o temperaturze 850 °C i oddaje ciepło do otoczenia o temperaturze 30 °C. Jego rzeczywista sprawność cieplna wynosi 40 %. Wyznacz jego sprawność egzergetyczną ηII, czyli stosunek sprawności rzeczywistej do sprawności Carnota.

TH = 850 °CTL = 30 °Cη = 40%
ηII = %
22.Silniki cieplne i chłodziarkiTrudny
TH​ = 55 °CTL​ = -5 °CHPQ̇H​ = 12 kWQ̇L​Ẇ = ?

Pompa ciepła dostarcza 12 kW ciepła do instalacji grzewczej o temperaturze 55 °C, pobierając ciepło ze źródła o temperaturze -5 °C. Jej rzeczywisty COP wynosi 50 % COP Carnota między tymi temperaturami. Wyznacz pobieraną przez nią moc elektryczną Ẇ.

TH = 55 °CTL = -5 °CηII = 50%Q̇H = 12 kW
Ẇ = kW
23.Obiegi silników cieplnych gazowychŚredni
vP1234qin​qout​

Idealny obieg Diesla ma stopień sprężania r = 16 i stopień obciążenia rc = v₃/v₂ = 2,2. Przyjmij założenia modelu zimnego powietrza: powietrze jest gazem doskonałym o stałych ciepłach właściwych, k = 1,4, cp = 1,005 kJ/(kg·K) i cv = 0,718 kJ/(kg·K). Wyznacz jego sprawność cieplną η w procentach.

r = 16rc = 2,2k = 1,4
η = %
24.Obiegi silników cieplnych gazowychTrudny
vP1234qin​qout​

Idealny obieg Braytona pracuje przy stosunku ciśnień w sprężarce 15. Przyjmij założenia modelu zimnego powietrza: powietrze jest gazem doskonałym o stałych ciepłach właściwych, k = 1,4, cp = 1,005 kJ/(kg·K) i cv = 0,718 kJ/(kg·K). Wyznacz jego sprawność cieplną η w procentach.

rp = 15k = 1,4
η = %
25.Ciągłość i BernoulliŚredni
12d1​ = 75 mm, p1​ = 390 kPad2​ = 25 mm, p2​ = ?

Woda płynie przez poziomą rurę, która zwęża się z 75 mm do 25 mm. W szerokim odcinku ciśnienie wynosi 390 kPa, a prędkość 1,25 m/s. Pomijając straty, wyznacz ciśnienie p2 w wąskim odcinku.

d1 = 75 mmd2 = 25 mmp1 = 390 kPav1 = 1,25 m/sρ = 1000 kg/m³
p2 = kPa
26.Ciągłość i BernoulliTrudny
pair​ = 150 kPav = ?d = 20 mmh = 4,5 m

Zamknięty zbiornik z wodą ma nad wodą sprężone powietrze o nadciśnieniu 150 kPa. Woda wypływa do atmosfery przez otwór o średnicy 20 mm, położony 4,5 m pod powierzchnią wody. Pomijając straty, wyznacz prędkość strumienia v. Przyjmij g = 9,81 m/s².

pair = 150 kPa (nadciśnienie)h = 4,5 md = 20 mmρ = 1000 kg/m³g = 9,81 m/s²
v = m/s
27.Tarcie w rurach i straty wysokościŚredni
Q = 4 L/sD = 50 mmL = 100 mf = 0,02

Woda płynie z natężeniem 4 L/s przez rurę o długości 100 m i średnicy wewnętrznej 50 mm. Przyjmij współczynnik tarcia Darcy’ego f = 0,02 i g = 9,81 m/s². Wyznacz stratę wysokości na tarcie hf.

Q = 4 L/sL = 100 mD = 50 mmf = 0,02g = 9,81 m/s²
hf = m
28.Tarcie w rurach i straty wysokościTrudny
K = 2Δz = 11 mL = 200 m, D = 200 mm, f = 0,0252× kolano (K = 0,3)

Dwa otwarte zbiorniki są połączone rurą o długości 200 m, średnicy wewnętrznej 200 mm i współczynniku tarcia Darcy’ego f = 0,025. Zwierciadło wody w górnym zbiorniku leży 11 m wyżej niż w dolnym. Straty miejscowe to: wlot o ostrej krawędzi (K = 0,5), kolana w liczbie 2 (K = 0,3 każde), zawór (K = 2) i wylot do dolnego zbiornika (K = 1,0). Przyjmij g = 9,81 m/s². Wyznacz natężenie przepływu Q w L/s.

Δz = 11 mL = 200 mD = 200 mmf = 0,025Kentrance = 0,5Kbend = 0,3 (×2)Kvalve = 2Kexit = 1,0
Q = L/s
29.Pompy i moc pompŚredni
PQ = 40 L/sH = 32 m

Pompa tłoczy 40 L/s wody (ρ = 1000 kg/m³) przy wysokości podnoszenia 32 m. Sprawność pompy wynosi 75%, a napędzającego ją silnika elektrycznego 90%. Przyjmij g = 9,81 m/s². Wyznacz moc elektryczną Pel pobieraną przez silnik, w kW.

Q = 40 L/sH = 32 mηpump = 75%ηmotor = 90%g = 9,81 m/s²
Pel = kW
30.Pompy i moc pompTrudny
PQ = 60 L/sΔz = 9 mL = 100 m, D = 200 mm, f = 0,015, ΣK = 8

Pompa przetłacza 60 L/s wody (ρ = 1000 kg/m³) z otwartego dolnego zbiornika do górnego zbiornika, którego zwierciadło wody leży 9 m wyżej. Rura ma długość 100 m, średnicę wewnętrzną 200 mm i współczynnik tarcia Darcy’ego f = 0,015; współczynniki strat miejscowych sumują się do ΣK = 8 (łącznie z wlotem i wylotem). Przyjmij g = 9,81 m/s². Wyznacz wysokość podnoszenia H, którą musi zapewnić pompa.

Q = 60 L/sΔz = 9 mL = 100 mD = 200 mmf = 0,015ΣK = 8
H = m
31.Przewodzenie ciepłaŚredni
Cegła, k = 0,72100 mmDrewno, k = 0,1725 mm24 °C3 °CQ̇ = ?

Ściana składa się z 100 mm (cegła) i 25 mm (drewno), a jej pole wynosi 6 m². Powierzchnie zewnętrzne mają 24 °C i 3 °C. Wyznacz strumień ciepła Q̇ przez ścianę.

L1 = 100 mm, k1 = 0,72 W/(m·K)L2 = 25 mm, k2 = 0,17 W/(m·K)A = 6 m²T1 = 24 °CT3 = 3 °C
Q̇ = W
32.Przewodzenie ciepłaTrudny
Styropian, k = 0,03550 mmCegła, k = 0,7250 mmpomieszczenie 21 °Cna zewnątrz -12 °CTs​ = ?Q̇

Ściana z 50 mm (styropian) i 50 mm (cegła) oddziela pomieszczenie o temperaturze 21 °C od powietrza zewnętrznego o temperaturze -12 °C. Współczynniki przejmowania ciepła wynoszą 7,7 W/(m²·K) wewnątrz i 25 W/(m²·K) na zewnątrz, a pole ściany 6 m². Wyznacz temperaturę Ts wewnętrznej powierzchni ściany.

L1 = 50 mm, k1 = 0,035 W/(m·K)L2 = 50 mm, k2 = 0,72 W/(m·K)A = 6 m²hi = 7,7 W/(m²·K)ho = 25 W/(m²·K)Troom = 21 °CTout = -12 °C
Ts = °C
33.Koła i przekładnie zębateŚredni
Nw​ = 48z = 1, 1450 rpm

Ślimak o liczbie zwojów 1 obraca się z prędkością 1450 rpm i napędza ślimacznicę o 48 zębach. Wyznacz prędkość obrotową nw ślimacznicy.

z = 1Nw = 48n = 1450 rpm
nw = rpm
34.Koła i przekładnie zębateTrudny
A = 45B = 60C = 16D = 25

Silnik o mocy 2,2 kW i prędkości 960 rpm napędza koło A (45 zębów), które współpracuje z kołem B (60 zębów). Koło C (16 zębów) na tym samym wale co B napędza koło D (25 zębów) na wale wyjściowym. Każda para kół ma sprawność 98 %. Wyznacz moment wyjściowy Tout.

P = 2,2 kWnin = 960 rpmA = 45B = 60C = 16D = 25ηmesh = 98%
Tout = N·m
35.Transformaty Laplace’aŚredni
ReIm−5−10−15−20−25−30−35

F(s) = (s + 31)/(s² + 14s + 48). Jej transformata odwrotna ma postać f(t) = A·e−6t + B·e−8t. Znajdź współczynnik A metodą rozkładu na ułamki proste.

F(s) = (s + 31)/(s² + 14s + 48)
A =
36.Transformaty Laplace’aTrudny
t (s)y(t)t = 0,8y(0,8) = ?

Układ spełnia równanie y″ + 7y′ + 6y = 15·u(t). Startuje ze stanu spoczynku (y(0) = 0, y′(0) = 0), a wejściem jest skok u(t) = 10 dla t ≥ 0. Za pomocą transformaty Laplace’a znajdź wyjście y(0,8).

y″ + 7y′ + 6y = 15·u(t)u(t) = 10y(0) = 0, y′(0) = 0t = 0,8 s
y(0.8) =
37.Schematy blokoweŚredni
R(s)+17s + 9G(s)C(s)0,25H−

Obiekt G(s) = 17/(s + 9) jest regulowany w pętli ujemnego sprzężenia zwrotnego, w której czujnik ma wzmocnienie H = 0,25. Na wejście R podano skok o wielkości 12. Znajdź wyjście w stanie ustalonym css.

G(s) = 17/(s + 9)H = 0,25R = 12
css =
38.Schematy blokoweTrudny
R(s)+64s(s + 9)G(s)C(s)2H−

Obiekt G(s) = 64/(s(s + 9)) jest w pętli ujemnego sprzężenia zwrotnego z czujnikiem o wzmocnieniu H = 2. Bieguny układu zamkniętego są zespolone. Znajdź ich część rzeczywistą σ.

G(s) = 64/(s(s + 9))H = 2
σ =
39.Odpowiedź skokowaŚredni
t (s)y(t)3tp​ = ?

Układ ma transmitancję G(s) = 138/(s² + 7s + 46) i otrzymuje jednostkowe wymuszenie skokowe. Znajdź czas szczytu tp odpowiedzi skokowej.

G(s) = 138/(s² + 7s + 46)
tp = s
40.Odpowiedź skokowaTrudny
+−12 VR = 78 ΩL = 200 mHC = 47 µFvC​

Szeregowy obwód RLC z R = 78 Ω, L = 200 mH i C = 47 µF jest dołączony do skoku 12 V w chwili t = 0. Kondensator jest nienaładowany i przed skokiem nie płynie prąd. Znajdź czas tp pierwszego szczytu napięcia, w ms.

R = 78 ΩL = 200 mHC = 47 µFV = 12 V
tp = ms
41.Uchyb ustalony i regulacja PIDŚredni
r(t) = 10 + 5t+5(s + 1)s(s + 1)(s + 10)C(s)−

Na pętlę ze sprzężeniem jednostkowym o G(s) = 5(s + 1)/(s(s + 1)(s + 10)) podano wejście r(t) = 10 + 5t (skok plus narastanie liniowe). Wyznacz uchyb ustalony ess.

G(s) = 5(s + 1)/(s(s + 1)(s + 10))r(t) = 10 + 5t
ess =
42.Uchyb ustalony i regulacja PIDTrudny
R(s)+Kp​ + Ki​/s2s + 1C(s)−

Regulator PI C(s) = Kp + Ki/s steruje obiektem G(s) = 2/(s + 1) w pętli ze sprzężeniem jednostkowym. Układ zamknięty ma mieć ωₙ = 3 rad/s i ζ = 1. Wyznacz wzmocnienie całkowania Ki.

G(s) = 2/(s + 1)ωₙ = 3 rad/sζ = 1
Ki =

3Roztwory

1.Reakcje podpór belek

Odpowiedź: MA = 720 kN·m

Równowaga wspornika: utwierdzenie może pchać w górę i również przeciwdziałać obrotowi, więc daje reakcję pionową RA i moment utwierdzenia MA. ΣFy = 0 daje RA, a momenty względem ściany dają MA = Σ(obciążenie × odległość od ściany). Obciążenie ciągłe działa jak jego wypadkowa w środku ciężkości wykresu obciążenia.

Wypadkowa obciążenia ciągłego

W = 10·8 = 80 kN, x̄ = 4 m

Moment względem ściany

MA = P·a + W·x̄ = 50·8 + 80·4 = 720 kN·m

Odpowiedź

MA = 720 kN·m
2.Reakcje podpór belek

Odpowiedź: MA = 158,3 kN·m

Równowaga wspornika: utwierdzenie może pchać w górę i również przeciwdziałać obrotowi, więc daje reakcję pionową RA i moment utwierdzenia MA. ΣFy = 0 daje RA, a momenty względem ściany dają MA = Σ(obciążenie × odległość od ściany). Obciążenie ciągłe działa jak jego wypadkowa w środku ciężkości wykresu obciążenia.

Wypadkowa obciążenia trójkątnego

W = ½·20·4 = 40 kN
x̄ = L/3 = 1,333 m

Moment względem ściany

MA = W·x̄ + P·a = 40·1,333 + 70·1,5 = 158,3 kN·m

Odpowiedź

MA = 158,3 kN·m
3.Naprężenie i odkształcenie osiowe

Odpowiedź: P = 22,28 kN

Dwa materiały: przy połączeniu szeregowym (jeden za drugim) każdy odcinek przenosi tę samą siłę, a wydłużenia się dodają: δ = P·Σ L/(E·A). Przy połączeniu równoległym (obok siebie między sztywnymi płytami) obie części skracają się o tyle samo, więc ich odkształcenia są równe: sztywniejsza część przejmuje większy udział, Pi = P·EiAi/ΣEA, a σi = Ei·ε. Pręt zamocowany między dwiema sztywnymi ścianami z obciążeniem w miejscu połączenia jest statycznie niewyznaczalny: dwa odcinki dzielą obciążenie proporcjonalnie do sztywności k = EA/L.

Pola

A₁ = 314,2 mm², A₂ = 1963 mm²

Podatności

L₁/(E₁A₁) = 500/(200000·314,2) = 7,958 × 10⁻⁶ mm/N
L₂/(E₂A₂) = 200/(100000·1963) = 1,019 × 10⁻⁶ mm/N

Obciążenie

P = 0,2/8,976 × 10⁻⁶ = 22281 N = 22,28 kN

Odpowiedź

P = 22,28 kN
4.Naprężenie i odkształcenie osiowe

Odpowiedź: δ = 2,199 mm

Pręty obciążone osiowo: każdy pryzmatyczny odcinek spełnia δ = PL/(AE). Dla odcinków połączonych szeregowo wydłużenia się sumują.

Ta sama siła w obu odcinkach

Swobodny koniec przenosi P, więc każdy przekrój pręta przenosi tę samą siłę osiową P. Zmieniają się tylko pola przekroju i materiały.

A1 = π·25²/4 = 490,9 mm², A2 = π·12²/4 = 113,1 mm²
σ1 = 30,56 MPa
σ2 = 132,6 MPa

Dodaj wydłużenia odcinków

δ1 = σ1L1/E1 = 30,56·2000/70000 = 0,8731 mm
δ2 = σ2L2/E2 = 132,6·1000/100000 = 1,326 mm
δ = δ1 + δ2 = 2,199 mm

Odpowiedź

δ = 2,199 mm
5.Skręcanie wałów

Odpowiedź: φ = 2,192 °

Skręcanie wałów okrągłych: naprężenie styczne rośnie liniowo od zera w środku do τmax = T·r/J na powierzchni, gdzie J = πd⁴/32 dla wału pełnego i J = π(do⁴ − di⁴)/32 dla drążonego. Kąt skręcenia to φ = TL/(GJ) w radianach. Moc i moment skręcający wiąże P = T·ω.

Biegunowy moment bezwładności

J = π·40⁴/32 = 251327 mm⁴

Kąt skręcenia

φ = TL/(GJ) = 750·1000·500/(39000·251327) = 0,03826 rad
φ = 0,03826·180/π = 2,192°

Odpowiedź

φ = 2,192°
6.Skręcanie wałów

Odpowiedź: do = 29,18 mm

Skręcanie wałów okrągłych: naprężenie styczne rośnie liniowo od zera w środku do τmax = T·r/J na powierzchni, gdzie J = πd⁴/32 dla wału pełnego i J = π(do⁴ − di⁴)/32 dla drążonego. Kąt skręcenia to φ = TL/(GJ) w radianach. Moc i moment skręcający wiąże P = T·ω.

Moment skręcający

ω = 2πN/60 = 2π·1500/60 = 157,1 rad/s
T = P/ω = 50000/157,1 = 318,3 N·m

Średnica zewnętrzna z naprężenia dopuszczalnego

T = 318,3 N·m = 318310 N·mm
1 − k⁴ = 1 − 0,6⁴ = 0,8704
do = ∛(16T/(π·τallow·(1 − k⁴))) = ∛(16·318310/(π·75·0,8704)) = 29,18 mm
di = 0,6·29,18 = 17,51 mm

Odpowiedź

do = 29,18 mm

Wał pełny wymagałby d = 27,86 mm; drążony jest nieco większy, ale znacznie lżejszy.

7.Siła tnąca i moment zginający

Odpowiedź: E (x = 8 m)

Największy moment zginający: nachylenie wykresu M zmienia się tam, gdzie zmienia się siła tnąca, więc jego maksima leżą tam, gdzie siła tnąca przechodzi przez zero: pod siłą skupioną, wewnątrz obciążenia ciągłego albo nad podporą przy wsporniku. W utwierdzeniu belki wspornikowej moment jest największy przy ścianie. Porównaj wartości M w tych miejscach, zarówno dodatnie (zginanie „w dół”), jak i ujemne.

Reakcje podpór

Zapisz momenty względem jednej podpory, a potem użyj bilansu sił pionowych.

R1 = 23,33 kN, R2 = 76,67 kN

Moment zginający w każdym punkcie

Zsumuj momenty wszystkich sił po jednej stronie przekroju. Zginanie „w dół” jest dodatnie, „w górę” ujemne.

MA = 0, MB = 46,67, MC = 93,33, MD = 140, ME = 186,7 kN·m, MF = 153,3, MG = 0

Odpowiedź

Punkt E: ME = 186,7 kN·m
8.Siła tnąca i moment zginający

Odpowiedź: Mmax = 94,5 kN·m

Metoda przekrojów: przetnij belkę tam, gdzie chcesz znać siły wewnętrzne, i zostaw jedną stronę. Siła tnąca V równoważy siły pionowe po tej stronie, a moment zginający M ich momenty względem przekroju. Moment maksymalny: ponieważ dM/dx = V, moment zginający osiąga ekstremum tam, gdzie siła tnąca przechodzi przez zero, albo w miejscu siły skupionej, albo wewnątrz obciążenia ciągłego.

Reakcje podpór

Równowaga całej belki.

Wypadkowa obciążenia ciągłego: W = w·ℓ = 4·5 = 20 kN, przyłożona w x̄ = 2,5 m
ΣMA = 0: RB·5 = 80·1,5 + 4·5·2,5
RB = 170/5 = 34 kN
ΣFy = 0: RA = 100 − 34 = 66 kN

Gdzie V = 0?

V(1,5⁻) = RA − w·1,5 = 60 kN
V(1,5⁺) = 60 − 80 = −20 kN
Siła tnąca zmienia znak w obrębie obciążenia, więc maksimum leży pod obciążeniem.

Moment maksymalny

Mmax = RA·1,5 − w·1,5²/2 = 66·1,5 − 4·1,5²/2 = 94,5 kN·m

Odpowiedź

Mmax = 94,5 kN·m

Występuje w x = 1,5 m od A.

9.Koło Mohra i transformacja naprężeń

Odpowiedź: τmax = 30,41 MPa

Transformacja płaskiego stanu naprężenia: na płaszczyźnie obróconej o kąt θ: σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ. Koło Mohra pokazuje wszystkie płaszczyzny naraz: jego środek leży w C = (σₓ + σᵧ)/2, a promień wynosi R = √(((σₓ − σᵧ)/2)² + τₓᵧ²). Naprężenia główne to C ± R, a największe naprężenie styczne w płaszczyźnie to R.

σxyCτmax​

Promień koła Mohra

τmax = R = √(((110 − 100)/2)² + 30²) = 30,41 MPa

Odpowiedź

τmax = 30,41 MPa
10.Koło Mohra i transformacja naprężeń

Odpowiedź: σ1 = 76,49 MPa

Transformacja płaskiego stanu naprężenia: na płaszczyźnie obróconej o kąt θ: σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ. Koło Mohra pokazuje wszystkie płaszczyzny naraz: jego środek leży w C = (σₓ + σᵧ)/2, a promień wynosi R = √(((σₓ − σᵧ)/2)² + τₓᵧ²). Naprężenia główne to C ± R, a największe naprężenie styczne w płaszczyźnie to R.

Naprężenia na powierzchni

σ = F/A = −80000/(π·20²/4) = −80000/314,2 = −254,6 MPa
τ = 16T/(πd³) = 16·250000/(π·20³) = 159,2 MPa
στxyCσ1​σ2​

Koło Mohra

C = σ/2 = −127,3 MPa
R = √((−127,3)² + 159,2²) = 203,8 MPa
σ1 = C + R = 76,49 MPa
σ2 = C − R = −331,1 MPa

Odpowiedź

σ1 = 76,49 MPa
11.Ugięcie belki

Odpowiedź: δmax = 19,05 mm

Ugięcie belek: ugięcia otrzymuje się, całkując dwukrotnie krzywiznę M/(EI). Typowe przypadki są stablicowane, na przykład PL³/(3EI) dla wspornika z siłą na końcu i 5wL⁴/(384EI) dla belki swobodnie podpartej z obciążeniem równomiernym na całej rozpiętości. Ugięcia od różnych obciążeń się sumują (superpozycja).

Spójne jednostki

Licz w N i mm: P w N, w w N/mm (kN/m = N/mm), L w mm, E w N/mm² (GPa × 1000), I w mm⁴.

w = 6 N/mm, L = 8000 mm
EI = 1,68 × 10¹³ N·mm²

Ugięcie

δmax = 5wL⁴/(384EI) = 5·6·8000⁴/(384·1,68 × 10¹³) = 19,05 mm

Odpowiedź

δmax = 19,05 mm
12.Ugięcie belki

Odpowiedź: RB = 6,25 kN

Ugięcie belek: ugięcia otrzymuje się, całkując dwukrotnie krzywiznę M/(EI). Typowe przypadki są stablicowane, na przykład PL³/(3EI) dla wspornika z siłą na końcu i 5wL⁴/(384EI) dla belki swobodnie podpartej z obciążeniem równomiernym na całej rozpiętości. Ugięcia od różnych obciążeń się sumują (superpozycja). W belce statycznie niewyznaczalnej usuń podporę nadliczbową i zażądaj, aby ugięcie w tym miejscu było równe zeru (warunek zgodności).

Usuń podporę przesuwną

Bez podpory przesuwnej belka jest wspornikiem. Jej swobodny koniec ugiąłby się pod obciążeniem.

δload = P(L/2)²(3L − L/2)/(6EI) = 5PL³/(48EI)

Warunek zgodności przemieszczeń na podporze

Sama reakcja podpory RB podnosi swobodny koniec o RBL³/(3EI). Podpora utrzymuje B przy zerowym ugięciu, więc oba ugięcia muszą się znosić. EI się skraca.

RB·L³/(3EI) = 5PL³/(48EI) → RB = 5P/16 = 5·20/16 = 6,25 kN

Odpowiedź

RB = 6,25 kN
13.Wyboczenie słupów

Odpowiedź: Pcr = 46,06 kN

Wyboczenie Eulera: smukły słup traci stateczność przy Pcr = π²EI/(KL)², gdzie KL to długość wyboczeniowa: K = 1 dla przegub–przegub, 2 dla utwierdzenie–wolny koniec, 0,7 dla utwierdzenie–przegub i 0,5 dla utwierdzenie–utwierdzenie. Wybocza się względem osi o najmniejszym I. Wzór dotyczy słupów smukłych, których naprężenie krytyczne jest mniejsze od granicy plastyczności.

Najmniejszy moment bezwładności

Prostokąt wybocza się względem słabszej osi, z mniejszym wymiarem do sześcianu.

Imin = h·b³/12 = 50·40³/12 = 266667 mm⁴

Długość wyboczeniowa

KL = 1·2000 = 2000 mm

Siła krytyczna

Pcr = π²EI/(KL)² = π²·70000·266667/2000² = 46058 N = 46,06 kN

Odpowiedź

Pcr = 46,06 kN
14.Wyboczenie słupów

Odpowiedź: σcr = 4,318 MPa

Wyboczenie Eulera: smukły słup traci stateczność przy Pcr = π²EI/(KL)², gdzie KL to długość wyboczeniowa: K = 1 dla przegub–przegub, 2 dla utwierdzenie–wolny koniec, 0,7 dla utwierdzenie–przegub i 0,5 dla utwierdzenie–utwierdzenie. Wybocza się względem osi o najmniejszym I. Wzór dotyczy słupów smukłych, których naprężenie krytyczne jest mniejsze od granicy plastyczności.

Najmniejszy moment bezwładności

Przekrój kołowy ma to samo I względem każdej osi.

I = πd⁴/64 = π·60⁴/64 = 636173 mm⁴

Promień bezwładności i smukłość

A = 2827 mm²
r = √(I/A) = √(636173/2827) = 15 mm
KL = 2·3000 = 6000 mm
KL/r = 6000/15 = 400

Naprężenie krytyczne

σcr = π²E/(KL/r)² = π²·70000/400² = 4,318 MPa

Odpowiedź

σcr = 4,318 MPa
15.Rzut ukośny

Odpowiedź: x = 145 m

Rzut ukośny: bez oporu powietrza x = vxt i y = h + vyt − ½gt². Oba kierunki są niezależne, łączy je tylko ten sam czas.

Rozłóż prędkość początkową

Prędkość pozioma jest stała; ruch pionowy ma stałe przyspieszenie g = 9,81 m/s² w dół.

vx = v0·cos θ = 38·cos 50° = 24,43 m/s
vy = v0·sin θ = 38·sin 50° = 29,11 m/s

Czas lotu

Start i lądowanie na tej samej wysokości: t = 2vy/g.

t = 2·29,11/9,81 = 5,935 s

Odległość pozioma

x = vx·t = 24,43·5,935 = 145 m

Odpowiedź

x = 145 m
16.Rzut ukośny

Odpowiedź: x = 42,01 m

Rzut ukośny: bez oporu powietrza x = vxt i y = h + vyt − ½gt². Oba kierunki są niezależne, łączy je tylko ten sam czas.

Rozłóż prędkość początkową

Prędkość pozioma jest stała; ruch pionowy ma stałe przyspieszenie g = 9,81 m/s² w dół.

vx = v0·cos θ = 14·cos 40° = 10,72 m/s
vy = v0·sin θ = 14·sin 40° = 8,999 m/s

Czas lotu

Rozwiąż y(t) = h + vy·t − ½gt² = 0, biorąc dodatni pierwiastek.

0 = 40 + 8,999·t − 4,905·t²
t = (vy + √(vy² + 2gh))/g = (8,999 + √(80,98 + 784,8))/9,81 = 3,917 s

Odległość pozioma

x = vx·t = 10,72·3,917 = 42,01 m

Odpowiedź

x = 42,01 m
17.Praca i energia

Odpowiedź: Wf = 51,23 J

Zasada pracy i energii: zmiana energii kinetycznej równa się pracy wszystkich sił. Zapisanie grawitacji i sprężyn jako energii potencjalnej (mgh i ½kx²) daje bilans: energia początkowa = energia końcowa + energia stracona na tarcie (μN·d).

Siła tarcia

N = mg·cos θ = 24,53·0,8192 = 20,09 N
F = μN = 0,15·20,09 = 3,013 N

Praca przeciw tarciu

Wf = F·d = 3,013·17 = 51,23 J

Odpowiedź

Wf = 51,23 J
18.Praca i energia

Odpowiedź: d = 1,214 m

Zasada pracy i energii: zmiana energii kinetycznej równa się pracy wszystkich sił. Zapisanie grawitacji i sprężyn jako energii potencjalnej (mgh i ½kx²) daje bilans: energia początkowa = energia końcowa + energia stracona na tarcie (μN·d).

Energia sprężyny

x = 5 cm = 0,05 m
½kx² = ½·3000·0,05² = 3,75 J

Bilans energii wzdłuż zbocza

m·g = 0,5·9,81 = 4,905 N
½kx² = mg·d·(sin θ + μ·cos θ)
d = 3,75/(4,905·(0,5 + 0,15·0,866)) = 1,214 m

Odpowiedź

d = 1,214 m
19.Drgania mechaniczne

Odpowiedź: ζ = 0,04992

Drgania swobodne układu masa–sprężyna: pulsacja drgań własnych nietłumionych to ωn = √(k/m) w rad/s, a fn = ωn/(2π) w Hz. Tłumik wiskotyczny c daje współczynnik tłumienia ζ = c/cc, przy tłumieniu krytycznym cc = 2·√(k·m) = 2·m·ωn. Układ z tłumieniem podkrytycznym (ζ < 1) drga z pulsacją ωd = ωn·√(1 − ζ²). Sztywności sprężyn połączonych równolegle się dodaje, k = k₁ + k₂; dla sprężyn szeregowych dodaje się odwrotności, 1/k = 1/k₁ + 1/k₂.

Logarytmiczny dekrement tłumienia

δ = (1/5)·ln(25/5,2) = 0,314

Współczynnik tłumienia

ζ = δ/√(4π² + δ²) = 0,314/√(39,48 + 0,09862) = 0,04992

Odpowiedź

ζ = 0,04992
20.Drgania mechaniczne

Odpowiedź: fn = 2,373 Hz

Drgania swobodne układu masa–sprężyna: pulsacja drgań własnych nietłumionych to ωn = √(k/m) w rad/s, a fn = ωn/(2π) w Hz. Tłumik wiskotyczny c daje współczynnik tłumienia ζ = c/cc, przy tłumieniu krytycznym cc = 2·√(k·m) = 2·m·ωn. Układ z tłumieniem podkrytycznym (ζ < 1) drga z pulsacją ωd = ωn·√(1 − ζ²). Sztywności sprężyn połączonych równolegle się dodaje, k = k₁ + k₂; dla sprężyn szeregowych dodaje się odwrotności, 1/k = 1/k₁ + 1/k₂.

Sztywność zastępcza

Tylko k1 i k2 są połączone równolegle. Ta para jest szeregowo z k3, więc połącz je przez odwrotności.

k12 = 5 + 20 = 25 kN/m
keq = 25·20/(25 + 20) = 11,11 kN/m

Częstotliwość drgań własnych

ωn = √(11111/50) = 14,91 rad/s
fn = 14,91/(2π) = 2,373 Hz

Odpowiedź

fn = 2,373 Hz
21.Silniki cieplne i chłodziarki

Odpowiedź: ηII = 54,79 %

Second-law limits: no cycle working between reservoirs at absolute temperatures TH and TL can beat the reversible (Carnot) cycle: ηmax = 1 − TL/TH, COPR,max = TL/(TH − TL) and COPHP,max = TH/(TH − TL). Energy is always conserved: QH = W + QL.

Temperatury bezwzględne

TH = 850 + 273,15 = 1123 K
TL = 30 + 273,15 = 303,1 K

Sprawność Carnota

ηCarnot = 1 − 303,1/1123 = 73,01 %

Sprawność egzergetyczna

ηII = η/ηCarnot = 40/73,01 = 54,79 %

Odpowiedź

ηII = 54,79 %
22.Silniki cieplne i chłodziarki

Odpowiedź: Ẇ = 4,388 kW

Second-law limits: no cycle working between reservoirs at absolute temperatures TH and TL can beat the reversible (Carnot) cycle: ηmax = 1 − TL/TH, COPR,max = TL/(TH − TL) and COPHP,max = TH/(TH − TL). Energy is always conserved: QH = W + QL.

Temperatury bezwzględne

TH = 55 + 273,15 = 328,1 K
TL = -5 + 273,15 = 268,1 K

COP Carnota i rzeczywisty

COPmax = 328,1/(328,1 − 268,1) = 5,469
COP = 0,5·5,469 = 2,735

Moc pobierana

Ẇ = Q̇H/COP = 12/2,735 = 4,388 kW

Odpowiedź

Ẇ = 4,388 kW
23.Obiegi silników cieplnych gazowych

Odpowiedź: η = 60,42 %

Obieg Diesla (silnik o zapłonie samoczynnym): 1→2 sprężanie izentropowe, 2→3 doprowadzenie ciepła przy stałym ciśnieniu, 3→4 rozprężanie izentropowe, 4→1 odprowadzenie ciepła przy stałej objętości. Przy stopniu sprężania r = v₁/v₂ i stopniu obciążenia rc = v₃/v₂: T₂ = T₁·rk−1, T₃ = rc·T₂, qin = cp·(T₃ − T₂) oraz η = 1 − (1/rk−1)·(rck − 1)/(k·(rc − 1)).

Człon stopnia obciążenia

(2,21,4 − 1)/(1,4·(2,2 − 1)) = 1,2

Sprawność cieplna

η = 1 − 0,3299·1,2 = 60,42%

Odpowiedź

η = 60,42%
24.Obiegi silników cieplnych gazowych

Odpowiedź: η = 53,87 %

Obieg Braytona (turbina gazowa): 1→2 sprężanie izentropowe, 2→3 doprowadzenie ciepła przy stałym ciśnieniu, 3→4 rozprężanie izentropowe w turbinie, 4→1 odprowadzenie ciepła przy stałym ciśnieniu. Przy stosunku ciśnień rp = P₂/P₁: T₂ = T₁·rp(k−1)/k, T₄ = T₃/rp(k−1)/k oraz η = 1 − 1/rp(k−1)/k. Praca netto w = cp·[(T₃ − T₄) − (T₂ − T₁)]. Przy sprawnościach izentropowych sprężarka wymaga ws/ηc, a turbina daje ηt·ws.

Sprawność cieplna

η = 1 − 1/150,2857 = 1 − 1/2,168 = 53,87%

Odpowiedź

η = 53,87%
25.Ciągłość i Bernoulli

Odpowiedź: p2 = 327,5 kPa

Ciągłość i Bernoulli: dla płynu nieściśliwego natężenie przepływu jest takie samo w każdym przekroju, Q = A₁v₁ = A₂v₂. Wzdłuż linii prądu bez strat p + ½ρv² + ρgz jest stałe, więc tam, gdzie płyn jest szybszy lub wyżej, ciśnienie jest niższe.

Ciągłość i Bernoulli

Licz w Pa, a potem podziel przez 1000, aby dostać kPa.

v2 = v1·(d1/d2)² = 1,25·(75/25)² = 11,25 m/s
½ρ(v1² − v2²) = ½·1000·(1,25² − 11,25²)/1000 = −62,5 kPa
p2 = p1 + ½ρ(v1² − v2²) = 390 + (−62,5) = 327,5 kPa

Odpowiedź

p2 = 327,5 kPa
26.Ciągłość i Bernoulli

Odpowiedź: v = 19,71 m/s

Continuity and Bernoulli: for an incompressible fluid the flow rate is the same at every section, Q = A₁v₁ = A₂v₂. Along a streamline without losses, p + ½ρv² + ρgz is constant, so where the fluid is faster or higher its pressure is lower.

Równanie Bernoulliego od powierzchni do strumienia

v = √(2·(150000/1000 + 9,81·4,5)) = 19,71 m/s

Odpowiedź

v = 19,71 m/s
27.Tarcie w rurach i straty wysokości

Odpowiedź: hf = 8,461 m

Przepływ w rurach: liczba Reynoldsa Re = V·D/ν = ρ·V·D/μ odróżnia przepływ laminarny (Re poniżej ok. 2000) od turbulentnego (powyżej ok. 4000). Tarcie wzdłuż rury powoduje stratę wysokości Darcy’ego–Weisbacha hf = f·(L/D)·V²/(2g). Kształtki, zawory, wlot i wylot dodają straty miejscowe hm = ΣK·V²/(2g). Między dwiema otwartymi powierzchniami zbiorników równanie energii daje Δz = hf + hm, z uwzględnioną stratą na wylocie K = 1.

Prędkość średnia

A = π·0,05²/4 = 0,001963 m²
V = Q/A = 0,004/0,001963 = 2,037 m/s

Wysokość prędkości

V²/(2g) = 2,037²/(2·9,81) = 0,2115 m

Darcy–Weisbach

hf = f·(L/D)·V²/(2g) = 0,02·(100/0,05)·2,037²/(2·9,81) = 8,461 m

Odpowiedź

hf = 8,461 m
28.Tarcie w rurach i straty wysokości

Odpowiedź: Q = 85,56 L/s

Przepływ w rurach: liczba Reynoldsa Re = V·D/ν = ρ·V·D/μ odróżnia przepływ laminarny (Re poniżej ok. 2000) od turbulentnego (powyżej ok. 4000). Tarcie wzdłuż rury powoduje stratę wysokości Darcy’ego–Weisbacha hf = f·(L/D)·V²/(2g). Kształtki, zawory, wlot i wylot dodają straty miejscowe hm = ΣK·V²/(2g). Między dwiema otwartymi powierzchniami zbiorników równanie energii daje Δz = hf + hm, z uwzględnioną stratą na wylocie K = 1.

Współczynniki strat

f·L/D = 0,025·200/0,2 = 25
ΣK = 0,5 + 2·0,3 + 2 + 1,0 = 4,1

Równanie energii

Obie powierzchnie są otwarte do atmosfery i praktycznie w spoczynku, więc całe Δz jest tracone w rurze.

11 = (25 + 4,1)·V²/(2·9,81)
V = √(2·9,81·11/29,1) = 2,723 m/s

Natężenie przepływu

Q = V·π·D²/4 = 2,723·0,03142 = 0,08556 m³/s = 85,56 L/s

Odpowiedź

Q = 85,56 L/s
29.Pompy i moc pomp

Odpowiedź: Pel = 18,6 kW

Moc pompy: moc przekazana wodzie to moc hydrauliczna Ph = ρ·g·Q·H = Q·Δp. Wał musi dostarczyć więcej, bo pompa ma straty: Pshaft = Ph/ηpump, a silnik pobiera Pel = Pshaft/ηmotor. Wysokość podnoszenia pompy H wynika z równania energii: geometryczna wysokość podnoszenia + wysokość ciśnienia + wszystkie straty tarcia i miejscowe.

Moc hydrauliczna

Ph = ρ·g·Q·H = 1000·9,81·0,04·32 = 12,56 kW

Moc na wale

Pshaft = 12,56/0,75 = 16,74 kW

Moc elektryczna

Pel = 16,74/0,9 = 18,6 kW

Odpowiedź

Pel = 18,6 kW
30.Pompy i moc pomp

Odpowiedź: H = 11,88 m

Moc pompy: moc przekazana wodzie to moc hydrauliczna Ph = ρ·g·Q·H = Q·Δp. Wał musi dostarczyć więcej, bo pompa ma straty: Pshaft = Ph/ηpump, a silnik pobiera Pel = Pshaft/ηmotor. Wysokość podnoszenia pompy H wynika z równania energii: geometryczna wysokość podnoszenia + wysokość ciśnienia + wszystkie straty tarcia i miejscowe.

Prędkość

V = Q/A = 0,06/0,03142 = 1,91 m/s
V²/(2g) = 0,1859 m

Straty wysokości

hL = (f·L/D + ΣK)·V²/(2g) = (7,5 + 8)·0,1859 = 2,882 m

Wysokość podnoszenia pompy

Oba zwierciadła są otwarte i praktycznie w spoczynku, więc człony ciśnienia i prędkości znikają.

H = Δz + hL = 9 + 2,882 = 11,88 m

Odpowiedź

H = 11,88 m
31.Przewodzenie ciepła

Odpowiedź: Q̇ = 440,6 W

Opór cieplny: ustalony przepływ ciepła przez warstwy połączone szeregowo zachowuje się jak prąd w rezystorach szeregowych: Q̇ = ΔT/ΣR, gdzie R = L/(kA) dla warstwy przewodzącej i R = 1/(hA) dla powierzchni z konwekcją. Spadek temperatury na każdej warstwie jest proporcjonalny do jej oporu.

Opory warstw

R1 = L1/(k1A) = (100/1000)/(0,72·6) = 0,02315 K/W
R2 = L2/(k2A) = (25/1000)/(0,17·6) = 0,02451 K/W

Strumień ciepła

Q̇ = (T1 − T3)/(R1 + R2) = 21/0,04766 = 440,6 W

Odpowiedź

Q̇ = 440,6 W
32.Przewodzenie ciepła

Odpowiedź: Ts = 18,43 °C

Opór cieplny: ustalony przepływ ciepła przez warstwy połączone szeregowo zachowuje się jak prąd w rezystorach szeregowych: Q̇ = ΔT/ΣR, gdzie R = L/(kA) dla warstwy przewodzącej i R = 1/(hA) dla powierzchni z konwekcją. Spadek temperatury na każdej warstwie jest proporcjonalny do jej oporu.

Opory szeregowo

Ri = 1/(hiA) = 1/(7,7·6) = 0,02165 K/W
R1 = L1/(k1A) = (50/1000)/(0,035·6) = 0,2381 K/W
R2 = L2/(k2A) = (50/1000)/(0,72·6) = 0,01157 K/W
Ro = 1/(hoA) = 1/(25·6) = 0,006667 K/W
ΣR = 0,278 K/W

Straty ciepła

Q̇ = (21 − (-12))/0,278 = 118,7 W
Ts = Troom − Q̇·Ri = 21 − 118,7·0,02165 = 18,43 °C

Odpowiedź

Ts = 18,43 °C
33.Koła i przekładnie zębate

Odpowiedź: nw = 30,21 rpm

Gear trains: meshing gears have the same speed at the teeth, so n₁·N₁ = n₂·N₂: the gear with fewer teeth turns faster. Gears on one shaft turn at the same speed. Power is passed on (minus losses), so torque goes up when speed goes down: T₂ = T₁·(N₂/N₁)·η, and P = T·ω with ω = 2π·n/60.

Prędkość ślimacznicy

nw = n·z/Nw = 1450·1/48 = 30,21 rpm

Odpowiedź

nw = 30,21 rpm
34.Koła i przekładnie zębate

Odpowiedź: Tout = 43,79 N·m

Gear trains: meshing gears have the same speed at the teeth, so n₁·N₁ = n₂·N₂: the gear with fewer teeth turns faster. Gears on one shaft turn at the same speed. Power is passed on (minus losses), so torque goes up when speed goes down: T₂ = T₁·(N₂/N₁)·η, and P = T·ω with ω = 2π·n/60.

Prędkość obrotowa wyjściowa

nout = 960·(45/60)·(16/25) = 460,8 rpm
ωout = 2π·460,8/60 = 48,25 rad/s

Moc i moment wyjściowy

Pout = 2200·0,98² = 2113 W
Tout = Pout/ωout = 2113/48,25 = 43,79 N·m

Odpowiedź

Tout = 43,79 N·m
35.Transformaty Laplace’a

Odpowiedź: A = 12,5

Transformata Laplace’a: F(s) = ∫₀∞ f(t)·e−st dt zamienia pochodne na mnożenie przez s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Podstawowe pary: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), a przesunięcie s → s + σ mnoży funkcję czasu przez e−σt. Aby odwrócić wymierną F(s), rozłóż ją na ułamki proste: residuum prostego bieguna p wynosi A = [(s − p)·F(s)] przy s = p. Twierdzenie o wartości końcowej: f(∞) = lims→0 s·F(s), słuszne tylko wtedy, gdy wszystkie bieguny s·F(s) leżą w lewej półpłaszczyźnie. Twierdzenie o wartości początkowej: f(0⁺) = lims→∞ s·F(s).

Rozłóż mianownik na czynniki

s² + 14s + 48 = (s + 6)(s + 8)
F(s) = A/(s + 6) + B/(s + 8)

Metoda zakrywania

Pomnóż przez (s + 6) i podstaw s = −6: składnik z B znika.

A = (s + 31)/(s + 8) |s = −6 = 25/2 = 12,5

Odpowiedź

A = 12,5
(B = −11,5, więc f(t) = 12,5·e−6t − 11,5·e−8t)
36.Transformaty Laplace’a

Odpowiedź: y(0.8) = 11,56

Transformata Laplace’a: F(s) = ∫₀∞ f(t)·e−st dt zamienia pochodne na mnożenie przez s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Podstawowe pary: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), a przesunięcie s → s + σ mnoży funkcję czasu przez e−σt. Aby odwrócić wymierną F(s), rozłóż ją na ułamki proste: residuum prostego bieguna p wynosi A = [(s − p)·F(s)] przy s = p. Twierdzenie o wartości końcowej: f(∞) = lims→0 s·F(s), słuszne tylko wtedy, gdy wszystkie bieguny s·F(s) leżą w lewej półpłaszczyźnie. Twierdzenie o wartości początkowej: f(0⁺) = lims→∞ s·F(s).

Przekształcenie równania

Y(s) = 150/(s(s + 1)(s + 6))

Ułamki proste

A = 150/(1·6) = 25
B = 150/((−1)·5) = −30
C = 150/((−6)·(−5)) = 5

Wydalone

y(t) = 25 − 30·e−t + 5·e−6t
y(0,8) = 11,56

Odpowiedź

y(0,8) = 11,56
37.Schematy blokowe

Odpowiedź: css = 15,4

Algebra schematów blokowych: bloki połączone szeregowo mnożą się (G₁·G₂), bloki połączone równolegle dodają się (G₁ + G₂), a pętla ze ścieżką prostą G i ścieżką sprzężenia zwrotnego H redukuje się do T = G/(1 + G·H) dla ujemnego sprzężenia zwrotnego i G/(1 − G·H) dla dodatniego. Najpierw redukuj pętle wewnętrzne, potem idź na zewnątrz. Bieguny układu zamkniętego są pierwiastkami równania charakterystycznego 1 + G(s)H(s) = 0. Dla stabilnej pętli wymuszenie skokowe o wielkości A daje w stanie ustalonym wyjście A·T(0); przy jednostkowym sprzężeniu zwrotnym uchyb ustalony przy skoku wynosi A/(1 + Kp), gdzie Kp = G(0).

Transmitancja układu zamkniętego

T(s) = G/(1 + G·H) = 17/(s + 9 + 4,25) = 17/(s + 13,25)

Wynik

css = R·T(0) = 12·17/13,25 = 15,4

Odpowiedź

css = 15,4
38.Schematy blokowe

Odpowiedź: σ = −4,5

Algebra schematów blokowych: bloki połączone szeregowo mnożą się (G₁·G₂), bloki połączone równolegle dodają się (G₁ + G₂), a pętla ze ścieżką prostą G i ścieżką sprzężenia zwrotnego H redukuje się do T = G/(1 + G·H) dla ujemnego sprzężenia zwrotnego i G/(1 − G·H) dla dodatniego. Najpierw redukuj pętle wewnętrzne, potem idź na zewnątrz. Bieguny układu zamkniętego są pierwiastkami równania charakterystycznego 1 + G(s)H(s) = 0. Dla stabilnej pętli wymuszenie skokowe o wielkości A daje w stanie ustalonym wyjście A·T(0); przy jednostkowym sprzężeniu zwrotnym uchyb ustalony przy skoku wynosi A/(1 + Kp), gdzie Kp = G(0).

Równanie charakterystyczne

s(s + 9) + 64·2 = 0 → s² + 9s + 128 = 0

Pierwiastki

s = −4,5 ± j·√(128 − 20,25) = −4,5 ± j10,38
σ = −4,5

Odpowiedź

σ = −4,5
39.Odpowiedź skokowa

Odpowiedź: tp = 0,5408 s

Odpowiedź skokowa: układ pierwszego rzędu G(s) = K/(τs + 1) odpowiada na skok o wielkości U przebiegiem y(t) = K·U·(1 − e−t/τ): po jednej stałej czasowej τ osiąga 63,2 % wartości końcowej K·U, a po 4τ 98 %. Układ drugiego rzędu G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) ma bieguny s = −ζωₙ ± jωd, gdzie ωd = ωₙ√(1 − ζ²). Dla 0 < ζ < 1 odpowiedź skokowa ma przeregulowanie: przeregulowanie procentowe %OS = 100·e−ζπ/√(1−ζ²), czas szczytu tp = π/ωd, a 2-procentowy czas ustalania ts ≈ 4/(ζωₙ). Porównaj mianownik z s² + 2ζωₙs + ωₙ², aby odczytać ωₙ i ζ. Dla układu masa–sprężyna–tłumik ωₙ = √(k/m) i ζ = c/(2√(km)); dla szeregowego obwodu RLC ωₙ = 1/√(LC) i ζ = (R/2)·√(C/L).

Porównaj z postacią standardową

Mianownik musi odpowiadać s² + 2ζωₙs + ωₙ².

ωₙ² = 46 → ωₙ = 6,782 rad/s
2ζωₙ = 7 → ζ = 7/(2·6,782) = 0,516

Wynik

ωd = 6,782·√(1 − 0,516²) = 5,809 rad/s → tp = π/ωd = 0,5408 s

Odpowiedź

tp = 0,5408 s
40.Odpowiedź skokowa

Odpowiedź: tp = 12,02 ms

Odpowiedź skokowa: układ pierwszego rzędu G(s) = K/(τs + 1) odpowiada na skok o wielkości U przebiegiem y(t) = K·U·(1 − e−t/τ): po jednej stałej czasowej τ osiąga 63,2 % wartości końcowej K·U, a po 4τ 98 %. Układ drugiego rzędu G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) ma bieguny s = −ζωₙ ± jωd, gdzie ωd = ωₙ√(1 − ζ²). Dla 0 < ζ < 1 odpowiedź skokowa ma przeregulowanie: przeregulowanie procentowe %OS = 100·e−ζπ/√(1−ζ²), czas szczytu tp = π/ωd, a 2-procentowy czas ustalania ts ≈ 4/(ζωₙ). Porównaj mianownik z s² + 2ζωₙs + ωₙ², aby odczytać ωₙ i ζ. Dla układu masa–sprężyna–tłumik ωₙ = √(k/m) i ζ = c/(2√(km)); dla szeregowego obwodu RLC ωₙ = 1/√(LC) i ζ = (R/2)·√(C/L).

Pulsacja własna i współczynnik tłumienia

ωₙ = 1/√(LC) = 1/√(0,2·4,7 × 10⁻⁵) = 326,2 rad/s
ζ = (R/2)·√(C/L) = (78/2)·√(4,7 × 10⁻⁵/0,2) = 0,5979

Wynik

ωd = 261,5 rad/s
tp = π/ωd = 12,02 ms

Odpowiedź

tp = 12,02 ms
41.Uchyb ustalony i regulacja PID

Odpowiedź: ess = 10

Uchyb ustalony: w pętli ze sprzężeniem jednostkowym o transmitancji układu otwartego G(s) uchyb wynosi E(s) = R(s)/(1 + G(s)), a z twierdzenia o wartości końcowej ess = lims→0 s·E(s). Typ układu to liczba biegunów G(s) w s = 0. Stałe uchybowe to Kp = lim G(s), Kv = lim s·G(s) oraz Ka = lim s²·G(s), wszystkie dla s → 0. Skok o wartości A daje ess = A/(1 + Kp), narastanie liniowe A·t daje A/Kv, a parabola A·t²/2 daje A/Ka. Pętla typu 0 nadąża za skokiem ze skończonym uchybem, pętla typu 1 nadąża za nim dokładnie, a za narastaniem liniowym ze skończonym uchybem.

Stała prędkościowa

Kv = lim s·G(s) = 0,5

Część skokowa

Kp = ∞ → estep = 0

Uchyb dla narastania liniowego

ess = 5/0,5 = 10

Odpowiedź

ess = 10
42.Uchyb ustalony i regulacja PID

Odpowiedź: Ki = 4,5

Regulacja P, PI i PID: regulator C(s) = Kp + Ki/s + Kd·s działa na uchyb. Wzmocnienie proporcjonalne przyspiesza układ i zmniejsza uchyb obiektu typu 0, ale nigdy go nie usuwa. Człon całkujący dodaje biegun w s = 0: typ układu rośnie o jeden, a skok jest śledzony bez uchybu. Człon różniczkujący zwiększa tłumienie. Aby dobrać nastawy, zapisz wielomian charakterystyczny układu zamkniętego i porównaj go z s² + 2ζωns + ωn². Wtedy %OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd, gdzie ωd = ωn√(1 − ζ²), oraz ts ≈ 4/(ζωn).

Wielomian charakterystyczny

1 + C(s)G(s) = 0, pomnożone przez s(s + 1).

s² + (1 + 2Kp)s + 2Ki

Wielomian docelowy

s² + 2·1·3s + 3² = s² + 6s + 9

Dopasuj współczynniki

2Ki = 9 → Ki = 4,5

Odpowiedź

Ki = 4,5

Człon całkujący sprawia, że pętla jest typu 1, więc skok jest śledzony bez uchybu ustalonego.

To niezależne materiały do ćwiczeń. ExamProblems nie jest powiązane z organizacją przeprowadzającą ten egzamin ani przez nią rekomendowane, a nazwa egzaminu służy wyłącznie do opisu treści. Egzamin FE exam prowadzi NCEES, organizacja niezwiązana z ExamProblems, która go nie popiera.

ExamProblems · examproblems.com/courses/fe-exam-mechanical