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
P1​ = 35 kNP2​ = 25 kNABRA​ = ?RB​5 m3 m2 m

Belka swobodnie podparta ma podporę przegubową nieprzesuwną w A i przesuwną w B. Wyznacz pionową reakcję podpory RA.

L = 10 mP1 = 35 kN w 5 mP2 = 25 kN w 8 m
RA = kN
2.Reakcje podpór belekTrudny
w = 5 kN/mP = 15 kNABRA​RB​ = ?1 m1,5 m1,5 m

Belka swobodnie podparta ma podporę przegubową nieprzesuwną w A i przesuwną w B. Wyznacz pionową reakcję podpory RB.

L = 4 mP = 15 kN w 2,5 mw = 5 kN/m od 0 do 1 m
RB = kN
3.Naprężenie i odkształcenie osioweŚredni
P = ?12600 mm300 mmØ25Ø50δ = 0,3 mm

Pręt składa się z odcinka (materiał: stal, średnica 25 mm, długość 600 mm, E = 200 GPa) połączonego końcami z odcinkiem (materiał: brąz, średnica 50 mm, długość 300 mm, E = 110 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,3 mm.

odcinek 1: stal, Ø25 mm, 600 mm, E = 200 GPaodcinek 2: brąz, Ø50 mm, 300 mm, E = 110 GPaδ = 0,3 mm
P = kN
4.Naprężenie i odkształcenie osioweTrudny
Ø 30 mmMosiądz2 mØ 20 mmMiedź1 mP = 30 kN

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

P = 30 kNodcinek 1: Ø30 mm, 2 m, mosiądz (E = 100 GPa)odcinek 2: Ø20 mm, 1 m, miedź (E = 117 GPa)
δ = mm
5.Skręcanie wałówŚredni
T = 400 N·md = 50 mmL = 1,2 mφ = ?

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

T = 400 N·md = 50 mmL = 1,2 mG = 26 GPa
φ = °
6.Skręcanie wałówTrudny
Tdo​ = 50 mm, di​ = 38 mm

Wał drążony przenosi 50 kW przy 1800 obr/min. Wyznacz maksymalne naprężenie styczne τmax w wale.

P = 50 kWN = 1800 rpmdo = 50 mmdi = 38 mm
τmax = MPa
7.Siła tnąca i moment zginającyŚredni
w = 8 kN/mABprzekrój3,5 m1,5 m

Wyznacz moment zginający Mx w przekroju x = 3,5 m od podpory A. Przyjmij dodatnie momenty rozciągające dolne włókna i dodatnią siłę tnącą skierowaną w górę na lewej ściance.

L = 5 mw = 8 kN/m od 0 do 5 mx = 3,5 m
Mx = kN·m
8.Siła tnąca i moment zginającyTrudny
25 kN60 kN20 kNABCDEFG1 m1 m1 m1 m1 m1 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.

25 kN @ 0 m60 kN @ 6 m20 kN @ 2 m
  1. 1A (x = 0 m)
  2. 2B (x = 1 m)
  3. 3C (x = 2 m)
  4. 4D (x = 3 m)
  5. 5E (x = 4 m)
  6. 6F (x = 5 m)
  7. 7G (x = 6 m)
9.Koło Mohra i transformacja naprężeńŚredni
σx​ = 100σy​ = 60τxy​ = 20MPa

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

σx = 100 MPaσy = 60 MPaτxy = 20 MPa
τmax = MPa
10.Koło Mohra i transformacja naprężeńTrudny
σx​ = 90σy​ = −80τxy​ = 50MPa

W pewnym punkcie obciążonej części występują σₓ = 90 MPa, σᵧ = −80 MPa i τₓᵧ = 50 MPa (rozciąganie dodatnie). Wyznacz naprężenie normalne σθ na płaszczyźnie, której normalna jest obrócona o 20° przeciwnie do ruchu wskazówek zegara względem osi x.

σx = 90 MPaσy = −80 MPaτxy = 50 MPaθ = 20°
σθ = MPa
11.Ugięcie belkiŚredni
P = 2 kNABδ2 m3 m

Wspornik o długości L jest utwierdzony w A i przenosi siłę skupioną P w odległości a = 2 m od A. Wyznacz ugięcie δP pod siłą.

P = 2 kNa = 2 mL = 5 mE = 70 GPa (aluminium)I = 12 × 10⁶ mm⁴
δP = mm
12.Ugięcie belkiTrudny
w = 12 kN/mP = 11 kNδL = 4 m

Belka swobodnie podparta przenosi obciążenie równomiernie rozłożone w i siłę skupioną P w środku rozpiętości. Korzystając z superpozycji, wyznacz ugięcie w środku δmax.

P = 11 kNw = 12 kN/mL = 4 mE = 200 GPa (stal)I = 160 × 10⁶ mm⁴
δmax = mm
13.Wyboczenie słupówŚredni
PL = 2,5 md = 50 mmprzekrój poprzeczny

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

L = 2,5 mE = 200 GPa (stal)d = 50 mm
Pcr = kN
14.Wyboczenie słupówTrudny
P = 15 kNL = 1,5 md = ?przekrój poprzeczny

Pełny słup okrągły o długości L z materiału (stal) jest utwierdzony u podstawy i przegubowo podparty u góry. Ma przenieść 15 kN przy współczynniku bezpieczeństwa 2 względem wyboczenia. Wyznacz najmniejszą średnicę d.

P = 15 kNn = 2L = 1,5 mE = 200 GPa (stal)
d = mm
15.Rzut ukośnyŚredni
v0​ = 28 m/s30°zasięg

Piłka zostaje wyrzucona pod kątem 30° do poziomu z płaskiego terenu. Pomiń opór powietrza. Wyznacz czas lotu t, zanim piłka uderzy w ziemię.

v0 = 28 m/sθ = 30°g = 9,81 m/s²
t = s
16.Rzut ukośnyTrudny
h = 5 mv0​ = 30 m/s35°x = ?

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

v0 = 30 m/sθ = 35°h = 5 mg = 9,81 m/s²
x = m
17.Praca i energiaŚredni
30 km/h

Samochód jadący z prędkością 30 km/h po poziomej drodze gwałtownie hamuje i z zablokowanymi kołami sunie aż do zatrzymania. Współczynnik tarcia kinetycznego między oponami a drogą wynosi 0,4. Wyznacz drogę hamowania d.

v = 30 km/hμ = 0,4g = 9,81 m/s²
d = m
18.Praca i energiaTrudny
h = 2 mL = 3 mμ = 0,25k = 800 N/m

Klocek o masie 1 kg puszczono bez prędkości początkowej z wysokości 2 m na gładkiej rampie. Na dole przejeżdża przez szorstki poziomy pas o długości 3 m (μ = 0,25), a potem trafia na sprężynę (k = 800 N/m) na gładkiej powierzchni. Wyznacz maksymalne ugięcie x sprężyny.

m = 1 kgh = 2 mL = 3 mμ = 0,25k = 800 N/mg = 9,81 m/s²
x = cm
19.Drgania mechaniczneŚredni
k = 40 kN/mζ = 0,6m = 2 kg

Masa 2 kg na sprężynie o sztywności 40 kN/m ma tłumik wiskotyczny dający współczynnik tłumienia 0,6. Wyznacz częstotliwość drgań swobodnych tłumionych fd w Hz.

m = 2 kgk = 40 kN/mζ = 0,6
fd = Hz
20.Drgania mechaniczneTrudny
k1​ = 20k2​ = 30m = 100 kg

Masa 100 kg wisi na dwóch sprężynach obok siebie (równolegle) o sztywnościach k₁ = 20 kN/m i k₂ = 30 kN/m. Masa pozostaje pozioma. Pomiń tłumienie i masy sprężyn. Wyznacz częstotliwość drgań własnych fn w kierunku pionowym w Hz.

m = 100 kgk1 = 20 kN/mk2 = 30 kN/m
fn = Hz
21.Silniki cieplne i chłodziarkiŚredni
TH​ = 750 °CTL​ = 15 °CEQH​QL​η = 43%

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

TH = 750 °CTL = 15 °Cη = 43%
ηII = %
22.Silniki cieplne i chłodziarkiTrudny
TH​ = 21 °CTL​ = 5 °CHPQ̇H​ = 6 kWQ̇L​Ẇ = ?

Pompa ciepła utrzymuje w domu 21 °C, gdy na zewnątrz jest 5 °C, dostarczając do domu 6 kW ciepła. Wyznacz minimalną moc napędową Ẇmin.

TH = 21 °CTL = 5 °CQ̇H = 6 kW
Ẇmin = kW
23.Obiegi silników cieplnych gazowychŚredni
vP1234qin​qout​

Idealny obieg Diesla ma stopień sprężania r = 14 i stopień obciążenia rc = v₃/v₂ = 2,5. 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 = 14rc = 2,5k = 1,4
η = %
24.Obiegi silników cieplnych gazowychTrudny
vP1234qin​qout​

Elektrownia z turbiną gazową pracuje w obiegu Braytona przy stosunku ciśnień w sprężarce 8. Powietrze wpływa do sprężarki w temperaturze 295 K, a do turbiny w temperaturze 1400 K. Obie maszyny są izentropowe. 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 pracę netto wnet na kilogram powietrza, w kJ/kg.

rp = 8T1 = 295 KT3 = 1400 Kk = 1,4cp = 1,005 kJ/(kg·K)
wnet = kJ/kg
25.Ciągłość i BernoulliŚredni
12d1​ = 100 mm, p1​ = 240 kPad2​ = 75 mm, p2​ = ?

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

d1 = 100 mmd2 = 75 mmp1 = 240 kPav1 = 1 m/sρ = 1000 kg/m³
p2 = kPa
26.Ciągłość i BernoulliTrudny
H = ?d1​ = 75 mmp1​ = 50 kPad2​ = 20 mm

Dysza fontanny jest skierowana pionowo w górę. Tuż pod dyszą rura zasilająca ma średnicę 75 mm i nadciśnienie 50 kPa. Średnica wylotu dyszy wynosi 20 mm, a różnica wysokości na długości dyszy jest pomijalna. Pomijając straty i opór powietrza, wyznacz wysokość H, na jaką strumień wznosi się nad dyszę. Przyjmij g = 9,81 m/s².

d1 = 75 mmp1 = 50 kPa (nadciśnienie)d2 = 20 mmρ = 1000 kg/m³g = 9,81 m/s²
H = m
27.Tarcie w rurach i straty wysokościŚredni
V = 2,2 m/sD = 75 mmL = 300 mf = 0,03

Woda płynie z prędkością średnią 2,2 m/s przez rurę o długości 300 m i średnicy wewnętrznej 75 mm. Przyjmij współczynnik tarcia Darcy’ego f = 0,03 i g = 9,81 m/s². Wyznacz stratę wysokości na tarcie hf.

V = 2,2 m/sL = 300 mD = 75 mmf = 0,03g = 9,81 m/s²
hf = m
28.Tarcie w rurach i straty wysokościTrudny
K = 0,2Δz = 22 mL = 200 m, D = 250 mm, f = 0,0223× kolano (K = 0,3)

Dwa otwarte zbiorniki są połączone rurą o długości 200 m, średnicy wewnętrznej 250 mm i współczynniku tarcia Darcy’ego f = 0,022. Zwierciadło wody w górnym zbiorniku leży 22 m wyżej niż w dolnym. Straty miejscowe to: wlot o ostrej krawędzi (K = 0,5), kolana w liczbie 3 (K = 0,3 każde), zawór (K = 0,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 = 22 mL = 200 mD = 250 mmf = 0,022Kentrance = 0,5Kbend = 0,3 (×3)Kvalve = 0,2Kexit = 1,0
Q = L/s
29.Pompy i moc pompŚredni
PQ = 30 L/sH = 32 m

Pompa o sprawności 85% tłoczy 30 L/s wody (ρ = 1000 kg/m³) przy wysokości podnoszenia 32 m. Przyjmij g = 9,81 m/s². Wyznacz potrzebną moc na wale Pshaft, w kW.

Q = 30 L/sH = 32 mη = 85%g = 9,81 m/s²
Pshaft = kW
30.Pompy i moc pompTrudny
PQ = 50 L/sΔz = 12 mL = 200 m, D = 150 mm, f = 0,018, ΣK = 5

Pompa przetłacza 50 L/s wody (ρ = 1000 kg/m³) z otwartego dolnego zbiornika do górnego zbiornika, którego zwierciadło wody leży 12 m wyżej. Rura ma długość 200 m, średnicę wewnętrzną 150 mm i współczynnik tarcia Darcy’ego f = 0,018; współczynniki strat miejscowych sumują się do ΣK = 5 (łącznie z wlotem i wylotem). Przyjmij g = 9,81 m/s². Sprawność pompy wynosi 82%. Wyznacz moc na wale Pshaft w kW.

Q = 50 L/sΔz = 12 mL = 200 mD = 150 mmf = 0,018ΣK = 5η = 82%
Pshaft = kW
31.Przewodzenie ciepłaŚredni
Wełna mineralna, k = 0,04100 mmBeton, k = 1,4150 mm17 °C7 °CTi​ = ?Q̇

Ściana składa się z 100 mm (wełna mineralna) i 150 mm (beton), a jej pole wynosi 4 m². Powierzchnie zewnętrzne mają 17 °C i 7 °C. Wyznacz temperaturę Ti na styku dwóch warstw.

L1 = 100 mm, k1 = 0,04 W/(m·K)L2 = 150 mm, k2 = 1,4 W/(m·K)A = 4 m²T1 = 17 °CT3 = 7 °C
Ti = °C
32.Przewodzenie ciepłaTrudny
Beton, k = 1,4250 mmStyropian, k = 0,03575 mmpomieszczenie 20 °Cna zewnątrz -9 °CTs​ = ?Q̇

Ściana z 250 mm (beton) i 75 mm (styropian) oddziela pomieszczenie o temperaturze 20 °C od powietrza zewnętrznego o temperaturze -9 °C. Współczynniki przejmowania ciepła wynoszą 8 W/(m²·K) wewnątrz i 25 W/(m²·K) na zewnątrz, a pole ściany 4 m². Wyznacz temperaturę Ts wewnętrznej powierzchni ściany.

L1 = 250 mm, k1 = 1,4 W/(m·K)L2 = 75 mm, k2 = 0,035 W/(m·K)A = 4 m²hi = 8 W/(m²·K)ho = 25 W/(m²·K)Troom = 20 °CTout = -9 °C
Ts = °C
33.Koła i przekładnie zębateŚredni
N1​ = 12N2​ = 24

Silnik o mocy 5,5 kW pracujący z prędkością 1450 rpm napędza zębnik o 12 zębach, który zazębia się z kołem o 24 zębach. Sprawność pary kół wynosi 98%. Wyznacz moment wyjściowy T2 na wale koła.

P = 5,5 kWn1 = 1450 rpmN1 = 12N2 = 24η = 98%
T2 = N·m
34.Koła i przekładnie zębateTrudny
A = 12B = 16C = 16D = 50

W przekładni dwustopniowej koło A (12 zębów) na wale wejściowym, obracające się z prędkością 2900 rpm, napędza koło B (16 zębów). Koło C (16 zębów) jest osadzone na tym samym wale co B i napędza koło D (50 zębów) na wale wyjściowym. Wyznacz prędkość obrotową na wyjściu nout.

nin = 2900 rpmA = 12B = 16C = 16D = 50
nout = rpm
35.Transformaty Laplace’aŚredni
ReIm−2−4−6−8−10−12

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

F(s) = (2s + 22)/(s² + 9s + 14)
A =
36.Transformaty Laplace’aTrudny
t (s)y(t)t = 0,1y(0,1) = ?

Układ spełnia równanie y″ + 6y′ + 8y = 32·u(t), przy y(0) = 2 i y′(0) = 0. W chwili t = 0 podano na wejście skok u(t) = 1. Za pomocą transformaty Laplace’a znajdź y(0,1).

y″ + 6y′ + 8y = 32·u(t)y(0) = 2y′(0) = 0u(t) = 1t = 0,1 s
y(0.1) =
37.Schematy blokoweŚredni
R(s)+K(s + 4)(s + 6)G(s)C(s)−

Obiekt G(s) = K/((s + 4)(s + 6)) jest w pętli jednostkowego ujemnego sprzężenia zwrotnego. Znajdź wzmocnienie K, dla którego układ zamknięty jest tłumiony krytycznie (dwa równe bieguny rzeczywiste).

G(s) = K/((s + 4)(s + 6))jednostkowe ujemne sprzężenie zwrotne
K =
38.Schematy blokoweTrudny
R(s)+68s(s + 5)G(s)C(s)0,5H−

Obiekt G(s) = 68/(s(s + 5)) jest w pętli ujemnego sprzężenia zwrotnego z czujnikiem o wzmocnieniu H = 0,5. Bieguny układu zamkniętego są zespolone. Znajdź ich część urojoną ωd (wartość dodatnią), w rad/s.

G(s) = 68/(s(s + 5))H = 0,5
ωd = rad/s
39.Odpowiedź skokowaŚredni
t (s)y(t)2

Układ ma transmitancję G(s) = 182/(s² + 3s + 91) i otrzymuje jednostkowe wymuszenie skokowe. Wyznacz współczynnik tłumienia ζ.

G(s) = 182/(s² + 3s + 91)
ζ =
40.Odpowiedź skokowaTrudny
t (s)y(t)110 %

Układ drugiego rzędu bez zer ma mieć odpowiedź skokową z przeregulowaniem 10 %. Znajdź współczynnik tłumienia ζ.

%OS = 10 %
ζ =
41.Uchyb ustalony i regulacja PIDŚredni
R(s)+K = 551s + 1C(s)−

Proces G(s) = 5/(1s + 1) jest sterowany wzmocnieniem proporcjonalnym K = 5 w pętli ze sprzężeniem jednostkowym. Wyznacz uchyb ustalony ess dla skoku, jako procent wartości skoku.

G(s) = 5/(1s + 1)K = 5
ess = %
42.Uchyb ustalony i regulacja PIDTrudny
R(s)+K = 510s(s + 10)C(s)−

Regulator proporcjonalny K = 5 steruje obiektem G(s) = 10/(s(s + 10)) w pętli ze sprzężeniem jednostkowym. Znajdź przeregulowanie procentowe %OS odpowiedzi skokowej.

G(s) = 10/(s(s + 10))K = 5
%OS = %

3Roztwory

1.Reakcje podpór belek

Odpowiedź: RA = 22,5 kN

Równowaga statyczna: ciało w spoczynku ma ΣFx = 0, ΣFy = 0 i ΣM = 0 względem dowolnego punktu. Moment względem jednej podpory usuwa jej reakcję z równania.

P1​ = 35 kNP2​ = 25 kNABRA​RB​5 m3 m2 m

Schemat sił

Zastąp podpory siłami reakcji. Podpora przesuwna daje jedną reakcję pionową; przy samych obciążeniach pionowych podpora nieprzesuwna też daje jedną reakcję pionową.

Równowaga

Momenty względem A, dodatnie przeciwnie do ruchu wskazówek zegara. Potem równowaga sił pionowych.

ΣMA = 0: RB·10 = 35·5 + 25·8
RB = 375/10 = 37,5 kN
ΣFy = 0: RA = 60 − 37,5 = 22,5 kN

Odpowiedź

RA = 22,5 kN
2.Reakcje podpór belek

Odpowiedź: RB = 10 kN

Równowaga statyczna: ciało w spoczynku ma ΣFx = 0, ΣFy = 0 i ΣM = 0 względem dowolnego punktu. Moment względem jednej podpory usuwa jej reakcję z równania.

w = 5 kN/mP = 15 kNABRA​RB​1 m1,5 m1,5 m

Schemat sił

Zastąp podpory siłami reakcji. Podpora przesuwna daje jedną reakcję pionową; przy samych obciążeniach pionowych podpora nieprzesuwna też daje jedną reakcję pionową.

Równowaga

Momenty względem A, dodatnie przeciwnie do ruchu wskazówek zegara. Potem równowaga sił pionowych.

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

Odpowiedź

RB = 10 kN
3.Naprężenie i odkształcenie osiowe

Odpowiedź: P = 40 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₁ = 490,9 mm², A₂ = 1963 mm²

Podatności

L₁/(E₁A₁) = 600/(200000·490,9) = 6,112 × 10⁻⁶ mm/N
L₂/(E₂A₂) = 300/(110000·1963) = 1,389 × 10⁻⁶ mm/N

Obciążenie

P = 0,3/7,501 × 10⁻⁶ = 39997 N = 40 kN

Odpowiedź

P = 40 kN
4.Naprężenie i odkształcenie osiowe

Odpowiedź: δ = 1,665 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 = π·30²/4 = 706,9 mm², A2 = π·20²/4 = 314,2 mm²
σ1 = 42,44 MPa
σ2 = 95,49 MPa

Dodaj wydłużenia odcinków

δ1 = σ1L1/E1 = 42,44·2000/100000 = 0,8488 mm
δ2 = σ2L2/E2 = 95,49·1000/117000 = 0,8162 mm
δ = δ1 + δ2 = 1,665 mm

Odpowiedź

δ = 1,665 mm
5.Skręcanie wałów

Odpowiedź: φ = 1,724 °

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 = π·50⁴/32 = 613592 mm⁴

Kąt skręcenia

φ = TL/(GJ) = 400·1000·1200/(26000·613592) = 0,03009 rad
φ = 0,03009·180/π = 1,724°

Odpowiedź

φ = 1,724°
6.Skręcanie wałów

Odpowiedź: τmax = 16,22 MPa

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π·1800/60 = 188,5 rad/s
T = P/ω = 50000/188,5 = 265,3 N·m

Biegunowy moment bezwładności

J = π(50⁴ − 38⁴)/32 = 408885 mm⁴

Maksymalne naprężenie styczne

τmax = T·(do/2)/J = 265258·25/408885 = 16,22 MPa

Odpowiedź

τmax = 16,22 MPa
7.Siła tnąca i moment zginający

Odpowiedź: Mx = 21 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.

Reakcje podpór

Równowaga całej belki.

Wypadkowa obciążenia ciągłego: W = w·ℓ = 8·5 = 40 kN, przyłożona w x̄ = 2,5 m
ΣMA = 0: RB·5 = 8·5·2,5
RB = 100/5 = 20 kN
ΣFy = 0: RA = 40 − 20 = 20 kN
w = 8 kN/mABRA​ = 20RB​ = 20przekrój3,5 m1,5 m

Przetnij w x = 3,5 m i zostaw lewą część

Uwzględniasz tylko siły na lewo od przekroju. Fragment obciążenia ciągłego działa jak jedna siła w środku tego fragmentu.

M = RA·x − (każde obciążenie na lewo od przekroju × jego odległość od przekroju)
Mx = 20·3,5 − 8·3,5·1,75 = 21 kN·m

Odpowiedź

Mx = 21 kN·m
8.Siła tnąca i moment zginający

Odpowiedź: F (x = 5 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 = 31,25 kN, R2 = 73,75 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 = −25, MC = −18,75, MD = −32,5, ME = −46,25, MF = −60 kN·m, MG = 0

Odpowiedź

Punkt F: MF = −60 kN·m
9.Koło Mohra i transformacja naprężeń

Odpowiedź: τmax = 28,28 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 = √(((100 − 60)/2)² + 20²) = 28,28 MPa

Odpowiedź

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

Odpowiedź: σθ = 102,3 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.

Równanie transformacji

σθ = 5 + 85·0,766 + 50·0,6428 = 102,3 MPa

Odpowiedź

σθ = 102,3 MPa
11.Ugięcie belki

Odpowiedź: δP = 6,349 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⁴.

P = 2000 N, a = 2000 mm, L = 5000 mm
EI = 8,4 × 10¹¹ N·mm²

Ugięcie pod siłą

δP = Pa³/(3EI) = 2000·2000³/(3·8,4 × 10¹¹) = 6,349 mm

Odpowiedź

δP = 6,349 mm
12.Ugięcie belki

Odpowiedź: δmax = 1,708 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⁴.

P = 11000 N, w = 12 N/mm, L = 4000 mm, EI = 3,2 × 10¹³ N·mm²

Każde obciążenie osobno

δP = PL³/(48EI) = 11000·4000³/(48·3,2 × 10¹³) = 0,4583 mm
δw = 5wL⁴/(384EI) = 5·12·4000⁴/(384·3,2 × 10¹³) = 1,25 mm

Dodaj

δmax = 0,4583 + 1,25 = 1,708 mm

Odpowiedź

δmax = 1,708 mm
13.Wyboczenie słupów

Odpowiedź: Pcr = 387,6 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

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

I = πd⁴/64 = π·50⁴/64 = 306796 mm⁴

Długość wyboczeniowa

KL = 0,5·2500 = 1250 mm

Siła krytyczna

Pcr = π²EI/(KL)² = π²·200000·306796/1250² = 387578 N = 387,6 kN

Odpowiedź

Pcr = 387,6 kN
14.Wyboczenie słupów

Odpowiedź: d = 24,17 mm

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.

Wymagana siła krytyczna

Pcr = n·P = 2·15000 = 30000 N
KL = 0,7·1500 = 1050 mm

Wyznacz średnicę

d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·30000·1050²/(π³·200000))^(1/4) = 24,17 mm

Odpowiedź

d = 24,17 mm

Sprawdzenie: smukłość KL/r wynosi 173,8, więc słup jest smukły i wzór Eulera obowiązuje.

15.Rzut ukośny

Odpowiedź: t = 2,854 s

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 θ = 28·cos 30° = 24,25 m/s
vy = v0·sin θ = 28·sin 30° = 14 m/s

Czas lotu

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

t = 2·14/9,81 = 2,854 s

Odpowiedź

t = 2,854 s
16.Rzut ukośny

Odpowiedź: x = 92,84 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 θ = 30·cos 35° = 24,57 m/s
vy = v0·sin θ = 30·sin 35° = 17,21 m/s

Czas lotu

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

0 = 5 + 17,21·t − 4,905·t²
t = (vy + √(vy² + 2gh))/g = (17,21 + √(296,1 + 98,1))/9,81 = 3,778 s

Odległość pozioma

x = vx·t = 24,57·3,778 = 92,84 m

Odpowiedź

x = 92,84 m
17.Praca i energia

Odpowiedź: d = 8,849 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).

Prędkość w jednostkach podstawowych

v = 30/3,6 = 8,333 m/s

Bilans energii

½mv² = μmg·d → d = v²/(2μg) = 8,333²/(2·0,4·9,81) = 8,849 m

Odpowiedź

d = 8,849 m

Podwojenie prędkości wydłuża drogę hamowania czterokrotnie.

18.Praca i energia

Odpowiedź: x = 17,51 cm

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 po szorstkim pasie

Rampa jest gładka, więc energię zabiera tylko pas.

mgh − μmgL = 9,81·(2 − 0,25·3) = 12,26 J

Maksymalne ugięcie

½kx² = 12,26 → x = √(2·12,26/800) = 0,1751 m = 17,51 cm

Odpowiedź

x = 17,51 cm
19.Drgania mechaniczne

Odpowiedź: fd = 18,01 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₂.

Częstotliwość drgań własnych nietłumionych

ωn = √(40000/2) = 141,4 rad/s
fn = 141,4/(2π) = 22,51 Hz

Częstotliwość drgań tłumionych

fd = fn·√(1 − ζ²) = 22,51·√(1 − 0,6²) = 18,01 Hz

Odpowiedź

fd = 18,01 Hz
20.Drgania mechaniczne

Odpowiedź: fn = 3,559 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

Sprężyny obok siebie mają to samo ugięcie, a ich siły się sumują, więc ich sztywności po prostu się dodaje: k = k1 + k2.

keq = 20 + 30 = 50 kN/m

Częstotliwość drgań własnych

ωn = √(50000/100) = 22,36 rad/s
fn = 22,36/(2π) = 3,559 Hz

Odpowiedź

fn = 3,559 Hz
21.Silniki cieplne i chłodziarki

Odpowiedź: ηII = 59,86 %

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 = 750 + 273,15 = 1023 K
TL = 15 + 273,15 = 288,1 K

Sprawność Carnota

ηCarnot = 1 − 288,1/1023 = 71,84 %

Sprawność egzergetyczna

ηII = η/ηCarnot = 43/71,84 = 59,86 %

Odpowiedź

ηII = 59,86 %
22.Silniki cieplne i chłodziarki

Odpowiedź: Ẇmin = 0,3264 kW

Ograniczenia z drugiej zasady: żaden obieg pracujący między zbiornikami o temperaturach bezwzględnych TH i TL nie przewyższy obiegu odwracalnego (Carnota): ηmax = 1 − TL/TH, COPR,max = TL/(TH − TL) i COPHP,max = TH/(TH − TL). Energia jest zawsze zachowana: QH = W + QL.

Temperatury bezwzględne

TH = 21 + 273,15 = 294,1 K
TL = 5 + 273,15 = 278,1 K

Najlepszy możliwy COP

COPHP,max = TH/(TH − TL) = 294,1/(294,1 − 278,1) = 18,38

Minimalna moc

Najmniej pracy potrzebuje maszyna odwracalna, która ma największy COP.

Ẇmin = Q̇H/COP = 6/18,38 = 0,3264 kW

Odpowiedź

Ẇmin = 0,3264 kW
23.Obiegi silników cieplnych gazowych

Odpowiedź: η = 56,81 %

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,51,4 − 1)/(1,4·(2,5 − 1)) = 1,241

