การออกแบบเครื่องจักรกล: ชุดสรุปเร่งรัด

สรุปสูตร ขั้นตอนวิธีทำที่สำคัญของทุกหัวข้อ และโจทย์ฝึก 12 ข้อ พร้อมเฉลยละเอียดไว้ท้ายเล่ม

โจทย์ใหม่กลับไปที่รายวิชา

ในหน้าต่างพิมพ์ ให้เลือก “บันทึกเป็น PDF” เพื่อเก็บสำเนาไว้ในอุปกรณ์ของคุณ

1สูตรและวิธีทำ

เฟืองและชุดเฟือง

เฟืองขบกัน
n1·N1 = n2·N2
แรงบิด
T2 = T1·(N2/N1)·η
กำลัง
P = T·ω, ω = 2π·n/60
ชุดเฟืองประกอบ
อัตราทดรวม = ผลคูณของอัตราทดแต่ละขั้น

การบิดของเพลา

ความเค้นเฉือน
τ = T·r/J
เพลาตัน
J = πd⁴/32
เพลากลวง
J = π(D⁴ − d⁴)/32
มุมบิด
φ = TL/(GJ)
กำลัง
P = T·ω = 2πnT/60

วิธีการ

  1. หาโมเมนต์ความเฉื่อยเชิงขั้ว: J = πd⁴/32 สำหรับเพลาตัน
  2. ความเค้นเฉือนที่ผิว: τ = T·r/J โดย r คือรัศมีภายนอก
  3. มุมบิด: φ = TL/(GJ) เป็นเรเดียน
  4. ถ้าโจทย์ให้กำลัง ให้หาแรงบิดก่อนด้วย P = T·ω, ω = 2πn/60

วงกลมของมอร์และการแปลงความเค้น

ศูนย์กลาง
C = (σx + σy)/2
รัศมี
R = √(((σx − σy)/2)² + τxy²)
ความเค้นหลัก
σ1,2 = C ± R
ทิศทางหลัก
tan 2θp = 2τxy/(σx − σy)
พื้นเอียง
σθ = C + (σx − σy)/2·cos 2θ + τxy·sin 2θ

วิธีการ

  1. คำนวณศูนย์กลางของวงกลมมอร์: C = (σx + σy)/2
  2. คำนวณรัศมี: R = √(((σx − σy)/2)² + τxy²)
  3. ความเค้นหลักคือ C ± R และความเค้นเฉือนในระนาบที่มากที่สุดคือ R
  4. มุมบนวงกลมเป็นสองเท่าของมุมบนชิ้นส่วน: tan 2θp = 2τxy/(σx − σy)

ความเค้นและความเครียดตามแนวแกน

ความเค้นแนวฉาก
σ = F/A
ความเครียด
ε = ΔL/L
กฎของฮุก
σ = E·ε
ความยืด
ΔL = FL/(EA)
ความเครียดจากความร้อน
εT = α·ΔT

การโก่งเดาะของเสา

น้ำหนักวิกฤตของออยเลอร์
Pcr = π²EI/(KL)²
ตัวประกอบความยาวประสิทธิผล
K = 1 หมุด–หมุด, 0.5 ยึดแน่น–ยึดแน่น, 0.7 ยึดแน่น–หมุด, 2 ยึดแน่น–อิสระ
ความเค้นวิกฤต
σcr = Pcr/A
ความเพรียว
λ = KL/r, r = √(I/A)

การสั่นสะเทือนเชิงกล

ความถี่ธรรมชาติ
ωn = √(k/m), fn = ωn/(2π)
จากระยะแอ่นสถิต
ωn = √(g/δst)
การหน่วงวิกฤต
cc = 2√(km) = 2mωn
อัตราส่วนการหน่วง
ζ = c/cc
ความถี่ที่ถูกหน่วง
ωd = ωn√(1 − ζ²)
ลอการิทึมดีเครเมนต์
δ = (1/n)·ln(x0/xn), ζ = δ/√(4π² + δ²)
สปริงต่อขนาน
k = k1 + k2
สปริงต่ออนุกรม
1/k = 1/k1 + 1/k2

2โจทย์ฝึก

1.เฟืองและชุดเฟืองปานกลาง
N1​ = 16N2​ = 25

