Tổng quan nghiên cứu

Trong ngành cơ khí chế tạo máy hiện đại, bộ truyền bánh răng chiếm hơn 85% tổng số cơ cấu truyền động cơ học nhờ hiệu suất truyền lực cao và tỷ số truyền ổn định. Tuy nhiên, khoảng 70% đến 80% chi phí và thời gian chu kỳ sản xuất sản phẩm mới phụ thuộc vào giai đoạn thiết kế và thử nghiệm kỹ thuật. Các phương pháp mô hình hóa bánh răng truyền thống chủ yếu dựa trên công thức giải tích 2D hoặc các phần mềm vẽ tĩnh, thường bỏ qua các hiện tượng động học phức tạp trong quá trình gia công thực tế như cắt lẹm chân răng, nhọn đỉnh răng và sai số bao hình. Tại Việt Nam, việc ứng dụng công nghệ mô phỏng số trực tiếp (Direct Digital Simulation - DDS) nhằm số hóa quá trình tạo hình bề mặt xoắn ốc thân khai 3D vẫn còn nhiều hạn chế, dẫn đến sự phụ thuộc vào các phần mềm thương mại đắt tiền nhưng thiếu khả năng tùy biến sâu.

Luận văn thạc sĩ kỹ thuật với đề tài "Lập trình mô phỏng động học việc tạo hình bánh răng trụ răng nghiêng dùng thanh răng sinh" do học viên Trần Ngọc Nhẫn thực hiện dưới sự hướng dẫn khoa học của Phó Giáo sư, Tiến sĩ Phạm Huy Hoàng tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh (thời gian thực hiện từ ngày 04/01/2012 đến ngày 02/12/2012) được xây dựng nhằm giải quyết triệt để bài toán trên. Mục tiêu trọng tâm của nghiên cứu là phát triển một chương trình mô phỏng tự động hóa trên nền tảng Autodesk Inventor VBA, tái hiện chính xác 100% quá trình bao hình không gian 3 chiều của bánh răng trụ răng nghiêng với góc nghiêng bất kỳ bằng dao thanh răng sinh. Nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm ước tính 40% thời gian thiết kế thử nghiệm, kiểm soát sai số động học biên dạng dưới mức 0,05 mm, đồng thời xuất ra các mô hình 3D kỹ thuật số tham số hóa hoàn chỉnh phục vụ công tác giảng dạy, nghiên cứu chuyên sâu và gia công CNC đa trục.

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 được xây dựng trên nền tảng tích hợp của 3 trụ cột lý thuyết cơ học và hình học vi phân chính xác:

Thứ nhất, lý thuyết ăn khớp bánh răng thân khai không gian. Biên dạng răng được tạo thành từ đường thân khai của vòng tròn cơ sở có đường kính $d_b = d \cdot \cos(\alpha_t)$, trong đó quỹ đạo điểm tiếp xúc được mô tả qua phương trình hàm thân khai $inv(\alpha) = \tan(\alpha) - \alpha$. Đối với bánh răng trụ răng nghiêng, các thông số được phân định rõ ràng trên mặt cắt pháp tuyến và mặt cắt tiếp tuyến, bao gồm mô đun pháp $m_n$, mô đun tiếp tuyến $m_t = m_n / \cos(\beta)$, góc áp lực tiêu chuẩn $\alpha = 20^\circ$, và góc nghiêng răng $\beta$ trong phạm vi từ $0^\circ$ đến $45^\circ$. Mức độ êm dịu khi truyền động được định lượng qua hệ số trùng khớp tổng cộng gồm hệ số trùng khớp ngang $\epsilon_\alpha$ và hệ số trùng khớp dọc $\epsilon_\beta = b \cdot \sin(\beta) / (\pi \cdot m_n)$.

Thứ hai, lý thuyết bao hình bề mặt bằng thanh răng sinh tương đương. Thanh răng sinh được xem là trường hợp đặc biệt của bánh răng có bán kính vòng chia tiến đến vô cùng. Các thông số thanh răng chuẩn hóa theo tiêu chuẩn công nghiệp gồm hệ số chiều cao đỉnh răng $h_a^* = 1,0$, hệ số khe hở chân răng $c^* = 0,25$, bán kính góc lượn mũi dao $r_c = 0,3 \cdot m_n$, chiều cao đỉnh dao $h_a = 1,25 \cdot m_n$, và chiều cao chân dao $h_f = 1,25 \cdot m_n$. Mối quan hệ động học tạo hình được xác lập bởi phương trình liên hệ giữa vận tốc tịnh tiến của dao $v$ và vận tốc góc của phôi $\omega$: $v = \omega \cdot d / 2 = \omega \cdot m_n \cdot Z / (2 \cdot \cos\beta)$.