Sprawność cieplna

η = 1 − 0,348·1,241 = 56,81%

Odpowiedź

η = 56,81%
24.Obiegi silników cieplnych gazowych

Odpowiedź: wnet = 389,7 kJ/kg

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.

Temperatury izentropowe

T2 = 295·80,2857 = 534,4 K
T4 = 1400/80,2857 = 772,9 K

Praca sprężarki i turbiny

wc = 1,005·(534,4 − 295) = 240,6 kJ/kg
wt = 1,005·(1400 − 772,9) = 630,3 kJ/kg

Praca netto

wnet = 630,3 − 240,6 = 389,7 kJ/kg

Odpowiedź

wnet = 389,7 kJ/kg

Wskaźnik pracy zwrotnej: wc/wt = 38,17%.

25.Ciągłość i Bernoulli

Odpowiedź: p2 = 238,9 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·(100/75)² = 1,778 m/s
½ρ(v1² − v2²) = ½·1000·(1² − 1,778²)/1000 = −1,08 kPa
p2 = p1 + ½ρ(v1² − v2²) = 240 + (−1,08) = 238,9 kPa

Odpowiedź

p2 = 238,9 kPa
26.Ciągłość i Bernoulli

Odpowiedź: H = 5,123 m

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.

Prędkość wylotowa

