Luyện thi FE Cơ khí: gói ôn nước rút

Bảng công thức, các bước phương pháp chính cho từng chủ đề và 42 bài tập, kèm lời giải đầy đủ ở cuối.

Bài mớiQuay lại môn học

Trong cửa sổ in, hãy chọn “Lưu dưới dạng PDF” để giữ một bản trên thiết bị của bạn.

1Công thức và phương pháp

Phản lực gối của dầm

Cân bằng
ΣFx = 0, ΣFy = 0, ΣM = 0
Tải phân bố đều
hợp lực w·L đặt tại giữa đoạn tải
Mômen
M = F·d

Phương pháp

  1. Vẽ sơ đồ vật thể tự do và thay mỗi gối đỡ bằng các phản lực của nó.
  2. Thay mọi tải phân bố bằng một lực duy nhất đặt tại trọng tâm của nó.
  3. Lấy mômen quanh một gối để tìm phản lực ở gối kia.
  4. Dùng ΣFy = 0 để tìm phản lực còn lại, và kiểm tra bằng một phương trình mômen thứ hai.

Ứng suất và biến dạng dọc trục

Ứng suất pháp
σ = F/A
Biến dạng
ε = ΔL/L
Định luật Hooke
σ = E·ε
Độ giãn dài
ΔL = FL/(EA)
Biến dạng nhiệt
εT = α·ΔT

Xoắn của trục

Ứng suất tiếp
τ = T·r/J
Trục đặc
J = πd⁴/32
Trục rỗng
J = π(D⁴ − d⁴)/32
Góc xoắn
φ = TL/(GJ)
Công suất
P = T·ω = 2πnT/60

Phương pháp

  1. Tìm mômen quán tính độc cực: J = πd⁴/32 cho trục đặc.
  2. Ứng suất tiếp ở bề mặt: τ = T·r/J, với r là bán kính ngoài.
  3. Góc xoắn: φ = TL/(GJ) theo radian.
  4. Nếu cho công suất, hãy tìm mômen xoắn trước bằng P = T·ω, ω = 2πn/60.

Lực cắt và mômen uốn

Lực cắt và mômen
V = ΣF ở một phía, M = ΣF·d ở một phía
Ứng suất uốn
σ = M·y/I
Hình chữ nhật
I = bh³/12
Mômen chống uốn của mặt cắt
W = I/ymax
Dầm đơn giản, tải phân bố đều
Mmax = wL²/8

Phương pháp

  1. Tìm phản lực gối từ cân bằng của cả dầm.
  2. Cắt dầm tại mặt cắt bạn cần và giữ lại một phía.
  3. Lực cắt: tổng các lực thẳng đứng ở phía đó. Mômen uốn: tổng mômen của chúng quanh mặt cắt.
  4. Để tìm mômen lớn nhất, hãy tìm nơi lực cắt bằng không.

Vòng tròn Mohr và biến đổi ứng suất

Tâm
C = (σx + σy)/2
Bán kính
R = √(((σx − σy)/2)² + τxy²)
Ứng suất chính
σ1,2 = C ± R
Phương chính
tan 2θp = 2τxy/(σx − σy)
Mặt phẳng nghiêng
σθ = C + (σx − σy)/2·cos 2θ + τxy·sin 2θ

Phương pháp

  1. Tính tâm của vòng tròn Mohr: C = (σx + σy)/2.
  2. Tính bán kính: R = √(((σx − σy)/2)² + τxy²).
  3. Các ứng suất chính là C ± R và ứng suất tiếp lớn nhất trong mặt phẳng là R.
  4. Các góc trên vòng tròn gấp đôi các góc trên phần tử: tan 2θp = 2τxy/(σx − σy).

Độ võng của dầm

Tải tập trung ở giữa nhịp
δ = PL³/(48EI)
Tải phân bố đều, dầm đơn giản
δ = 5wL⁴/(384EI)
Công xôn, tải ở đầu
δ = PL³/(3EI)
Công xôn, tải phân bố đều
δ = wL⁴/(8EI)
Nguyên lý xếp chồng
δ = δ1 + δ2

Mất ổn định cột

Tải Euler
Pcr = π²EI/(KL)²
Hệ số chiều dài tính toán
K = 1 khớp–khớp, 0,5 ngàm–ngàm, 0,7 ngàm–khớp, 2 ngàm–tự do
Ứng suất tới hạn
σcr = Pcr/A
Độ mảnh
λ = KL/r, r = √(I/A)

Chuyển động ném

Các thành phần
vx = v0cos θ, vy = v0sin θ
Vị trí
x = vxt, y = h + vyt − ½gt²
Độ cao cực đại
H = h + vy²/(2g)
Tầm xa (mặt đất bằng phẳng)
R = v0²·sin 2θ/g

Phương pháp

  1. Phân tích vận tốc ban đầu thành các thành phần: vx = v₀·cos θ và vy = v₀·sin θ.
  2. Chuyển động ngang có tốc độ không đổi; chuyển động thẳng đứng có gia tốc −g.
  3. Tìm thời gian bay từ chuyển động thẳng đứng.
  4. Dùng thời gian đó trong chuyển động ngang để tìm tầm xa.

Công và năng lượng

Động năng
Ek = ½mv²
Thế năng
Ep = mgh
Cân bằng năng lượng
Ek1 + Ep1 + W = Ek2 + Ep2 + losses
Công của lực ma sát
W = μ·N·d
Công suất
P = W/t

Dao động cơ học

Tần số riêng
ωn = √(k/m), fn = ωn/(2π)
Từ độ võng tĩnh
ωn = √(g/δst)
Cản tới hạn
cc = 2√(km) = 2mωn
Hệ số cản
ζ = c/cc
Tần số có cản
ωd = ωn√(1 − ζ²)
Giảm lượng logarit
δ = (1/n)·ln(x0/xn), ζ = δ/√(4π² + δ²)
Lò xo song song
k = k1 + k2
Lò xo nối tiếp
1/k = 1/k1 + 1/k2

Động cơ nhiệt và máy lạnh

Hiệu suất
η = W/QH = 1 − QC/QH
Hiệu suất Carnot
ηC = 1 − TC/TH (kelvin)
COP của tủ lạnh
COPR = QC/W
COP của bơm nhiệt
COPHP = QH/W = COPR + 1
Cân bằng năng lượng
QH = W + QC

Chu trình công suất khí

Đẳng entropi, tỉ số thể tích
T2/T1 = rk−1
Đẳng entropi, tỉ số áp suất
T2/T1 = rp(k−1)/k
Hiệu suất Otto
η = 1 − 1/rk−1
Hiệu suất Diesel
η = 1 − (1/rk−1)·(rck − 1)/(k(rc − 1))
Hiệu suất Brayton
η = 1 − 1/rp(k−1)/k
Nhiệt cấp vào
Otto: cv(T3 − T2) Diesel, Brayton: cp(T3 − T2)
Công ròng của Brayton
w = cp[(T3 − T4) − (T2 − T1)]
Hiệu suất đẳng entropi
wc = wc,s/ηc, wt = ηt·wt,s
Không khí lạnh
k = 1,4, cp = 1,005, cv = 0,718 kJ/(kg·K)

Phương trình liên tục và Bernoulli

Phương trình liên tục
A1v1 = A2v2
Bernoulli
p1 + ½ρv1² + ρgz1 = p2 + ½ρv2² + ρgz2
Lưu lượng
Q = A·v
Diện tích hình tròn
A = πd²/4

Ma sát đường ống và tổn thất cột áp

Số Reynolds
Re = V·D/ν = ρ·V·D/μ
Chế độ chảy
chảy tầng dưới Re ≈ 2000, chảy rối trên Re ≈ 4000
Vận tốc trung bình
V = Q/A, A = π·D²/4
Darcy–Weisbach
hf = f·(L/D)·V²/(2g)
Tổn thất cục bộ
hm = ΣK·V²/(2g)
Giá trị K thường gặp
lối vào 0,5, lối ra 1,0, khúc cua 90° ≈ 0,3
Độ giảm áp suất
Δp = ρ·g·hf
Hai bể chứa hở
Δz = (f·L/D + ΣK)·V²/(2g)
Cửa ra tự do
H = (1 + f·L/D + ΣK)·V²/(2g)

Máy bơm và công suất bơm

Công suất thủy lực
Ph = ρ·g·Q·H = Q·Δp
Công suất trên trục
Pshaft = Ph/ηpump
Công suất điện
Pel = Pshaft/ηmotor
Cột áp của bơm
H = Δz + Δp/(ρ·g) + hL
Tổn thất của hệ thống
hL = (f·L/D + ΣK)·V²/(2g)

Dẫn nhiệt

Định luật Fourier
Q̇ = kAΔT/L
Nhiệt trở dẫn nhiệt
R = L/(kA)
Nhiệt trở đối lưu
R = 1/(hA)
Các nhiệt trở nối tiếp
Q̇ = ΔTtotal/ΣR
Hệ số U
U = 1/(ΣR·A)

Bánh răng và hệ bánh răng

Bánh răng ăn khớp
n1·N1 = n2·N2
Mômen xoắn
T2 = T1·(N2/N1)·η
Công suất
P = T·ω, ω = 2π·n/60
Hệ bánh răng ghép
tỉ số truyền tổng = tích các tỉ số truyền từng cấp

