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
Ø 16 mmThép3 mP = 90 kN

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

P = 90 kNd = 16 mmL = 3 mE = 200 GPa
δ = mm
2.Ứng suất và biến dạng dọc trụcKhó
Ø 25 mmĐồng thau1,5 mØ 16 mmThép1 mP = 55 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 = 55 kNđoạn 1: Ø25 mm, 1,5 m, đồng thau (E = 100 GPa)đoạn 2: Ø16 mm, 1 m, thép (E = 200 GPa)
σmax = MPa
3.Xoắn của trụcTrung bình
do​ = 30, di​ = 18 mmT = 200 N·mL = 3 m

Một trục rỗng nhôm (đường kính ngoài 30 mm, đường kính trong 18 mm, chiều dài 3 m) được ngàm một đầu và bị xoắn bởi mô-men xoắn 200 N·m ở đầu kia. Tìm góc xoắn φ theo độ.

T = 200 N·mdo = 30 mmdi = 18 mmL = 3 mG = 26 GPa
φ = °
4.Xoắn của trụcKhó
d = ?10 kW, 300 rpm

Một trục đặc nhôm truyền 10 kW ở 300 rpm. Ứng suất cắt không được vượt quá 75 MPa và góc xoắn không được vượt quá 0,5° trên mỗi m chiều dài. Tìm đường kính nhỏ nhất d thỏa mãn cả hai giới hạn.

P = 10 kWN = 300 rpmτallow = 75 MPaθallow = 0,5°/mG = 26 GPa
d = mm
5.Lực cắt và mômen uốnTrung bình
w = 8 kN/mABmặt cắt2,5 m1,5 m

Tìm mômen uốn Mx tại mặt cắt x = 2,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 mw = 8 kN/m từ 0 đến 4 mx = 2,5 m
Mx = kN·m
6.Lực cắt và mômen uốnKhó
w = 12 kN/mP = 45 kNABmặt cắt7 m3 m

Tìm mômen uốn Mx tại mặt cắt x = 7 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 = 45 kN tại 7 mw = 12 kN/m từ 0 đến 10 mx = 7 m
Mx = kN·m
7.Trọng tâm và mômen quán tính mặt cắtTrung bình
3001516020xkích thước tính bằng mm

Tìm mômen quán tính Ix của mặt cắt chữ T đối với trục ngang x đi qua trọng tâm.

bf = 300 mmtf = 15 mmtw = 20 mmhw = 160 mm
Ix = × 10⁶ mm⁴
8.Trọng tâm và mômen quán tính mặt cắtKhó
120 mm150 mm60 mm80 × 50xIx​ = ?

Một tấm 120 mm × 150 mm có một lỗ khoét hình chữ nhật rộng 80 mm và cao 50 mm, đặt giữa đường tâm thẳng đứng với tâm cách đáy 60 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 = 120 mmh = 150 mmlỗ khoét: 80 × 50 mmtâm lỗ: y = 60 mm
Ix = × 10⁶ mm⁴
9.Vòng tròn Mohr và biến đổi ứng suấtTrung bình
σx​ = −30σy​ = 100τxy​ = 40MPa

Một điểm trong chi tiết chịu tải có σₓ = −30 MPa, σᵧ = 100 MPa và τₓᵧ = 40 MPa (kéo là dương). Tìm góc θp từ trục x đến phương của ứng suất chính lớn nhất, trong khoảng −90° đến 90°.

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

Một điểm trong chi tiết chịu tải có σₓ = −20 MPa, σᵧ = 80 MPa và τₓᵧ = 40 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 25° ngược chiều kim đồng hồ so với trục x (dương cùng chiều với τₓᵧ).

σx = −20 MPaσy = 80 MPaτxy = 40 MPaθ = 25°
τθ = MPa
11.Độ võng của dầmTrung bình
w = 4 kN/mδL = 5 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 = 4 kN/mL = 5 mE = 200 GPa (thép)I = 23 × 10⁶ mm⁴
δ = mm
12.Độ võng của dầmKhó
w = 20 kN/mδL = 1,5 m

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

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

Một cột dài L làm bằng thép đượ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 = 3 mE = 200 GPa (thép)b × h = 25 × 40 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à khớp ở đỉ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ố: δ = 6,714 mm

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

Ứng suất

A = π·16²/4 = 201,1 mm²
σ = P/A = 90·1000/201,1 = 447,6 MPa

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

E for thép = 200 GPa = 200000 MPa.

ε = σ/E = 447,6/200000 = 0,002238 = 2238 µε

Độ giãn dài

δ = ε·L = 0,002238·3000 mm = 6,714 mm

Đáp số

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

Đáp số: σmax = 273,5 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 = π·25²/4 = 490,9 mm², A2 = π·16²/4 = 201,1 mm²
σ1 = 112 MPa
σ2 = 273,5 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 = 273,5 MPa

Đáp số

σmax = 273,5 MPa
3.Xoắn của trục

Đáp số: φ = 19,1 °

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⁴ − 18⁴)/32 = 69216 mm⁴

Góc xoắn

T = 200 N·m = 200000 N·mm
L = 3 m = 3000 mm
G = 26 GPa = 26000 MPa
φ = TL/(GJ) = 200000·3000/(26000·69216) = 0,3334 rad
φ = 0,3334·180/π = 19,1°

