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
w0​ = 15 kN/mP = 25 kNABRA​ = ?1 m9 m

Dầm đơn giản (gối cố định tại A, gối con lăn tại B) chịu tải thay đổi tuyến tính tăng từ không tại A đến 15 kN/m tại B. Một lực tập trung P cũng tác dụng tại vị trí cách A 1 m. Tìm phản lực thẳng đứng RA tại gối.

L = 10 mw0 = 15 kN/mP = 25 kN tại 1 m
RA = kN
2.Phản lực gối của dầmKhó
P = 55 kN, 50°Q = 20 kNABRA​ = ?3 m4 m3 m

Dầm đơn giản có gối cố định tại A và gối con lăn tại B. Một lực P tác dụng tại vị trí cách A 3 m, nghiêng 50° so với phương ngang, và một lực thẳng đứng Q tác dụng tại vị trí cách A 7 m. Tìm độ lớn của tổng phản lực RA tại gối cố định A.

L = 10 mP = 55 kN tại 3 m, θ = 50°Q = 20 kN tại 7 m
RA = kN
3.Ứng suất và biến dạng dọc trụcTrung bình
P = ?12600 mm300 mmØ25Ø30δ = 0,4 mm

Một thanh gồm một đoạn thép (đường kính 25 mm, chiều dài 600 mm, E = 200 GPa) nối đầu với đầu với một đoạn làm bằng nhôm (đường kính 30 mm, chiều dài 300 mm, E = 70 GPa). Nó được ngàm một đầu và bị kéo ở đầu kia. Tìm tải dọc trục P làm nó giãn tổng cộng 0,4 mm.

đoạn 1: thép, Ø25 mm, 600 mm, E = 200 GPađoạn 2: nhôm, Ø30 mm, 300 mm, E = 70 GPaδ = 0,4 mm
P = kN
4.Ứng suất và biến dạng dọc trụcKhó
Ø 40 mmThép2 mØ 12 mmNhôm0,5 mP = 50 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 = 50 kNđoạn 1: Ø40 mm, 2 m, thép (E = 200 GPa)đoạn 2: Ø12 mm, 0,5 m, nhôm (E = 70 GPa)
δ = mm
5.Xoắn của trụcTrung bình
T = 1000 N·md = 40 mmL = 1 mφ = ?

Một trục đặc đồng thau được ngàm vào tường và bị xoắn bởi mômen xoắn T ở đầu tự do. Tìm góc xoắn φ của đầu tự do, theo độ.

T = 1000 N·md = 40 mmL = 1 mG = 39 GPa
φ = °
6.Xoắn của trụcKhó
d = 60 mmT0​ = 1000 N·mACBa = 0,5 mb = 2,5 m

Một trục đặc đồng thau đường kính 60 mm được ngàm vào các tường cứng ở cả hai đầu A và B. Một mô-men xoắn 1000 N·m đặt tại C, cách A 0,5 m và cách B 2,5 m. Tìm góc quay φC của tiết diện C, theo độ.

d = 60 mmT0 = 1000 N·ma = 0,5 mb = 2,5 mG = 39 GPa
φC = °
7.Lực cắt và mômen uốnTrung bình
P1​ = 35 kNP2​ = 10 kNAB1,5 m0,5 m2 m

Dầm đơn giản chịu hai lực tập trung. Tìm mô-men uốn lớn nhất Mmax trong dầm.

L = 4 mP1 = 35 kN tại 1,5 mP2 = 10 kN tại 2 m
Mmax = kN·m
8.Lực cắt và mômen uốnKhó
w = 10 kN/mP = 30 kNABmặt cắt7,5 m1,5 m1 m

Tìm mômen uốn Mx tại mặt cắt x = 7,5 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 = 10 mP = 30 kN tại 9 mw = 10 kN/m từ 0 đến 10 mx = 7,5 m
Mx = kN·m
9.Vòng tròn Mohr và biến đổi ứng suấtTrung bình
σx​ = 50σy​ = 50τxy​ = −50MPa

Một điểm trong chi tiết chịu tải có σₓ = 50 MPa, σᵧ = 50 MPa và τₓᵧ = −50 MPa (kéo là dương). Tìm ứng suất tiếp lớn nhất trong mặt phẳng τmax.

σx = 50 MPaσy = 50 MPaτxy = −50 MPa
τmax = MPa
10.Vòng tròn Mohr và biến đổi ứng suấtKhó
σx​ = 120σy​ = 120τxy​ = 20MPa

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

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

Một dầm công xôn chịu tải phân bố đều w trên toàn chiều dài. Tìm độ võng δ tại đầu tự do.

w = 5 kN/mL = 4 mE = 200 GPa (thép)I = 120 × 10⁶ mm⁴
δ = mm
12.Độ võng của dầmKhó
w = 15 kN/mP = 23 kNδL = 2,5 m

