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
P = 55 kNM = 50 kN·mABRA​ = ?2 m3 m3 m

Belka swobodnie podparta ma podporę nieprzesuwną w A i przesuwną w B. Przenosi siłę skupioną P w odległości 2 m od A i moment skupiony M działający przeciwnie do ruchu wskazówek zegara, przyłożony w odległości 5 m od A. Wyznacz reakcję pionową RA, przyjmując zwrot w górę jako dodatni.

L = 8 mP = 55 kN w 2 mM = 50 kN·m (przeciwnie do ruchu wskazówek zegara) w 5 m
RA = kN
2.Reakcje podpór belekTrudny
w = 2 kN/mP = 60 kNABCRB​ = ?1,5 m1 m1,5 m1 m

Belka ma podporę nieprzesuwną w A i przesuwną w B, w odległości 4 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 2,5 m do swobodnego końca w 5 m. Wyznacz reakcję pionową RB, przyjmując zwrot w górę jako dodatni.

AB = 4 mAC = 5 mP = 60 kN w 1,5 mw = 2 kN/m od 2,5 do 5 m
RB = kN
3.Naprężenie i odkształcenie osioweŚredni
Ø 12 mmMosiądz1,5 mP = 85 kN

Wyznacz wydłużenie δ pręta (mosiądz) pod obciążeniem osiowym.

P = 85 kNd = 12 mmL = 1,5 mE = 100 GPa
δ = mm
4.Naprężenie i odkształcenie osioweTrudny
Ø 25 mmAluminium1 mØ 20 mmAluminium1 mP = 60 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 = 60 kNodcinek 1: Ø25 mm, 1 m, aluminium (E = 70 GPa)odcinek 2: Ø20 mm, 1 m, aluminium (E = 70 GPa)
σmax = MPa
5.Skręcanie wałówŚredni
d = ?T = 2000 N·mL = 0,5 m

Pełny wał (aluminium) o długości 0,5 m przenosi moment skręcający 2000 N·m. Kąt skręcenia między jego końcami nie może przekroczyć 1°. Wyznacz najmniejszą średnicę d, która spełnia ten warunek.

T = 2000 N·mL = 0,5 mG = 26 GPaφallow = 1°
d = mm
6.Skręcanie wałówTrudny
do​ = ?20 kW, 1000 rpm

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

P = 20 kWN = 1000 rpmdi/do = 0,8τallow = 40 MPa
do = mm
7.Siła tnąca i moment zginającyŚredni
10 kN35 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.

10 kN @ 2 m35 kN @ 10 m
  1. 1A (x = 0 m)
  2. 2B (x = 2 m)
  3. 3C (x = 4 m)
  4. 4D (x = 6 m)
  5. 5E (x = 8 m)
  6. 6F (x = 10 m)
  7. 7G (x = 12 m)
8.Siła tnąca i moment zginającyTrudny
w = 3 kN/mP = 15 kNAB7 m1 m

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

L = 8 mP = 15 kN w 7 mw = 3 kN/m od 0 do 8 m
Mmax = kN·m
9.Koło Mohra i transformacja naprężeńŚredni
σx​ = 20σy​ = −20τxy​ = 60MPa

W pewnym punkcie obciążonej części występują σₓ = 20 MPa, σᵧ = −20 MPa i τₓᵧ = 60 MPa (rozciąganie dodatnie). Wyznacz kąt θp od osi x do kierunku największego naprężenia głównego, z przedziału od −90° do 90°.

σx = 20 MPaσy = −20 MPaτxy = 60 MPa
θp = °
10.Koło Mohra i transformacja naprężeńTrudny
σx​τxy​górna powierzchnia wału

Pełny wał o średnicy 20 mm przenosi moment zginający 200 N·m wraz z momentem skręcającym 400 N·m. Na górnej powierzchni, gdzie naprężenie od zginania jest rozciągające i największe, wyznacz największe naprężenie główne σ1.

d = 20 mmM = 200 N·mT = 400 N·m
σ1 = MPa
11.Ugięcie belkiŚredni
w = 8 kN/mδL = 1,5 m