Thứ ba, lý thuyết dịch chỉnh biên dạng và điều kiện biên tránh cắt chân răng. Hiện tượng cắt chân răng xảy ra khi đỉnh của dao thanh răng cắt lẹm vào phần gốc răng dưới vòng tròn cơ sở. Số răng tối thiểu không bị cắt chân răng khi gia công không dịch chỉnh được tính theo công thức $Z_{min} = 2 / \sin^2(\alpha_t)$. Đối với bánh răng thẳng tiêu chuẩn ($\alpha = 20^\circ, \beta = 0^\circ$), $Z_{min} = 18$ răng; khi $\alpha = 15^\circ$, $Z_{min} = 30$ răng; và khi tăng góc nghiêng $\beta = 45^\circ$, $Z_{min}$ giảm xuống còn 13 răng. Luận văn áp dụng hệ số dịch chỉnh $X$ với điều kiện giới hạn $X_{min} = 1 - (Z \cdot \sin^2\alpha_t) / 2$ nhằm triệt tiêu hiện tượng cắt chân răng khi số răng $Z < 17$, đồng thời tính toán để ngăn ngừa hiện tượng nhọn đỉnh răng khi dịch chỉnh dương quá mức.

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

Nghiên cứu thu thập và thiết lập bộ dữ liệu thực nghiệm gồm 24 mô hình bánh răng với các dải thông số biến thiên đại diện: mô đun $m_n$ từ 1,0 mm đến 6,0 mm, số răng $Z$ từ 10 đến 50 răng, góc nghiêng $\beta$ từ $0^\circ$ đến $45^\circ$, và hệ số dịch chỉnh $X$ từ -0,5 đến +0,8. Phương pháp chọn mẫu là chọn mẫu phân tầng có chủ đích theo tiêu chuẩn chế tạo máy JIS B 1701-1973 của Nhật Bản, nhằm bao quát đầy đủ 3 trạng thái hình học: bánh răng bị cắt chân răng nghiêm trọng ($Z = 12, X = 0$), bánh răng tiêu chuẩn chuẩn tắc ($Z = 20, X = 0$), và bánh răng dịch chỉnh khắc phục cắt lẹm ($Z = 12, X = +0,4$).

Phương pháp phân tích cốt lõi là mô phỏng số trực tiếp (Direct Digital Simulation - DDS) kết hợp lập trình giao diện ứng dụng (API) thông qua công nghệ ActiveX Automation trong môi trường Autodesk Inventor VBA. Lý do lựa chọn giải pháp này thay vì các phần mềm như AutoCAD VBA, SolidWorks hay OpenGL C++ là vì Autodesk Inventor sở hữu hệ thống đối tượng 3D chuyên sâu (PartDocument, TransientBRep, SolidBody, Sketch3D), cho phép thực hiện các phép toán Boolean cắt gọt thể tích rắn theo thời gian thực với độ chính xác hình học cao. Toàn bộ tiến trình mô phỏng được điều khiển thông qua bộ điều khiển thời gian IeTimer.ocx, chia hành trình cắt thành chuỗi $n$ bước rời rạc để tính toán chính xác sai số động học. Timeline nghiên cứu được thực hiện xuyên suốt trong 11 tháng liên tục từ khâu xây dựng thuật toán toán học đến hoàn thiện bộ cài đặt và kiểm chứng sai số.

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

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

Quá trình mô phỏng và khảo sát thực nghiệm trên mô hình số đã mang lại 4 phát hiện kỹ thuật mang tính định lượng rõ rệt:

Thứ nhất, xây dựng thành công thuật toán đồng bộ động học 3 chiều, chuyển đổi quá trình cắt gọt liên tục thành chuỗi $n$ bước vi phân thời gian điều khiển bằng IeTimer. Với dải bước chia từ 60 đến 360 bước trên một chu kỳ ăn khớp, chương trình đã tái hiện chính xác quỹ đạo chuyển động tịnh tiến của thanh răng sinh kết hợp chuyển động quay của phôi bánh răng nghiêng theo đúng góc xoắn $\beta$.

