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 = 2 kN/mP = 25 kNABRA​ = ?4 m2 m4 m

Wspornik jest utwierdzony w ścianie w A. Przenosi obciążenie równomiernie rozłożone w od ściany do 6 m i siłę skupioną P w odległości 4 m od ściany. Wyznacz reakcję pionową RA w ścianie.

L = 10 mP = 25 kN w 4 mw = 2 kN/m od 0 do 6 m
RA = kN
2.Reakcje podpór belekTrudny
w = 8 kN/mP = 45 kNABCRA​ = ?1,5 m1 m1,5 m

Belka ma podporę nieprzesuwną w A i przesuwną w B, w odległości 2,5 m od A, i wystaje poza B aż do swobodnego końca C. Przenosi siłę skupioną P w odległości 1,5 m od A oraz obciążenie równomiernie rozłożone w od 1,5 m do swobodnego końca w 4 m. Wyznacz reakcję pionową RA, przyjmując zwrot w górę jako dodatni.

AB = 2,5 mAC = 4 mP = 45 kN w 1,5 mw = 8 kN/m od 1,5 do 4 m
RA = kN
3.Naprężenie i odkształcenie osioweŚredni
P = ?12300 mm500 mmØ20Ø30δ = 0,25 mm

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

odcinek 1: stal, Ø20 mm, 300 mm, E = 200 GPaodcinek 2: aluminium, Ø30 mm, 500 mm, E = 70 GPaδ = 0,25 mm
P = kN
4.Naprężenie i odkształcenie osioweTrudny
Ø 40 mmAluminium0,5 mØ 12 mmMiedź0,5 mP = 50 kN

Pręt schodkowy jest utwierdzony w ścianie i obciążony na swobodnym końcu. Wyznacz największe naprężenie normalne σmax w pręcie.

P = 50 kNodcinek 1: Ø40 mm, 0,5 m, aluminium (E = 70 GPa)odcinek 2: Ø12 mm, 0,5 m, miedź (E = 117 GPa)
σmax = MPa
5.Skręcanie wałówŚredni
T = 200 N·md = 25 mmL = 2 mφ = ?

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

T = 200 N·md = 25 mmL = 2 mG = 80 GPa
φ = °
6.Skręcanie wałówTrudny
Td = ?

Pełny wał ma przenosić 50 kW przy 1000 obr/min bez przekroczenia naprężenia stycznego 80 MPa. Wyznacz najmniejszą średnicę d, która wystarczy.

P = 50 kWN = 1000 rpmτallow = 80 MPa
d = mm
7.Siła tnąca i moment zginającyŚredni
8 kN/m45 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.

45 kN @ 6 m8 kN/m, 4–12 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
40 kN35 kN30 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.

40 kN @ 0 m35 kN @ 6 m30 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​ = −50σy​ = 70τxy​ = −20MPa

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

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

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

σx = −20 MPaσy = −40 MPaτxy = −20 MPaθ = 30°
σθ = MPa
11.Ugięcie belkiŚredni
P = 18 kNABδ5 m1 m

Belka swobodnie podparta przenosi siłę skupioną P w odległości a = 5 m od podpory A. Wyznacz ugięcie δP pod siłą.

P = 18 kNa = 5 mL = 6 mE = 11 GPa (drewno konstrukcyjne)I = 36 × 10⁶ mm⁴
δP = mm
12.Ugięcie belkiTrudny
w = 6 kN/mδL = 4 m

Belka swobodnie podparta przenosi obciążenie równomiernie rozłożone w. Jej ugięcie nie może przekroczyć rozpiętość/360. Wyznacz najmniejszy moment bezwładności Ireq, który spełnia ten warunek.

w = 6 kN/mL = 4 mE = 11 GPa (drewno konstrukcyjne)
Ireq = × 10⁶ mm⁴
13.Wyboczenie słupówŚredni
PL = 1 m25 × 80 mmprzekrój poprzeczny

Słup o długości L z materiału (aluminium) jest utwierdzony u podstawy i wolny u góry. Korzystając ze wzoru Eulera, wyznacz jego siłę krytyczną Pcr.

L = 1 mE = 70 GPa (aluminium)b × h = 25 × 80 mm
Pcr = kN
14.Wyboczenie słupówTrudny
PL = 2,5 md = 20 mmprzekrój poprzeczny

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