β = 20/75 = 0,2667
v2 = √(2p1/(ρ(1 − β⁴))) = √(2·50000/(1000·0,9949)) = 10,03 m/s

Wysokość strumienia

Równanie Bernoulliego wzdłuż strumienia: w najwyższym punkcie prędkość jest równa zeru, a ciśnienie znów jest atmosferyczne.

H = v2²/(2g) = 10,03²/(2·9,81) = 5,123 m

Odpowiedź

H = 5,123 m
27.Tarcie w rurach i straty wysokości

Odpowiedź: hf = 29,6 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.

Wysokość prędkości

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

Darcy–Weisbach

hf = f·(L/D)·V²/(2g) = 0,03·(300/0,075)·2,2²/(2·9,81) = 29,6 m

Odpowiedź

hf = 29,6 m
28.Tarcie w rurach i straty wysokości

Odpowiedź: Q = 226,9 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,022·200/0,25 = 17,6
ΣK = 0,5 + 3·0,3 + 0,2 + 1,0 = 2,6

Równanie energii

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

22 = (17,6 + 2,6)·V²/(2·9,81)
V = √(2·9,81·22/20,2) = 4,623 m/s

Natężenie przepływu

Q = V·π·D²/4 = 4,623·0,04909 = 0,2269 m³/s = 226,9 L/s