Sai số xác lập và điều khiển PID

Loại của hệ
số cực của G(s) tại s = 0
Hằng số sai số
Kp = lim G(s), Kv = lim s·G(s), Ka = lim s²·G(s) (s → 0)
Bước nhảy A
ess = A/(1 + Kp)
Tín hiệu dốc A·t
ess = A/Kv
Parabol A·t²/2
ess = A/Ka
Bộ điều khiển PID
C(s) = Kp + Ki/s + Kd·s
Mục tiêu bậc hai
s² + 2ζωns + ωn²
Độ vượt quá và thời gian
%OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd, ts ≈ 4/(ζωn)

Phương pháp

  1. Đếm các cực của vòng hở C(s)·G(s) tại s = 0: đó là loại của hệ.
  2. Tính hằng số sai số khớp với đầu vào: Kp cho bước nhảy, Kv cho tín hiệu dốc, Ka cho parabol.
  3. Áp dụng ess = A/(1 + Kp), A/Kv hay A/Ka, hoặc giải nó theo độ lợi bạn cần.
  4. Với yêu cầu về độ cản hoặc tốc độ, viết đa thức đặc trưng vòng kín và so sánh nó với s² + 2ζω_n·s + ω_n².

2Bài tập

1.Phản lực gối của dầmTrung bình
P1​ = 55 kNP2​ = 70 kNABCRA​ = ?2 m2 m2 m

Dầm có gối cố định tại A và gối con lăn tại B, cách A 4 m, và đoạn công-xôn nhô ra khỏi B tới đầu tự do C. Tìm phản lực thẳng đứng RA, lấy chiều hướng lên là dương.

AB = 4 mAC = 6 mP1 = 55 kN tại 2 mP2 = 70 kN tại 6 m
RA = kN
2.Phản lực gối của dầmKhó
w = 2 kN/mP = 10 kNADCRD​ = ?6 m2 m2 m

Một dầm ghép được ngàm vào tường tại A, có một khớp bản lề bên trong tại C (cách A 6 m) và tựa trên một gối con lăn tại D (cách A 10 m). Dầm chịu một lực tập trung P tại vị trí cách A 8 m, một tải phân bố đều w trên toàn bộ đoạn CD. Tìm phản lực RD tại gối con lăn.

AC = 6 mAD = 10 mP = 10 kN tại 8 mw = 2 kN/m từ 6 đến 10 m
RD = kN
3.Ứng suất và biến dạng dọc trụcTrung bình
Ø 16 mmĐồng thau2 mP = 100 kN

Tìm độ giãn dài δ của thanh đồng thau dưới tải trọng dọc trục.

P = 100 kNd = 16 mmL = 2 mE = 100 GPa
δ = mm
4.Ứng suất và biến dạng dọc trụcKhó
Ø 25 mmNhôm1 mØ 16 mmThép1 mP = 100 kN

Thanh bậc được ngàm vào tường và chịu tải ở đầu tự do. Tìm tổng độ giãn dài δ.

P = 100 kNđoạn 1: Ø25 mm, 1 m, nhôm (E = 70 GPa)đoạn 2: Ø16 mm, 1 m, thép (E = 200 GPa)
δ = mm
5.Xoắn của trụcTrung bình
d = ?T = 200 N·mL = 1,5 m

Một trục thép đặc dài 1,5 m chịu mô-men xoắn 200 N·m. Góc xoắn giữa hai đầu không được vượt quá 2°. Tìm đường kính nhỏ nhất d đáp ứng giới hạn này.

T = 200 N·mL = 1,5 mG = 80 GPaφallow = 2°
d = mm
6.Xoắn của trụcKhó
do​ = ?10 kW, 500 rpm

Một trục rỗng phải truyền 10 kW ở 500 rpm. Đường kính trong của nó bằng 0,75 lần đường kính ngoài, và ứng suất cắt không được vượt quá 40 MPa. Tìm đường kính ngoài nhỏ nhất do.

P = 10 kWN = 500 rpmdi/do = 0,75τallow = 40 MPa
do = mm
7.Lực cắt và mômen uốnTrung bình
15 kN15 kNABCDEFG2 m2 m2 m2 m2 m2 m

Các điểm đánh dấu từ A đến G chia dầm thành sáu phần bằng nhau. Tại điểm nào mômen uốn có độ lớn lớn nhất, dù võng hay vồng? Nhấp vào nó trong sơ đồ.

15 kN @ 4 m15 kN @ 10 m
  1. 1A (x = 0 m)
  2. 2B (x = 2 m)
  3. 3C (x = 4 m)
  4. 4D (x = 6 m)
  5. 5E (x = 8 m)
  6. 6F (x = 10 m)
  7. 7G (x = 12 m)
8.Lực cắt và mômen uốnKhó
w = 8 kN/mP = 50 kNABmặt cắt4 m4 m

Tìm mômen uốn Mx tại mặt cắt x = 4 m tính từ gối A. Coi mômen làm dầm võng xuống là dương và lực cắt hướng lên trên mặt trái là dương.

L = 8 mP = 50 kN tại 4 mw = 8 kN/m từ 0 đến 8 mx = 4 m
Mx = kN·m
9.Vòng tròn Mohr và biến đổi ứng suấtTrung bình
σx​ = 20σy​ = 60MPaτxy​ = ?

Một điểm có σₓ = 20 MPa và σᵧ = 60 MPa. Ứng suất cắt τxy (độ lớn) có thể lớn đến mức nào trước khi ứng suất chính lớn nhất đạt 70 MPa?

σx = 20 MPaσy = 60 MPaσ1 = 70 MPa
τxy = MPa
10.Vòng tròn Mohr và biến đổi ứng suấtKhó
σx​ = −80σy​ = −30τxy​ = 60MPa

Một điểm trong chi tiết chịu tải có σₓ = −80 MPa, σᵧ = −30 MPa và τₓᵧ = 60 MPa (kéo là dương). Tìm ứng suất pháp σθ trên một mặt có pháp tuyến xoay 30° ngược chiều kim đồng hồ so với trục x.

σx = −80 MPaσy = −30 MPaτxy = 60 MPaθ = 30°
σθ = MPa
11.Độ võng của dầmTrung bình
w = 5 kN/mδL = 3 m

Một dầm đơn giản chịu tải phân bố đều w trên toàn nhịp. Tìm độ võng lớn nhất δmax, tại giữa nhịp.

w = 5 kN/mL = 3 mE = 11 GPa (gỗ xẻ)I = 23 × 10⁶ mm⁴
δmax = mm
12.Độ võng của dầmKhó
w = 15 kN/mδL = 1,5 m

Một dầm đơn giản chịu tải phân bố đều w. Độ võng của nó không được vượt quá nhịp/360. Tìm mômen quán tính nhỏ nhất của mặt cắt Ireq thỏa mãn giới hạn này.

w = 15 kN/mL = 1,5 mE = 70 GPa (nhôm)
Ireq = × 10⁶ mm⁴
13.Mất ổn định cộtTrung bình
PL = 6 m20 × 80 mmmặt cắt ngang

Một cột dài L làm bằng thép được ngàm ở chân và khớp ở đỉnh. Dùng công thức Euler, tìm tải trọng oằn tới hạn Pcr của nó.

L = 6 mE = 200 GPa (thép)b × h = 20 × 80 mm
Pcr = kN
14.Mất ổn định cộtKhó
P = 95 kNL = 2,5 md = ?mặt cắt ngang

Một cột tròn đặc dài L làm bằng thép được ngàm cả hai đầu. Nó phải chịu 95 kN với hệ số an toàn 2,5 chống oằn. Tìm đường kính nhỏ nhất d.

P = 95 kNn = 2,5L = 2,5 mE = 200 GPa (thép)
d = mm
15.Chuyển động némTrung bình
v0​ = 34 m/sθ = ?R = 80 m

Một quả bóng được phóng từ mặt đất bằng phẳng với tốc độ 34 m/s và phải rơi cách 80 m. Bỏ qua sức cản không khí và dùng g = 9,81 m/s². Tìm góc nhỏ hơn trong hai góc phóng có thể θ.

v₀ = 34 m/sR = 80 mg = 9,81 m/s²
θ = °
16.Chuyển động némKhó
h = 35 mv0​ = 23 m/s50°tầm xa

Một quả bóng được ném với góc 50° so với phương ngang từ đỉnh một vách đá cao 35 m. Bỏ qua lực cản không khí. Tìm thời gian bay t cho đến khi quả bóng chạm đất.

v0 = 23 m/sθ = 50°h = 35 mg = 9,81 m/s²
t = s
17.Công và năng lượngTrung bình
0,05 kgk = 400 N/mh = ?

Một bộ phóng đồ chơi có lò xo thẳng đứng độ cứng 400 N/m. Một quả cầu 0,05 kg được đặt lên nó và lò xo bị nén 4 cm, rồi được thả ra. Bỏ qua sức cản không khí và khối lượng của lò xo, tìm độ cao h mà quả cầu bay lên, đo từ vị trí của nó khi lò xo bị nén.

m = 0,05 kgk = 400 N/mx = 4 cmg = 9,81 m/s²
h = m
18.Công và năng lượngKhó
h = 12 mR = 4 mN = ?