L = 2,5 mE = 200 GPa (stal)d = 20 mm
σcr = MPa
15.Rzut ukośnyŚredni
v0​ = ?70°R = 30 m

Piłkę wyrzucono z poziomego terenu pod kątem 70° do poziomu i upadła 30 m dalej. Pomiń opór powietrza i przyjmij g = 9,81 m/s². Wyznacz prędkość początkową v0.

θ = 70°R = 30 mg = 9,81 m/s²
v0 = m/s
16.Rzut ukośnyTrudny
h = 35 mv0​ = 27 m/s35°x = ?

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

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

Klocek zsuwa się bez prędkości początkowej po zboczu nachylonym pod kątem 40° do poziomu. Współczynnik tarcia kinetycznego wynosi 0,4. Wyznacz jego prędkość v po przebyciu 20 m wzdłuż zbocza.

θ = 40°μ = 0,4d = 20 mg = 9,81 m/s²
v = m/s
18.Praca i energiaTrudny
h = 11,3 mR = 2,5 mv = ?

Wózek o masie 10 kg zostaje puszczony bez prędkości początkowej z wysokości 11,3 m na gładkim torze, który przechodzi w pionową pętlę o promieniu 2,5 m. Wyznacz prędkość v wózka w najwyższym punkcie pętli.

m = 10 kgh = 11,3 mR = 2,5 mg = 9,81 m/s²
v = m/s
19.Drgania mechaniczneŚredni
k = 40 kN/mc = 310 N·s/mm = 15 kg

Masa 15 kg jest podparta na sprężynie o sztywności 40 kN/m i tłumiku wiskotycznym o współczynniku tłumienia 310 N·s/m. Wyznacz współczynnik tłumienia względnego ζ.

m = 15 kgk = 40 kN/mc = 310 N·s/m
ζ =
20.Drgania mechaniczneTrudny
k1​ = 30k2​ = 30m = 100 kg

Masa 100 kg wisi na dwóch sprężynach obok siebie (równolegle) o sztywnościach k₁ = 30 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 = 30 kN/mk2 = 30 kN/m
fn = Hz
21.Silniki cieplne i chłodziarkiŚredni
TH​ = 700 °CTL​ = 25 °CEQH​ = 2000 kJQL​W = ?

Silnik cieplny pracuje między zbiornikiem gorącym o temperaturze 700 °C a zbiornikiem zimnym o temperaturze 25 °C. Pobiera 2000 kJ ciepła. Wyznacz największą pracę Wmax, jaką mógłby wykonać.

TH = 700 °CTL = 25 °CQH = 2000 kJ
Wmax = kJ
22.Silniki cieplne i chłodziarkiTrudny
TH​ = 750 °CT0​ = 20 °CTL​ = -25 °CERWQ̇H​ = 30 kWQ̇L​ = ?

Silnik cieplny Carnota pobiera 30 kW ciepła ze źródła o temperaturze 750 °C i oddaje ciepło do otoczenia o temperaturze 20 °C. Cała jego praca napędza chłodziarkę Carnota, która odbiera ciepło z zamrażarki o temperaturze -25 °C i również oddaje ciepło do otoczenia o temperaturze 20 °C. Wyznacz strumień ciepła Q̇L odbieranego z zamrażarki.

TH = 750 °CT0 = 20 °CTL = -25 °CQ̇H = 30 kW
Q̇L = kW
23.Obiegi silników cieplnych gazowychŚredni
vP1234qin​qout​

Idealny obieg Diesla ma stopień sprężania r = 15 i stopień obciążenia rc = v₃/v₂ = 2,5. Na początku sprężania powietrze ma temperaturę 290 K. 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 temperaturę T3 na końcu doprowadzania ciepła, w K.

r = 15rc = 2,5T1 = 290 Kk = 1,4
T3 = K
24.Obiegi silników cieplnych gazowychTrudny
vP1234qin​qout​

Idealny obieg Braytona pracuje przy stosunku ciśnień w sprężarce 8. 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 = 8k = 1,4
η = %
25.Ciągłość i BernoulliŚredni
d1​ = 150 mmd2​ = 75 mm12p1​ − p2​ = 5 kPaQ = ?

Woda płynie przez poziomą zwężkę Venturiego o średnicy wlotowej 150 mm i średnicy gardzieli 75 mm. Ciśnienie na wlocie jest o 5 kPa wyższe niż w gardzieli. Pomijając straty, wyznacz natężenie przepływu Q.

