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​ = 10 kNP2​ = 35 kNABRA​RB​ = ?4 m3 m3 m

Dầm đơn giản có gối cố định tại A và gối di động tại B. Tìm phản lực thẳng đứng RB.

L = 10 mP1 = 10 kN tại 4 mP2 = 35 kN tại 7 m
RB = kN
2.Phản lực gối của dầmKhó
w0​ = 8 kN/mP = 40 kNABMA​ = ?6 m2 m

Dầm công-xôn được ngàm vào tường tại A. Nó chịu tải thay đổi tuyến tính bằng 8 kN/m tại tường và giảm về không tại đầu tự do B, cùng một lực tập trung P tại vị trí cách tường 6 m. Tìm độ lớn của mô-men phản lực MA tại tường.

L = 8 mw0 = 8 kN/mP = 40 kN tại 6 m
MA = kN·m
3.Ứng suất và biến dạng dọc trụcTrung bình
Ø 25 mmNhôm1 mP = 75 kN

Tìm biến dạng dọc trục ε trong thanh nhôm, tính bằng microstrain (µε = 10⁻⁶).

P = 75 kNd = 25 mmL = 1 mE = 70 GPa
ε = µε
4.Ứng suất và biến dạng dọc trụcKhó
80 kN12300 mm200 mmØ25Ø50σ2​ = ?

Một đoạn thép (đường kính 25 mm, chiều dài 300 mm, E = 200 GPa) và một đoạn làm bằng đồng thanh (đường kính 50 mm, chiều dài 200 mm, E = 110 GPa) được nối đầu với đầu và ngàm giữa hai vách cứng. Một tải dọc trục 80 kN tác dụng tại chỗ nối, hướng về phía đoạn 2. Tìm độ lớn ứng suất σ2 trong đoạn 2.

đoạn 1: thép, Ø25 mm, 300 mm, E = 200 GPađoạn 2: đồng thanh, Ø50 mm, 200 mm, E = 110 GPaP = 80 kN
σ2 = MPa
5.Xoắn của trụcTrung bình
T = 400 N·md = 30 mmL = 1,2 mφ = ?

Một trục đặc thép đượ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 = 400 N·md = 30 mmL = 1,2 mG = 80 GPa
φ = °
6.Xoắn của trụcKhó
Td = ?

Một trục đặc phải truyền 50 kW ở 750 rpm mà ứng suất tiếp không vượt quá 60 MPa. Tìm đường kính nhỏ nhất d thỏa mãn.

P = 50 kWN = 750 rpmτallow = 60 MPa
d = mm
7.Lực cắt và mômen uốnTrung bình
4 kN/m20 kNABCDEFG1 m1 m1 m1 m1 m1 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ơ đồ.

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

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

Một điểm trong một tấm mỏng ở trạng thái ứng suất phẳng: σₓ = 160 MPa, σᵧ = 70 MPa và τₓᵧ = −60 MPa, với σz = 0. Tìm ứng suất cắt lớn nhất tuyệt đối τabs tại điểm đó (xét mọi mặt trong không gian ba chiều).

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

Một điểm trong chi tiết chịu tải có σₓ = 30 MPa, σᵧ = 50 MPa và τₓᵧ = −10 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 20° theo chiều kim đồng hồ so với trục x.

σx = 30 MPaσy = 50 MPaτxy = −10 MPaθ = -20°
σθ = MPa
11.Độ võng của dầmTrung bình
w = 15 kN/mABθ8 m

Một dầm đơn giản chịu tải phân bố đều w trên toàn bộ nhịp. Tìm góc xoay θA của dầm tại gối A, theo radian.

w = 15 kN/mL = 8 mE = 200 GPa (thép)I = 36 × 10⁶ mm⁴
θA = rad
12.Độ võng của dầmKhó
w = 12 kN/mAB5 m

Một dầm công-xôn có gối chống ngàm tại A và tựa đơn (chống) tại B. Dầm chịu tải phân bố đều w trên toàn bộ chiều dài. EI không đổi. Tìm phản lực RB tại gối chống.