Thứ hai, định lượng chính xác hiện tượng cắt chân răng và hiệu quả của dịch chỉnh biên dạng. Khi khảo sát bánh răng có $Z = 12, \alpha = 20^\circ, \beta = 0^\circ$ không dịch chỉnh ($X = 0$), thể tích vật liệu chân răng bị cắt lẹm làm giảm chiều dày chân răng $S_f$ khoảng 18,5% so với biên dạng lý thuyết, gây suy giảm nghiêm trọng độ bền uốn. Khi áp dụng hệ số dịch chỉnh dương $X = +0,38$, hiện tượng cắt lẹm bị triệt tiêu hoàn toàn, chiều dày chân răng tăng 22,4%, khôi phục toàn bộ đoạn thân khai hữu ích gần vòng cơ sở.

Thứ ba, làm rõ mối quan hệ giữa góc nghiêng $\beta$ và việc giảm thiểu kích thước cơ cấu. Khi tăng góc nghiêng $\beta$ từ $0^\circ$ lên lần lượt $15^\circ$, $30^\circ$, và $45^\circ$ với góc áp lực $\alpha = 20^\circ$, số răng tối thiểu $Z_{min}$ không bị cắt chân răng giảm từ 18 răng xuống tương ứng là 16, 14, và 13 răng (giảm tương đương 27,8%). Phát hiện này chứng minh rằng việc tăng góc nghiêng cho phép thiết kế các bánh răng có số răng rất nhỏ mà không cần dịch chỉnh quá mức, giúp giảm kích thước hướng tâm của hộp giảm tốc khoảng 22%.

Thứ tư, xác định quy luật biến thiên của sai số động học biên dạng răng ($H$). Kết quả đo đạc trên mô hình số cho thấy khi tăng mật độ bước chia mô phỏng $n$ từ 60 bước lên 240 bước, sai số biên dạng cực đại $H$ giảm từ 0,082 mm xuống còn 0,022 mm (giảm 73,2%), đạt độ chính xác bề mặt tương đương cấp chính xác 7 theo tiêu chuẩn gia công cơ khí quốc tế.

Thảo luận kết quả

Nguyên nhân cốt lõi tạo nên sai số động học $H$ trong phương pháp DDS là do việc xấp xỉ bề mặt cong thân khai liên tục bằng một tập hợp hữu hạn các mặt phẳng cắt từ các vết gọt của thanh răng sinh. Số bước cắt $n$ càng lớn thì các vết gọt càng mịn, làm cho sai số $H$ tiệm cận về 0 nhưng sẽ làm tăng thời gian tính toán của CPU từ vài giây lên hàng chục phút. So với các nghiên cứu trước đây như nghiên cứu mô phỏng trên nền Flash MX của Đào Ngọc Biên (vốn chỉ mang tính trình diễn 2D) hay mô phỏng 2D trên AutoCAD VBA của Diệp Bảo Duy (chưa xử lý được mặt xoắn ốc 3D), giải pháp trên Autodesk Inventor VBA của luận văn này đã vượt trội hoàn toàn về khả năng tạo hình khối rắn 3D tham số hóa.

Dữ liệu thực nghiệm về mối tương quan đa biến giữa sai số động học $H$, hệ số dịch chỉnh $X$, góc nghiêng $\beta$ và số bước chia $n$ có thể được trực quan hóa một cách tối ưu thông qua biểu đồ bề mặt 3D (3D Surface Plot) hoặc bảng ma trận so sánh 5 cột bao gồm: Mô đun ($m_n$), Góc xoắn ($\beta$), Hệ số dịch chỉnh ($X$), Số bước cắt ($n$), và Sai số đo lường ($H$). Việc xuất thành công các mô hình 3D kỹ thuật số sang các định dạng chuẩn công nghiệp như STEP (.stp), ACIS (.sat), DWG và CATPart mang ý nghĩa đặc biệt quan trọng, cho phép liên kết trực tiếp với các phần mềm CAE (Ansys, Nastran) để phân tích độ bền ứng suất tiếp xúc và phần mềm CAM (Mastercam, PowerMill) để lập trình gia công CNC 5 trục mà không cần dựng lại hình học từ đầu.

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

