Tổng quan nghiên cứu

Trong xu hướng phát triển mạnh mẽ của cuộc Cách mạng Công nghiệp lần thứ tư, công nghệ tạo mẫu nhanh với hơn 30 phương pháp thương mại hóa đã mở ra bước đột phá toàn diện cho ngành công nghệ chế tạo máy. Trong các hệ thống sản xuất bồi đắp hiện hành, công nghệ mô hình hóa lắng đọng nóng chảy FDM (Fused Deposition Modelling) giữ vị thế chiến lược nhờ ưu điểm vận hành linh hoạt, không sử dụng nguồn năng lượng laser đắt đỏ và có khả năng gia công trực tiếp từ mô hình CAD 3D sang sản phẩm vật lý. Tuy nhiên, rào cản lớn nhất ngăn FDM chuyển đổi từ giai đoạn tạo mẫu thẩm mỹ sang chế tạo trực tiếp các chi tiết kỹ thuật chịu tải là chất lượng đầu ra chưa ổn định, trong đó sai số kích thước hình học có thể vượt quá 3.40% và độ nhám bề mặt Ra thường dao động lớn từ 1.98 µm đến 38.60 µm tùy thuộc vào chế độ cài đặt.

Mục tiêu trọng tâm của luận văn thạc sĩ chuyên ngành Công nghệ Chế tạo Máy là nghiên cứu định lượng mức độ ảnh hưởng của ba thông số công nghệ cốt lõi gồm độ dày lớp in từ 0.116 mm đến 0.284 mm, nhiệt độ đầu đùn trong dải 177°C đến 193°C và tốc độ điền đầy từ 22 mm/s đến 38 mm/s. Từ đó, tác giả xây dựng hàm mục tiêu toán học nhằm tối ưu hóa đồng thời độ chính xác kích thước không gian ba chiều và giảm thiểu độ nhám bề mặt cho sản phẩm chế tạo từ vật liệu nhựa nhiệt dẻo PLA.

Được thực hiện tại Phòng Cơ học máy thuộc Viện Cơ học và Tin học ứng dụng trong khoảng thời gian từ ngày 21/01/2013 đến ngày 21/06/2013, công trình đã xác lập bộ thông số công nghệ chuẩn xác, giúp giảm thiểu độ nhám bề mặt xuống dưới 2.50 µm và duy trì sai số kích thước trong ngưỡng 1.20%. Kết quả này không chỉ khẳng định tính cấp thiết về mặt học thuật mà còn mang lại giá trị kinh tế trực tiếp, giúp rút ngắn hơn 50% thời gian thiết kế thử nghiệm và giảm thiểu tỷ lệ phế phẩm trong sản xuất khuôn mẫu và chế tạo chi tiết cơ khí chính xá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 dựa trên sự tích hợp của các hệ thống lý thuyết nền tảng trong cơ học vật liệu polyme và toán học ứng dụng:

  • Lý thuyết nhiệt động học dòng chảy polyme phi Newton: Giải thích hành vi hóa dẻo và biến dạng chảy của sợi nhựa nhiệt dẻo PLA có công thức phân tử (C3H4O2)n khi đi qua buồng gia nhiệt và miệng vòi đùn đường kính 0.375 mm. Dải nhiệt độ nóng chảy tiêu chuẩn của vật liệu từ 173°C đến 178°C quyết định trực tiếp đến độ nhớt động học, áp lực đẩy của hệ thống con lăn và khả năng duy trì hình dạng sợi nhựa sau khi rời vòi đùn.
  • Cơ chế bồi đắp và liên kết liên lớp (Layer-by-Layer Deposition Mechanics): Quá trình hình thành liên kết phân tử và độ bền bám dính giữa các lớp in liền kề phụ thuộc chặt chẽ vào hiện tượng khuếch tán nhiệt và ứng suất dư co ngót. Nếu ứng suất nhiệt cục bộ vượt quá giới hạn chịu đựng, chi tiết sẽ xuất hiện hiện tượng tách lớp, cong vênh góc hoặc tạo rãnh phân cách làm suy giảm độ bền kéo đạt chuẩn 69.90 MPa của vật liệu.
  • Phương pháp bề mặt đáp ứng (Response Surface Methodology - RSM): Cung cấp cơ sở toán học để thiết lập phương trình hồi quy đa thức bậc hai, mô tả mối quan hệ phi tuyến phức tạp giữa các biến công nghệ đầu vào với các chỉ tiêu chất lượng đầu ra.
  • Các khái niệm kỹ thuật cốt lõi:
    • Độ dày mỗi lớp (Layer height): Khoảng cách dịch chuyển tịnh tiến của bàn máy theo trục Z sau mỗi chu trình quét biên dạng.
    • Tốc độ điền đầy (Speed infill): Vận tốc di chuyển tương đối của đầu đùn trong quá trình quét vật liệu bên trong đường bao với mật độ đặc 100%.
    • Độ nhám bề mặt (Surface roughness - Ra): Sai lệch số học trung bình của biên dạng nhám đo đạc tại mặt trên (RaT) và mặt bên (RaS) của mẫu.

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

