Tổng quan nghiên cứu

Công nghệ tạo mẫu nhanh FDM (Fused Deposition Modeling) ra đời từ năm 1998 và đã nhanh chóng trở thành một trong những phương thức sản xuất đắp dần cốt lõi trong ngành cơ khí chính xác và chế tạo thử nghiệm. Mặc dù các hệ thống thương mại của các tập đoàn quốc tế như Stratasys có thể đạt độ chính xác từ 0.09mm đến 0.0015mm, phần lớn các thiết bị FDM tự nghiên cứu và chế tạo trong nước vẫn gặp nhiều hạn chế nghiêm trọng về độ chính xác kích thước và độ nhám bề mặt. Vấn đề cốt lõi bắt nguồn từ sai số tích lũy qua các khâu: sai số động học cơ khí, sai số chuyển đổi dữ liệu file STL, sai số cắt lớp và sự biến dạng vật liệu trong quá trình đùn ép nhiệt.

Nghiên cứu này được thực hiện nhằm giải quyết bài toán nâng cao độ chính xác tạo hình cho hệ thống máy Vina FDM 2015 được thiết kế và lắp ráp tại Phòng thí nghiệm Đo lường thuộc Khoa Cơ khí, Trường Đại học Bách Khoa - Đại học Quốc gia Thành phố Hồ Chí Minh. Mục tiêu cụ thể của đề tài là khảo sát toàn diện các nguồn sai số, thực hiện hiệu chỉnh cơ khí hệ thống 3 trục tọa độ và tối ưu hóa các thông số công nghệ chủ chốt bao gồm chiều dày lớp cắt, đường kính đầu đùn, tốc độ đùn và tốc độ chạy máy.

Phạm vi thực nghiệm được triển khai từ ngày 19/01/2015 đến ngày 20/11/2015, tập trung đánh giá mẫu thử tiêu chuẩn có kích thước danh nghĩa 25mm x 25mm trên vật liệu nhựa ABS. Kết quả nghiên cứu có ý nghĩa quyết định trong việc thiết lập quy trình calip chuẩn cho máy in 3D nội địa, kiểm soát sai số kích thước dưới ngưỡng mục tiêu 0.05mm, đồng thời giảm thiểu độ nhám bề mặt Ra, tạo nền tảng vững chắc cho ứng dụng chế tạo khuôn mẫu và y học kỹ thuật cao tại Việt Nam.

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 xây dựng trên nền tảng lý thuyết tạo mẫu nhanh FDM theo nguyên lý bồi đắp từng lớp (layer-by-layer) từ vật liệu nhiệt dẻo nóng chảy kết hợp với lý thuyết sai số động học máy công cụ 3 trục. Hệ thống động học chuyển động thẳng trong không gian 3 chiều chịu sự chi phối của 21 biến số sai số độc lập, bao gồm 18 thành phần sai số tịnh tiến và sai số góc (Roll, Pitch, Yaw) cùng 3 sai số vuông góc giữa các trục X, Y và Z. Khi các trục chuyển động mang theo đầu đùn dọc hành trình 100mm, các sai số tích lũy này sẽ trực tiếp làm sai lệch vị trí lắng đọng vật liệu.

Bên cạnh đó, đề tài ứng dụng lý thuyết xử lý dữ liệu hình học máy tính để phân tích chuỗi sai số chuyển đổi từ mô hình CAD sang định dạng lưới tam giác STL. Các khái niệm chính bao gồm:

  • Chiều cao dây cung (Chordal Tolerance): Khoảng cách lớn nhất giữa bề mặt cong thực tế và bề mặt đa giác xấp xỉ, quyết định độ mịn của dữ liệu hình học.
  • Góc điều khiển (Angle Control): Thông số kiểm soát góc lệch của các vector chỉ phương pháp tuyến tam giác nội suy.
  • Hiệu ứng bậc thang (Staircase Effect): Sai số hình học tất yếu theo phương thẳng đứng Z do quá trình phân lớp phẳng gây ra.
  • Hiện tượng biến dạng sợi đùn (Offset đường biên): Sự giãn nở và co bẹp của sợi nhựa nóng chảy đường kính 2mm thành kích thước thực tế sau khi ép qua lỗ vòi phun dưới tác động của áp lực đùn.
  • Độ nhám bề mặt Ra: Chỉ tiêu đo lường biên độ nhấp nhô tế vi của bề mặt chi tiết sau khi hoàn thành bồi đắp.

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 từ các chuỗi gia công mẫu thử trên máy Vina FDM 2015 sử dụng vật liệu sợi nhựa ABS chuyên dụng. Để đánh giá và tối ưu hóa đa thông số một cách khoa học, nghiên cứu áp dụng phương pháp quy hoạch thực nghiệm theo mô hình bề mặt đáp ứng RSM (Response Surface Methodology) với cấu trúc ma trận Box-Behnken bậc hai.