Wspornik przenosi obciążenie równomiernie rozłożone w na całej długości. Wyznacz ugięcie δ swobodnego końca.

w = 8 kN/mL = 1,5 mE = 11 GPa (drewno konstrukcyjne)I = 36 × 10⁶ mm⁴
δ = mm
12.Ugięcie belkiTrudny
w = 4 kN/mABδ6 m

Belka drewniana swobodnie podparta o przekroju prostokątnym i szerokości 150 mm przenosi obciążenie równomiernie rozłożone w na całej rozpiętości. Jej ugięcie w środku rozpiętości nie może przekroczyć rozpiętość/300. Wyznacz najmniejszą wysokość hreq przekroju.

w = 4 kN/mL = 6 mE = 11 GPa (drewno konstrukcyjne)b = 150 mm
hreq = mm
13.Wyboczenie słupówŚredni
PL = 1,5 md = 40 mmprzekrój poprzeczny

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

L = 1,5 mE = 70 GPa (aluminium)d = 40 mm
Pcr = kN
14.Wyboczenie słupówTrudny
P = 40 kNL = 4 md = ?przekrój poprzeczny

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

P = 40 kNn = 2L = 4 mE = 70 GPa (aluminium)
d = mm
15.Rzut ukośnyŚredni
25 my = ?v0​ = 36 m/s55°

Piłkę kopnięto z poziomego terenu z prędkością 36 m/s pod kątem 55° do poziomu w stronę muru odległego o 25 m. Pomiń opór powietrza i przyjmij g = 9,81 m/s². Na jakiej wysokości y piłka mija mur?

v₀ = 36 m/sθ = 55°x = 25 mg = 9,81 m/s²
y = m
16.Rzut ukośnyTrudny
h = 40 mv0​ = 22 m/s50°zasięg

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

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

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

θ = 45°μ = 0,35d = 14 mg = 9,81 m/s²
v = m/s
18.Praca i energiaTrudny
h = 2,5 mL = 4,5 mμ = 0,15k = 1500 N/m

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

m = 0,5 kgh = 2,5 mL = 4,5 mμ = 0,15k = 1500 N/mg = 9,81 m/s²
x = cm
19.Drgania mechaniczneŚredni
k = 200 kN/mζ = 0,7m = 5 kg

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

m = 5 kgk = 200 kN/mζ = 0,7
fd = Hz
20.Drgania mechaniczneTrudny
k1​ = 10k2​ = 80m = 40 kg

Masa 40 kg wisi na dwóch sprężynach obok siebie (równolegle) o sztywnościach k₁ = 10 kN/m i k₂ = 80 kN/m. Masa pozostaje pozioma. Wyznacz sztywność zastępczą keq układu sprężyn w kN/m.

m = 40 kgk1 = 10 kN/mk2 = 80 kN/m
keq = kN/m
21.Silniki cieplne i chłodziarkiŚredni
TH​ = 28 °CTL​ = 0 °CRQ̇H​ = ?Q̇L​ = 2 kWẆ

Odwracalna chłodziarka utrzymuje w chłodni temperaturę 0 °C, odbierając 2 kW ciepła, i oddaje ciepło do otoczenia o temperaturze 28 °C. Wyznacz strumień ciepła oddawanego Q̇H.

TL = 0 °CTH = 28 °CQ̇L = 2 kW
Q̇H = kW
22.Silniki cieplne i chłodziarkiTrudny
TH​ = 300 °CTL​ = 10 °CEQH​ = 2500 kJQL​ = 1625 kJW = 875 kJ

Silnik cieplny w każdym cyklu pobiera 2500 kJ ciepła ze zbiornika o temperaturze 300 °C, wykonuje 875 kJ pracy, a resztę oddaje do zbiornika o temperaturze 10 °C. Wyznacz całkowitą entropię generowaną Sgen w jednym cyklu.

TH = 300 °CTL = 10 °CQH = 2500 kJW = 875 kJ
Sgen = kJ/K
23.Obiegi silników cieplnych gazowychŚredni
vP1234qin​qout​

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

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