Nguồn dữ liệu thực nghiệm được thu thập trực tiếp thông qua quy trình chế tạo mẫu thử nghiệm trên hệ thống máy in 3D ba chiều tại Phòng Cơ học máy. Mẫu thực nghiệm được thiết kế dạng khối hình học tiêu chuẩn trên phần mềm 3D SolidWorks, sau đó chuyển đổi định dạng dữ liệu sang tệp STL và xử lý mã lệnh G-code thông qua phần mềm Slic3r. Vật liệu sử dụng là sợi nhựa PLA đường kính 3.00 mm với mô đun đàn hồi 1968 MPa và độ bền uốn 94.00 MPa. Kích thước hình học ba chiều của mẫu được đo kiểm bằng thước cặp điện tử Mitutoyo có sai số chỉ 0.01 mm. Độ nhám bề mặt được xác định bằng máy đo độ nhám tiếp xúc Mitutoyo của Nhật Bản với độ phân giải 0.01 µm và độ phóng đại dọc lên đến 500,000X.

Về quy mô và phương pháp chọn mẫu, nghiên cứu áp dụng quy hoạch thực nghiệm bậc hai tâm xoay hỗn hợp quay đều Box-Hunter (Central Composite Design - CCD) với tổng cỡ mẫu gồm N = 20 thí nghiệm độc lập. Cấu trúc mẫu bao gồm 8 thí nghiệm tại các đỉnh của mô hình trực giao toàn phần 2^3, 6 thí nghiệm tại các điểm sao với khoảng cách cánh tay đòn sao alpha = 1.682, và 6 thí nghiệm lặp tại tâm quy hoạch (n0 = 6).

Phương pháp quy hoạch Box-Hunter được lựa chọn vì tính ưu việt vượt trội: cho phép giảm thiểu số lần thí nghiệm xuống mức tối ưu nhưng vẫn đảm bảo tính độc lập thống kê, xây dựng được bề mặt hồi quy liên tục và ước lượng chính xác sai số tái hiện thực nghiệm. Các hệ số hồi quy được giải bằng phương pháp bình phương bé nhất, sau đó kiểm tra mức ý nghĩa bằng tiêu chuẩn Student với tb = 2.57 tại mức ý nghĩa q = 0.05 và kiểm tra tính tương thích của mô hình thông qua tiêu chuẩn Fisher với Fb = 4.88. Toàn bộ thuật toán tối ưu hóa đa mục tiêu được lập trình và giải trên phần mềm Matlab trong suốt thời gian 5 tháng triển khai đề tài.

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

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

