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 ma podporę nieprzesuwną w A i przesuwną w B, w odległości 3 m od A, i wystaje poza B aż do swobodnego końca C. Wyznacz reakcję pionową RA, przyjmując zwrot w górę jako dodatni.
AB = 3 mAC = 4 mP1 = 10 kN w 0,5 mP2 = 55 kN w 4 m
RA = kN
2.Reakcje podpór belekTrudny
Belka ma podporę nieprzesuwną w A i przesuwną w B, w odległości 5 m od A, i wystaje poza B aż do swobodnego końca C. Przenosi siłę skupioną P w odległości 4 m od A oraz obciążenie równomiernie rozłożone w od 5 m do swobodnego końca w 8 m. Wyznacz reakcję pionową RA, przyjmując zwrot w górę jako dodatni.
AB = 5 mAC = 8 mP = 20 kN w 4 mw = 2 kN/m od 5 do 8 m
RA = kN
3.Naprężenie i odkształcenie osioweŚredni
Wyznacz odkształcenie osiowe ε w pręcie (aluminium), w mikroodkształceniach (µε = 10⁻⁶).
P = 110 kNd = 10 mmL = 0,5 mE = 70 GPa
ε = µε
4.Naprężenie i odkształcenie osioweTrudny
Pręt (materiał: stal, średnica 40 mm, E = 200 GPa, naprężenie dopuszczalne 250 MPa) znajduje się wewnątrz rury (materiał: miedź, średnica zewnętrzna 120 mm, średnica wewnętrzna 60 mm, E = 120 GPa, naprężenie dopuszczalne 50 MPa). Są razem ściskane między sztywnymi płytami. Wyznacz największe obciążenie Pmax.
pręt: Ø40 mm, E = 200 GPa, σallow = 250 MParura: Ø120/Ø60 mm, E = 120 GPa, σallow = 50 MPa
Pmax = kN
5.Skręcanie wałówŚredni
Pełny wał (mosiądz) jest utwierdzony w ścianie i skręcany momentem T na swobodnym końcu. Wyznacz kąt skręcenia φ swobodnego końca w stopniach.
T = 500 N·md = 50 mmL = 0,5 mG = 39 GPa
φ = °
6.Skręcanie wałówTrudny
Pełny wał ma przenosić 50 kW przy 1000 obr/min bez przekroczenia naprężenia stycznego 50 MPa. Wyznacz najmniejszą średnicę d, która wystarczy.
P = 50 kWN = 1000 rpmτallow = 50 MPa
d = mm
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.
60 kN @ 2 m8 kN/m, 0–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
Wyznacz moment zginający Mx w przekroju x = 5 m od podpory A. Przyjmij dodatnie momenty rozciągające dolne włókna i dodatnią siłę tnącą skierowaną w górę na lewej ściance.
L = 6 mP = 40 kN w 3,5 mw = 10 kN/m od 0 do 6 mx = 5 m
Mx = kN·m
9.Koło Mohra i transformacja naprężeńŚredni
Punkt cienkiej płyty jest w płaskim stanie naprężenia: σₓ = 90 MPa, σᵧ = 50 MPa i τₓᵧ = −60 MPa, przy σz = 0. Wyznacz bezwzględne maksymalne naprężenie styczne τabs w tym punkcie (z uwzględnieniem wszystkich płaszczyzn w trzech wymiarach).
σx = 90 MPaσy = 50 MPaτxy = −60 MPaσz = 0
τabs = MPa
10.Koło Mohra i transformacja naprężeńTrudny
W pewnym punkcie obciążonej części występują σₓ = −70 MPa, σᵧ = 10 MPa i τₓᵧ = −40 MPa (rozciąganie dodatnie). Wyznacz naprężenie styczne τθ na płaszczyźnie, której normalna jest obrócona o 20° przeciwnie do ruchu wskazówek zegara względem osi x (dodatnie w tym samym sensie co τₓᵧ).
σx = −70 MPaσy = 10 MPaτxy = −40 MPaθ = 20°
τθ = MPa
11.Ugięcie belkiŚredni
Belka swobodnie podparta przenosi obciążenie równomiernie rozłożone w na całej rozpiętości. Wyznacz maksymalne ugięcie δmax w środku rozpiętości.
w = 20 kN/mL = 5 mE = 70 GPa (aluminium)I = 160 × 10⁶ mm⁴
δmax = mm
12.Ugięcie belkiTrudny
Belka swobodnie podparta przenosi obciążenie równomiernie rozłożone w i siłę skupioną P w środku rozpiętości. Korzystając z superpozycji, wyznacz ugięcie w środku δmax.
P = 18 kNw = 12 kN/mL = 4 mE = 70 GPa (aluminium)I = 84 × 10⁶ mm⁴
δmax = 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 = 2 mE = 200 GPa (stal)d = 60 mm
σcr = MPa
15.Rzut ukośnyŚredni
Piłka zostaje wyrzucona pod kątem 25° do poziomu z płaskiego terenu. Pomiń opór powietrza. Wyznacz największą wysokość H, jaką osiąga piłka.
v0 = 16 m/sθ = 25°g = 9,81 m/s²
H = m
16.Rzut ukośnyTrudny
Piłkę rzucono z poziomu ziemi z prędkością 32 m/s pod kątem 65° do poziomu. Ląduje, już opadając, na płaskim dachu budynku o wysokości 9 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₀ = 32 m/sθ = 65°hroof = 9 mg = 9,81 m/s²
x = m
17.Praca i energiaŚredni
Klocek o masie 2,5 kg zsuwa się na odcinku 19 m po zboczu nachylonym pod kątem 25° do poziomu, przy współczynniku tarcia kinetycznego 0,15. Wyznacz pracę Wf wykonaną przeciw tarciu.
m = 2,5 kgθ = 25°μ = 0,15d = 19 mg = 9,81 m/s²
Wf = J
18.Praca i energiaTrudny
Klocek o masie 2 kg leży na chropowatym poziomym stole (μ = 0,25). Lekka nić biegnie od niego przez gładki bloczek na krawędzi stołu do wiszącego klocka o masie 6 kg. Układ zostaje puszczony ze spoczynku. Wyznacz prędkość v klocków po opadnięciu wiszącego klocka o 3 m.
m1 = 2 kgm2 = 6 kgμ = 0,25d = 3 mg = 9,81 m/s²
v = m/s
19.Drgania mechaniczneŚredni
Masa 2 kg jest podparta na sprężynie o sztywności 12 kN/m i tłumiku wiskotycznym o współczynniku tłumienia 60 N·s/m. Wyznacz współczynnik tłumienia względnego ζ.
m = 2 kgk = 12 kN/mc = 60 N·s/m
ζ =
20.Drgania mechaniczneTrudny
Masa 25 kg wisi na sprężynie k₃ = 60 kN/m. Górny koniec k₃ jest przymocowany do sztywnej belki, która wisi na dwóch sprężynach obok siebie, k₁ = 10 kN/m i k₂ = 50 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 = 10 kN/mk2 = 50 kN/mk3 = 60 kN/m
fn = Hz
21.Silniki cieplne i chłodziarkiŚredni
Silnik cieplny pobiera ciepło ze źródła o temperaturze 700 °C i oddaje ciepło do otoczenia o temperaturze 20 °C. Jego rzeczywista sprawność cieplna wynosi 28 %. Wyznacz jego sprawność egzergetyczną ηII, czyli stosunek sprawności rzeczywistej do sprawności Carnota.
TH = 700 °CTL = 20 °Cη = 28%
ηII = %
22.Silniki cieplne i chłodziarkiTrudny
Lodówka utrzymuje w komorze chłodniczej -20 °C w kuchni o temperaturze 28 °C, odbierając ciepło z komory w ilości 3 kW. Wyznacz minimalną moc napędową Ẇmin.
TH = 28 °CTL = -20 °CQ̇L = 3 kW
Ẇmin = kW
23.Obiegi silników cieplnych gazowychŚredni
Idealny obieg Diesla ma stopień sprężania r = 17 i stopień obciążenia rc = v₃/v₂ = 3. Przyjmij założenia modelu zimnego powietrza: powietrze jest gazem doskonałym o stałych ciepłach właściwych, k = 1,4, cp = 1,005 kJ/(kg·K) i cv = 0,718 kJ/(kg·K). Wyznacz jego sprawność cieplną η w procentach.
r = 17rc = 3k = 1,4
η = %
24.Obiegi silników cieplnych gazowychTrudny
Idealny obieg Braytona pracuje przy stosunku ciśnień w sprężarce 20. Przyjmij założenia modelu zimnego powietrza: powietrze jest gazem doskonałym o stałych ciepłach właściwych, k = 1,4, cp = 1,005 kJ/(kg·K) i cv = 0,718 kJ/(kg·K). Wyznacz jego sprawność cieplną η w procentach.
rp = 20k = 1,4
η = %
25.Ciągłość i BernoulliŚredni
Woda wypływa z dużego otwartego zbiornika przez otwór o ostrej krawędzi i średnicy 30 mm. Poziom wody utrzymuje się 4 m nad otworem, a współczynnik wypływu wynosi 0,62. Wyznacz natężenie przepływu Q. Przyjmij g = 9,81 m/s².
d = 30 mmh = 4 mCd = 0,62g = 9,81 m/s²
Q = L/s
26.Ciągłość i BernoulliTrudny
Woda płynie w górę przez rurę, która zwęża się z 100 mm do 40 mm i wznosi się o 7 m. W przekroju 1 ciśnienie wynosi 150 kPa, a prędkość 0,75 m/s. Pomijając straty, wyznacz ciśnienie p2 w przekroju 2.
Woda płynie z prędkością średnią 2,8 m/s przez rurę o długości 100 m i średnicy wewnętrznej 200 mm. Przyjmij współczynnik tarcia Darcy’ego f = 0,028 i g = 9,81 m/s². Wyznacz stratę wysokości na tarcie hf.
Dwa otwarte zbiorniki są połączone rurą o długości 200 m, średnicy wewnętrznej 250 mm i współczynniku tarcia Darcy’ego f = 0,02. Zwierciadło wody w górnym zbiorniku leży 15 m wyżej niż w dolnym. Straty miejscowe to: wlot o ostrej krawędzi (K = 0,5), kolana w liczbie 3 (K = 0,3 każde), zawór (K = 10) i wylot do dolnego zbiornika (K = 1,0). Przyjmij g = 9,81 m/s². Wyznacz natężenie przepływu Q w L/s.
Podczas próby pompa tłoczy 10 L/s wody (ρ = 1000 kg/m³) przy wysokości podnoszenia 68 m, a jej wał pobiera 9,3 kW. Przyjmij g = 9,81 m/s². Wyznacz sprawność pompy η w procentach.
Q = 10 L/sH = 68 mPshaft = 9,3 kWg = 9,81 m/s²
η = %
30.Pompy i moc pompTrudny
Pompa przetłacza 60 L/s wody (ρ = 1000 kg/m³) z otwartego dolnego zbiornika do górnego zbiornika, którego zwierciadło wody leży 11 m wyżej. Rura ma długość 100 m, średnicę wewnętrzną 200 mm i współczynnik tarcia Darcy’ego f = 0,022; współczynniki strat miejscowych sumują się do ΣK = 4 (łącznie z wlotem i wylotem). Przyjmij g = 9,81 m/s². Górny zbiornik jest zamknięty, a powietrze nad wodą ma nadciśnienie 250 kPa. Wyznacz wysokość podnoszenia H, którą musi zapewnić pompa.
Ściana składa się z 100 mm (płyta gipsowo-kartonowa) i 25 mm (drewno), a jej pole wynosi 12 m². Powierzchnie zewnętrzne mają 22 °C i -3 °C. Wyznacz strumień ciepła Q̇ przez ścianę.
L1 = 100 mm, k1 = 0,17 W/(m·K)L2 = 25 mm, k2 = 0,17 W/(m·K)A = 12 m²T1 = 22 °CT3 = -3 °C
Q̇ = W
32.Przewodzenie ciepłaTrudny
Ściana z 100 mm (wełna mineralna) i 100 mm (styropian) oddziela pomieszczenie o temperaturze 20 °C od powietrza zewnętrznego o temperaturze -5 °C. Współczynniki przejmowania ciepła wynoszą 6 W/(m²·K) wewnątrz i 20 W/(m²·K) na zewnątrz, a pole ściany 10 m². Wyznacz straty ciepła Q̇ przez ścianę.
L1 = 100 mm, k1 = 0,04 W/(m·K)L2 = 100 mm, k2 = 0,035 W/(m·K)A = 10 m²hi = 6 W/(m²·K)ho = 20 W/(m²·K)Troom = 20 °CTout = -5 °C
Q̇ = W
33.Koła i przekładnie zębateŚredni
Silnik o mocy 3 kW pracujący z prędkością 1440 rpm napędza zębnik o 15 zębach, który zazębia się z kołem o 72 zębach. Sprawność pary kół wynosi 97%. Wyznacz moment wyjściowy T2 na wale koła.
P = 3 kWn1 = 1440 rpmN1 = 15N2 = 72η = 97%
T2 = N·m
34.Koła i przekładnie zębateTrudny
W przekładni planetarnej koło słoneczne ma 16 zębów, a każdy satelita 18 zębów. Jarzmo jest unieruchomione, a koło słoneczne napędzane z prędkością 3000 rpm. Wyznacz prędkość obrotową nr koła koronowego (wartość bezwzględną).
Ns = 16Np = 18ns = 3000 rpm
nr = rpm
35.Transformaty Laplace’aŚredni
Znajdź f(t) dla F(s) = (2s + 39)/(s² + 6s + 9) i oblicz f(0,25).
F(s) = (2s + 39)/(s² + 6s + 9)t = 0,25 s
f(0.25) =
36.Transformaty Laplace’aTrudny
Układ opisuje równanie różniczkowe y″ + 2y′ + 5y = 4u′ + 7u, z wejściem u(t) i wyjściem y(t). Znajdź jego transmitancję G(s) = Y(s)/U(s) (zerowe warunki początkowe) i oblicz odpowiedź impulsową g(t) dla t = 1,5 s. Znajdź g(1,5).
y″ + 2y′ + 5y = 4u′ + 7ut = 1,5 s
g(1.5) =
37.Schematy blokoweŚredni
Obiekt G(s) = 6/(s + 6) ma wewnętrzną (podrzędną) pętlę ujemnego sprzężenia zwrotnego ze wzmocnieniem H₁ = 1. Wzmocnienie regulatora K = 9 i zewnętrzna pętla jednostkowego ujemnego sprzężenia zwrotnego zamykają układ. Znajdź biegun układu zamkniętego p.
G(s) = 6/(s + 6)H₁ = 1K = 9
p =
38.Schematy blokoweTrudny
Regulator proporcjonalny K steruje obiektem G(s) = 1/(s + 7) w pętli jednostkowego ujemnego sprzężenia zwrotnego. Skok zakłócenia o wielkości 10 wchodzi na wejście obiektu, gdy R = 0. Znajdź najmniejsze wzmocnienie K, przy którym wyjście w stanie ustalonym nie przekracza 0,02857.
G(s) = 1/(s + 7)D = 10css ≤ 0,02857
K =
39.Odpowiedź skokowaŚredni
Układ ma transmitancję G(s) = 50/(s² + 3s + 100) i otrzymuje jednostkowe wymuszenie skokowe. Znajdź pulsację własną ωn w rad/s.
G(s) = 50/(s² + 3s + 100)
ωn = rad/s
40.Odpowiedź skokowaTrudny
Szeregowy obwód RLC z R = 175 Ω, L = 100 mH i C = 4,7 µF jest dołączony do skoku 12 V w chwili t = 0. Kondensator jest nienaładowany i przed skokiem nie płynie prąd. Znajdź szczytowe napięcie na kondensatorze vC,max.
R = 175 ΩL = 100 mHC = 4,7 µFV = 12 V
vC,max = V
41.Uchyb ustalony i regulacja PIDŚredni
Proces G(s) = 2/(0,5s + 1) jest sterowany wzmocnieniem proporcjonalnym K = 9 w pętli ze sprzężeniem jednostkowym. Wyznacz uchyb ustalony ess dla skoku, jako procent wartości skoku.
G(s) = 2/(0,5s + 1)K = 9
ess = %
42.Uchyb ustalony i regulacja PIDTrudny
Regulator proporcjonalny K = 5 steruje obiektem G(s) = 1/(s(s + 2)) w pętli ze sprzężeniem jednostkowym. Znajdź 2-procentowy czas ustalania ts odpowiedzi skokowej, korzystając z ts ≈ 4/(ζωₙ).
G(s) = 1/(s(s + 2))K = 5
ts = s
3Roztwory
1.Reakcje podpór belek
Odpowiedź: RA = −10 kN
Równowaga statyczna: ciało w spoczynku spełnia ΣFx = 0, ΣFy = 0 i ΣM = 0 względem dowolnego punktu. Zastąp każde obciążenie ciągłe jego wypadkową (polem wykresu obciążenia) przyłożoną w środku ciężkości wykresu obciążenia i licz momenty względem podpory, aby usunąć z równania reakcję tej podpory.
Momenty względem A
Przeciwnie do ruchu wskazówek zegara dodatnio.
ΣMA = 0: RB·3 = 10·0,5 + 55·4
RB = 225/3 = 75 kN
Równowaga sił pionowych
RA = 10 + 55 − 75 = −10 kN
Odpowiedź
RA = −10 kN
RA jest ujemne: podpora w A musi przytrzymywać belkę.
2.Reakcje podpór belek
Odpowiedź: RA = 2,2 kN
Równowaga statyczna: ciało w spoczynku spełnia ΣFx = 0, ΣFy = 0 i ΣM = 0 względem dowolnego punktu. Zastąp każde obciążenie ciągłe jego wypadkową (polem wykresu obciążenia) przyłożoną w środku ciężkości wykresu obciążenia i licz momenty względem podpory, aby usunąć z równania reakcję tej podpory.
Wypadkowa obciążenia ciągłego
W = 2·3 = 6 kN, x̄ = (5 + 8)/2 = 6,5 m
Równowaga
Momenty względem A, dodatnie przeciwnie do ruchu wskazówek zegara. Potem równowaga sił pionowych.
ΣMA = 0: RB·5 = 20·4 + 6·6,5
RB = 119/5 = 23,8 kN
RA = 26 − 23,8 = 2,2 kN
Odpowiedź
RA = 2,2 kN
3.Naprężenie i odkształcenie osiowe
Odpowiedź: ε = 20008 µε
Prawo Hooke’a: w zakresie sprężystym σ = E·ε. Łącząc to z σ = P/A i ε = δ/L, dostajemy δ = PL/(AE).
Naprężenie
A = π·10²/4 = 78,54 mm²
σ = P/A = 110·1000/78,54 = 1401 MPa
Odkształcenie (prawo Hooke’a)
E for aluminium = 70 GPa = 70000 MPa.
ε = σ/E = 1401/70000 = 0,02001 = 20008 µε
Odpowiedź
ε = 20008 µε
4.Naprężenie i odkształcenie osiowe
Odpowiedź: Pmax = 528,8 kN
Dwa materiały: przy połączeniu szeregowym (jeden za drugim) każdy odcinek przenosi tę samą siłę, a wydłużenia się dodają: δ = P·Σ L/(E·A). Przy połączeniu równoległym (obok siebie między sztywnymi płytami) obie części skracają się o tyle samo, więc ich odkształcenia są równe: sztywniejsza część przejmuje większy udział, Pi = P·EiAi/ΣEA, a σi = Ei·ε. Pręt zamocowany między dwiema sztywnymi ścianami z obciążeniem w miejscu połączenia jest statycznie niewyznaczalny: dwa odcinki dzielą obciążenie proporcjonalnie do sztywności k = EA/L.
Skręcanie wałów okrągłych: naprężenie styczne rośnie liniowo od zera w środku do τmax = T·r/J na powierzchni, gdzie J = πd⁴/32 dla wału pełnego i J = π(do⁴ − di⁴)/32 dla drążonego. Kąt skręcenia to φ = TL/(GJ) w radianach. Moc i moment skręcający wiąże P = T·ω.
Biegunowy moment bezwładności
J = π·50⁴/32 = 613592 mm⁴
Kąt skręcenia
φ = TL/(GJ) = 500·1000·500/(39000·613592) = 0,01045 rad
φ = 0,01045·180/π = 0,5986°
Odpowiedź
φ = 0,5986°
6.Skręcanie wałów
Odpowiedź: d = 36,5 mm
Skręcanie wałów okrągłych: naprężenie styczne rośnie liniowo od zera w środku do τmax = T·r/J na powierzchni, gdzie J = πd⁴/32 dla wału pełnego i J = π(do⁴ − di⁴)/32 dla drążonego. Kąt skręcenia to φ = TL/(GJ) w radianach. Moc i moment skręcający wiąże P = T·ω.
Moment skręcający
ω = 2πN/60 = 2π·1000/60 = 104,7 rad/s
T = P/ω = 50000/104,7 = 477,5 N·m
Średnica z naprężenia dopuszczalnego
d = ∛(16T/(π·τallow)) = ∛(16·477465/(π·50)) = 36,5 mm
Odpowiedź
d = 36,5 mm
W praktyce zaokrągla się w górę do najbliższego wymiaru znormalizowanego.
7.Siła tnąca i moment zginający
Odpowiedź: C (x = 2 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 = 64 kN, R2 = 44 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 = 60, MC = 112 kN·m, MD = 96, ME = 72, MF = 40, MG = 0
Odpowiedź
Punkt C: MC = 112 kN·m
8.Siła tnąca i moment zginający
Odpowiedź: Mx = 48,33 kN·m
Metoda przekrojów: przetnij belkę tam, gdzie chcesz znać siły wewnętrzne, i zostaw jedną stronę. Siła tnąca V równoważy siły pionowe po tej stronie, a moment zginający M ich momenty względem przekroju.
Reakcje podpór
Równowaga całej belki.
Wypadkowa obciążenia ciągłego: W = w·ℓ = 10·6 = 60 kN, przyłożona w x̄ = 3 m
ΣMA = 0: RB·6 = 40·3,5 + 10·6·3
RB = 320/6 = 53,33 kN
ΣFy = 0: RA = 100 − 53,33 = 46,67 kN
Przetnij w x = 5 m i zostaw lewą część
Uwzględniasz tylko siły na lewo od przekroju. Fragment obciążenia ciągłego działa jak jedna siła w środku tego fragmentu.
M = RA·x − (każde obciążenie na lewo od przekroju × jego odległość od przekroju)
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. W trzech wymiarach bezwzględne maksymalne naprężenie styczne wynosi (σmax − σmin)/2 ze wszystkich trzech naprężeń głównych; w płaskim stanie naprężenia trzecie jest równe zeru.
Tutaj decyduje koło spoza płaszczyzny: maksimum bezwzględne jest większe niż maksimum w płaszczyźnie R.
10.Koło Mohra i transformacja naprężeń
Odpowiedź: τθ = −4,93 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
τθ = −−40·0,6428 + −40·0,766 = −4,93 MPa
Odpowiedź
τθ = −4,93 MPa
11.Ugięcie belki
Odpowiedź: δmax = 14,53 mm
Ugięcie belek: ugięcia otrzymuje się, całkując dwukrotnie krzywiznę M/(EI). Typowe przypadki są stablicowane, na przykład PL³/(3EI) dla wspornika z siłą na końcu i 5wL⁴/(384EI) dla belki swobodnie podpartej z obciążeniem równomiernym na całej rozpiętości. Ugięcia od różnych obciążeń się sumują (superpozycja).
Spójne jednostki
Licz w N i mm: P w N, w w N/mm (kN/m = N/mm), L w mm, E w N/mm² (GPa × 1000), I w mm⁴.
w = 20 N/mm, L = 5000 mm
EI = 1,12 × 10¹³ N·mm²
Ugięcie
δmax = 5wL⁴/(384EI) = 5·20·5000⁴/(384·1,12 × 10¹³) = 14,53 mm
Odpowiedź
δmax = 14,53 mm
12.Ugięcie belki
Odpowiedź: δmax = 10,88 mm
Ugięcie belek: ugięcia otrzymuje się, całkując dwukrotnie krzywiznę M/(EI). Typowe przypadki są stablicowane, na przykład PL³/(3EI) dla wspornika z siłą na końcu i 5wL⁴/(384EI) dla belki swobodnie podpartej z obciążeniem równomiernym na całej rozpiętości. Ugięcia od różnych obciążeń się sumują (superpozycja).
Spójne jednostki
Licz w N i mm: P w N, w w N/mm (kN/m = N/mm), L w mm, E w N/mm² (GPa × 1000), I w mm⁴.
P = 18000 N, w = 12 N/mm, L = 4000 mm, EI = 5,88 × 10¹² N·mm²
Każde obciążenie osobno
δP = PL³/(48EI) = 18000·4000³/(48·5,88 × 10¹²) = 4,082 mm
δw = 5wL⁴/(384EI) = 5·12·4000⁴/(384·5,88 × 10¹²) = 6,803 mm
Dodaj
δmax = 4,082 + 6,803 = 10,88 mm
Odpowiedź
δmax = 10,88 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 = π·60⁴/64 = 636173 mm⁴
Promień bezwładności i smukłość
A = 2827 mm²
r = √(I/A) = √(636173/2827) = 15 mm
KL = 2·2000 = 4000 mm
KL/r = 4000/15 = 266,7
Naprężenie krytyczne
σcr = π²E/(KL/r)² = π²·200000/266,7² = 27,76 MPa
Odpowiedź
σcr = 27,76 MPa
15.Rzut ukośny
Odpowiedź: H = 2,33 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 θ = 16·cos 25° = 14,5 m/s
vy = v0·sin θ = 16·sin 25° = 6,762 m/s
Najwyższy punkt
Na szczycie prędkość pionowa jest zerowa: vy² = 2gH.
H = vy²/(2g) = 6,762²/(2·9,81) = 2,33 m
Odpowiedź
H = 2,33 m
16.Rzut ukośny
Odpowiedź: x = 75,52 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 = 13,52 m/s, vy = 29 m/s
Czas do lądowania na dachu
Większy pierwiastek: piłka opada.
9 = 29·t − 4,905·t²
t = (29 + √(29² − 2·9,81·9))/9,81 = 5,584 s
Odległość pozioma
x = vx·t = 13,52·5,584 = 75,52 m
Odpowiedź
x = 75,52 m
17.Praca i energia
Odpowiedź: Wf = 63,35 J
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).
Siła tarcia
N = mg·cos θ = 24,53·0,9063 = 22,23 N
F = μN = 0,15·22,23 = 3,334 N
Praca przeciw tarciu
Wf = F·d = 3,334·19 = 63,35 J
Odpowiedź
Wf = 63,35 J
18.Praca i energia
Odpowiedź: v = 6,361 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 uwolniona i stracona
m₂g·d = 58,86·3 = 176,6 J
μ·m₁g·d = 0,25·19,62·3 = 14,71 J
Energia kinetyczna obu klocków
½(m₁ + m₂)v² = 161,9 → v = √(2·161,9/8) = 6,361 m/s
Odpowiedź
v = 6,361 m/s
19.Drgania mechaniczne
Odpowiedź: ζ = 0,1936
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·√(12000·2) = 309,8 N·s/m
Współczynnik tłumienia
ζ = c/cc = 60/309,8 = 0,1936
Odpowiedź
ζ = 0,1936
ζ < 1, więc tłumienie jest podkrytyczne, a układ drga z malejącą amplitudą.
20.Drgania mechaniczne
Odpowiedź: fn = 5,513 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 = 10 + 50 = 60 kN/m
keq = 60·60/(60 + 60) = 30 kN/m
Częstotliwość drgań własnych
ωn = √(30000/25) = 34,64 rad/s
fn = 34,64/(2π) = 5,513 Hz
Odpowiedź
fn = 5,513 Hz
21.Silniki cieplne i chłodziarki
Odpowiedź: ηII = 40,07 %
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
TL = 20 + 273,15 = 293,1 K
Sprawność Carnota
ηCarnot = 1 − 293,1/973,1 = 69,88 %
Sprawność egzergetyczna
ηII = η/ηCarnot = 28/69,88 = 40,07 %
Odpowiedź
ηII = 40,07 %
22.Silniki cieplne i chłodziarki
Odpowiedź: Ẇmin = 0,5688 kW
Ograniczenia z drugiej zasady: żaden obieg pracujący między zbiornikami o temperaturach bezwzględnych TH i TL nie przewyższy obiegu odwracalnego (Carnota): ηmax = 1 − TL/TH, COPR,max = TL/(TH − TL) i COPHP,max = TH/(TH − TL). Energia jest zawsze zachowana: QH = W + QL.
Najmniej pracy potrzebuje maszyna odwracalna, która ma największy COP.
Ẇmin = Q̇L/COP = 3/5,274 = 0,5688 kW
Odpowiedź
Ẇmin = 0,5688 kW
23.Obiegi silników cieplnych gazowych
Odpowiedź: η = 57,96 %
Obieg Diesla (silnik o zapłonie samoczynnym): 1→2 sprężanie izentropowe, 2→3 doprowadzenie ciepła przy stałym ciśnieniu, 3→4 rozprężanie izentropowe, 4→1 odprowadzenie ciepła przy stałej objętości. Przy stopniu sprężania r = v₁/v₂ i stopniu obciążenia rc = v₃/v₂: T₂ = T₁·rk−1, T₃ = rc·T₂, qin = cp·(T₃ − T₂) oraz η = 1 − (1/rk−1)·(rck − 1)/(k·(rc − 1)).
Człon stopnia obciążenia
(31,4 − 1)/(1,4·(3 − 1)) = 1,306
Sprawność cieplna
η = 1 − 0,322·1,306 = 57,96%
Odpowiedź
η = 57,96%
24.Obiegi silników cieplnych gazowych
Odpowiedź: η = 57,51 %
Obieg Braytona (turbina gazowa): 1→2 sprężanie izentropowe, 2→3 doprowadzenie ciepła przy stałym ciśnieniu, 3→4 rozprężanie izentropowe w turbinie, 4→1 odprowadzenie ciepła przy stałym ciśnieniu. Przy stosunku ciśnień rp = P₂/P₁: T₂ = T₁·rp(k−1)/k, T₄ = T₃/rp(k−1)/k oraz η = 1 − 1/rp(k−1)/k. Praca netto w = cp·[(T₃ − T₄) − (T₂ − T₁)]. Przy sprawnościach izentropowych sprężarka wymaga ws/ηc, a turbina daje ηt·ws.
Sprawność cieplna
η = 1 − 1/200,2857 = 1 − 1/2,354 = 57,51%
Odpowiedź
η = 57,51%
25.Ciągłość i Bernoulli
Odpowiedź: Q = 3,882 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.
Idealna prędkość strumienia
√(2gh) = √(2·9,81·4) = 8,859 m/s
Rzeczywiste natężenie przepływu
Współczynnik wypływu uwzględnia kontrakcję strumienia i tarcie na krawędzi.
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.
Przepływ w rurach: liczba Reynoldsa Re = V·D/ν = ρ·V·D/μ odróżnia przepływ laminarny (Re poniżej ok. 2000) od turbulentnego (powyżej ok. 4000). Tarcie wzdłuż rury powoduje stratę wysokości Darcy’ego–Weisbacha hf = f·(L/D)·V²/(2g). Kształtki, zawory, wlot i wylot dodają straty miejscowe hm = ΣK·V²/(2g). Między dwiema otwartymi powierzchniami zbiorników równanie energii daje Δz = hf + hm, z uwzględnioną stratą na wylocie K = 1.
Wysokość prędkości
V²/(2g) = 2,8²/(2·9,81) = 0,3996 m
Darcy–Weisbach
hf = f·(L/D)·V²/(2g) = 0,028·(100/0,2)·2,8²/(2·9,81) = 5,594 m
Odpowiedź
hf = 5,594 m
28.Tarcie w rurach i straty wysokości
Odpowiedź: Q = 158 L/s
Przepływ w rurach: liczba Reynoldsa Re = V·D/ν = ρ·V·D/μ odróżnia przepływ laminarny (Re poniżej ok. 2000) od turbulentnego (powyżej ok. 4000). Tarcie wzdłuż rury powoduje stratę wysokości Darcy’ego–Weisbacha hf = f·(L/D)·V²/(2g). Kształtki, zawory, wlot i wylot dodają straty miejscowe hm = ΣK·V²/(2g). Między dwiema otwartymi powierzchniami zbiorników równanie energii daje Δz = hf + hm, z uwzględnioną stratą na wylocie K = 1.
Współczynniki strat
f·L/D = 0,02·200/0,25 = 16
ΣK = 0,5 + 3·0,3 + 10 + 1,0 = 12,4
Równanie energii
Obie powierzchnie są otwarte do atmosfery i praktycznie w spoczynku, więc całe Δz jest tracone w rurze.
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,01·68 = 6,671 kW
Sprawność
η = Ph/Pshaft = 6,671/9,3 = 71,73%
Odpowiedź
η = 71,73%
30.Pompy i moc pomp
Odpowiedź: H = 39,27 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ą praktycznie w spoczynku, więc nie ma członu wysokości prędkości.
H = Δz + p/(ρ·g) + hL = 11 + 25,48 + 2,789 = 39,27 m
Odpowiedź
H = 39,27 m
31.Przewodzenie ciepła
Odpowiedź: Q̇ = 408 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.
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.
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π·1440/60 = 150,8 rad/s
T1 = P/ω₁ = 3000/150,8 = 19,89 N·m
Moment wyjściowy
T2 = 19,89·72/15·0,97 = 92,63 N·m
Odpowiedź
T2 = 92,63 N·m
34.Koła i przekładnie zębate
Odpowiedź: nr = 923,1 rpm
Przekładnie planetarne: prędkości koła słonecznego, koronowego i jarzma wiąże wzór Willisa (ns − nc)/(nr − nc) = −Nr/Ns: widziana z jarzma przekładnia jest zwykłą przekładnią z satelitami jako kołami pośredniczącymi. Unieruchomienie jednego członu zamienia ją w reduktor o jednym przełożeniu.
Liczba zębów koła koronowego
Nr = Ns + 2·Np = 16 + 2·18 = 52
Wzór Willisa przy nieruchomym jarzmie
Przy nc = 0 koło słoneczne napędza koło koronowe przez satelity, które działają jak koła pośredniczące.
nr = −ns·Ns/Nr = −3000·16/52 = −923,1 rpm
|nr| = 923,1 rpm
Odpowiedź
nr = 923,1 rpm
Koło koronowe obraca się przeciwnie do koła słonecznego.
35.Transformaty Laplace’a
Odpowiedź: f(0.25) = 4,842
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² + 6s + 9 = (s + 3)²
Ułamki proste
2s + 39 = 2(s + 3) + 33
F(s) = 2/(s + 3) + 33/(s + 3)²
Transformata odwrotna
f(t) = 2·e−3t + 33·t·e−3t
f(0,25) = (2 + 33·0,25)·e−3·0,25 = 4,842
Odpowiedź
f(0,25) = 4,842
36.Transformaty Laplace’a
Odpowiedź: g(1.5) = −0,8364
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).
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)) = 10·1/(7 + 1K)
Rozwiąż względem K
7 + 1K = 10·1/0,02857 = 350
K = (350 − 7)/1 = 343
Odpowiedź
K = 343
39.Odpowiedź skokowa
Odpowiedź: ωn = 10 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 + ωₙ².
ωₙ² = 100 → ωₙ = 10 rad/s
2ζωₙ = 3 → ζ = 3/(2·10) = 0,15
Wynik
ωₙ = 10 rad/s
Odpowiedź
ωn = 10 rad/s
40.Odpowiedź skokowa
Odpowiedź: vC,max = 13,14 V
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).
Uchyb ustalony: w pętli ze sprzężeniem jednostkowym o transmitancji układu otwartego G(s) uchyb wynosi E(s) = R(s)/(1 + G(s)), a z twierdzenia o wartości końcowej ess = lims→0 s·E(s). Typ układu to liczba biegunów G(s) w s = 0. Stałe uchybowe to Kp = lim G(s), Kv = lim s·G(s) oraz Ka = lim s²·G(s), wszystkie dla s → 0. Skok o wartości A daje ess = A/(1 + Kp), narastanie liniowe A·t daje A/Kv, a parabola A·t²/2 daje A/Ka. Pętla typu 0 nadąża za skokiem ze skończonym uchybem, pętla typu 1 nadąża za nim dokładnie, a za narastaniem liniowym ze skończonym uchybem.
Stała położeniowa
Kp = 9·2 = 18
Uchyb skokowy
ess = 100/(1 + 18) = 5,263 %
Odpowiedź
ess = 5,263 %
42.Uchyb ustalony i regulacja PID
Odpowiedź: ts = 4 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 + 2) + 5·1 = s² + 2s + 5
ωₙ i ζ
ωₙ = √5 = 2,236 rad/s
ζ = 2/(2·2,236) = 0,4472
Czas ustalania
ζωₙ = 2/2 = 1
ts ≈ 4/1 = 4 s
Odpowiedź
ts = 4 s
Czas ustalania nie zależy tu od K: zwiększenie K tylko przyspiesza oscylacje odpowiedzi.
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.