Dựa trên kết quả đạt được, luận văn đưa ra 4 khuyến nghị và giải pháp hành động cụ thể nhằm nâng cao hiệu quả ứng dụng thực tiễn:

Thứ nhất, nâng cấp thuật toán điều khiển thời gian sang cơ chế bước chia thích ứng (Adaptive Step Control). Nhóm nghiên cứu cần tối ưu hóa mã lệnh để chương trình tự động tăng mật độ vết cắt tại vùng chân răng và giảm mật độ tại vùng ít cong, giúp rút ngắn 35% thời gian tính toán mô phỏng trong vòng 6 tháng tới. Chủ thể thực hiện: Nhóm phát triển thuật toán và lập trình CAD/CAM.

Thứ hai, mở rộng phạm vi mô phỏng sang các công cụ cắt bao hình phức tạp khác như dao phay lăn trục vít (Hobbing) và dao xọc răng dạng bánh răng (Fellows). Mục tiêu hướng đến việc kiểm soát sai số động học dưới ngưỡng 0,01 mm cho cả bánh răng côn xoắn và bánh răng chữ V trong lộ trình 12 tháng. Chủ thể thực hiện: Phòng thí nghiệm Cơ khí chính xác và CAD/CAM tại các trường đại học kỹ thuật.

Thứ ba, tích hợp toàn diện bộ phần mềm mô phỏng DDS vào chương trình giảng dạy thực nghiệm số cho 100% sinh viên chuyên ngành Cơ khí chế tạo máy và Kỹ thuật công nghiệp trong thời gian 18 tháng. Giải pháp này giúp sinh viên quan sát trực quan hiện tượng dịch chỉnh và cắt chân răng, nâng cao tỷ lệ tiếp thu kiến thức thiết kế máy trên 90%. Chủ thể thực hiện: Bộ môn Thiết kế máy và Khoa Cơ khí.

Thứ tư, tái cấu trúc mã nguồn từ VBA sang dạng Add-in chuyên nghiệp sử dụng ngôn ngữ C# và nền tảng .NET Framework. Việc này sẽ tận dụng khả năng xử lý đa luồng (Multi-threading) của CPU hiện đại, tăng tốc độ xử lý đồ họa lên 50% và hỗ trợ xuất dữ liệu sang chuẩn STEP AP242 trong vòng 9 tháng. Chủ thể thực hiện: Các kỹ sư phần mềm công nghiệp và các viện nghiên cứu thiết kế máy.

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

Nội dung luận văn mang lại giá trị thực tiễn và học thuật sâu sắc cho 4 nhóm đối tượng chính:

  1. Giảng viên và nghiên cứu viên ngành Kỹ thuật Cơ khí: Sử dụng tài liệu và phần mềm như một giáo cụ trực quan sinh động để giảng dạy các môn học Nguyên lý chi tiết máy, Thiết kế CAD/CAM/CAE và Công nghệ chế tạo phụ tùng máy.

  2. Học viên cao học và sinh viên đại học chuyên ngành Chế tạo máy, Cơ điện tử: Tham khảo phương pháp xây dựng mô hình toán học đường thân khai 3D, kỹ thuật lập trình tương tác ActiveX API trên phần mềm thiết kế chuyên dụng và phương pháp phân tích sai số hình học.

  3. Kỹ sư thiết kế R&D tại các nhà máy sản xuất cơ khí và hộp giảm tốc: Ứng dụng công cụ để thiết kế tối ưu các bộ truyền bánh răng phi tiêu chuẩn, tính toán chính xác hệ số dịch chỉnh nhằm tránh hiện tượng cắt lẹm hoặc nhọn đỉnh răng trước khi đưa vào sản xuất hàng loạt.

  4. Kỹ sư lập trình CAM và vận hành máy CNC đa trục: Sử dụng các tệp tin 3D kỹ thuật số (STEP, SAT) được xuất trực tiếp từ mô phỏng để phân tích đường chạy dao, kiểm tra va chạm động học và lập trình gia công bánh răng chính xác cao.

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

Phương pháp Direct Digital Simulation (DDS) khác gì so với phương pháp dựng hình CAD truyền thống?

Phương pháp CAD truyền thống tạo biên dạng răng bằng cách nội suy toán học các đường cong lý thuyết tĩnh. Ngược lại, phương pháp DDS tái hiện chính xác các chuyển động cắt gọt động học tương đối giữa dao thanh răng và phôi, tạo ra biên dạng răng thực tế bao gồm cả các sai số hình học và hiện tượng cắt chân răng.