Cỡ mẫu được thiết lập gồm 27 mẫu thực nghiệm tương ứng với các tổ hợp biến thiên của 4 nhân tố công nghệ độc lập: chiều dày lớp cắt, đường kính đầu đùn, tốc độ đùn và tốc độ chạy máy. Phương pháp chọn mẫu trực giao ngẫu nhiên hóa hoàn toàn giúp loại trừ các tác động gây nhiễu từ nhiệt độ môi trường và sự dao động nguồn cấp.

Quá trình đo kiểm được thực hiện bằng hệ thống thiết bị đo chuẩn cấp chính xác cao: thước Eke vuông góc chuẩn có độ vuông góc cạnh ngoài 0.012mm/100mm và cạnh trong 0.009mm/100mm, kết hợp đồng hồ so điện tử Mitutoyo có độ phân giải 0.001mm để kiểm tra độ song song và vuông góc cơ khí. Độ nhám bề mặt Ra được kiểm tra bằng máy đo biên dạng tiếp xúc chuyên dụng.

Phần mềm Minitab 16 được lựa chọn để phân tích hồi quy phi tuyến tính và giải bài toán tối ưu hóa đa mục tiêu. Phương pháp này cho phép nhận diện định lượng mức độ ảnh hưởng của từng nhân tố đơn lẻ, đánh giá tương tác bậc cao giữa các thông số và loại bỏ các hệ số không có ý nghĩa thống kê ở mức độ tin cậy 95%, đảm bảo mô hình toán học đạt độ tương thích cao nhất.

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

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

Quá trình nghiên cứu và thực nghiệm đã mang lại 4 phát hiện khoa học quan trọng:

Thứ nhất, việc hiệu chỉnh cơ khí giúp loại bỏ phần lớn sai số lắp ráp ban đầu của hệ thống 3 trục. Sau 10 lần đo kiểm và căn chỉnh liên tục bằng đồng hồ so điện tử trên cự ly 100mm, độ vuông góc giữa trục X và trục Y đã giảm mạnh từ mức ban đầu 0.029mm/100mm xuống còn 0.009mm/100mm (mức cải thiện đạt khoảng 68.9%). Độ vuông góc giữa trục X và trục Z đạt 0.02mm/100mm, trong khi độ vuông góc giữa trục Y và trục Z được kiểm soát ở mức 0.07mm/100mm.

Thứ hai, chiều dày lớp cắt và đường kính đầu đùn là hai thông số có tác động chi phối mạnh mẽ nhất đến chất lượng chi tiết. Chiều dày lớp cắt đóng góp hơn 50% vào sự biến thiên của độ nhám bề mặt Ra. Khi chiều dày lớp cắt giảm xuống, độ nhám Ra giảm đáng kể nhưng thời gian gia công tăng lên theo tỷ lệ phi tuyến tính.

Thứ ba, sự tương tác giữa tốc độ đùn vật liệu và tốc độ chạy máy tạo ra sự biến dạng hình học không đối xứng của sợi nhựa. Khi vật liệu ABS bị ép qua đầu đùn, sợi nhựa bị co theo trục Z và phình nở theo mặt phẳng X-Y, dẫn đến kích thước chi tiết thực tế thường bị dương so với thiết kế CAD nếu không cài đặt giá trị offset đường biên chính xác.

Thứ tư, mô hình hồi quy bậc hai xây dựng trên Minitab 16 đã xác lập thành công bộ thông số công nghệ tối ưu hóa đồng thời cho 3 hàm mục tiêu: kích thước phương X tiến tới 25.00mm, kích thước phương Y tiến tới 25.00mm và độ nhám bề mặt Ra đạt giá trị cực tiểu.

Thảo luận kết quả

Nguyên nhân cốt lõi của sự sai lệch kích thước và độ nhám bề mặt bắt nguồn từ bản chất nhiệt cơ học của polymer nóng chảy khi kết hợp với chuyển động cơ học nhiều bậc tự do. Hiện tượng phình sợi đùn xuất hiện do ứng suất cắt dư và sự giải phóng áp lực khi nhựa rời khỏi vòi phun, làm thay đổi bề rộng đường bồi đắp thực tế so với đường kính lỗ đùn danh nghĩa.

