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 = 8 kN/mP = 70 kNABRA​ = ?2,5 m2 m0,5 m

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

L = 5 mP = 70 kN w 4,5 mw = 8 kN/m od 0 do 2,5 m
RA = kN
2.Reakcje podpór belekTrudny
w = 3 kN/mP = 55 kNABRA​RB​ = ?2 m0,5 m1,5 m2 m

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

L = 6 mP = 55 kN w 2,5 mw = 3 kN/m od 2 do 4 m
RB = kN
3.Naprężenie i odkształcenie osioweŚredni
P = ?12300 mm500 mmØ30Ø50δ = 0,5 mm

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

odcinek 1: stal, Ø30 mm, 300 mm, E = 200 GPaodcinek 2: aluminium, Ø50 mm, 500 mm, E = 70 GPaδ = 0,5 mm
P = kN
4.Naprężenie i odkształcenie osioweTrudny
Ø 35 mmStal0,5 mØ 12 mmMiedź0,5 mP = 30 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 = 30 kNodcinek 1: Ø35 mm, 0,5 m, stal (E = 200 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,5 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,5 mG = 80 GPa
φ = °
6.Skręcanie wałówTrudny
Td = ?

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

P = 75 kWN = 1450 rpmτallow = 60 MPa
d = mm
7.Siła tnąca i moment zginającyŚredni
w = 8 kN/mP = 75 kNAB3 m

Wspornik jest utwierdzony w A. Przenosi obciążenie równomiernie rozłożone w na całej długości oraz siłę skupioną P w odległości 3 m od A. Wyznacz moment zginający MA w utwierdzeniu A. Przyjmij dodatnie momenty rozciągające dolne włókna i dodatnią siłę tnącą skierowaną w górę na lewej ściance.

L = 3 mP = 75 kN w 3 mw = 8 kN/m od 0 do 3 m
MA = kN·m
8.Siła tnąca i moment zginającyTrudny
w = 5 kN/mP = 60 kNAB1,5 m2,5 m

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

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

W punkcie σₓ = 140 MPa, σᵧ = −70 MPa i τₓᵧ = 20 MPa (rozciąganie dodatnie). Element obrócono o 35° przeciwnie do ruchu wskazówek zegara. Wyznacz naprężenie normalne σy′ na ściankach obróconego elementu prostopadłych do obróconej osi y′.

σx = 140 MPaσy = −70 MPaτxy = 20 MPaθ = 35°
σy′ = MPa
10.Koło Mohra i transformacja naprężeńTrudny
35°D = 1000 mm, t = 8 mm, p = 0,8 MPap

Cienkościenny zbiornik cylindryczny o średnicy wewnętrznej 1000 mm i grubości ścianki 8 mm zawiera gaz pod nadciśnieniem 0,8 MPa. Wykonano go z blachy spawanej wzdłuż spoiny śrubowej nachylonej pod kątem 35° do osi walca. Wyznacz wartość naprężenia stycznego |τn| działającego wzdłuż spoiny.

D = 1000 mmt = 8 mmp = 0,8 MPaα = 35°
|τn| = MPa
11.Ugięcie belkiŚredni
w = 5 kN/mP = 27 kNABδ1,5 m

Wspornik przenosi obciążenie równomiernie rozłożone w na całej długości i siłę skupioną P na swobodnym końcu. Korzystając z superpozycji, wyznacz ugięcie δ swobodnego końca.

P = 27 kNw = 5 kN/mL = 1,5 mE = 200 GPa (stal)I = 160 × 10⁶ mm⁴
δ = mm
12.Ugięcie belkiTrudny
w = 20 kN/mδL = 8 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 = 20 kN/mL = 8 mE = 70 GPa (aluminium)
Ireq = × 10⁶ mm⁴
13.Wyboczenie słupówŚredni
PL = 3 md = 40 mmprzekrój poprzeczny

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

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

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

P = 135 kNn = 2,5L = 1,5 mE = 200 GPa (stal)
d = mm
15.Rzut ukośnyŚredni
v0​ = ?45°R = 120 m

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

θ = 45°R = 120 mg = 9,81 m/s²
v0 = m/s
16.Rzut ukośnyTrudny
10 mv0​ = 26 m/s40°x = ?

Piłkę rzucono z poziomu ziemi z prędkością 26 m/s pod kątem 40° do poziomu. Ląduje, już opadając, na płaskim dachu budynku o wysokości 10 m. Pomiń opór powietrza i przyjmij g = 9,81 m/s². Wyznacz odległość poziomą x od punktu wyrzutu do miejsca lądowania.

v₀ = 26 m/sθ = 40°hroof = 10 mg = 9,81 m/s²
x = m
17.Praca i energiaŚredni
80 km/h

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

v = 80 km/hμ = 0,6g = 9,81 m/s²
d = m
18.Praca i energiaTrudny
h = 5 mL = 2,5 mμ = 0,25k = 4000 N/m

Klocek o masie 4 kg puszczono bez prędkości początkowej z wysokości 5 m na gładkiej rampie. Na dole przejeżdża przez szorstki poziomy pas o długości 2,5 m (μ = 0,25), a potem trafia na sprężynę (k = 4000 N/m) na gładkiej powierzchni. Wyznacz prędkość v klocka, gdy dociera do sprężyny.

m = 4 kgh = 5 mL = 2,5 mμ = 0,25k = 4000 N/mg = 9,81 m/s²
v = m/s
19.Drgania mechaniczneŚredni
k = 100 kN/mc = 320 N·s/mm = 40 kg

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

m = 40 kgk = 100 kN/mc = 320 N·s/m
ζ =
20.Drgania mechaniczneTrudny
k1​ = 12k2​ = 40k3​ = 50m = 50 kg

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

m = 50 kgk1 = 12 kN/mk2 = 40 kN/mk3 = 50 kN/m
fn = Hz
21.Silniki cieplne i chłodziarkiŚredni
TH​ = 300 °CTL​ = 5 °CEQH​QL​W

Silnik cieplny pracuje między zbiornikiem gorącym o temperaturze 300 °C a zbiornikiem zimnym o temperaturze 5 °C. Wyznacz największą możliwą sprawność cieplną ηmax.

TH = 300 °CTL = 5 °C
ηmax = %
22.Silniki cieplne i chłodziarkiTrudny
TH​ = 700 °CT0​ = 27 °CTL​ = 0 °CERWQ̇H​ = 50 kWQ̇L​ = ?

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

TH = 700 °CT0 = 27 °CTL = 0 °CQ̇H = 50 kW
Q̇L = kW
23.Obiegi silników cieplnych gazowychŚredni
vP1234qin​qout​

Idealny obieg Diesla ma stopień sprężania r = 19 i stopień obciążenia rc = v₃/v₂ = 1,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 = 19rc = 1,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 16. 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 = 16k = 1,4
η = %
25.Ciągłość i BernoulliŚredni
12d1​ = 150 mm, p1​ = 280 kPad2​ = 100 mm, p2​ = ?

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

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

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

d1 = 150 mmd2 = 50 mmp1 = 250 kPav1 = 1 m/sΔz = 5 mρ = 1000 kg/m³
p2 = kPa
27.Tarcie w rurach i straty wysokościŚredni
V = 1,6 m/sD = 100 mmL = 50 mf = 0,03

Woda (ρ = 1000 kg/m³) płynie z prędkością 1,6 m/s przez poziomą rurę o długości 50 m i średnicy wewnętrznej 100 mm. Przyjmij współczynnik tarcia Darcy’ego f = 0,03 i g = 9,81 m/s². Wyznacz spadek ciśnienia Δp na rurze w kPa.

V = 1,6 m/sL = 50 mD = 100 mmf = 0,03ρ = 1000 kg/m³
Δp = kPa
28.Tarcie w rurach i straty wysokościTrudny
K = 10Δz = ?L = 200 m, D = 150 mm, f = 0,0183× kolano (K = 0,3)

Dwa otwarte zbiorniki są połączone rurą o długości 200 m, średnicy wewnętrznej 150 mm i współczynniku tarcia Darcy’ego f = 0,018. Straty miejscowe to: wlot o ostrej krawędzi (K = 0,5), kolana w liczbie 3 (K = 0,3 każde), zawór (K = 10) i wylot do dolnego zbiornika (K = 1,0). Przyjmij g = 9,81 m/s². Jaka różnica poziomów wody Δz między zbiornikami daje natężenie przepływu 50 L/s?

Q = 50 L/sL = 200 mD = 150 mmf = 0,018Kentrance = 0,5Kbend = 0,3 (×3)Kvalve = 10Kexit = 1,0
Δz = m
29.Pompy i moc pompŚredni
PQ = 50 L/sH = 30 m

Podczas próby pompa tłoczy 50 L/s wody (ρ = 1000 kg/m³) przy wysokości podnoszenia 30 m, a jej wał pobiera 18,4 kW. Przyjmij g = 9,81 m/s². Wyznacz sprawność pompy η w procentach.

Q = 50 L/sH = 30 mPshaft = 18,4 kWg = 9,81 m/s²
η = %
30.Pompy i moc pompTrudny
300 kPa (nadciśnienie)PQ = 80 L/sΔz = 22 mL = 150 m, D = 200 mm, f = 0,015, ΣK = 5

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

Q = 80 L/sΔz = 22 mL = 150 mD = 200 mmf = 0,015ΣK = 5ptank = 300 kPa
H = m
31.Przewodzenie ciepłaŚredni
Wełna mineralna, k = 0,0450 mmCegła, k = 0,7225 mm19 °C-13 °CTi​ = ?Q̇

Ściana składa się z 50 mm (wełna mineralna) i 25 mm (cegła), a jej pole wynosi 10 m². Powierzchnie zewnętrzne mają 19 °C i -13 °C. Wyznacz temperaturę Ti na styku dwóch warstw.

L1 = 50 mm, k1 = 0,04 W/(m·K)L2 = 25 mm, k2 = 0,72 W/(m·K)A = 10 m²T1 = 19 °CT3 = -13 °C
Ti = °C
32.Przewodzenie ciepłaTrudny
Cegła, k = 0,7250 mmSzkło, k = 0,96100 mmpomieszczenie 20 °Cna zewnątrz -14 °CQ̇

Ściana z 50 mm (cegła) i 100 mm (szkło) oddziela pomieszczenie o temperaturze 20 °C od powietrza zewnętrznego o temperaturze -14 °C. Współczynniki przejmowania ciepła wynoszą 7,7 W/(m²·K) wewnątrz i 30 W/(m²·K) na zewnątrz, a pole ściany 20 m². Wyznacz współczynnik przenikania ciepła U ściany.

L1 = 50 mm, k1 = 0,72 W/(m·K)L2 = 100 mm, k2 = 0,96 W/(m·K)A = 20 m²hi = 7,7 W/(m²·K)ho = 30 W/(m²·K)Troom = 20 °CTout = -14 °C
U = W/(m²·K)
33.Koła i przekładnie zębateŚredni
N1​N2​ = ?

Dwa równoległe wały są oddalone o 200 mm. Mają być połączone parą kół zębatych walcowych o module 5 mm i przełożeniu redukującym 4 : 1. Ile zębów N2 musi mieć koło?

a = 200 mmm = 5 mmi = 4
N2 =
34.Koła i przekładnie zębateTrudny
A = 18B = 24C = 20D = 32

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

P = 3 kWnin = 1450 rpmA = 18B = 24C = 20D = 32ηmesh = 98%
Tout = N·m
35.Transformaty Laplace’aŚredni
t (s)f(t)t = 0,3f(0,3) = ?

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

F(s) = 50/(s(s + 1)(s + 4))t = 0,3 s
f(0.3) =
36.Transformaty Laplace’aTrudny
t (s)y(t)t = 0,25y(0,25) = ?

Rozwiąż y″ + 10y′ + 24y = 0 dla y(0) = 10 i y′(0) = −5 za pomocą transformaty Laplace’a i znajdź y(0,25).

y″ + 10y′ + 24y = 0y(0) = 10y′(0) = −5t = 0,25 s
y(0.25) =
37.Schematy blokoweŚredni
R(s)+5K+7s + 6G(s)C(s)5H₁−−

Obiekt G(s) = 7/(s + 6) ma wewnętrzną (podrzędną) pętlę ujemnego sprzężenia zwrotnego ze wzmocnieniem H₁ = 5. Wzmocnienie regulatora K = 5 i zewnętrzna pętla jednostkowego ujemnego sprzężenia zwrotnego zamykają układ. Znajdź biegun układu zamkniętego p.

G(s) = 7/(s + 6)H₁ = 5K = 5
p =
38.Schematy blokoweTrudny
R(s)+15regulator+D(s)+9s + 9obiektC(s)−

Regulator proporcjonalny K = 15 steruje obiektem G(s) = 9/(s + 9) w pętli jednostkowego ujemnego sprzężenia zwrotnego. Zakłócenie D wchodzi na wejście obiektu (dodaje się do wyjścia regulatora). Przy R = 0 podano skok zakłócenia o wielkości 1. Znajdź wyjście w stanie ustalonym css.

K = 15G(s) = 9/(s + 9)D = 1R = 0
css =
39.Odpowiedź skokowaŚredni
t (s)y(t)3

Układ ma transmitancję G(s) = 72/(s² + 5s + 24) i otrzymuje jednostkowe wymuszenie skokowe. Bieguny to s = σ ± jωd. Znajdź część rzeczywistą σ.

G(s) = 72/(s² + 5s + 24)
σ =
40.Odpowiedź skokowaTrudny
+−5 VR = 134 ΩL = 100 mHC = 2 µFvC​

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

R = 134 ΩL = 100 mHC = 2 µFV = 5 V
tp = ms
41.Uchyb ustalony i regulacja PIDŚredni
R(s)+K(s + 5)(s + 10)C(s)−

Pętla ze sprzężeniem jednostkowym ma G(s) = K/((s + 5)(s + 10)). Wyznacz wzmocnienie K, które daje uchyb ustalony 25 % dla wejścia skokowego.

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

Regulator proporcjonalny K = 10 steruje obiektem G(s) = 4/(s(s + 8)) w pętli ze sprzężeniem jednostkowym. Znajdź czas szczytu tp odpowiedzi skokowej.

G(s) = 4/(s(s + 8))K = 10
tp = s

3Roztwory

1.Reakcje podpór belek

Odpowiedź: RA = 90 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 = 8·2,5 = 20 kN, x̄ = 1,25 m

Równowaga sił pionowych

RA = P + W = 70 + 20 = 90 kN

Odpowiedź

RA = 90 kN
2.Reakcje podpór belek

Odpowiedź: RB = 25,92 kN

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

w = 3 kN/mP = 55 kNABRA​RB​2 m0,5 m1,5 m2 m

Schemat sił

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

Równowaga

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

Wypadkowa obciążenia ciągłego: W = w·ℓ = 3·2 = 6 kN, przyłożona w x̄ = 3 m
ΣMA = 0: RB·6 = 55·2,5 + 3·2·3
RB = 155,5/6 = 25,92 kN
ΣFy = 0: RA = 61 − 25,92 = 35,08 kN

Odpowiedź

RB = 25,92 kN
3.Naprężenie i odkształcenie osiowe

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

Podatności

L₁/(E₁A₁) = 300/(200000·706,9) = 2,122 × 10⁻⁶ mm/N
L₂/(E₂A₂) = 500/(70000·1963) = 3,638 × 10⁻⁶ mm/N

Obciążenie

P = 0,5/5,76 × 10⁻⁶ = 86807 N = 86,81 kN

Odpowiedź

P = 86,81 kN
4.Naprężenie i odkształcenie osiowe

Odpowiedź: σmax = 265,3 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 = π·35²/4 = 962,1 mm², A2 = π·12²/4 = 113,1 mm²
σ1 = 31,18 MPa
σ2 = 265,3 MPa

Największe naprężenie

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

σmax = 265,3 MPa

Odpowiedź

σmax = 265,3 MPa
5.Skręcanie wałów

Odpowiedź: φ = 9,338 °

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·2500/(80000·38350) = 0,163 rad
φ = 0,163·180/π = 9,338°

Odpowiedź

φ = 9,338°
6.Skręcanie wałów

Odpowiedź: d = 34,74 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π·1450/60 = 151,8 rad/s
T = P/ω = 75000/151,8 = 493,9 N·m

Średnica z naprężenia dopuszczalnego

d = ∛(16T/(π·τallow)) = ∛(16·493929/(π·60)) = 34,74 mm

Odpowiedź

d = 34,74 mm

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

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

Odpowiedź: MA = −261 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.

Zostaw część ze swobodnym końcem

P: −75·3 = −225 kN·m
w: −8·3²/2 = −36 kN·m
MA = −225 − 36 = −261 kN·m

Odpowiedź

MA = −261 kN·m

Wartość ujemna oznacza moment ujemny: rozciągane są górne włókna.

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

Odpowiedź: Mmax = 65,63 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·ℓ = 5·4 = 20 kN, przyłożona w x̄ = 2 m
ΣMA = 0: RB·4 = 60·1,5 + 5·4·2
RB = 130/4 = 32,5 kN
ΣFy = 0: RA = 80 − 32,5 = 47,5 kN

Gdzie V = 0?

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

Moment maksymalny

Mmax = RA·1,5 − w·1,5²/2 = 47,5·1,5 − 5·1,5²/2 = 65,63 kN·m

Odpowiedź

Mmax = 65,63 kN·m

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

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

Odpowiedź: σy′ = −19,71 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 dla ścianki y′

σy′ = 35 − 105·0,342 − 20·0,9397 = −19,71 MPa

Sprawdzenie

Suma naprężeń normalnych nie zmienia się przy obrocie.

σx′ = 89,71 MPa, σx′ + σy′ = 70 = σx + σy

Odpowiedź

σy′ = −19,71 MPa
10.Koło Mohra i transformacja naprężeń

Odpowiedź: |τn| = 11,75 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. Cylinder cienkościenny: naprężenie obwodowe σh = pr/t i naprężenie osiowe σa = pr/(2t); są to naprężenia główne w ściance.

Naprężenia w ściance

r = 1000/2 = 500 mm
σh = pr/t = 0,8·500/8 = 50 MPa
σa = pr/(2t) = 25 MPa

Transformacja na kierunek spoiny

Przyjmij oś x wzdłuż osi walca (σa), a oś y po obwodzie (σh). Normalna spoiny tworzy z osią kąt 90° − α.

|τn| = (50 − 25)·sin 35°·cos 35° = 11,75 MPa

Odpowiedź

|τn| = 11,75 MPa
11.Ugięcie belki

Odpowiedź: δ = 1,048 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 = 27000 N, w = 5 N/mm, L = 1500 mm
EI = 3,2 × 10¹³ N·mm²

Każde obciążenie osobno

δP = PL³/(3EI) = 27000·1500³/(3·3,2 × 10¹³) = 0,9492 mm
δw = wL⁴/(8EI) = 5·1500⁴/(8·3,2 × 10¹³) = 0,09888 mm

Dodaj

δ = 0,9492 + 0,09888 = 1,048 mm

Odpowiedź

δ = 1,048 mm
12.Ugięcie belki

Odpowiedź: Ireq = 685,7 × 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 = 8000/360 = 22,22 mm

Wymagany moment bezwładności

I = 5wL⁴/(384·E·δallow) = 5·20·8000⁴/(384·70000·22,22) = 6,857 × 10⁸ mm⁴
Ireq = 685,7 × 10⁶ mm⁴

Odpowiedź

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

Odpowiedź: Pcr = 110,2 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 = 0,5·3000 = 1500 mm

Siła krytyczna

Pcr = π²EI/(KL)² = π²·200000·125664/1500² = 110245 N = 110,2 kN

Odpowiedź

Pcr = 110,2 kN
14.Wyboczenie słupów

Odpowiedź: d = 74,83 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,5·135000 = 337500 N
KL = 2·1500 = 3000 mm

Wyznacz średnicę

d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·337500·3000²/(π³·200000))^(1/4) = 74,83 mm