มอเตอร์ 3 kW หมุนที่ 1440 rpm ขับเฟืองเล็กที่มี 16 ฟัน ซึ่งขบกับเฟืองที่มี 25 ฟัน คู่เฟืองมีประสิทธิภาพ 98% หาแรงบิดขาออก T2 บนเพลาของเฟือง

P = 3 kWn1 = 1440 rpmN1 = 16N2 = 25η = 98%
T2 = N·m
2.เฟืองและชุดเฟืองยาก
A = 20B = 36C = 36D = 40

ในชุดเฟืองประกอบ เฟือง A (20 ฟัน) บนเพลาขาเข้าหมุนที่ 1450 rpm ขับเฟือง B (36 ฟัน) เฟือง C (36 ฟัน) ติดอยู่กับเพลาเดียวกับ B และขับเฟือง D (40 ฟัน) บนเพลาขาออก หาความเร็วรอบขาออก nout

nin = 1450 rpmA = 20B = 36C = 36D = 40
nout = rpm
3.การบิดของเพลาปานกลาง
T = 750 N·md = 30 mmL = 1 mφ = ?

เพลาตันอะลูมิเนียมถูกยึดติดกับผนังและถูกบิดด้วยแรงบิด T ที่ปลายอิสระ หามุมบิด φ ของปลายอิสระ ในหน่วยองศา

T = 750 N·md = 30 mmL = 1 mG = 26 GPa
φ = °
4.การบิดของเพลายาก
d = ?20 kW, 1450 rpm

เพลา อะลูมิเนียม ตันส่งกำลัง 20 kW ที่ 1450 rpm ความเค้นเฉือนต้องไม่เกิน 75 MPa และมุมบิดต้องไม่เกิน 0.5° ต่อความยาว m จงหาเส้นผ่านศูนย์กลางที่น้อยที่สุด d ที่ตรงตามทั้งสองขีดจำกัด

P = 20 kWN = 1450 rpmτallow = 75 MPaθallow = 0.5°/mG = 26 GPa
d = mm
5.วงกลมของมอร์และการแปลงความเค้นปานกลาง
σx​ = −10σy​ = 60τxy​ = 10MPa

จุดหนึ่งในชิ้นส่วนที่รับน้ำหนักมี σₓ = −10 MPa, σᵧ = 60 MPa และ τₓᵧ = 10 MPa (ความเค้นดึงเป็นบวก) หามุม θp จากแกน x ไปยังทิศทางของความเค้นหลักที่มากที่สุด ระหว่าง −90° ถึง 90°

σx = −10 MPaσy = 60 MPaτxy = 10 MPa
θp = °
6.วงกลมของมอร์และการแปลงความเค้นยาก
60°D = 800 mm, t = 12 mm, p = 1 MPap

ถังทรงกระบอกผนังบางเส้นผ่านศูนย์กลางภายใน 800 mm ความหนาผนัง 12 mm บรรจุแก๊สที่ความดันเกจ 1 MPa ทำจากแผ่นเชื่อมตามรอยเชื่อมเกลียวที่ทำมุม 60° กับแกนของทรงกระบอก จงหาขนาดของความเค้นเฉือน |τn| ที่กระทำตามรอยเชื่อม

D = 800 mmt = 12 mmp = 1 MPaα = 60°
|τn| = MPa
7.ความเค้นและความเครียดตามแนวแกนปานกลาง
P = ?12500 mm500 mmØ40Ø30δ = 0.2 mm

แท่งหนึ่งประกอบด้วยท่อนเหล็กกล้า (เส้นผ่านศูนย์กลาง 40 mm, ยาว 500 mm, E = 200 GPa) ต่อปลายชนปลายกับท่อนที่ทำจากบรอนซ์ (เส้นผ่านศูนย์กลาง 30 mm, ยาว 500 mm, E = 110 GPa) ยึดปลายหนึ่งและดึงอีกปลาย จงหาโหลดตามแนวแกน P ที่ทำให้มันยืดทั้งหมด 0.2 mm

ท่อนที่ 1: เหล็กกล้า, Ø40 mm, 500 mm, E = 200 GPaท่อนที่ 2: บรอนซ์, Ø30 mm, 500 mm, E = 110 GPaδ = 0.2 mm
P = kN
8.ความเค้นและความเครียดตามแนวแกนยาก
Ø 35 mmทองเหลือง2 mØ 20 mmอะลูมิเนียม2 mP = 15 kN