Một dầm đơn giản chịu tải phân bố đều w và một tải tập trung P ở giữa nhịp. Dùng nguyên lý xếp chồng, tìm độ võng giữa nhịp δmax.

P = 23 kNw = 15 kN/mL = 2,5 mE = 11 GPa (gỗ xẻ)I = 57 × 10⁶ mm⁴
δmax = mm
13.Mất ổn định cộtTrung bình
PL = 2 md = 50 mmmặt cắt ngang

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

L = 2 mE = 70 GPa (nhôm)d = 50 mm
Pcr = kN
14.Mất ổn định cộtKhó
PL = 2 md = 40 mmmặt cắt ngang

Một cột dài L làm bằng thép được ngàm cả hai đầu. Tìm độ mảnh của nó và ứng suất oằn tới hạn σcr theo công thức Euler.

L = 2 mE = 200 GPa (thép)d = 40 mm
σcr = MPa
15.Chuyển động némTrung bình
v0​ = 14 m/s60°tầm xaH = ?

Một quả bóng được ném với góc 60° so với phương ngang từ mặt đất bằng phẳng. Bỏ qua lực cản không khí. Tìm độ cao cực đại H mà quả bóng đạt được.

v0 = 14 m/sθ = 60°g = 9,81 m/s²
H = m
16.Chuyển động némKhó
h = 15 mv0​ = 35 m/s45°x = ?

Một quả bóng được ném với góc 45° so với phương ngang từ đỉnh một vách đá cao 15 m. Bỏ qua lực cản không khí. Tìm quãng đường ngang x mà quả bóng đi được trước khi chạm đất.

v0 = 35 m/sθ = 45°h = 15 mg = 9,81 m/s²
x = m
17.Công và năng lượngTrung bình
25°d = 12 mμ = 0,15

Một vật 5 kg trượt 12 m xuống một mặt dốc nghiêng 25° so với phương ngang, với hệ số ma sát trượt là 0,15. Tìm công Wf chống lại ma sát.

m = 5 kgθ = 25°μ = 0,15d = 12 mg = 9,81 m/s²
Wf = J
18.Công và năng lượngKhó
h = 3 mL = 2 mμ = 0,25k = 2000 N/m

Một vật 5 kg được thả từ trạng thái nghỉ ở độ cao 3 m trên một mặt dốc nhẵn. Ở chân dốc, nó đi qua một đoạn ngang nhám dài 2 m (μ = 0,25) rồi chạy vào một lò xo (k = 2000 N/m) trên mặt nhẵn. Tìm độ nén cực đại x của lò xo.

m = 5 kgh = 3 mL = 2 mμ = 0,25k = 2000 N/mg = 9,81 m/s²
x = cm
19.Dao động cơ họcTrung bình
k = 150 kN/mc = 1370 N·s/mm = 50 kg

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

m = 50 kgk = 150 kN/mc = 1370 N·s/m
ζ =
20.Dao động cơ họcKhó
k1​ = 25k2​ = 25k3​ = 80m = 100 kg

Một vật 100 kg treo vào lò xo k₃ = 80 kN/m. Đầu trên của k₃ được gắn vào một thanh cứng treo bằng hai lò xo đặt cạnh nhau, k₁ = 25 kN/m và k₂ = 25 kN/m; thanh giữ nằm ngang. Bỏ qua lực cản và khối lượng của các lò xo và thanh. Tìm tần số riêng fn của dao động thẳng đứng theo Hz.

m = 100 kgk1 = 25 kN/mk2 = 25 kN/mk3 = 80 kN/m
fn = Hz
21.Động cơ nhiệt và máy lạnhTrung bình
TH​ = 850 °CTL​ = 15 °CEQH​QL​η = 41%

Một động cơ nhiệt nhận nhiệt từ nguồn ở 850 °C và thải nhiệt ra môi trường ở 15 °C. Hiệu suất nhiệt thực tế của nó là 41%. Tìm hiệu suất theo định luật hai ηII, là tỉ số giữa hiệu suất thực tế và hiệu suất Carnot.

TH = 850 °CTL = 15 °Cη = 41%
ηII = %
22.Động cơ nhiệt và máy lạnhKhó
TH​ = 20 °CTL​ = 2 °CHPQ̇H​ = 8 kWQ̇L​Ẇ = ?

Một bơm nhiệt giữ ngôi nhà ở 20 °C trong khi bên ngoài là 2 °C, cung cấp 8 kW nhiệt cho ngôi nhà. Tìm công suất đầu vào tối thiểu Ẇmin.

TH = 20 °CTL = 2 °CQ̇H = 8 kW
Ẇmin = 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 = 22 và tỉ số cắt nhiên liệu rc = v₃/v₂ = 2. Khi bắt đầu nén không khí ở 305 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 độ T3 ở cuối quá trình cấp nhiệt, theo K.