L = 5 mw = 12 kN/m
RB = kN
13.Mất ổn định cộtTrung bình
PL = 6 m40 × 100 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 = 40 × 100 mm
Pcr = kN
14.Mất ổn định cộtKhó
PL = 1 md = 25 mmmặt cắt ngang

Một cột dài L làm bằng thép được ngàm ở chân và tự do ở đỉnh. 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 = 1 mE = 200 GPa (thép)d = 25 mm
σcr = MPa
15.Chuyển động némTrung bình
v0​ = 32 m/s60°tầm xa

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 thời gian bay t cho đến khi quả bóng chạm đất.

v0 = 32 m/sθ = 60°g = 9,81 m/s²
t = s
16.Chuyển động némKhó
h = 50 mv0​ = 12 m/s20°x = ?

Một quả bóng được ném với góc 20° so với phương ngang từ đỉnh một vách đá cao 50 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 = 12 m/sθ = 20°h = 50 mg = 9,81 m/s²
x = m
17.Công và năng lượngTrung bình
30 km/h

Một ô tô chạy với vận tốc 30 km/h trên đường nằm ngang phanh gấp và trượt cho đến khi dừng với các bánh xe bị khóa. Hệ số ma sát trượt giữa lốp và mặt đường là 0,7. Tìm quãng đường dừng d.

v = 30 km/hμ = 0,7g = 9,81 m/s²
d = m
18.Công và năng lượngKhó
5 kg5 kgμ = 0,2d = 2,5 m

Một khối 5 kg nằm trên mặt bàn nằm ngang nhám (μ = 0,2). Một sợi dây nhẹ đi từ nó qua một ròng rọc nhẵn ở mép bàn tới một khối 5 kg treo. Hệ được thả từ trạng thái nghỉ. Tìm vận tốc v của các khối sau khi khối treo hạ xuống 2,5 m.

m1 = 5 kgm2 = 5 kgμ = 0,2d = 2,5 mg = 9,81 m/s²
v = m/s
19.Dao động cơ họcTrung bình
k = 50 kN/mc = 1000 N·s/mm = 20 kg

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

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

Một vật 100 kg treo vào lò xo k₃ = 25 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₁ = 5 kN/m và k₂ = 20 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 = 5 kN/mk2 = 20 kN/mk3 = 25 kN/m
fn = Hz
21.Động cơ nhiệt và máy lạnhTrung bình
TH​ = 600 °CTL​ = ?EQH​ = 3000 kJQL​ = 1093 kJW

Một động cơ nhiệt thuận nghịch (Carnot) nhận 3000 kJ nhiệt từ một nguồn ở 600 °C và thải 1093 kJ mỗi chu trình. Tìm nhiệt độ TL của nguồn lạnh theo °C.

TH = 600 °CQH = 3000 kJQL = 1093 kJ
TL = °C
22.Động cơ nhiệt và máy lạnhKhó
TH​ = 22 °CTL​ = 7 °CHPQ̇H​ = 4 kWQ̇L​Ẇ = ?

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

TH = 22 °CTL = 7 °CQ̇H = 4 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 = 17 và tỉ số cắt nhiên liệu rc = v₃/v₂ = 3. Khi bắt đầu nén không khí ở 295 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 = 17rc = 3T1 = 295 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 chu trình Brayton lý tưởng hoạt động với tỉ số áp suất máy nén 15. 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 = 15k = 1,4
η = %
25.Phương trình liên tục và BernoulliTrung bình
Q = ?d = 30 mm, Cd​ = 0,63h = 3 m

Nước chảy ra từ một bể hở lớn qua một lỗ có cạnh sắc đường kính 30 mm. Mực nước được giữ cao hơn lỗ 3 m và hệ số lưu lượng là 0,63. Tìm lưu lượng Q. Dùng g = 9,81 m/s².

