Tổng quan nghiên cứu

Trong kỹ thuật cơ khí ô tô, cơ cấu piston - trục khuỷu - thanh truyền đóng vai trò trung tâm trong việc chuyển đổi năng lượng nhiệt thành cơ năng quay để dẫn động phương tiện. Nghiên cứu tập trung vào động cơ xăng B6 dung tích 1,6 lít (1598 cc) trang bị trên dòng xe Mazda 323, sở hữu công suất cực đại 62 kW tại 5200 vòng/phút và mô-men xoắn cực đại 110 N.m tại 3600 vòng/phút. Vấn đề cốt lõi đặt ra là cơ cấu này liên tục vận hành dưới tải trọng cơ học và tải trọng nhiệt khắc nghiệt, chịu áp suất khí cháy đỉnh lên tới 5,3 MN/m² kết hợp cùng lực quán tính chuyển động tịnh tiến và chuyển động quay tuần hoàn. Tình trạng này dễ dẫn tới biến dạng hình học, tập trung ứng suất uốn - xoắn phức tạp và làm suy giảm tuổi thọ làm việc của cụm chi tiết.

Mục tiêu cụ thể của nghiên cứu là xây dựng quy trình tính toán nhiệt động lực học hoàn chỉnh, xác định chính xác các vector tải trọng tác dụng, thiết kế mô hình không gian 3D tham số hóa và ứng dụng phương pháp phần tử hữu hạn (FEM) trên phần mềm CATIA V5R21 để phân tích biến dạng cũng như phân bố ứng suất của trục khuỷu. Đề tài được triển khai thực nghiệm tính toán tại Viện Cơ khí thuộc Trường Đại học Giao thông Vận tải Thành phố Hồ Chí Minh trong giai đoạn năm 2023. Ý nghĩa thực tiễn của công trình thể hiện qua việc tối ưu hóa độ bền kết cấu, giảm thiểu từ 15% đến 20% nguy cơ phá hủy mỏi và rút ngắn khoảng 30% chu kỳ thiết kế - kiểm nghiệm so với các phương pháp tính toán giải tích truyền thống.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu vận dụng đồng thời hai nền tảng lý thuyết chủ đạo: lý thuyết nhiệt động lực học động cơ đốt trong và cơ học biến dạng vật rắn thông qua phương pháp phần tử hữu hạn. Quá trình nhiệt động được mô hình hóa dựa trên chu trình công tác của động cơ 4 kỳ nạp khí tự nhiên với chỉ số nén đa biến n1 bằng 1,35 và chỉ số giãn nở đa biến n2 bằng 1,25. Phương pháp phân tích động lực học phẳng xác định sự biến thiên của lực khí thể, lực quán tính chuyển động tịnh tiến và lực tiếp tuyến - pháp tuyến tác dụng lên cơ cấu theo góc quay trục khuỷu từ 0 độ đến 720 độ.

Về mặt cơ học kết cấu, nguyên lý cực tiểu hóa thế năng toàn phần và tiêu chuẩn ứng suất tương đương Von-Mises được sử dụng làm cơ sở đánh giá độ bền. Các khái niệm then chốt được phát triển xuyên suốt bao gồm: đồ thị công chỉ thị P-V, đồ thị gia tốc dịch chuyển Briks, hệ số kết cấu động cơ λ bằng 0,25, điều kiện biên liên kết tiếp xúc (Contact Constraint) và ngàm giữ cố định (Clamp Constraint) trong không gian 3D.

Phương pháp nghiên cứu

Nghiên cứu sử dụng nguồn dữ liệu kỹ thuật từ 1 mẫu động cơ Mazda B6 tiêu chuẩn nguyên bản với 4 xilanh thẳng hàng. Cỡ mẫu nghiên cứu bao gồm bộ thông số hình học và cơ tính hoàn chỉnh của 1 trục khuỷu nguyên khối gồm 5 cổ trục chính (đường kính 50 mm) và 4 chốt khuỷu (đường kính 45 mm), kết hợp cùng 4 bộ piston đường kính 78 mm và 4 thanh truyền tiêu chuẩn có chiều dài tương ứng. Phương pháp chọn mẫu là chọn mẫu có chủ đích điển hình nhằm đại diện chính xác cho kết cấu động cơ xăng I4 phổ thông.