Một xe đẩy 0,5 kg được thả từ trạng thái nghỉ ở độ cao 12 m trên một đường ray nhẵn dẫn vào một vòng tròn thẳng đứng bán kính 4 m. Tìm phản lực pháp tuyến N mà đường ray tác dụng lên xe đẩy tại đỉnh vòng.

m = 0,5 kgh = 12 mR = 4 mg = 9,81 m/s²
N = N
19.Dao động cơ họcTrung bình
k = 12 kN/mc = 240 N·s/mm = 20 kg

Một vật 20 kg được đỡ bởi một lò xo có độ cứng 12 kN/m và một bộ cản nhớt có hệ số cản 240 N·s/m. Tìm hệ số cản ζ.

m = 20 kgk = 12 kN/mc = 240 N·s/m
ζ =
20.Dao động cơ họcKhó
k1​ = 40k2​ = 50m = 10 kg

Một vật 10 kg treo vào hai lò xo nối đầu với đầu (nối tiếp): lò xo trên có k₁ = 40 kN/m và lò xo dưới có k₂ = 50 kN/m. Tìm độ cứng tương đương keq của hệ lò xo theo kN/m.

m = 10 kgk1 = 40 kN/mk2 = 50 kN/m
keq = kN/m
21.Động cơ nhiệt và máy lạnhTrung bình
TH​ = 225 °CTL​ = 10 °CEQH​QL​W

Một động cơ nhiệt làm việc giữa nguồn nóng ở 225 °C và nguồn lạnh ở 10 °C. Tìm hiệu suất nhiệt lớn nhất có thể ηmax.

TH = 225 °CTL = 10 °C
ηmax = %
22.Động cơ nhiệt và máy lạnhKhó
TH​ = 650 °CT0​ = 30 °CTL​ = 0 °CERWQ̇H​ = 120 kWQ̇L​ = ?

Một động cơ nhiệt Carnot nhận 120 kW nhiệt từ nguồn ở 650 °C và thải nhiệt ra môi trường ở 30 °C. Toàn bộ công của nó chạy một tủ lạnh Carnot lấy nhiệt từ một tủ đông ở 0 °C và cũng thải nhiệt ra môi trường ở 30 °C. Tìm tốc độ Q̇L mà nhiệt được lấy ra khỏi tủ đông.

TH = 650 °CT0 = 30 °CTL = 0 °CQ̇H = 120 kW
Q̇L = kW
23.Chu trình công suất khíTrung bình
vP1234qin​qout​

Một chu trình Diesel lý tưởng có tỉ số nén r = 16 và tỉ số cắt nhiên liệu rc = v₃/v₂ = 2,2. Khi bắt đầu nén không khí ở 320 K. Dùng các giả thiết không khí lạnh tiêu chuẩn: không khí là khí lý tưởng có nhiệt dung riêng không đổi, k = 1,4, cp = 1,005 kJ/(kg·K) và cv = 0,718 kJ/(kg·K). Tìm nhiệt cấp vào qin trên mỗi kilogam không khí, theo kJ/kg.

r = 16rc = 2,2T1 = 320 Kk = 1,4cp = 1,005 kJ/(kg·K)
qin = kJ/kg
24.Chu trình công suất khíKhó
vP1234qin​qout​

Một nhà máy tuabin khí chạy theo chu trình Brayton với tỉ số áp suất máy nén 10. Không khí vào máy nén ở 310 K và vào tuabin ở 1300 K. Cả hai máy đều là đẳng entropi. Dùng các giả thiết không khí lạnh tiêu chuẩn: không khí là khí lý tưởng có nhiệt dung riêng không đổi, k = 1,4, cp = 1,005 kJ/(kg·K) và cv = 0,718 kJ/(kg·K). Tìm công ròng sinh ra wnet trên mỗi kilogam không khí, theo kJ/kg.

rp = 10T1 = 310 KT3 = 1300 Kk = 1,4cp = 1,005 kJ/(kg·K)
wnet = kJ/kg
25.Phương trình liên tục và BernoulliTrung bình
12d1​ = 75 mm, p1​ = 110 kPad2​ = 40 mm, p2​ = ?

Nước chảy qua một ống nằm ngang thu hẹp từ 75 mm xuống 40 mm. Ở đoạn rộng, áp suất là 110 kPa và vận tốc là 1,5 m/s. Bỏ qua tổn thất, tìm áp suất p2 ở đoạn hẹp.

d1 = 75 mmd2 = 40 mmp1 = 110 kPav1 = 1,5 m/sρ = 1000 kg/m³
p2 = kPa
26.Phương trình liên tục và BernoulliKhó
12d1​ = 200 mm, p1​ = 400 kPad2​ = 75 mm, p2​ = ?Δz = 8 m

Nước chảy lên qua một ống thu hẹp từ 200 mm xuống 75 mm trong khi dâng cao 8 m. Tại tiết diện 1, áp suất là 400 kPa và vận tốc là 1,75 m/s. Bỏ qua tổn thất, tìm áp suất p2 tại tiết diện 2.

d1 = 200 mmd2 = 75 mmp1 = 400 kPav1 = 1,75 m/sΔz = 8 mρ = 1000 kg/m³
p2 = kPa
27.Ma sát đường ống và tổn thất cột ápTrung bình
V = 2,3 m/sD = 100 mmL = 50 mf = 0,015

Nước chảy với vận tốc trung bình 2,3 m/s qua một ống dài 50 m, đường kính trong 100 mm. Dùng hệ số ma sát Darcy f = 0,015 và g = 9,81 m/s². Tìm tổn thất cột áp do ma sát hf.

V = 2,3 m/sL = 50 mD = 100 mmf = 0,015g = 9,81 m/s²
hf = m
28.Ma sát đường ống và tổn thất cột ápKhó
K = 5cửa xả tự doH = ?L = 150 m, D = 200 mm, f = 0,0153× khúc cua (K = 0,3)

Một ống dài 150 m, đường kính trong 200 mm (hệ số ma sát Darcy f = 0,015) chạy từ một bể hở và xả thành tia tự do vào không khí. Các tổn thất cục bộ là lối vào cạnh vuông (K = 0,5), 3 khúc cua (K = 0,3 mỗi khúc) và một van (K = 5). Dùng g = 9,81 m/s². Mặt bể phải cao hơn cửa ra của ống bao nhiêu H để có lưu lượng 75 L/s?

Q = 75 L/sL = 150 mD = 200 mmf = 0,015Kentrance = 0,5Kbend = 0,3 (×3)Kvalve = 5
H = m
29.Máy bơm và công suất bơmTrung bình
PQ = 80 L/sH = 52 m

Trong một thử nghiệm, một máy bơm cấp 80 L/s nước (ρ = 1000 kg/m³) ngược với cột áp 52 m trong khi trục của nó nhận 56,7 kW. Dùng g = 9,81 m/s². Tìm hiệu suất bơm η theo phần trăm.

Q = 80 L/sH = 52 mPshaft = 56,7 kWg = 9,81 m/s²
η = %
30.Máy bơm và công suất bơmKhó
PQ = 25 L/sΔz = 8 mL = 100 m, D = 150 mm, f = 0,018, ΣK = 8

Một máy bơm chuyển 25 L/s nước (ρ = 1000 kg/m³) từ một bể hở phía dưới lên một bồn phía trên có mặt nước cao hơn 8 m. Ống dài 100 m, đường kính trong 150 mm và hệ số ma sát Darcy f = 0,018; các hệ số tổn thất cục bộ cộng lại là ΣK = 8 (đã gồm lối vào và lối ra). Dùng g = 9,81 m/s². Tìm cột áp H mà bơm phải bổ sung.

Q = 25 L/sΔz = 8 mL = 100 mD = 150 mmf = 0,018ΣK = 8
H = m
31.Dẫn nhiệtTrung bình
Bông khoáng, k = 0,04250 mmGỗ, k = 0,1725 mm15 °C6 °CTi​ = ?Q̇

Một bức tường gồm 250 mm bông khoáng và 25 mm gỗ, có diện tích 15 m². Hai mặt ngoài ở 15 °C và 6 °C. Tìm nhiệt độ Ti tại mặt tiếp giáp giữa hai lớp.

L1 = 250 mm, k1 = 0,04 W/(m·K)L2 = 25 mm, k2 = 0,17 W/(m·K)A = 15 m²T1 = 15 °CT3 = 6 °C
Ti = °C
32.Dẫn nhiệtKhó
Gạch, k = 0,7250 mmPolystyrene giãn nở, k = 0,03575 mmphòng 20 °Cbên ngoài 5 °CQ̇

Một bức tường gồm 50 mm gạch và 75 mm polystyrene giãn nở ngăn cách một căn phòng ở 20 °C với không khí bên ngoài ở 5 °C. Hệ số đối lưu là 10 W/(m²·K) bên trong và 20 W/(m²·K) bên ngoài, và diện tích tường là 12 m². Tìm hệ số truyền nhiệt tổng thể U của tường.

L1 = 50 mm, k1 = 0,72 W/(m·K)L2 = 75 mm, k2 = 0,035 W/(m·K)A = 12 m²hi = 10 W/(m²·K)ho = 20 W/(m²·K)Troom = 20 °CTout = 5 °C
U = W/(m²·K)
33.Bánh răng và hệ bánh răngTrung bình
Tw​ = ?z = 1, 2900 rpm