Odpowiedź

d = 74,83 mm

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

15.Rzut ukośny

Odpowiedź: v0 = 34,31 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·120/sin 90°) = 34,31 m/s

Odpowiedź

v0 = 34,31 m/s
16.Rzut ukośny

Odpowiedź: x = 52,44 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.

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

vx = 19,92 m/s, vy = 16,71 m/s

Czas do lądowania na dachu

Większy pierwiastek: piłka opada.

10 = 16,71·t − 4,905·t²
t = (16,71 + √(16,71² − 2·9,81·10))/9,81 = 2,633 s

Odległość pozioma

x = vx·t = 19,92·2,633 = 52,44 m

Odpowiedź

x = 52,44 m
17.Praca i energia

Odpowiedź: d = 41,95 m

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

Prędkość w jednostkach podstawowych

v = 80/3,6 = 22,22 m/s

Bilans energii

½mv² = μmg·d → d = v²/(2μg) = 22,22²/(2·0,6·9,81) = 41,95 m

Odpowiedź

d = 41,95 m

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

18.Praca i energia

Odpowiedź: v = 9,265 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).

Energia po szorstkim pasie

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

mgh − μmgL = 39,24·(5 − 0,25·2,5) = 171,7 J

Prędkość przy sprężynie

½mv² = mg(h − μL) → v = √(2g(h − μL)) = √(2·9,81·4,375) = 9,265 m/s