d1 = 150 mmd2 = 75 mmp1 − p2 = 5 kPaρ = 1000 kg/m³
Q = L/s
26.Ciągłość i BernoulliTrudny
12d1​ = 75 mm, p1​ = 290 kPad2​ = 25 mm, p2​ = ?Δz = 10 m

Woda płynie w górę przez rurę, która zwęża się z 75 mm do 25 mm i wznosi się o 10 m. W przekroju 1 ciśnienie wynosi 290 kPa, a prędkość 1,25 m/s. Pomijając straty, wyznacz ciśnienie p2 w przekroju 2.

d1 = 75 mmd2 = 25 mmp1 = 290 kPav1 = 1,25 m/sΔz = 10 mρ = 1000 kg/m³
p2 = kPa
27.Tarcie w rurach i straty wysokościŚredni
V = 0,9 m/sD = 150 mmL = 300 mf = 0,025

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

V = 0,9 m/sL = 300 mD = 150 mmf = 0,025g = 9,81 m/s²
hf = m
28.Tarcie w rurach i straty wysokościTrudny
K = 5swobodny wypływH = ?L = 100 m, D = 200 mm, f = 0,0184× kolano (K = 0,3)

Rura o długości 100 m i średnicy wewnętrznej 200 mm (współczynnik tarcia Darcy’ego f = 0,018) wychodzi z otwartego zbiornika i wypływa swobodnym strumieniem do atmosfery. Straty miejscowe to: wlot o ostrej krawędzi (K = 0,5), kolana w liczbie 4 (K = 0,3 każde) i zawór (K = 5). Przyjmij g = 9,81 m/s². O ile H zwierciadło wody w zbiorniku musi leżeć wyżej niż wylot rury, aby natężenie przepływu wynosiło 80 L/s?

Q = 80 L/sL = 100 mD = 200 mmf = 0,018Kentrance = 0,5Kbend = 0,3 (×4)Kvalve = 5
H = m
29.Pompy i moc pompŚredni
PQ = 80 L/sH = 14 m

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

Q = 80 L/sH = 14 mηpump = 70%ηmotor = 90%g = 9,81 m/s²
Pel = kW
30.Pompy i moc pompTrudny
PQ = 60 L/sΔz = 6 mL = 200 m, D = 250 mm, f = 0,02, Σ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 6 m wyżej. Rura ma długość 200 m, średnicę wewnętrzną 250 mm i współczynnik tarcia Darcy’ego f = 0,02; współczynniki strat miejscowych sumują się do ΣK = 8 (łącznie z wlotem i wylotem). Przyjmij g = 9,81 m/s². Sprawność pompy wynosi 82%. Wyznacz moc na wale Pshaft w kW.

Q = 60 L/sΔz = 6 mL = 200 mD = 250 mmf = 0,02ΣK = 8η = 82%
Pshaft = kW
31.Przewodzenie ciepłaŚredni
Szkło, k = 0,96150 mmPłyta gipsowo-kartonowa, k = 0,1725 mm23 °C-12 °CQ̇ = ?

Ściana składa się z 150 mm (szkło) i 25 mm (płyta gipsowo-kartonowa), a jej pole wynosi 8 m². Powierzchnie zewnętrzne mają 23 °C i -12 °C. Wyznacz strumień ciepła Q̇ przez ścianę.

L1 = 150 mm, k1 = 0,96 W/(m·K)L2 = 25 mm, k2 = 0,17 W/(m·K)A = 8 m²T1 = 23 °CT3 = -12 °C
Q̇ = W
32.Przewodzenie ciepłaTrudny
Drewno, k = 0,17200 mmBeton, k = 1,4150 mmpomieszczenie 22 °Cna zewnątrz -12 °CQ̇

Ściana z 200 mm (drewno) i 150 mm (beton) oddziela pomieszczenie o temperaturze 22 °C od powietrza zewnętrznego o temperaturze -12 °C. Współczynniki przejmowania ciepła wynoszą 10 W/(m²·K) wewnątrz i 30 W/(m²·K) na zewnątrz, a pole ściany 6 m². Wyznacz współczynnik przenikania ciepła U ściany.