r = 22rc = 2T1 = 305 Kk = 1,4
T3 = K
24.Chu trình công suất khíKhó
vP1234qin​qout​

Một chu trình Brayton lý tưởng hoạt động với tỉ số áp suất máy nén 10. 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 hiệu suất nhiệt η của nó theo phần trăm.

rp = 10k = 1,4
η = %
25.Phương trình liên tục và BernoulliTrung bình
12d1​ = 200 mm, p1​ = 170 kPad2​ = 100 mm, p2​ = ?

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

d1 = 200 mmd2 = 100 mmp1 = 170 kPav1 = 1 m/sρ = 1000 kg/m³
p2 = kPa
26.Phương trình liên tục và BernoulliKhó
C213 m4,5 md = 40 mmQ = ?

Một ống mềm có đường kính trong 40 mm hút nước xi-phông từ một bể hở. Điểm cao nhất C của ống cao hơn mặt nước 3 m, và cửa ra (thông với không khí) thấp hơn mặt nước 4,5 m. Bỏ qua tổn thất và lấy g = 9,81 m/s². Tìm lưu lượng Q.

d = 40 mmzC = +3 mz2 = −4,5 mg = 9,81 m/s²
Q = L/s
27.Ma sát đường ống và tổn thất cột ápTrung bình
V = 2,4 m/sD = 50 mmL = 200 mf = 0,03

Nước (ρ = 1000 kg/m³) chảy với 2,4 m/s qua một ống nằm ngang dài 200 m, đường kính trong 50 mm. Dùng hệ số ma sát Darcy f = 0,03 và g = 9,81 m/s². Tìm độ giảm áp suất Δp trên ống theo kPa.

V = 2,4 m/sL = 200 mD = 50 mmf = 0,03ρ = 1000 kg/m³
Δp = kPa
28.Ma sát đường ống và tổn thất cột ápKhó
K = 10Δz = ?L = 150 m, D = 300 mm, f = 0,0253× khúc cua (K = 0,3)

Hai bể chứa hở được nối bằng một ống dài 150 m, đường kính trong 300 mm và hệ số ma sát Darcy f = 0,025. 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), một van (K = 10) và lối ra vào bể thấp hơn (K = 1,0). Dùng g = 9,81 m/s². Chênh lệch mực nước Δz giữa các bể là bao nhiêu để có lưu lượng 100 L/s?

Q = 100 L/sL = 150 mD = 300 mmf = 0,025Kentrance = 0,5Kbend = 0,3 (×3)Kvalve = 10Kexit = 1,0
Δz = m
29.Máy bơm và công suất bơmTrung bình
PQ = 15 L/sH = 48 m

Một máy bơm cấp 15 L/s nước (ρ = 1000 kg/m³) ngược với cột áp 48 m. Hiệu suất bơm là 82% và động cơ điện dẫn động nó có hiệu suất 88%. Dùng g = 9,81 m/s². Tìm công suất điện Pel mà động cơ tiêu thụ, theo kW.

Q = 15 L/sH = 48 mηpump = 82%ηmotor = 88%g = 9,81 m/s²
Pel = kW
30.Máy bơm và công suất bơmKhó
PQ = 40 L/sΔz = 5 mL = 600 m, D = 150 mm, f = 0,025, ΣK = 8

Một máy bơm chuyển 40 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 5 m. Ống dài 600 m, đường kính trong 150 mm và hệ số ma sát Darcy f = 0,025; 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 = 40 L/sΔz = 5 mL = 600 mD = 150 mmf = 0,025ΣK = 8
H = m
31.Dẫn nhiệtTrung bình
Kính, k = 0,96100 mmBông khoáng, k = 0,0475 mm16 °C-7 °CTi​ = ?Q̇

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

L1 = 100 mm, k1 = 0,96 W/(m·K)L2 = 75 mm, k2 = 0,04 W/(m·K)A = 15 m²T1 = 16 °CT3 = -7 °C
Ti = °C
32.Dẫn nhiệtKhó
Polystyrene giãn nở, k = 0,035200 mmKính, k = 0,96100 mmphòng 20 °Cbên ngoài -5 °CQ̇ = ?

Một bức tường gồm 200 mm polystyrene giãn nở và 100 mm kính 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à 7,7 W/(m²·K) bên trong và 30 W/(m²·K) bên ngoài, và diện tích tường là 4 m². Tìm lượng nhiệt thất thoát Q̇ qua tường.

L1 = 200 mm, k1 = 0,035 W/(m·K)L2 = 100 mm, k2 = 0,96 W/(m·K)A = 4 m²hi = 7,7 W/(m²·K)ho = 30 W/(m²·K)Troom = 20 °CTout = -5 °C
Q̇ = W
33.Bánh răng và hệ bánh răngTrung bình
N1​ = 16N2​ = 32