Idealny obieg Braytona pracuje przy stosunku ciśnień w sprężarce 20. 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 = 20k = 1,4
η = %
25.Ciągłość i BernoulliŚredni
12d1​ = 150 mm, p1​ = 150 kPad2​ = 50 mm, p2​ = ?

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

d1 = 150 mmd2 = 50 mmp1 = 150 kPav1 = 1,25 m/sρ = 1000 kg/m³
p2 = kPa
26.Ciągłość i BernoulliTrudny
12d1​ = 150 mm, p1​ = 180 kPad2​ = 100 mm, p2​ = ?Δz = 9 m

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

d1 = 150 mmd2 = 100 mmp1 = 180 kPav1 = 1 m/sΔz = 9 mρ = 1000 kg/m³
p2 = kPa
27.Tarcie w rurach i straty wysokościŚredni
Q = 120 L/sD = 300 mmL = 800 mf = 0,02

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

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

Rura o długości 150 m i średnicy wewnętrznej 150 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 2 (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 25 L/s?

Q = 25 L/sL = 150 mD = 150 mmf = 0,018Kentrance = 0,5Kbend = 0,3 (×2)Kvalve = 5
H = m
29.Pompy i moc pompŚredni
PQ = 15 L/sH = 38 m

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

Q = 15 L/sH = 38 mη = 70%g = 9,81 m/s²
Pshaft = kW
30.Pompy i moc pompTrudny
200 kPa (nadciśnienie)PQ = 20 L/sΔz = 8 mL = 100 m, D = 100 mm, f = 0,025, ΣK = 4

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

Q = 20 L/sΔz = 8 mL = 100 mD = 100 mmf = 0,025ΣK = 4ptank = 200 kPa
H = m
31.Przewodzenie ciepłaŚredni
Płyta gipsowo-kartonowa, k = 0,17150 mmDrewno, k = 0,1775 mm22 °C-3 °CQ̇ = ?

Ściana składa się z 150 mm (płyta gipsowo-kartonowa) i 75 mm (drewno), a jej pole wynosi 6 m². Powierzchnie zewnętrzne mają 22 °C i -3 °C. Wyznacz strumień ciepła Q̇ przez ścianę.

L1 = 150 mm, k1 = 0,17 W/(m·K)L2 = 75 mm, k2 = 0,17 W/(m·K)A = 6 m²T1 = 22 °CT3 = -3 °C
Q̇ = W
32.Przewodzenie ciepłaTrudny
Drewno, k = 0,1750 mmPłyta gipsowo-kartonowa, k = 0,1725 mmpomieszczenie 21 °Cna zewnątrz -12 °CQ̇

Ściana z 50 mm (drewno) i 25 mm (płyta gipsowo-kartonowa) oddziela pomieszczenie o temperaturze 21 °C od powietrza zewnętrznego o temperaturze -12 °C. Współczynniki przejmowania ciepła wynoszą 10 W/(m²·K) wewnątrz i 20 W/(m²·K) na zewnątrz, a pole ściany 20 m². Wyznacz współczynnik przenikania ciepła U ściany.

L1 = 50 mm, k1 = 0,17 W/(m·K)L2 = 25 mm, k2 = 0,17 W/(m·K)A = 20 m²hi = 10 W/(m²·K)ho = 20 W/(m²·K)Troom = 21 °CTout = -12 °C
U = W/(m²·K)
33.Koła i przekładnie zębateŚredni
Nw​ = 45z = 3, 1450 rpm

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

z = 3Nw = 45n = 1450 rpm
nw = rpm
34.Koła i przekładnie zębateTrudny
d1​ = 160 mmd2​ = 315 mm1440 rpmC = 800 mmθ = ?

Pas otwarty łączy koło 160 mm i koło 315 mm, których środki są oddalone o 800 mm. Wyznacz kąt opasania θ na mniejszym kole, w stopniach.

d1 = 160 mmd2 = 315 mmC = 800 mm
θ = °
35.Transformaty Laplace’aŚredni
t (s)f(t)t = 0,25f(0,25) = ?

Znajdź f(t) dla F(s) = 55/(s(s + 4)(s + 6)) metodą rozkładu na ułamki proste i oblicz f(0,25).

F(s) = 55/(s(s + 4)(s + 6))t = 0,25 s
f(0.25) =
36.Transformaty Laplace’aTrudny
t (s)y(t)t = 0,6y(0,6) = ?

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

y″ + 6y′ + 8y = 16·u(t)u(t) = 5y(0) = 0, y′(0) = 0t = 0,6 s
y(0.6) =
37.Schematy blokoweŚredni
R(s)+10(s + 8)(s + 2)(s + 7)G(s)C(s)−

Pętla jednostkowego ujemnego sprzężenia zwrotnego ma G(s) = 10(s + 8)/((s + 2)(s + 7)). Na wartość zadaną podano skok o wielkości 1. Znajdź uchyb ustalony ess.

G(s) = 10(s + 8)/((s + 2)(s + 7))R = 1
ess =
38.Schematy blokoweTrudny
R(s)+Kregulator+D(s)+3s + 10obiektC(s)−

Regulator proporcjonalny K steruje obiektem G(s) = 3/(s + 10) w pętli jednostkowego ujemnego sprzężenia zwrotnego. Skok zakłócenia o wielkości 10 wchodzi na wejście obiektu, gdy R = 0. Znajdź najmniejsze wzmocnienie K, przy którym wyjście w stanie ustalonym nie przekracza 0,15.

G(s) = 3/(s + 10)D = 10css ≤ 0,15
K =
39.Odpowiedź skokowaŚredni
t (s)y(t)2ts​ = ?

Układ ma transmitancję G(s) = 104/(s² + 2s + 52) i otrzymuje jednostkowe wymuszenie skokowe. Znajdź 2-procentowy czas ustalania ts odpowiedzi skokowej, korzystając z ts ≈ 4/(ζωₙ).

G(s) = 104/(s² + 2s + 52)
ts = s
40.Odpowiedź skokowaTrudny
k = 5000 N/mc = 63 N·s/mm = 5 kgF = 10 Nx(t)

Masa m = 5 kg wisi na sprężynie k = 5000 N/m z tłumikiem lepkim c = 63 N·s/m (m·x″ + c·x′ + k·x = F). W spoczynku w położeniu równowagi nagle przyłożono stałą siłę F = 10 N. Znajdź największe wychylenie xmax w mm.

m = 5 kgk = 5000 N/mc = 63 N·s/mF = 10 N
xmax = mm
41.Uchyb ustalony i regulacja PIDŚredni
r(t) = 2·t²/2+2(s + 1)s²(s + 1)C(s)−

Na pętlę ze sprzężeniem jednostkowym o G(s) = 2(s + 1)/(s²(s + 1)) podano wejście paraboliczne r(t) = 2·t²/2. Wyznacz uchyb ustalony ess.

G(s) = 2(s + 1)/(s²(s + 1))r(t) = 2·t²/2
ess =
42.Uchyb ustalony i regulacja PIDTrudny
R(s)+K = 252s(s + 10)C(s)−

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

G(s) = 2/(s(s + 10))K = 25
%OS = %

3Roztwory

1.Reakcje podpór belek

Odpowiedź: RA = 47,5 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.

Momenty względem A

Przeciwnie do ruchu wskazówek zegara dodatnio.

ΣMA = 0: RB·8 − 55·2 + 50 = 0
RB = 60/8 = 7,5 kN

Równowaga sił pionowych

RA = 55 − 7,5 = 47,5 kN

Odpowiedź

RA = 47,5 kN

Ujemna reakcja działa w dół.

2.Reakcje podpór belek

Odpowiedź: RB = 27,19 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 = 2·2,5 = 5 kN, x̄ = (2,5 + 5)/2 = 3,75 m

Równowaga

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

ΣMA = 0: RB·4 = 60·1,5 + 5·3,75
RB = 108,8/4 = 27,19 kN
RA = 65 − 27,19 = 37,81 kN

Odpowiedź

RB = 27,19 kN
3.Naprężenie i odkształcenie osiowe

Odpowiedź: δ = 11,27 mm

Prawo Hooke’a: w zakresie sprężystym σ = E·ε. Łącząc to z σ = P/A i ε = δ/L, dostajemy δ = PL/(AE).

Naprężenie

A = π·12²/4 = 113,1 mm²
σ = P/A = 85·1000/113,1 = 751,6 MPa

Odkształcenie (prawo Hooke’a)

E for mosiądz = 100 GPa = 100000 MPa.

ε = σ/E = 751,6/100000 = 0,007516 = 7516 µε

Wydłużenie

δ = ε·L = 0,007516·1500 mm = 11,27 mm

Odpowiedź

δ = 11,27 mm
4.Naprężenie i odkształcenie osiowe

Odpowiedź: σmax = 191 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 = π·25²/4 = 490,9 mm², A2 = π·20²/4 = 314,2 mm²
σ1 = 122,2 MPa
σ2 = 191 MPa

Największe naprężenie

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

σmax = 191 MPa

Odpowiedź

σmax = 191 MPa
5.Skręcanie wałów

Odpowiedź: d = 68,83 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·ω.

Wymagany biegunowy moment bezwładności

T = 2000 N·m = 2 × 10⁶ N·mm
L = 0,5 m = 500 mm
G = 26 GPa = 26000 MPa
φ = 1·π/180 = 0,01745 rad
J = TL/(Gφ) = 2 × 10⁶·500/(26000·0,01745) = 2,204 × 10⁶ mm⁴

Średnica

d = (32J/π)1/4 = 68,83 mm

Odpowiedź

d = 68,83 mm

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

6.Skręcanie wałów

Odpowiedź: do = 34,53 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/ω = 20000/104,7 = 191 N·m

Średnica zewnętrzna z naprężenia dopuszczalnego

T = 191 N·m = 190986 N·mm
1 − k⁴ = 1 − 0,8⁴ = 0,5904
do = ∛(16T/(π·τallow·(1 − k⁴))) = ∛(16·190986/(π·40·0,5904)) = 34,53 mm
di = 0,8·34,53 = 27,63 mm

Odpowiedź

do = 34,53 mm

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

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

Odpowiedź: F (x = 10 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 = 14,17 kN, R2 = 30,83 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 = 28,33, MC = 36,67, MD = 45, ME = 53,33, MF = 61,67 kN·m, MG = 0

Odpowiedź

Punkt F: MF = 61,67 kN·m
8.Siła tnąca i moment zginający

Odpowiedź: Mmax = 32,09 kN·m

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

Reakcje podpór

Równowaga całej belki.

Wypadkowa obciążenia ciągłego: W = w·ℓ = 3·8 = 24 kN, przyłożona w x̄ = 4 m
ΣMA = 0: RB·8 = 15·7 + 3·8·4
RB = 201/8 = 25,13 kN
ΣFy = 0: RA = 39 − 25,13 = 13,88 kN

Gdzie V = 0?

RA − w·x0 = 0 → x0 = 13,88/3 = 4,625 m

Moment maksymalny

Mmax = RA·x0 − w·x0²/2 = 13,88·4,625 − 3·4,625²/2 = 32,09 kN·m

Odpowiedź

Mmax = 32,09 kN·m

Występuje w x = 4,625 m od A.

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

Odpowiedź: θp = 35,78 °

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σ1​σ2​

Kąt na kole

Na kole Mohra punkt x leży pod kątem 2θp od osi σ; ćwiartka wynika ze znaków σₓ − σᵧ i τₓᵧ.

2θp = atan2(2·60, 20 − −20) = 71,57°

Kąt na elemencie

θp = 35,78°

Odpowiedź

θp = 35,78°

Kąty dodatnie mierzy się przeciwnie do ruchu wskazówek zegara.

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

Odpowiedź: σ1 = 412 MPa

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

Naprężenia na powierzchni

σ = 32M/(πd³) = 32·200000/(π·20³) = 254,6 MPa
τ = 16T/(πd³) = 16·400000/(π·20³) = 254,6 MPa
στxyCσ1​σ2​

Koło Mohra

C = σ/2 = 127,3 MPa
R = √(127,3² + 254,6²) = 284,7 MPa
σ1 = C + R = 412 MPa
σ2 = C − R = −157,4 MPa

Odpowiedź

σ1 = 412 MPa
11.Ugięcie belki

Odpowiedź: δ = 12,78 mm

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

Spójne jednostki

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

w = 8 N/mm, L = 1500 mm
EI = 3,96 × 10¹¹ N·mm²

Ugięcie

δ = wL⁴/(8EI) = 8·1500⁴/(8·3,96 × 10¹¹) = 12,78 mm

Odpowiedź

δ = 12,78 mm
12.Ugięcie belki

Odpowiedź: hreq = 290,6 mm

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

Spójne jednostki

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

w = 4 N/mm, L = 6000 mm, E = 11000 N/mm²

Warunek ugięcia

δallow = L/300 = 6000/300 = 20 mm

Wymagany moment bezwładności

Ireq = 5wL⁴/(384·E·δallow) = 5·4·6000⁴/(384·11000·20) = 3,068 × 10⁸ mm⁴

Wymagana wysokość

bh³/12 = Ireq → h = (12·3,068 × 10⁸/150)1/3 = 290,6 mm

Odpowiedź

hreq = 290,6 mm

W praktyce wybierz najbliższą większą wysokość znormalizowaną.

13.Wyboczenie słupów

Odpowiedź: Pcr = 38,59 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 = π·40⁴/64 = 125664 mm⁴

Długość wyboczeniowa

KL = 1·1500 = 1500 mm

Siła krytyczna

Pcr = π²EI/(KL)² = π²·70000·125664/1500² = 38586 N = 38,59 kN

Odpowiedź

Pcr = 38,59 kN
14.Wyboczenie słupów

Odpowiedź: d = 110,8 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·40000 = 80000 N
KL = 2·4000 = 8000 mm

Wyznacz średnicę

d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·80000·8000²/(π³·70000))^(1/4) = 110,8 mm

Odpowiedź

d = 110,8 mm

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

15.Rzut ukośny

Odpowiedź: y = 28,51 m

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.

Czas dotarcia do muru

vx = 36·cos 55° = 20,65 m/s
t = x/vx = 25/20,65 = 1,211 s

Wysokość w tej chwili

y = 29,49·1,211 − ½·9,81·1,211² = 28,51 m

Odpowiedź

y = 28,51 m

Przy murze piłka jeszcze się wznosi.

16.Rzut ukośny

Odpowiedź: t = 5,051 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 θ = 22·cos 50° = 14,14 m/s
vy = v0·sin θ = 22·sin 50° = 16,85 m/s

Czas lotu

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

0 = 40 + 16,85·t − 4,905·t²
t = (vy + √(vy² + 2gh))/g = (16,85 + √(284 + 784,8))/9,81 = 5,051 s

Odpowiedź

t = 5,051 s
17.Praca i energia

Odpowiedź: v = 11,24 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·14·(0,7071 − 0,35·0,7071)) = 11,24 m/s