L1 = 200 mm, k1 = 0,17 W/(m·K)L2 = 150 mm, k2 = 1,4 W/(m·K)A = 6 m²hi = 10 W/(m²·K)ho = 30 W/(m²·K)Troom = 22 °CTout = -12 °C
U = W/(m²·K)
33.Koła i przekładnie zębateŚredni
d1​ = 80 mmd2​ = 300 mm960 rpmn2​ = ?

Koło pasowe 80 mm obracające się z prędkością 960 obr/min napędza przez pas koło 300 mm z poślizgiem 2 %. Wyznacz prędkość obrotową n2 koła napędzanego.

n1 = 960 rpmd1 = 80 mmd2 = 300 mms = 2 %
n2 = rpm
34.Koła i przekładnie zębateTrudny
A = 18B = 45C = 40D = 60

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

P = 7,5 kWnin = 1440 rpmA = 18B = 45C = 40D = 60ηmesh = 95%
Tout = N·m
35.Transformaty Laplace’aŚredni
ReIm−2−4−6−8−10−122

Znajdź f(t) dla F(s) = (3s + 32)/(s² + 4s + 8) metodą dopełnienia do kwadratu i oblicz f(0,2).

F(s) = (3s + 32)/(s² + 4s + 8)t = 0,2 s
f(0.2) =
36.Transformaty Laplace’aTrudny
t (s)y(t)t = 0,25y(0,25) = ?

Rozwiąż y″ + 7y′ + 12y = 0 dla y(0) = −4 i y′(0) = −8 za pomocą transformaty Laplace’a i znajdź y(0,25).

y″ + 7y′ + 12y = 0y(0) = −4y′(0) = −8t = 0,25 s
y(0.25) =
37.Schematy blokoweŚredni
R(s)+K(s + 6)(s + 11)G(s)C(s)−

Obiekt G(s) = K/((s + 6)(s + 11)) 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 + 6)(s + 11))jednostkowe ujemne sprzężenie zwrotne
K =
38.Schematy blokoweTrudny
R(s)+Ks(s + 1)(s + 6)G(s)C(s)−

Obiekt G(s) = K/(s(s + 1)(s + 6)) jest w pętli jednostkowego ujemnego sprzężenia zwrotnego. Przy największym stabilnym wzmocnieniu układ oscyluje. Znajdź częstość drgań ω w rad/s.

G(s) = K/(s(s + 1)(s + 6))jednostkowe ujemne sprzężenie zwrotne
ω = rad/s
39.Odpowiedź skokowaŚredni
t (s)y(t)1

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

G(s) = 52/(s² + 6s + 52)
ζ =
40.Odpowiedź skokowaTrudny
t (s)y(t)1tp​ = 0,8 s20 %

Odpowiedź skokowa układu drugiego rzędu bez zer wykazuje przeregulowanie 20 %, a pierwszy szczyt występuje w chwili t = 0,8 s. Znajdź pulsację własną ωn.

%OS = 20 %tp = 0,8 s
ωn = rad/s
41.Uchyb ustalony i regulacja PIDŚredni
R(s)+K = 941s + 1C(s)−

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

G(s) = 4/(1s + 1)K = 9
ess = %
42.Uchyb ustalony i regulacja PIDTrudny
R(s)+Kp​ + Ki​/s4s + 1C(s)−

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

G(s) = 4/(s + 1)ωₙ = 5 rad/sζ = 0,5
Kp =

3Roztwory

1.Reakcje podpór belek

Odpowiedź: RA = 37 kN

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 = 2·6 = 12 kN, x̄ = 3 m

Równowaga sił pionowych

RA = P + W = 25 + 12 = 37 kN

Odpowiedź

RA = 37 kN
2.Reakcje podpór belek

Odpowiedź: RA = 16 kN

Równowaga statyczna: ciało w spoczynku spełnia ΣFx = 0, ΣFy = 0 i ΣM = 0 względem dowolnego punktu. Zastąp każde obciążenie ciągłe jego wypadkową (polem wykresu obciążenia) przyłożoną w środku ciężkości wykresu obciążenia i licz momenty względem podpory, aby usunąć z równania reakcję tej podpory.

Wypadkowa obciążenia ciągłego

W = 8·2,5 = 20 kN, x̄ = (1,5 + 4)/2 = 2,75 m

Równowaga

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

ΣMA = 0: RB·2,5 = 45·1,5 + 20·2,75
RB = 122,5/2,5 = 49 kN
RA = 65 − 49 = 16 kN

Odpowiedź

RA = 16 kN
3.Naprężenie i odkształcenie osiowe