Quá trình phân tích số liệu từ 20 mẫu thực nghiệm đã xác định rõ mức độ chi phối của từng biến số công nghệ đến các chỉ tiêu chất lượng của mẫu in:

  • Sai số kích thước chiều dài (%L): Nhiệt độ đầu đùn là nhân tố ảnh hưởng mạnh mẽ nhất với hệ số tuyến tính dương b2 = 0.2126, tiếp theo là tác động tương hỗ giữa độ dày lớp và tốc độ điền đầy (b13 = 0.1750). Sai số tương đối chiều dài dao động từ mức thấp nhất 1.18% (kích thước thực 50.59 mm tại chế độ 0.15 mm, 180°C, 35 mm/s) lên đến mức cao nhất 2.98% (kích thước thực 51.49 mm tại chế độ 0.20 mm, 193°C, 30 mm/s), tương đương mức biến thiên sai số tăng thêm 152.54%.
  • Sai số kích thước chiều rộng (%W) và chiều cao (%T): Độ dày lớp và nhiệt độ đầu đùn tác động đáng kể đến sự giãn nở theo phương ngang. Sai số chiều rộng ghi nhận giá trị nhỏ nhất là 0.85% tại điểm sao độ dày lớp 0.116 mm và tăng lên cực đại 2.45% khi tăng đồng thời nhiệt độ lên 190°C và tốc độ lên 35 mm/s. Đối với chiều cao, sai số dao động trong dải từ 0.60% đến 3.40% (kích thước thực đạt từ 4.83 mm đến 4.98 mm so với thiết kế 5.00 mm).
  • Đặc tính độ nhám bề mặt:
    • Độ nhám mặt trên (RaT): Dao động trong khoảng rộng từ 19.65 µm đến 38.60 µm. Giá trị nhám tăng vọt thêm 96.44% khi độ dày lớp tăng từ 0.15 mm lên 0.25 mm dưới tốc độ điền đầy 35 mm/s.
    • Độ nhám mặt bên (RaS): Đạt chất lượng tối ưu ở mức 1.98 µm khi vận hành tại độ dày lớp cực tiểu 0.116 mm và tăng lên đến 4.80 µm khi sử dụng độ dày lớp 0.25 mm ở nhiệt độ 190°C (mức độ nhám tăng hơn 142.42%).

Thảo luận kết quả

Cơ chế vật lý giải thích hiện tượng tăng sai số kích thước khi nhiệt độ đầu đùn vượt ngưỡng 185°C bắt nguồn từ sự suy giảm độ nhớt động học của nhựa PLA. Dưới tác dụng của áp lực đùn liên tục, dòng chất lỏng polyme nóng chảy có xu hướng lan tỏa tự do sang hai bên trước khi kịp đông đặc hoàn toàn dưới nhiệt độ môi trường, dẫn đến sự mở rộng biên dạng thực tế so với bản vẽ CAD. Ngược lại, khi độ dày lớp in tăng lên, hiệu ứng bậc thang hình thành rõ nét hơn tại các cạnh tiếp giáp giữa các lớp, tạo nên các khe rãnh tế vi làm tăng mạnh chỉ số độ nhám bề mặt.

Toàn bộ dữ liệu thực nghiệm được trực quan hóa hiệu quả thông qua hệ thống biểu đồ bề mặt đáp ứng 3D (Response Surface Plot) và biểu đồ đẳng trị (Contour Plot) xuất từ phần mềm Matlab. Biểu đồ 3D thể hiện rõ vùng trũng tối ưu của sai số kích thước và độ nhám tập trung quanh điểm tâm quy hoạch với dải nhiệt độ 180°C - 185°C và tốc độ 25 mm/s - 30 mm/s. Đồng thời, bảng ma trận phân tích phương sai ANOVA chứng minh phương trình hồi quy có hệ số xác định cao, đạt độ tương thích chuẩn xác theo phân phối Fisher với Ftt < 4.88.

So sánh với các công bố quốc tế, kết quả nghiên cứu này tương đồng với nhận định của Rajan Bansal (2011) và Mohamad Alhubail (2012) về vai trò then chốt của độ dày lớp đối với độ nhám bề mặt. Đồng thời, luận văn đóng góp thêm phát hiện mới khi chứng minh rằng đối với dòng máy in 3D sử dụng nhựa nhiệt dẻo PLA, tương tác chéo giữa độ dày lớp và tốc độ điền đầy (X1X3) có mức độ ảnh hưởng đến sai số kích thước lớn hơn nhiều so với ảnh hưởng đơn lẻ của tốc độ đùn.

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

