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ę nieprzesuwną w A i przesuwną w B. Przenosi siłę skupioną P w odległości 4,5 m od A i moment skupiony M działający zgodnie z ruchem wskazówek zegara, przyłożony w odległości 0,5 m od A. Wyznacz reakcję pionową RB, przyjmując zwrot w górę jako dodatni.
L = 5 mP = 60 kN w 4,5 mM = 10 kN·m (zgodnie z ruchem wskazówek zegara) w 0,5 m
RB = kN
2.Reakcje podpór belekTrudny
Belka złożona jest utwierdzona w ścianie w A, ma przegub wewnętrzny w C (2 m od A) i opiera się na podporze przesuwnej w D (6 m od A). Przenosi siłę skupioną P w odległości 4,5 m od A. Wyznacz reakcję RD na podporze przesuwnej.
AC = 2 mAD = 6 mP = 65 kN w 4,5 m
RD = kN
3.Naprężenie i odkształcenie osioweŚredni
Wyznacz wydłużenie δ pręta (stal) pod obciążeniem osiowym.
P = 55 kNd = 12 mmL = 2,5 mE = 200 GPa
δ = mm
4.Naprężenie i odkształcenie osioweTrudny
Pręt (materiał: stal, średnica 30 mm, E = 200 GPa, naprężenie dopuszczalne 250 MPa) znajduje się wewnątrz rury (materiał: mosiądz, średnica zewnętrzna 80 mm, średnica wewnętrzna 50 mm, E = 100 GPa, naprężenie dopuszczalne 50 MPa). Są razem ściskane między sztywnymi płytami. Wyznacz największe obciążenie Pmax.
pręt: Ø30 mm, E = 200 GPa, σallow = 250 MParura: Ø80/Ø50 mm, E = 100 GPa, σallow = 50 MPa
Pmax = kN
5.Skręcanie wałówŚredni
Pełny wał schodkowy (stal) jest utwierdzony w A. Odcinek AB ma średnicę 40 mm i długość 2 m; odcinek BC ma średnicę 30 mm i długość 3 m. Na swobodnym końcu C działa moment skręcający 500 N·m. Wyznacz maksymalne naprężenie styczne τmax w wale.
T = 500 N·md1 = 40 mm, L1 = 2 md2 = 30 mm, L2 = 3 m
τmax = MPa
6.Skręcanie wałówTrudny
Wał drążony przenosi 15 kW przy 1800 obr/min. Wyznacz maksymalne naprężenie styczne τmax w wale.
P = 15 kWN = 1800 rpmdo = 40 mmdi = 20 mm
τmax = MPa
7.Siła tnąca i moment zginającyŚredni
Wyznacz moment zginający Mx w przekroju x = 6,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 = 8 mw = 10 kN/m od 0 do 8 mx = 6,5 m
Mx = kN·m
8.Siła tnąca i moment zginającyTrudny
Belka opiera się na podporze nieprzesuwnej w A i przesuwnej w B i wystaje poza B o 1 m aż do swobodnego końca C. Na całej długości A–C działa obciążenie równomiernie rozłożone w. Wyznacz maksymalny dodatni moment zginający Mmax między podporami.
L = 7 mw = 2 kN/m od 0 do 7 mpodpory w x = 0 i x = 6 m
Mmax = kN·m
9.Koło Mohra i transformacja naprężeńŚredni
W pewnym punkcie obciążonej części występują σₓ = 0 MPa, σᵧ = −30 MPa i τₓᵧ = 10 MPa (rozciąganie dodatnie). Wyznacz maksymalne naprężenie styczne w płaszczyźnie τmax.
σx = 0 MPaσy = −30 MPaτxy = 10 MPa
τmax = MPa
10.Koło Mohra i transformacja naprężeńTrudny
Cienkościenny zbiornik cylindryczny o średnicy wewnętrznej 400 mm i grubości ścianki 6 mm zawiera gaz pod nadciśnieniem 3 MPa. Wykonano go z blachy spawanej wzdłuż spoiny śrubowej nachylonej pod kątem 50° do osi walca. Wyznacz naprężenie normalne σn działające w poprzek spoiny.
D = 400 mmt = 6 mmp = 3 MPaα = 50°
σn = MPa
11.Ugięcie belkiŚredni
Wspornik o długości L jest utwierdzony w A i przenosi siłę skupioną P w odległości a = 2 m od A. Wyznacz ugięcie δB swobodnego końca B.
P = 4 kNa = 2 mL = 5 mE = 200 GPa (stal)I = 57 × 10⁶ mm⁴
δB = mm
12.Ugięcie belkiTrudny
Belka drewniana swobodnie podparta o przekroju prostokątnym i szerokości 75 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ść/250. Wyznacz najmniejszą wysokość hreq przekroju.
w = 1 kN/mL = 6 mE = 11 GPa (drewno konstrukcyjne)b = 75 mm
hreq = mm
13.Wyboczenie słupówŚredni
Słup o długości L z materiału (aluminium) jest przegubowo podparty na obu końcach. Korzystając ze wzoru Eulera, wyznacz jego siłę krytyczną Pcr.
L = 3 mE = 70 GPa (aluminium)d = 40 mm
Pcr = kN
14.Wyboczenie słupówTrudny
Pełny słup okrągły o długości L z materiału (stal) jest przegubowo podparty na obu końcach. Ma przenieść 30 kN przy współczynniku bezpieczeństwa 3 względem wyboczenia. Wyznacz najmniejszą średnicę d.
P = 30 kNn = 3L = 5 mE = 200 GPa (stal)
d = mm
15.Rzut ukośnyŚredni
Piłka zostaje wyrzucona pod kątem 65° do poziomu z płaskiego terenu. Pomiń opór powietrza. Wyznacz odległość poziomą x pokonaną, zanim piłka uderzy w ziemię.
v0 = 30 m/sθ = 65°g = 9,81 m/s²
x = m
16.Rzut ukośnyTrudny
Piłkę rzucono z poziomu ziemi z prędkością 26 m/s pod kątem 55° do poziomu. Ląduje, już opadając, na płaskim dachu budynku o wysokości 15 m. Pomiń opór powietrza i przyjmij g = 9,81 m/s². Wyznacz odległość poziomą x od punktu wyrzutu do miejsca lądowania.
v₀ = 26 m/sθ = 55°hroof = 15 mg = 9,81 m/s²
x = m
17.Praca i energiaŚredni
Samochód jadący z prędkością 40 km/h po poziomej drodze gwałtownie hamuje i z zablokowanymi kołami sunie aż do zatrzymania. Współczynnik tarcia kinetycznego między oponami a drogą wynosi 0,6. Wyznacz drogę hamowania d.
v = 40 km/hμ = 0,6g = 9,81 m/s²
d = m
18.Praca i energiaTrudny
Klocek o masie 20 kg rusza ze spoczynku u podstawy zbocza nachylonego pod kątem 15° i jest ciągnięty w górę zbocza stałą siłą 115 N równoległą do zbocza. Współczynnik tarcia kinetycznego wynosi 0,1. Wyznacz prędkość v klocka po przebyciu 4 m w górę zbocza.
Masa 20 kg jest podparta na sprężynie o sztywności 100 kN/m i tłumiku wiskotycznym o współczynniku tłumienia 570 N·s/m. Wyznacz współczynnik tłumienia względnego ζ.
m = 20 kgk = 100 kN/mc = 570 N·s/m
ζ =
20.Drgania mechaniczneTrudny
Masa 40 kg wisi na dwóch sprężynach obok siebie (równolegle) o sztywnościach k₁ = 25 kN/m i k₂ = 10 kN/m. Masa pozostaje pozioma. Wyznacz sztywność zastępczą keq układu sprężyn w kN/m.
m = 40 kgk1 = 25 kN/mk2 = 10 kN/m
keq = kN/m
21.Silniki cieplne i chłodziarkiŚredni
Silnik cieplny pracuje między zbiornikiem gorącym o temperaturze 425 °C a zbiornikiem zimnym o temperaturze 20 °C. Wyznacz największą możliwą sprawność cieplną ηmax.
TH = 425 °CTL = 20 °C
ηmax = %
22.Silniki cieplne i chłodziarkiTrudny
Pompa ciepła utrzymuje w domu 20 °C, gdy na zewnątrz jest -10 °C, dostarczając do domu 13 kW ciepła. Wyznacz minimalną moc napędową Ẇmin.
TH = 20 °CTL = -10 °CQ̇H = 13 kW
Ẇmin = kW
23.Obiegi silników cieplnych gazowychŚredni
Idealny obieg Diesla ma stopień sprężania r = 18 i stopień obciążenia rc = v₃/v₂ = 1,5. Na początku sprężania powietrze ma temperaturę 300 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 = 18rc = 1,5T1 = 300 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 6. Powietrze wpływa do sprężarki w temperaturze 288 K, a do turbiny w temperaturze 1100 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 = 6T1 = 288 KT3 = 1100 Kk = 1,4cp = 1,005 kJ/(kg·K)
wnet = kJ/kg
25.Ciągłość i BernoulliŚredni
Woda płynie przez poziomą rurę, która zwęża się z 75 mm do 25 mm. W szerokim odcinku ciśnienie wynosi 110 kPa, a prędkość 1 m/s. Pomijając straty, wyznacz ciśnienie p2 w wąskim odcinku.
Pozioma zwężka Venturiego (wlot 250 mm, gardziel 125 mm) przepuszcza wodę. Rtęciowy manometr U-rurkowy podłączony między wlotem a gardzielą pokazuje różnicę poziomów 140 mm. Współczynnik przepływu wynosi 0,97. Wyznacz natężenie przepływu Q. Przyjmij ρ = 1000 kg/m³ dla wody, 13 600 kg/m³ dla rtęci i g = 9,81 m/s².
Woda płynie z natężeniem 5 L/s przez rurę o długości 300 m i średnicy wewnętrznej 50 mm. Przyjmij współczynnik tarcia Darcy’ego f = 0,03 i g = 9,81 m/s². Wyznacz stratę wysokości na tarcie hf.
Q = 5 L/sL = 300 mD = 50 mmf = 0,03g = 9,81 m/s²
hf = m
28.Tarcie w rurach i straty wysokościTrudny
Dwa otwarte zbiorniki są połączone rurą o długości 150 m, średnicy wewnętrznej 250 mm i współczynniku tarcia Darcy’ego f = 0,018. Zwierciadło wody w górnym zbiorniku leży 16 m wyżej niż w dolnym. Straty miejscowe to: wlot o ostrej krawędzi (K = 0,5), kolana w liczbie 2 (K = 0,3 każde), zawór (K = 0,2) i wylot do dolnego zbiornika (K = 1,0). Przyjmij g = 9,81 m/s². Wyznacz natężenie przepływu Q w L/s.
Pompa o sprawności 60% tłoczy 30 L/s wody (ρ = 1000 kg/m³) przy wysokości podnoszenia 54 m. Przyjmij g = 9,81 m/s². Wyznacz potrzebną moc na wale Pshaft, w kW.
Q = 30 L/sH = 54 mη = 60%g = 9,81 m/s²
Pshaft = kW
30.Pompy i moc pompTrudny
Pompa przetłacza 80 L/s wody (ρ = 1000 kg/m³) z otwartego dolnego zbiornika do górnego zbiornika, którego zwierciadło wody leży 7 m wyżej. Rura ma długość 200 m, średnicę wewnętrzną 200 mm i współczynnik tarcia Darcy’ego f = 0,015; współczynniki strat miejscowych sumują się do ΣK = 6 (łącznie z wlotem i wylotem). Przyjmij g = 9,81 m/s². Górny zbiornik jest zamknięty, a powietrze nad wodą ma nadciśnienie 350 kPa. Wyznacz wysokość podnoszenia H, którą musi zapewnić pompa.
Ściana składa się z 250 mm (drewno) i 150 mm (styropian), a jej pole wynosi 8 m². Powierzchnie zewnętrzne mają 19 °C i -11 °C. Wyznacz strumień ciepła Q̇ przez ścianę.
L1 = 250 mm, k1 = 0,17 W/(m·K)L2 = 150 mm, k2 = 0,035 W/(m·K)A = 8 m²T1 = 19 °CT3 = -11 °C
Q̇ = W
32.Przewodzenie ciepłaTrudny
Ściana z 200 mm (drewno) i 100 mm (szkło) oddziela pomieszczenie o temperaturze 20 °C od powietrza zewnętrznego o temperaturze -8 °C. Współczynniki przejmowania ciepła wynoszą 7,7 W/(m²·K) wewnątrz i 25 W/(m²·K) na zewnątrz, a pole ściany 15 m². Wyznacz straty ciepła Q̇ przez ścianę.
L1 = 200 mm, k1 = 0,17 W/(m·K)L2 = 100 mm, k2 = 0,96 W/(m·K)A = 15 m²hi = 7,7 W/(m²·K)ho = 25 W/(m²·K)Troom = 20 °CTout = -8 °C
Q̇ = W
33.Koła i przekładnie zębateŚredni
Silnik o mocy 1,5 kW pracujący z prędkością 1500 rpm napędza zębnik o 16 zębach, który zazębia się z kołem o 80 zębach. Sprawność pary kół wynosi 96%. Wyznacz moment wyjściowy T2 na wale koła.
P = 1,5 kWn1 = 1500 rpmN1 = 16N2 = 80η = 96%
T2 = N·m
34.Koła i przekładnie zębateTrudny
Silnik o mocy 15 kW i prędkości 2900 rpm napędza koło A (20 zębów), które współpracuje z kołem B (40 zębów). Koło C (16 zębów) na tym samym wale co B napędza koło D (25 zębów) na wale wyjściowym. Każda para kół ma sprawność 96 %. Wyznacz moment wyjściowy Tout.
Znajdź f(t) dla F(s) = 35/(s(s + 3)(s + 4)) metodą rozkładu na ułamki proste i oblicz f(0,1).
F(s) = 35/(s(s + 3)(s + 4))t = 0,1 s
f(0.1) =
36.Transformaty Laplace’aTrudny
Układ opisuje równanie różniczkowe y″ + 4y′ + 13y = u′ + 11u, 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,5 s. Znajdź g(0,5).
y″ + 4y′ + 13y = u′ + 11ut = 0,5 s
g(0.5) =
37.Schematy blokoweŚredni
Obiekt G(s) = K/(s(s + 3)) jest w pętli jednostkowego ujemnego sprzężenia zwrotnego. Znajdź wzmocnienie K, dla którego układ zamknięty jest tłumiony krytycznie (dwa równe bieguny rzeczywiste).
Regulator proporcjonalny K steruje obiektem G(s) = 10/(s + 10) w pętli jednostkowego ujemnego sprzężenia zwrotnego. Skok zakłócenia o wielkości 5 wchodzi na wejście obiektu, gdy R = 0. Znajdź najmniejsze wzmocnienie K, przy którym wyjście w stanie ustalonym nie przekracza 0,1.
G(s) = 10/(s + 10)D = 5css ≤ 0,1
K =
39.Odpowiedź skokowaŚredni
Układ ma transmitancję G(s) = 3,5/(s² + 3s + 7) i otrzymuje jednostkowe wymuszenie skokowe. Znajdź przeregulowanie procentowe %OS odpowiedzi skokowej.
G(s) = 3,5/(s² + 3s + 7)
%OS = %
40.Odpowiedź skokowaTrudny
Pętla jednostkowego ujemnego sprzężenia zwrotnego ma G(s) = 95/(s(s + 2)). Znajdź 2-procentowy czas ustalania ts odpowiedzi skokowej układu zamkniętego, korzystając z ts ≈ 4/(ζωₙ).
G(s) = 95/(s(s + 2))
ts = s
41.Uchyb ustalony i regulacja PIDŚredni
Na pętlę ze sprzężeniem jednostkowym o G(s) = 20(s + 2)/(s(s + 10)) podano wejście r(t) = 10 + 1t (skok plus narastanie liniowe). Wyznacz uchyb ustalony ess.
G(s) = 20(s + 2)/(s(s + 10))r(t) = 10 + 1t
ess =
42.Uchyb ustalony i regulacja PIDTrudny
Regulator PI C(s) = 1 + 1/s steruje obiektem G(s) = 2/((s + 2)(s + 5)) w pętli ze sprzężeniem jednostkowym. Wyznacz uchyb ustalony ess dla wejścia liniowo narastającego r(t) = 5t.
C(s) = 1 + 1/sG(s) = 2/((s + 2)(s + 5))r(t) = 5t
ess =
3Roztwory
1.Reakcje podpór belek
Odpowiedź: RB = 56 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·5 − 60·4,5 − 10 = 0
RB = 280/5 = 56 kN
Równowaga sił pionowych
RA = 60 − 56 = 4 kN
Odpowiedź
RB = 56 kN
Ujemna reakcja działa w dół.
2.Reakcje podpór belek
Odpowiedź: RD = 40,63 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ć.
ΣMC = 0: RD·4 = 65·2,5
RD = 40,63 kN
VC = 65 − 40,63 = 24,38 kN
Odpowiedź
RD = 40,63 kN
3.Naprężenie i odkształcenie osiowe
Odpowiedź: δ = 6,079 mm
Prawo Hooke’a: w zakresie sprężystym σ = E·ε. Łącząc to z σ = P/A i ε = δ/L, dostajemy δ = PL/(AE).
Naprężenie
A = π·12²/4 = 113,1 mm²
σ = P/A = 55·1000/113,1 = 486,3 MPa
Odkształcenie (prawo Hooke’a)
E for stal = 200 GPa = 200000 MPa.
ε = σ/E = 486,3/200000 = 0,002432 = 2432 µε
Wydłużenie
δ = ε·L = 0,002432·2500 mm = 6,079 mm
Odpowiedź
δ = 6,079 mm
4.Naprężenie i odkształcenie osiowe
Odpowiedź: Pmax = 223,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.
Pola
Ar = π·30²/4 = 706,9 mm²
At = π(80² − 50²)/4 = 3063 mm²
Łączna sztywność osiowa
ΣEA = 200·10³·706,9 + 100·10³·3063 = 4,477 × 10⁸ N
Decydujące odkształcenie
ε = min(0,00125, 0,0005) = 0,0005 (mosiądz)
Największe obciążenie
Pmax = 0,0005·4,477 × 10⁸ = 223838 N = 223,8 kN
Odpowiedź
Pmax = 223,8 kN
5.Skręcanie wałów
Odpowiedź: τmax = 94,31 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·ω.
Naprężenie w każdym odcinku
T = 500 N·m = 500000 N·mm
τAB = 16·500000/(π·40³) = 39,79 MPa
τBC = 16·500000/(π·30³) = 94,31 MPa
Odpowiedź
τmax = 94,31 MPa
6.Skręcanie wałów
Odpowiedź: τmax = 6,755 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π·1800/60 = 188,5 rad/s
T = P/ω = 15000/188,5 = 79,58 N·m
Biegunowy moment bezwładności
J = π(40⁴ − 20⁴)/32 = 235619 mm⁴
Maksymalne naprężenie styczne
τmax = T·(do/2)/J = 79577·20/235619 = 6,755 MPa
Odpowiedź
τmax = 6,755 MPa
7.Siła tnąca i moment zginający
Odpowiedź: Mx = 48,75 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·8 = 80 kN, przyłożona w x̄ = 4 m
ΣMA = 0: RB·8 = 10·8·4
RB = 320/8 = 40 kN
ΣFy = 0: RA = 80 − 40 = 40 kN
Przetnij w x = 6,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)
Mx = 40·6,5 − 10·6,5·3,25 = 48,75 kN·m
Odpowiedź
Mx = 48,75 kN·m
8.Siła tnąca i moment zginający
Odpowiedź: Mmax = 8,507 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, łącznie z obciążeniem na wsporniku.
Wypadkowa obciążenia ciągłego: W = w·ℓ = 2·7 = 14 kN, przyłożona w x̄ = 3,5 m
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 = √(((0 − −30)/2)² + 10²) = 18,03 MPa
Odpowiedź
τmax = 18,03 MPa
10.Koło Mohra i transformacja naprężeń
Odpowiedź: σn = 70,66 MPa
Transformacja płaskiego stanu naprężenia: na płaszczyźnie obróconej o kąt θ: σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ. Koło Mohra pokazuje wszystkie płaszczyzny naraz: jego środek leży w C = (σₓ + σᵧ)/2, a promień wynosi R = √(((σₓ − σᵧ)/2)² + τₓᵧ²). Naprężenia główne to C ± R, a największe naprężenie styczne w płaszczyźnie to R. Cylinder cienkościenny: naprężenie obwodowe σh = pr/t i naprężenie osiowe σa = pr/(2t); są to naprężenia główne w ściance.
Naprężenia w ściance
r = 400/2 = 200 mm
σh = pr/t = 3·200/6 = 100 MPa
σa = pr/(2t) = 50 MPa
Transformacja na kierunek spoiny
Przyjmij oś x wzdłuż osi walca (σa), a oś y po obwodzie (σh). Normalna spoiny tworzy z osią kąt 90° − α.
σn = 50·sin²50° + 100·cos²50° = 70,66 MPa
Odpowiedź
σn = 70,66 MPa
11.Ugięcie belki
Odpowiedź: δB = 3,041 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 = 4000 N, a = 2000 mm, L = 5000 mm
EI = 1,14 × 10¹³ N·mm²
Ugięcie swobodnego końca
Ugięcie pod siłą plus kąt obrotu pod siłą razy pozostała długość.
δP = Pa³/(3EI) = 0,9357 mm
θP = Pa²/(2EI) = 0,0007018 rad
δB = δP + θP·(L − a) = Pa²(3L − a)/(6EI) = 3,041 mm
Odpowiedź
δB = 3,041 mm
12.Ugięcie belki
Odpowiedź: hreq = 217,1 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·6,392 × 10⁷/75)1/3 = 217,1 mm
Odpowiedź
hreq = 217,1 mm
W praktyce wybierz najbliższą większą wysokość znormalizowaną.
13.Wyboczenie słupów
Odpowiedź: Pcr = 9,646 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.
Wymagana siła krytyczna
Pcr = n·P = 3·30000 = 90000 N
KL = 1·5000 = 5000 mm
Wyznacz średnicę
d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·90000·5000²/(π³·200000))^(1/4) = 69,42 mm
Odpowiedź
d = 69,42 mm
Sprawdzenie: smukłość KL/r wynosi 288,1, więc słup jest smukły i wzór Eulera obowiązuje.
15.Rzut ukośny
Odpowiedź: x = 70,28 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 65° = 12,68 m/s
vy = v0·sin θ = 30·sin 65° = 27,19 m/s
Czas lotu
Start i lądowanie na tej samej wysokości: t = 2vy/g.
t = 2·27,19/9,81 = 5,543 s
Odległość pozioma
x = vx·t = 12,68·5,543 = 70,28 m
Odpowiedź
x = 70,28 m
16.Rzut ukośny
Odpowiedź: x = 51,56 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 = 14,91 m/s, vy = 21,3 m/s
Czas do lądowania na dachu
Większy pierwiastek: piłka opada.
15 = 21,3·t − 4,905·t²
t = (21,3 + √(21,3² − 2·9,81·15))/9,81 = 3,458 s
Odległość pozioma
x = vx·t = 14,91·3,458 = 51,56 m
Odpowiedź
x = 51,56 m
17.Praca i energia
Odpowiedź: d = 10,49 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 = 40/3,6 = 11,11 m/s
Bilans energii
½mv² = μmg·d → d = v²/(2μg) = 11,11²/(2·0,6·9,81) = 10,49 m
Odpowiedź
d = 10,49 m
Podwojenie prędkości wydłuża drogę hamowania czterokrotnie.
18.Praca i energia
Odpowiedź: v = 4,255 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).
Siły wzdłuż równi
m·g = 20·9,81 = 196,2 N
mg·sin θ = 50,78 N
μ·mg·cos θ = 18,95 N
Praca i energia
½mv² = (115 − 69,73)·4 = 181,1 J
v = √(2·181,1/20) = 4,255 m/s
Odpowiedź
v = 4,255 m/s
19.Drgania mechaniczne
Odpowiedź: ζ = 0,2015
Drgania swobodne układu masa–sprężyna: pulsacja drgań własnych nietłumionych to ωn = √(k/m) w rad/s, a fn = ωn/(2π) w Hz. Tłumik wiskotyczny c daje współczynnik tłumienia ζ = c/cc, przy tłumieniu krytycznym cc = 2·√(k·m) = 2·m·ωn. Układ z tłumieniem podkrytycznym (ζ < 1) drga z pulsacją ωd = ωn·√(1 − ζ²). Sztywności sprężyn połączonych równolegle się dodaje, k = k₁ + k₂; dla sprężyn szeregowych dodaje się odwrotności, 1/k = 1/k₁ + 1/k₂.
Tłumienie krytyczne
cc = 2·√(k·m) = 2·√(100000·20) = 2828 N·s/m
Współczynnik tłumienia
ζ = c/cc = 570/2828 = 0,2015
Odpowiedź
ζ = 0,2015
ζ < 1, więc tłumienie jest podkrytyczne, a układ drga z malejącą amplitudą.
20.Drgania mechaniczne
Odpowiedź: keq = 35 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 = 25 + 10 = 35 kN/m
Odpowiedź
keq = 35 kN/m
21.Silniki cieplne i chłodziarki
Odpowiedź: ηmax = 58,01 %
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 = 425 + 273,15 = 698,1 K
TL = 20 + 273,15 = 293,1 K
Sprawność Carnota
ηmax = 1 − TL/TH = 1 − 293,1/698,1 = 58,01 %
Odpowiedź
ηmax = 58,01 %
22.Silniki cieplne i chłodziarki
Odpowiedź: Ẇmin = 1,33 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 = 13/9,772 = 1,33 kW
Odpowiedź
Ẇmin = 1,33 kW
23.Obiegi silników cieplnych gazowych
Odpowiedź: T3 = 1430 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 = 300·180,4 = 953,3 K
Doprowadzenie ciepła przy stałym ciśnieniu
T3 = rc·T2 = 1,5·953,3 = 1430 K
Odpowiedź
T3 = 1430 K
24.Obiegi silników cieplnych gazowych
Odpowiedź: wnet = 249,4 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·60,2857 = 480,5 K
T4 = 1100/60,2857 = 659,3 K
Praca sprężarki i turbiny
wc = 1,005·(480,5 − 288) = 193,5 kJ/kg
wt = 1,005·(1100 − 659,3) = 442,9 kJ/kg
Praca netto
wnet = 442,9 − 193,5 = 249,4 kJ/kg
Odpowiedź
wnet = 249,4 kJ/kg
Wskaźnik pracy zwrotnej: wc/wt = 43,68%.
25.Ciągłość i Bernoulli
Odpowiedź: p2 = 70 kPa
Ciągłość i Bernoulli: dla płynu nieściśliwego natężenie przepływu jest takie samo w każdym przekroju, Q = A₁v₁ = A₂v₂. Wzdłuż linii prądu bez strat p + ½ρv² + ρgz jest stałe, więc tam, gdzie płyn jest szybszy lub wyżej, ciśnienie jest niższe.
Ciągłość i Bernoulli
Licz w Pa, a potem podziel przez 1000, aby dostać kPa.
v2 = v1·(d1/d2)² = 1·(75/25)² = 9 m/s
½ρ(v1² − v2²) = ½·1000·(1² − 9²)/1000 = −40 kPa
p2 = p1 + ½ρ(v1² − v2²) = 110 + (−40) = 70 kPa
Odpowiedź
p2 = 70 kPa
26.Ciągłość i Bernoulli
Odpowiedź: Q = 72,33 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.
Różnica ciśnień z manometru
Δp = (13600 − 1000)·9,81·0,14 = 17305 Pa
Prędkość w gardzieli
Równanie ciągłości i równanie Bernoulliego razem dają idealną prędkość w gardzieli.
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,05²/4 = 0,001963 m²
V = Q/A = 0,005/0,001963 = 2,546 m/s
Wysokość prędkości
V²/(2g) = 2,546²/(2·9,81) = 0,3305 m
Darcy–Weisbach
hf = f·(L/D)·V²/(2g) = 0,03·(300/0,05)·2,546²/(2·9,81) = 59,49 m
Odpowiedź
hf = 59,49 m
28.Tarcie w rurach i straty wysokości
Odpowiedź: Q = 240,3 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,018·150/0,25 = 10,8
ΣK = 0,5 + 2·0,3 + 0,2 + 1,0 = 2,3
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,03·54 = 15,89 kW
Moc na wale
Pshaft = Ph/η = 15,89/0,6 = 26,49 kW
Odpowiedź
Pshaft = 26,49 kW
30.Pompy i moc pomp
Odpowiedź: H = 49,62 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 = 7 + 35,68 + 6,941 = 49,62 m
Odpowiedź
H = 49,62 m
31.Przewodzenie ciepła
Odpowiedź: Q̇ = 41,69 W
Opór cieplny: ustalony przepływ ciepła przez warstwy połączone szeregowo zachowuje się jak prąd w rezystorach szeregowych: Q̇ = ΔT/ΣR, gdzie R = L/(kA) dla warstwy przewodzącej i R = 1/(hA) dla powierzchni z konwekcją. Spadek temperatury na każdej warstwie jest proporcjonalny do jej oporu.
Opory warstw
R1 = L1/(k1A) = (250/1000)/(0,17·8) = 0,1838 K/W
R2 = L2/(k2A) = (150/1000)/(0,035·8) = 0,5357 K/W
Strumień ciepła
Q̇ = (T1 − T3)/(R1 + R2) = 30/0,7195 = 41,69 W
Odpowiedź
Q̇ = 41,69 W
32.Przewodzenie ciepła
Odpowiedź: Q̇ = 289,6 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.
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π·1500/60 = 157,1 rad/s
T1 = P/ω₁ = 1500/157,1 = 9,549 N·m
Moment wyjściowy
T2 = 9,549·80/16·0,96 = 45,84 N·m
Odpowiedź
T2 = 45,84 N·m
34.Koła i przekładnie zębate
Odpowiedź: Tout = 142,3 N·m
Gear trains: meshing gears have the same speed at the teeth, so n₁·N₁ = n₂·N₂: the gear with fewer teeth turns faster. Gears on one shaft turn at the same speed. Power is passed on (minus losses), so torque goes up when speed goes down: T₂ = T₁·(N₂/N₁)·η, and P = T·ω with ω = 2π·n/60.
Prędkość obrotowa wyjściowa
nout = 2900·(20/40)·(16/25) = 928 rpm
ωout = 2π·928/60 = 97,18 rad/s
Moc i moment wyjściowy
Pout = 15000·0,96² = 13824 W
Tout = Pout/ωout = 13824/97,18 = 142,3 N·m
Odpowiedź
Tout = 142,3 N·m
35.Transformaty Laplace’a
Odpowiedź: f(0.1) = 0,1391
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).
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).
Równanie charakterystyczne
1 + K/(s(s + 3)) = 0 → s² + 3s + K = 0
Pierwiastek podwójny
Trójmian s² + βs + γ ma pierwiastek podwójny, gdy β² = 4γ.
3² = 4K → K = 9/4 = 2,25
Odpowiedź
K = 2,25
Oba bieguny układu zamkniętego leżą wtedy w s = −1,5.
38.Schematy blokowe
Odpowiedź: K = 49
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)) = 5·10/(10 + 10K)
Rozwiąż względem K
10 + 10K = 5·10/0,1 = 500
K = (500 − 10)/10 = 49
Odpowiedź
K = 49
39.Odpowiedź skokowa
Odpowiedź: %OS = 11,51 %
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 + ωₙ².
ωₙ² = 7 → ωₙ = 2,646 rad/s
2ζωₙ = 3 → ζ = 3/(2·2,646) = 0,5669
Wynik
%OS = 100·e−ζπ/√(1−ζ²) = 100·e−2,162 = 11,51 %
Odpowiedź
%OS = 11,51 %
40.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).
Dla tej pętli ζωₙ = a/2 niezależnie od K: zwiększenie K zwiększa przeregulowanie, ale nie skraca czasu ustalania.
41.Uchyb ustalony i regulacja PID
Odpowiedź: ess = 0,25
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) = 4
Część skokowa
Kp = ∞ → estep = 0
Uchyb dla narastania liniowego
ess = 1/4 = 0,25
Odpowiedź
ess = 0,25
42.Uchyb ustalony i regulacja PID
Odpowiedź: ess = 25
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).
Stała prędkościowa
s·C(s) = 1s + 1 → 1 przy s → 0.
Kv = 1·2/10 = 0,2
Uchyb dla narastania liniowego
ess = 5/0,2 = 25
Odpowiedź
ess = 25
Regulator PI całkowicie usuwa uchyb skokowy; narastanie liniowe jest śledzone ze stałym opóźnieniem.
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.