Một động cơ 3 kW chạy ở 1440 rpm dẫn động một bánh răng nhỏ có 16 răng, ăn khớp với một bánh răng có 32 răng. Cặp bánh răng có hiệu suất 96%. Tìm mômen xoắn T1 trên trục động cơ.

P = 3 kWn1 = 1440 rpmN1 = 16N2 = 32η = 96%
T1 = N·m
34.Bánh răng và hệ bánh răngKhó
d1​ = 80 mmd2​ = 315 mm2880 rpmC = 500 mmθ = ?

Một dây đai hở nối puli 80 mm và puli 315 mm có tâm cách nhau 500 mm. Tìm góc ôm θ trên puli nhỏ hơn, theo độ.

d1 = 80 mmd2 = 315 mmC = 500 mm
θ = °
35.Biến đổi LaplaceTrung bình
ReIm−5−105

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

F(s) = (4s + 40)/(s² + 6s + 34)t = 0,6 s
f(0.6) =
36.Biến đổi LaplaceKhó
ReIm−224

Một hệ thống được mô tả bởi phương trình vi phân y″ + 2y′ + 17y = 4u′ + 8u, với đầu vào u(t) và đầu ra y(t). Tìm hàm truyền G(s) = Y(s)/U(s) của nó (điều kiện ban đầu bằng không) và tính đáp ứng xung g(t) tại t = 0,3 s. Tìm g(0,3).

y″ + 2y′ + 17y = 4u′ + 8ut = 0,3 s
g(0.3) =
37.Sơ đồ khốiTrung bình
R(s)+11K+8s + 3G(s)C(s)1H₁−−

Đối tượng G(s) = 8/(s + 3) có một vòng hồi tiếp âm trong (phụ) với độ lợi H₁ = 1. Độ lợi bộ điều khiển K = 11 và một vòng hồi tiếp âm đơn vị ngoài khép kín hệ thống. Tìm cực vòng kín p.

G(s) = 8/(s + 3)H₁ = 1K = 11
p =
38.Sơ đồ khốiKhó
R(s)+K(s + 1)(s + 4)(s + 8)G(s)C(s)−

Đối tượng G(s) = K/((s + 1)(s + 4)(s + 8)) nằm trong vòng hồi tiếp âm đơn vị. Dùng tiêu chuẩn Routh–Hurwitz để tìm độ lợi lớn nhất Kmax mà vòng kín còn ổn định.

G(s) = K/((s + 1)(s + 4)(s + 8))hồi tiếp âm đơn vị
Kmax =
39.Đáp ứng bướcTrung bình
t (s)y(t)1

Một hệ có hàm truyền G(s) = 43/(s² + 5s + 43) và nhận tín hiệu vào là bước nhảy đơn vị. Tìm hệ số cản ζ.

G(s) = 43/(s² + 5s + 43)
ζ =
40.Đáp ứng bướcKhó
t (s)y(t)1tp​ = 1 s25 %

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

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

Một vòng phản hồi đơn vị với G(s) = K/(s(s + 2)) phải bám theo tín hiệu dốc r(t) = 1t với sai số xác lập 0,01. Tìm độ lợi K.

G(s) = K/(s(s + 2))r(t) = 1tess = 0,01
K =
42.Sai số xác lập và điều khiển PIDKhó
R(s)+K = 501s(s + 5)C(s)−

Một bộ điều khiển tỉ lệ K = 50 điều khiển đối tượng G(s) = 1/(s(s + 5)) trong một vòng phản hồi đơn vị. Tìm thời gian đỉnh tp của đáp ứng bước.

G(s) = 1/(s(s + 5))K = 50
tp = s

3Dung dịch

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

Đáp số: RA = 47,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.

Hợp lực của tải tam giác

W = ½·15·10 = 75 kN
x̄ = 2L/3 = 6,667 m

Cân bằng

Mômen quanh A, chiều ngược kim đồng hồ là dương. Sau đó cân bằng lực theo phương thẳng đứng.

ΣMA = 0: RB·10 = 75·6,667 + 25·1
RB = 52,5 kN
RA = 100 − 52,5 = 47,5 kN

Đáp số

RA = 47,5 kN
2.Phản lực gối của dầm

Đáp số: RA = 50,1 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.

Các thành phần của P

Pv = 55·sin 50° = 42,13 kN
Ph = 55·cos 50° = 35,35 kN

Các phản lực

ΣMA = 0: RB·10 = 42,13·3 + 20·7 → RB = 26,64 kN
Ay = 42,13 + 20 − 26,64 = 35,49 kN
Ax = Ph = 35,35 kN

Hợp lực tại A

RA = √(35,35² + 35,49²) = 50,1 kN

Đáp số

RA = 50,1 kN
3.Ứng suất và biến dạng dọc trục