Lý do lựa chọn phương pháp phân tích kết hợp giữa giải tích toán học truyền thống và mô phỏng số 3D CAE là nhằm khắc phục nhược điểm của các công thức giải tích phẳng vốn bỏ qua sự tập trung ứng suất tại các góc lượn chuyển tiếp phức tạp. Quy trình phân tích được thực hiện tuần tự qua các module Part Design, Assembly Design và Generative Structural Analysis của phần mềm CATIA V5R21, chia lưới phần tử tứ diện bậc cao với kích thước mắt lưới từ 1,5 mm đến 3,0 mm. Toàn bộ quá trình thu thập thông số, thiết lập mô hình toán và giải bài toán mô phỏng được hoàn thành trong khung thời gian 12 tháng của năm học 2022 - 2023.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình tính toán nhiệt động và mô phỏng số đã mang lại 4 kết quả định lượng quan trọng:

  • Xác lập chính xác các thông số nhiệt động của động cơ B6: Thể tích buồng cháy Vc đạt 0,048 lít, thể tích công tác một xilanh Vh đạt 0,399 lít và thể tích toàn phần Va đạt 0,447 lít, tương ứng với tỷ số nén danh định 9,3:1.
  • Xác định phổ tải trọng cực hạn trong chu trình công tác: Lực khí thể cực đại đạt 5,3 MN/m² tại thời điểm sau điểm chết trên, tạo ra lực pháp tuyến Zmax và lực tiếp tuyến Tmax tác động trực tiếp lên hệ thống chốt khuỷu và má khuỷu.
  • Nhận diện vùng tập trung ứng suất nguy hiểm: Kết quả mô phỏng 3D chỉ ra ứng suất Von-Mises cao nhất xuất hiện tại vị trí góc lượn chuyển tiếp giữa chốt khuỷu và má khuỷu, với biên độ ứng suất tăng hơn 35% so với giá trị trung bình trên thân trục.
  • Đánh giá chuyển vị biến dạng kết cấu: Độ biến dạng đàn hồi tổng thể của trục khuỷu dưới tác dụng của lực tổng hợp đạt giá trị cực đại trong giới hạn cho phép (dưới 0,08 mm), đảm bảo màng dầu bôi trơn duy trì ổn định không bị phá vỡ.

Thảo luận kết quả

Nguyên nhân dẫn đến sự tập trung ứng suất cục bộ tại các góc lượn là do sự thay đổi đột ngột về hình học và diện tích tiết diện giữa chốt khuỷu đường kính 45 mm và má khuỷu phẳng. Khi so sánh với các nghiên cứu động cơ đốt trong cùng phân khúc 1,6 lít, kết cấu trục khuỷu nguyên 5 cổ trục của động cơ B6 thể hiện độ cứng vững xoắn vượt trội hơn khoảng 18% so với kết cấu trục khuỷu thiếu cổ. Dữ liệu động học và phân bố lực có thể được trình bày trực quan qua đồ thị công P-V, đồ thị mô-men xoắn tổng hợp theo góc quay trục khuỷu và bảng tổng hợp ứng suất đa trục tại từng chế độ tải.

Mô hình hóa màu sắc 3D trong CATIA cho phép quan sát trực tiếp các dải đẳng trị ứng suất, giúp đối chiếu sai số giữa phương pháp giải tích lý thuyết và mô phỏng phần tử hữu hạn ở mức dưới 7,5%. Điều này chứng minh tính chuẩn xác và độ tin cậy cao của quy trình thiết kế mô phỏng số.

Đề xuất và khuyến nghị