Odpowiedź

v = 9,265 m/s
19.Drgania mechaniczne

Odpowiedź: ζ = 0,08

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·√(100000·40) = 4000 N·s/m

Współczynnik tłumienia

ζ = c/cc = 320/4000 = 0,08

Odpowiedź

ζ = 0,08

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

20.Drgania mechaniczne

Odpowiedź: fn = 3,594 Hz

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

Sztywność zastępcza

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

k12 = 12 + 40 = 52 kN/m
keq = 52·50/(52 + 50) = 25,49 kN/m

Częstotliwość drgań własnych

ωn = √(25490/50) = 22,58 rad/s
fn = 22,58/(2π) = 3,594 Hz

Odpowiedź

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

Odpowiedź: ηmax = 51,47 %

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 = 300 + 273,15 = 573,1 K
TL = 5 + 273,15 = 278,1 K

Sprawność Carnota

ηmax = 1 − TL/TH = 1 − 278,1/573,1 = 51,47 %

Odpowiedź

ηmax = 51,47 %
22.Silniki cieplne i chłodziarki

Odpowiedź: Q̇L = 349,8 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 = 700 + 273,15 = 973,1 K
T0 = 27 + 273,15 = 300,1 K
TL = 0 + 273,15 = 273,1 K

Praca silnika