Một động cơ 0,75 kW ở 2900 rpm dẫn động một trục vít 1 đầu mối. Bánh vít có 45 răng và bộ truyền có hiệu suất 85%. Tìm mô-men đầu ra Tw trên trục bánh vít.

z = 1Nw = 45n = 2900 rpmP = 0,75 kWη = 85%
Tw = N·m
34.Bánh răng và hệ bánh răngKhó
A = 12B = 60C = 24D = 25

Trong một hệ bánh răng ghép, bánh răng A (12 răng) trên trục vào quay 3000 rpm dẫn động bánh răng B (60 răng). Bánh răng C (24 răng) được cố định trên cùng trục với B và dẫn động bánh răng D (25 răng) trên trục ra. Tìm tốc độ ra nout.

nin = 3000 rpmA = 12B = 60C = 24D = 25
nout = rpm
35.Biến đổi LaplaceTrung bình
ReIm−2−4−6−82

Tìm f(t) cho F(s) = (4s + 29)/(s² + 6s + 18) bằng cách hoàn thành bình phương, và tính f(0,2).

F(s) = (4s + 29)/(s² + 6s + 18)t = 0,2 s
f(0.2) =
36.Biến đổi LaplaceKhó
t (s)y(t)t = 0,25y(0,25) = ?

Giải y″ + 9y′ + 18y = 0 với y(0) = −4 và y′(0) = −13 bằng biến đổi Laplace, và tìm y(0,25).

y″ + 9y′ + 18y = 0y(0) = −4y′(0) = −13t = 0,25 s
y(0.25) =
37.Sơ đồ khốiTrung bình
R(s)+Ks(s + 15)G(s)C(s)−

Đối tượng G(s) = K/(s(s + 15)) nằm trong vòng hồi tiếp âm đơn vị. Tìm độ lợi K để vòng kín tắt dần tới hạn (hai cực thực bằng nhau).

G(s) = K/(s(s + 15))hồi tiếp âm đơn vị
K =
38.Sơ đồ khốiKhó
R(s)+10bộ điều khiển+D(s)+2s + 8đối tượngC(s)−

Một bộ điều khiển tỷ lệ K = 10 điều khiển đối tượng G(s) = 2/(s + 8) trong vòng hồi tiếp âm đơn vị. Một nhiễu D đi vào đầu vào của đối tượng (được cộng vào đầu ra bộ điều khiển). Một bậc thang tín hiệu chuẩn có độ lớn 1 và một bậc thang nhiễu có độ lớn 5 được đặt cùng lúc. Tìm đầu ra xác lập css.

K = 10G(s) = 2/(s + 8)D = 5R = 1
css =
39.Đáp ứng bướcTrung bình
t (s)y(t)0,5?

Một hệ có hàm truyền G(s) = 30,5/(s² + 7s + 61) và nhận tín hiệu vào là bước nhảy đơn vị. Tìm độ vọt lố phần trăm %OS của đáp ứng bước.

G(s) = 30,5/(s² + 7s + 61)
%OS = %
40.Đáp ứng bướcKhó
t (s)y(t)1tp​ = 2 s10 %

Đáp ứng bước của một hệ bậc hai không có điểm không có độ vọt lố 10 %, và đỉnh đầu tiên xuất hiện tại t = 2 s. Tìm tần số riêng ωn.

%OS = 10 %tp = 2 s
ωn = rad/s
41.Sai số xác lập và điều khiển PIDTrung bình
r(t) = 1 + 0,5t+10(s + 2)s(s + 4)(s + 10)C(s)−

Đầu vào r(t) = 1 + 0,5t (một bước nhảy cộng một tín hiệu dốc) được đưa vào một vòng phản hồi đơn vị với G(s) = 10(s + 2)/(s(s + 4)(s + 10)). Tìm sai số xác lập ess.

G(s) = 10(s + 2)/(s(s + 4)(s + 10))r(t) = 1 + 0,5t
ess =
42.Sai số xác lập và điều khiển PIDKhó
R(s)+Kp​ + Ki​/s4s + 5C(s)−

Một bộ điều khiển PI C(s) = Kp + Ki/s điều khiển đối tượng G(s) = 4/(s + 5) trong một vòng phản hồi đơn vị. Vòng kín phải có ωₙ = 6 rad/s và ζ = 0,5. Tìm độ lợi tỉ lệ Kp.

G(s) = 4/(s + 5)ωₙ = 6 rad/sζ = 0,5
Kp =

3Dung dịch

1.Phản lực gối của dầm

Đáp số: RA = −7,5 kN

Cân bằng tĩnh: một vật ở trạng thái nghỉ có ΣFx = 0, ΣFy = 0 và ΣM = 0 đối với mọi điểm. Thay mỗi tải phân bố bằng hợp lực của nó (diện tích) đặt tại trọng tâm của biểu đồ tải, và lấy mô-men đối với một gối để khử phản lực của gối đó khỏi phương trình.

Mômen quanh A

Ngược chiều kim đồng hồ là dương.

ΣMA = 0: RB·4 = 55·2 + 70·6
RB = 530/4 = 132,5 kN

Cân bằng theo phương thẳng đứng

RA = 55 + 70 − 132,5 = −7,5 kN

Đáp số

RA = −7,5 kN

RA âm: gối cố định tại A phải giữ dầm xuống.

2.Phản lực gối của dầm

Đáp số: RD = 9 kN

Cân bằng tĩnh: một vật ở trạng thái nghỉ có ΣFx = 0, ΣFy = 0 và ΣM = 0 đối với mọi điểm. Thay mỗi tải phân bố bằng hợp lực của nó (diện tích) đặt tại trọng tâm của biểu đồ tải, và lấy mô-men đối với một gối để khử phản lực của gối đó khỏi phương trình. Một khớp bản lề bên trong không truyền được mô-men uốn, điều này cho thêm một phương trình: mô-men đối với khớp của các lực ở mỗi phía bằng không.

Phần CD

Khớp không có mô-men, nên mô-men đối với C chỉ trên CD phải cân bằng.

W = 2·4 = 8 kN tại 2 m từ C
ΣMC = 0: RD·4 = 10·2 + 8·2
RD = 9 kN
VC = 18 − 9 = 9 kN

Đáp số

RD = 9 kN
3.Ứng suất và biến dạng dọc trục

Đáp số: δ = 9,947 mm

Định luật Hooke: trong giới hạn đàn hồi σ = E·ε. Kết hợp với σ = P/A và ε = δ/L ta được δ = PL/(AE).

Ứng suất

A = π·16²/4 = 201,1 mm²
σ = P/A = 100·1000/201,1 = 497,4 MPa

Biến dạng (định luật Hooke)

E for đồng thau = 100 GPa = 100000 MPa.

ε = σ/E = 497,4/100000 = 0,004974 = 4974 µε

Độ giãn dài

δ = ε·L = 0,004974·2000 mm = 9,947 mm

Đáp số

δ = 9,947 mm
4.Ứng suất và biến dạng dọc trục

Đáp số: δ = 5,397 mm

Cấu kiện chịu tải dọc trục: mỗi đoạn lăng trụ tuân theo δ = PL/(AE). Với các cấu kiện nối tiếp, độ giãn dài cộng lại.

Cùng một lực trong cả hai đoạn

Đầu tự do chịu lực P, nên mọi mặt cắt của thanh đều chịu cùng một lực dọc trục P. Chỉ có diện tích và vật liệu thay đổi.

A1 = π·25²/4 = 490,9 mm², A2 = π·16²/4 = 201,1 mm²
σ1 = 203,7 MPa
σ2 = 497,4 MPa

Cộng độ giãn dài của các đoạn

δ1 = σ1L1/E1 = 203,7·1000/70000 = 2,91 mm
δ2 = σ2L2/E2 = 497,4·1000/200000 = 2,487 mm
δ = δ1 + δ2 = 5,397 mm

Đáp số

δ = 5,397 mm
5.Xoắn của trục

Đáp số: d = 32,34 mm

Xoắn thanh tròn: ứng suất tiếp tăng tuyến tính từ không ở tâm đến τmax = T·r/J ở bề mặt, với J = πd⁴/32 cho trục đặc và J = π(do⁴ − di⁴)/32 cho trục rỗng. Góc xoắn là φ = TL/(GJ) tính bằng radian. Công suất và mômen xoắn liên hệ với nhau qua P = T·ω.

Mô-men quán tính độc cực cần thiết

T = 200 N·m = 200000 N·mm
L = 1,5 m = 1500 mm
G = 80 GPa = 80000 MPa
φ = 2·π/180 = 0,03491 rad
J = TL/(Gφ) = 200000·1500/(80000·0,03491) = 107430 mm⁴

Đường kính

d = (32J/π)1/4 = 32,34 mm

Đáp số

d = 32,34 mm

Trong thực tế bạn sẽ làm tròn lên kích thước tiêu chuẩn kế tiếp.

6.Xoắn của trục

Đáp số: do = 32,89 mm

Xoắn thanh tròn: ứng suất tiếp tăng tuyến tính từ không ở tâm đến τmax = T·r/J ở bề mặt, với J = πd⁴/32 cho trục đặc và J = π(do⁴ − di⁴)/32 cho trục rỗng. Góc xoắn là φ = TL/(GJ) tính bằng radian. Công suất và mômen xoắn liên hệ với nhau qua P = T·ω.