แท่งขั้นบันไดถูกยึดติดกับผนังและรับแรงที่ปลายอิสระ หาความเค้นแนวฉากสูงสุด σmax ในแท่ง

P = 15 kNช่วงที่ 1: Ø35 mm, 2 m, ทองเหลือง (E = 100 GPa)ช่วงที่ 2: Ø20 mm, 2 m, อะลูมิเนียม (E = 70 GPa)
σmax = MPa
9.การโก่งเดาะของเสาปานกลาง
PL = 1 md = 20 mmหน้าตัด

เสายาว L ทำจากเหล็กกล้า โดยยึดด้วยหมุดทั้งสองปลาย ใช้สูตรของออยเลอร์หาน้ำหนักวิกฤตที่ทำให้โก่งเดาะ Pcr

L = 1 mE = 200 GPa (เหล็กกล้า)d = 20 mm
Pcr = kN
10.การโก่งเดาะของเสายาก
P = 5 kNL = 1 md = ?หน้าตัด

เสากลมตันยาว L ทำจากเหล็กกล้า โดยยึดแน่นที่ฐานและยึดด้วยหมุดที่ปลายบน ต้องรับ 5 kN ด้วยตัวประกอบความปลอดภัย 2.5 ต่อการโก่งเดาะ หาเส้นผ่านศูนย์กลางที่เล็กที่สุด d

P = 5 kNn = 2.5L = 1 mE = 200 GPa (เหล็กกล้า)
d = mm
11.การสั่นสะเทือนเชิงกลปานกลาง
k = 15 kN/mc = 310 N·s/mm = 10 kg

มวล 10 kg ถูกรองรับด้วยสปริงที่มีค่าความแข็ง 15 kN/m และตัวหน่วงหนืดที่มีสัมประสิทธิ์การหน่วง 310 N·s/m หาอัตราส่วนการหน่วง ζ

m = 10 kgk = 15 kN/mc = 310 N·s/m
ζ =
12.การสั่นสะเทือนเชิงกลยาก
k1​ = 15k2​ = 5m = 10 kg

มวล 10 kg แขวนจากสปริงสองตัวที่วางเคียงกัน (ขนาน) ที่มีค่าความแข็ง k₁ = 15 kN/m และ k₂ = 5 kN/m มวลอยู่ในแนวระดับ หาค่าความแข็งสมมูล keq ของระบบสปริงใน kN/m

m = 10 kgk1 = 15 kN/mk2 = 5 kN/m
keq = kN/m

3สารละลาย

1.เฟืองและชุดเฟือง

คำตอบ: T2 = 30.46 N·m

ชุดเฟือง: เฟืองที่ขบกันมีความเร็วที่ฟันเท่ากัน ดังนั้น n₁·N₁ = n₂·N₂ เฟืองที่มีฟันน้อยกว่าหมุนเร็วกว่า เฟืองบนเพลาเดียวกันหมุนด้วยความเร็วเท่ากัน กำลังถูกส่งต่อ (หักการสูญเสีย) แรงบิดจึงเพิ่มเมื่อความเร็วลดลง: T₂ = T₁·(N₂/N₁)·η และ P = T·ω โดย ω = 2π·n/60

แรงบิดของมอเตอร์

ω₁ = 2π·1440/60 = 150.8 rad/s
T1 = P/ω₁ = 3000/150.8 = 19.89 N·m

แรงบิดขาออก

T2 = 19.89·25/16·0.98 = 30.46 N·m

คำตอบ

T2 = 30.46 N·m
2.เฟืองและชุดเฟือง

คำตอบ: nout = 725 rpm

ชุดเฟือง: เฟืองที่ขบกันมีความเร็วที่ฟันเท่ากัน ดังนั้น n₁·N₁ = n₂·N₂ เฟืองที่มีฟันน้อยกว่าหมุนเร็วกว่า เฟืองบนเพลาเดียวกันหมุนด้วยความเร็วเท่ากัน กำลังถูกส่งต่อ (หักการสูญเสีย) แรงบิดจึงเพิ่มเมื่อความเร็วลดลง: T₂ = T₁·(N₂/N₁)·η และ P = T·ω โดย ω = 2π·n/60