Odpowiedź

Q = 226,9 L/s
29.Pompy i moc pomp

Odpowiedź: Pshaft = 11,08 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,03·32 = 9,418 kW

Moc na wale

Pshaft = Ph/η = 9,418/0,85 = 11,08 kW

Odpowiedź

Pshaft = 11,08 kW
30.Pompy i moc pomp

Odpowiedź: Pshaft = 14,26 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.

Prędkość

V = Q/A = 0,05/0,01767 = 2,829 m/s
V²/(2g) = 0,408 m

Straty wysokości

hL = (f·L/D + ΣK)·V²/(2g) = (24 + 5)·0,408 = 11,83 m

Wysokość podnoszenia pompy

H = Δz + hL = 12 + 11,83 = 23,83 m

Moc na wale

Pshaft = ρ·g·Q·H/η = 1000·9,81·0,05·23,83/0,82 = 14,26 kW

Odpowiedź

Pshaft = 14,26 kW
31.Przewodzenie ciepła

Odpowiedź: Ti = 7,411 °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 warstw

R1 = L1/(k1A) = (100/1000)/(0,04·4) = 0,625 K/W
R2 = L2/(k2A) = (150/1000)/(1,4·4) = 0,02679 K/W

Strumień ciepła

Q̇ = (T1 − T3)/(R1 + R2) = 10/0,6518 = 15,34 W