Mômen xoắn

ω = 2πN/60 = 2π·500/60 = 52,36 rad/s
T = P/ω = 10000/52,36 = 191 N·m

Đường kính ngoài từ ứng suất cho phép

T = 191 N·m = 190986 N·mm
1 − k⁴ = 1 − 0,75⁴ = 0,6836
do = ∛(16T/(π·τallow·(1 − k⁴))) = ∛(16·190986/(π·40·0,6836)) = 32,89 mm
di = 0,75·32,89 = 24,67 mm

Đáp số

do = 32,89 mm

Một trục đặc sẽ cần d = 28,97 mm; trục rỗng lớn hơn một chút nhưng nhẹ hơn nhiều.

7.Lực cắt và mômen uốn

Đáp số: C (x = 4 m)

Mômen uốn lớn nhất: M đổi độ dốc ở bất cứ chỗ nào lực cắt thay đổi, nên các đỉnh của nó nằm nơi lực cắt đi qua không: dưới một tải tập trung, bên trong một tải phân bố, hoặc trên một gối của phần công-xon. Tại ngàm của dầm công-xon, mômen lớn nhất ở chân tường. So sánh độ lớn của M tại những chỗ đó, dù là mômen làm võng hay làm vồng.

Phản lực gối

Lấy mômen quanh một gối tựa, rồi dùng cân bằng lực thẳng đứng.

R1 = 12,5 kN, R2 = 17,5 kN

Mômen uốn tại mỗi điểm đánh dấu

Cộng các mômen của mọi lực ở một phía của mặt cắt. Võng là dương, vồng là âm.

MA = 0, MB = 25, MC = 50 kN·m, MD = 45, ME = 40, MF = 35, MG = 0

Đáp số

Điểm đánh dấu C: MC = 50 kN·m
8.Lực cắt và mômen uốn

Đáp số: Mx = 164 kN·m

Phương pháp mặt cắt: cắt dầm tại vị trí cần tìm nội lực và giữ lại một phía. Lực cắt V cân bằng với các lực thẳng đứng ở phía đó, và mômen uốn M cân bằng với mômen của chúng quanh mặt cắt.

Phản lực gối

Cân bằng của cả dầm.

Hợp lực của tải phân bố đều: W = w·ℓ = 8·8 = 64 kN, đặt tại x̄ = 4 m
ΣMA = 0: RB·8 = 50·4 + 8·8·4
RB = 456/8 = 57 kN
ΣFy = 0: RA = 114 − 57 = 57 kN
w = 8 kN/mP = 50 kNABRA​ = 57RB​ = 57mặt cắt4 m4 m

Cắt tại x = 4 m và giữ phần bên trái

Chỉ các lực bên trái mặt cắt xuất hiện. Phần tải phân bố đều được thay bằng một lực đặt tại giữa phần đó.

M = RA·x − (mỗi tải bên trái mặt cắt × cánh tay đòn đến mặt cắt)
Mx = 57·4 − 8·4·2 = 164 kN·m

Đáp số

Mx = 164 kN·m
9.Vòng tròn Mohr và biến đổi ứng suất

Đáp số: τxy = 22,36 MPa

Biến đổi ứng suất phẳng: trên một mặt xoay góc θ, σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ. Vòng tròn Mohr cho thấy tất cả các mặt cùng lúc: tâm của nó ở C = (σₓ + σᵧ)/2 và bán kính là R = √(((σₓ − σᵧ)/2)² + τₓᵧ²). Các ứng suất chính là C ± R, và ứng suất tiếp lớn nhất trong mặt phẳng là R.

Bán kính cần thiết

C = (20 + 60)/2 = 40 MPa
R = σ1 − C = 70 − 40 = 30 MPa
στxyCσ1​σ2​

Ứng suất tiếp

τxy = √(R² − ((σₓ − σᵧ)/2)²) = √(30² − 20²) = 22,36 MPa

Đáp số

τxy = 22,36 MPa
10.Vòng tròn Mohr và biến đổi ứng suất

Đáp số: σθ = −15,54 MPa

Biến đổi ứng suất phẳng: trên một mặt xoay góc θ, σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ. Vòng tròn Mohr cho thấy tất cả các mặt cùng lúc: tâm của nó ở C = (σₓ + σᵧ)/2 và bán kính là R = √(((σₓ − σᵧ)/2)² + τₓᵧ²). Các ứng suất chính là C ± R, và ứng suất tiếp lớn nhất trong mặt phẳng là R.

Phương trình biến đổi

σθ = −55 + −25·0,5 + 60·0,866 = −15,54 MPa

Đáp số

σθ = −15,54 MPa
11.Độ võng của dầm

Đáp số: δmax = 20,84 mm

Độ võng của dầm: độ võng có được bằng cách tích phân độ cong M/(EI) hai lần. Các trường hợp tiêu chuẩn đã được lập bảng, ví dụ PL³/(3EI) cho dầm công xôn chịu tải ở đầu tự do và 5wL⁴/(384EI) cho dầm đơn giản chịu tải phân bố đều toàn nhịp. Độ võng do các tải riêng rẽ cộng lại (xếp chồng).

Đơn vị nhất quán

Làm việc theo N và mm: P theo N, w theo N/mm (kN/m = N/mm), L theo mm, E theo N/mm² (GPa × 1000), I theo mm⁴.

w = 5 N/mm, L = 3000 mm
EI = 2,53 × 10¹¹ N·mm²

Độ võng

δmax = 5wL⁴/(384EI) = 5·5·3000⁴/(384·2,53 × 10¹¹) = 20,84 mm

Đáp số

δmax = 20,84 mm
12.Độ võng của dầm

Đáp số: Ireq = 3,39 × 10⁶ mm⁴

Độ võng của dầm: độ võng có được bằng cách tích phân độ cong M/(EI) hai lần. Các trường hợp tiêu chuẩn đã được lập bảng, ví dụ PL³/(3EI) cho dầm công xôn chịu tải ở đầu tự do và 5wL⁴/(384EI) cho dầm đơn giản chịu tải phân bố đều toàn nhịp. Độ võng do các tải riêng rẽ cộng lại (xếp chồng).

Giới hạn độ võng

δallow = L/360 = 1500/360 = 4,167 mm

Mômen quán tính mặt cắt yêu cầu

I = 5wL⁴/(384·E·δallow) = 5·15·1500⁴/(384·70000·4,167) = 3,39 × 10⁶ mm⁴
Ireq = 3,39 × 10⁶ mm⁴

Đáp số

Ireq = 3,39 × 10⁶ mm⁴
13.Mất ổn định cột

Đáp số: Pcr = 5,968 kN

Mất ổn định Euler: một cột mảnh mất ổn định tại Pcr = π²EI/(KL)², trong đó KL là chiều dài tính toán: K = 1 cho khớp–khớp, 2 cho ngàm–tự do, 0,7 cho ngàm–khớp và 0,5 cho ngàm–ngàm. Cột oằn quanh trục có I nhỏ nhất. Công thức đúng với cột mảnh, có ứng suất tới hạn thấp hơn ứng suất chảy.

Mômen quán tính nhỏ nhất của mặt cắt

Hình chữ nhật oằn quanh trục yếu, dùng kích thước nhỏ hơn lập phương.

Imin = h·b³/12 = 80·20³/12 = 53333 mm⁴

Chiều dài tính toán

KL = 0,7·6000 = 4200 mm

Tải trọng tới hạn

Pcr = π²EI/(KL)² = π²·200000·53333/4200² = 5968 N = 5,968 kN

Đáp số

Pcr = 5,968 kN
14.Mất ổn định cột

Đáp số: d = 44,24 mm

Mất ổn định Euler: một cột mảnh mất ổn định tại Pcr = π²EI/(KL)², trong đó KL là chiều dài tính toán: K = 1 cho khớp–khớp, 2 cho ngàm–tự do, 0,7 cho ngàm–khớp và 0,5 cho ngàm–ngàm. Cột oằn quanh trục có I nhỏ nhất. Công thức đúng với cột mảnh, có ứng suất tới hạn thấp hơn ứng suất chảy.

Tải trọng tới hạn yêu cầu

Pcr = n·P = 2,5·95000 = 237500 N
KL = 0,5·2500 = 1250 mm

Giải tìm đường kính

d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·237500·1250²/(π³·200000))^(1/4) = 44,24 mm

Đáp số

d = 44,24 mm

Kiểm tra: độ mảnh KL/r là 113, nên cột là mảnh và công thức Euler áp dụng được.

15.Chuyển động ném

Đáp số: θ = 21,38 °

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.

Công thức tầm xa

R = v₀²·sin 2θ/g → sin 2θ = 9,81·80/34² = 0,6789

Góc

2θ = 42,76° → θ = 21,38°

Đáp số

θ = 21,38°

Nghiệm còn lại là 90° − θ = 68,62°.

16.Chuyển động ném

Đáp số: t = 5,015 s

Chuyển động ném: khi bỏ qua lực cản không khí, x = vxt và y = h + vyt − ½gt². Hai phương độc lập với nhau, chỉ chung thời gian.

Phân tích vận tốc ban đầu

Vận tốc ngang không đổi; chuyển động thẳng đứng có gia tốc không đổi g = 9,81 m/s² hướng xuống.

