Phương pháp thiết kế thí nghiệm cho phát triển quy trình in 3D: từ sàng lọc trực giao, mô hình hóa bề mặt đáp ứng đến hiệu chuẩn cửa sổ quy trình

Giới thiệu: Vì sao tối ưu tham số luôn giậm chân tại chỗ

Trong quá trình hiệu chỉnh tại hiện trường đối với nung chảy bột kim loại bằng laser (LPBF) và thiêu kết polymer, tình huống phổ biến nhất là: kỹ sư tăng công suất laser từ 250 W lên 280 W, độ đặc tăng 0,3%, nên tiếp tục tăng; đến 300 W xuất hiện hiện tượng lỗ khóa, rồi lại phải điều chỉnh giảm xuống. Sau một vòng thực hiện hơn chục mẫu thử, kết quả chỉ là một đường cong kinh nghiệm rời rạc; khi thay máy hoặc thay lô bột, mọi việc phải làm lại từ đầu. Vấn đề không phải kỹ sư thiếu nỗ lực, mà là phương pháp thay đổi luân phiên từng yếu tố vốn không thể xử lý hiệu ứng tương tác — công suất và tốc độ quét cùng quyết định mật độ năng lượng thể tích, còn chiều dày lớp rải bột lại thay đổi điều kiện tản nhiệt của bể nóng chảy; ba yếu tố liên kết với nhau, nên thay đổi từng biến một sẽ không bao giờ cho thấy toàn cảnh.

Thiết kế thí nghiệm (Design of Experiments, DoE) chính là con đường kỹ thuật để giải quyết các vấn đề ghép nối đa yếu tố này: với số lượng mẫu ít nhất, có thể đồng thời ước lượng hiệu ứng chính, hiệu ứng tương tác và bề mặt đáp ứng của các yếu tố, đồng thời xác định cửa sổ quy trình có thể phê duyệt kèm khoảng tin cậy. Bài viết này hướng đến bối cảnh phát triển quy trình in 3D, trình bày phương pháp hoàn chỉnh và có thể triển khai từ sàng lọc trực giao, mô hình hóa bề mặt đáp ứng đến xác minh cửa sổ quy trình.

Một, trước hết cần phân biệt ba loại thí nghiệm: sàng lọc, mô hình hóa và xác minh

Phát triển quy trình không thể hoàn tất chỉ bằng một thí nghiệm mà thường được triển khai theo ba tầng, trong đó mỗi tầng có mục đích, số lượng yếu tố và số lượng mẫu khác nhau.

Tầng thứ nhất là thí nghiệm sàng lọc, nhằm chọn ra 2–4 yếu tố thực sự có ý nghĩa trong số 6–10 yếu tố ứng viên. Các yếu tố ứng viên điển hình gồm công suất laser, tốc độ quét, khoảng cách quét, chiều dày lớp, nhiệt độ gia nhiệt trước, chiến lược quét (bề rộng dải), góc xoay giữa các lớp và lưu lượng khí bảo vệ. Giai đoạn sàng lọc sử dụng thiết kế nhân tố phân đoạn hoặc thiết kế Plackett-Burman; chỉ cần 8–12 lượt thí nghiệm, với hai chỉ tiêu đáp ứng thô là độ đặc và độ nhám bề mặt.

Tầng thứ hai là thí nghiệm mô hình hóa, nhằm xây dựng mô hình bề mặt đáp ứng cho 2–3 yếu tố then chốt đã được sàng lọc. Thiết kế hỗn hợp trung tâm (CCD) hoặc thiết kế Box-Behnken thường được sử dụng. Một CCD ba yếu tố gồm 6 điểm trục, 8 điểm lập phương và 6 lần lặp tại tâm, tổng cộng 20 lượt thí nghiệm. Các chỉ tiêu đáp ứng được mở rộng thành bốn hạng mục: độ đặc, độ bền kéo, sai lệch kích thước và độ nhám bề mặt Ra.

Tầng thứ ba là thí nghiệm xác minh, trong vùng tối ưu do mô hình dự báo, lựa chọn 3–5 điểm xác nhận để in lặp lại và kiểm tra xem sai lệch giữa giá trị dự báo của mô hình với giá trị đo thực tế có nằm trong ±5% hay không. Sau khi xác minh đạt yêu cầu, thu hẹp biên cửa sổ 10% để xác định phạm vi tham số cho phiếu quy trình sản xuất hàng loạt.

Hai, thí nghiệm trực giao: bao phủ không gian lớn nhất với số mẫu ít nhất