η = 1 − 300,1/973,1 = 0,6916
Ẇ = 0,6916·50 = 34,58 kW

Ciepło odbierane przez chłodziarkę

COPR = 273,1/(300,1 − 273,1) = 10,12
Q̇L = 10,12·34,58 = 349,8 kW

Odpowiedź

Q̇L = 349,8 kW
23.Obiegi silników cieplnych gazowych

Odpowiedź: T3 = 1413 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·190,4 = 941,7 K

Doprowadzenie ciepła przy stałym ciśnieniu

T3 = rc·T2 = 1,5·941,7 = 1413 K

Odpowiedź

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

Odpowiedź: η = 54,71 %

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/160,2857 = 1 − 1/2,208 = 54,71%

Odpowiedź

η = 54,71%
25.Ciągłość i Bernoulli

Odpowiedź: p2 = 267,3 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)² = 2,5·(150/100)² = 5,625 m/s
½ρ(v1² − v2²) = ½·1000·(2,5² − 5,625²)/1000 = −12,7 kPa
p2 = p1 + ½ρ(v1² − v2²) = 280 + (−12,7) = 267,3 kPa

Odpowiedź

p2 = 267,3 kPa
26.Ciągłość i Bernoulli

Odpowiedź: p2 = 160,9 kPa

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