Temperatura na styku

Na pierwszej warstwie temperatura spada o Q̇·R₁.

Ti = T1 − Q̇·R1 = 17 − 15,34·0,625 = 7,411 °C

Odpowiedź

Ti = 7,411 °C
32.Przewodzenie ciepła

Odpowiedź: Ts = 18,54 °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/(8·4) = 0,03125 K/W
R1 = L1/(k1A) = (250/1000)/(1,4·4) = 0,04464 K/W
R2 = L2/(k2A) = (75/1000)/(0,035·4) = 0,5357 K/W
Ro = 1/(hoA) = 1/(25·4) = 0,01 K/W
ΣR = 0,6216 K/W

Straty ciepła

Q̇ = (20 − (-9))/0,6216 = 46,65 W
Ts = Troom − Q̇·Ri = 20 − 46,65·0,03125 = 18,54 °C

Odpowiedź

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

Odpowiedź: T2 = 70,99 N·m

Przekładnie zębate: zazębione koła mają tę samą prędkość na zębach, więc n₁·N₁ = n₂·N₂: koło o mniejszej liczbie zębów obraca się szybciej. Koła na jednym wale obracają się z tą samą prędkością obrotową. Moc jest przekazywana dalej (minus straty), więc moment obrotowy rośnie, gdy prędkość obrotowa maleje: T₂ = T₁·(N₂/N₁)·η, a P = T·ω, gdzie ω = 2π·n/60.