Đáp số: P = 32,86 kN

Hai vật liệu: khi nối tiếp (cái này sau cái kia) mỗi đoạn chịu cùng một lực và các độ giãn dài cộng lại: δ = P·Σ L/(E·A). Khi song song (cạnh nhau giữa các tấm cứng) cả hai phần co lại cùng một lượng, nên biến dạng của chúng bằng nhau: phần cứng hơn nhận phần tải lớn hơn, Pi = P·EiAi/ΣEA, và σi = Ei·ε. Một thanh bị giữ giữa hai vách cứng với tải đặt tại chỗ nối là siêu tĩnh: hai đoạn chia tải theo tỉ lệ độ cứng k = EA/L của chúng.

Diện tích

A₁ = 490,9 mm², A₂ = 706,9 mm²

Các độ mềm

L₁/(E₁A₁) = 600/(200000·490,9) = 6,112 × 10⁻⁶ mm/N
L₂/(E₂A₂) = 300/(70000·706,9) = 6,063 × 10⁻⁶ mm/N

Tải

P = 0,4/1,217 × 10⁻⁵ = 32855 N = 32,86 kN

Đáp số

P = 32,86 kN
4.Ứng suất và biến dạng dọc trục

Đáp số: δ = 3,556 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 = π·40²/4 = 1257 mm², A2 = π·12²/4 = 113,1 mm²
σ1 = 39,79 MPa
σ2 = 442,1 MPa

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

δ1 = σ1L1/E1 = 39,79·2000/200000 = 0,3979 mm
δ2 = σ2L2/E2 = 442,1·500/70000 = 3,158 mm
δ = δ1 + δ2 = 3,556 mm

Đáp số

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

Đáp số: φ = 5,845 °

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ủa mặt cắt

J = π·40⁴/32 = 251327 mm⁴

Góc xoắn

φ = TL/(GJ) = 1000·1000·1000/(39000·251327) = 0,102 rad
φ = 0,102·180/π = 5,845°

Đáp số

φ = 5,845°
6.Xoắn của trục

Đáp số: φC = 0,4811 °

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·ω. Khi cả hai đầu bị ngàm, tĩnh học chỉ cho TA + TB = T₀; phương trình thứ hai là điều kiện tương thích (góc xoắn của hai phần bằng nhau tại tiết diện chịu tải).

Cân bằng và tương thích

TA + TB = 1000
TA·a = TB·b
TA = 1000·2,5/3 = 833,3 N·m
TB = 1000·0,5/3 = 166,7 N·m

Góc quay của C

Dùng đoạn AC: nó được ngàm tại A và chịu TA.

J = π·60⁴/32 = 1,272 × 10⁶ mm⁴
G = 39 GPa = 39000 MPa
φC = TA·a/(GJ) = 833333·500/(39000·1,272 × 10⁶) = 0,008397 rad = 0,4811°

Đáp số

φC = 0,4811°
7.Lực cắt và mômen uốn

Đáp số: Mmax = 40,31 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. Mô-men lớn nhất: vì dM/dx = V, mô-men uốn đạt cực trị tại nơi lực cắt đi qua không, hoặc tại một lực tập trung hoặc bên trong một tải phân bố.

Phản lực gối

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

ΣMA = 0: RB·4 = 35·1,5 + 10·2
RB = 72,5/4 = 18,13 kN
ΣFy = 0: RA = 45 − 18,13 = 26,88 kN

Mô-men dưới mỗi tải

Giữ phần bên trái của mỗi tải.

dưới P1: M = RA·1,5 = 26,88·1,5 = 40,31 kN·m
dưới P2: M = RB·2 = 18,13·2 = 36,25 kN·m

Giá trị lớn nhất

Mmax = 40,31 kN·m (dưới P1)

Đáp số

Mmax = 40,31 kN·m
8.Lực cắt và mômen uốn

Đáp số: Mx = 116,3 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·ℓ = 10·10 = 100 kN, đặt tại x̄ = 5 m
ΣMA = 0: RB·10 = 30·9 + 10·10·5
RB = 770/10 = 77 kN
ΣFy = 0: RA = 130 − 77 = 53 kN
w = 10 kN/mP = 30 kNABRA​ = 53RB​ = 77mặt cắt7,5 m1,5 m1 m

Cắt tại x = 7,5 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 = 53·7,5 − 10·7,5·3,75 = 116,3 kN·m

Đáp số

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

Đáp số: τmax = 50 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.

στxyCτmax​

Bán kính vòng tròn Mohr

τmax = R = √(((50 − 50)/2)² + −50²) = 50 MPa

Đáp số

τmax = 50 MPa
10.Vòng tròn Mohr và biến đổi ứng suất

Đáp số: τθ = −10 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

τθ = −0·0,866 + 20·−0,5 = −10 MPa