Odpowiedź

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

Odpowiedź: x = 10,92 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 = 4,905·(2,5 − 0,15·4,5) = 8,952 J

Maksymalne ugięcie

½kx² = 8,952 → x = √(2·8,952/1500) = 0,1092 m = 10,92 cm

Odpowiedź

x = 10,92 cm
19.Drgania mechaniczne

Odpowiedź: fd = 22,73 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 = √(200000/5) = 200 rad/s
fn = 200/(2π) = 31,83 Hz

Częstotliwość drgań tłumionych

fd = fn·√(1 − ζ²) = 31,83·√(1 − 0,7²) = 22,73 Hz

Odpowiedź

fd = 22,73 Hz
20.Drgania mechaniczne

Odpowiedź: keq = 90 kN/m

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 = 10 + 80 = 90 kN/m

Odpowiedź

keq = 90 kN/m
21.Silniki cieplne i chłodziarki

Odpowiedź: Q̇H = 2,205 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 = 28 + 273,15 = 301,1 K
TL = 0 + 273,15 = 273,1 K

Stosunek ciepeł w obiegu odwracalnym

Q̇H = Q̇L·TH/TL = 2·301,1/273,1 = 2,205 kW

Odpowiedź

Q̇H = 2,205 kW

Moc pobierana: Ẇ = 0,205 kW.