Moment silnika

ω₁ = 2π·1450/60 = 151,8 rad/s
T1 = P/ω₁ = 5500/151,8 = 36,22 N·m

Moment wyjściowy

T2 = 36,22·24/12·0,98 = 70,99 N·m

Odpowiedź

T2 = 70,99 N·m
34.Koła i przekładnie zębate

Odpowiedź: nout = 696 rpm

Przekładnie zębate: zazębione koła mają tę samą prędkość na zębach, więc n₁·N₁ = n₂·N₂: koło o mniejszej liczbie zębów obraca się szybciej. Koła na jednym wale obracają się z tą samą prędkością obrotową. Moc jest przekazywana dalej (minus straty), więc moment obrotowy rośnie, gdy prędkość obrotowa maleje: T₂ = T₁·(N₂/N₁)·η, a P = T·ω, gdzie ω = 2π·n/60.

Pierwszy stopień

nB = 2900·12/16 = 2175 rpm = nC

Drugi stopień

nout = 2175·16/50 = 696 rpm

Odpowiedź

nout = 696 rpm

przełożenie całkowite = 4,167 : 1

35.Transformaty Laplace’a

Odpowiedź: A = 3,6

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² + 9s + 14 = (s + 2)(s + 7)
F(s) = A/(s + 2) + B/(s + 7)