vx = v0·cos θ = 23·cos 50° = 14,78 m/s
vy = v0·sin θ = 23·sin 50° = 17,62 m/s

Thời gian trong không khí

Giải y(t) = h + vy·t − ½gt² = 0 và lấy nghiệm dương.

0 = 35 + 17,62·t − 4,905·t²
t = (vy + √(vy² + 2gh))/g = (17,62 + √(310,4 + 686,7))/9,81 = 5,015 s

Đáp số

t = 5,015 s
17.Công và năng lượng

Đáp số: h = 0,6524 m

Nguyên lý công–năng: độ biến thiên động năng bằng công của mọi lực. Viết trọng lực và lò xo thành thế năng (mgh và ½kx²) biến nó thành một phương trình cân bằng: năng lượng ban đầu = năng lượng cuối + năng lượng mất do ma sát (μN·d).

Thế năng đàn hồi

x = 4 cm = 0,04 m
½kx² = ½·400·0,04² = 0,32 J

Độ cao

m·g = 0,05·9,81 = 0,4905 N
h = ½kx²/(mg) = 0,32/0,4905 = 0,6524 m

Đáp số

h = 0,6524 m
18.Công và năng lượng

Đáp số: N = 4,905 N

Nguyên lý công–năng: độ biến thiên động năng bằng công của mọi lực. Viết trọng lực và lò xo thành thế năng (mgh và ½kx²) biến nó thành một phương trình cân bằng: năng lượng ban đầu = năng lượng cuối + năng lượng mất do ma sát (μN·d). Trên quỹ đạo tròn, hợp lực hướng về tâm phải bằng mv²/r.

Vận tốc tại đỉnh

v² = 2g(h − 2R) = 2·9,81·(12 − 8) = 78,48 m²/s²

Các lực tại đỉnh

m·g = 0,5·9,81 = 4,905 N
N = mv²/R − mg = 0,5·78,48/4 − 4,905 = 4,905 N

Đáp số

N = 4,905 N

Xe đẩy vừa đủ đi hết vòng khi h = 2,5R; khi đó N = 0 tại đỉnh.

19.Dao động cơ học

Đáp số: ζ = 0,2449

Dao động tự do của hệ vật–lò xo: tần số riêng không cản là ωn = √(k/m) theo rad/s, và fn = ωn/(2π) theo Hz. Một bộ cản nhớt c cho hệ số cản ζ = c/cc, với cản tới hạn cc = 2·√(k·m) = 2·m·ωn. Hệ cản dưới tới hạn (ζ < 1) dao động ở ωd = ωn·√(1 − ζ²). Các lò xo song song cộng lại, k = k₁ + k₂; các lò xo nối tiếp cộng theo nghịch đảo, 1/k = 1/k₁ + 1/k₂.

Cản tới hạn

cc = 2·√(k·m) = 2·√(12000·20) = 979,8 N·s/m

Hệ số cản

ζ = c/cc = 240/979,8 = 0,2449

Đáp số

ζ = 0,2449

ζ < 1, nên hệ cản dưới tới hạn và dao động khi tắt dần.

20.Dao động cơ học

Đáp số: keq = 22,22 kN/m

Dao động tự do của hệ vật–lò xo: tần số riêng không cản là ωn = √(k/m) theo rad/s, và fn = ωn/(2π) theo Hz. Một bộ cản nhớt c cho hệ số cản ζ = c/cc, với cản tới hạn cc = 2·√(k·m) = 2·m·ωn. Hệ cản dưới tới hạn (ζ < 1) dao động ở ωd = ωn·√(1 − ζ²). Các lò xo song song cộng lại, k = k₁ + k₂; các lò xo nối tiếp cộng theo nghịch đảo, 1/k = 1/k₁ + 1/k₂.

Độ cứng tương đương

Các lò xo nối tiếp chịu cùng một lực và cộng độ võng, nên chúng gộp theo nghịch đảo: 1/k = 1/k1 + 1/k2.

1/keq = 1/40 + 1/50 → keq = 40·50/(40 + 50) = 22,22 kN/m

Đáp số

keq = 22,22 kN/m
21.Động cơ nhiệt và máy lạnh

Đáp số: ηmax = 43,16 %

Giới hạn của định luật hai: không chu trình nào làm việc giữa các nguồn nhiệt có nhiệt độ tuyệt đối TH và TL có thể vượt qua chu trình thuận nghịch (Carnot): ηmax = 1 − TL/TH, COPR,max = TL/(TH − TL) và COPHP,max = TH/(TH − TL). Năng lượng luôn được bảo toàn: QH = W + QL.

Nhiệt độ tuyệt đối

TH = 225 + 273,15 = 498,1 K
TL = 10 + 273,15 = 283,1 K

Hiệu suất Carnot

ηmax = 1 − TL/TH = 1 − 283,1/498,1 = 43,16 %

Đáp số

ηmax = 43,16 %
22.Động cơ nhiệt và máy lạnh

Đáp số: Q̇L = 733,8 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.

Nhiệt độ tuyệt đối

TH = 650 + 273,15 = 923,1 K
T0 = 30 + 273,15 = 303,1 K
TL = 0 + 273,15 = 273,1 K

Công của động cơ

η = 1 − 303,1/923,1 = 0,6716
Ẇ = 0,6716·120 = 80,59 kW

Nhiệt lấy đi bởi tủ lạnh

COPR = 273,1/(303,1 − 273,1) = 9,105
Q̇L = 9,105·80,59 = 733,8 kW

Đáp số

Q̇L = 733,8 kW
23.Chu trình công suất khí

Đáp số: qin = 1170 kJ/kg

Chu trình Diesel (động cơ tự cháy do nén): 1→2 nén đẳng entropi, 2→3 cấp nhiệt đẳng áp, 3→4 giãn nở đẳng entropi, 4→1 thải nhiệt đẳng tích. Với tỉ số nén r = v₁/v₂ và tỉ số cắt nhiên liệu rc = v₃/v₂: T₂ = T₁·rk−1, T₃ = rc·T₂, qin = cp·(T₃ − T₂) và η = 1 − (1/rk−1)·(rck − 1)/(k·(rc − 1)).

Các nhiệt độ

T2 = 320·160,4 = 970,1 K
T3 = 2,2·970,1 = 2134 K

Nhiệt cấp vào đẳng áp

qin = 1,005·(2134 − 970,1) = 1170 kJ/kg

Đáp số

qin = 1170 kJ/kg
24.Chu trình công suất khí

Đáp số: wnet = 339,8 kJ/kg

Chu trình Brayton (tuabin khí): 1→2 nén đẳng entropi, 2→3 cấp nhiệt đẳng áp, 3→4 giãn nở đẳng entropi trong tuabin, 4→1 thải nhiệt đẳng áp. Với tỉ số áp suất rp = P₂/P₁: T₂ = T₁·rp(k−1)/k, T₄ = T₃/rp(k−1)/k và η = 1 − 1/rp(k−1)/k. Công ròng w = cp·[(T₃ − T₄) − (T₂ − T₁)]. Với hiệu suất đẳng entropi, máy nén cần ws/ηc và tuabin cho ηt·ws.

Các nhiệt độ đẳng entropi

T2 = 310·100,2857 = 598,5 K
T4 = 1300/100,2857 = 673,3 K

Công của máy nén và tuabin

wc = 1,005·(598,5 − 310) = 290 kJ/kg
wt = 1,005·(1300 − 673,3) = 629,8 kJ/kg

Công ròng

wnet = 629,8 − 290 = 339,8 kJ/kg

Đáp số

wnet = 339,8 kJ/kg

Tỉ số công ngược: wc/wt = 46,04%.

25.Phương trình liên tục và Bernoulli

Đáp số: p2 = 97,22 kPa

Phương trình liên tục và Bernoulli: với chất lỏng không nén được, lưu lượng ở mọi tiết diện đều như nhau, Q = A₁v₁ = A₂v₂. Dọc theo một đường dòng không có tổn thất, p + ½ρv² + ρgz là hằng số, nên ở nơi chất lỏng nhanh hơn hoặc cao hơn thì áp suất thấp hơn.

Phương trình liên tục và Bernoulli

Làm việc theo Pa, rồi chia cho 1000 để có kPa.

v2 = v1·(d1/d2)² = 1,5·(75/40)² = 5,273 m/s
½ρ(v1² − v2²) = ½·1000·(1,5² − 5,273²)/1000 = −12,78 kPa
p2 = p1 + ½ρ(v1² − v2²) = 110 + (−12,78) = 97,22 kPa

Đáp số

p2 = 97,22 kPa
26.Phương trình liên tục và Bernoulli

Đáp số: p2 = 245,6 kPa

Phương trình liên tục và Bernoulli: với chất lỏng không nén được, lưu lượng ở mọi tiết diện đều như nhau, Q = A₁v₁ = A₂v₂. Dọc theo một đường dòng không có tổn thất, p + ½ρv² + ρgz là hằng số, nên ở nơi chất lỏng nhanh hơn hoặc cao hơn thì áp suất thấp hơn.

Phương trình liên tục và Bernoulli

Làm việc theo Pa, rồi chia cho 1000 để có kPa.

