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 RB.
L = 10 mP1 = 20 kN w 5 mP2 = 60 kN w 9 m
RB = kN
2.Reakcje podpór belekTrudny
Belka złożona jest utwierdzona w ścianie w A, ma przegub wewnętrzny w C (5 m od A) i opiera się na podporze przesuwnej w D (10 m od A). Przenosi siłę skupioną P w odległości 6 m od A i siłę skupioną Q w odległości 1 m od A. Wyznacz reakcję RD na podporze przesuwnej.
AC = 5 mAD = 10 mP = 25 kN w 6 mQ = 60 kN w 1 m
RD = kN
3.Naprężenie i odkształcenie osioweŚredni
Pręt składa się z odcinka (materiał: stal, średnica 40 mm, długość 400 mm, E = 200 GPa) połączonego końcami z odcinkiem (materiał: miedź, średnica 60 mm, długość 400 mm, E = 120 GPa). Jest zamocowany na jednym końcu i rozciągany na drugim. Wyznacz obciążenie osiowe P, które wydłuża go łącznie o 0,25 mm.
odcinek 1: stal, Ø40 mm, 400 mm, E = 200 GPaodcinek 2: miedź, Ø60 mm, 400 mm, E = 120 GPaδ = 0,25 mm
P = kN
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 = 50 kNodcinek 1: Ø30 mm, 2 m, mosiądz (E = 100 GPa)odcinek 2: Ø20 mm, 0,5 m, miedź (E = 117 GPa)
σmax = MPa
5.Skręcanie wałówŚredni
Pełny wał (stal) o długości 1,5 m przenosi moment skręcający 500 N·m. Kąt skręcenia między jego końcami nie może przekroczyć 0,5°. Wyznacz najmniejszą średnicę d, która spełnia ten warunek.
T = 500 N·mL = 1,5 mG = 80 GPaφallow = 0,5°
d = mm
6.Skręcanie wałówTrudny
Pełny wał ma przenosić 75 kW przy 3000 obr/min bez przekroczenia naprężenia stycznego 50 MPa. Wyznacz najmniejszą średnicę d, która wystarczy.
P = 75 kWN = 3000 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.
15 kN @ 2 m50 kN @ 8 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)
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.
45 kN @ 0 m60 kN @ 6 m50 kN @ 3 m10 kN/m, 1–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)
9.Koło Mohra i transformacja naprężeńŚredni
W pewnym punkcie obciążonej części występują σₓ = 110 MPa, σᵧ = −10 MPa i τₓᵧ = −20 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 = 110 MPaσy = −10 MPaτxy = −20 MPa
θp = °
10.Koło Mohra i transformacja naprężeńTrudny
Punkt cienkiej płyty jest w płaskim stanie naprężenia: σₓ = 10 MPa, σᵧ = −50 MPa i τₓᵧ = −30 MPa (σz = 0). Granica plastyczności materiału wynosi 300 MPa. Wyznacz współczynnik bezpieczeństwa n względem uplastycznienia według hipotezy Hubera-Misesa (energii odkształcenia postaciowego).
σx = 10 MPaσy = −50 MPaτxy = −30 MPaσY = 300 MPa
n =
11.Ugięcie belkiŚredni
Wspornik przenosi obciążenie równomiernie rozłożone w na całej długości i siłę skupioną P na swobodnym końcu. Korzystając z superpozycji, wyznacz ugięcie δ swobodnego końca.
P = 10 kNw = 12 kN/mL = 1,5 mE = 70 GPa (aluminium)I = 8,5 × 10⁶ mm⁴
δ = mm
12.Ugięcie belkiTrudny
Belka swobodnie podparta przenosi obciążenie równomiernie rozłożone w. Jej ugięcie nie może przekroczyć rozpiętość/360. Wyznacz najmniejszy moment bezwładności Ireq, który spełnia ten warunek.
w = 8 kN/mL = 5 mE = 70 GPa (aluminium)
Ireq = × 10⁶ mm⁴
13.Wyboczenie słupówŚredni
Słup o długości L z materiału (stal) jest przegubowo podparty na obu końcach. Korzystając ze wzoru Eulera, wyznacz jego siłę krytyczną Pcr.
L = 6 mE = 200 GPa (stal)b × h = 20 × 100 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 utwierdzony u podstawy i przegubowo podparty u góry. Ma przenieść 55 kN przy współczynniku bezpieczeństwa 2,5 względem wyboczenia. Wyznacz najmniejszą średnicę d.
P = 55 kNn = 2,5L = 4 mE = 200 GPa (stal)
d = mm
15.Rzut ukośnyŚredni
Piłka zostaje wyrzucona pod kątem 55° do poziomu z płaskiego terenu. Pomiń opór powietrza. Wyznacz czas lotu t, zanim piłka uderzy w ziemię.
v0 = 36 m/sθ = 55°g = 9,81 m/s²
t = s
16.Rzut ukośnyTrudny
Piłkę wyrzucono z prędkością 33 m/s pod kątem 25° do poziomu ze szczytu klifu o wysokości 35 m. Pomiń opór powietrza i przyjmij g = 9,81 m/s². Wyznacz największą wysokość H nad ziemią u podnóża klifu.
v₀ = 33 m/sθ = 25°h = 35 mg = 9,81 m/s²
H = m
17.Praca i energiaŚredni
Klocek o masie 5 kg zsuwa się na odcinku 18 m po zboczu nachylonym pod kątem 45° do poziomu, przy współczynniku tarcia kinetycznego 0,2. Wyznacz pracę Wf wykonaną przeciw tarciu.
m = 5 kgθ = 45°μ = 0,2d = 18 mg = 9,81 m/s²
Wf = J
18.Praca i energiaTrudny
Klocek o masie 4 kg leży na chropowatym poziomym stole (μ = 0,4). Lekka nić biegnie od niego przez gładki bloczek na krawędzi stołu do wiszącego klocka o masie 5 kg. Układ zostaje puszczony ze spoczynku. Wyznacz prędkość v klocków po opadnięciu wiszącego klocka o 3 m.
m1 = 4 kgm2 = 5 kgμ = 0,4d = 3 mg = 9,81 m/s²
v = m/s
19.Drgania mechaniczneŚredni
Masa 20 kg na sprężynie o sztywności 2 kN/m ma tłumik wiskotyczny dający współczynnik tłumienia 0,3. Wyznacz częstotliwość drgań swobodnych tłumionych fd w Hz.
m = 20 kgk = 2 kN/mζ = 0,3
fd = Hz
20.Drgania mechaniczneTrudny
Masa 50 kg wisi na dwóch sprężynach połączonych jedna pod drugą (szeregowo): górna sprężyna ma k₁ = 50 kN/m, a dolna k₂ = 25 kN/m. Wyznacz sztywność zastępczą keq układu sprężyn w kN/m.
m = 50 kgk1 = 50 kN/mk2 = 25 kN/m
keq = kN/m
21.Silniki cieplne i chłodziarkiŚredni
Silnik cieplny pracuje między zbiornikiem gorącym o temperaturze 225 °C a zbiornikiem zimnym o temperaturze 30 °C. Wyznacz największą możliwą sprawność cieplną ηmax.
TH = 225 °CTL = 30 °C
ηmax = %
22.Silniki cieplne i chłodziarkiTrudny
Pompa ciepła dostarcza 9 kW ciepła do instalacji grzewczej o temperaturze 55 °C, pobierając ciepło ze źródła o temperaturze 10 °C. Jej rzeczywisty COP wynosi 60 % COP Carnota między tymi temperaturami. Wyznacz pobieraną przez nią moc elektryczną Ẇ.
TH = 55 °CTL = 10 °CηII = 60%Q̇H = 9 kW
Ẇ = kW
23.Obiegi silników cieplnych gazowychŚredni
Idealny obieg Diesla ma stopień sprężania r = 16 i stopień obciążenia rc = v₃/v₂ = 1,5. 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 = 16rc = 1,5k = 1,4
η = %
24.Obiegi silników cieplnych gazowychTrudny
Elektrownia z turbiną gazową pracuje w obiegu Braytona przy stosunku ciśnień w sprężarce 8. Powietrze wpływa do sprężarki w temperaturze 290 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 = 8T1 = 290 KT3 = 1100 Kk = 1,4cp = 1,005 kJ/(kg·K)
wnet = kJ/kg
25.Ciągłość i BernoulliŚredni
Rurka Prandtla w rurociągu wodnym wskazuje różnicę 15 kPa między ciśnieniem spiętrzenia a ciśnieniem statycznym. Wyznacz prędkość wody v przy końcówce.
p0 − p = 15 kPaρ = 1000 kg/m³
v = m/s
26.Ciągłość i BernoulliTrudny
Woda płynie w górę przez rurę, która zwęża się z 150 mm do 75 mm i wznosi się o 10 m. W przekroju 1 ciśnienie wynosi 330 kPa, a prędkość 1,5 m/s. Pomijając straty, wyznacz ciśnienie p2 w przekroju 2.
Woda (ρ = 1000 kg/m³) płynie z prędkością 0,8 m/s przez poziomą rurę o długości 800 m i średnicy wewnętrznej 150 mm. Przyjmij współczynnik tarcia Darcy’ego f = 0,03 i g = 9,81 m/s². Wyznacz spadek ciśnienia Δp na rurze w kPa.
Rura o długości 150 m i średnicy wewnętrznej 150 mm (współczynnik tarcia Darcy’ego f = 0,022) wychodzi z otwartego zbiornika i wypływa swobodnym strumieniem do atmosfery. Straty miejscowe to: wlot o ostrej krawędzi (K = 0,5), kolana w liczbie 4 (K = 0,3 każde) i zawór (K = 5). Przyjmij g = 9,81 m/s². O ile H zwierciadło wody w zbiorniku musi leżeć wyżej niż wylot rury, aby natężenie przepływu wynosiło 60 L/s?
Podczas próby pompa tłoczy 60 L/s wody (ρ = 1000 kg/m³) przy wysokości podnoszenia 36 m, a jej wał pobiera 24,9 kW. Przyjmij g = 9,81 m/s². Wyznacz sprawność pompy η w procentach.
Q = 60 L/sH = 36 mPshaft = 24,9 kWg = 9,81 m/s²
η = %
30.Pompy i moc pompTrudny
Pompa przetłacza 120 L/s wody (ρ = 1000 kg/m³) z otwartego dolnego zbiornika do górnego zbiornika, którego zwierciadło wody leży 14 m wyżej. Rura ma długość 400 m, średnicę wewnętrzną 250 mm i współczynnik tarcia Darcy’ego f = 0,02; współczynniki strat miejscowych sumują się do ΣK = 5 (łącznie z wlotem i wylotem). Przyjmij g = 9,81 m/s². Górny zbiornik jest zamknięty, a powietrze nad wodą ma nadciśnienie 250 kPa. Wyznacz wysokość podnoszenia H, którą musi zapewnić pompa.
Ściana składa się z 50 mm (szkło) i 100 mm (beton), a jej pole wynosi 12 m². Powierzchnie zewnętrzne mają 21 °C i -11 °C. Wyznacz strumień ciepła Q̇ przez ścianę.
L1 = 50 mm, k1 = 0,96 W/(m·K)L2 = 100 mm, k2 = 1,4 W/(m·K)A = 12 m²T1 = 21 °CT3 = -11 °C
Q̇ = W
32.Przewodzenie ciepłaTrudny
Ściana z 50 mm (szkło) i 50 mm (płyta gipsowo-kartonowa) oddziela pomieszczenie o temperaturze 21 °C od powietrza zewnętrznego o temperaturze -9 °C. Współczynniki przejmowania ciepła wynoszą 10 W/(m²·K) wewnątrz i 20 W/(m²·K) na zewnątrz, a pole ściany 6 m². Wyznacz straty ciepła Q̇ przez ścianę.
L1 = 50 mm, k1 = 0,96 W/(m·K)L2 = 50 mm, k2 = 0,17 W/(m·K)A = 6 m²hi = 10 W/(m²·K)ho = 20 W/(m²·K)Troom = 21 °CTout = -9 °C
Q̇ = W
33.Koła i przekładnie zębateŚredni
Silnik o mocy 0,75 kW pracujący z prędkością 1440 rpm napędza zębnik o 20 zębach, który zazębia się z kołem o 60 zębach. Pomiń straty. Wyznacz moment wyjściowy T2 na wale koła.
P = 0,75 kWn1 = 1440 rpmN1 = 20N2 = 60η = 100%
T2 = N·m
34.Koła i przekładnie zębateTrudny
Pas otwarty łączy koło napędzające 80 mm o prędkości 1450 obr/min z kołem 200 mm; środki są oddalone o 500 mm. Współczynnik tarcia wynosi 0,3, a siła w cięgnie czynnym nie może przekroczyć 800 N. Wyznacz największą moc P, jaką pas może przenieść bez poślizgu, w kW (pomiń efekty odśrodkowe).
F(s) = (4s + 22)/(s² + 11s + 28). Jej transformata odwrotna ma postać f(t) = A·e−4t + B·e−7t. Znajdź współczynnik A metodą rozkładu na ułamki proste.
F(s) = (4s + 22)/(s² + 11s + 28)
A =
36.Transformaty Laplace’aTrudny
Rozwiąż y″ + 7y′ + 10y = 0 dla y(0) = 9 i y′(0) = 12 za pomocą transformaty Laplace’a i znajdź y(0,5).
y″ + 7y′ + 10y = 0y(0) = 9y′(0) = 12t = 0,5 s
y(0.5) =
37.Schematy blokoweŚredni
Obiekt G(s) = 16/(s + 4) jest regulowany w pętli ujemnego sprzężenia zwrotnego, w której czujnik ma wzmocnienie H = 0,1. Na wejście R podano skok o wielkości 1. Znajdź wyjście w stanie ustalonym css.
G(s) = 16/(s + 4)H = 0,1R = 1
css =
38.Schematy blokoweTrudny
Regulator proporcjonalny K = 25 steruje obiektem G(s) = 3/(s + 3) w pętli jednostkowego ujemnego sprzężenia zwrotnego. Zakłócenie D wchodzi na wejście obiektu (dodaje się do wyjścia regulatora). Przy R = 0 podano skok zakłócenia o wielkości 0,5. Znajdź wyjście w stanie ustalonym css.
K = 25G(s) = 3/(s + 3)D = 0,5R = 0
css =
39.Odpowiedź skokowaŚredni
Układ ma transmitancję G(s) = 186/(s² + 10s + 93) i otrzymuje jednostkowe wymuszenie skokowe. Bieguny to s = σ ± jωd. Znajdź część rzeczywistą σ.
G(s) = 186/(s² + 10s + 93)
σ =
40.Odpowiedź skokowaTrudny
Pętla jednostkowego ujemnego sprzężenia zwrotnego ma G(s) = K/(s(s + 8)). Znajdź wzmocnienie K, które daje odpowiedź skokową z przeregulowaniem 20 %.
G(s) = K/(s(s + 8))%OS = 20 %
K =
41.Uchyb ustalony i regulacja PIDŚredni
Pętla ze sprzężeniem jednostkowym ma G(s) = K/((s + 1)(s + 10)). Wyznacz wzmocnienie K, które daje uchyb ustalony 4 % dla wejścia skokowego.
G(s) = K/((s + 1)(s + 10))ess = 4 %
K =
42.Uchyb ustalony i regulacja PIDTrudny
Regulator proporcjonalny K = 25 steruje obiektem G(s) = 2/(s(s + 8)) w pętli ze sprzężeniem jednostkowym. Znajdź czas szczytu tp odpowiedzi skokowej.
G(s) = 2/(s(s + 8))K = 25
tp = s
3Roztwory
1.Reakcje podpór belek
Odpowiedź: RB = 64 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·10 = 20·5 + 60·9
RB = 640/10 = 64 kN
ΣFy = 0: RA = 80 − 64 = 16 kN
Odpowiedź
RB = 64 kN
2.Reakcje podpór belek
Odpowiedź: RD = 5 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·5 = 25·1
RD = 5 kN
VC = 25 − 5 = 20 kN
Odpowiedź
RD = 5 kN
3.Naprężenie i odkształcenie osiowe
Odpowiedź: P = 90,24 kN
Dwa materiały: przy połączeniu szeregowym (jeden za drugim) każdy odcinek przenosi tę samą siłę, a wydłużenia się dodają: δ = P·Σ L/(E·A). Przy połączeniu równoległym (obok siebie między sztywnymi płytami) obie części skracają się o tyle samo, więc ich odkształcenia są równe: sztywniejsza część przejmuje większy udział, Pi = P·EiAi/ΣEA, a σi = Ei·ε. Pręt zamocowany między dwiema sztywnymi ścianami z obciążeniem w miejscu połączenia jest statycznie niewyznaczalny: dwa odcinki dzielą obciążenie proporcjonalnie do sztywności k = EA/L.
Pola
A₁ = 1257 mm², A₂ = 2827 mm²
Podatności
L₁/(E₁A₁) = 400/(200000·1257) = 1,592 × 10⁻⁶ mm/N
L₂/(E₂A₂) = 400/(120000·2827) = 1,179 × 10⁻⁶ mm/N
Obciążenie
P = 0,25/2,77 × 10⁻⁶ = 90237 N = 90,24 kN
Odpowiedź
P = 90,24 kN
4.Naprężenie i odkształcenie osiowe
Odpowiedź: σmax = 159,2 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.
Cieńszy odcinek ma mniejsze pole, więc przenosi większe naprężenie.
σmax = 159,2 MPa
Odpowiedź
σmax = 159,2 MPa
5.Skręcanie wałów
Odpowiedź: d = 57,51 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·ω.
W praktyce zaokrągla się w górę do najbliższego wymiaru znormalizowanego.
6.Skręcanie wałów
Odpowiedź: d = 28,97 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π·3000/60 = 314,2 rad/s
T = P/ω = 75000/314,2 = 238,7 N·m
Średnica z naprężenia dopuszczalnego
d = ∛(16T/(π·τallow)) = ∛(16·238732/(π·50)) = 28,97 mm
Odpowiedź
d = 28,97 mm
W praktyce zaokrągla się w górę do najbliższego wymiaru znormalizowanego.
7.Siła tnąca i moment zginający
Odpowiedź: E (x = 8 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 = 29,17 kN, R2 = 35,83 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 = 58,33, MC = 86,67, MD = 115, ME = 143,3 kN·m, MF = 71,67, MG = 0
Odpowiedź
Punkt E: ME = 143,3 kN·m
8.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 = 51,67 kN, R2 = 133,3 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 = −45, MC = −43,33, MD = −51,67, ME = −120 kN·m, MF = −60, MG = 0
Odpowiedź
Punkt E: ME = −120 kN·m
9.Koło Mohra i transformacja naprężeń
Odpowiedź: θp = −9,217 °
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·−20, 110 − −10) = −18,43°
Kąt na elemencie
θp = −9,217°
Odpowiedź
θp = −9,217°
Kąty dodatnie mierzy się przeciwnie do ruchu wskazówek zegara.
10.Koło Mohra i transformacja naprężeń
Odpowiedź: n = 3,939
Transformacja płaskiego stanu naprężenia: na płaszczyźnie obróconej o kąt θ: σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ. Koło Mohra pokazuje wszystkie płaszczyzny naraz: jego środek leży w C = (σₓ + σᵧ)/2, a promień wynosi R = √(((σₓ − σᵧ)/2)² + τₓᵧ²). Naprężenia główne to C ± R, a największe naprężenie styczne w płaszczyźnie to R. Hipotezy wytężeniowe: według hipotezy Hubera-Misesa uplastycznienie następuje, gdy σvM = √(σ₁² − σ₁σ₂ + σ₂²) osiąga σY; według hipotezy Tresci, gdy bezwzględne maksymalne naprężenie styczne osiąga σY/2.
n > 1 oznacza, że punkt się nie uplastycznia; hipoteza Tresci jest nieco bardziej zachowawcza niż hipoteza Hubera-Misesa.
11.Ugięcie belki
Odpowiedź: δ = 31,67 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 = 10000 N, w = 12 N/mm, L = 1500 mm
EI = 5,95 × 10¹¹ N·mm²
Każde obciążenie osobno
δP = PL³/(3EI) = 10000·1500³/(3·5,95 × 10¹¹) = 18,91 mm
δw = wL⁴/(8EI) = 12·1500⁴/(8·5,95 × 10¹¹) = 12,76 mm
Dodaj
δ = 18,91 + 12,76 = 31,67 mm
Odpowiedź
δ = 31,67 mm
12.Ugięcie belki
Odpowiedź: Ireq = 66,96 × 10⁶ mm⁴
Ugięcie belek: ugięcia otrzymuje się, całkując dwukrotnie krzywiznę M/(EI). Typowe przypadki są stablicowane, na przykład PL³/(3EI) dla wspornika z siłą na końcu i 5wL⁴/(384EI) dla belki swobodnie podpartej z obciążeniem równomiernym na całej rozpiętości. Ugięcia od różnych obciążeń się sumują (superpozycja).
Warunek ugięcia
δallow = L/360 = 5000/360 = 13,89 mm
Wymagany moment bezwładności
I = 5wL⁴/(384·E·δallow) = 5·8·5000⁴/(384·70000·13,89) = 6,696 × 10⁷ mm⁴
Ireq = 66,96 × 10⁶ mm⁴
Odpowiedź
Ireq = 66,96 × 10⁶ mm⁴
13.Wyboczenie słupów
Odpowiedź: Pcr = 3,655 kN
Wyboczenie Eulera: smukły słup traci stateczność przy Pcr = π²EI/(KL)², gdzie KL to długość wyboczeniowa: K = 1 dla przegub–przegub, 2 dla utwierdzenie–wolny koniec, 0,7 dla utwierdzenie–przegub i 0,5 dla utwierdzenie–utwierdzenie. Wybocza się względem osi o najmniejszym I. Wzór dotyczy słupów smukłych, których naprężenie krytyczne jest mniejsze od granicy plastyczności.
Najmniejszy moment bezwładności
Prostokąt wybocza się względem słabszej osi, z mniejszym wymiarem do sześcianu.
Wyboczenie Eulera: smukły słup traci stateczność przy Pcr = π²EI/(KL)², gdzie KL to długość wyboczeniowa: K = 1 dla przegub–przegub, 2 dla utwierdzenie–wolny koniec, 0,7 dla utwierdzenie–przegub i 0,5 dla utwierdzenie–utwierdzenie. Wybocza się względem osi o najmniejszym I. Wzór dotyczy słupów smukłych, których naprężenie krytyczne jest mniejsze od granicy plastyczności.
Wymagana siła krytyczna
Pcr = n·P = 2,5·55000 = 137500 N
KL = 0,7·4000 = 2800 mm
Wyznacz średnicę
d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·137500·2800²/(π³·200000))^(1/4) = 57,75 mm
Odpowiedź
d = 57,75 mm
Sprawdzenie: smukłość KL/r wynosi 193,9, więc słup jest smukły i wzór Eulera obowiązuje.
15.Rzut ukośny
Odpowiedź: t = 6,012 s
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 θ = 36·cos 55° = 20,65 m/s
vy = v0·sin θ = 36·sin 55° = 29,49 m/s
Czas lotu
Start i lądowanie na tej samej wysokości: t = 2vy/g.
t = 2·29,49/9,81 = 6,012 s
Odpowiedź
t = 6,012 s
16.Rzut ukośny
Odpowiedź: H = 44,91 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.
Pionowa prędkość początkowa
vy = 33·sin 25° = 13,95 m/s
Wzniesienie ponad klif
Δh = vy²/(2g) = 13,95²/(2·9,81) = 9,913 m
Wysokość nad ziemią
H = 35 + 9,913 = 44,91 m
Odpowiedź
H = 44,91 m
17.Praca i energia
Odpowiedź: Wf = 124,9 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,7071 = 34,68 N
F = μN = 0,2·34,68 = 6,937 N
Praca przeciw tarciu
Wf = F·d = 6,937·18 = 124,9 J
Odpowiedź
Wf = 124,9 J
18.Praca i energia
Odpowiedź: v = 4,716 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 = 49,05·3 = 147,2 J
μ·m₁g·d = 0,4·39,24·3 = 47,09 J
Energia kinetyczna obu klocków
½(m₁ + m₂)v² = 100,1 → v = √(2·100,1/9) = 4,716 m/s
Odpowiedź
v = 4,716 m/s
19.Drgania mechaniczne
Odpowiedź: fd = 1,518 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₂.
Częstotliwość drgań własnych nietłumionych
ωn = √(2000/20) = 10 rad/s
fn = 10/(2π) = 1,592 Hz
Częstotliwość drgań tłumionych
fd = fn·√(1 − ζ²) = 1,592·√(1 − 0,3²) = 1,518 Hz
Odpowiedź
fd = 1,518 Hz
20.Drgania mechaniczne
Odpowiedź: keq = 16,67 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 szeregowe przenoszą tę samą siłę, a ich ugięcia się sumują, więc łączy się je przez odwrotności: 1/k = 1/k1 + 1/k2.
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 = 225 + 273,15 = 498,1 K
TL = 30 + 273,15 = 303,1 K
Sprawność Carnota
ηmax = 1 − TL/TH = 1 − 303,1/498,1 = 39,14 %
Odpowiedź
ηmax = 39,14 %
22.Silniki cieplne i chłodziarki
Odpowiedź: Ẇ = 2,057 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 = 55 + 273,15 = 328,1 K
TL = 10 + 273,15 = 283,1 K
COP Carnota i rzeczywisty
COPmax = 328,1/(328,1 − 283,1) = 7,292
COP = 0,6·7,292 = 4,375
Moc pobierana
Ẇ = Q̇H/COP = 9/4,375 = 2,057 kW
Odpowiedź
Ẇ = 2,057 kW
23.Obiegi silników cieplnych gazowych
Odpowiedź: η = 63,99 %
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
(1,51,4 − 1)/(1,4·(1,5 − 1)) = 1,092
Sprawność cieplna
η = 1 − 0,3299·1,092 = 63,99%
Odpowiedź
η = 63,99%
24.Obiegi silników cieplnych gazowych
Odpowiedź: wnet = 258,7 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 = 290·80,2857 = 525,3 K
T4 = 1100/80,2857 = 607,2 K
Praca sprężarki i turbiny
wc = 1,005·(525,3 − 290) = 236,5 kJ/kg
wt = 1,005·(1100 − 607,2) = 495,2 kJ/kg
Praca netto
wnet = 495,2 − 236,5 = 258,7 kJ/kg
Odpowiedź
wnet = 258,7 kJ/kg
Wskaźnik pracy zwrotnej: wc/wt = 47,76%.
25.Ciągłość i Bernoulli
Odpowiedź: v = 5,477 m/s
Continuity and Bernoulli: for an incompressible fluid the flow rate is the same at every section, Q = A₁v₁ = A₂v₂. Along a streamline without losses, p + ½ρv² + ρgz is constant, so where the fluid is faster or higher its pressure is lower.
Równanie Bernoulliego do punktu spiętrzenia
Przepływ zatrzymuje się przy końcówce, więc jego energia kinetyczna ujawnia się jako dodatkowe ciśnienie.
p0 − p = ½ρv²
v = √(2·15000/1000) = 5,477 m/s
Odpowiedź
v = 5,477 m/s
26.Ciągłość i Bernoulli
Odpowiedź: p2 = 215 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.
Przepływ w rurach: liczba Reynoldsa Re = V·D/ν = ρ·V·D/μ odróżnia przepływ laminarny (Re poniżej ok. 2000) od turbulentnego (powyżej ok. 4000). Tarcie wzdłuż rury powoduje stratę wysokości Darcy’ego–Weisbacha hf = f·(L/D)·V²/(2g). Kształtki, zawory, wlot i wylot dodają straty miejscowe hm = ΣK·V²/(2g). Między dwiema otwartymi powierzchniami zbiorników równanie energii daje Δz = hf + hm, z uwzględnioną stratą na wylocie K = 1.
Strata wysokości
hf = f·(L/D)·V²/(2g) = 0,03·(800/0,15)·0,8²/(2·9,81) = 5,219 m
Spadek ciśnienia
Rura jest pozioma, a prędkość się nie zmienia, więc równanie energii upraszcza się do Δp = ρ·g·hf.
Przepływ w rurach: liczba Reynoldsa Re = V·D/ν = ρ·V·D/μ odróżnia przepływ laminarny (Re poniżej ok. 2000) od turbulentnego (powyżej ok. 4000). Tarcie wzdłuż rury powoduje stratę wysokości Darcy’ego–Weisbacha hf = f·(L/D)·V²/(2g). Kształtki, zawory, wlot i wylot dodają straty miejscowe hm = ΣK·V²/(2g). Między dwiema otwartymi powierzchniami zbiorników równanie energii daje Δz = hf + hm, z uwzględnioną stratą na wylocie K = 1.
Prędkość
V = Q/A = 0,06/0,01767 = 3,395 m/s
V²/(2g) = 0,5876 m
Współczynniki strat
f·L/D = 0,022·150/0,15 = 22
ΣK = 0,5 + 4·0,3 + 5 = 6,7
Równanie energii
Strumień na wylocie nadal niesie swoją wysokość prędkości; to nie strata, ale trzeba ją dostarczyć.
H = (1 + 22 + 6,7)·0,5876 = 17,45 m
Odpowiedź
H = 17,45 m
29.Pompy i moc pomp
Odpowiedź: η = 85,1 %
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,06·36 = 21,19 kW
Sprawność
η = Ph/Pshaft = 21,19/24,9 = 85,1%
Odpowiedź
η = 85,1%
30.Pompy i moc pomp
Odpowiedź: H = 50,75 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 = 14 + 25,48 + 11,27 = 50,75 m
Odpowiedź
H = 50,75 m
31.Przewodzenie ciepła
Odpowiedź: Q̇ = 3109 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.
Opory szeregowo
Ri = 1/(hiA) = 1/(10·6) = 0,01667 K/W
R1 = L1/(k1A) = (50/1000)/(0,96·6) = 0,008681 K/W
R2 = L2/(k2A) = (50/1000)/(0,17·6) = 0,04902 K/W
Ro = 1/(hoA) = 1/(20·6) = 0,008333 K/W
ΣR = 0,0827 K/W
Straty ciepła
Q̇ = (21 − (-9))/0,0827 = 362,8 W
Odpowiedź
Q̇ = 362,8 W
33.Koła i przekładnie zębate
Odpowiedź: T2 = 14,92 N·m
Przekładnie zębate: zazębione koła mają tę samą prędkość na zębach, więc n₁·N₁ = n₂·N₂: koło o mniejszej liczbie zębów obraca się szybciej. Koła na jednym wale obracają się z tą samą prędkością obrotową. Moc jest przekazywana dalej (minus straty), więc moment obrotowy rośnie, gdy prędkość obrotowa maleje: T₂ = T₁·(N₂/N₁)·η, a P = T·ω, gdzie ω = 2π·n/60.
Moment silnika
ω₁ = 2π·1440/60 = 150,8 rad/s
T1 = P/ω₁ = 750/150,8 = 4,974 N·m
Moment wyjściowy
T2 = 4,974·60/20·1 = 14,92 N·m
Odpowiedź
T2 = 14,92 N·m
34.Koła i przekładnie zębate
Odpowiedź: P = 2,824 kW
Przekładnie pasowe: bez poślizgu oba koła pasowe mają tę samą prędkość obwodową v = π·d·n/60, więc n₁·d₁ = n₂·d₂, a przełożenie wynosi d₂/d₁. Poślizg s sprawia, że koło napędzane obraca się nieco wolniej: n₂ = n₁·(d₁/d₂)·(1 − s). Przenoszona moc to P = (T₁ − T₂)·v, gdzie T₁ to siła w cięgnie czynnym, a T₂ w biernym. Tarcie ogranicza stosunek sił: T₁/T₂ = eμθ, gdzie θ jest kątem opasania w radianach na mniejszym kole; dla pasa otwartego θ = π − 2·asin((d₂ − d₁)/(2C)) przy rozstawie osi C.
Kąt opasania
θ = 2,901 rad
Siła w cięgnie biernym
eμθ = e0,8703 = 2,388, T₂ = 800/2,388 = 335,1 N
Moc
P = (800 − 335,1)·6,074 = 2,824 kW
Odpowiedź
P = 2,824 kW
35.Transformaty Laplace’a
Odpowiedź: A = 2
Transformata Laplace’a: F(s) = ∫₀∞ f(t)·e−st dt zamienia pochodne na mnożenie przez s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Podstawowe pary: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), a przesunięcie s → s + σ mnoży funkcję czasu przez e−σt. Aby odwrócić wymierną F(s), rozłóż ją na ułamki proste: residuum prostego bieguna p wynosi A = [(s − p)·F(s)] przy s = p. Twierdzenie o wartości końcowej: f(∞) = lims→0 s·F(s), słuszne tylko wtedy, gdy wszystkie bieguny s·F(s) leżą w lewej półpłaszczyźnie. Twierdzenie o wartości początkowej: f(0⁺) = lims→∞ s·F(s).
Rozłóż mianownik na czynniki
s² + 11s + 28 = (s + 4)(s + 7)
F(s) = A/(s + 4) + B/(s + 7)
Metoda zakrywania
Pomnóż przez (s + 4) i podstaw s = −4: składnik z B znika.
A = (4s + 22)/(s + 7) |s = −4 = 6/3 = 2
Odpowiedź
A = 2
(B = 2, więc f(t) = 2·e−4t + 2·e−7t)
36.Transformaty Laplace’a
Odpowiedź: y(0.5) = 6,169
Transformata Laplace’a: F(s) = ∫₀∞ f(t)·e−st dt zamienia pochodne na mnożenie przez s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Podstawowe pary: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), a przesunięcie s → s + σ mnoży funkcję czasu przez e−σt. Aby odwrócić wymierną F(s), rozłóż ją na ułamki proste: residuum prostego bieguna p wynosi A = [(s − p)·F(s)] przy s = p. Twierdzenie o wartości końcowej: f(∞) = lims→0 s·F(s), słuszne tylko wtedy, gdy wszystkie bieguny s·F(s) leżą w lewej półpłaszczyźnie. Twierdzenie o wartości początkowej: f(0⁺) = lims→∞ s·F(s).
Przekształcenie równania
s²Y − 9s − 12 + 7(sY − 9) + 10Y = 0
Y(s) = (9s + 75)/((s + 2)(s + 5))
Ułamki proste
A = 57/3 = 19
B = 30/(−3) = −10
Wynik
y(t) = 19·e−2t − 10·e−5t
y(0,5) = 6,169
Odpowiedź
y(0,5) = 6,169
Sprawdzenie: y(0) = A + B = 9.
37.Schematy blokowe
Odpowiedź: css = 2,857
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).
Transmitancje
C/D = G/(1 + K·G) = 3/(s + 3 + 75)
Stan ustalony
css = 0,5·3/78 = 0,01923
Odpowiedź
css = 0,01923
Większe K zmniejsza wpływ zakłócenia, ale przy regulatorze proporcjonalnym nigdy do zera.
39.Odpowiedź skokowa
Odpowiedź: σ = −5
Odpowiedź skokowa: układ pierwszego rzędu G(s) = K/(τs + 1) odpowiada na skok o wielkości U przebiegiem y(t) = K·U·(1 − e−t/τ): po jednej stałej czasowej τ osiąga 63,2 % wartości końcowej K·U, a po 4τ 98 %. Układ drugiego rzędu G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) ma bieguny s = −ζωₙ ± jωd, gdzie ωd = ωₙ√(1 − ζ²). Dla 0 < ζ < 1 odpowiedź skokowa ma przeregulowanie: przeregulowanie procentowe %OS = 100·e−ζπ/√(1−ζ²), czas szczytu tp = π/ωd, a 2-procentowy czas ustalania ts ≈ 4/(ζωₙ). Porównaj mianownik z s² + 2ζωₙs + ωₙ², aby odczytać ωₙ i ζ. Dla układu masa–sprężyna–tłumik ωₙ = √(k/m) i ζ = c/(2√(km)); dla szeregowego obwodu RLC ωₙ = 1/√(LC) i ζ = (R/2)·√(C/L).
Porównaj z postacią standardową
Mianownik musi odpowiadać s² + 2ζωₙs + ωₙ².
ωₙ² = 93 → ωₙ = 9,644 rad/s
2ζωₙ = 10 → ζ = 10/(2·9,644) = 0,5185
Wynik
σ = −ζωₙ = −10/2 = −5
Odpowiedź
σ = −5
40.Odpowiedź skokowa
Odpowiedź: K = 76,96
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² + 8s + K)
Współczynnik tłumienia z przeregulowania
ζ = −ln(0,2)/√(π² + ln²(0,2)) = 0,4559
Wzmocnienie
ωₙ = 8/(2ζ) = 8,773 rad/s
K = ωₙ² = 76,96
Odpowiedź
K = 76,96
41.Uchyb ustalony i regulacja PID
Odpowiedź: K = 240
Uchyb ustalony: w pętli ze sprzężeniem jednostkowym o transmitancji układu otwartego G(s) uchyb wynosi E(s) = R(s)/(1 + G(s)), a z twierdzenia o wartości końcowej ess = lims→0 s·E(s). Typ układu to liczba biegunów G(s) w s = 0. Stałe uchybowe to Kp = lim G(s), Kv = lim s·G(s) oraz Ka = lim s²·G(s), wszystkie dla s → 0. Skok o wartości A daje ess = A/(1 + Kp), narastanie liniowe A·t daje A/Kv, a parabola A·t²/2 daje A/Ka. Pętla typu 0 nadąża za skokiem ze skończonym uchybem, pętla typu 1 nadąża za nim dokładnie, a za narastaniem liniowym ze skończonym uchybem.
Wymagana stała uchybowa
1/(1 + Kp) = 0,04 → Kp = 24
Wzmocnienie
K/10 = 24 → K = 240
Odpowiedź
K = 240
Sprawdź też stabilność układu zamkniętego: dla tej pętli drugiego rzędu każde K > 0 jest stabilne.
42.Uchyb ustalony i regulacja PID
Odpowiedź: tp = 0,5388 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) + 25·2 = s² + 8s + 50
ωₙ i ζ
ωₙ = √50 = 7,071 rad/s
ζ = 8/(2·7,071) = 0,5657
Czas szczytu
ωd = 7,071·√(1 − 0,5657²) = 5,831 rad/s
tp = π/5,831 = 0,5388 s
Odpowiedź
tp = 0,5388 s
To niezależne materiały do ćwiczeń. ExamProblems nie jest powiązane z organizacją przeprowadzającą ten egzamin ani przez nią rekomendowane, a nazwa egzaminu służy wyłącznie do opisu treści. Egzamin FE exam prowadzi NCEES, organizacja niezwiązana z ExamProblems, która go nie popiera.