22.Silniki cieplne i chłodziarki

Odpowiedź: Sgen = 1,377 kJ/K

Generacja entropii: w pełnym cyklu czynnik roboczy nie zmienia swojej entropii. Zbiornik gorący traci QH/TH, a zbiornik chłodny zyskuje QL/TL. Ich suma to entropia generowana, równa zeru dla obiegu odwracalnego i dodatnia dla każdego rzeczywistego.

Ciepło oddane

QL = QH − W = 2500 − 875 = 1625 kJ

Temperatury bezwzględne

TH = 300 + 273,15 = 573,1 K
TL = 10 + 273,15 = 283,1 K

Zmiany entropii zbiorników

ΔSH = −2500/573,1 = −4,362 kJ/K
ΔSL = 1625/283,1 = 5,739 kJ/K
Sgen = 5,739 − 4,362 = 1,377 kJ/K

Odpowiedź

Sgen = 1,377 kJ/K

Sgen > 0, więc silnik jest możliwy, ale nieodwracalny.

23.Obiegi silników cieplnych gazowych

Odpowiedź: η = 61,37 %

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

Człon stopnia obciążenia

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

Sprawność cieplna

η = 1 − 0,322·1,2 = 61,37%

Odpowiedź

η = 61,37%
24.Obiegi silników cieplnych gazowych