Metoda zakrywania

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

A = (2s + 22)/(s + 7) |s = −2 = 18/5 = 3,6

Odpowiedź

A = 3,6
(B = −1,6, więc f(t) = 3,6·e−2t − 1,6·e−7t)
36.Transformaty Laplace’a

Odpowiedź: y(0.1) = 2,066

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

(s² + 6s + 8)·Y(s) = 2s + 12 + 32/s
Y(s) = (2s² + 12s + 32)/(s(s + 2)(s + 4))

Ułamki proste

A = 32/(2·4) = 4
B = 16/((−2)·2) = −4
C = 16/((−4)·(−2)) = 2

Wynik

y(t) = 4 − 4·e−2t + 2·e−4t
y(0,1) = 2,066

Odpowiedź

y(0,1) = 2,066

Sprawdzenie: y(0) = A + B + C = 2.

37.Schematy blokowe

Odpowiedź: K = 1

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

1 + K/((s + 4)(s + 6)) = 0 → s² + 10s + 24 + K = 0

Pierwiastek podwójny

Trójmian s² + βs + γ ma pierwiastek podwójny, gdy β² = 4γ.

10² = 4(24 + K) → K = 100/4 − 24 = 1

Odpowiedź

K = 1
Oba bieguny układu zamkniętego leżą wtedy w s = −5.
38.Schematy blokowe