ขั้นแรก

nB = 1450·20/36 = 805.6 rpm = nC

ขั้นที่สอง

nout = 805.6·36/40 = 725 rpm

คำตอบ

nout = 725 rpm

อัตราทดรวม = 2 : 1

3.การบิดของเพลา

คำตอบ: φ = 20.78 °

การบิดของเพลากลม: ความเค้นเฉือนเพิ่มขึ้นเป็นเส้นตรงจากศูนย์ที่ศูนย์กลางถึง τmax = T·r/J ที่ผิว โดย J = πd⁴/32 สำหรับเพลาตัน และ J = π(do⁴ − di⁴)/32 สำหรับเพลากลวง มุมบิดคือ φ = TL/(GJ) เป็นเรเดียน กำลังและแรงบิดสัมพันธ์กันด้วย P = T·ω

โมเมนต์ความเฉื่อยเชิงขั้วของพื้นที่

J = π·30⁴/32 = 79522 mm⁴

มุมบิด

φ = TL/(GJ) = 750·1000·1000/(26000·79522) = 0.3627 rad
φ = 0.3627·180/π = 20.78°

คำตอบ

φ = 20.78°
4.การบิดของเพลา

คำตอบ: d = 49.31 mm

การบิดของเพลากลม: ความเค้นเฉือนเพิ่มขึ้นเป็นเส้นตรงจากศูนย์ที่ศูนย์กลางถึง τmax = T·r/J ที่ผิว โดย J = πd⁴/32 สำหรับเพลาตัน และ J = π(do⁴ − di⁴)/32 สำหรับเพลากลวง มุมบิดคือ φ = TL/(GJ) เป็นเรเดียน กำลังและแรงบิดสัมพันธ์กันด้วย P = T·ω

แรงบิด

ω = 2πN/60 = 2π·1450/60 = 151.8 rad/s
T = P/ω = 20000/151.8 = 131.7 N·m

เส้นผ่านศูนย์กลางจากความแข็งแรง

T = 131.7 N·m = 131714 N·mm
dτ = ∛(16·131714/(π·75)) = 20.76 mm

เส้นผ่านศูนย์กลางจากความแข็งเกร็ง

θ = 0.5·π/180/1000 = 8.727 × 10⁻⁶ rad/mm
G = 26 GPa = 26000 MPa
dφ = (32·131714/(π·26000·8.727 × 10⁻⁶))1/4 = 49.31 mm

เส้นผ่านศูนย์กลางที่เป็นตัวกำหนด

ทั้งสองขีดจำกัดต้องเป็นจริง ดังนั้นเส้นผ่านศูนย์กลางที่มากกว่าเป็นตัวกำหนด

d = max(20.76, 49.31) = 49.31 mm

คำตอบ

d = 49.31 mm

ในที่นี้ความแข็งเกร็งเป็นตัวกำหนด

5.วงกลมของมอร์และการแปลงความเค้น

คำตอบ: θp = 82.03 °

การแปลงความเค้นระนาบ: บนระนาบที่หมุนไป θ σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ วงกลมของมอร์แสดงทุกระนาบพร้อมกัน ศูนย์กลางอยู่ที่ C = (σₓ + σᵧ)/2 และรัศมีคือ R = √(((σₓ − σᵧ)/2)² + τₓᵧ²) ความเค้นหลักคือ C ± R และความเค้นเฉือนในระนาบที่มากที่สุดคือ R

στxyCσ1​σ2​

มุมบนวงกลม

จุด x อยู่ที่มุม 2θp จากแกน σ บนวงกลมของมอร์ ควอดรันต์ได้จากเครื่องหมายของ σₓ − σᵧ และ τₓᵧ

2θp = atan2(2·10, −10 − 60) = 164.1°

มุมบนชิ้นส่วน

θp = 82.03°

คำตอบ

θp = 82.03°

มุมบวกคือทวนเข็มนาฬิกา

6.วงกลมของมอร์และการแปลงความเค้น

คำตอบ: |τn| = 7.217 MPa