Odpowiedź: η = 57,51 %

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/200,2857 = 1 − 1/2,354 = 57,51%

Odpowiedź

η = 57,51%
25.Ciągłość i Bernoulli

Odpowiedź: p2 = 87,5 kPa

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

Ciągłość i Bernoulli

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

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

Odpowiedź

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

Odpowiedź: p2 = 89,68 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·(150/100)² = 2,25 m/s
½ρ(v1² − v2²) = ½·1000·(1² − 2,25²)/1000 = −2,031 kPa
ρgΔz = 1000·9,81·9/1000 = 88,29 kPa
p2 = p1 + ½ρ(v1² − v2²) − ρgΔz = 180 + (−2,031) − 88,29 = 89,68 kPa

Odpowiedź

p2 = 89,68 kPa
27.Tarcie w rurach i straty wysokości

Odpowiedź: hf = 7,834 m

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

Prędkość średnia

A = π·0,3²/4 = 0,07069 m²
V = Q/A = 0,12/0,07069 = 1,698 m/s

Wysokość prędkości

V²/(2g) = 1,698²/(2·9,81) = 0,1469 m

Darcy–Weisbach

hf = f·(L/D)·V²/(2g) = 0,02·(800/0,3)·1,698²/(2·9,81) = 7,834 m

Odpowiedź