Odpowiedź: P = 16,8 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₂ = 706,9 mm²

Podatności

L₁/(E₁A₁) = 300/(200000·314,2) = 4,775 × 10⁻⁶ mm/N
L₂/(E₂A₂) = 500/(70000·706,9) = 1,011 × 10⁻⁵ mm/N

Obciążenie

P = 0,25/1,488 × 10⁻⁵ = 16801 N = 16,8 kN

Odpowiedź

P = 16,8 kN
4.Naprężenie i odkształcenie osiowe

Odpowiedź: σmax = 442,1 MPa

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 = π·40²/4 = 1257 mm², A2 = π·12²/4 = 113,1 mm²
σ1 = 39,79 MPa
σ2 = 442,1 MPa

Największe naprężenie

Cieńszy odcinek ma mniejsze pole, więc przenosi większe naprężenie.

σmax = 442,1 MPa

Odpowiedź

σmax = 442,1 MPa
5.Skręcanie wałów

Odpowiedź: φ = 7,47 °

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 = π·25⁴/32 = 38350 mm⁴

Kąt skręcenia

φ = TL/(GJ) = 200·1000·2000/(80000·38350) = 0,1304 rad
φ = 0,1304·180/π = 7,47°

Odpowiedź

φ = 7,47°
6.Skręcanie wałów

Odpowiedź: d = 31,21 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π·1000/60 = 104,7 rad/s
T = P/ω = 50000/104,7 = 477,5 N·m

Średnica z naprężenia dopuszczalnego

d = ∛(16T/(π·τallow)) = ∛(16·477465/(π·80)) = 31,21 mm

Odpowiedź

d = 31,21 mm

W praktyce zaokrągla się w górę do najbliższego wymiaru znormalizowanego.

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

Odpowiedź: D (x = 6 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 = 43,83 kN, R2 = 65,17 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 = 87,67, MC = 175,3, MD = 247 kN·m, ME = 196,7, MF = 114,3, MG = 0

Odpowiedź

Punkt D: MD = 247 kN·m
8.Siła tnąca i moment zginający

Odpowiedź: E (x = 4 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 = 50 kN, R2 = 55 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 = −40, MC = −30, MD = −50, ME = −70 kN·m, MF = −35, MG = 0

Odpowiedź

Punkt E: ME = −70 kN·m
9.Koło Mohra i transformacja naprężeń

Odpowiedź: τmax = 63,25 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 = √(((−50 − 70)/2)² + −20²) = 63,25 MPa

Odpowiedź

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

Odpowiedź: σθ = −42,32 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

σθ = −30 + 10·0,5 + −20·0,866 = −42,32 MPa

Odpowiedź

σθ = −42,32 MPa
11.Ugięcie belki

Odpowiedź: δP = 63,13 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 = 18000 N, a = 5000 mm, b = 1000 mm, L = 6000 mm
EI = 3,96 × 10¹¹ N·mm²

Ugięcie pod siłą

δP = Pa²b²/(3EIL) = 18000·5000²·1000²/(3·3,96 × 10¹¹·6000) = 63,13 mm

Odpowiedź

δP = 63,13 mm
12.Ugięcie belki

Odpowiedź: Ireq = 163,6 × 10⁶ 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).

Warunek ugięcia

δallow = L/360 = 4000/360 = 11,11 mm

Wymagany moment bezwładności

I = 5wL⁴/(384·E·δallow) = 5·6·4000⁴/(384·11000·11,11) = 1,636 × 10⁸ mm⁴
Ireq = 163,6 × 10⁶ mm⁴

Odpowiedź

Ireq = 163,6 × 10⁶ mm⁴
13.Wyboczenie słupów

Odpowiedź: Pcr = 17,99 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 = 80·25³/12 = 104167 mm⁴

Długość wyboczeniowa

KL = 2·1000 = 2000 mm

Siła krytyczna

Pcr = π²EI/(KL)² = π²·70000·104167/2000² = 17991 N = 17,99 kN

Odpowiedź

Pcr = 17,99 kN
14.Wyboczenie słupów

Odpowiedź: σcr = 16,11 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 = π·20⁴/64 = 7854 mm⁴

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

A = 314,2 mm²
r = √(I/A) = √(7854/314,2) = 5 mm
KL = 0,7·2500 = 1750 mm
KL/r = 1750/5 = 350

Naprężenie krytyczne

σcr = π²E/(KL/r)² = π²·200000/350² = 16,11 MPa

Odpowiedź

σcr = 16,11 MPa
15.Rzut ukośny

Odpowiedź: v0 = 21,4 m/s

Projectile motion: with no air resistance, x = vxt and y = h + vyt − ½gt². The two directions are independent except that they share the same time.

Wyznacz v₀ ze wzoru na zasięg

v₀ = √(g·R/sin 2θ) = √(9,81·30/sin 140°) = 21,4 m/s

Odpowiedź

v0 = 21,4 m/s
16.Rzut ukośny

Odpowiedź: x = 103,5 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 θ = 27·cos 35° = 22,12 m/s
vy = v0·sin θ = 27·sin 35° = 15,49 m/s

Czas lotu

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

0 = 35 + 15,49·t − 4,905·t²
t = (vy + √(vy² + 2gh))/g = (15,49 + √(239,8 + 686,7))/9,81 = 4,682 s

Odległość pozioma

x = vx·t = 22,12·4,682 = 103,5 m

Odpowiedź

x = 103,5 m
17.Praca i energia

Odpowiedź: v = 11,49 m/s

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

Bilans energii wzdłuż zbocza

Klocek obniża się o d·sin θ, a tarcie zabiera μmg·cos θ·d. Masa się skraca.

½v² = g·d·(sin θ − μ·cos θ)
v = √(2·9,81·20·(0,6428 − 0,4·0,766)) = 11,49 m/s

Odpowiedź

v = 11,49 m/s
18.Praca i energia

Odpowiedź: v = 11,12 m/s

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). Na torze kołowym siła wypadkowa skierowana do środka musi być równa mv²/r.