Khi so sánh với các nghiên cứu quốc tế tiêu biểu của Sood và cộng sự (2009) hay Anitha và cộng sự (2001), kết quả của luận văn hoàn toàn đồng nhất về quy luật ảnh hưởng của chiều dày lớp cắt đến độ nhám Ra và tính phi tuyến tính của tốc độ đùn. Tuy nhiên, điểm cải tiến vượt trội của đề tài là đã tích hợp thành công việc bù trừ sai số cơ khí thực nghiệm trực tiếp với thuật toán tối ưu hóa phần mềm, giúp máy Vina FDM 2015 đạt được cấp chính xác tương đương các thiết bị thương mại cùng phân khúc.

Dữ liệu thực nghiệm được trực quan hóa rõ nét thông qua biểu đồ bề mặt đáp ứng 3D (Surface Plots) và biểu đồ đường mức (Contour Plots) trong Minitab 16. Các biểu đồ này mô tả sinh động vùng làm việc tối ưu, nơi các đường đồng mức của kích thước X và Y giao nhau tại điểm hội tụ sai số bằng 0, đồng thời phản ánh rõ nét mối quan hệ đánh đổi giữa tốc độ in và độ bóng bề mặt chi tiết.

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

Để duy trì và nâng cao hơn nữa độ chính xác của hệ thống máy tạo mẫu nhanh FDM trong sản xuất thực tế, luận văn đề xuất 4 nhóm giải pháp kỹ thuật cụ thể:

  1. Chuẩn hóa quy trình calip động học cơ khí định kỳ: Áp dụng phương pháp gá Eke chuẩn và đồng hồ so điện tử kiểm tra độ vuông góc các trục định kỳ 3 tháng một lần; mục tiêu duy trì sai số vuông góc trục X-Y dưới 0.010mm/100mm và trục Z dưới 0.025mm/100mm; do đội ngũ kỹ sư vận hành và bảo trì thực hiện xuyên suốt vòng đời máy.

  2. Ứng dụng giải thuật cắt lớp thích nghi (Adaptive Slicing): Nghiên cứu tích hợp thuật toán phân lớp tự động điều chỉnh chiều dày lớp cắt theo độ dốc biên dạng của chi tiết vào phần mềm CAM; mục tiêu giảm 25% tổng thời gian tạo mẫu mà vẫn đảm bảo độ nhám Ra tại các bề mặt nghiêng đạt dưới 3.2 micromet; do nhóm kỹ sư phần mềm điều khiển triển khai trong vòng 6 tháng.

  3. Tích hợp bảng thông số công nghệ chuẩn cho từng loại vật liệu: Thiết lập cơ sở dữ liệu thông số tối ưu cho vật liệu ABS và PLA dựa trên mô hình bề mặt đáp ứng Minitab; kiểm soát chặt chẽ tỷ lệ đùn và offset đường biên nhằm duy trì dung sai kích thước sản phẩm trong khoảng cộng trừ 0.03mm; do kỹ thuật viên phòng thí nghiệm áp dụng ngay trong các quy trình in thực tế.

  4. Nâng cấp hệ thống kiểm soát môi trường tạo mẫu: Trang bị buồng in gia nhiệt kín duy trì nhiệt độ ổn định ở mức 50 độ C và lắp đặt tủ sấy chuyên dụng để kiểm soát độ ẩm của cuộn dây nhựa ABS dưới 0.05%; mục tiêu triệt tiêu hiện tượng co ngót nhiệt và rỗ khí sợi đùn; do ban quản trị xưởng chế tạo đầu tư thực hiện trong quý tới.

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

Công trình nghiên cứu cung cấp khối lượng dữ liệu thực nghiệm và giải pháp kỹ thuật phong phú, đặc biệt hữu ích cho 4 nhóm đối tượng sau:

  1. Học viên cao học và nghiên cứu sinh chuyên ngành Cơ khí - Tự động hóa: Luận văn cung cấp phương pháp luận chặt chẽ về phân tích sai số động học 21 bậc tự do, phương pháp quy hoạch thực nghiệm Box-Behnken và cách xử lý dữ liệu tối ưu hóa bề mặt đáp ứng trên phần mềm Minitab.

  2. Kỹ sư thiết kế và chế tạo máy in 3D: Nắm vững quy trình căn chỉnh cơ khí thực tế bằng dụng cụ đo tiêu chuẩn, từ đó nâng cao độ chính xác lắp ráp khung máy và tối ưu hóa hệ dẫn động trục X-Y-Z với chi phí gia công thấp nhất.

  3. Doanh nghiệp sản xuất mẫu nhanh và công nghệ khuôn mẫu: Ứng dụng bộ thông số công nghệ tối ưu để rút ngắn 20% thời gian thử nghiệm sản phẩm mới, giảm tỷ lệ phế phẩm kích thước xuống dưới 2% và nâng cao chất lượng thẩm mỹ của chi tiết nhựa mẫu.

  4. Giảng viên và cán bộ giảng dạy tại các trường đại học kỹ thuật: Sử dụng tài liệu như một nghiên cứu điển hình (case study) mẫu mực phục vụ giảng dạy các môn học Đo lường kỹ thuật, Công nghệ CAD/CAM/CNC và Thiết kế thực nghiệm công nghiệp.

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

