Sức bền vật liệu: 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à 14 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

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

Trọng tâm và mômen quán tính mặt cắt

Trọng tâm
ȳ = ΣAiyi/ΣAi
Hình chữ nhật
I = bh³/12
Định lý trục song song
I = Ic + A·d²
Mặt cắt ghép
I = Σ(Ic,i + Aidi²)

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)

2Bài tập

1.Ứng suất và biến dạng dọc trụcTrung bình
Ø 30 mmĐồng thau2,5 mP = 90 kN

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

P = 90 kNd = 30 mmL = 2,5 mE = 100 GPa
δ = mm
2.Ứng suất và biến dạng dọc trụcKhó
Ø 30 mmNhôm0,5 mØ 12 mmĐồng1 mP = 75 kN

Thanh bậc được ngàm vào tường và chịu tải ở đầu tự do. Tìm ứng suất pháp lớn nhất σmax trong thanh.

P = 75 kNđoạn 1: Ø30 mm, 0,5 m, nhôm (E = 70 GPa)đoạn 2: Ø12 mm, 1 m, đồng (E = 117 GPa)
σmax = MPa
3.Xoắn của trụcTrung bình
d = 20 mmT = 100 N·m400 mmφ = 3,74°

Trong một thí nghiệm xoắn, một mẫu đặc đường kính 20 mm bị xoắn bởi mô-men xoắn 100 N·m. Hai đầu của một đoạn đo dài 400 mm quay 3,74° so với nhau (vật liệu vẫn còn đàn hồi). Tìm mô-đun cắt G của vật liệu.

T = 100 N·md = 20 mmL = 400 mmφ = 3,74°
G = GPa
4.Xoắn của trụcKhó
Td = ?

Một trục đặc phải truyền 10 kW ở 300 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 = 10 kWN = 300 rpmτallow = 60 MPa
d = mm
5.Lực cắt và mômen uốnTrung bình
15 kN45 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ơ đồ.

15 kN @ 3 m45 kN @ 2 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)
6.Lực cắt và mômen uốnKhó
w = 12 kN/mP = 80 kNAB4 m6 m

Dầm đơn giản chịu tải phân bố đều w trên toàn bộ nhịp và một lực tập trung P. Tìm mô-men uốn lớn nhất Mmax trong dầm.

L = 10 mP = 80 kN tại 4 mw = 12 kN/m từ 0 đến 10 m
Mmax = kN·m
7.Trọng tâm và mômen quán tính mặt cắtTrung bình
150 mm300 mm210 mmØ40xIx​ = ?

Một tấm 150 mm × 300 mm có một lỗ tròn đường kính 40 mm trên đường tâm thẳng đứng, tâm lỗ cách đáy 210 mm. Tìm mômen quán tính diện tích Ix quanh trục nằm ngang đi qua trọng tâm.

b = 150 mmh = 300 mmd = 40 mmtâm lỗ: y = 210 mm
Ix = × 10⁶ mm⁴
8.Trọng tâm và mômen quán tính mặt cắtKhó
2103602210xM = 80 kN·mkích thước tính bằng mm

Một tiết diện chữ I đối xứng có các cánh rộng 210 mm và dày 22 mm, chiều cao tổng thể 360 mm và bụng dày 10 mm. Nó chịu mô-men uốn 80 kN·m quanh trục x. Tìm độ lớn của ứng suất uốn σj tại chỗ nối giữa bụng và cánh trên.

B = 210 mmH = 360 mmtf = 22 mmtw = 10 mmM = 80 kN·m
σj = MPa
9.Vòng tròn Mohr và biến đổi ứng suấtTrung bình
σx​ = 100σy​ = −20MPaτxy​ = ?

Một điểm có σₓ = 100 MPa và σᵧ = −20 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 130 MPa?

σx = 100 MPaσy = −20 MPaσ1 = 130 MPa
τxy = MPa
10.Vòng tròn Mohr và biến đổi ứng suấtKhó
σx​ = 140σy​ = 40τxy​ = 20MPa

Một điểm trong một tấm mỏng ở trạng thái ứng suất phẳng: σₓ = 140 MPa, σᵧ = 40 MPa và τₓᵧ = 20 MPa (σz = 0). Vật liệu có giới hạn chảy 275 MPa. Tìm hệ số an toàn n chống chảy theo tiêu chuẩn Tresca (ứng suất cắt lớn nhất).

