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 ma podporę przegubową nieprzesuwną w A i przesuwną w B. Wyznacz pionową reakcję podpory RA.
L = 4 mP1 = 40 kN w 1,5 mP2 = 25 kN w 2,5 m
RA = kN
2.Reakcje podpór belekTrudny
Belka złożona jest utwierdzona w ścianie w A, ma przegub wewnętrzny w C (3 m od A) i opiera się na podporze przesuwnej w D (6 m od A). Przenosi siłę skupioną P w odległości 5 m od A, obciążenie równomiernie rozłożone w na całej części CD. Wyznacz reakcję RD na podporze przesuwnej.
AC = 3 mAD = 6 mP = 25 kN w 5 mw = 3 kN/m od 3 do 6 m
RD = 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ł: miedź, średnica zewnętrzna 80 mm, średnica wewnętrzna 50 mm, E = 120 GPa). Oba mają długość 800 mm i są razem ściskane między sztywnymi płytami siłą 50 kN. Wyznacz skrócenie δ.
P = 50 kNpręt: Ø40 mm, E = 200 GParura: Ø80/Ø50 mm, E = 120 GPaL = 800 mm
δ = mm
4.Naprężenie i odkształcenie osioweTrudny
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 = 20 kNodcinek 1: Ø40 mm, 2 m, aluminium (E = 70 GPa)odcinek 2: Ø20 mm, 2 m, aluminium (E = 70 GPa)
σmax = MPa
5.Skręcanie wałówŚredni
Pełny wał schodkowy (stal) jest utwierdzony w A. Odcinek AB ma średnicę 40 mm i długość 1 m; odcinek BC ma średnicę 25 mm i długość 0,5 m. Na swobodnym końcu C działa moment skręcający 1500 N·m. Wyznacz kąt skręcenia φC końca C w stopniach.
T = 1500 N·md1 = 40 mm, L1 = 1 md2 = 25 mm, L2 = 0,5 mG = 80 GPa
φC = °
6.Skręcanie wałówTrudny
Pełny wał (mosiądz) przenosi moc 20 kW przy 300 rpm. Naprężenie styczne nie może przekroczyć 50 MPa, a kąt skręcenia 0,75° na m długości. Wyznacz najmniejszą średnicę d, która spełnia oba warunki.
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.
20 kN @ 5 m50 kN @ 4 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
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.
20 kN @ 12 m25 kN @ 2 m
1A (x = 0 m)
2B (x = 2 m)
3C (x = 4 m)
4D (x = 6 m)
5E (x = 8 m)
6F (x = 10 m)
7G (x = 12 m)
9.Koło Mohra i transformacja naprężeńŚredni
W pewnym punkcie obciążonej części występują σₓ = 100 MPa, σᵧ = 0 MPa i τₓᵧ = 10 MPa (rozciąganie dodatnie). Wyznacz kąt θp od osi x do kierunku największego naprężenia głównego, z przedziału od −90° do 90°.
σx = 100 MPaσy = 0 MPaτxy = 10 MPa
θp = °
10.Koło Mohra i transformacja naprężeńTrudny
W pewnym punkcie obciążonej części występują σₓ = 100 MPa, σᵧ = 80 MPa i τₓᵧ = 20 MPa (rozciąganie dodatnie). Wyznacz naprężenie styczne τθ na płaszczyźnie, której normalna jest obrócona o 25° przeciwnie do ruchu wskazówek zegara względem osi x (dodatnie w tym samym sensie co τₓᵧ).
σx = 100 MPaσy = 80 MPaτxy = 20 MPaθ = 25°
τθ = 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 kąt obrotu θA belki na podporze A, w radianach.
w = 12 kN/mL = 8 mE = 200 GPa (stal)I = 120 × 10⁶ mm⁴
θA = rad
12.Ugięcie belkiTrudny
Belka drewniana swobodnie podparta o przekroju prostokątnym i szerokości 125 mm przenosi obciążenie równomiernie rozłożone w na całej rozpiętości. Jej ugięcie w środku rozpiętości nie może przekroczyć rozpiętość/360. Wyznacz najmniejszą wysokość hreq przekroju.
w = 2 kN/mL = 4,5 mE = 11 GPa (drewno konstrukcyjne)b = 125 mm
hreq = mm
13.Wyboczenie słupówŚredni
Słup o długości L z materiału (stal) jest przegubowo podparty na obu końcach. Korzystając ze wzoru Eulera, wyznacz jego siłę krytyczną Pcr.
L = 4 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 = 4 mE = 200 GPa (stal)d = 80 mm
σcr = MPa
15.Rzut ukośnyŚredni
Piłka zostaje wyrzucona pod kątem 45° do poziomu z płaskiego terenu. Pomiń opór powietrza. Wyznacz największą wysokość H, jaką osiąga piłka.
v0 = 14 m/sθ = 45°g = 9,81 m/s²
H = m
16.Rzut ukośnyTrudny
Piłkę wyrzucono z prędkością 20 m/s pod kątem 15° do poziomu ze szczytu klifu o wysokości 55 m. Pomiń opór powietrza i przyjmij g = 9,81 m/s². Pod jakim kątem α poniżej poziomu uderza w ziemię?
v₀ = 20 m/sθ = 15°h = 55 mg = 9,81 m/s²
α = °
17.Praca i energiaŚredni
Klocek o masie 5 kg zsuwa się na odcinku 3 m po zboczu nachylonym pod kątem 30° do poziomu, przy współczynniku tarcia kinetycznego 0,3. Wyznacz pracę Wf wykonaną przeciw tarciu.
m = 5 kgθ = 30°μ = 0,3d = 3 mg = 9,81 m/s²
Wf = J
18.Praca i energiaTrudny
Klocek o masie 5 kg puszczono bez prędkości początkowej z wysokości 4,5 m na gładkiej rampie. Na dole przejeżdża przez szorstki poziomy pas o długości 2,5 m (μ = 0,15), a potem trafia na sprężynę (k = 1000 N/m) na gładkiej powierzchni. Wyznacz maksymalne ugięcie x sprężyny.
m = 5 kgh = 4,5 mL = 2,5 mμ = 0,15k = 1000 N/mg = 9,81 m/s²
x = cm
19.Drgania mechaniczneŚredni
Masa 200 kg jest podparta na sprężynie o sztywności 20 kN/m i tłumiku wiskotycznym o współczynniku tłumienia 800 N·s/m. Wyznacz współczynnik tłumienia względnego ζ.
m = 200 kgk = 20 kN/mc = 800 N·s/m
ζ =
20.Drgania mechaniczneTrudny
Masa 50 kg wisi na dwóch sprężynach obok siebie (równolegle) o sztywnościach k₁ = 12 kN/m i k₂ = 20 kN/m. Masa pozostaje pozioma. Wyznacz sztywność zastępczą keq układu sprężyn w kN/m.
m = 50 kgk1 = 12 kN/mk2 = 20 kN/m
keq = kN/m
21.Silniki cieplne i chłodziarkiŚredni
Odwracalna chłodziarka utrzymuje w chłodni temperaturę -5 °C, odbierając 8,5 kW ciepła, i oddaje ciepło do otoczenia o temperaturze 30 °C. Wyznacz strumień ciepła oddawanego Q̇H.
TL = -5 °CTH = 30 °CQ̇L = 8,5 kW
Q̇H = kW
22.Silniki cieplne i chłodziarkiTrudny
Pompa ciepła utrzymuje w domu 20 °C, gdy na zewnątrz jest 5 °C, dostarczając do domu 7 kW ciepła. Wyznacz minimalną moc napędową Ẇmin.
TH = 20 °CTL = 5 °CQ̇H = 7 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₂ = 2,2. 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 = 17rc = 2,2T1 = 290 Kk = 1,4
T3 = K
24.Obiegi silników cieplnych gazowychTrudny
Elektrownia z turbiną gazową pracuje w obiegu Braytona przy stosunku ciśnień w sprężarce 15. Powietrze wpływa do sprężarki w temperaturze 288 K, a do turbiny w temperaturze 1400 K. Obie maszyny są izentropowe. 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 pracę netto wnet na kilogram powietrza, w kJ/kg.
rp = 15T1 = 288 KT3 = 1400 Kk = 1,4cp = 1,005 kJ/(kg·K)
wnet = kJ/kg
25.Ciągłość i BernoulliŚredni
Woda wypływa z dużego otwartego zbiornika przez otwór o ostrej krawędzi i średnicy 50 mm. Poziom wody utrzymuje się 1 m nad otworem, a współczynnik wypływu wynosi 0,63. Wyznacz natężenie przepływu Q. Przyjmij g = 9,81 m/s².
d = 50 mmh = 1 mCd = 0,63g = 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 200 mm do 100 mm i wznosi się o 5 m. W przekroju 1 ciśnienie wynosi 390 kPa, a prędkość 2,75 m/s. Pomijając straty, wyznacz ciśnienie p2 w przekroju 2.
Woda płynie z natężeniem 25 L/s przez rurę o długości 200 m i średnicy wewnętrznej 200 mm. Przyjmij współczynnik tarcia Darcy’ego f = 0,02 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 200 mm i współczynniku tarcia Darcy’ego f = 0,015. 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 40 L/s?
Pompa o sprawności 60% tłoczy 50 L/s wody (ρ = 1000 kg/m³) przy wysokości podnoszenia 60 m. Przyjmij g = 9,81 m/s². Wyznacz potrzebną moc na wale Pshaft, w kW.
Q = 50 L/sH = 60 mη = 60%g = 9,81 m/s²
Pshaft = kW
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 13 m wyżej. Rura ma długość 600 m, średnicę wewnętrzną 200 mm i współczynnik tarcia Darcy’ego f = 0,018; 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 (styropian) i 50 mm (beton), a jej pole wynosi 20 m². Powierzchnie zewnętrzne mają 21 °C i -1 °C. Wyznacz temperaturę Ti na styku dwóch warstw.
L1 = 50 mm, k1 = 0,035 W/(m·K)L2 = 50 mm, k2 = 1,4 W/(m·K)A = 20 m²T1 = 21 °CT3 = -1 °C
Ti = °C
32.Przewodzenie ciepłaTrudny
Ściana z 150 mm (styropian) i 75 mm (wełna mineralna) oddziela pomieszczenie o temperaturze 21 °C od powietrza zewnętrznego o temperaturze 0 °C. Współczynniki przejmowania ciepła wynoszą 8 W/(m²·K) wewnątrz i 30 W/(m²·K) na zewnątrz, a pole ściany 12 m². Wyznacz współczynnik przenikania ciepła U ściany.
L1 = 150 mm, k1 = 0,035 W/(m·K)L2 = 75 mm, k2 = 0,04 W/(m·K)A = 12 m²hi = 8 W/(m²·K)ho = 30 W/(m²·K)Troom = 21 °CTout = 0 °C
U = W/(m²·K)
33.Koła i przekładnie zębateŚredni
Silnik o mocy 2,2 kW pracujący z prędkością 1450 rpm napędza zębnik o 25 zębach i module 5 mm. Wyznacz siłę obwodową Ft działającą na zęby zębnika.
P = 2,2 kWn = 1450 rpmN1 = 25m = 5 mm
Ft = N
34.Koła i przekładnie zębateTrudny
W przekładni dwustopniowej koło A (15 zębów) na wale wejściowym, obracające się z prędkością 1440 rpm, napędza koło B (25 zębów). Koło C (16 zębów) jest osadzone na tym samym wale co B i napędza koło D (30 zębów) na wale wyjściowym. Wyznacz prędkość obrotową na wyjściu nout.
nin = 1440 rpmA = 15B = 25C = 16D = 30
nout = rpm
35.Transformaty Laplace’aŚredni
Znajdź f(t) dla F(s) = (6s + 5)/(s² + 8s + 16) i oblicz f(0,4).
F(s) = (6s + 5)/(s² + 8s + 16)t = 0,4 s
f(0.4) =
36.Transformaty Laplace’aTrudny
Układ opisuje równanie różniczkowe y″ + 2y′ + 26y = 19u, 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 = 0,8 s. Znajdź g(0,8).
y″ + 2y′ + 26y = 19ut = 0,8 s
g(0.8) =
37.Schematy blokoweŚredni
Pętla jednostkowego ujemnego sprzężenia zwrotnego ma G(s) = 12(s + 3)/((s + 2)(s + 6)). Na wartość zadaną podano skok o wielkości 1. Znajdź uchyb ustalony ess.
G(s) = 12(s + 3)/((s + 2)(s + 6))R = 1
ess =
38.Schematy blokoweTrudny
Obiekt G(s) = K/(s(s + 2)(s + 7)) jest w pętli jednostkowego ujemnego sprzężenia zwrotnego. Za pomocą kryterium Routha–Hurwitza znajdź największe wzmocnienie Kmax, przy którym układ zamknięty jest stabilny.
Układ ma transmitancję G(s) = 54/(s² + 2s + 54) i otrzymuje jednostkowe wymuszenie skokowe. Znajdź 2-procentowy czas ustalania ts odpowiedzi skokowej, korzystając z ts ≈ 4/(ζωₙ).
G(s) = 54/(s² + 2s + 54)
ts = s
40.Odpowiedź skokowaTrudny
Pętla jednostkowego ujemnego sprzężenia zwrotnego ma G(s) = K/(s(s + 10)). Znajdź wzmocnienie K, które daje odpowiedź skokową z przeregulowaniem 15 %.
G(s) = K/(s(s + 10))%OS = 15 %
K =
41.Uchyb ustalony i regulacja PIDŚredni
Na pętlę ze sprzężeniem jednostkowym o G(s) = 4(s + 1)/(s(s + 4)) podano wejście r(t) = 5 + 2t (skok plus narastanie liniowe). Wyznacz uchyb ustalony ess.
G(s) = 4(s + 1)/(s(s + 4))r(t) = 5 + 2t
ess =
42.Uchyb ustalony i regulacja PIDTrudny
Regulator proporcjonalny K = 20 steruje obiektem G(s) = 10/(s(s + 8)) w pętli ze sprzężeniem jednostkowym. Znajdź 2-procentowy czas ustalania ts odpowiedzi skokowej, korzystając z ts ≈ 4/(ζωₙ).
G(s) = 10/(s(s + 8))K = 20
ts = s
3Roztwory
1.Reakcje podpór belek
Odpowiedź: RA = 34,38 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.
ΣMA = 0: RB·4 = 40·1,5 + 25·2,5
RB = 122,5/4 = 30,63 kN
ΣFy = 0: RA = 65 − 30,63 = 34,38 kN
Odpowiedź
RA = 34,38 kN
2.Reakcje podpór belek
Odpowiedź: RD = 21,17 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. Przegub wewnętrzny nie może przenieść momentu zginającego, co daje jedno dodatkowe równanie: moment względem przegubu sił po jednej stronie jest równy zeru.
Część CD
W przegubie nie ma momentu, więc momenty względem C dla samej części CD muszą się równoważyć.
W = 3·3 = 9 kN w 1,5 m od C
ΣMC = 0: RD·3 = 25·2 + 9·1,5
RD = 21,17 kN
VC = 34 − 21,17 = 12,83 kN
Odpowiedź
RD = 21,17 kN
3.Naprężenie i odkształcenie osiowe
Odpowiedź: δ = 0,06463 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 = π(80² − 50²)/4 = 3063 mm²
Łączna sztywność osiowa
ΣEA = 200·10³·1257 + 120·10³·3063 = 6,189 × 10⁸ N
Skrócenie
δ = 50000·800/6,189 × 10⁸ = 0,06463 mm
Odpowiedź
δ = 0,06463 mm
4.Naprężenie i odkształcenie osiowe
Odpowiedź: σmax = 63,66 MPa
Pręty obciążone osiowo: każdy pryzmatyczny odcinek spełnia δ = PL/(AE). Dla odcinków połączonych szeregowo wydłużenia się sumują.
Ta sama siła w obu odcinkach
Swobodny koniec przenosi P, więc każdy przekrój pręta przenosi tę samą siłę osiową P. Zmieniają się tylko pola przekroju i materiały.
A1 = π·40²/4 = 1257 mm², A2 = π·20²/4 = 314,2 mm²
σ1 = 15,92 MPa
σ2 = 63,66 MPa
Największe naprężenie
Cieńszy odcinek ma mniejsze pole, więc przenosi większe naprężenie.
σmax = 63,66 MPa
Odpowiedź
σmax = 63,66 MPa
5.Skręcanie wałów
Odpowiedź: φC = 18,28 °
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·ω.
Biegunowe momenty bezwładności
J1 = π·40⁴/32 = 251327 mm⁴
J2 = π·25⁴/32 = 38350 mm⁴
Kąt skręcenia każdego odcinka
T = 1500 N·m = 1,5 × 10⁶ N·mm
G = 80 GPa = 80000 MPa
φAB = 1,5 × 10⁶·1000/(80000·251327) = 0,0746 rad
φBC = 1,5 × 10⁶·500/(80000·38350) = 0,2445 rad
Całkowity kąt skręcenia
φC = 0,0746 + 0,2445 = 0,3191 rad = 18,28°
Odpowiedź
φC = 18,28°
6.Skręcanie wałów
Odpowiedź: d = 59,7 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π·300/60 = 31,42 rad/s
T = P/ω = 20000/31,42 = 636,6 N·m
Średnica z warunku wytrzymałości
T = 636,6 N·m = 636620 N·mm
dτ = ∛(16·636620/(π·50)) = 40,18 mm
Średnica z warunku sztywności
θ = 0,75·π/180/1000 = 1,309 × 10⁻⁵ rad/mm
G = 39 GPa = 39000 MPa
dφ = (32·636620/(π·39000·1,309 × 10⁻⁵))1/4 = 59,7 mm
Średnica decydująca
Oba warunki muszą być spełnione, więc decyduje większa średnica.
d = max(40,18, 59,7) = 59,7 mm
Odpowiedź
d = 59,7 mm
Tutaj decyduje sztywność.
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 = 20 kN, R2 = 50 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 = 20, MC = 40, MD = 60, ME = 80 kN·m, MF = 50, MG = 0
Odpowiedź
Punkt E: ME = 80 kN·m
8.Siła tnąca i moment zginający
Odpowiedź: F (x = 10 m)
Największy moment zginający: nachylenie wykresu M zmienia się tam, gdzie zmienia się siła tnąca, więc jego maksima leżą tam, gdzie siła tnąca przechodzi przez zero: pod siłą skupioną, wewnątrz obciążenia ciągłego albo nad podporą przy wsporniku. W utwierdzeniu belki wspornikowej moment jest największy przy ścianie. Porównaj wartości M w tych miejscach, zarówno dodatnie (zginanie „w dół”), jak i ujemne.
Reakcje podpór
Zapisz momenty względem jednej podpory, a potem użyj bilansu sił pionowych.
R1 = 16 kN, R2 = 29 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 = 32, MC = 14, MD = −4, ME = −22, MF = −40 kN·m, MG = 0
Odpowiedź
Punkt F: MF = −40 kN·m
9.Koło Mohra i transformacja naprężeń
Odpowiedź: θp = 5,655 °
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.
Kąt na kole
Na kole Mohra punkt x leży pod kątem 2θp od osi σ; ćwiartka wynika ze znaków σₓ − σᵧ i τₓᵧ.
2θp = atan2(2·10, 100 − 0) = 11,31°
Kąt na elemencie
θp = 5,655°
Odpowiedź
θp = 5,655°
Kąty dodatnie mierzy się przeciwnie do ruchu wskazówek zegara.
10.Koło Mohra i transformacja naprężeń
Odpowiedź: τθ = 5,195 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
τθ = −10·0,766 + 20·0,6428 = 5,195 MPa
Odpowiedź
τθ = 5,195 MPa
11.Ugięcie belki
Odpowiedź: θA = 0,01067 rad
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). Kąt obrotu θ = dv/dx to pierwsza całka z M/(EI); w spójnych jednostkach wychodzi w radianach.
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 = 12 N/mm, L = 8000 mm
EI = 2,4 × 10¹³ N·mm²
Kąt obrotu końca
θA = wL³/(24EI) = 12·8000³/(24·2,4 × 10¹³) = 0,01067 rad
Odpowiedź
θA = 0,01067 rad
To 0,6112°.
12.Ugięcie belki
Odpowiedź: hreq = 195,4 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⁴.
bh³/12 = Ireq → h = (12·7,766 × 10⁷/125)1/3 = 195,4 mm
Odpowiedź
hreq = 195,4 mm
W praktyce wybierz najbliższą większą wysokość znormalizowaną.
13.Wyboczenie słupów
Odpowiedź: Pcr = 15,5 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 = π·80⁴/64 = 2,011 × 10⁶ mm⁴
Promień bezwładności i smukłość
A = 5027 mm²
r = √(I/A) = √(2,011 × 10⁶/5027) = 20 mm
KL = 2·4000 = 8000 mm
KL/r = 8000/20 = 400
Naprężenie krytyczne
σcr = π²E/(KL/r)² = π²·200000/400² = 12,34 MPa
Odpowiedź
σcr = 12,34 MPa
15.Rzut ukośny
Odpowiedź: H = 4,995 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 θ = 14·cos 45° = 9,899 m/s
vy = v0·sin θ = 14·sin 45° = 9,899 m/s
Najwyższy punkt
Na szczycie prędkość pionowa jest zerowa: vy² = 2gH.
H = vy²/(2g) = 9,899²/(2·9,81) = 4,995 m
Odpowiedź
H = 4,995 m
16.Rzut ukośny
Odpowiedź: α = 59,85 °
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 = 20·cos 15° = 19,32 m/s
vy0 = 20·sin 15° = 5,176 m/s
Prędkość pionowa przy uderzeniu
Energia albo v² = v₀² + 2as w kierunku pionowym.
vy = √(5,176² + 2·9,81·55) = 33,25 m/s
Kąt uderzenia
tan α = vy/vx = 33,25/19,32 → α = 59,85°
Odpowiedź
α = 59,85°
17.Praca i energia
Odpowiedź: Wf = 38,23 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 θ = 49,05·0,866 = 42,48 N
F = μN = 0,3·42,48 = 12,74 N
Praca przeciw tarciu
Wf = F·d = 12,74·3 = 38,23 J
Odpowiedź
Wf = 38,23 J
18.Praca i energia
Odpowiedź: x = 63,61 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 = 49,05·(4,5 − 0,15·2,5) = 202,3 J
Maksymalne ugięcie
½kx² = 202,3 → x = √(2·202,3/1000) = 0,6361 m = 63,61 cm
Odpowiedź
x = 63,61 cm
19.Drgania mechaniczne
Odpowiedź: ζ = 0,2
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·√(20000·200) = 4000 N·s/m
Współczynnik tłumienia
ζ = c/cc = 800/4000 = 0,2
Odpowiedź
ζ = 0,2
ζ < 1, więc tłumienie jest podkrytyczne, a układ drga z malejącą amplitudą.
20.Drgania mechaniczne
Odpowiedź: keq = 32 kN/m
Drgania swobodne układu masa–sprężyna: pulsacja drgań własnych nietłumionych to ωn = √(k/m) w rad/s, a fn = ωn/(2π) w Hz. Tłumik wiskotyczny c daje współczynnik tłumienia ζ = c/cc, przy tłumieniu krytycznym cc = 2·√(k·m) = 2·m·ωn. Układ z tłumieniem podkrytycznym (ζ < 1) drga z pulsacją ωd = ωn·√(1 − ζ²). Sztywności sprężyn połączonych równolegle się dodaje, k = k₁ + k₂; dla sprężyn szeregowych dodaje się odwrotności, 1/k = 1/k₁ + 1/k₂.
Sztywność zastępcza
Sprężyny obok siebie mają to samo ugięcie, a ich siły się sumują, więc ich sztywności po prostu się dodaje: k = k1 + k2.
keq = 12 + 20 = 32 kN/m
Odpowiedź
keq = 32 kN/m
21.Silniki cieplne i chłodziarki
Odpowiedź: Q̇H = 9,609 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 = 30 + 273,15 = 303,1 K
TL = -5 + 273,15 = 268,1 K
Stosunek ciepeł w obiegu odwracalnym
Q̇H = Q̇L·TH/TL = 8,5·303,1/268,1 = 9,609 kW
Odpowiedź
Q̇H = 9,609 kW
Moc pobierana: Ẇ = 1,109 kW.
22.Silniki cieplne i chłodziarki
Odpowiedź: Ẇmin = 0,3582 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̇H/COP = 7/19,54 = 0,3582 kW
Odpowiedź
Ẇmin = 0,3582 kW
23.Obiegi silników cieplnych gazowych
Odpowiedź: T3 = 1982 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·170,4 = 900,7 K
Doprowadzenie ciepła przy stałym ciśnieniu
T3 = rc·T2 = 2,2·900,7 = 1982 K
Odpowiedź
T3 = 1982 K
24.Obiegi silników cieplnych gazowych
Odpowiedź: wnet = 419,9 kJ/kg
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.
Temperatury izentropowe
T2 = 288·150,2857 = 624,3 K
T4 = 1400/150,2857 = 645,8 K
Praca sprężarki i turbiny
wc = 1,005·(624,3 − 288) = 338 kJ/kg
wt = 1,005·(1400 − 645,8) = 758 kJ/kg
Praca netto
wnet = 758 − 338 = 419,9 kJ/kg
Odpowiedź
wnet = 419,9 kJ/kg
Wskaźnik pracy zwrotnej: wc/wt = 44,6%.
25.Ciągłość i Bernoulli
Odpowiedź: Q = 5,479 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·1) = 4,429 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.
Prędkość średnia
A = π·0,2²/4 = 0,03142 m²
V = Q/A = 0,025/0,03142 = 0,7958 m/s
Wysokość prędkości
V²/(2g) = 0,7958²/(2·9,81) = 0,03228 m
Darcy–Weisbach
hf = f·(L/D)·V²/(2g) = 0,02·(200/0,2)·0,7958²/(2·9,81) = 0,6455 m
Odpowiedź
hf = 0,6455 m
28.Tarcie w rurach i straty wysokości
Odpowiedź: Δz = 2,264 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,04/0,03142 = 1,273 m/s
V²/(2g) = 0,08263 m
Współczynniki strat
f·L/D = 0,015·200/0,2 = 15
ΣK = 0,5 + 3·0,3 + 10 + 1,0 = 12,4
Całkowita strata wysokości
hf = 15·0,08263 = 1,239 m
hm = 12,4·0,08263 = 1,025 m
Δz = hf + hm = 2,264 m
Odpowiedź
Δz = 2,264 m
29.Pompy i moc pomp
Odpowiedź: Pshaft = 49,05 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,05·60 = 29,43 kW
Moc na wale
Pshaft = Ph/η = 29,43/0,6 = 49,05 kW
Odpowiedź
Pshaft = 49,05 kW
30.Pompy i moc pomp
Odpowiedź: H = 24,53 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 = 13 + 11,53 = 24,53 m
Odpowiedź
H = 24,53 m
31.Przewodzenie ciepła
Odpowiedź: Ti = −0,4634 °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.
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.
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.
Moment na zębniku
ω = 2π·1450/60 = 151,8 rad/s
T = P/ω = 2200/151,8 = 14,49 N·m
Średnica podziałowa
d = m·N = 5·25 = 125 mm
Siła obwodowa
Ft = 2T/d = 2·14,49/0,125 = 231,8 N
Odpowiedź
Ft = 231,8 N
34.Koła i przekładnie zębate
Odpowiedź: nout = 460,8 rpm
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.
Pierwszy stopień
nB = 1440·15/25 = 864 rpm = nC
Drugi stopień
nout = 864·16/30 = 460,8 rpm
Odpowiedź
nout = 460,8 rpm
przełożenie całkowite = 3,125 : 1
35.Transformaty Laplace’a
Odpowiedź: f(0.4) = −0,323
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² + 8s + 16 = (s + 4)²
Ułamki proste
6s + 5 = 6(s + 4) + −19
F(s) = 6/(s + 4) − 19/(s + 4)²
Transformata odwrotna
f(t) = 6·e−4t − 19·t·e−4t
f(0,4) = (6 − 19·0,4)·e−4·0,4 = −0,323
Odpowiedź
f(0,4) = −0,323
36.Transformaty Laplace’a
Odpowiedź: g(0.8) = −1,292
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).
Transmitancja
G(s) = Y(s)/U(s) = 19/(s² + 2s + 26)
Ułamki proste
s² + 2s + 26 = (s + 1)² + 5²
G(s) = (19/5)·5/((s + 1)² + 5²)
Odpowiedź impulsowa
g(t) = e−t·(3,8·sin 5t)
g(0,8) = −1,292
Odpowiedź
g(0,8) = −1,292
37.Schematy blokowe
Odpowiedź: ess = 0,25
Algebra schematów blokowych: bloki połączone szeregowo mnożą się (G₁·G₂), bloki połączone równolegle dodają się (G₁ + G₂), a pętla ze ścieżką prostą G i ścieżką sprzężenia zwrotnego H redukuje się do T = G/(1 + G·H) dla ujemnego sprzężenia zwrotnego i G/(1 − G·H) dla dodatniego. Najpierw redukuj pętle wewnętrzne, potem idź na zewnątrz. Bieguny układu zamkniętego są pierwiastkami równania charakterystycznego 1 + G(s)H(s) = 0. Dla stabilnej pętli wymuszenie skokowe o wielkości A daje w stanie ustalonym wyjście A·T(0); przy jednostkowym sprzężeniu zwrotnym uchyb ustalony przy skoku wynosi A/(1 + Kp), gdzie Kp = G(0).
Stała uchybu położeniowego
Kₚ = G(0) = 12·3/(2·6) = 3
Uchyb ustalony
ess = A/(1 + Kₚ) = 1/(1 + 3) = 0,25
Odpowiedź
ess = 0,25
38.Schematy blokowe
Odpowiedź: Kmax = 126
Algebra schematów blokowych: bloki połączone szeregowo mnożą się (G₁·G₂), bloki połączone równolegle dodają się (G₁ + G₂), a pętla ze ścieżką prostą G i ścieżką sprzężenia zwrotnego H redukuje się do T = G/(1 + G·H) dla ujemnego sprzężenia zwrotnego i G/(1 − G·H) dla dodatniego. Najpierw redukuj pętle wewnętrzne, potem idź na zewnątrz. Bieguny układu zamkniętego są pierwiastkami równania charakterystycznego 1 + G(s)H(s) = 0. Dla stabilnej pętli wymuszenie skokowe o wielkości A daje w stanie ustalonym wyjście A·T(0); przy jednostkowym sprzężeniu zwrotnym uchyb ustalony przy skoku wynosi A/(1 + Kp), gdzie Kp = G(0).
Równanie charakterystyczne
s³ + 9s² + 14s + K = 0
Tablica Routha
s³ | 1 14
s² | 9 K
s¹ | (9·14 − K)/9
s⁰ | K
Granica stabilności
Wszystkie elementy pierwszej kolumny muszą być dodatnie.
Kmax = 9·14 = 126
Odpowiedź
Kmax = 126
39.Odpowiedź skokowa
Odpowiedź: ts = 4 s
Odpowiedź skokowa: układ pierwszego rzędu G(s) = K/(τs + 1) odpowiada na skok o wielkości U przebiegiem y(t) = K·U·(1 − e−t/τ): po jednej stałej czasowej τ osiąga 63,2 % wartości końcowej K·U, a po 4τ 98 %. Układ drugiego rzędu G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) ma bieguny s = −ζωₙ ± jωd, gdzie ωd = ωₙ√(1 − ζ²). Dla 0 < ζ < 1 odpowiedź skokowa ma przeregulowanie: przeregulowanie procentowe %OS = 100·e−ζπ/√(1−ζ²), czas szczytu tp = π/ωd, a 2-procentowy czas ustalania ts ≈ 4/(ζωₙ). Porównaj mianownik z s² + 2ζωₙs + ωₙ², aby odczytać ωₙ i ζ. Dla układu masa–sprężyna–tłumik ωₙ = √(k/m) i ζ = c/(2√(km)); dla szeregowego obwodu RLC ωₙ = 1/√(LC) i ζ = (R/2)·√(C/L).
Porównaj z postacią standardową
Mianownik musi odpowiadać s² + 2ζωₙs + ωₙ².
ωₙ² = 54 → ωₙ = 7,348 rad/s
2ζωₙ = 2 → ζ = 2/(2·7,348) = 0,1361
Wynik
ζωₙ = 1 → ts ≈ 4/1 = 4 s
Odpowiedź
ts = 4 s
40.Odpowiedź skokowa
Odpowiedź: K = 93,56
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).
Transmitancja układu zamkniętego
T(s) = K/(s² + 10s + K)
Współczynnik tłumienia z przeregulowania
ζ = −ln(0,15)/√(π² + ln²(0,15)) = 0,5169
Wzmocnienie
ωₙ = 10/(2ζ) = 9,672 rad/s
K = ωₙ² = 93,56
Odpowiedź
K = 93,56
41.Uchyb ustalony i regulacja PID
Odpowiedź: ess = 2
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 prędkościowa
Kv = lim s·G(s) = 1
Część skokowa
Kp = ∞ → estep = 0
Uchyb dla narastania liniowego
ess = 2/1 = 2
Odpowiedź
ess = 2
42.Uchyb ustalony i regulacja PID
Odpowiedź: ts = 1 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) + 20·10 = s² + 8s + 200
ωₙ i ζ
ωₙ = √200 = 14,14 rad/s
ζ = 8/(2·14,14) = 0,2828
Czas ustalania
ζωₙ = 8/2 = 4
ts ≈ 4/4 = 1 s
Odpowiedź
ts = 1 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.