v2 = v1·(d1/d2)² = 1,75·(200/75)² = 12,44 m/s
½ρ(v1² − v2²) = ½·1000·(1,75² − 12,44²)/1000 = −75,9 kPa
ρgΔz = 1000·9,81·8/1000 = 78,48 kPa
p2 = p1 + ½ρ(v1² − v2²) − ρgΔz = 400 + (−75,9) − 78,48 = 245,6 kPa

Đáp số

p2 = 245,6 kPa
27.Ma sát đường ống và tổn thất cột áp

Đáp số: hf = 2,022 m

Dòng chảy trong ống: số Reynolds Re = V·D/ν = ρ·V·D/μ phân biệt chảy tầng (Re dưới khoảng 2000) với chảy rối (trên khoảng 4000). Ma sát dọc ống gây tổn thất cột áp Darcy–Weisbach hf = f·(L/D)·V²/(2g). Các phụ kiện, van, lối vào và lối ra gây thêm tổn thất cục bộ hm = ΣK·V²/(2g). Giữa hai mặt thoáng của bể hở, phương trình năng lượng cho Δz = hf + hm, đã tính cả tổn thất lối ra K = 1.

Cột áp vận tốc

V²/(2g) = 2,3²/(2·9,81) = 0,2696 m

Darcy–Weisbach

hf = f·(L/D)·V²/(2g) = 0,015·(50/0,1)·2,3²/(2·9,81) = 2,022 m

Đáp số

hf = 2,022 m
28.Ma sát đường ống và tổn thất cột áp

Đáp số: H = 5,418 m

Dòng chảy trong ống: số Reynolds Re = V·D/ν = ρ·V·D/μ phân biệt chảy tầng (Re dưới khoảng 2000) với chảy rối (trên khoảng 4000). Ma sát dọc ống gây tổn thất cột áp Darcy–Weisbach hf = f·(L/D)·V²/(2g). Các phụ kiện, van, lối vào và lối ra gây thêm tổn thất cục bộ hm = ΣK·V²/(2g). Giữa hai mặt thoáng của bể hở, phương trình năng lượng cho Δz = hf + hm, đã tính cả tổn thất lối ra K = 1.

Vận tốc

V = Q/A = 0,075/0,03142 = 2,387 m/s
V²/(2g) = 0,2905 m

Các hệ số tổn thất

f·L/D = 0,015·150/0,2 = 11,25
ΣK = 0,5 + 3·0,3 + 5 = 6,4

Phương trình năng lượng

Tia ra vẫn mang cột áp vận tốc của nó; đó không phải là tổn thất nhưng phải được cung cấp.

H = (1 + 11,25 + 6,4)·0,2905 = 5,418 m

Đáp số

H = 5,418 m
29.Máy bơm và công suất bơm

Đáp số: η = 71,97 %

Công suất bơm: công suất truyền cho nước là công suất thủy lực Ph = ρ·g·Q·H = Q·Δp. Trục phải cung cấp nhiều hơn vì bơm có tổn hao: Pshaft = Ph/ηpump, và động cơ tiêu thụ Pel = Pshaft/ηmotor. Cột áp bơm H suy ra từ phương trình năng lượng: chiều cao nâng tĩnh + cột áp áp suất + mọi tổn thất ma sát và cục bộ.

Công suất thủy lực

Ph = ρ·g·Q·H = 1000·9,81·0,08·52 = 40,81 kW

Hiệu suất

η = Ph/Pshaft = 40,81/56,7 = 71,97%

Đáp số

η = 71,97%
30.Máy bơm và công suất bơm

Đáp số: H = 10,04 m

Công suất bơm: công suất truyền cho nước là công suất thủy lực Ph = ρ·g·Q·H = Q·Δp. Trục phải cung cấp nhiều hơn vì bơm có tổn hao: Pshaft = Ph/ηpump, và động cơ tiêu thụ Pel = Pshaft/ηmotor. Cột áp bơm H suy ra từ phương trình năng lượng: chiều cao nâng tĩnh + cột áp áp suất + mọi tổn thất ma sát và cục bộ.

Vận tốc

V = Q/A = 0,025/0,01767 = 1,415 m/s
V²/(2g) = 0,102 m

Các tổn thất cột áp

hL = (f·L/D + ΣK)·V²/(2g) = (12 + 8)·0,102 = 2,04 m

Cột áp của bơm

Cả hai mặt đều hở và thực tế đứng yên, nên các số hạng áp suất và vận tốc triệt tiêu.

H = Δz + hL = 8 + 2,04 = 10,04 m

Đáp số

H = 10,04 m
31.Dẫn nhiệt

Đáp số: Ti = 6,207 °C

Nhiệt trở: dòng nhiệt ổn định qua các lớp nối tiếp hoạt động giống như dòng điện qua các điện trở nối tiếp: Q̇ = ΔT/ΣR, với R = L/(kA) cho lớp dẫn nhiệt và R = 1/(hA) cho bề mặt đối lưu. Độ giảm nhiệt độ qua mỗi lớp tỉ lệ với nhiệt trở của nó.

Nhiệt trở của các lớp

R1 = L1/(k1A) = (250/1000)/(0,04·15) = 0,4167 K/W
R2 = L2/(k2A) = (25/1000)/(0,17·15) = 0,009804 K/W

Dòng nhiệt

Q̇ = (T1 − T3)/(R1 + R2) = 9/0,4265 = 21,1 W

Nhiệt độ tại mặt tiếp giáp

Nhiệt độ giảm qua lớp thứ nhất một lượng Q̇·R₁.

Ti = T1 − Q̇·R1 = 15 − 21,1·0,4167 = 6,207 °C

Đáp số

Ti = 6,207 °C
32.Dẫn nhiệt

Đáp số: U = 0,4233 W/(m²·K)

Nhiệt trở: dòng nhiệt ổn định qua các lớp nối tiếp hoạt động giống như dòng điện qua các điện trở nối tiếp: Q̇ = ΔT/ΣR, với R = L/(kA) cho lớp dẫn nhiệt và R = 1/(hA) cho bề mặt đối lưu. Độ giảm nhiệt độ qua mỗi lớp tỉ lệ với nhiệt trở của nó.

Các nhiệt trở nối tiếp

Ri = 1/(hiA) = 1/(10·12) = 0,008333 K/W
R1 = L1/(k1A) = (50/1000)/(0,72·12) = 0,005787 K/W
R2 = L2/(k2A) = (75/1000)/(0,035·12) = 0,1786 K/W
Ro = 1/(hoA) = 1/(20·12) = 0,004167 K/W
ΣR = 0,1969 K/W

Hệ số tổng thể

U = 1/(A·ΣR) = 1/(12·0,1969) = 0,4233 W/(m²·K)

Đáp số

U = 0,4233 W/(m²·K)
33.Bánh răng và hệ bánh răng

Đáp số: Tw = 94,46 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.

Tốc độ bánh vít

nw = n·z/Nw = 2900·1/45 = 64,44 rpm

Mômen xoắn ra

ωw = 2π·64,44/60 = 6,749 rad/s
Tw = η·P/ωw = 0,85·750/6,749 = 94,46 N·m

Đáp số

Tw = 94,46 N·m
34.Bánh răng và hệ bánh răng

Đáp số: nout = 576 rpm

Hệ bánh răng: các bánh răng ăn khớp có cùng tốc độ tại răng, nên n₁·N₁ = n₂·N₂: bánh răng có ít răng hơn quay nhanh hơn. Các bánh răng trên cùng một trục quay cùng tốc độ. Công suất được truyền đi (trừ tổn thất), nên mômen xoắn tăng khi tốc độ giảm: T₂ = T₁·(N₂/N₁)·η, và P = T·ω với ω = 2π·n/60.

Cấp thứ nhất

nB = 3000·12/60 = 600 rpm = nC

Cấp thứ hai

nout = 600·24/25 = 576 rpm

Đáp số

nout = 576 rpm

tỉ số truyền tổng = 5,208 : 1

35.Biến đổi Laplace

Đáp số: f(0.2) = 3,568

Biến đổi Laplace: F(s) = ∫₀∞ f(t)·e−st dt biến đạo hàm thành phép nhân với s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Các cặp chuẩn: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), và phép dịch s → s + σ nhân hàm theo thời gian với e−σt. Để biến đổi ngược một F(s) hữu tỉ, hãy tách thành các phân thức riêng: thặng dư của cực đơn p là A = [(s − p)·F(s)] tại s = p. Định lý giá trị cuối: f(∞) = lims→0 s·F(s), chỉ đúng khi mọi cực của s·F(s) nằm trong nửa mặt phẳng trái. Định lý giá trị đầu: f(0⁺) = lims→∞ s·F(s).

Hoàn thành bình phương

s² + 6s + 18 = (s + 3)² + 9 = (s + 3)² + 3²

Ghép tử số

F(s) = 4·(s + 3)/((s + 3)² + 3²) + (17/3)·3/((s + 3)² + 3²)

Biến đổi ngược

f(t) = e−3t·(4·cos 3t + 5,667·sin 3t)
f(0,2) = e−3·0,2·(4·cos(3·0,2) + 5,667·sin(3·0,2)) = 3,568

Đáp số

f(0,2) = 3,568
36.Biến đổi Laplace

Đáp số: y(0.25) = −3,966