d = 30 mmh = 3 mCd = 0,63g = 9,81 m/s²
Q = L/s
26.Phương trình liên tục và BernoulliKhó
12d1​ = 75 mm, p1​ = 230 kPad2​ = 50 mm, p2​ = ?Δz = 4 m

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

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

Nước chảy với vận tốc trung bình 2,4 m/s qua một ống dài 150 m, đường kính trong 75 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,4 m/sL = 150 mD = 75 mmf = 0,015g = 9,81 m/s²
hf = m
28.Ma sát đường ống và tổn thất cột ápKhó
K = 2cửa xả tự doH = ?L = 150 m, D = 150 mm, f = 0,0153× khúc cua (K = 0,3)

Một ống dài 150 m, đường kính trong 150 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 = 2). 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 25 L/s?

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

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

Q = 50 L/sH = 40 mPshaft = 26,2 kWg = 9,81 m/s²
η = %
30.Máy bơm và công suất bơmKhó
100 kPa (áp suất dư)PQ = 50 L/sΔz = 6 mL = 600 m, D = 200 mm, f = 0,025, ΣK = 8

Một máy bơm chuyển 50 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 6 m. Ống dài 600 m, đường kính trong 200 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². Bồn phía trên đóng kín và không khí phía trên nước ở áp suất dư 100 kPa. Tìm cột áp H mà bơm phải bổ sung.

Q = 50 L/sΔz = 6 mL = 600 mD = 200 mmf = 0,025ΣK = 8ptank = 100 kPa
H = m
31.Dẫn nhiệtTrung bình
Gỗ, k = 0,1750 mmTấm thạch cao, k = 0,1775 mm20 °C0 °CQ̇ = ?

Một bức tường gồm 50 mm gỗ và 75 mm tấm thạch cao, có diện tích 15 m². Hai mặt ngoài ở 20 °C và 0 °C. Tìm dòng nhiệt Q̇ qua tường.

L1 = 50 mm, k1 = 0,17 W/(m·K)L2 = 75 mm, k2 = 0,17 W/(m·K)A = 15 m²T1 = 20 °CT3 = 0 °C
Q̇ = W
32.Dẫn nhiệtKhó
Tấm thạch cao, k = 0,17250 mmKính, k = 0,96100 mmphòng 21 °Cbên ngoài -15 °CQ̇ = ?

Một bức tường gồm 250 mm tấm thạch cao và 100 mm kính ngăn cách một căn phòng ở 21 °C với không khí bên ngoài ở -15 °C. Hệ số đối lưu là 6 W/(m²·K) bên trong và 30 W/(m²·K) bên ngoài, và diện tích tường là 12 m². Tìm lượng nhiệt thất thoát Q̇ qua tường.

L1 = 250 mm, k1 = 0,17 W/(m·K)L2 = 100 mm, k2 = 0,96 W/(m·K)A = 12 m²hi = 6 W/(m²·K)ho = 30 W/(m²·K)Troom = 21 °CTout = -15 °C
Q̇ = W
33.Bánh răng và hệ bánh răngTrung bình
N1​ = 25N2​ = 48

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

P = 0,75 kWn1 = 1500 rpmN1 = 25N2 = 48η = 97%
T1 = N·m
34.Bánh răng và hệ bánh răngKhó
d1​ = 200 mmd2​ = 250 mm1500 rpmC = 1200 mmT1​ = 800 NT2​P = ?

Một dây đai hở nối puli dẫn 200 mm ở 1500 vòng/phút với puli 250 mm; tâm cách nhau 1200 mm. Hệ số ma sát là 0,35 và lực căng nhánh căng không được vượt quá 800 N. Tìm công suất lớn nhất P mà dây đai có thể truyền mà không trượt, theo kW (bỏ qua hiệu ứng ly tâm).

d1 = 200 mmn1 = 1500 rpmd2 = 250 mmC = 1200 mmμ = 0,35T1 = 800 N
P = kW
35.Biến đổi LaplaceTrung bình
ReIm−2−4−6−8−10−12