Đáp số

τθ = −10 MPa
11.Độ võng của dầm

Đáp số: δ = 6,667 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 = 4000 mm
EI = 2,4 × 10¹³ N·mm²

Độ võng

δ = wL⁴/(8EI) = 5·4000⁴/(8·2,4 × 10¹³) = 6,667 mm

Đáp số

δ = 6,667 mm
12.Độ võng của dầm

Đáp số: δmax = 24,11 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⁴.

P = 23000 N, w = 15 N/mm, L = 2500 mm, EI = 6,27 × 10¹¹ N·mm²

Từng tải riêng rẽ

δP = PL³/(48EI) = 23000·2500³/(48·6,27 × 10¹¹) = 11,94 mm
δw = 5wL⁴/(384EI) = 5·15·2500⁴/(384·6,27 × 10¹¹) = 12,17 mm

Cộng lại

δmax = 11,94 + 12,17 = 24,11 mm

Đáp số

δmax = 24,11 mm
13.Mất ổn định cột

Đáp số: Pcr = 52,99 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

Mặt cắt tròn có cùng I đối với mọi trục.

I = πd⁴/64 = π·50⁴/64 = 306796 mm⁴

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

KL = 1·2000 = 2000 mm

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

Pcr = π²EI/(KL)² = π²·70000·306796/2000² = 52989 N = 52,99 kN

Đáp số

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

Đáp số: σcr = 197,4 MPa

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

Mặt cắt tròn có cùng I đối với mọi trục.

I = πd⁴/64 = π·40⁴/64 = 125664 mm⁴

Bán kính quán tính và độ mảnh

A = 1257 mm²
r = √(I/A) = √(125664/1257) = 10 mm
KL = 0,5·2000 = 1000 mm
KL/r = 1000/10 = 100

Ứng suất tới hạn

σcr = π²E/(KL/r)² = π²·200000/100² = 197,4 MPa

Đáp số

σcr = 197,4 MPa
15.Chuyển động ném

Đáp số: H = 7,492 m

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 θ = 14·cos 60° = 7 m/s
vy = v0·sin θ = 14·sin 60° = 12,12 m/s

Điểm cao nhất

Ở đỉnh, vận tốc thẳng đứng bằng không: vy² = 2gH.

H = vy²/(2g) = 12,12²/(2·9,81) = 7,492 m

Đáp số

H = 7,492 m
16.Chuyển động ném

Đáp số: x = 138,4 m

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 θ = 35·cos 45° = 24,75 m/s
vy = v0·sin θ = 35·sin 45° = 24,75 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 = 15 + 24,75·t − 4,905·t²
t = (vy + √(vy² + 2gh))/g = (24,75 + √(612,5 + 294,3))/9,81 = 5,592 s

Khoảng cách theo phương ngang

x = vx·t = 24,75·5,592 = 138,4 m

Đáp số

x = 138,4 m
17.Công và năng lượng

Đáp số: Wf = 80,02 J

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).

Lực ma sát

N = mg·cos θ = 49,05·0,9063 = 44,45 N
F = μN = 0,15·44,45 = 6,668 N

Công chống ma sát

Wf = F·d = 6,668·12 = 80,02 J

Đáp số

Wf = 80,02 J
18.Công và năng lượng

Đáp số: x = 35,02 cm

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).

Năng lượng còn lại sau đoạn nhám

Mặt dốc nhẵn, nên chỉ có đoạn nhám lấy đi năng lượng.

mgh − μmgL = 49,05·(3 − 0,25·2) = 122,6 J

Độ nén cực đại

½kx² = 122,6 → x = √(2·122,6/2000) = 0,3502 m = 35,02 cm

Đáp số

x = 35,02 cm
19.Dao động cơ học

Đáp số: ζ = 0,2501

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·√(150000·50) = 5477 N·s/m

Hệ số cản

ζ = c/cc = 1370/5477 = 0,2501

Đáp số

ζ = 0,2501

ζ < 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ố: fn = 2,792 Hz

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

Chỉ k1 và k2 song song. Cặp đó nối tiếp với k3, nên gộp chúng theo nghịch đảo.

k12 = 25 + 25 = 50 kN/m
keq = 50·80/(50 + 80) = 30,77 kN/m

Tần số riêng

ωn = √(30769/100) = 17,54 rad/s
fn = 17,54/(2π) = 2,792 Hz

Đáp số

fn = 2,792 Hz
21.Động cơ nhiệt và máy lạnh

Đáp số: ηII = 55,15 %

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 = 850 + 273,15 = 1123 K
TL = 15 + 273,15 = 288,1 K

Hiệu suất Carnot

ηCarnot = 1 − 288,1/1123 = 74,34 %

Hiệu suất theo định luật hai

ηII = η/ηCarnot = 41/74,34 = 55,15 %

Đáp số