การแปลงความเค้นระนาบ: บนระนาบที่หมุนไป θ σ = (σₓ + σᵧ)/2 + (σₓ − σᵧ)/2·cos 2θ + τₓᵧ·sin 2θ วงกลมของมอร์แสดงทุกระนาบพร้อมกัน ศูนย์กลางอยู่ที่ C = (σₓ + σᵧ)/2 และรัศมีคือ R = √(((σₓ − σᵧ)/2)² + τₓᵧ²) ความเค้นหลักคือ C ± R และความเค้นเฉือนในระนาบที่มากที่สุดคือ R ทรงกระบอกผนังบาง: ความเค้นตามเส้นรอบวง σh = pr/t และความเค้นตามแนวแกน σa = pr/(2t) ทั้งสองคือความเค้นหลักในผนัง

ความเค้นในผนัง

r = 800/2 = 400 mm
σh = pr/t = 1·400/12 = 33.33 MPa
σa = pr/(2t) = 16.67 MPa

แปลงไปยังรอยเชื่อม

ให้ x ตามแนวแกน (σa) และ y รอบเส้นรอบวง (σh) เส้นแนวฉากของรอยเชื่อมทำมุม 90° − α กับแกน

|τn| = (33.33 − 16.67)·sin 60°·cos 60° = 7.217 MPa

คำตอบ

|τn| = 7.217 MPa
7.ความเค้นและความเครียดตามแนวแกน

คำตอบ: P = 23.75 kN

วัสดุสองชนิด: แบบอนุกรม (ต่อกันทีละส่วน) ทุกส่วนรับแรงเท่ากันและการยืดตัวบวกกัน: δ = P·Σ L/(E·A) แบบขนาน (วางเคียงกันระหว่างแผ่นแข็ง) ทั้งสองส่วนหดเท่ากัน ดังนั้นความเครียดเท่ากัน: ส่วนที่แข็งกว่ารับสัดส่วนมากกว่า Pi = P·EiAi/ΣEA และ σi = Ei·ε แท่งที่ยึดระหว่างผนังแข็งสองด้านโดยมีโหลดที่รอยต่อเป็นโครงสร้างอินดีเทอร์มิเนต: สองส่วนแบ่งโหลดตามสัดส่วนความแข็ง k = EA/L

พื้นที่

A₁ = 1257 mm², A₂ = 706.9 mm²

ความอ่อนตัว

L₁/(E₁A₁) = 500/(200000·1257) = 1.989 × 10⁻⁶ mm/N
L₂/(E₂A₂) = 500/(110000·706.9) = 6.431 × 10⁻⁶ mm/N

โหลด

P = 0.2/8.42 × 10⁻⁶ = 23753 N = 23.75 kN

คำตอบ

P = 23.75 kN
8.ความเค้นและความเครียดตามแนวแกน

คำตอบ: σmax = 47.75 MPa

ชิ้นส่วนรับแรงตามแนวแกน: แต่ละช่วงที่มีหน้าตัดคงที่เป็นไปตาม δ = PL/(AE) สำหรับชิ้นส่วนที่ต่ออนุกรมกัน ความยืดจะบวกกัน

แรงเท่ากันทั้งสองช่วง

ปลายอิสระรับแรง P ดังนั้นทุกส่วนของแท่งมีแรงตามแนวแกน P เท่ากัน มีเพียงพื้นที่และวัสดุที่เปลี่ยนไป

A1 = π·35²/4 = 962.1 mm², A2 = π·20²/4 = 314.2 mm²
σ1 = 15.59 MPa
σ2 = 47.75 MPa

ความเค้นสูงสุด

ช่วงที่บางกว่ามีพื้นที่น้อยกว่า จึงรับความเค้นมากกว่า

σmax = 47.75 MPa

คำตอบ

σmax = 47.75 MPa
9.การโก่งเดาะของเสา

คำตอบ: Pcr = 15.5 kN

การโก่งเดาะของออยเลอร์: เสาเพรียวจะไม่เสถียรที่ Pcr = π²EI/(KL)² โดย KL คือความยาวประสิทธิผล: K = 1 สำหรับหมุด–หมุด, 2 สำหรับยึดแน่น–อิสระ, 0.7 สำหรับยึดแน่น–หมุด และ 0.5 สำหรับยึดแน่น–ยึดแน่น เสาจะโก่งเดาะรอบแกนที่มี I น้อยที่สุด สูตรนี้ใช้ได้กับเสาเพรียวซึ่งความเค้นวิกฤตยังต่ำกว่าความเค้นครากของวัสดุ