σx = 140 MPaσy = 40 MPaτxy = 20 MPaσY = 275 MPa
n =
11.Độ võng của dầmTrung bình
w = 10 kN/mδL = 2,5 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 = 10 kN/mL = 2,5 mE = 200 GPa (thép)I = 23 × 10⁶ mm⁴
δmax = mm
12.Độ võng của dầmKhó
w = 4 kN/mδL = 2,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 = 4 kN/mL = 2,5 mE = 70 GPa (nhôm)
Ireq = × 10⁶ mm⁴
13.Mất ổn định cộtTrung bình
PL = 4 md = 80 mmmặt cắt ngang

Một cột dài L làm bằng nhôm đượ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 = 4 mE = 70 GPa (nhôm)d = 80 mm
Pcr = kN
14.Mất ổn định cộtKhó
PL = 5 md = 50 mmmặt cắt ngang

Một cột dài L làm bằng nhôm đượ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 = 5 mE = 70 GPa (nhôm)d = 50 mm
σcr = MPa

3Dung dịch

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

Đáp số: δ = 3,183 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 = π·30²/4 = 706,9 mm²
σ = P/A = 90·1000/706,9 = 127,3 MPa

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

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

ε = σ/E = 127,3/100000 = 0,001273 = 1273 µε

Độ giãn dài

δ = ε·L = 0,001273·2500 mm = 3,183 mm

Đáp số

δ = 3,183 mm
2.Ứng suất và biến dạng dọc trục

Đáp số: σmax = 663,1 MPa

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 = π·30²/4 = 706,9 mm², A2 = π·12²/4 = 113,1 mm²
σ1 = 106,1 MPa
σ2 = 663,1 MPa

Ứng suất lớn nhất

Đoạn mảnh hơn có diện tích nhỏ hơn, nên chịu ứng suất lớn hơn.

σmax = 663,1 MPa

Đáp số

σmax = 663,1 MPa
3.Xoắn của trục

Đáp số: G = 39,01 GPa

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 và góc

J = π·20⁴/32 = 15708 mm⁴
φ = 3,74·π/180 = 0,06528 rad

Mô đun trượt

T = 100 N·m = 100000 N·mm
G = TL/(Jφ) = 100000·400/(15708·0,06528) = 39011 MPa
G = 39,01 GPa

Đáp số

G = 39,01 GPa

So sánh: thép ≈ 80 GPa, đồng thau ≈ 39 GPa, nhôm ≈ 26 GPa.

4.Xoắn của trục

Đáp số: d = 30,01 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π·300/60 = 31,42 rad/s
T = P/ω = 10000/31,42 = 318,3 N·m

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

d = ∛(16T/(π·τallow)) = ∛(16·318310/(π·60)) = 30,01 mm

Đáp số

d = 30,01 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.

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