Đáp số

φ = 19,1°
4.Xoắn của trục

Đáp số: d = 61,48 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 theo điều kiện bền

T = 318,3 N·m = 318310 N·mm
dτ = ∛(16·318310/(π·75)) = 27,86 mm

Đường kính theo điều kiện cứng

θ = 0,5·π/180/1000 = 8,727 × 10⁻⁶ rad/mm
G = 26 GPa = 26000 MPa
dφ = (32·318310/(π·26000·8,727 × 10⁻⁶))1/4 = 61,48 mm

Đường kính quyết định

Cả hai giới hạn đều phải thỏa, nên đường kính lớn hơn quyết định.

d = max(27,86, 61,48) = 61,48 mm

Đáp số

d = 61,48 mm

Ở đây điều kiện cứng quyết định.

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

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

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

Phản lực gối

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

Hợp lực của tải phân bố đều: W = w·ℓ = 8·4 = 32 kN, đặt tại x̄ = 2 m
ΣMA = 0: RB·4 = 8·4·2
RB = 64/4 = 16 kN
ΣFy = 0: RA = 32 − 16 = 16 kN
w = 8 kN/mABRA​ = 16RB​ = 16mặt cắt2,5 m1,5 m

Cắt tại x = 2,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 = 16·2,5 − 8·2,5·1,25 = 15 kN·m

Đáp số

Mx = 15 kN·m
6.Lực cắt và mômen uốn

Đáp số: Mx = 220,5 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·10 = 120 kN, đặt tại x̄ = 5 m
ΣMA = 0: RB·10 = 45·7 + 12·10·5
RB = 915/10 = 91,5 kN
ΣFy = 0: RA = 165 − 91,5 = 73,5 kN
w = 12 kN/mP = 45 kNABRA​ = 73,5RB​ = 91,5mặt cắt7 m3 m

Cắt tại x = 7 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 = 73,5·7 − 12·7·3,5 = 220,5 kN·m

Đáp số

Mx = 220,5 kN·m
7.Trọng tâm và mômen quán tính mặt cắt

Đáp số: Ix = 21,23 × 10⁶ mm⁴

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.

Chia thành hai hình chữ nhật

Đo mọi y từ mép dưới.

cánh: Af = 300·15 = 4500 mm², yf = 160 + 15/2 = 167,5 mm
bụng: Aw = 20·160 = 3200 mm², yw = 160/2 = 80 mm

Trọng tâm

ȳ = (Afyf + Awyw)/(Af + Aw) = (4500·167,5 + 3200·80)/7700 = 131,1 mm

Định lý trục song song

cánh: 300·15³/12 + 4500·(167,5 − 131,1)² = 6,035 × 10⁶ mm⁴
bụng: 20·160³/12 + 3200·(80 − 131,1)² = 1,519 × 10⁷ mm⁴
Ix = 6,035 × 10⁶ + 1,519 × 10⁷ = 2,123 × 10⁷ mm⁴ = 21,23 × 10⁶ mm⁴

Đáp số

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

Đáp số: Ix = 31,76 × 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

ȳ = (18000·75 − 4000·60)/(18000 − 4000) = 79,29 mm

Tấm đầy đủ quanh ȳ

120·150³/12 + 18000·−4,286² = 34,08 × 10⁶ mm⁴

Lỗ quanh ȳ

80·50³/12 + 4000·−19,29² = 2,321 × 10⁶ mm⁴

Trừ

Ix = 31,76 × 10⁶ mm⁴

Đáp số

Ix = 31,76 × 10⁶ mm⁴
9.Vòng tròn Mohr và biến đổi ứng suất

Đáp số: θp = 74,2 °

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σ1​σ2​

Góc trên vòng tròn

Điểm x nằm ở góc 2θp so với trục σ trên vòng tròn Mohr; góc phần tư suy ra từ dấu của σₓ − σᵧ và τₓᵧ.

2θp = atan2(2·40, −30 − 100) = 148,4°

Góc trên phần tử

θp = 74,2°

Đáp số

θp = 74,2°

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

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

Đáp số: τθ = 64,01 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

τθ = −−50·0,766 + 40·0,6428 = 64,01 MPa

Đáp số

τθ = 64,01 MPa
11.Độ võng của dầm

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

Độ võng

δ = wL⁴/(8EI) = 4·5000⁴/(8·4,6 × 10¹²) = 67,93 mm

Đáp số

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

Đáp số: Ireq = 4,52 × 10⁶ mm⁴

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

Giới hạn độ võng

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

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

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

Đáp số

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

Đáp số: Pcr = 11,42 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 = 40·25³/12 = 52083 mm⁴

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

KL = 1·3000 = 3000 mm

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

Pcr = π²EI/(KL)² = π²·200000·52083/3000² = 11423 N = 11,42 kN

Đáp số

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

Đáp số: σcr = 8,812 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 = 0,7·5000 = 3500 mm
KL/r = 3500/12,5 = 280

Ứng suất tới hạn

σcr = π²E/(KL/r)² = π²·70000/280² = 8,812 MPa

Đáp số

σcr = 8,812 MPa

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