Nhằm nâng cao chất lượng và độ chính xác của sản phẩm tạo mẫu bằng công nghệ FDM, các giải pháp kỹ thuật cụ thể được đề xuất như sau:

  • Chuẩn hóa thông số công nghệ tối ưu cho quy trình gia công nhựa PLA: Thiết lập chế độ in tiêu chuẩn bao gồm độ dày lớp X1 = 0.15 mm, nhiệt độ đầu đùn X2 = 180°C và tốc độ điền đầy X3 = 25 mm/s. Chế độ này hướng tới mục tiêu duy trì sai số kích thước chiều dài dưới 1.20% và đưa độ nhám mặt bên RaS đạt ngưỡng dưới 2.50 µm. Kế hoạch áp dụng ngay trong quý đầu tiên do kỹ sư vận hành xưởng in 3D chịu trách nhiệm thực thi.
  • Ứng dụng thuật toán bù trừ co ngót nhiệt trên phần mềm CAM: Tích hợp hệ số bù trừ kích thước từ 1.50% đến 1.80% theo phương X-Y trực tiếp vào mô hình hình học trước khi chuyển đổi sang mã G-code trên phần mềm Slic3r. Giải pháp này giúp triệt tiêu sai lệch do giãn nở nhiệt, đưa dung sai chi tiết hoàn thiện về mức tiệm cận 0.05 mm. Thời gian triển khai trong vòng 3 tháng do nhóm kỹ sư thiết kế R&D phụ trách.
  • Nâng cấp hệ thống điều khiển nhiệt độ đầu phun bằng thuật toán PID chính xác cao: Thay thế bộ điều nhiệt thông thường bằng hệ thống điều khiển nhiệt độ vòng kín PID có khả năng kiểm soát độ trôi nhiệt trong phạm vi hẹp ±0.50°C xung quanh mức cài đặt 185°C. Giải pháp giúp loại bỏ hiện tượng đùn thừa hoặc tắc nghẽn vòi phun, hoàn thành lắp đặt sau 6 tháng dưới sự giám sát của bộ phận kỹ thuật bảo trì.
  • Mở rộng mô hình tối ưu hóa đa mục tiêu trên vật liệu kỹ thuật cao: Tiếp tục nghiên cứu ứng dụng quy hoạch thực nghiệm trên các họ vật liệu phức hợp như ABS M30i, PETG và nhựa gia cường sợi carbon. Mục tiêu đặt ra là tối ưu hóa đồng thời ba tiêu chí: độ chính xác kích thước dưới 0.80%, độ nhám Ra dưới 1.50 µm và gia tăng độ bền kéo vượt mức 40.00 MPa. Lộ trình thực hiện trong 12 tháng do nhóm nghiên cứu sau đại học tại các viện chuyên ngành chủ trì.

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

Công trình luận văn mang lại giá trị tham khảo khoa học và ứng dụng thực tiễn sâu rộng cho các nhóm đối tượng sau:

  • Học viên cao học và nghiên cứu sinh ngành Cơ khí chế tạo: Cung cấp tài liệu tham khảo chuẩn mực về phương pháp luận quy hoạch thực nghiệm tâm xoay Box-Hunter, quy trình xử lý số liệu bằng kiểm định thống kê Student - Fisher và kỹ thuật tối ưu hóa đa mục tiêu bằng Matlab.
  • Kỹ sư R&D và chuyên viên thiết kế sản phẩm: Tiếp cận hướng dẫn chi tiết về cách thức điều chỉnh thông số cắt lớp Slic3r nhằm chế tạo các mẫu thử nghiệm chức năng có độ chính xác cao, phục vụ hiệu quả cho công tác kiểm tra lắp ráp trong ngành công nghiệp ô tô và thiết bị y tế.
  • Nhà quản lý xưởng in 3D và các doanh nghiệp dịch vụ tạo mẫu nhanh: Ứng dụng trực tiếp bộ thông số công nghệ tối ưu để giảm thiểu tỷ lệ lỗi phế phẩm từ 3.40% xuống dưới 1.00%, tối ưu hóa chu kỳ sản xuất và tiết kiệm chi phí sợi nhựa PLA trong quá trình vận hành dịch vụ thương mại.
  • Giảng viên và cán bộ nghiên cứu tại các trường đại học kỹ thuật: Sử dụng làm giáo trình mẫu trong các học phần Công nghệ CAD/CAM/CNC, Kỹ thuật đo lường kiểm tra và Công nghệ vật liệu mới, giúp sinh viên nắm vững nguyên lý bồi đắp vật liệu hiện đại.

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