F(s) = (s + 12)/(s² + 10s + 24). Biến đổi ngược của nó có dạng f(t) = A·e−4t + B·e−6t. Tìm hệ số A bằng phân thức riêng.

F(s) = (s + 12)/(s² + 10s + 24)
A =
36.Biến đổi LaplaceKhó
ReIm−0,5−1−1,5−2−2,5−3−3,5−4−4,50,51

Một hệ thống được mô tả bởi phương trình vi phân y″ + 6y′ + 8y = 3u′ + 11u, 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,4 s. Tìm g(0,4).

y″ + 6y′ + 8y = 3u′ + 11ut = 0,4 s
g(0.4) =
37.Sơ đồ khốiTrung bình
R(s)+30s + 5G(s)C(s)0,5H−

Một đối tượng G(s) = 30/(s + 5) được điều khiển trong vòng hồi tiếp âm mà cảm biến có độ lợi H = 0,5. Một bậc thang có độ lớn 5 được đặt vào R. Tìm đầu ra xác lập css.

G(s) = 30/(s + 5)H = 0,5R = 5
css =
38.Sơ đồ khốiKhó
R(s)+29bộ điều khiển+D(s)+9s + 8đối tượngC(s)−

Một bộ điều khiển tỷ lệ K = 29 điều khiển đối tượng G(s) = 9/(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 5 và một bậc thang nhiễu có độ lớn 2 được đặt cùng lúc. Tìm đầu ra xác lập css.

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

Một hệ có hàm truyền G(s) = 49/(s² + 8s + 49) và nhận tín hiệu vào là bước nhảy đơn vị. Tìm tần số riêng ωn theo rad/s.

G(s) = 49/(s² + 8s + 49)
ωn = rad/s
40.Đáp ứng bướcKhó
t (s)y(t)1tp​ = 0,25 s30 %

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

%OS = 30 %tp = 0,25 s
ωn = rad/s
41.Sai số xác lập và điều khiển PIDTrung bình
R(s)+K = 1025s + 1C(s)−

Một quá trình G(s) = 2/(5s + 1) được điều khiển bởi độ lợi tỉ lệ K = 10 trong một vòng phản hồi đơn vị. Tìm sai số xác lập ess cho một bước nhảy, theo phần trăm độ lớn bước nhảy.

G(s) = 2/(5s + 1)K = 10
ess = %
42.Sai số xác lập và điều khiển PIDKhó
R(s)+20 + Kd​s10s(s + 4)C(s)−

Một bộ điều khiển PD C(s) = 20 + Kd·s điều khiển đối tượng G(s) = 10/(s(s + 4)) trong một vòng phản hồi đơn vị. Tìm độ lợi vi phân Kd cho vòng kín có tỷ số cản ζ = 0,9.

G(s) = 10/(s(s + 4))Kp = 20ζ = 0,9
Kd = s

3Dung dịch

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

Đáp số: RB = 28,5 kN

Cân bằng tĩnh: một vật đứng yên có ΣFx = 0, ΣFy = 0 và ΣM = 0 quanh mọi điểm. Lấy mômen quanh một gối sẽ loại phản lực của gối đó khỏi phương trình.

P1​ = 10 kNP2​ = 35 kNABRA​RB​4 m3 m3 m

Sơ đồ vật thể tự do

Thay các gối đỡ bằng phản lực. Gối di động cho một phản lực thẳng đứng; khi chỉ có tải thẳng đứng, gối cố định cũng cho một phản lực thẳng đứng.

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 = 10·4 + 35·7
RB = 285/10 = 28,5 kN
ΣFy = 0: RA = 45 − 28,5 = 16,5 kN

Đáp số

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

Đáp số: MA = 325,3 kN·m

Cân bằng của dầm công-xôn: một gối ngàm có thể đẩy lên và còn chống lại sự quay, nên nó cho một phản lực thẳng đứng RA và một mô-men phản lực MA. ΣFy = 0 cho RA, và mô-men đối với tường cho MA = Σ(tải × khoảng cách đến tường). Tải phân bố tác dụng như hợp lực của nó tại trọng tâm của biểu đồ tải.

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

W = ½·8·8 = 32 kN
x̄ = L/3 = 2,667 m

Mô-men đối với tường

MA = W·x̄ + P·a = 32·2,667 + 40·6 = 325,3 kN·m

Đáp số

MA = 325,3 kN·m
3.Ứng suất và biến dạng dọc trục

Đáp số: ε = 2183 µε

Đị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 = π·25²/4 = 490,9 mm²
σ = P/A = 75·1000/490,9 = 152,8 MPa

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

E for nhôm = 70 GPa = 70000 MPa.

ε = σ/E = 152,8/70000 = 0,002183 = 2183 µε

Đáp số

ε = 2183 µε
4.Ứng suất và biến dạng dọc trục

Đáp số: σ2 = 31,27 MPa

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.

Các độ cứng

k₁ = 200000·490,9/300 = 327249 N/mm
k₂ = 110000·1963/200 = 1,08 × 10⁶ N/mm

Phần của đoạn 2

F2 = 80000·1,08 × 10⁶/1,407 × 10⁶ = 61395 N

Ứng suất

σ2 = 61395/1963 = 31,27 MPa (nén)

Đáp số

σ2 = 31,27 MPa
5.Xoắn của trục

Đáp số: φ = 4,323 °

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 = π·30⁴/32 = 79522 mm⁴

Góc xoắn

φ = TL/(GJ) = 400·1000·1200/(80000·79522) = 0,07545 rad
φ = 0,07545·180/π = 4,323°

Đáp số

φ = 4,323°
6.Xoắn của trục

Đáp số: d = 37,81 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π·750/60 = 78,54 rad/s
T = P/ω = 50000/78,54 = 636,6 N·m

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

d = ∛(16T/(π·τallow)) = ∛(16·636620/(π·60)) = 37,81 mm

Đáp số

d = 37,81 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.

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

Đáp số: D (x = 3 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 = 22 kN, R2 = 22 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 = 20, MC = 36, MD = 48 kN·m, ME = 36, MF = 20, MG = 0

Đáp số

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

Đáp số: Mx = 30 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·ℓ = 12·4 = 48 kN, đặt tại x̄ = 2 m
ΣMA = 0: RB·4 = 20·3 + 12·4·2
RB = 156/4 = 39 kN
ΣFy = 0: RA = 68 − 39 = 29 kN
w = 12 kN/mP = 20 kNABRA​ = 29RB​ = 39mặt cắt1,5 m1,5 m1 m

Cắt tại x = 1,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 = 29·1,5 − 12·1,5·0,75 = 30 kN·m

Đáp số

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

Đáp số: τabs = 95 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. Trong không gian ba chiều, ứng suất cắt lớn nhất tuyệt đối là (σmax − σmin)/2 trên cả ba ứng suất chính; trong ứng suất phẳng, ứng suất thứ ba bằng không.

στxyCσ1​σ2​

Các ứng suất chính trong mặt phẳng

C = 115 MPa, R = √(45² + (−60)²) = 75 MPa
σ1 = 190 MPa, σ2 = 40 MPa, σz = 0

Vòng tròn lớn nhất trong ba vòng tròn

τ = max(|σ1 − σ2|, |σ1|, |σ2|)/2 = max(150, 190, 40)/2 = 95 MPa

Đáp số

τabs = 95 MPa

Ở đây vòng tròn ngoài mặt phẳng quyết định: giá trị lớn nhất tuyệt đối lớn hơn giá trị lớn nhất trong mặt phẳng R.

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

Đáp số: σθ = 38,77 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

σθ = 40 + −10·0,766 + −10·−0,6428 = 38,77 MPa

Đáp số

σθ = 38,77 MPa
11.Độ võng của dầm

Đáp số: θA = 0,04444 rad

Độ 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). Góc xoay θ = dv/dx là tích phân bậc một của M/(EI); với các đơn vị nhất quán, nó có đơn vị radian.

Đơ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 = 15 N/mm, L = 8000 mm
EI = 7,2 × 10¹² N·mm²

Góc xoay ở đầu

θA = wL³/(24EI) = 15·8000³/(24·7,2 × 10¹²) = 0,04444 rad

Đáp số

θA = 0,04444 rad

Đó là 2,546°.

12.Độ võng của dầm

Đáp số: RB = 22,5 kN

Độ 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). Với dầm siêu tĩnh, hãy bỏ một liên kết thừa và đòi hỏi độ võng tại đó bằng không (điều kiện tương thích).