hf = 7,834 m
28.Tarcie w rurach i straty wysokości

Odpowiedź: H = 2,56 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,025/0,01767 = 1,415 m/s
V²/(2g) = 0,102 m

Współczynniki strat

f·L/D = 0,018·150/0,15 = 18
ΣK = 0,5 + 2·0,3 + 5 = 6,1

Równanie energii

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

H = (1 + 18 + 6,1)·0,102 = 2,56 m

Odpowiedź

H = 2,56 m
29.Pompy i moc pomp

Odpowiedź: Pshaft = 7,988 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,015·38 = 5,592 kW

Moc na wale

Pshaft = Ph/η = 5,592/0,7 = 7,988 kW

Odpowiedź

Pshaft = 7,988 kW
30.Pompy i moc pomp

Odpowiedź: H = 37,97 m

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

Prędkość

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

Straty wysokości

hL = (f·L/D + ΣK)·V²/(2g) = (25 + 4)·0,3305 = 9,585 m

Wysokość ciśnienia

p/(ρ·g) = 200000/(1000·9,81) = 20,39 m

Wysokość podnoszenia pompy

Oba zwierciadła są praktycznie w spoczynku, więc nie ma członu wysokości prędkości.

H = Δz + p/(ρ·g) + hL = 8 + 20,39 + 9,585 = 37,97 m

Odpowiedź

H = 37,97 m
31.Przewodzenie ciepła