Odpowiedź: ωd = 5,268 rad/s

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 + 5) + 68·0,5 = 0 → s² + 5s + 34 = 0

Pierwiastki

s = −2,5 ± j·√(34 − 6,25) = −2,5 ± j5,268
ωd = 5,268 rad/s

Odpowiedź

ωd = 5,268 rad/s
39.Odpowiedź skokowa

Odpowiedź: ζ = 0,1572

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 + ωₙ².

ωₙ² = 91 → ωₙ = 9,539 rad/s
2ζωₙ = 3 → ζ = 3/(2·9,539) = 0,1572

Wynik

ζ = 0,1572

Odpowiedź

ζ = 0,1572
40.Odpowiedź skokowa

Odpowiedź: ζ = 0,5912

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).

Współczynnik tłumienia

ζ = −ln(0,1)/√(π² + ln²(0,1)) = 2,303/√(9,87 + 5,302) = 0,5912

Odpowiedź

ζ = 0,5912
41.Uchyb ustalony i regulacja PID

Odpowiedź: ess = 3,846 %

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 położeniowa

Kp = 5·5 = 25

Uchyb skokowy

ess = 100/(1 + 25) = 3,846 %

Odpowiedź

ess = 3,846 %
42.Uchyb ustalony i regulacja PID

Odpowiedź: %OS = 4,321 %

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

s(s + 10) + 5·10 = s² + 10s + 50

ωₙ i ζ

ωₙ = √50 = 7,071 rad/s
ζ = 10/(2·7,071) = 0,7071

Przeregulowanie

%OS = 100·e−0,7071π/√(1 − 0,7071²) = 4,321 %

Odpowiedź

%OS = 4,321 %

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