Vật liệu nhựa PLA có những đặc tính nổi bật nào để ứng dụng trong công nghệ FDM?
Nhựa sinh học PLA sở hữu độ bền kéo đứt đạt 69.90 MPa, mô đun đàn hồi 1968 MPa và nhiệt độ nóng chảy thấp trong dải 173°C - 178°C. Đặc tính co ngót nhiệt thấp giúp PLA ít bị biến dạng cong vênh hơn đáng kể so với nhựa ABS, rất phù hợp cho tạo mẫu chính xác.

Thông số công nghệ nào quyết định lớn nhất đến độ nhám bề mặt chi tiết in 3D?
Độ dày lớp in là nhân tố chi phối hàng đầu. Kết quả thực nghiệm chứng minh khi giảm độ dày lớp từ 0.25 mm xuống 0.116 mm, độ nhám mặt bên RaS giảm mạnh từ 4.80 µm xuống 1.98 µm (cải thiện hơn 58.75% độ mịn bề mặt) nhờ giảm thiểu hiệu ứng bậc thang.

Vì sao nhiệt độ đầu đùn quá cao lại làm gia tăng sai số kích thước sản phẩm?
Khi nhiệt độ vượt quá 190°C, độ nhớt của PLA giảm mạnh làm tăng tính lưu động của dòng chảy polyme. Dưới áp lực đùn, sợi nhựa bị tràn sang hai biên trước khi đông kết, làm sai số kích thước chiều dài %L tăng vọt từ 1.18% lên 2.98%.

Phương pháp quy hoạch thực nghiệm Box-Hunter mang lại lợi ích gì cho đề tài?
Phương pháp Box-Hunter cho phép khảo sát toàn diện 3 nhân tố với đầy đủ tương tác phi tuyến chỉ qua 20 lần thử nghiệm thay vì hàng trăm mẫu đơn biến. Quy hoạch này đảm bảo tính quay đều, giúp dự báo chính xác bề mặt đáp ứng với độ tin cậy thống kê 95%.

Kết quả nghiên cứu này có thể chuyển giao trực tiếp cho sản xuất công nghiệp không?
Hoàn toàn khả thi. Doanh nghiệp sử dụng máy in 3D FDM có đường kính vòi phun 0.375 mm có thể cài đặt ngay bộ thông số: độ dày lớp 0.15 mm, nhiệt độ 180°C và tốc độ 25 mm/s để chế tạo đồ gá, khuôn sáp đúc chính xác với sai số dưới 0.10 mm.

Kết luận

  • Xây dựng thành công phương trình hồi quy bậc hai mô tả chính xác tương tác phi tuyến giữa độ dày lớp, nhiệt độ đầu đùn và tốc độ điền đầy đến sai số kích thước và độ nhám bề mặt sản phẩm FDM.
  • Xác định bộ thông số công nghệ tối ưu: độ dày lớp 0.15 mm, nhiệt độ vòi phun 180°C và tốc độ điền đầy 25 mm/s, giúp kiểm soát sai số kích thước dưới 1.20% và độ nhám mặt bên đạt mức 2.50 µm.
  • Chứng minh tính hiệu quả vượt bậc của quy hoạch thực nghiệm Box-Hunter với 20 mẫu chuẩn, tiết kiệm hơn 65% chi phí và thời gian thực nghiệm so với phương pháp truyền thống.
  • Kiểm chứng độ tin cậy của mô hình toán học thông qua hệ thống tiêu chuẩn thống kê Student (tb = 2.57) và Fisher (Fb = 4.88) đạt mức ý nghĩa 0.05.
  • Tạo nền tảng kỹ thuật vững chắc để thúc đẩy công nghệ in 3D FDM chuyển từ vai trò tạo mẫu thử nghiệm sang sản xuất chi tiết kỹ thuật hoàn chỉnh trong ngành chế tạo máy.

Kế hoạch tiếp theo là mở rộng bài toán tối ưu hóa đa mục tiêu trên các vật liệu kỹ thuật cao trong giai đoạn 6 đến 12 tháng tới. Quý độc giả và các nhà nghiên cứu quan tâm vui lòng tham khảo toàn văn luận văn thạc sĩ để khai thác chi tiết ma trận thực nghiệm và thuật toán tối ưu phục vụ nghiên cứu và sản xuất!