Bỏ gối chống đi

Khi không có gối chống, dầm là một dầm công-xôn. Đầu tự do của nó sẽ võng dưới tải.

δload = wL⁴/(8EI)

Điều kiện tương thích tại gối chống

Chỉ riêng lực chống RB nâng đầu tự do lên một đoạn RBL³/(3EI). Gối chống giữ B có độ võng bằng không, nên hai độ võng phải triệt tiêu nhau. EI được khử.

RB·L³/(3EI) = wL⁴/(8EI) → RB = 3wL/8 = 3·12·5/8 = 22,5 kN

Đáp số

RB = 22,5 kN
13.Mất ổn định cột

Đáp số: Pcr = 59,68 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 = 100·40³/12 = 533333 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·533333/4200² = 59680 N = 59,68 kN

Đáp số

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

Đáp số: σcr = 19,28 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 = π·25⁴/64 = 19175 mm⁴

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

A = 490,9 mm²
r = √(I/A) = √(19175/490,9) = 6,25 mm
KL = 2·1000 = 2000 mm
KL/r = 2000/6,25 = 320

Ứng suất tới hạn

σcr = π²E/(KL/r)² = π²·200000/320² = 19,28 MPa

Đáp số

σcr = 19,28 MPa
15.Chuyển động ném

Đáp số: t = 5,65 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 θ = 32·cos 60° = 16 m/s
vy = v0·sin θ = 32·sin 60° = 27,71 m/s