ηII = 55,15 %
22.Động cơ nhiệt và máy lạnh

Đáp số: Ẇmin = 0,4912 kW

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 = 20 + 273,15 = 293,1 K
TL = 2 + 273,15 = 275,1 K

COP tốt nhất có thể

COPHP,max = TH/(TH − TL) = 293,1/(293,1 − 275,1) = 16,29

Công suất tối thiểu

Máy thuận nghịch cần công ít nhất, và nó có COP cao nhất.

Ẇmin = Q̇H/COP = 8/16,29 = 0,4912 kW

Đáp số

Ẇmin = 0,4912 kW
23.Chu trình công suất khí

Đáp số: T3 = 2100 K

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)).

Cuối kỳ nén

T2 = 305·220,4 = 1050 K

Cấp nhiệt đẳng áp

T3 = rc·T2 = 2·1050 = 2100 K

Đáp số

T3 = 2100 K
24.Chu trình công suất khí

Đáp số: η = 48,21 %

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.

Hiệu suất nhiệt

η = 1 − 1/100,2857 = 1 − 1/1,931 = 48,21%

Đáp số

η = 48,21%
25.Phương trình liên tục và Bernoulli

Đáp số: p2 = 162,5 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·(200/100)² = 4 m/s
½ρ(v1² − v2²) = ½·1000·(1² − 4²)/1000 = −7,5 kPa
p2 = p1 + ½ρ(v1² − v2²) = 170 + (−7,5) = 162,5 kPa

Đáp số

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

Đáp số: Q = 11,81 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.

Vận tốc trong ống

Mặt thoáng và cửa ra đều thông với không khí, nên các số hạng áp suất triệt tiêu.

v = √(2g·4,5) = √(2·9,81·4,5) = 9,396 m/s

Lưu lượng

Q = A·v = 0,001257·9,396 = 0,01181 m³/s = 11,81 L/s

Đáp số

Q = 11,81 L/s
27.Ma sát đường ống và tổn thất cột áp

Đáp số: Δp = 345,6 kPa

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.

Tổn thất cột áp

hf = f·(L/D)·V²/(2g) = 0,03·(200/0,05)·2,4²/(2·9,81) = 35,23 m

Độ giảm áp suất

Ống nằm ngang và vận tốc không đổi, nên phương trình năng lượng rút gọn thành Δp = ρ·g·hf.

Δp = ρ·g·hf = 1000·9,81·35,23 = 345600 Pa = 345,6 kPa

Đáp số

Δp = 345,6 kPa
28.Ma sát đường ống và tổn thất cột áp

Đáp số: Δz = 2,54 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,1/0,07069 = 1,415 m/s
V²/(2g) = 0,102 m

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

f·L/D = 0,025·150/0,3 = 12,5
ΣK = 0,5 + 3·0,3 + 10 + 1,0 = 12,4

Tổng tổn thất cột áp

hf = 12,5·0,102 = 1,275 m
hm = 12,4·0,102 = 1,265 m
Δz = hf + hm = 2,54 m

Đáp số

Δz = 2,54 m
29.Máy bơm và công suất bơm

Đáp số: Pel = 9,788 kW

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,015·48 = 7,063 kW

Công suất trên trục

Pshaft = 7,063/0,82 = 8,614 kW

Công suất điện

Pel = 8,614/0,88 = 9,788 kW

Đáp số

Pel = 9,788 kW
30.Máy bơm và công suất bơm

Đáp số: H = 33,2 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,04/0,01767 = 2,264 m/s
V²/(2g) = 0,2611 m

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

hL = (f·L/D + ΣK)·V²/(2g) = (100 + 8)·0,2611 = 28,2 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 = 5 + 28,2 = 33,2 m

Đáp số

H = 33,2 m
31.Dẫn nhiệt

Đáp số: Ti = 14,79 °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) = (100/1000)/(0,96·15) = 0,006944 K/W
R2 = L2/(k2A) = (75/1000)/(0,04·15) = 0,125 K/W

Dòng nhiệt

Q̇ = (T1 − T3)/(R1 + R2) = 23/0,1319 = 174,3 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 = 16 − 174,3·0,006944 = 14,79 °C

Đáp số

Ti = 14,79 °C
32.Dẫn nhiệt

Đáp số: Q̇ = 16,72 W

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/(7,7·4) = 0,03247 K/W
R1 = L1/(k1A) = (200/1000)/(0,035·4) = 1,429 K/W
R2 = L2/(k2A) = (100/1000)/(0,96·4) = 0,02604 K/W
Ro = 1/(hoA) = 1/(30·4) = 0,008333 K/W
ΣR = 1,495 K/W

Nhiệt thất thoát

Q̇ = (20 − (-5))/1,495 = 16,72 W

Đáp số

Q̇ = 16,72 W
33.Bánh răng và hệ bánh răng