Đáp số: C (x = 2 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 = 37,5 kN, R2 = 22,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 = 37,5, MC = 75 kN·m, MD = 67,5, ME = 45, MF = 22,5, MG = 0

Đáp số

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

Đáp số: Mmax = 336 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.

Hợp lực của tải phân bố đều: W = w·ℓ = 12·10 = 120 kN, đặt tại x̄ = 5 m
ΣMA = 0: RB·10 = 80·4 + 12·10·5
RB = 920/10 = 92 kN
ΣFy = 0: RA = 200 − 92 = 108 kN

V = 0 ở đâu?

V(4⁻) = RA − w·4 = 60 kN
V(4⁺) = 60 − 80 = −20 kN
Lực cắt đổi dấu tại tải, nên giá trị lớn nhất nằm dưới tải.

Mô-men lớn nhất

Mmax = RA·4 − w·4²/2 = 108·4 − 12·4²/2 = 336 kN·m

Đáp số

Mmax = 336 kN·m

Nó xảy ra tại x = 4 m tính từ A.

7.Trọng tâm và mômen quán tính mặt cắt

Đáp số: Ix = 332,7 × 10⁶ mm⁴

Tiết diện có lỗ: coi lỗ là một diện tích âm. Trọng tâm là ȳ = (A₁y₁ − A₂y₂)/(A₁ − A₂), với A₁, y₁ cho tấm đầy đủ và A₂, y₂ cho lỗ. Với mômen quán tính diện tích quanh trục đi qua trọng tâm, dùng định lý trục song song cho từng phần rồi trừ lỗ: I = (Ic1 + A₁d₁²) − (Ic2 + A₂d₂²), trong đó d là khoảng cách từ trọng tâm riêng của mỗi phần đến trục đi qua ȳ. Hình chữ nhật có Ic = bh³/12 và hình tròn có Ic = πd⁴/64.

Trọng tâm

ȳ = (45000·150 − 1257·210)/(45000 − 1257) = 148,3 mm

Tấm đầy đủ quanh ȳ

150·300³/12 + 45000·1,724² = 337,6 × 10⁶ mm⁴

Lỗ quanh ȳ

π·40⁴/64 + 1257·61,72² = 4,913 × 10⁶ mm⁴

Trừ

Ix = 332,7 × 10⁶ mm⁴

Đáp số

Ix = 332,7 × 10⁶ mm⁴
8.Trọng tâm và mômen quán tính mặt cắt

Đáp số: σj = 43,5 MPa

Diện tích ghép: chia mặt cắt thành các hình chữ nhật. Trọng tâm là trung bình có trọng số theo diện tích ȳ = ΣAᵢyᵢ/ΣAᵢ. Mômen quán tính của mặt cắt đối với trục đi qua trọng tâm bằng tổng bh³/12 của từng phần và số hạng trục song song A·d², trong đó d là khoảng cách từ trọng tâm của phần đó đến trọng tâm chung. Ứng suất uốn là σ = M·y/I.

Hình chữ nhật ngoài trừ hai phần khuyết hai bên

Hai phần khuyết cộng lại rộng (B − tw) và cao hw = H − 2tf, và có chung trục x.

hw = 360 − 2·22 = 316 mm
Ix = (B·H³ − (B − tw)·hw³)/12 = (210·360³ − 200·316³)/12
Ix = 2,906 × 10⁸ mm⁴

Ứng suất uốn

y = 180 − 22 = 158 mm
σ = M·y/I = 8 × 10⁷·158/2,906 × 10⁸ = 43,5 MPa

Đáp số

σj = 43,5 MPa
9.Vòng tròn Mohr và biến đổi ứng suất

Đáp số: τxy = 67,08 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 = (100 + (−20))/2 = 40 MPa
R = σ1 − C = 130 − 40 = 90 MPa
στxyCσ1​σ2​

Ứng suất tiếp

τxy = √(R² − ((σₓ − σᵧ)/2)²) = √(90² − 60²) = 67,08 MPa

Đáp số

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

Đáp số: n = 1,912

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. Tiêu chuẩn chảy: von Mises dự đoán vật liệu chảy khi σvM = √(σ₁² − σ₁σ₂ + σ₂²) đạt σY; Tresca khi ứng suất cắt lớn nhất tuyệt đối đạt σY/2.

στxyCσ1​σ2​

Ứng suất chính

C = 90 MPa, R = 53,85 MPa
σ1 = 143,9 MPa, σ2 = 36,15 MPa, σ3 = 0

Tresca

σmax − σmin = 143,9 MPa
n = 275/143,9 = 1,912

Đáp số

n = 1,912

n > 1 nghĩa là điểm không bị chảy; Tresca thận trọng hơn von Mises một chút.

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

Đáp số: δmax = 1,106 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 = 10 N/mm, L = 2500 mm
EI = 4,6 × 10¹² N·mm²

Độ võng

δmax = 5wL⁴/(384EI) = 5·10·2500⁴/(384·4,6 × 10¹²) = 1,106 mm

Đáp số

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

Đáp số: Ireq = 4,185 × 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 = 2500/360 = 6,944 mm

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

I = 5wL⁴/(384·E·δallow) = 5·4·2500⁴/(384·70000·6,944) = 4,185 × 10⁶ mm⁴
Ireq = 4,185 × 10⁶ mm⁴

Đáp số

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

Đáp số: Pcr = 177,2 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 = π·80⁴/64 = 2,011 × 10⁶ mm⁴

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

KL = 0,7·4000 = 2800 mm

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

Pcr = π²EI/(KL)² = π²·70000·2,011 × 10⁶/2800² = 177179 N = 177,2 kN

Đáp số

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

Đáp số: σcr = 1,079 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 = π·50⁴/64 = 306796 mm⁴

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

A = 1963 mm²
r = √(I/A) = √(306796/1963) = 12,5 mm
KL = 2·5000 = 10000 mm
KL/r = 10000/12,5 = 800

Ứng suất tới hạn

σcr = π²E/(KL/r)² = π²·70000/800² = 1,079 MPa

Đáp số

σcr = 1,079 MPa

ExamProblems · examproblems.com/courses/mechanics-of-materials