Odpowiedź: Q̇ = 113,3 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,17·6) = 0,1471 K/W
R2 = L2/(k2A) = (75/1000)/(0,17·6) = 0,07353 K/W

Strumień ciepła

Q̇ = (T1 − T3)/(R1 + R2) = 25/0,2206 = 113,3 W

Odpowiedź

Q̇ = 113,3 W
32.Przewodzenie ciepła

Odpowiedź: U = 1,692 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·20) = 0,005 K/W
R1 = L1/(k1A) = (50/1000)/(0,17·20) = 0,01471 K/W
R2 = L2/(k2A) = (25/1000)/(0,17·20) = 0,007353 K/W
Ro = 1/(hoA) = 1/(20·20) = 0,0025 K/W
ΣR = 0,02956 K/W

Współczynnik przenikania

U = 1/(A·ΣR) = 1/(20·0,02956) = 1,692 W/(m²·K)

Odpowiedź

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

Odpowiedź: nw = 96,67 rpm

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

Prędkość ślimacznicy

nw = n·z/Nw = 1450·3/45 = 96,67 rpm

Odpowiedź

nw = 96,67 rpm
34.Koła i przekładnie zębate

Odpowiedź: θ = 168,9 °

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.

Kąt pasa

asin(0,09688) = 5,559°

Kąt opasania

θ = 180 − 2·5,559 = 168,9°

Odpowiedź

θ = 168,9°
35.Transformaty Laplace’a

Odpowiedź: f(0.25) = 0,7852

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

Ułamki proste

F(s) = A/s + B/(s + 4) + C/(s + 6)
A = 55/(4·6) = 2,292
B = 55/((−4)·2) = −6,875
C = 55/((−6)·(−2)) = 4,583

Transformata odwrotna

f(t) = 2,292 − 6,875·e−4t + 4,583·e−6t

Oblicz

f(0,25) = 2,292 − 6,875·e−4·0,25 + 4,583·e−6·0,25 = 0,7852

Odpowiedź

f(0,25) = 0,7852
36.Transformaty Laplace’a

Odpowiedź: y(0.6) = 4,883

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

Przekształcenie równania

Y(s) = 80/(s(s + 2)(s + 4))

Ułamki proste

A = 80/(2·4) = 10
B = 80/((−2)·2) = −20
C = 80/((−4)·(−2)) = 10

Wydalone

y(t) = 10 − 20·e−2t + 10·e−4t
y(0,6) = 4,883

Odpowiedź

y(0,6) = 4,883
37.Schematy blokowe

Odpowiedź: ess = 0,1489

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

Stała uchybu położeniowego

Kₚ = G(0) = 10·8/(2·7) = 5,714

Uchyb ustalony

ess = A/(1 + Kₚ) = 1/(1 + 5,714) = 0,1489

Odpowiedź

ess = 0,1489
38.Schematy blokowe

Odpowiedź: K = 63,33

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

Odpowiedź na zakłócenie

css = D·G(0)/(1 + K·G(0)) = 10·3/(10 + 3K)

Rozwiąż względem K

10 + 3K = 10·3/0,15 = 200
K = (200 − 10)/3 = 63,33

Odpowiedź

K = 63,33
39.Odpowiedź skokowa

Odpowiedź: ts = 4 s

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

Porównaj z postacią standardową

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

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

Wynik

ζωₙ = 1 → ts ≈ 4/1 = 4 s

Odpowiedź

ts = 4 s
40.Odpowiedź skokowa

Odpowiedź: xmax = 3,056 mm

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

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

ωₙ = √(k/m) = √(5000/5) = 31,62 rad/s
ζ = c/(2√(km)) = 63/(2·√(5000·5)) = 0,1992

Wynik

xst = F/k = 10/5000 m = 2 mm
%OS = 52,8 %
xmax = 2·(1 + 0,528) = 3,056 mm

Odpowiedź

xmax = 3,056 mm
41.Uchyb ustalony i regulacja PID

Odpowiedź: ess = 1

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 przyspieszeniowa

Ka = 2·1/1 = 2

Uchyb dla paraboli

ess = 2/2 = 1

Odpowiedź

ess = 1
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) + 25·2 = 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