Prędkość w najwyższym punkcie

v² = 2g(h − 2R) = 2·9,81·(11,3 − 5) = 123,6 m²/s²
v = 11,12 m/s

Odpowiedź

v = 11,12 m/s

Wózek ledwo pokonuje pętlę, gdy h = 2,5R; wtedy N = 0 w najwyższym punkcie.

19.Drgania mechaniczne

Odpowiedź: ζ = 0,2001

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

Tłumienie krytyczne

cc = 2·√(k·m) = 2·√(40000·15) = 1549 N·s/m

Współczynnik tłumienia

ζ = c/cc = 310/1549 = 0,2001

Odpowiedź

ζ = 0,2001

ζ < 1, więc tłumienie jest podkrytyczne, a układ drga z malejącą amplitudą.

20.Drgania mechaniczne

Odpowiedź: fn = 3,898 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 = 30 + 30 = 60 kN/m

Częstotliwość drgań własnych

ωn = √(60000/100) = 24,49 rad/s
fn = 24,49/(2π) = 3,898 Hz

Odpowiedź

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

Odpowiedź: Wmax = 1387 kJ

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 = 700 + 273,15 = 973,1 K
TL = 25 + 273,15 = 298,1 K

Sprawność Carnota

ηmax = 1 − TL/TH = 1 − 298,1/973,1 = 69,36 %

Maksymalna praca

Wmax = ηmax·QH = 0,6936·2000 = 1387 kJ

Odpowiedź

Wmax = 1387 kJ
22.Silniki cieplne i chłodziarki

Odpowiedź: Q̇L = 118 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 = 750 + 273,15 = 1023 K
T0 = 20 + 273,15 = 293,1 K
TL = -25 + 273,15 = 248,1 K

Praca silnika

η = 1 − 293,1/1023 = 0,7135
Ẇ = 0,7135·30 = 21,4 kW

Ciepło odbierane przez chłodziarkę

COPR = 248,1/(293,1 − 248,1) = 5,514
Q̇L = 5,514·21,4 = 118 kW

Odpowiedź

Q̇L = 118 kW
23.Obiegi silników cieplnych gazowych

Odpowiedź: T3 = 2142 K

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

Koniec sprężania

T2 = 290·150,4 = 856,7 K

Doprowadzenie ciepła przy stałym ciśnieniu

T3 = rc·T2 = 2,5·856,7 = 2142 K

Odpowiedź

T3 = 2142 K
24.Obiegi silników cieplnych gazowych

Odpowiedź: η = 44,8 %

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/80,2857 = 1 − 1/1,811 = 44,8%

Odpowiedź

η = 44,8%
25.Ciągłość i Bernoulli

Odpowiedź: Q = 14,43 L/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.