Sau quá trình tối ưu hóa, độ chính xác kích thước của máy Vina FDM 2015 đạt được là bao nhiêu? Sau khi hoàn thành cân chỉnh cơ khí và cài đặt bộ thông số công nghệ tối ưu từ mô hình Minitab 16, máy Vina FDM 2015 đạt độ sai lệch kích thước dưới 0.05mm trên khối mẫu chuẩn 25mm x 25mm, tiệm cận tiêu chuẩn chính xác của các dòng máy FDM thương mại quốc tế.

Tại sao phải hiệu chuẩn độ vuông góc giữa các trục chuyển động trước khi tối ưu hóa thông số in? Hệ thống chuyển động 3 trục có đến 21 thành phần sai số động học. Nếu các trục X, Y, Z không vuông góc chuẩn, sai số góc sẽ tích lũy liên tục theo hành trình 100mm, gây méo biên dạng hình học của sản phẩm mà không có phần mềm cắt lớp nào có thể bù trừ hoàn toàn được.

Chiều dày lớp cắt ảnh hưởng như thế nào đến chất lượng chi tiết in 3D? Chiều dày lớp cắt là yếu tố quyết định hơn 50% độ nhám bề mặt Ra do hiệu ứng bậc thang trên bề mặt cong. Giảm chiều dày lớp cắt giúp bề mặt chi tiết mịn màng và biên dạng sắc nét hơn nhưng sẽ làm tăng đáng kể thời gian bồi đắp và chi phí vận hành máy.

Nguyên nhân nào khiến sợi vật liệu đùn ra bị thay đổi kích thước so với đường kính đầu phun? Khi sợi nhựa ABS đường kính 2mm được nung chảy và đẩy qua lỗ vòi phun dưới áp lực đùn, vật liệu chịu lực nén cơ học khi bồi đắp lên lớp trước. Hiện tượng này làm sợi nhựa co theo trục Z và phình rộng theo mặt phẳng X-Y, đòi hỏi phải cài đặt thông số offset đường biên phù hợp.

Mô hình thực nghiệm Box-Behnken mang lại lợi ích gì so với phương pháp thử - sai truyền thống? Phương pháp Box-Behnken cho phép khảo sát toàn diện 4 yếu tố công nghệ chỉ với 27 lần chạy thực nghiệm thay vì hàng trăm lần thử nghiệm đơn lẻ. Mô hình thiết lập phương trình hồi quy bậc hai chính xác, giúp phát hiện các tương tác phi tuyến và xác định điểm tối ưu toàn cục nhanh chóng.

Kết luận

• Luận văn đã phân tích toàn diện 21 thành phần sai số động học và xây dựng thành công quy trình hiệu chuẩn cơ khí cho máy Vina FDM 2015, đưa độ vuông góc trục X-Y đạt mức chuẩn 0.009mm/100mm. • Khảo sát chi tiết chuỗi sai số thiết lập file STL và cơ chế biến dạng hình học của sợi đùn polymer nhiệt dẻo trong quá trình tạo hình từng lớp. • Ứng dụng xuất sắc ma trận quy hoạch thực nghiệm Box-Behnken với 27 mẫu thử để xác lập các phương trình hồi quy mô tả kích thước X, Y và độ nhám Ra. • Tối ưu hóa thành công tổ hợp thông số công nghệ trên phần mềm Minitab 16, kiểm soát sai số kích thước sản phẩm dưới 0.05mm và cực tiểu hóa độ nhám bề mặt. • Khẳng định tính khả thi của việc chế tạo và làm chủ máy in 3D công nghệ cao trong nước với chi phí thấp nhưng đạt độ chính xác tương đương thiết bị ngoại nhập.

Đóng góp lớn nhất của công trình là cung cấp một phương pháp luận chuẩn xác kết hợp giữa hiệu chuẩn phần cứng cơ khí và tối ưu hóa phần mềm điều khiển. Trong lộ trình 12 tháng tiếp theo, nhóm nghiên cứu khuyến nghị tích hợp thuật toán cắt lớp thích nghi tự động và hoàn thiện hệ thống cảm biến nhiệt buồng kín. Hãy khai thác ngay bộ thông số và quy trình tối ưu này để nâng tầm chất lượng cho các dự án in 3D và sản xuất mẫu nhanh của bạn!