Ciągłość i Bernoulli

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

v2 = v1·(d1/d2)² = 1·(150/50)² = 9 m/s
½ρ(v1² − v2²) = ½·1000·(1² − 9²)/1000 = −40 kPa
ρgΔz = 1000·9,81·5/1000 = 49,05 kPa
p2 = p1 + ½ρ(v1² − v2²) − ρgΔz = 250 + (−40) − 49,05 = 160,9 kPa

Odpowiedź

p2 = 160,9 kPa
27.Tarcie w rurach i straty wysokości

Odpowiedź: Δp = 19,2 kPa

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.

Strata wysokości

hf = f·(L/D)·V²/(2g) = 0,03·(50/0,1)·1,6²/(2·9,81) = 1,957 m

Spadek ciśnienia

Rura jest pozioma, a prędkość się nie zmienia, więc równanie energii upraszcza się do Δp = ρ·g·hf.

Δp = ρ·g·hf = 1000·9,81·1,957 = 19200 Pa = 19,2 kPa

Odpowiedź

Δp = 19,2 kPa
28.Tarcie w rurach i straty wysokości

Odpowiedź: Δz = 14,85 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,05/0,01767 = 2,829 m/s
V²/(2g) = 0,408 m

Współczynniki strat

f·L/D = 0,018·200/0,15 = 24
ΣK = 0,5 + 3·0,3 + 10 + 1,0 = 12,4