Thời gian trong không khí

Xuất phát và rơi xuống cùng độ cao: t = 2vy/g.

t = 2·27,71/9,81 = 5,65 s

Đáp số

t = 5,65 s
16.Chuyển động ném

Đáp số: x = 41,03 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 θ = 12·cos 20° = 11,28 m/s
vy = v0·sin θ = 12·sin 20° = 4,104 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 = 50 + 4,104·t − 4,905·t²
t = (vy + √(vy² + 2gh))/g = (4,104 + √(16,84 + 981))/9,81 = 3,638 s

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

x = vx·t = 11,28·3,638 = 41,03 m

Đáp số

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

Đáp số: d = 5,056 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).

Vận tốc theo đơn vị cơ bản

v = 30/3,6 = 8,333 m/s

Cân bằng năng lượng

½mv² = μmg·d → d = v²/(2μg) = 8,333²/(2·0,7·9,81) = 5,056 m

Đáp số

d = 5,056 m

Tăng gấp đôi vận tốc làm quãng đường dừng dài gấp bốn lần.

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

Đáp số: v = 4,429 m/s

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 giải phóng và mất đi

m₂g·d = 49,05·2,5 = 122,6 J
μ·m₁g·d = 0,2·49,05·2,5 = 24,53 J

Động năng của cả hai khối

½(m₁ + m₂)v² = 98,1 → v = √(2·98,1/10) = 4,429 m/s

Đáp số

v = 4,429 m/s
19.Dao động cơ học

Đáp số: ζ = 0,5

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·√(50000·20) = 2000 N·s/m

Hệ số cản

ζ = c/cc = 1000/2000 = 0,5

Đáp số

ζ = 0,5