Ciągłość

v1 = v2·(75/150)² = 0,25·v2

Równanie Bernoulliego między wlotem a gardzielą

Rura jest pozioma, więc zmieniają się tylko ciśnienie i prędkość.

p1 − p2 = ½ρv2²(1 − β⁴)
v2 = √(2·5000/(1000·(1 − 0,0625))) = 3,266 m/s
v1 = 0,8165 m/s

Natężenie przepływu

Q = A2·v2 = 0,004418·3,266 = 0,01443 m³/s = 14,43 L/s

Odpowiedź

Q = 14,43 L/s
26.Ciągłość i Bernoulli

Odpowiedź: p2 = 129,4 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
ρgΔz = 1000·9,81·10/1000 = 98,1 kPa
p2 = p1 + ½ρ(v1² − v2²) − ρgΔz = 290 + (−62,5) − 98,1 = 129,4 kPa

Odpowiedź

p2 = 129,4 kPa
27.Tarcie w rurach i straty wysokości

Odpowiedź: hf = 2,064 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) = 0,9²/(2·9,81) = 0,04128 m

Darcy–Weisbach

hf = f·(L/D)·V²/(2g) = 0,025·(300/0,15)·0,9²/(2·9,81) = 2,064 m

Odpowiedź

hf = 2,064 m
28.Tarcie w rurach i straty wysokości

Odpowiedź: H = 5,519 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ść

V = Q/A = 0,08/0,03142 = 2,546 m/s
V²/(2g) = 0,3305 m

Współczynniki strat

f·L/D = 0,018·100/0,2 = 9
ΣK = 0,5 + 4·0,3 + 5 = 6,7

Równanie energii

Strumień na wylocie nadal niesie swoją wysokość prędkości; to nie strata, ale trzeba ją dostarczyć.

H = (1 + 9 + 6,7)·0,3305 = 5,519 m

Odpowiedź

H = 5,519 m
29.Pompy i moc pomp

Odpowiedź: Pel = 17,44 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,08·14 = 10,99 kW

Moc na wale

Pshaft = 10,99/0,7 = 15,7 kW

Moc elektryczna

Pel = 15,7/0,9 = 17,44 kW

Odpowiedź

Pel = 17,44 kW
30.Pompy i moc pomp

Odpowiedź: Pshaft = 5,619 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,06/0,04909 = 1,222 m/s
V²/(2g) = 0,07615 m

Straty wysokości

hL = (f·L/D + ΣK)·V²/(2g) = (16 + 8)·0,07615 = 1,828 m

Wysokość podnoszenia pompy

H = Δz + hL = 6 + 1,828 = 7,828 m

Moc na wale

Pshaft = ρ·g·Q·H/η = 1000·9,81·0,06·7,828/0,82 = 5,619 kW

Odpowiedź

Pshaft = 5,619 kW
31.Przewodzenie ciepła

Odpowiedź: Q̇ = 923,2 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) = (150/1000)/(0,96·8) = 0,01953 K/W
R2 = L2/(k2A) = (25/1000)/(0,17·8) = 0,01838 K/W

Strumień ciepła

Q̇ = (T1 − T3)/(R1 + R2) = 35/0,03791 = 923,2 W

Odpowiedź

Q̇ = 923,2 W
32.Przewodzenie ciepła

Odpowiedź: U = 0,7057 W/(m²·K)

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/(10·6) = 0,01667 K/W
R1 = L1/(k1A) = (200/1000)/(0,17·6) = 0,1961 K/W
R2 = L2/(k2A) = (150/1000)/(1,4·6) = 0,01786 K/W
Ro = 1/(hoA) = 1/(30·6) = 0,005556 K/W
ΣR = 0,2362 K/W

Współczynnik przenikania

U = 1/(A·ΣR) = 1/(6·0,2362) = 0,7057 W/(m²·K)

Odpowiedź

U = 0,7057 W/(m²·K)
33.Koła i przekładnie zębate

Odpowiedź: n2 = 250,9 rpm