Całkowita strata wysokości

hf = 24·0,408 = 9,793 m
hm = 12,4·0,408 = 5,06 m
Δz = hf + hm = 14,85 m

Odpowiedź

Δz = 14,85 m
29.Pompy i moc pomp

Odpowiedź: η = 79,97 %

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,05·30 = 14,71 kW

Sprawność

η = Ph/Pshaft = 14,71/18,4 = 79,97%

Odpowiedź

η = 79,97%
30.Pompy i moc pomp

Odpowiedź: H = 57,95 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,08/0,03142 = 2,546 m/s
V²/(2g) = 0,3305 m

Straty wysokości

hL = (f·L/D + ΣK)·V²/(2g) = (11,25 + 5)·0,3305 = 5,371 m

Wysokość ciśnienia

p/(ρ·g) = 300000/(1000·9,81) = 30,58 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 = 22 + 30,58 + 5,371 = 57,95 m

Odpowiedź

H = 57,95 m
31.Przewodzenie ciepła

Odpowiedź: Ti = −12,14 °C

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

Opory warstw

R1 = L1/(k1A) = (50/1000)/(0,04·10) = 0,125 K/W
R2 = L2/(k2A) = (25/1000)/(0,72·10) = 0,003472 K/W

Strumień ciepła

Q̇ = (T1 − T3)/(R1 + R2) = 32/0,1285 = 249,1 W