Khi số lượng yếu tố không nhiều và số mức của các yếu tố bằng nhau, bảng trực giao là lựa chọn kinh tế nhất. Lấy LPBF thép không gỉ 316L làm ví dụ, xét ba yếu tố ba mức gồm công suất laser P, tốc độ quét v và khoảng cách quét h. Thiết kế toàn nhân tố cần 27 lượt thí nghiệm, trong khi sử dụng bảng trực giao L9(3^4) chỉ cần 9 lượt.

Thiết kế cụ thể: P là 200/240/280 W, v là 800/1000/1200 mm/s, h là 0,09/0,11/0,13 mm. Sau khi bố trí thí nghiệm theo L9, in một nhóm mẫu lập phương 10×10×10 mm cho mỗi tổ hợp mức, đo độ đặc bằng phương pháp Archimedes và đồng thời ghi nhận độ nhám bề mặt Ra.

Có thể nhanh chóng xếp hạng bằng phân tích khoảng biến thiên: tính đáp ứng trung bình của mỗi yếu tố tại từng mức, với khoảng biến thiên R = max(đáp ứng trung bình) − min(đáp ứng trung bình); R càng lớn thì ảnh hưởng của yếu tố càng đáng kể. Nếu thu được R(P)=2,1%, R(v)=1,4%, R(h)=0,6%, thứ tự ảnh hưởng là công suất > tốc độ > khoảng cách. Tiếp tục thực hiện phân tích phương sai (ANOVA), sử dụng kiểm định F để đánh giá ý nghĩa: các yếu tố có giá trị F lớn hơn giá trị tới hạn và p nhỏ hơn 0,05 sẽ được đưa vào vòng mô hình hóa tiếp theo.

Hạn chế của thí nghiệm trực giao là giả định không có tương tác giữa các yếu tố. Nếu nghi ngờ công suất và tốc độ có ghép nối mạnh với nhau (mật độ năng lượng thể tích E = P/(v·h·t) vốn là quan hệ tích), nên chuyển sang thiết kế bề mặt đáp ứng trong giai đoạn mô hình hóa thay vì tiếp tục áp dụng cứng nhắc bảng trực giao.

Ba, phương pháp bề mặt đáp ứng: từ các điểm rời rạc đến bề mặt liên tục

Cốt lõi của phương pháp bề mặt đáp ứng (RSM) là dùng đa thức bậc hai để khớp quan hệ giữa đáp ứng và các yếu tố:

Y = β0 + Σβi·xi + Σβii·xi² + Σβij·xi·xj + ε

Trong đó, hạng chéo βij là biểu diễn định lượng của hiệu ứng tương tác. Lấy thiết kế CCD làm ví dụ, tham số then chốt của phương án ba yếu tố là khoảng cách trục α. Với thiết kế có tính quay, α = 2^(k/4) = 2^(3/4) ≈ 1,682 (k là số yếu tố). Các điểm trục nằm tại ±1,682 trong không gian mã hóa, điểm lập phương nằm tại ±1, còn điểm tâm được lặp 6 lần để ước lượng sai số thuần.

Sau khi khớp mô hình, cần thực hiện ba kiểm tra: thứ nhất là kiểm định thiếu phù hợp (lack-of-fit); p lớn hơn 0,05 cho thấy dạng mô hình không sai khác đáng kể so với dữ liệu. Thứ hai là hệ số xác định R² và R² hiệu chỉnh; chênh lệch giữa hai giá trị nên nhỏ hơn 0,2, nếu không có thể xảy ra quá khớp. Thứ ba là kiểm tra tính chuẩn của phần dư và tính đồng nhất phương sai. Trong thực tế, trường hợp thường gặp là R² đạt 0,95 nhưng R² hiệu chỉnh chỉ là 0,86; khi đó nên loại bỏ các hạng bậc ba không có ý nghĩa rồi khớp lại mô hình.

Tối ưu hóa đa đáp ứng là một thách thức trong in 3D: độ đặc cần tối đa, độ nhám cần tối thiểu, còn sai lệch kích thước cần gần bằng không; ba mục tiêu này thường xung đột với nhau. Trong kỹ thuật, hàm mong muốn (desirability function) được dùng để ánh xạ từng đáp ứng vào khoảng 0–1, sau đó lấy trung bình hình học để có độ mong muốn tổng hợp D và xác định tổ hợp tham số tại vị trí D lớn nhất. Trọng số được phân bổ theo mức độ ưu tiên trong nghiệm thu của khách hàng; ví dụ, với chi tiết hàng không, có thể đặt trọng số độ đặc là 0,5, độ nhám là 0,3 và sai lệch kích thước là 0,2.

Bốn, hiệu chuẩn cửa sổ quy trình: biến điểm tối ưu thành khoảng có thể phê duyệt