Tại sao tăng góc nghiêng $\beta$ lại giúp giảm nguy cơ cắt chân răng ở bánh răng có ít răng?

Khi tăng góc nghiêng $\beta$, số răng tương đương của bánh răng tăng lên theo công thức $Z_v = Z / \cos^3(\beta)$. Khi góc áp lực pháp tuyến $\alpha = 20^\circ$ và $\beta = 45^\circ$, số răng tối thiểu $Z_{min}$ để không xảy ra cắt chân răng giảm từ 18 răng xuống chỉ còn 13 răng, giúp gia công bánh răng nhỏ an toàn.

Việc lựa chọn Autodesk Inventor VBA mang lại ưu thế gì so với ngôn ngữ lập trình độc lập như C++ OpenGL?

Autodesk Inventor VBA cung cấp sẵn hệ thống thư viện đối tượng 3D API hoàn chỉnh và công cụ tính toán hình học khối rắn mạnh mẽ. Người nghiên cứu không cần phải tự viết lại các thuật toán đồ họa phức tạp từ đầu như trong OpenGL, đồng thời có thể trích xuất trực tiếp các đặc tính vật lý và tệp CAD chuẩn.

Hệ số dịch chỉnh $X$ có giá trị bao nhiêu là tối ưu để loại bỏ hiện tượng cắt chân răng khi $Z = 12$?

Theo tính toán lý thuyết và kiểm chứng mô phỏng với $\alpha = 20^\circ$ và $\beta = 0^\circ$, hệ số dịch chỉnh tối thiểu là $X_{min} = 1 - (12 \cdot \sin^2 20^\circ) / 2 \approx +0,30$. Giá trị tối ưu được khuyến nghị là $X = +0,35$ đến $+0,40$, vừa triệt tiêu hoàn toàn cắt chân răng vừa tránh làm nhọn đỉnh răng.

Sai số động học $H$ trong mô phỏng bị chi phối nhiều nhất bởi yếu tố nào?

Sai số động học $H$ phụ thuộc chủ yếu vào số lượng bước chia thời gian $n$ của chu trình gia công và bước dịch chuyển tịnh tiến của thanh răng sinh. Khi tăng số bước chia $n$ từ 60 lên 240 bước, sai số $H$ giảm tới 73,2%, giúp bề mặt thân khai đạt độ mịn và độ chính xác hình học cao.

Kết luận

  • Xây dựng thành công bộ chương trình mô phỏng động học 3D quá trình tạo hình bánh răng trụ răng nghiêng dùng thanh răng sinh trên nền tảng Autodesk Inventor VBA.
  • Giải quyết triệt để và trực quan hóa các bài toán phi tuyến hình học phức tạp bao gồm hiện tượng cắt chân răng, nhọn đỉnh răng và dịch chỉnh góc ăn khớp.
  • Thiết lập mô hình toán học và phương pháp đo kiểm sai số động học biên dạng răng ($H$), chứng minh khả năng giảm sai số trên 70% khi tối ưu hóa bước chia thời gian.
  • Tạo lập hoàn chỉnh nguồn mô hình 3D kỹ thuật số tham số hóa đa định dạng (STEP, SAT, DWG), phục vụ liên thông trực tiếp với hệ sinh thái phần mềm CAE/CAM.
  • Khẳng định tính khả thi của việc ứng dụng công nghệ mô phỏng số trực tiếp (DDS) nhằm nâng cao chất lượng đào tạo đại học và tối ưu hóa quy trình thiết kế cơ khí tại Việt Nam.

Nghiên cứu đã đóng góp một công cụ khoa học có giá trị cao, kết nối chặt chẽ giữa lý thuyết giải tích bánh răng và thực tiễn sản xuất công nghiệp. Trong lộ trình 6 đến 12 tháng tiếp theo, nhóm nghiên cứu sẽ tập trung nâng cấp mã nguồn sang chuẩn C# Add-in đa luồng và mở rộng thuật toán cho các loại dao cắt định hình phức tạp. Các tổ chức giáo dục, viện nghiên cứu và doanh nghiệp cơ khí chế tạo hãy áp dụng ngay giải pháp mô phỏng số này vào quy trình R&D để rút ngắn thời gian phát triển sản phẩm và tối ưu hóa chi phí sản xuất.