Temperatura na styku

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

Ti = T1 − Q̇·R1 = 19 − 249,1·0,125 = −12,14 °C

Odpowiedź

Ti = −12,14 °C
32.Przewodzenie ciepła

Odpowiedź: U = 2,969 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/(7,7·20) = 0,006494 K/W
R1 = L1/(k1A) = (50/1000)/(0,72·20) = 0,003472 K/W
R2 = L2/(k2A) = (100/1000)/(0,96·20) = 0,005208 K/W
Ro = 1/(hoA) = 1/(30·20) = 0,001667 K/W
ΣR = 0,01684 K/W

Współczynnik przenikania

U = 1/(A·ΣR) = 1/(20·0,01684) = 2,969 W/(m²·K)

Odpowiedź

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

Odpowiedź: N2 = 64

Gear geometry: the module m sets the tooth size. A gear with N teeth has a pitch diameter d = m·N, two meshing gears have the same module, and their centre distance is a = (d₁ + d₂)/2 = m·(N₁ + N₂)/2. The pitch-line speed is v = π·d·n/60, and the tangential tooth force follows from the torque: Ft = 2T/d.

Łączna liczba zębów

N1 + N2 = 2a/m = 2·200/5 = 80

Podział według przełożenia

N1 = 80/(1 + 4) = 16
N2 = 4·16 = 64