Biến đổi Laplace: F(s) = ∫₀∞ f(t)·e−st dt biến đạo hàm thành phép nhân với s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Các cặp chuẩn: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), và phép dịch s → s + σ nhân hàm theo thời gian với e−σt. Để biến đổi ngược một F(s) hữu tỉ, hãy tách thành các phân thức riêng: thặng dư của cực đơn p là A = [(s − p)·F(s)] tại s = p. Định lý giá trị cuối: f(∞) = lims→0 s·F(s), chỉ đúng khi mọi cực của s·F(s) nằm trong nửa mặt phẳng trái. Định lý giá trị đầu: f(0⁺) = lims→∞ s·F(s).

Biến đổi phương trình

s²Y − (−4)s − (−13) + 9(sY − (−4)) + 18Y = 0
Y(s) = (−4s − 49)/((s + 3)(s + 6))

Phân thức riêng

A = −37/3 = −12,33
B = −25/(−3) = 8,333

Kết quả

y(t) = −12,33·e−3t + 8,333·e−6t
y(0,25) = −3,966

Đáp số

y(0,25) = −3,966

Kiểm tra: y(0) = A + B = −4.

37.Sơ đồ khối

Đáp số: K = 56,25

Đại số sơ đồ khối: các khối nối tiếp thì nhân (G₁·G₂), các khối song song thì cộng (G₁ + G₂), và một vòng có đường thuận G và đường hồi tiếp H rút gọn thành T = G/(1 + G·H) với hồi tiếp âm và G/(1 − G·H) với hồi tiếp dương. Rút gọn các vòng bên trong trước, rồi làm dần ra ngoài. Các cực vòng kín là nghiệm của phương trình đặc trưng 1 + G(s)H(s) = 0. Với vòng ổn định, đầu vào bậc thang có độ lớn A cho đầu ra xác lập A·T(0); với hồi tiếp đơn vị, sai số xác lập đối với bậc thang là A/(1 + Kp) với Kp = G(0).

Phương trình đặc trưng

1 + K/(s(s + 15)) = 0 → s² + 15s + K = 0

Nghiệm kép

Tam thức s² + βs + γ có nghiệm kép khi β² = 4γ.

15² = 4K → K = 225/4 = 56,25

Đáp số

K = 56,25
Khi đó cả hai cực vòng kín đều tại s = −7,5.
38.Sơ đồ khối

Đáp số: css = 1,071

Đại số sơ đồ khối: các khối nối tiếp thì nhân (G₁·G₂), các khối song song thì cộng (G₁ + G₂), và một vòng có đường thuận G và đường hồi tiếp H rút gọn thành T = G/(1 + G·H) với hồi tiếp âm và G/(1 − G·H) với hồi tiếp dương. Rút gọn các vòng bên trong trước, rồi làm dần ra ngoài. Các cực vòng kín là nghiệm của phương trình đặc trưng 1 + G(s)H(s) = 0. Với vòng ổn định, đầu vào bậc thang có độ lớn A cho đầu ra xác lập A·T(0); với hồi tiếp đơn vị, sai số xác lập đối với bậc thang là A/(1 + Kp) với Kp = G(0).

Hàm truyền

C/D = G/(1 + K·G) = 2/(s + 8 + 20)
C/R = K·G/(1 + K·G) = 20/(s + 8 + 20)

Trạng thái ổn định

css = 1·20/28 + 5·2/28 = 1,071

Đáp số

css = 1,071

K lớn hơn làm giảm ảnh hưởng của nhiễu, nhưng với bộ điều khiển tỷ lệ thì không bao giờ về không.

39.Đáp ứng bước

Đáp số: %OS = 20,7 %

Đáp ứng bước: một hệ bậc nhất G(s) = K/(τs + 1) đáp ứng với một bước nhảy có biên độ U là y(t) = K·U·(1 − e−t/τ): nó đạt 63,2 % giá trị cuối K·U sau một hằng số thời gian τ và 98 % sau 4τ. Một hệ bậc hai G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) có các cực s = −ζωₙ ± jωd với ωd = ωₙ√(1 − ζ²). Với 0 < ζ < 1, đáp ứng bước có độ vọt lố: độ vọt lố phần trăm %OS = 100·e−ζπ/√(1−ζ²), thời gian đỉnh tp = π/ωd, và thời gian xác lập 2 % ts ≈ 4/(ζωₙ). So sánh mẫu số với s² + 2ζωₙs + ωₙ² để suy ra ωₙ và ζ. Với hệ khối lượng–lò xo–giảm chấn ωₙ = √(k/m) và ζ = c/(2√(km)); với mạch RLC nối tiếp ωₙ = 1/√(LC) và ζ = (R/2)·√(C/L).

So sánh với dạng chuẩn

Mẫu số phải khớp với s² + 2ζωₙs + ωₙ².

ωₙ² = 61 → ωₙ = 7,81 rad/s
2ζωₙ = 7 → ζ = 7/(2·7,81) = 0,4481

Kết quả

%OS = 100·e−ζπ/√(1−ζ²) = 100·e−1,575 = 20,7 %

Đáp số

%OS = 20,7 %
40.Đáp ứng bước

Đáp số: ωn = 1,948 rad/s

Đáp ứng bước: một hệ bậc nhất G(s) = K/(τs + 1) đáp ứng với một bước nhảy có biên độ U là y(t) = K·U·(1 − e−t/τ): nó đạt 63,2 % giá trị cuối K·U sau một hằng số thời gian τ và 98 % sau 4τ. Một hệ bậc hai G(s) = K·ωₙ²/(s² + 2ζωₙs + ωₙ²) có các cực s = −ζωₙ ± jωd với ωd = ωₙ√(1 − ζ²). Với 0 < ζ < 1, đáp ứng bước có độ vọt lố: độ vọt lố phần trăm %OS = 100·e−ζπ/√(1−ζ²), thời gian đỉnh tp = π/ωd, và thời gian xác lập 2 % ts ≈ 4/(ζωₙ). So sánh mẫu số với s² + 2ζωₙs + ωₙ² để suy ra ωₙ và ζ. Với hệ khối lượng–lò xo–giảm chấn ωₙ = √(k/m) và ζ = c/(2√(km)); với mạch RLC nối tiếp ωₙ = 1/√(LC) và ζ = (R/2)·√(C/L).

Hệ số cản

ζ = −ln(0,1)/√(π² + ln²(0,1)) = 2,303/√(9,87 + 5,302) = 0,5912

Tần số riêng

ωd = π/tp = π/2 = 1,571 rad/s
ωₙ = ωd/√(1 − ζ²) = 1,571/0,8066 = 1,948 rad/s

Đáp số

ωn = 1,948 rad/s
41.Sai số xác lập và điều khiển PID

Đáp số: ess = 1

Sai số xác lập: trong một vòng phản hồi đơn vị với hàm truyền vòng hở G(s), sai số là E(s) = R(s)/(1 + G(s)), và theo định lý giá trị cuối ess = lims→0 s·E(s). Loại của hệ là số cực của G(s) tại s = 0. Các hằng số sai số là Kp = lim G(s), Kv = lim s·G(s) và Ka = lim s²·G(s), đều với s → 0. Một bước nhảy có độ lớn A cho ess = A/(1 + Kp), một tín hiệu dốc A·t cho A/Kv và một parabol A·t²/2 cho A/Ka. Một vòng loại 0 bám theo bước nhảy với sai số hữu hạn, một vòng loại 1 bám theo nó chính xác và bám theo tín hiệu dốc với sai số hữu hạn.

Hằng số vận tốc

Kv = lim s·G(s) = 0,5

Phần bước nhảy

Kp = ∞ → estep = 0

Sai số tín hiệu dốc

ess = 0,5/0,5 = 1

Đáp số

ess = 1
42.Sai số xác lập và điều khiển PID

Đáp số: Kp = 0,25

Điều khiển P, PI và PID: bộ điều khiển C(s) = Kp + Ki/s + Kd·s tác động lên sai số. Một độ lợi tỉ lệ làm vòng nhanh hơn và giảm sai số của đối tượng loại 0, nhưng không bao giờ loại bỏ nó. Thành phần tích phân thêm một cực tại s = 0: loại của hệ tăng thêm một, và một bước nhảy được bám theo không có sai số. Thành phần vi phân thêm sự cản. Để chọn các độ lợi, viết đa thức đặc trưng vòng kín và so sánh nó với s² + 2ζωns + ωn². Khi đó %OS = 100·e−ζπ/√(1−ζ²), tp = π/ωd với ωd = ωn√(1 − ζ²), và ts ≈ 4/(ζωn).

Đa thức đặc trưng

1 + C(s)G(s) = 0, nhân với s(s + 5).

s² + (5 + 4Kp)s + 4Ki

Đa thức mục tiêu

s² + 2·0,5·6s + 6² = s² + 6s + 36

Cân bằng các hệ số

5 + 4Kp = 6 → Kp = 0,25

Đáp số

Kp = 0,25

Thành phần tích phân làm vòng thành loại 1, nên một bước nhảy được bám theo mà không có sai số xác lập.

Đây là tài liệu luyện tập độc lập. ExamProblems không liên kết với hay được xác nhận bởi tổ chức tổ chức kỳ thi này, và tên kỳ thi chỉ được dùng để mô tả nội dung. Kỳ thi FE exam do NCEES tổ chức; tổ chức này không liên kết với ExamProblems và không bảo trợ nó.

ExamProblems · examproblems.com/courses/fe-exam-mechanical