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
Narysuj schemat sił i zastąp każdą podporę jej reakcjami.
Zastąp każde obciążenie ciągłe jedną siłą przyłożoną w jego środku ciężkości.
Napisz równanie momentów względem jednej podpory, aby wyznaczyć reakcję w drugiej.
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
Wyznacz biegunowy moment bezwładności: J = πd⁴/32 dla wału pełnego.
Naprężenie styczne na powierzchni: τ = T·r/J, gdzie r to promień zewnętrzny.
Kąt skręcenia: φ = TL/(GJ) w radianach.
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
Wyznacz reakcje podpór z równowagi całej belki.
Przetnij belkę w potrzebnym przekroju i zostaw jedną stronę.
Siła tnąca: suma sił pionowych po tej stronie. Moment zginający: suma ich momentów względem przekroju.
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
Oblicz środek koła Mohra: C = (σx + σy)/2.
Oblicz promień: R = √(((σx − σy)/2)² + τxy²).
Naprężenia główne to C ± R, a największe naprężenie styczne w płaszczyźnie to R.
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
Rozłóż prędkość początkową na składowe: vx = v₀·cos θ i vy = v₀·sin θ.
Ruch poziomy ma stałą prędkość; ruch pionowy ma przyspieszenie −g.
Wyznacz czas lotu z ruchu pionowego.
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
Policz bieguny układu otwartego C(s)·G(s) w s = 0: to typ układu.
Wyznacz stałą uchybową pasującą do wejścia: Kp dla skoku, Kv dla narastania liniowego, Ka dla paraboli.
Zastosuj ess = A/(1 + Kp), A/Kv lub A/Ka albo rozwiąż to względem potrzebnego wzmocnienia.
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
Belka swobodnie podparta (podpora nieprzesuwna w A, przesuwna w B) przenosi obciążenie liniowo zmienne, rosnące od zera w A do 15 kN/m w B. W odległości 4,5 m od A działa też siła skupiona P. Wyznacz pionową reakcję podpory RB.
L = 6 mw0 = 15 kN/mP = 55 kN w 4,5 m
RB = kN
2.Reakcje podpór belekTrudny
Belka swobodnie podparta ma podporę przegubową nieprzesuwną w A i przesuwną w B. Wyznacz pionową reakcję podpory RB.
L = 4 mP = 55 kN w 1 mw = 2 kN/m od 0 do 4 m
RB = kN
3.Naprężenie i odkształcenie osioweŚredni
Pręt (materiał: stal, średnica 40 mm, E = 200 GPa) znajduje się wewnątrz rury (materiał: mosiądz, średnica zewnętrzna 100 mm, średnica wewnętrzna 70 mm, E = 100 GPa). Oba mają długość 600 mm i są razem ściskane między sztywnymi płytami siłą 250 kN. Wyznacz skrócenie δ.
P = 250 kNpręt: Ø40 mm, E = 200 GParura: Ø100/Ø70 mm, E = 100 GPaL = 600 mm
δ = mm
4.Naprężenie i odkształcenie osioweTrudny
Pręt (materiał: stal, średnica 25 mm, E = 200 GPa, naprężenie dopuszczalne 250 MPa) znajduje się wewnątrz rury (materiał: mosiądz, średnica zewnętrzna 100 mm, średnica wewnętrzna 60 mm, E = 100 GPa, naprężenie dopuszczalne 80 MPa). Są razem ściskane między sztywnymi płytami. Wyznacz największe obciążenie Pmax.
pręt: Ø25 mm, E = 200 GPa, σallow = 250 MParura: Ø100/Ø60 mm, E = 100 GPa, σallow = 80 MPa
Pmax = kN
5.Skręcanie wałówŚredni
W próbie skręcania pełna próbka o średnicy 20 mm jest skręcana momentem 20 N·m. Końce odcinka pomiarowego o długości 300 mm obracają się względem siebie o 0,56° (materiał jest jeszcze w zakresie sprężystym). Wyznacz moduł sprężystości postaciowej G materiału.
T = 20 N·md = 20 mmL = 300 mmφ = 0,56°
G = GPa
6.Skręcanie wałówTrudny
Wał drążony przenosi 50 kW przy 1450 obr/min. Wyznacz maksymalne naprężenie styczne τmax w wale.
P = 50 kWN = 1450 rpmdo = 40 mmdi = 30 mm
τmax = MPa
7.Siła tnąca i moment zginającyŚredni
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.
30 kN @ 4 m12 kN/m, 2–6 m
1A (x = 0 m)
2B (x = 1 m)
3C (x = 2 m)
4D (x = 3 m)
5E (x = 4 m)
6F (x = 5 m)
7G (x = 6 m)
8.Siła tnąca i moment zginającyTrudny
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 = 20 kN w 2 mw = 3 kN/m od 0 do 4 m
Mmax = kN·m
9.Koło Mohra i transformacja naprężeńŚredni
W pewnym punkcie obciążonej części występują σₓ = −70 MPa, σᵧ = −30 MPa i τₓᵧ = −40 MPa (rozciąganie dodatnie). Wyznacz maksymalne naprężenie styczne w płaszczyźnie τmax.
σx = −70 MPaσy = −30 MPaτxy = −40 MPa
τmax = MPa
10.Koło Mohra i transformacja naprężeńTrudny
Punkt cienkiej płyty jest w płaskim stanie naprężenia: σₓ = −80 MPa, σᵧ = −30 MPa i τₓᵧ = 10 MPa (σz = 0). Granica plastyczności materiału wynosi 300 MPa. Wyznacz współczynnik bezpieczeństwa n względem uplastycznienia według hipotezy Hubera-Misesa (energii odkształcenia postaciowego).
σx = −80 MPaσy = −30 MPaτxy = 10 MPaσY = 300 MPa
n =
11.Ugięcie belkiŚredni
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 = 14 kNw = 12 kN/mL = 3 mE = 200 GPa (stal)I = 36 × 10⁶ mm⁴
δ = mm
12.Ugięcie belkiTrudny
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 = 10 kN/mL = 1,5 mE = 70 GPa (aluminium)
Ireq = × 10⁶ mm⁴
13.Wyboczenie słupówŚredni
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 = 5 mE = 200 GPa (stal)d = 40 mm
Pcr = kN
14.Wyboczenie słupówTrudny
Słup o długości L z materiału (stal) jest utwierdzony u podstawy i wolny u góry. Wyznacz jego smukłość i naprężenie krytyczne σcr ze wzoru Eulera.
L = 5 mE = 200 GPa (stal)d = 50 mm
σcr = MPa
15.Rzut ukośnyŚredni
Piłka zostaje wyrzucona pod kątem 60° do poziomu z płaskiego terenu. Pomiń opór powietrza. Wyznacz największą wysokość H, jaką osiąga piłka.
v0 = 30 m/sθ = 60°g = 9,81 m/s²
H = m
16.Rzut ukośnyTrudny
Piłka zostaje wyrzucona pod kątem 50° do poziomu ze szczytu klifu o wysokości 10 m. Pomiń opór powietrza. Wyznacz odległość poziomą x pokonaną, zanim piłka uderzy w ziemię.
v0 = 19 m/sθ = 50°h = 10 mg = 9,81 m/s²
x = m
17.Praca i energiaŚredni
Samochód jadący z prędkością 100 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,8. Wyznacz drogę hamowania d.
v = 100 km/hμ = 0,8g = 9,81 m/s²
d = m
18.Praca i energiaTrudny
Klocek o masie 10 kg puszczono bez prędkości początkowej z wysokości 3,5 m na gładkiej rampie. Na dole przejeżdża przez szorstki poziomy pas o długości 3 m (μ = 0,3), a potem trafia na sprężynę (k = 2500 N/m) na gładkiej powierzchni. Wyznacz maksymalne ugięcie x sprężyny.
m = 10 kgh = 3,5 mL = 3 mμ = 0,3k = 2500 N/mg = 9,81 m/s²
x = cm
19.Drgania mechaniczneŚredni
W próbie drgań swobodnych tłumionego układu masa–sprężyna po maksimum przemieszczenia 50 mm po 3 cyklach następuje maksimum 16,1 mm. Wyznacz współczynnik tłumienia ζ.
x0 = 50 mmx3 = 16,1 mmn = 3
ζ =
20.Drgania mechaniczneTrudny
Masa 25 kg wisi na sprężynie k₃ = 12 kN/m. Górny koniec k₃ jest przymocowany do sztywnej belki, która wisi na dwóch sprężynach obok siebie, k₁ = 30 kN/m i k₂ = 20 kN/m; belka pozostaje pozioma. Pomiń tłumienie oraz masy sprężyn i belki. Wyznacz częstotliwość drgań własnych fn w kierunku pionowym w Hz.
m = 25 kgk1 = 30 kN/mk2 = 20 kN/mk3 = 12 kN/m
fn = Hz
21.Silniki cieplne i chłodziarkiŚredni
Silnik cieplny pracuje między zbiornikiem gorącym o temperaturze 800 °C a zbiornikiem zimnym o temperaturze 15 °C. Wyznacz największą możliwą sprawność cieplną ηmax.
TH = 800 °CTL = 15 °C
ηmax = %
22.Silniki cieplne i chłodziarkiTrudny
Silnik cieplny Carnota pobiera 80 kW ciepła ze źródła o temperaturze 650 °C i oddaje ciepło do otoczenia o temperaturze 20 °C. Cała jego praca napędza chłodziarkę Carnota, która odbiera ciepło z zamrażarki o temperaturze 0 °C i również oddaje ciepło do otoczenia o temperaturze 20 °C. Wyznacz strumień ciepła Q̇L odbieranego z zamrażarki.
TH = 650 °CT0 = 20 °CTL = 0 °CQ̇H = 80 kW
Q̇L = kW
23.Obiegi silników cieplnych gazowychŚredni
Idealny obieg Diesla ma stopień sprężania r = 14 i stopień obciążenia rc = v₃/v₂ = 2,5. Na początku sprężania powietrze ma temperaturę 310 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 doprowadzone ciepło qin na kilogram powietrza, w kJ/kg.
r = 14rc = 2,5T1 = 310 Kk = 1,4cp = 1,005 kJ/(kg·K)
qin = kJ/kg
24.Obiegi silników cieplnych gazowychTrudny
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
Rurka Prandtla w tunelu aerodynamicznym wskazuje różnicę 50 Pa między ciśnieniem spiętrzenia a ciśnieniem statycznym. Gęstość powietrza wynosi 1,2 kg/m³. Wyznacz prędkość powietrza v.
p0 − p = 50 Paρ = 1,2 kg/m³
v = m/s
26.Ciągłość i BernoulliTrudny
Wąż o średnicy wewnętrznej 32 mm działa jak lewar i odprowadza wodę z otwartego zbiornika. Najwyższy punkt węża C leży 1 m nad powierzchnią wody, a wylot (otwarty do atmosfery) 2 m pod powierzchnią wody. Pomiń straty i przyjmij g = 9,81 m/s². Wyznacz natężenie przepływu Q.
d = 32 mmzC = +1 mz2 = −2 mg = 9,81 m/s²
Q = L/s
27.Tarcie w rurach i straty wysokościŚredni
Woda (ρ = 1000 kg/m³) płynie z prędkością 2,2 m/s przez poziomą rurę o długości 800 m i średnicy wewnętrznej 300 mm. Przyjmij współczynnik tarcia Darcy’ego f = 0,018 i g = 9,81 m/s². Wyznacz spadek ciśnienia Δp na rurze w kPa.
Rura o długości 100 m i średnicy wewnętrznej 200 mm (współczynnik tarcia Darcy’ego f = 0,02) 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 = 10). 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 60 L/s?
Pompa tłoczy 100 L/s wody (ρ = 1000 kg/m³) przy wysokości podnoszenia 80 m. Sprawność pompy wynosi 78%, a napędzającego ją silnika elektrycznego 88%. Przyjmij g = 9,81 m/s². Wyznacz moc elektryczną Pel pobieraną przez silnik, w kW.
Pompa przetłacza 40 L/s wody (ρ = 1000 kg/m³) z otwartego dolnego zbiornika do górnego zbiornika, którego zwierciadło wody leży 20 m wyżej. Rura ma długość 300 m, średnicę wewnętrzną 150 mm i współczynnik tarcia Darcy’ego f = 0,022; współczynniki strat miejscowych sumują się do ΣK = 8 (łącznie z wlotem i wylotem). Przyjmij g = 9,81 m/s². Wyznacz wysokość podnoszenia H, którą musi zapewnić pompa.
Ściana składa się z 50 mm (szkło) i 25 mm (płyta gipsowo-kartonowa), a jej pole wynosi 4 m². Powierzchnie zewnętrzne mają 20 °C i -1 °C. Wyznacz temperaturę Ti na styku dwóch warstw.
L1 = 50 mm, k1 = 0,96 W/(m·K)L2 = 25 mm, k2 = 0,17 W/(m·K)A = 4 m²T1 = 20 °CT3 = -1 °C
Ti = °C
32.Przewodzenie ciepłaTrudny
Ściana z 250 mm (styropian) i 75 mm (drewno) oddziela pomieszczenie o temperaturze 21 °C od powietrza zewnętrznego o temperaturze 3 °C. Współczynniki przejmowania ciepła wynoszą 8 W/(m²·K) wewnątrz i 30 W/(m²·K) na zewnątrz, a pole ściany 4 m². Wyznacz straty ciepła Q̇ przez ścianę.
L1 = 250 mm, k1 = 0,035 W/(m·K)L2 = 75 mm, k2 = 0,17 W/(m·K)A = 4 m²hi = 8 W/(m²·K)ho = 30 W/(m²·K)Troom = 21 °CTout = 3 °C
Q̇ = W
33.Koła i przekładnie zębateŚredni
Silnik o mocy 1,5 kW pracujący z prędkością 1450 rpm napędza zębnik o 25 zębach, który zazębia się z kołem o 80 zębach. Sprawność pary kół wynosi 98%. Wyznacz moment obrotowy T1 na wale silnika.
P = 1,5 kWn1 = 1450 rpmN1 = 25N2 = 80η = 98%
T1 = N·m
34.Koła i przekładnie zębateTrudny
Pas otwarty łączy koło 100 mm i koło 250 mm, których środki są oddalone o 500 mm. Wyznacz kąt opasania θ na mniejszym kole, w stopniach.
d1 = 100 mmd2 = 250 mmC = 500 mm
θ = °
35.Transformaty Laplace’aŚredni
F(s) = (5s + 35)/(s² + 11s + 30). Jej transformata odwrotna ma postać f(t) = A·e−5t + B·e−6t. Znajdź współczynnik A metodą rozkładu na ułamki proste.
F(s) = (5s + 35)/(s² + 11s + 30)
A =
36.Transformaty Laplace’aTrudny
Układ spełnia równanie y″ + 3y′ + 2y = 8·u(t), przy y(0) = 5 i y′(0) = 0. W chwili t = 0 podano na wejście skok u(t) = 2. Za pomocą transformaty Laplace’a znajdź y(0,25).
Obiekt G(s) = 19/(s + 1) jest regulowany w pętli ujemnego sprzężenia zwrotnego, w której czujnik ma wzmocnienie H = 0,25. Znajdź stałą czasową τ układu zamkniętego.
G(s) = 19/(s + 1)H = 0,25
τ = s
38.Schematy blokoweTrudny
Regulator proporcjonalny K steruje obiektem G(s) = 2/(s + 6) w pętli jednostkowego ujemnego sprzężenia zwrotnego. Skok zakłócenia o wielkości 1 wchodzi na wejście obiektu, gdy R = 0. Znajdź najmniejsze wzmocnienie K, przy którym wyjście w stanie ustalonym nie przekracza 0,06667.
G(s) = 2/(s + 6)D = 1css ≤ 0,06667
K =
39.Odpowiedź skokowaŚredni
Układ ma transmitancję G(s) = 6/(s² + 3s + 12) i otrzymuje jednostkowe wymuszenie skokowe. Znajdź pulsację własną ωn w rad/s.
G(s) = 6/(s² + 3s + 12)
ωn = rad/s
40.Odpowiedź skokowaTrudny
Masa m = 5 kg wisi na sprężynie k = 5000 N/m z tłumikiem lepkim c = 32 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 = 25 N. Znajdź największe wychylenie xmax w mm.
m = 5 kgk = 5000 N/mc = 32 N·s/mF = 25 N
xmax = mm
41.Uchyb ustalony i regulacja PIDŚredni
Na pętlę ze sprzężeniem jednostkowym o G(s) = 4(s + 3)/(s²(s + 5)) podano wejście paraboliczne r(t) = 2·t²/2. Wyznacz uchyb ustalony ess.
G(s) = 4(s + 3)/(s²(s + 5))r(t) = 2·t²/2
ess =
42.Uchyb ustalony i regulacja PIDTrudny
Regulator PD C(s) = 20 + Kd·s steruje obiektem G(s) = 4/(s(s + 4)) w pętli ze sprzężeniem jednostkowym. Wyznacz wzmocnienie członu różniczkującego Kd, przy którym układ zamknięty ma współczynnik tłumienia ζ = 0,9.
G(s) = 4/(s(s + 4))Kp = 20ζ = 0,9
Kd = s
3Roztwory
1.Reakcje podpór belek
Odpowiedź: RB = 71,25 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 trójkątnego
W = ½·15·6 = 45 kN
x̄ = 2L/3 = 4 m
Równowaga
Momenty względem A, dodatnie przeciwnie do ruchu wskazówek zegara. Potem równowaga sił pionowych.
ΣMA = 0: RB·6 = 45·4 + 55·4,5
RB = 71,25 kN
RA = 100 − 71,25 = 28,75 kN
Odpowiedź
RB = 71,25 kN
2.Reakcje podpór belek
Odpowiedź: RB = 17,75 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.
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·ℓ = 2·4 = 8 kN, przyłożona w x̄ = 2 m
ΣMA = 0: RB·4 = 55·1 + 2·4·2
RB = 71/4 = 17,75 kN
ΣFy = 0: RA = 63 − 17,75 = 45,25 kN
Odpowiedź
RB = 17,75 kN
3.Naprężenie i odkształcenie osiowe
Odpowiedź: δ = 0,2301 mm
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
Ar = π·40²/4 = 1257 mm²
At = π(100² − 70²)/4 = 4006 mm²
Łączna sztywność osiowa
ΣEA = 200·10³·1257 + 100·10³·4006 = 6,519 × 10⁸ N
Skrócenie
δ = 250000·600/6,519 × 10⁸ = 0,2301 mm
Odpowiedź
δ = 0,2301 mm
4.Naprężenie i odkształcenie osiowe
Odpowiedź: Pmax = 480,7 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
Ar = π·25²/4 = 490,9 mm²
At = π(100² − 60²)/4 = 5027 mm²
Łączna sztywność osiowa
ΣEA = 200·10³·490,9 + 100·10³·5027 = 6,008 × 10⁸ N
Decydujące odkształcenie
ε = min(0,00125, 0,0008) = 0,0008 (mosiądz)
Największe obciążenie
Pmax = 0,0008·6,008 × 10⁸ = 480664 N = 480,7 kN
Odpowiedź
Pmax = 480,7 kN
5.Skręcanie wałów
Odpowiedź: G = 39,08 GPa
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 i kąt
J = π·20⁴/32 = 15708 mm⁴
φ = 0,56·π/180 = 0,009774 rad
Moduł sprężystości poprzecznej
T = 20 N·m = 20000 N·mm
G = TL/(Jφ) = 20000·300/(15708·0,009774) = 39081 MPa
Skręcanie wałów okrągłych: naprężenie styczne rośnie liniowo od zera w środku do τmax = T·r/J na powierzchni, gdzie J = πd⁴/32 dla wału pełnego i J = π(do⁴ − di⁴)/32 dla drążonego. Kąt skręcenia to φ = TL/(GJ) w radianach. Moc i moment skręcający wiąże P = T·ω.
Moment skręcający
ω = 2πN/60 = 2π·1450/60 = 151,8 rad/s
T = P/ω = 50000/151,8 = 329,3 N·m
Biegunowy moment bezwładności
J = π(40⁴ − 30⁴)/32 = 171806 mm⁴
Maksymalne naprężenie styczne
τmax = T·(do/2)/J = 329286·20/171806 = 38,33 MPa
Odpowiedź
τmax = 38,33 MPa
7.Siła tnąca i moment zginający
Odpowiedź: E (x = 4 m)
Największy moment zginający: nachylenie wykresu M zmienia się tam, gdzie zmienia się siła tnąca, więc jego maksima leżą tam, gdzie siła tnąca przechodzi przez zero: pod siłą skupioną, wewnątrz obciążenia ciągłego albo nad podporą przy wsporniku. W utwierdzeniu belki wspornikowej moment jest największy przy ścianie. Porównaj wartości M w tych miejscach, zarówno dodatnie (zginanie „w dół”), jak i ujemne.
Reakcje podpór
Zapisz momenty względem jednej podpory, a potem użyj bilansu sił pionowych.
R1 = 26 kN, R2 = 52 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 = 26, MC = 52, MD = 72, ME = 80 kN·m, MF = 46, MG = 0
Odpowiedź
Punkt E: ME = 80 kN·m
8.Siła tnąca i moment zginający
Odpowiedź: Mmax = 26 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·4 = 12 kN, przyłożona w x̄ = 2 m
ΣMA = 0: RB·4 = 20·2 + 3·4·2
RB = 64/4 = 16 kN
ΣFy = 0: RA = 32 − 16 = 16 kN
Gdzie V = 0?
V(2⁻) = RA − w·2 = 10 kN
V(2⁺) = 10 − 20 = −10 kN
Siła tnąca zmienia znak w obrębie obciążenia, więc maksimum leży pod obciążeniem.
Moment maksymalny
Mmax = RA·2 − w·2²/2 = 16·2 − 3·2²/2 = 26 kN·m
Odpowiedź
Mmax = 26 kN·m
Występuje w x = 2 m od A.
9.Koło Mohra i transformacja naprężeń
Odpowiedź: τmax = 44,72 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.
Promień koła Mohra
τmax = R = √(((−70 − −30)/2)² + −40²) = 44,72 MPa
Odpowiedź
τmax = 44,72 MPa
10.Koło Mohra i transformacja naprężeń
Odpowiedź: n = 4,16
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. Hipotezy wytężeniowe: według hipotezy Hubera-Misesa uplastycznienie następuje, gdy σvM = √(σ₁² − σ₁σ₂ + σ₂²) osiąga σY; według hipotezy Tresci, gdy bezwzględne maksymalne naprężenie styczne osiąga σY/2.
n > 1 oznacza, że punkt się nie uplastycznia; hipoteza Tresci jest nieco bardziej zachowawcza niż hipoteza Hubera-Misesa.
11.Ugięcie belki
Odpowiedź: δ = 34,38 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 = 14000 N, w = 12 N/mm, L = 3000 mm
EI = 7,2 × 10¹² N·mm²
Każde obciążenie osobno
δP = PL³/(3EI) = 14000·3000³/(3·7,2 × 10¹²) = 17,5 mm
δw = wL⁴/(8EI) = 12·3000⁴/(8·7,2 × 10¹²) = 16,88 mm
Dodaj
δ = 17,5 + 16,88 = 34,38 mm
Odpowiedź
δ = 34,38 mm
12.Ugięcie belki
Odpowiedź: Ireq = 2,26 × 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 = 1500/360 = 4,167 mm
Wymagany moment bezwładności
I = 5wL⁴/(384·E·δallow) = 5·10·1500⁴/(384·70000·4,167) = 2,26 × 10⁶ mm⁴
Ireq = 2,26 × 10⁶ mm⁴
Odpowiedź
Ireq = 2,26 × 10⁶ mm⁴
13.Wyboczenie słupów
Odpowiedź: Pcr = 39,69 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.
Wyboczenie Eulera: smukły słup traci stateczność przy Pcr = π²EI/(KL)², gdzie KL to długość wyboczeniowa: K = 1 dla przegub–przegub, 2 dla utwierdzenie–wolny koniec, 0,7 dla utwierdzenie–przegub i 0,5 dla utwierdzenie–utwierdzenie. Wybocza się względem osi o najmniejszym I. Wzór dotyczy słupów smukłych, których naprężenie krytyczne jest mniejsze od granicy plastyczności.
Najmniejszy moment bezwładności
Przekrój kołowy ma to samo I względem każdej osi.
I = πd⁴/64 = π·50⁴/64 = 306796 mm⁴
Promień bezwładności i smukłość
A = 1963 mm²
r = √(I/A) = √(306796/1963) = 12,5 mm
KL = 2·5000 = 10000 mm
KL/r = 10000/12,5 = 800
Naprężenie krytyczne
σcr = π²E/(KL/r)² = π²·200000/800² = 3,084 MPa
Odpowiedź
σcr = 3,084 MPa
15.Rzut ukośny
Odpowiedź: H = 34,4 m
Rzut ukośny: bez oporu powietrza x = vxt i y = h + vyt − ½gt². Oba kierunki są niezależne, łączy je tylko ten sam czas.
Rozłóż prędkość początkową
Prędkość pozioma jest stała; ruch pionowy ma stałe przyspieszenie g = 9,81 m/s² w dół.
vx = v0·cos θ = 30·cos 60° = 15 m/s
vy = v0·sin θ = 30·sin 60° = 25,98 m/s
Najwyższy punkt
Na szczycie prędkość pionowa jest zerowa: vy² = 2gH.
H = vy²/(2g) = 25,98²/(2·9,81) = 34,4 m
Odpowiedź
H = 34,4 m
16.Rzut ukośny
Odpowiedź: x = 43,27 m
Rzut ukośny: bez oporu powietrza x = vxt i y = h + vyt − ½gt². Oba kierunki są niezależne, łączy je tylko ten sam czas.
Rozłóż prędkość początkową
Prędkość pozioma jest stała; ruch pionowy ma stałe przyspieszenie g = 9,81 m/s² w dół.
t = (vy + √(vy² + 2gh))/g = (14,55 + √(211,8 + 196,2))/9,81 = 3,543 s
Odległość pozioma
x = vx·t = 12,21·3,543 = 43,27 m
Odpowiedź
x = 43,27 m
17.Praca i energia
Odpowiedź: d = 49,16 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 = 100/3,6 = 27,78 m/s
Bilans energii
½mv² = μmg·d → d = v²/(2μg) = 27,78²/(2·0,8·9,81) = 49,16 m
Odpowiedź
d = 49,16 m
Podwojenie prędkości wydłuża drogę hamowania czterokrotnie.
18.Praca i energia
Odpowiedź: x = 45,17 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 = 98,1·(3,5 − 0,3·3) = 255,1 J
Maksymalne ugięcie
½kx² = 255,1 → x = √(2·255,1/2500) = 0,4517 m = 45,17 cm
Odpowiedź
x = 45,17 cm
19.Drgania mechaniczne
Odpowiedź: ζ = 0,06001
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₂.
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 = 30 + 20 = 50 kN/m
keq = 50·12/(50 + 12) = 9,677 kN/m
Częstotliwość drgań własnych
ωn = √(9677/25) = 19,67 rad/s
fn = 19,67/(2π) = 3,131 Hz
Odpowiedź
fn = 3,131 Hz
21.Silniki cieplne i chłodziarki
Odpowiedź: ηmax = 73,15 %
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 = 800 + 273,15 = 1073 K
TL = 15 + 273,15 = 288,1 K
Sprawność Carnota
ηmax = 1 − TL/TH = 1 − 288,1/1073 = 73,15 %
Odpowiedź
ηmax = 73,15 %
22.Silniki cieplne i chłodziarki
Odpowiedź: Q̇L = 745,6 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 = 650 + 273,15 = 923,1 K
T0 = 20 + 273,15 = 293,1 K
TL = 0 + 273,15 = 273,1 K
Praca silnika
η = 1 − 293,1/923,1 = 0,6824
Ẇ = 0,6824·80 = 54,6 kW
Ciepło odbierane przez chłodziarkę
COPR = 273,1/(293,1 − 273,1) = 13,66
Q̇L = 13,66·54,6 = 745,6 kW
Odpowiedź
Q̇L = 745,6 kW
23.Obiegi silników cieplnych gazowych
Odpowiedź: qin = 1343 kJ/kg
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)).
Temperatury
T2 = 310·140,4 = 890,9 K
T3 = 2,5·890,9 = 2227 K
Ciepło doprowadzone przy stałym ciśnieniu
qin = 1,005·(2227 − 890,9) = 1343 kJ/kg
Odpowiedź
qin = 1343 kJ/kg
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ź: v = 9,129 m/s
Continuity and Bernoulli: for an incompressible fluid the flow rate is the same at every section, Q = A₁v₁ = A₂v₂. Along a streamline without losses, p + ½ρv² + ρgz is constant, so where the fluid is faster or higher its pressure is lower.
Równanie Bernoulliego do punktu spiętrzenia
Przepływ zatrzymuje się przy końcówce, więc jego energia kinetyczna ujawnia się jako dodatkowe ciśnienie.
p0 − p = ½ρv²
v = √(2·50/1,2) = 9,129 m/s
Odpowiedź
v = 9,129 m/s
26.Ciągłość i Bernoulli
Odpowiedź: Q = 5,038 L/s
Continuity and Bernoulli: for an incompressible fluid the flow rate is the same at every section, Q = A₁v₁ = A₂v₂. Along a streamline without losses, p + ½ρv² + ρgz is constant, so where the fluid is faster or higher its pressure is lower.
Prędkość w wężu
Powierzchnia i wylot są otwarte do atmosfery, więc wyrazy ciśnieniowe się znoszą.
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,018·(800/0,3)·2,2²/(2·9,81) = 11,84 m
Spadek ciśnienia
Rura jest pozioma, a prędkość się nie zmienia, więc równanie energii upraszcza się do Δp = ρ·g·hf.
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,06/0,03142 = 1,91 m/s
V²/(2g) = 0,1859 m
Współczynniki strat
f·L/D = 0,02·100/0,2 = 10
ΣK = 0,5 + 2·0,3 + 10 = 11,1
Równanie energii
Strumień na wylocie nadal niesie swoją wysokość prędkości; to nie strata, ale trzeba ją dostarczyć.
H = (1 + 10 + 11,1)·0,1859 = 4,109 m
Odpowiedź
H = 4,109 m
29.Pompy i moc pomp
Odpowiedź: Pel = 114,3 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,1·80 = 78,48 kW
Moc na wale
Pshaft = 78,48/0,78 = 100,6 kW
Moc elektryczna
Pel = 100,6/0,88 = 114,3 kW
Odpowiedź
Pel = 114,3 kW
30.Pompy i moc pomp
Odpowiedź: H = 33,58 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.
Oba zwierciadła są otwarte i praktycznie w spoczynku, więc człony ciśnienia i prędkości znikają.
H = Δz + hL = 20 + 13,58 = 33,58 m
Odpowiedź
H = 33,58 m
31.Przewodzenie ciepła
Odpowiedź: Ti = 14,51 °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,96·4) = 0,01302 K/W
R2 = L2/(k2A) = (25/1000)/(0,17·4) = 0,03676 K/W
Strumień ciepła
Q̇ = (T1 − T3)/(R1 + R2) = 21/0,04979 = 421,8 W
Temperatura na styku
Na pierwszej warstwie temperatura spada o Q̇·R₁.
Ti = T1 − Q̇·R1 = 20 − 421,8·0,01302 = 14,51 °C
Odpowiedź
Ti = 14,51 °C
32.Przewodzenie ciepła
Odpowiedź: Q̇ = 9,299 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 szeregowo
Ri = 1/(hiA) = 1/(8·4) = 0,03125 K/W
R1 = L1/(k1A) = (250/1000)/(0,035·4) = 1,786 K/W
R2 = L2/(k2A) = (75/1000)/(0,17·4) = 0,1103 K/W
Ro = 1/(hoA) = 1/(30·4) = 0,008333 K/W
ΣR = 1,936 K/W
Straty ciepła
Q̇ = (21 − (3))/1,936 = 9,299 W
Odpowiedź
Q̇ = 9,299 W
33.Koła i przekładnie zębate
Odpowiedź: T1 = 9,879 N·m
Przekładnie zębate: zazębione koła mają tę samą prędkość na zębach, więc n₁·N₁ = n₂·N₂: koło o mniejszej liczbie zębów obraca się szybciej. Koła na jednym wale obracają się z tą samą prędkością obrotową. Moc jest przekazywana dalej (minus straty), więc moment obrotowy rośnie, gdy prędkość obrotowa maleje: T₂ = T₁·(N₂/N₁)·η, a P = T·ω, gdzie ω = 2π·n/60.
Moment silnika
ω₁ = 2π·1450/60 = 151,8 rad/s
T1 = P/ω₁ = 1500/151,8 = 9,879 N·m
Odpowiedź
T1 = 9,879 N·m
34.Koła i przekładnie zębate
Odpowiedź: θ = 162,7 °
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,15) = 8,627°
Kąt opasania
θ = 180 − 2·8,627 = 162,7°
Odpowiedź
θ = 162,7°
35.Transformaty Laplace’a
Odpowiedź: A = 10
Transformata Laplace’a: F(s) = ∫₀∞ f(t)·e−st dt zamienia pochodne na mnożenie przez s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Podstawowe pary: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), a przesunięcie s → s + σ mnoży funkcję czasu przez e−σt. Aby odwrócić wymierną F(s), rozłóż ją na ułamki proste: residuum prostego bieguna p wynosi A = [(s − p)·F(s)] przy s = p. Twierdzenie o wartości końcowej: f(∞) = lims→0 s·F(s), słuszne tylko wtedy, gdy wszystkie bieguny s·F(s) leżą w lewej półpłaszczyźnie. Twierdzenie o wartości początkowej: f(0⁺) = lims→∞ s·F(s).
Rozłóż mianownik na czynniki
s² + 11s + 30 = (s + 5)(s + 6)
F(s) = A/(s + 5) + B/(s + 6)
Metoda zakrywania
Pomnóż przez (s + 5) i podstaw s = −5: składnik z B znika.
A = (5s + 35)/(s + 6) |s = −5 = 10/1 = 10
Odpowiedź
A = 10
(B = −5, więc f(t) = 10·e−5t − 5·e−6t)
36.Transformaty Laplace’a
Odpowiedź: y(0.25) = 5,147
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² + 3s + 2)·Y(s) = 5s + 15 + 16/s
Y(s) = (5s² + 15s + 16)/(s(s + 1)(s + 2))
Ułamki proste
A = 16/(1·2) = 8
B = 6/((−1)·1) = −6
C = 6/((−2)·(−1)) = 3
Wynik
y(t) = 8 − 6·e−t + 3·e−2t
y(0,25) = 5,147
Odpowiedź
y(0,25) = 5,147
Sprawdzenie: y(0) = A + B + C = 5.
37.Schematy blokowe
Odpowiedź: τ = 0,1739 s
Algebra schematów blokowych: bloki połączone szeregowo mnożą się (G₁·G₂), bloki połączone równolegle dodają się (G₁ + G₂), a pętla ze ścieżką prostą G i ścieżką sprzężenia zwrotnego H redukuje się do T = G/(1 + G·H) dla ujemnego sprzężenia zwrotnego i G/(1 − G·H) dla dodatniego. Najpierw redukuj pętle wewnętrzne, potem idź na zewnątrz. Bieguny układu zamkniętego są pierwiastkami równania charakterystycznego 1 + G(s)H(s) = 0. Dla stabilnej pętli wymuszenie skokowe o wielkości A daje w stanie ustalonym wyjście A·T(0); przy jednostkowym sprzężeniu zwrotnym uchyb ustalony przy skoku wynosi A/(1 + Kp), gdzie Kp = G(0).
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)) = 1·2/(6 + 2K)
Rozwiąż względem K
6 + 2K = 1·2/0,06667 = 30
K = (30 − 6)/2 = 12
Odpowiedź
K = 12
39.Odpowiedź skokowa
Odpowiedź: ωn = 3,464 rad/s
Odpowiedź skokowa: układ pierwszego rzędu G(s) = K/(τs + 1) odpowiada na skok o wielkości U przebiegiem y(t) = K·U·(1 − e−t/τ): po jednej stałej czasowej τ osiąga 63,2 % wartości końcowej K·U, a po 4τ 98 %. Układ drugiego rzędu G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) ma bieguny s = −ζωₙ ± jωd, gdzie ωd = ωₙ√(1 − ζ²). Dla 0 < ζ < 1 odpowiedź skokowa ma przeregulowanie: przeregulowanie procentowe %OS = 100·e−ζπ/√(1−ζ²), czas szczytu tp = π/ωd, a 2-procentowy czas ustalania ts ≈ 4/(ζωₙ). Porównaj mianownik z s² + 2ζωₙs + ωₙ², aby odczytać ωₙ i ζ. Dla układu masa–sprężyna–tłumik ωₙ = √(k/m) i ζ = c/(2√(km)); dla szeregowego obwodu RLC ωₙ = 1/√(LC) i ζ = (R/2)·√(C/L).
Porównaj z postacią standardową
Mianownik musi odpowiadać s² + 2ζωₙs + ωₙ².
ωₙ² = 12 → ωₙ = 3,464 rad/s
2ζωₙ = 3 → ζ = 3/(2·3,464) = 0,433
Wynik
ωₙ = 3,464 rad/s
Odpowiedź
ωn = 3,464 rad/s
40.Odpowiedź skokowa
Odpowiedź: xmax = 8,632 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)) = 32/(2·√(5000·5)) = 0,1012
Wynik
xst = F/k = 25/5000 m = 5 mm
%OS = 72,65 %
xmax = 5·(1 + 0,7265) = 8,632 mm
Odpowiedź
xmax = 8,632 mm
41.Uchyb ustalony i regulacja PID
Odpowiedź: ess = 0,8333
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 = 4·3/5 = 2,4
Uchyb dla paraboli
ess = 2/2,4 = 0,8333
Odpowiedź
ess = 0,8333
42.Uchyb ustalony i regulacja PID
Odpowiedź: Kd = 3,025 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 + 4) + 4(20 + Kds) = s² + (4 + 4Kd)s + 80
Częstotliwość drgań własnych
ωₙ = √80 = 8,944 rad/s
Wzmocnienie członu różniczkującego
4 + 4Kd = 2·0,9·8,944 = 16,1
Kd = 3,025 s
Odpowiedź
Kd = 3,025 s
Człon różniczkujący zwiększa tłumienie bez zmiany ωₙ.
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.