ζ < 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 = 1,779 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 = 5 + 20 = 25 kN/m
keq = 25·25/(25 + 25) = 12,5 kN/m

Tần số riêng

ωn = √(12500/100) = 11,18 rad/s
fn = 11,18/(2π) = 1,779 Hz

Đáp số

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

Đáp số: TL = 44,97 °C

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 của nguồn

TH = 600 + 273,15 = 873,1 K

Tỉ số nhiệt thuận nghịch

TL = TH·QL/QH = 873,1·1093/3000 = 318,1 K = 44,97 °C

Đáp số

TL = 44,97 °C
22.Động cơ nhiệt và máy lạnh

Đáp số: Ẇmin = 0,2033 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 = 22 + 273,15 = 295,1 K
TL = 7 + 273,15 = 280,1 K

COP tốt nhất có thể

COPHP,max = TH/(TH − TL) = 295,1/(295,1 − 280,1) = 19,68

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 = 4/19,68 = 0,2033 kW

Đáp số

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

Đáp số: qin = 1842 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 = 295·170,4 = 916,2 K
T3 = 3·916,2 = 2749 K

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

qin = 1,005·(2749 − 916,2) = 1842 kJ/kg

Đáp số

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

Đáp số: η = 53,87 %

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/150,2857 = 1 − 1/2,168 = 53,87%

Đáp số

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

Đáp số: Q = 3,417 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 tia lý tưởng

√(2gh) = √(2·9,81·3) = 7,672 m/s

Lưu lượng thực

Hệ số lưu lượng tính đến sự co hẹp của tia và ma sát ở mép lỗ.

Q = Cd·A·√(2gh) = 0,63·0,0007069·7,672 = 0,003417 m³/s = 3,417 L/s

Đáp số

Q = 3,417 L/s
26.Phương trình liên tục và Bernoulli

Đáp số: p2 = 188,7 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/50)² = 2,25 m/s
½ρ(v1² − v2²) = ½·1000·(1² − 2,25²)/1000 = −2,031 kPa
ρgΔz = 1000·9,81·4/1000 = 39,24 kPa
p2 = p1 + ½ρ(v1² − v2²) − ρgΔz = 230 + (−2,031) − 39,24 = 188,7 kPa

Đáp số

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

Đáp số: hf = 8,807 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,4²/(2·9,81) = 0,2936 m

Darcy–Weisbach

hf = f·(L/D)·V²/(2g) = 0,015·(150/0,075)·2,4²/(2·9,81) = 8,807 m

Đáp số

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

Đáp số: H = 1,979 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,025/0,01767 = 1,415 m/s
V²/(2g) = 0,102 m

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

f·L/D = 0,015·150/0,15 = 15
ΣK = 0,5 + 3·0,3 + 2 = 3,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 + 15 + 3,4)·0,102 = 1,979 m

Đáp số

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

Đáp số: η = 74,89 %

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,05·40 = 19,62 kW

Hiệu suất

η = Ph/Pshaft = 19,62/26,2 = 74,89%

Đáp số

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

Đáp số: H = 26,91 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,05/0,03142 = 1,592 m/s
V²/(2g) = 0,1291 m

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

hL = (f·L/D + ΣK)·V²/(2g) = (75 + 8)·0,1291 = 10,72 m

Cột áp áp suất

p/(ρ·g) = 100000/(1000·9,81) = 10,19 m

Cột áp của bơm

Cả hai mặt nước thực tế đều đứng yên, nên không có số hạng cột áp vận tốc.

H = Δz + p/(ρ·g) + hL = 6 + 10,19 + 10,72 = 26,91 m

Đáp số

H = 26,91 m
31.Dẫn nhiệt

Đáp số: Q̇ = 408 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ó.

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

R1 = L1/(k1A) = (50/1000)/(0,17·15) = 0,01961 K/W
R2 = L2/(k2A) = (75/1000)/(0,17·15) = 0,02941 K/W

Dòng nhiệt