Odpowiedź

N2 = 64
34.Koła i przekładnie zębate

Odpowiedź: Tout = 40,48 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 = 1450·(18/24)·(20/32) = 679,7 rpm
ωout = 2π·679,7/60 = 71,18 rad/s

Moc i moment wyjściowy

Pout = 3000·0,98² = 2881 W
Tout = Pout/ωout = 2881/71,18 = 40,48 N·m

Odpowiedź

Tout = 40,48 N·m
35.Transformaty Laplace’a

Odpowiedź: f(0.3) = 1,408

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 + 1) + C/(s + 4)
A = 50/(1·4) = 12,5
B = 50/((−1)·3) = −16,67
C = 50/((−4)·(−3)) = 4,167

Transformata odwrotna

f(t) = 12,5 − 16,67·e−t + 4,167·e−4t

Oblicz

f(0,3) = 12,5 − 16,67·e−0,3 + 4,167·e−4·0,3 = 1,408

Odpowiedź

f(0,3) = 1,408
36.Transformaty Laplace’a

Odpowiedź: y(0.25) = 6,212

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 − 10s − (−5) + 10(sY − 10) + 24Y = 0
Y(s) = (10s + 95)/((s + 4)(s + 6))

Ułamki proste

A = 55/2 = 27,5
B = 35/(−2) = −17,5

Wynik

y(t) = 27,5·e−4t − 17,5·e−6t
y(0,25) = 6,212

Odpowiedź

y(0,25) = 6,212

Sprawdzenie: y(0) = A + B = 10.

37.Schematy blokowe

Odpowiedź: p = −76

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

Pętla wewnętrzna

Gᵢ = G/(1 + G·H₁) = 7/(s + 6 + 7·5) = 7/(s + 41)

Pętla zewnętrzna

T(s) = K·Gᵢ/(1 + K·Gᵢ) = 35/(s + 41 + 35) = 35/(s + 76)

Wynik

p = −76

Odpowiedź

p = −76
38.Schematy blokowe

Odpowiedź: css = 0,0625

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

Transmitancje

C/D = G/(1 + K·G) = 9/(s + 9 + 135)

Stan ustalony

css = 1·9/144 = 0,0625

Odpowiedź

css = 0,0625

Większe K zmniejsza wpływ zakłócenia, ale przy regulatorze proporcjonalnym nigdy do zera.

39.Odpowiedź skokowa

Odpowiedź: σ = −2,5

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

ωₙ² = 24 → ωₙ = 4,899 rad/s
2ζωₙ = 5 → ζ = 5/(2·4,899) = 0,5103

Wynik

σ = −ζωₙ = −5/2 = −2,5

Odpowiedź

σ = −2,5
40.Odpowiedź skokowa

Odpowiedź: tp = 1,473 ms

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

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

ωₙ = 1/√(LC) = 1/√(0,1·2 × 10⁻⁶) = 2236 rad/s
ζ = (R/2)·√(C/L) = (134/2)·√(2 × 10⁻⁶/0,1) = 0,2996

Wynik

ωd = 2133 rad/s
tp = π/ωd = 1,473 ms

Odpowiedź

tp = 1,473 ms
41.Uchyb ustalony i regulacja PID

Odpowiedź: K = 150

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.

Wymagana stała uchybowa

1/(1 + Kp) = 0,25 → Kp = 3

Wzmocnienie

K/50 = 3 → K = 150

Odpowiedź

K = 150

Sprawdź też stabilność układu zamkniętego: dla tej pętli drugiego rzędu każde K > 0 jest stabilne.

42.Uchyb ustalony i regulacja PID

Odpowiedź: tp = 0,6413 s

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 + 8) + 10·4 = s² + 8s + 40

ωₙ i ζ

ωₙ = √40 = 6,325 rad/s
ζ = 8/(2·6,325) = 0,6325

Czas szczytu

ωd = 6,325·√(1 − 0,6325²) = 4,899 rad/s
tp = π/4,899 = 0,6413 s

Odpowiedź

tp = 0,6413 s

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