Przekładnie pasowe: bez poślizgu oba koła pasowe mają tę samą prędkość obwodową v = π·d·n/60, więc n₁·d₁ = n₂·d₂, a przełożenie wynosi d₂/d₁. Poślizg s sprawia, że koło napędzane obraca się nieco wolniej: n₂ = n₁·(d₁/d₂)·(1 − s). Przenoszona moc to P = (T₁ − T₂)·v, gdzie T₁ to siła w cięgnie czynnym, a T₂ w biernym. Tarcie ogranicza stosunek sił: T₁/T₂ = eμθ, gdzie θ jest kątem opasania w radianach na mniejszym kole; dla pasa otwartego θ = π − 2·asin((d₂ − d₁)/(2C)) przy rozstawie osi C.

Prędkość bez poślizgu

960·80/300 = 256 rpm

Z poślizgiem

n2 = 256·0,98 = 250,9 rpm

Odpowiedź

n2 = 250,9 rpm
34.Koła i przekładnie zębate

Odpowiedź: Tout = 168,3 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 = 1440·(18/45)·(40/60) = 384 rpm
ωout = 2π·384/60 = 40,21 rad/s

Moc i moment wyjściowy

Pout = 7500·0,95² = 6769 W
Tout = Pout/ωout = 6769/40,21 = 168,3 N·m

Odpowiedź

Tout = 168,3 N·m
35.Transformaty Laplace’a

Odpowiedź: f(0.2) = 5,246

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

Dopełnienie do kwadratu

s² + 4s + 8 = (s + 2)² + 4 = (s + 2)² + 2²

Dopasuj licznik

F(s) = 3·(s + 2)/((s + 2)² + 2²) + (26/2)·2/((s + 2)² + 2²)

Transformata odwrotna

f(t) = e−2t·(3·cos 2t + 13·sin 2t)
f(0,2) = e−2·0,2·(3·cos(2·0,2) + 13·sin(2·0,2)) = 5,246

Odpowiedź

f(0,2) = 5,246
36.Transformaty Laplace’a

Odpowiedź: y(0.25) = −3,979

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²Y − (−4)s − (−8) + 7(sY − (−4)) + 12Y = 0
Y(s) = (−4s − 36)/((s + 3)(s + 4))

Ułamki proste

A = −24/1 = −24
B = −20/(−1) = 20

Wynik

y(t) = −24·e−3t + 20·e−4t
y(0,25) = −3,979

Odpowiedź

y(0,25) = −3,979

Sprawdzenie: y(0) = A + B = −4.

37.Schematy blokowe

Odpowiedź: K = 6,25

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 + 6)(s + 11)) = 0 → s² + 17s + 66 + K = 0

Pierwiastek podwójny

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

17² = 4(66 + K) → K = 289/4 − 66 = 6,25

Odpowiedź

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

Odpowiedź: ω = 2,449 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³ + 7s² + 6s + K = 0

Tablica Routha

s³ | 1 6
s² | 7 K
s¹ | (7·6 − K)/7
s⁰ | K

Granica stabilności

Wszystkie elementy pierwszej kolumny muszą być dodatnie.

Kmax = 7·6 = 42

Częstość drgań

7·(jω)² + 42 = 0 → ω² = 42/7 = 6
ω = 2,449 rad/s

Odpowiedź

ω = 2,449 rad/s
39.Odpowiedź skokowa

Odpowiedź: ζ = 0,416

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

ωₙ² = 52 → ωₙ = 7,211 rad/s
2ζωₙ = 6 → ζ = 6/(2·7,211) = 0,416

Wynik

ζ = 0,416

Odpowiedź

ζ = 0,416
40.Odpowiedź skokowa

Odpowiedź: ωn = 4,412 rad/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).

Współczynnik tłumienia

ζ = −ln(0,2)/√(π² + ln²(0,2)) = 1,609/√(9,87 + 2,59) = 0,4559

Częstotliwość drgań własnych

ωd = π/tp = π/0,8 = 3,927 rad/s
ωₙ = ωd/√(1 − ζ²) = 3,927/0,89 = 4,412 rad/s

Odpowiedź

ωn = 4,412 rad/s
41.Uchyb ustalony i regulacja PID

Odpowiedź: ess = 2,703 %

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 = 9·4 = 36

Uchyb skokowy

ess = 100/(1 + 36) = 2,703 %

Odpowiedź

ess = 2,703 %
42.Uchyb ustalony i regulacja PID

Odpowiedź: Kp = 1

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 + 4Kp)s + 4Ki

Wielomian docelowy

s² + 2·0,5·5s + 5² = s² + 5s + 25

Dopasuj współczynniki

1 + 4Kp = 5 → Kp = 1

Odpowiedź

Kp = 1

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