Dựa trên kết quả phân tích động lực học và biến dạng kết cấu, 4 nhóm giải pháp kỹ thuật cụ thể được đề xuất nhằm nâng cao độ tin cậy vận hành:

  • Tối ưu hóa bán kính góc lượn chuyển tiếp giữa cổ trục, chốt khuỷu và má khuỷu: Tăng bán kính góc lượn từ 2,5 mm lên 3,2 - 3,5 mm nhằm giảm thiểu từ 12% đến 15% hệ số tập trung ứng suất; giải pháp do đội ngũ kỹ sư thiết kế R&D thực hiện trong vòng 6 tháng thử nghiệm mẫu.
  • Nâng cao chất lượng cơ tính bề mặt làm việc: Ứng dụng công nghệ thấm nitơ hoặc tôi cao tần cho các bề mặt làm việc của cổ trục và chốt khuỷu để đạt độ cứng từ 55 đến 60 HRC; mục tiêu giảm 25% tốc độ mài mòn cơ học, triển khai tại các xưởng gia công chế tạo trong lộ trình 12 tháng.
  • Chuẩn hóa quy trình mô phỏng CAE trong thiết kế động cơ: Tích hợp phần mềm CATIA và các công cụ mô phỏng phần tử hữu hạn chuyên sâu vào quy trình kiểm nghiệm thiết kế chi tiết máy tại các cơ sở đào tạo và viện nghiên cứu trong giai đoạn 2024 - 2025; giúp cắt giảm 35% chi phí chế tạo mẫu thử vật lý.
  • Kiểm soát và nâng cấp hệ thống bôi trơn cưỡng bức: Duy trì áp suất dầu bôi trơn ổn định trong dải 0,3 - 0,5 MPa tại các bạc lót đầu to thanh truyền và cổ trục chính; bộ phận kỹ thuật bảo dưỡng ô tô cần định kỳ kiểm tra khe hở dầu nhằm tăng tuổi thọ làm việc của cụm chi tiết thêm ít nhất 20%.

Đối tượng nên tham khảo luận văn

Nội dung và phương pháp luận của luận văn mang lại giá trị thực tiễn cho 4 nhóm đối tượng:

  • Kỹ sư thiết kế động cơ và cơ khí ô tô: Tham khảo quy trình tham số hóa 3D và các thông số tải trọng động lực học của động cơ 1,6L để ứng dụng trực tiếp vào quá trình thiết kế, cải tiến các cụm chi tiết máy tương đương.
  • Giảng viên, học viên cao học và sinh viên kỹ thuật động lực: Sử dụng làm tài liệu học thuật mẫu về phương pháp tính toán động lực học động cơ đốt trong kết hợp mô phỏng phần tử hữu hạn CAE, giúp chuẩn hóa nội dung giảng dạy và nghiên cứu chuyên sâu.
  • Kỹ thuật viên chẩn đoán, đại tu và bảo dưỡng ô tô: Nắm vững các vùng xung yếu dễ nứt vỡ và biến dạng của trục khuỷu, thanh truyền để tối ưu hóa quy trình kiểm tra không phá hủy (NDT) và phục hồi chi tiết chính xác.
  • Chuyên gia phân tích mô phỏng kết cấu (CAE): Tiếp cận phương pháp thiết lập điều kiện biên, ràng buộc lắp ráp và kỹ thuật chia lưới phần tử hữu hạn trên phần mềm CATIA V5R21 ứng dụng cho các bài toán cơ học phi tuyến phức tạp.

Câu hỏi thường gặp

Cơ cấu trục khuỷu động cơ Mazda B6 có đặc điểm kết cấu nổi bật nào?
Trục khuỷu động cơ B6 thuộc loại trục khuỷu nguyên khối với 4 chốt khuỷu và 5 cổ trục chính đường kính 50 mm. Kết cấu đủ cổ này giúp phân bổ tải trọng uốn đồng đều, giảm độ võng của trục khi chịu áp suất cháy cực đại 5,3 MN/m² và tăng cường độ cứng vững xoắn trong suốt chu trình làm việc.