Q̇ = (T1 − T3)/(R1 + R2) = 20/0,04902 = 408 W

Đáp số

Q̇ = 408 W
32.Dẫn nhiệt

Đáp số: Q̇ = 243,4 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/(6·12) = 0,01389 K/W
R1 = L1/(k1A) = (250/1000)/(0,17·12) = 0,1225 K/W
R2 = L2/(k2A) = (100/1000)/(0,96·12) = 0,008681 K/W
Ro = 1/(hoA) = 1/(30·12) = 0,002778 K/W
ΣR = 0,1479 K/W

Nhiệt thất thoát

Q̇ = (21 − (-15))/0,1479 = 243,4 W

Đáp số

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

Đáp số: T1 = 4,775 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π·1500/60 = 157,1 rad/s
T1 = P/ω₁ = 750/157,1 = 4,775 N·m

Đáp số

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

Đáp số: P = 8,32 kW

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 ôm

θ = 3,1 rad

Lực căng nhánh chùng

eμθ = e1,085 = 2,959, T₂ = 800/2,959 = 270,3 N

Công suất

P = (800 − 270,3)·15,71 = 8,32 kW

Đáp số

P = 8,32 kW
35.Biến đổi Laplace

Đáp số: A = 4

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

Phân tích mẫu số thành nhân tử

s² + 10s + 24 = (s + 4)(s + 6)
F(s) = A/(s + 4) + B/(s + 6)

Quy tắc che

Nhân với (s + 4) và đặt s = −4: số hạng B biến mất.

A = (s + 12)/(s + 6) |s = −4 = 8/2 = 4

Đáp số

A = 4
(B = −3, nên f(t) = 4·e−4t − 3·e−6t)
36.Biến đổi Laplace

Đáp số: g(0.4) = 1,224

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) = (3s + 11)/(s² + 6s + 8)

Phân thức riêng

s² + 6s + 8 = (s + 2)(s + 4)
A = 2,5, B = 0,5

Đáp ứng xung

g(t) = 2,5·e−2t + 0,5·e−4t
g(0,4) = 1,224

Đáp số

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

Đáp số: css = 7,5

Đạ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 vòng kín

T(s) = G/(1 + G·H) = 30/(s + 5 + 15) = 30/(s + 20)

Kết quả

css = R·T(0) = 5·30/20 = 7,5

Đáp số

css = 7,5
38.Sơ đồ khối

Đáp số: css = 4,918

Đạ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) = 9/(s + 8 + 261)
C/R = K·G/(1 + K·G) = 261/(s + 8 + 261)

Trạng thái ổn định

css = 5·261/269 + 2·9/269 = 4,918

Đáp số

css = 4,918

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ố: ωn = 7 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).

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

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

ωₙ² = 49 → ωₙ = 7 rad/s
2ζωₙ = 8 → ζ = 8/(2·7) = 0,5714

Kết quả

ωₙ = 7 rad/s

Đáp số

ωn = 7 rad/s
40.Đáp ứng bước

Đáp số: ωn = 13,46 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,3)/√(π² + ln²(0,3)) = 1,204/√(9,87 + 1,45) = 0,3579

Tần số riêng

ωd = π/tp = π/0,25 = 12,57 rad/s
ωₙ = ωd/√(1 − ζ²) = 12,57/0,9338 = 13,46 rad/s

Đáp số

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

Đáp số: ess = 4,762 %

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ị trí

Kp = 10·2 = 20

Sai số bước nhảy

ess = 100/(1 + 20) = 4,762 %

Đáp số

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

Đáp số: Kd = 2,146 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 + 4) + 10(20 + Kds) = s² + (4 + 10Kd)s + 200

Tần số riêng

ωₙ = √200 = 14,14 rad/s

Độ lợi vi phân

4 + 10Kd = 2·0,9·14,14 = 25,46
Kd = 2,146 s

Đáp số

Kd = 2,146 s

Thành phần vi phân thêm sự cản mà không làm đổi ωₙ.

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