Đáp số: T1 = 19,89 N·m

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.

Mômen xoắn động cơ

ω₁ = 2π·1440/60 = 150,8 rad/s
T1 = P/ω₁ = 3000/150,8 = 19,89 N·m

Đáp số

T1 = 19,89 N·m
34.Bánh răng và hệ bánh răng

Đáp số: θ = 152,8 °

Truyền động đai: khi không trượt, cả hai puli có cùng tốc độ mép v = π·d·n/60, nên n₁·d₁ = n₂·d₂ và tỉ số tốc độ là d₂/d₁. Độ trượt s làm puli bị dẫn chậm hơn một chút: n₂ = n₁·(d₁/d₂)·(1 − s). Công suất truyền là P = (T₁ − T₂)·v, với T₁ là lực căng ở nhánh căng và T₂ ở nhánh chùng. Ma sát giới hạn tỉ số các lực căng: T₁/T₂ = eμθ, trong đó θ là góc ôm tính theo radian trên puli nhỏ hơn; với đai hở θ = π − 2·asin((d₂ − d₁)/(2C)) với khoảng cách tâm C.

Góc của dây đai

asin(0,235) = 13,59°

Góc ôm

θ = 180 − 2·13,59 = 152,8°

Đáp số

θ = 152,8°
35.Biến đổi Laplace

Đáp số: f(0.6) = −0,5239

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 + 34 = (s + 3)² + 25 = (s + 3)² + 5²

Ghép tử số

F(s) = 4·(s + 3)/((s + 3)² + 5²) + (28/5)·5/((s + 3)² + 5²)

Biến đổi ngược

f(t) = e−3t·(4·cos 5t + 5,6·sin 5t)
f(0,6) = e−3·0,6·(4·cos(5·0,6) + 5,6·sin(5·0,6)) = −0,5239

Đáp số

f(0,6) = −0,5239
36.Biến đổi Laplace

Đáp số: g(0.3) = 1,764

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).

Hàm truyền

G(s) = Y(s)/U(s) = (4s + 8)/(s² + 2s + 17)

Phân thức riêng

s² + 2s + 17 = (s + 1)² + 4²
G(s) = 4·(s + 1)/((s + 1)² + 4²) + (4/4)·4/((s + 1)² + 4²)

Đáp ứng xung

g(t) = e−t·(4·cos 4t + sin 4t)
g(0,3) = 1,764

Đáp số

g(0,3) = 1,764
37.Sơ đồ khối

Đáp số: p = −99

Đạ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).

Vòng trong

Gᵢ = G/(1 + G·H₁) = 8/(s + 3 + 8·1) = 8/(s + 11)

Vòng ngoài

T(s) = K·Gᵢ/(1 + K·Gᵢ) = 88/(s + 11 + 88) = 88/(s + 99)

Kết quả

p = −99

Đáp số

p = −99
38.Sơ đồ khối

Đáp số: Kmax = 540

Đạ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

s³ + 13s² + 44s + 32 + K = 0

Bảng Routh

s³ | 1 44
s² | 13 32 + K
s¹ | (13·44 − 32 − K)/13
s⁰ | 32 + K

Giới hạn ổn định

Mọi phần tử của cột đầu tiên phải dương.

Kmax = 13·44 − 32 = 540

Đáp số

Kmax = 540
39.Đáp ứng bước

Đáp số: ζ = 0,3812

Đá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 + ωₙ².

ωₙ² = 43 → ωₙ = 6,557 rad/s
2ζωₙ = 5 → ζ = 5/(2·6,557) = 0,3812

Kết quả

ζ = 0,3812

Đáp số

ζ = 0,3812
40.Đáp ứng bước

Đáp số: ωn = 3,434 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,25)/√(π² + ln²(0,25)) = 1,386/√(9,87 + 1,922) = 0,4037

Tần số riêng

ωd = π/tp = π/1 = 3,142 rad/s
ωₙ = ωd/√(1 − ζ²) = 3,142/0,9149 = 3,434 rad/s

Đáp số

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

Đáp số: K = 200

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 yêu cầu

Kv = A/ess = 1/0,01 = 100

Độ khuếch đại

K = 2·Kv = 200

Đáp số

K = 200

Một độ lợi lớn làm vòng này cản rất yếu: ζ = 0,07071. Mâu thuẫn đó là lý do dùng các bộ bù.

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

Đáp số: tp = 0,475 s

Đ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

s(s + 5) + 50·1 = s² + 5s + 50

ωₙ và ζ

ωₙ = √50 = 7,071 rad/s
ζ = 5/(2·7,071) = 0,3536

Thời gian đỉnh

ωd = 7,071·√(1 − 0,3536²) = 6,614 rad/s
tp = π/6,614 = 0,475 s

Đáp số

tp = 0,475 s

Đâ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