โมเมนต์ความเฉื่อยของพื้นที่ที่น้อยที่สุด

หน้าตัดกลมมี I เท่ากันรอบทุกแกน

I = πd⁴/64 = π·20⁴/64 = 7854 mm⁴

ความยาวประสิทธิผล

KL = 1·1000 = 1000 mm

น้ำหนักวิกฤต

Pcr = π²EI/(KL)² = π²·200000·7854/1000² = 15503 N = 15.5 kN

คำตอบ

Pcr = 15.5 kN
10.การโก่งเดาะของเสา

คำตอบ: d = 15.86 mm

การโก่งเดาะของออยเลอร์: เสาเพรียวจะไม่เสถียรที่ Pcr = π²EI/(KL)² โดย KL คือความยาวประสิทธิผล: K = 1 สำหรับหมุด–หมุด, 2 สำหรับยึดแน่น–อิสระ, 0.7 สำหรับยึดแน่น–หมุด และ 0.5 สำหรับยึดแน่น–ยึดแน่น เสาจะโก่งเดาะรอบแกนที่มี I น้อยที่สุด สูตรนี้ใช้ได้กับเสาเพรียวซึ่งความเค้นวิกฤตยังต่ำกว่าความเค้นครากของวัสดุ

น้ำหนักวิกฤตที่ต้องการ

Pcr = n·P = 2.5·5000 = 12500 N
KL = 0.7·1000 = 700 mm

แก้หาเส้นผ่านศูนย์กลาง

d = (64·n·P·(KL)²/(π³E))^(1/4) = (64·12500·700²/(π³·200000))^(1/4) = 15.86 mm

คำตอบ

d = 15.86 mm

ตรวจสอบ: ค่าความเพรียว KL/r คือ 176.6 ดังนั้นเสาเป็นเสาเพรียวและใช้สูตรออยเลอร์ได้

11.การสั่นสะเทือนเชิงกล

คำตอบ: ζ = 0.4002

การสั่นอิสระของระบบมวล–สปริง: ความถี่ธรรมชาติแบบไม่มีการหน่วงคือ ωn = √(k/m) ใน rad/s และ fn = ωn/(2π) ใน Hz ตัวหน่วงหนืด c ให้อัตราส่วนการหน่วง ζ = c/cc โดยการหน่วงวิกฤต cc = 2·√(k·m) = 2·m·ωn ระบบที่หน่วงน้อยกว่าวิกฤต (ζ < 1) สั่นที่ ωd = ωn·√(1 − ζ²) สปริงต่อขนานบวกกัน k = k₁ + k₂ สปริงต่ออนุกรมบวกกันแบบส่วนกลับ 1/k = 1/k₁ + 1/k₂

การหน่วงวิกฤต

cc = 2·√(k·m) = 2·√(15000·10) = 774.6 N·s/m

อัตราส่วนการหน่วง

ζ = c/cc = 310/774.6 = 0.4002

คำตอบ

ζ = 0.4002

ζ < 1 ระบบจึงหน่วงน้อยกว่าวิกฤตและสั่นขณะที่ลดลง

12.การสั่นสะเทือนเชิงกล

คำตอบ: keq = 20 kN/m

การสั่นอิสระของระบบมวล–สปริง: ความถี่ธรรมชาติแบบไม่มีการหน่วงคือ ωn = √(k/m) ใน rad/s และ fn = ωn/(2π) ใน Hz ตัวหน่วงหนืด c ให้อัตราส่วนการหน่วง ζ = c/cc โดยการหน่วงวิกฤต cc = 2·√(k·m) = 2·m·ωn ระบบที่หน่วงน้อยกว่าวิกฤต (ζ < 1) สั่นที่ ωd = ωn·√(1 − ζ²) สปริงต่อขนานบวกกัน k = k₁ + k₂ สปริงต่ออนุกรมบวกกันแบบส่วนกลับ 1/k = 1/k₁ + 1/k₂

ค่าความแข็งสมมูล

สปริงที่วางเคียงกันมีระยะแอ่นเท่ากันและแรงบวกกัน ค่าความแข็งจึงบวกกันตรง ๆ: k = k1 + k2

keq = 15 + 5 = 20 kN/m

คำตอบ

keq = 20 kN/m

ExamProblems · examproblems.com/courses/machine-design