Tìm được điểm tối ưu mới chỉ là bước đầu; sản xuất hàng loạt cần một khoảng vận hành. Cách thực hiện phù hợp là dùng biểu đồ đường đồng mức chồng lấp (overlay contour plot): vẽ đường đồng mức cho ngưỡng đạt yêu cầu của từng đáp ứng (độ đặc ≥ 99,5%, Ra ≤ 12 μm, sai lệch kích thước ≤ ±0,1 mm), rồi lấy vùng giao nhau của tất cả các đường đồng mức làm cửa sổ quy trình.

Sau khi xác định cửa sổ, bắt buộc phải chừa biên độ an toàn. Khuyến nghị thu hẹp biên cửa sổ vào trong 10%; ví dụ, nếu mô hình cho khoảng công suất khả thi là 235–265 W thì phiếu quy trình ghi 238–262 W. Lý do là biến động giữa các lô bột, sự trôi độ tinh khiết khí bảo vệ và suy hao công suất của hệ quang học đều có thể khiến điểm quy trình thực tế dịch chuyển dần; chừa biên độ giúp nâng chỉ số năng lực quy trình Cpk từ 1,0 lên trên 1,33.

Xác minh cửa sổ áp dụng ba lô xác nhận: mỗi lô 5 chi tiết, lần lượt in một nhóm tại tâm cửa sổ, sát phía trong biên trên và sát phía trong biên dưới. Nếu cả ba lô liên tiếp đều đạt yêu cầu và độ trôi của giá trị trung bình giữa các lô nhỏ hơn 1/3 độ lệch chuẩn của đáp ứng, có thể xác định cửa sổ ổn định. Dữ liệu xác minh cũng phải được nhập vào sổ theo dõi quy trình làm đường cơ sở cho việc truy xuất sau này.

Năm, những sai lầm thường gặp và danh sách kiểm tra

Sai lầm thứ nhất: coi số lượng mẫu thử là toàn bộ chi phí. Trên thực tế, 9 lượt thí nghiệm trực giao cộng với 20 lượt CCD có tổng chi phí mẫu thấp hơn nhiều so với 30 lượt điều chỉnh từng yếu tố một cách mù quáng, đồng thời cho ra một mô hình có thể tái sử dụng.

Sai lầm thứ hai: bỏ qua hiệu ứng khối. Nếu 20 lượt thí nghiệm CCD được thực hiện trong ba ngày, những thay đổi về trạng thái thiết bị, lô bột, nhiệt độ và độ ẩm môi trường sẽ lẫn vào phần dư. Cách làm đúng là đưa ngày thực hiện vào thiết kế như một yếu tố khối, hoặc bố trí lặp điểm tâm trong từng khối để theo dõi độ trôi.

Sai lầm thứ ba: chỉ xem độ đặc. Điểm quy trình đạt độ đặc 99,5% chưa chắc đạt yêu cầu về tính chất cơ học; khuyết tật thiếu nóng chảy có thể bị đánh giá thấp trong phép đo mật độ. Khuyến nghị tối thiểu phải đồng thời kiểm tra độ bền kéo và độ giãn dài; khi cần, bổ sung kiểm tra kim tương vi mô và kiểm tra lấy mẫu bằng CT công nghiệp.

Danh sách kiểm tra triển khai: các mức yếu tố có bao phủ phạm vi vận hành ổn định của thiết bị hay không; có bố trí lặp điểm tâm để ước lượng sai số thuần hay không; đã thực hiện kiểm định thiếu phù hợp hay chưa; trọng số đa đáp ứng có nhất quán với tiêu chí nghiệm thu của khách hàng hay không; cửa sổ có chừa biên độ 10% hay không; các lô xác minh có bao phủ biên cửa sổ hay không.

Tổng kết

Thiết kế thí nghiệm chuyển hoạt động phát triển quy trình in 3D từ thử-sai dựa trên kinh nghiệm thành hoạt động kỹ thuật có thể định lượng và tái lập: thí nghiệm trực giao hoàn thành sàng lọc yếu tố trong 9–12 lượt thử, phương pháp bề mặt đáp ứng xây dựng mô hình bậc hai có hiệu ứng tương tác với khoảng 20 lượt thử, hàm mong muốn giải quyết xung đột đa mục tiêu, còn đường đồng mức chồng lấp đưa ra cửa sổ quy trình có thể phê duyệt với biên độ an toàn. Đối với các nền tảng quản lý vòng đời toàn diện như lantu3D, vốn cần đi từ thiết kế bản vẽ đến bàn giao sản phẩm thực, giá trị của DoE không chỉ nằm ở điều chỉnh tham số — nó tích lũy tri thức quy trình thành các mô hình tham số có thể truy xuất, tái sử dụng và kiểm chứng, để mỗi lần phát triển quy trình trở thành điểm xuất phát cho lần tiếp theo thay vì lại bắt đầu từ con số không.