Tham số kết cấu λ = 0,25 có ảnh hưởng gì đến động học của thanh truyền?
Hệ số λ là tỷ số giữa bán kính quay của khuỷu trục (R = 41,8 mm) và chiều dài thanh truyền (L = 167,2 mm). Giá trị λ = 0,25 giúp cân bằng tối ưu giữa chiều cao tổng thể của động cơ và độ lớn của lực ngang tác dụng lên thành xilanh, giảm thiểu ma sát mài mòn piston khoảng 15%.

Quy trình mô phỏng biến dạng trên phần mềm CATIA V5R21 được tiến hành qua các bước nào?
Quy trình gồm 4 bước chính: Dựng hình 3D tham số trong module Part Design, lắp ráp cụm cơ cấu trong Assembly Design, chuyển sang Generative Structural Analysis để chia lưới phần tử tứ diện bậc cao, áp đặt điều kiện biên (ngàm, liên kết tiếp xúc) và giải bài toán ứng suất Von-Mises dưới tải trọng cực đại.

Vị trí nào trên cụm cơ cấu chịu ứng suất tập trung lớn nhất?
Kết quả phân tích phần tử hữu hạn chỉ ra vị trí góc lượn chuyển tiếp giữa chốt khuỷu và má khuỷu là nơi chịu ứng suất Von-Mises cao nhất, vượt hơn 35% so với mức ứng suất trung bình thân trục. Đây là khu vực có nguy cơ mỏi và nứt gãy cao nhất cần được gia công góc lượn chính xác.

Độ tin cậy của phương pháp mô phỏng số FEM so với tính toán giải tích được đánh giá ra sao?
Kết quả đối chiếu giữa phương pháp giải tích cổ điển và mô phỏng số 3D FEM cho thấy mức độ sai lệch chỉ dưới 7,5%. Mô phỏng FEM phản ánh chân thực hơn sự phân bố ứng suất cục bộ đa chiều và biến dạng đàn hồi tại các vị trí hình học không gian phức tạp mà công thức giải tích phẳng chưa bao quát hết.

Kết luận

  • Hoàn thành trọn vẹn bộ tính toán nhiệt động và động lực học cho động cơ xăng B6 1,6L trên xe Mazda 323 với áp suất cực đại 5,3 MN/m² và tỷ số nén 9,3:1.
  • Thiết lập thành công mô hình 3D chính xác của cụm cơ cấu piston, thanh truyền và trục khuỷu 5 cổ chính trên nền tảng phần mềm CATIA V5R21.
  • Ứng dụng hiệu quả phương pháp phần tử hữu hạn (FEM) để xác định chính xác các điểm tập trung ứng suất Von-Mises và dải biến dạng chuyển vị của trục khuỷu.
  • Chứng minh tính ưu việt của giải pháp tích hợp CAD/CAE trong cơ khí ô tô, giúp kiểm soát sai số thiết kế dưới 7,5% và tối ưu hóa kết cấu bền vững.
  • Đề xuất các giải pháp kỹ thuật khả thi về bán kính góc lượn và tôi cứng bề mặt đạt 55 - 60 HRC nhằm kéo dài tuổi thọ động cơ trong giai đoạn vận hành 2024 - 2025.

Đóng góp chính của luận văn là xây dựng khung phương pháp luận kết hợp chuẩn xác giữa tính toán lý thuyết và mô phỏng số hóa 3D cho động cơ đốt trong. Để tiếp tục phát triển đề tài, các viện nghiên cứu và doanh nghiệp cơ khí nên mở rộng mô phỏng phân tích dao động xoắn và ứng suất nhiệt đa trường vật lý. Hãy liên hệ và khai thác toàn văn tài liệu kỹ thuật để chuẩn hóa quy trình thiết kế động cơ ngay hôm nay.