Tổng quan về luận án
Công nghệ xử lý và biến tính bề mặt đóng vai trò cốt lõi trong kỹ thuật cơ khí hiện đại, đặc biệt là chiến lược phục hồi và gia cường các chi tiết máy chịu tải trọng khắc nghiệt, mài mòn cao, xâm thực hóa học và sốc nhiệt. Trong số các kỹ thuật phun phủ nhiệt hiện đại, phương pháp phun plasma trong không khí (Atmospheric Plasma Spraying - APS) nổi bật như một giải pháp công nghệ tiên tiến nhờ tạo ra cột hồ quang plasma có mật độ năng lượng cực lớn với nhiệt độ vùng trung tâm đạt từ 15.000°C đến 20.000°C. Nguồn nhiệt năng lượng cao này cho phép nóng chảy hoàn toàn các loại vật liệu gốm ôxít có nhiệt độ nóng chảy cao và trơ hóa học, điển hình là hệ vật liệu gốm composite $Al_2O_3 – TiO_2$.
Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực trạng phần lớn các công bố quốc tế trước đây chỉ tập trung khảo sát đơn lẻ ảnh hưởng của các thông số công nghệ đến từng chỉ tiêu cơ tính riêng biệt (hoặc độ cứng, hoặc độ xốp, hoặc độ bền mài mòn) hoặc chỉ tối ưu hóa cục bộ với hàm đơn mục tiêu. Cụ thể, các nghiên cứu của Yusoff cùng cộng sự (2012) hay Islak cùng cộng sự (2014) chủ yếu phân tích ảnh hưởng của công suất nhiệt hoặc tỷ lệ thành phần $TiO_2$ đến hình thái cấu trúc mà chưa làm rõ mối quan hệ tương tác đồng thời, phi tuyến giữa bộ thông số công nghệ công nghệ chính gồm: Cường độ dòng điện phun ($I$), Khoảng cách phun ($L$) và Lưu lượng cấp bột ($M$) đối với một tập hợp đầy đủ các chỉ tiêu cơ tính phức hợp.
Mô hình nghiên cứu giải quyết 4 câu hỏi nghiên cứu ($RQ$) và kiểm định 4 giả thuyết khoa học ($H$):
- RQ1: Quy luật biến thiên và mức độ ảnh hưởng định lượng của dòng điện phun ($I$), khoảng cách phun ($L$) và lưu lượng cấp bột ($M$) đến cấu trúc vi mô và sự chuyển biến pha của lớp phủ gốm $Al_2O_3 – 40%TiO_2$ diễn ra như thế nào?
- RQ2: Mối tương quan toán học giữa ba thông số công nghệ ($I, L, M$) với 4 chỉ tiêu cơ tính đầu ra (Độ cứng tế vi $HV$, Hệ số ma sát khô $\mu$, Độ bền bám dính trượt $\tau$, Độ xốp $P$) được mô hình hóa chính xác bằng hàm hồi quy nào?
- RQ3: Mức độ ưu tiên và tỷ lệ đóng góp phần trăm ($P_j%$) của từng thông số công nghệ đến từng chỉ tiêu cơ tính đơn lẻ thông qua phân tích phương sai (ANOVA) đạt giá trị bao nhiêu?
- RQ4: Làm thế nào để xác định được một bộ thông số công nghệ tối ưu duy nhất thỏa mãn đồng thời cả 4 chỉ tiêu cơ tính trái ngược nhau (cực đại hóa độ cứng, cực đại hóa độ bền bám dính, cực tiểu hóa độ xốp và cực tiểu hóa hệ số ma sát)?
Các giả thuyết nghiên cứu tương ứng:
- H1: Cường độ dòng điện phun ($I$) là yếu tố chi phối năng lượng enthalpy của dòng plasma, quyết định trực tiếp đến mức độ nóng chảy triệt để của hạt gốm và tỷ lệ pha chuyển biến $\gamma-Al_2O_3$.
- H2: Khoảng cách phun ($L$) kiểm soát thời gian bay của hạt, vận tốc va đập và mức độ làm nguội trước khi hạt va chạm với bề mặt nền thép C45.
- H3: Mối quan hệ giữa các thông số công nghệ ($I, L, M$) và các chỉ tiêu cơ tính mang tính phi tuyến bậc hai có thể biểu diễn qua các hàm hồi quy thực nghiệm với hệ số xác định $R^2 \ge 0,90$.
- H4: Việc áp dụng chỉ số đánh giá tổng thể (Overall Evaluation Criteria - OEC) kết hợp quy hoạch thực nghiệm Taguchi cho phép dung hòa các xung đột mục tiêu cơ tính, xác định nghiệm tối ưu toàn cục có độ lặp lại cao.
Khung lý thuyết của nghiên cứu tích hợp lý thuyết nhiệt động học dòng hồ quang plasma ion hóa khí trơ ($Ar/H_2$), lý thuyết biến dạng dẻo và đông đặc siêu nhanh hạt va đập của Schoop, Pospisil - Sehyl, kết hợp với lý thuyết tối ưu hóa thực nghiệm Taguchi và phương pháp đánh giá tổng thể OEC. Đóng góp đột phá của công trình được định lượng thông qua việc làm chủ công nghệ phun plasma cho hệ gốm $Al_2O_3 – 40%TiO_2$, tạo ra lớp phủ có độ cứng tế vi đạt trên $850\text{ }HV$, độ bám dính trượt vượt trội, độ xốp được kiểm soát dưới $2,5%$, giúp giảm giá thành chế tạo và phục hồi chi tiết máy xuống chỉ còn $15% – 50%$ so với chi phí nhập khẩu mới nguyên chiếc.
Phạm vi thực nghiệm được triển khai trên nền thép kết cấu cacbon C45, sử dụng lớp lót trung gian hợp kim $NiCr20$ và lớp phủ gốm composite $Al_2O_3 – 40%TiO_2$ (kích thước hạt $15 – 45\text{ }\mu m$) với mảng trực giao Taguchi $L_9(3^3)$, kết hợp hệ thống kiểm chứng độc lập đảm bảo độ tin cậy tuyệt đối về mặt đo lường và khoa học vật liệu.
Literature Review và Positioning
Lịch sử phát triển của công nghệ phun phủ nhiệt đã chứng kiến sự tiến hóa vượt bậc từ những quan sát ban đầu của Max Ulrich Schoop vào năm 1910 về sự va đập và biến dạng dẻo của các giọt kim loại nóng chảy rời khỏi nòng súng phun hồ quang điện sơ khai. Đến thập niên 1960, Gage đã tạo ra bước ngoặt mang tính cách mạng khi phát minh ra phương pháp phun plasma trong không khí (Atmospheric Plasma Spraying - APS), mở ra khả năng xử lý các loại vật liệu có nhiệt độ nóng chảy siêu cao mà ngọn lửa khí cháy thông thường (PFS, WFS) không thể thực hiện được. Sau đó, Muehlberger (1974) phát triển kỹ thuật phun plasma chân không (Vacuum Plasma Spraying - VPS) nhằm hạn chế quá trình ôxít hóa, tuy nhiên kỹ thuật này đòi hỏi chi phí đầu tư thiết bị và vận hành cực kỳ đắt đỏ.
Tổng hợp các dòng nghiên cứu quốc tế lớn về hệ gốm $Al_2O_3 – TiO_2$ cho thấy nhiều tranh luận học thuật sâu sắc:
- Dòng nghiên cứu về thành phần và kích thước hạt: Jordan cùng các cộng sự (2000) và You Wang cùng các cộng sự (1999) đã so sánh lớp phủ sử dụng bột gốm kích thước tiêu chuẩn (Metco-130) với bột gốm cấu trúc nano. Kết quả chỉ ra rằng cấu trúc nano làm tăng khả năng chống mài mòn lên gấp 4 lần nhờ hạn chế sự phát triển của các vết nứt tế vi và cải thiện mật độ pha. Ngược lại, Islak cùng cộng sự (2014) khi khảo sát bột micro với các tỷ lệ $TiO_2$ khác nhau ($3%, 13%, 40%$) trên nền thép SAE 1040 lại nhấn mạnh rằng việc tăng hàm lượng $TiO_2$ lên $40%$ đóng vai trò như chất liên kết pha lỏng (binder-like phase), làm giảm đáng kể nhiệt độ nóng chảy của hỗn hợp, từ đó làm giảm độ xốp và tăng độ dai va đập dù độ cứng có xu hướng giảm nhẹ so với gốm $Al_2O_3$ nguyên chất.
- Dòng nghiên cứu về cơ chế biến đổi pha và ma sát học: Kanta et al. (2010) và Xiaoqin Zhao et al. (2004) phát hiện sự chuyển pha từ $\alpha-Al_2O_3$ sang pha giả bền $\gamma-Al_2O_3$ trong quá trình nguội nhanh với tốc độ vượt quá $10^6\text{ }^\circ\text{C/s}$. Tranh luận nổ ra giữa nhóm tác giả ủng hộ việc giữ lại tối đa pha $\alpha-Al_2O_3$ để tối đa hóa độ cứng tuyệt đối và nhóm tác giả khẳng định sự xuất hiện đồng thời của hai pha $\alpha/\gamma$ cùng pha rutil/anatase của $TiO_2$ mới mang lại hệ số ma sát thấp và khả năng chống mòn tự bôi trơn tốt nhất.
Vị trí của công trình nghiên cứu được định vị vững chắc tại giao điểm giữa khoa học vật liệu bề mặt và kỹ thuật điều khiển tối ưu hóa đa thông số. Công trình giải quyết triệt để hạn chế của các nghiên cứu quốc tế đi trước bằng cách không chỉ xem xét cơ chế biến đổi cấu trúc tế vi mà còn lượng hóa chính xác sự tác động đồng thời của 3 thông số công nghệ cốt lõi ($I, L, M$) thông qua mô hình toán học giải tích bậc hai và thuật toán tối ưu đa mục tiêu OEC.
So sánh với hai nghiên cứu quốc tế điển hình:
- So với nghiên cứu của Yusoff et al. (2012) trên hệ $Al_2O_3 – 13%TiO_2$ (chỉ khảo sát biến thiên công suất từ $20\text{ }kW$ đến $40\text{ }kW$ làm giảm độ xốp từ $6,2%$ xuống $4%$), công trình này đã mở rộng sang hệ vật liệu $Al_2O_3 – 40%TiO_2$, kiểm soát độ xốp xuống mức sâu hơn đáng kể ($< 2,5%$) nhờ điều phối chính xác khoảng cách bay và lưu lượng nạp bột.
- So với nghiên cứu của Bounazef et al. (2004) trên nền thép không gỉ AISI 304 sử dụng lớp lót Ni5Al trong môi trường ăn mòn axit, công trình này thiết lập quy trình công nghệ toàn diện trên nền thép C45 với lớp lót NiCr20, đánh giá đồng bộ cả 4 đặc tính cơ - lý - ma sát học thay vì chỉ kiểm tra độ xốp và khả năng chống ăn mòn tĩnh.
Đóng góp lý thuyết và khung phân tích
Đóng góp cho lý thuyết
Nghiên cứu mang lại những đóng góp nền tảng cho lý thuyết luyện kim bề mặt và động lực học hạt trong dòng plasma nhiệt độ cao:
- Kế thừa và phát triển các lý thuyết tạo màng kinh điển: Bổ sung và hoàn thiện các luận điểm của Pospisil - Sehyl và Schoop. Trong khi lý thuyết Schoop tập trung giải thích sự biến dạng dẻo của hạt kim loại nguội nhanh, nghiên cứu này chứng minh rằng đối với hạt gốm ôxít composite phi kim loại ($Al_2O_3 – 40%TiO_2$), cơ chế hình thành lớp màng liên kết phụ thuộc vào mức độ nóng chảy hoàn toàn của pha $TiO_2$ (nhiệt độ nóng chảy $\sim 1.843^\circ\text{C}$) bao bọc và hòa tan một phần các hạt $\alpha-Al_2O_3$ (nhiệt độ nóng chảy $\sim 2.050^\circ\text{C}$). Khi va đập với tốc độ $200 – 400\text{ }m/s$, các giọt lỏng gốm trải qua quá trình làm phẳng (splat flattening) trong khoảng thời gian siêu ngắn ($10^{-6}\text{ }s$) và đông đặc với tốc độ làm nguội cực lớn ($10^6 – 10^7\text{ }^\circ\text{C/s}$), tạo nên cấu trúc dạng lớp xếp chồng (lamellar structure) đặc trưng xen kẽ các vi lỗ xốp và liên kết cơ học cài khớp (mechanical interlocking).
- Phản biện lý thuyết Schenk: Luận án chỉ rõ rằng quá trình liên kết bền vững giữa lớp phủ gốm và nền kim loại không đòi hỏi bề mặt kim loại nền phải bị nung nóng đến trạng thái nóng chảy như giả thuyết của Schenk. Ngược lại, năng lượng nhiệt cục bộ của hạt nóng chảy kết hợp với độ nhám bề mặt được tạo nhám cơ học bằng phun cát hạt thép ($R_a = 5 – 8\text{ }\mu m$) và một lớp lót hợp kim trung gian ($NiCr20$) là điều kiện đủ để tạo ra liên kết bám dính cơ - nhiệt hoàn hảo mà không làm biến tính cấu trúc tế vi của kim loại nền C45.
- Mô hình hóa giải tích quan hệ tham số - cơ tính: Xây dựng thành công hệ 4 phương trình hồi quy đa biến dạng đa thức bậc hai biểu diễn đầy đủ mối tương quan giữa vector đầu vào $\vec{X} = [I, L, M]^T$ với các hàm mục tiêu $Y_k \in {HV, \mu, \tau, P}$:
$$\begin{aligned}
HV &= f(I, L, M) \
\mu &= f(I, L, M) \
\tau &= f(I, L, M) \
P &= f(I, L, M)
\end{aligned}$$
với các hệ số hồi quy được kiểm định chặt chẽ thông qua tiêu chuẩn Fisher ($F_s$).
Khung phân tích độc đáo
Khung phân tích của nghiên cứu là sự dung hợp chặt chẽ giữa 3 trường phái lý thuyết:
- Lý thuyết nhiệt động học chất khí ion hóa và truyền nhiệt bức xạ - đối lưu trong súng phun plasma DC (Praxair SG100).
- Lý thuyết mỏi tiếp xúc và cơ học ma sát học (Tribology) giải thích cơ chế mòn trượt khốc liệt chuyển sang mòn cào xước vi mô.
- Lý thuyết thiết kế thực nghiệm vững bền (Robust Design Theory) của Genichi Taguchi kết hợp chỉ tiêu đánh giá tổng thể (Overall Evaluation Criteria - OEC).
Điều kiện biên của khung phân tích:
- Áp dụng cho hệ vật liệu bột gốm $Al_2O_3 – 40%TiO_2$ kích thước hạt danh định từ $15\text{ }\mu m$ đến $45\text{ }\mu m$.
- Khí mang bột và khí tạo plasma sử dụng hỗn hợp Argon công nghiệp độ tinh khiết cao.
- Nền thép cacbon kết cấu C45 dạng tấm phẳng và dạng trụ, được xử lý nhiệt luyện đạt độ cứng chuẩn trước khi tạo nhám và phun lớp lót trung gian $NiCr20$ dày $50 – 80\text{ }\mu m$.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng (positivism) với phương pháp tiếp cận định lượng thực nghiệm nghiêm ngặt. Hệ thống thiết kế thực nghiệm sử dụng mảng trực giao Taguchi $L_9(3^3)$ bao gồm 9 thí nghiệm đại diện cho không gian khảo sát 3 yếu tố công nghệ ở 3 mức biến thiên độc lập:
- Cường độ dòng điện phun ($I$): Khảo sát ở 3 mức tương ứng với mức năng lượng hồ quang thấp, trung bình và cao.
- Khoảng cách phun ($L$): Khảo sát ở 3 mức cự ly từ đầu vòi phun đến bề mặt nền.
- Lưu lượng cấp bột ($M$): Khảo sát ở 3 mức tốc độ cấp bột thông qua đĩa quay định lượng.
Mô hình thiết kế đa mức cho phép phân tách rành mạch tác động chính của từng biến độc lập, đồng thời giảm thiểu tối đa số lượng mẫu thử nghiệm mà vẫn đảm bảo độ bao phủ toàn diện của không gian trạng thái tham số.
Quy trình nghiên cứu rigorous
Quy trình thực nghiệm được chuẩn hóa qua các bước khép kín và đo lường độc lập:
- Chuẩn bị mẫu và tiền xử lý: Mẫu thép C45 kích thước tiêu chuẩn được gia công cơ khí, tẩy dầu mỡ bằng sóng siêu âm trong dung môi axeton, làm khô và tạo nhám bề mặt bằng thiết bị phun cát áp lực cao TM-R-6F (Italia) sử dụng hạt oxit nhôm/hạt thép để đạt độ nhám bề mặt $R_a = 5 – 8\text{ }\mu m$, kiểm tra bằng máy đo độ nhám tiếp xúc Mitutoyo SJ-410.
- Quy trình phun phủ plasma: Thực hiện trên hệ thống thiết bị công nghiệp tự động hóa Praxair Model 3710 (Hoa Kỳ) kết hợp súng phun plasma DC Model SG100 gắn trên robot điều khiển tọa độ nhằm đảm bảo tính đồng nhất tuyệt đối về quỹ đạo quét và tốc độ dịch chuyển súng phun ($v_s$). Lớp lót $NiCr20$ được phun trước nhằm giảm chênh lệch hệ số giãn nở nhiệt giữa thép và gốm, tiếp theo là lớp phủ gốm chính $Al_2O_3 – 40%TiO_2$ đạt chiều dày danh nghĩa $200 – 300\text{ }\mu m$, theo dõi bằng thiết bị đo chiều dày điện từ Mitutoyo DIGI-DERM 745.
- Thu thập dữ liệu và thiết bị phân tích:
- Độ cứng tế vi ($HV$): Đo theo tiêu chuẩn quốc tế ASTM E384-17 trên máy đo vi độ cứng Buehler IndentaMet 1106 với mũi đâm kim cương Vickers, tải trọng tác dụng $300\text{ }g$ ($HV_{0,3}$), thời gian giữ tải $15\text{ }s$. Mỗi mẫu đo ngẫu nhiên ít nhất 10 điểm trên mặt cắt ngang để lấy giá trị trung bình thống kê.
- Hệ số ma sát khô ($\mu$): Đo trên thiết bị đo ma sát xoay tròn TRIBO-Tester, bi trượt đối tiếp gốm $Si_3N_4$ đường kính $6\text{ }mm$, tải trọng pháp tuyến $N = 5\text{ }N$, vận tốc trượt $0,1\text{ }m/s$, quãng đường trượt tổng cộng $100\text{ }m$.
- Độ bền bám dính trượt ($\tau$): Đánh giá theo tiêu chuẩn công nghiệp Nhật Bản JIS H 8666-1980 trên máy kéo nén vạn năng thủy lực WEW-1000B với đồ gá và khuôn trượt chuyên dụng do tác giả tự thiết kế và chế tạo chính xác.
- Độ xốp ($P$): Đánh giá theo tiêu chuẩn ASTM B276-21 thông qua phương pháp phân tích ảnh kim tương kỹ thuật số trên mặt cắt ngang ở độ phóng đại $200$ lần, xử lý ngưỡng ảnh bằng phần mềm phân tích chuyên dụng.
- Hình thái cấu trúc và phân tích pha: Sử dụng kính hiển vi điện tử quét phát xạ trường SEM/EDX Jeol JSM-6490 (Nhật Bản) để quan sát sự dàn mỏng của các hạt splat, vi nứt và sự phân bố nguyên tố; phân tích thành phần pha tinh thể bằng máy nhiễu xạ tia X D8-Advance (Bruker, Đức) với bức xạ $Cu-K_\alpha$ ($\lambda = 1,5406\text{ }\text{Å}$), góc quét $2\theta$ từ $20^\circ$ đến $80^\circ$.
Data và phân tích
Phân tích số liệu áp dụng tỷ số tín hiệu trên nhiễu ($S/N$) theo hai dạng hàm mục tiêu:
- Dạng Càng lớn càng tốt (Larger-the-Better) áp dụng cho Độ cứng tế vi ($HV$) và Độ bền bám dính trượt ($\tau$):
$$(S/N){LTB} = -10 \log{10} \left( \frac{1}{n} \sum_{i=1}^n \frac{1}{y_i^2} \right)$$
- Dạng Càng nhỏ càng tốt (Smaller-the-Better) áp dụng cho Hệ số ma sát ($\mu$) và Độ xốp ($P$):
$$(S/N){STB} = -10 \log{10} \left( \frac{1}{n} \sum_{i=1}^n y_i^2 \right)$$
Phân tích phương sai (ANOVA) xác định bậc tự do ($f_j$), tổng bình phương độ lệch ($SS$), phương sai ($MS$), chuẩn định Fisher ($F_s$) và tỷ lệ đóng góp phần trăm ảnh hưởng ($P_j%$):
$$P_j = \frac{SS_j - f_j \cdot MS_{error}}{SS_{total}} \times 100%$$
Phương pháp tối ưu hóa đa mục tiêu OEC tổng hợp 4 chỉ tiêu đơn lẻ thành một chỉ số vô lượng thứ duy nhất:
$$OEC = \sum_{k=1}^4 W_k \cdot \left( \frac{Y_k - Y_{k,min}}{Y_{k,max} - Y_{k,min}} \right)$$
với $W_k$ là vector trọng số thỏa mãn $\sum W_k = 1$, phản ánh mức độ quan trọng kỹ thuật của từng chỉ tiêu đối với điều kiện làm việc chống mài mòn.
Phát hiện đột phá và implications
Những phát hiện then chốt
-
Kiểm soát vi cấu trúc và chuyển biến pha nhiệt độ cao: Kết quả phân tích nhiễu xạ tia X (XRD) và ảnh SEM mặt cắt ngang khẳng định bột phun nguyên liệu ban đầu gồm pha bền nhiệt $\alpha-Al_2O_3$ và $TiO_2$ dạng rutil sau khi đi qua dòng plasma đã chuyển hóa một phần lớn thành pha giả bền $\gamma-Al_2O_3$. Sự hình thành dung dịch rắn vô định hình giữa $Al_2O_3$ và $TiO_2$ tại các vùng biên hạt nóng chảy hoàn toàn đóng vai trò triệt tiêu các vết nứt tế vi xuyên lớp, nâng cao mật độ sắp xếp hình học của các phiến dẹt splat.
-
Phát hiện quy luật chi phối của Cường độ dòng điện phun ($I$): Dòng điện phun $I$ là thông số có tỷ lệ đóng góp chi phối cao nhất đến độ cứng tế vi ($HV$) và độ xốp ($P$) của lớp phủ, với tỷ lệ đóng góp $P_j%$ trong phân tích ANOVA vượt trên $50%$. Khi tăng dòng điện $I$, nhiệt độ cột hồ quang tăng giúp các hạt gốm $Al_2O_3 – 40%TiO_2$ nóng chảy triệt để, giảm tỷ lệ hạt chưa nóng chảy (unmelted particles), do đó độ cứng tế vi trung bình tăng mạnh từ $680\text{ }HV$ lên trên $860\text{ }HV$, đồng thời độ xốp giảm từ $5,8%$ xuống dưới $2,2%$.
-
Hiện tượng phi tuyến của Khoảng cách phun ($L$): Khác với quan niệm tuyến tính thông thường, khoảng cách phun $L$ thể hiện mối quan hệ phi tuyến rõ rệt dạng parabol đối với độ bền bám dính trượt ($\tau$) và hệ số ma sát ($\mu$). Nếu khoảng cách $L$ quá ngắn ($< 90\text{ }mm$), hiệu ứng thổi nhiệt quá mức từ ngọn lửa plasma làm quá nhiệt chất nền và gia tăng ứng suất dư nhiệt gây suy giảm liên kết bám dính. Ngược lại, nếu khoảng cách $L$ quá lớn ($> 130\text{ }mm$), các hạt gốm bị mất vận tốc và nguội sớm trong không khí trước khi va đập, dẫn đến tăng độ xốp và giảm độ bền bám dính. Khoảng cách tối ưu được xác lập chặt chẽ trong khoảng $100 – 110\text{ }mm$.
-
Ảnh hưởng tương hỗ của Lưu lượng cấp bột ($M$): Lưu lượng nạp bột $M$ chi phối trực tiếp đến mật độ hạt bên trong thể tích ngọn lửa plasma. Khi $M$ vượt quá ngưỡng giới hạn, nhiệt lượng của dòng plasma bị phân tán không đủ để làm nóng chảy đồng đều toàn bộ lượng bột, làm tăng đột biến số lượng hạt bán nóng chảy, khiến hệ số ma sát khô $\mu$ tăng từ $0,28$ lên $0,45$ và độ cứng giảm sút.
-
Xác lập nghiệm tối ưu toàn cục OEC: Bằng cách giải bài toán quy hoạch đa mục tiêu trên hàm OEC, nghiên cứu đã xác định được bộ thông số công nghệ tối ưu đồng thời cả 4 chỉ tiêu cơ tính:
- Dòng điện phun: $I = 550\text{ }A$
- Khoảng cách phun: $L = 100\text{ }mm$
- Lưu lượng cấp bột: $M = 35\text{ }g/min$
Thực nghiệm kiểm chứng độc lập tại bộ thông số này cho kết quả thực nghiệm hoàn toàn trùng khớp với giá trị lý thuyết dự báo với sai số thống kê $< 4,5%$: Độ cứng tế vi đạt $872\text{ }HV$, độ bền bám dính trượt đạt $38,6\text{ }MPa$, độ xốp đạt mức tối thiểu $2,1%$, và hệ số ma sát khô duy trì ổn định ở mức $\mu = 0,29$.
========================================================================================
BẢNG TỔNG HỢP KẾT QUẢ THỰC NGHIỆM, MÔ HÌNH DỰ BÁO VÀ KIỂM CHỨNG TỐI ƯU OEC
========================================================================================
Chỉ tiêu cơ tính Mô hình hồi quy toán học (R^2) Yếu tố chi phối lớn nhất (Pj %) Giá trị tối ưu kiểm chứng
----------------------------------------------------------------------------------------
Độ cứng tế vi (HV) HV = f(I, L, M) (R^2 = 0.94) Dòng điện phun I (58.4%) 872 HV
Độ bền bám dính (MPa) τ = f(I, L, M) (R^2 = 0.91) Khoảng cách phun L (46.2%) 38.6 MPa
Hệ số ma sát khô (μ) μ = f(I, L, M) (R^2 = 0.92) Lưu lượng cấp bột M (42.8%) 0.29
Độ xốp thể tích (%) P = f(I, L, M) (R^2 = 0.95) Dòng điện phun I (54.1%) 2.1%
----------------------------------------------------------------------------------------
Chỉ số tổng hợp OEC: Đạt cực đại tại I = 550 A, L = 100 mm, M = 35 g/min (Sai số < 4.5%)
========================================================================================
Implications đa chiều
- Về mặt lý thuyết: Mở rộng lý thuyết nhiệt động lực học biến dạng hạt gốm trong môi trường plasma không khí; làm sáng tỏ cơ chế tự bôi trơn và chống mài mòn của pha composite $\gamma-Al_2O_3 / TiO_2$ khi chịu tải trọng trượt khô.
- Về mặt phương pháp luận: Thiết lập một khung quy trình thực nghiệm - toán học mẫu mực kết hợp Taguchi, ANOVA, mô hình hồi quy giải tích đa biến và chỉ tiêu đánh giá tổng thể OEC, có khả năng chuyển giao áp dụng cho việc nghiên cứu các hệ lớp phủ gốm kim loại phức tạp khác như $WC-Co$, $Cr_3C_2-NiCr$, $ZrO_2-Y_2O_3$.
- Về mặt ứng dụng thực tiễn: Cung cấp bộ thông số công nghệ và quy trình công nghệ hoàn chỉnh cho các viện nghiên cứu, nhà máy cơ khí chế tạo để phục hồi trực tiếp các chi tiết máy làm việc trong điều kiện mài mòn xâm thực nặng như: trục bơm nước biển, xi lanh thủy lực máy công trình, cánh tuabin nhà máy thủy điện, trục cán thép, thanh dẫn hướng máy công cụ.
- Về mặt chính sách và kinh tế: Cung cấp luận cứ khoa học thực chứng phục vụ chương trình nội địa hóa sản phẩm cơ khí trọng điểm của Bộ Công Thương, hỗ trợ doanh nghiệp giảm thiểu việc nhập khẩu phụ tùng thay thế đắt đỏ từ nước ngoài, tiết kiệm $50 – 80%$ ngoại tệ cho ngân sách quốc gia.
Limitations và Future Research
Nghiên cứu duy trì tính khách quan khoa học khi thẳng thắn thừa nhận một số giới hạn thực nghiệm:
- Phạm vi khảo sát vật liệu: Nghiên cứu giới hạn tỷ lệ phối trộn cố định $Al_2O_3 – 40%TiO_2$ theo khối lượng trên nền thép C45, chưa mở rộng khảo sát dải biến thiên liên tục của hàm lượng $TiO_2$ từ $0%$ đến $50%$ trong cùng một điều kiện thực nghiệm.
- Môi trường phun phủ: Thực nghiệm thực hiện bằng công nghệ APS trong môi trường không khí tự nhiên, dù có chi phí thấp và tính khả thi công nghiệp cao nhưng vẫn tồn tại lượng nhỏ tạp chất oxit và lỗ xốp tế vi không thể triệt tiêu hoàn toàn như trong môi trường buồng chân không (VPS) hay khí quyển có kiểm soát (CAPS).
- Số lượng thông số đầu vào: Mô hình quy hoạch tập trung vào 3 thông số chính ($I, L, M$), trong khi các thông số phụ trợ như tỷ lệ lưu lượng khí sơ cấp/thứ cấp ($Ar/H_2$), góc nghiêng súng phun và nhiệt độ sấy nóng sơ bộ chất nền được giữ cố định theo khuyến nghị của nhà sản xuất súng SG100.
- Phạm vi thử nghiệm mài mòn: Các phép thử ma sát học được tiến hành ở nhiệt độ phòng ($25^\circ\text{C}$) dưới tải trọng trượt khô tiếp xúc điểm, chưa mô phỏng hết môi trường mài mòn - xói mòn kết hợp ở nhiệt độ cao ($> 600^\circ\text{C}$) hoặc trong môi trường dung dịch axit/kiềm đậm đặc.
Chương trình nghiên cứu tiếp theo (Future Research Agenda):
- Phát triển mô hình mạng nơ-ron nhân tạo (Artificial Neural Network - ANN) kết hợp thuật toán tối ưu hóa bầy đàn (PSO) hoặc thuật toán di truyền (GA) để tối ưu hóa không gian thực nghiệm phi tuyến bậc cao với số lượng thông số đầu vào mở rộng ($\ge 6$ thông số).
- Nghiên cứu xử lý sau phun (post-treatment) bằng kỹ thuật tẩm bịt dung dịch muối nhôm phốt phát ($AlPO_4$) kết hợp ủ nhiệt biến tính bằng tia laser để triệt tiêu hoàn toàn độ xốp còn lại, đưa độ xốp về mức $< 0,5%$.
- Đánh giá độ bền mỏi nhiệt và độ bền ăn mòn điện hóa (EIS) dài hạn của lớp phủ $Al_2O_3 – 40%TiO_2$ trong môi trường nước biển nhân tạo mô phỏng giàn khoan dầu khí ngoài khơi.
- Khảo sát ứng dụng bột gốm cấu trúc nano tổng hợp để so sánh hiệu năng cơ tính và ma sát học với bột micro truyền thống.
Tác động và ảnh hưởng
- Tác động học thuật: Công trình đóng góp hệ thống dữ liệu thực nghiệm chuẩn mực về hệ lớp phủ gốm plasma $Al_2O_3 – 40%TiO_2$ trên nền thép kết cấu, dự kiến thu hút trích dẫn khoa học cao trong các tạp chí chuyên ngành uy tín thuộc danh mục ISI/Scopus như Journal of Thermal Spray Technology (ASM International), Surface and Coatings Technology (Elsevier), và Wear.
- Chuyển đổi ngành công nghiệp cơ khí: Thúc đẩy các doanh nghiệp cơ khí chế tạo, sửa chữa tại Việt Nam (như các đơn vị thuộc Tổng công ty Máy và Thiết bị công nghiệp MIE, Tập đoàn Than - Khoáng sản Việt Nam TKV, Viện Nghiên cứu Cơ khí NARIME) chuyển dịch từ phương pháp hàn đắp hoặc mạ crôm cứng truyền thống gây ô nhiễm môi trường sang công nghệ phun phủ plasma xanh, hiệu suất cao và có độ bền vượt trội.
- Tác động xã hội và môi trường: Giảm thiểu phát thải hóa chất độc hại (như axit cromic trong công nghệ mạ crôm cứng), gia tăng tuổi thọ làm việc của máy móc thiết bị công nghiệp lên gấp 2 đến 3 lần, giảm thiểu lượng phế thải kim loại thải ra môi trường và nâng cao an toàn vận hành trong các ngành công nghiệp nặng.
Đối tượng hưởng lợi
- Nghiên cứu sinh và học viên cao học ngành Kỹ thuật cơ khí / Khoa học vật liệu: Tiếp cận một mô hình nghiên cứu mẫu mực, kết hợp nhuần nhuyễn giữa cơ sở lý thuyết vật lý plasma, phương pháp quy hoạch thực nghiệm Taguchi và công cụ toán học tối ưu hóa đa mục tiêu OEC.
- Các nhà khoa học và chuyên gia bề mặt: Kế thừa hệ thống phương trình hồi quy và các dữ liệu phân tích pha XRD, ảnh SEM phân giải cao để phát triển các đề tài nghiên cứu chuyên sâu về composite gốm tiên tiến.
- Kỹ sư R&D và cán bộ kỹ thuật tại các nhà máy cơ khí: Áp dụng trực tiếp bộ thông số công nghệ tối ưu ($I = 550\text{ }A, L = 100\text{ }mm, M = 35\text{ }g/min$) vào quy trình sản xuất thực tế trên thiết bị phun plasma Praxair 3710 / SG100 mà không cần tốn kém chi phí thử nghiệm thăm dò.
- Các nhà hoạch định chính sách công nghiệp: Có cơ sở khoa học vững chắc để xây dựng các tiêu chuẩn kỹ thuật quốc gia về phun phủ nhiệt và ban hành các chính sách khuyến khích nội địa hóa công nghệ xử lý bề mặt cao định hướng công nghiệp 4.0.
Câu hỏi chuyên sâu
1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?
Đóng góp lý thuyết độc đáo nhất là việc mở rộng Lý thuyết hình thành lớp màng va đập của Max Ulrich Schoop và Pospisil - Sehyl từ vật liệu kim loại dẻo sang hệ vật liệu gốm composite oxit đa pha ($Al_2O_3 – 40%TiO_2$). Luận án chứng minh rằng cơ chế liên kết bám dính không chỉ dựa vào sự biến dạng dẻo thuần túy cơ học mà dựa trên sự nóng chảy chọn lọc của pha $TiO_2$ hoạt động như một chất trợ dung nhiệt độ cao, bao bọc các hạt $\alpha-Al_2O_3$, đồng thời sự làm nguội siêu nhanh ($10^6 – 10^7\text{ }^\circ\text{C/s}$) kích thích sự chuyển biến pha thành $\gamma-Al_2O_3$, tạo nên mạng lưới liên kết cài khớp vi mô có độ dẻo dai cao và giảm ứng suất cắt tại giao diện liên kết với lớp lót $NiCr20$.
2. Tính đổi mới trong phương pháp luận của luận án khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đây?
So với nghiên cứu của Islak et al. (2014) và Yusoff et al. (2012) (chỉ sử dụng phương pháp khảo sát đơn biến truyền thống từng yếu tố riêng lẻ hoặc phân tích thống kê mô tả), luận án đã thực hiện một bước tiến phương pháp luận vượt bậc: Thiết lập mô hình quy hoạch thực nghiệm vững bền Taguchi $L_9(3^3)$ kết hợp phân tích phương sai ANOVA đa biến, xây dựng hệ phương trình hồi quy giải tích phi tuyến bậc hai xác thực cho 4 chỉ tiêu chất lượng đồng thời, và lần đầu tiên tích hợp Chỉ số đánh giá tổng thể (Overall Evaluation Criteria - OEC) để giải quyết trọn vẹn bài toán tối ưu hóa đa mục tiêu nghịch biến trong công nghệ phun plasma.
3. Phát hiện bất ngờ nhất về mặt số liệu thực nghiệm là gì và giải thích cơ chế vật lý?
Phát hiện bất ngờ nhất là sự phi tuyến của khoảng cách phun ($L$) đối với hệ số ma sát khô ($\mu$) và độ bền bám dính ($\tau$). Khi tăng dòng điện $I$ lên $550\text{ }A$ nhưng nếu giữ khoảng cách phun ngắn ($L = 80\text{ }mm$), độ bám dính không tăng mà lại suy giảm nghiêm trọng ($< 25\text{ }MPa$) và độ xốp gia tăng cục bộ. Cơ chế vật lý được giải thích là do mật độ năng lượng nhiệt plasma quá lớn ở cự ly gần đã gây ra hiện tượng quá nhiệt bề mặt nền, làm bốc hơi một phần pha $TiO_2$ và tạo ra gradien nhiệt độ cực lớn giữa lớp phủ ngoài và nền kim loại C45, dẫn đến tích tụ ứng suất dư nhiệt gây nứt vi mô tức thời trong quá trình đông đặc.
4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) hoàn chỉnh không?
Quy trình thực nghiệm được mô tả chi tiết với độ chính xác cao: từ thông số chuẩn bị mẫu (thép C45, độ nhám hạt cát $R_a = 5 – 8\text{ }\mu m$), cấu hình súng phun Praxair SG100 (đường kính vòi phun, góc cấp bột $90^\circ$, lưu lượng khí chính Ar và khí phụ $H_2$), tốc độ quét của robot, quy chuẩn lớp lót $NiCr20$ ($50 – 80\text{ }\mu m$) đến kích thước khuôn trượt kiểm tra độ bám dính theo tiêu chuẩn JIS H 8666-1980 và phương pháp chuẩn bị mẫu kim tương hiển vi điện tử quét. Bất kỳ phòng thí nghiệm vật liệu bề mặt nào có trang bị hệ thống phun plasma công nghiệp đều có thể tái lập chính xác 100% các kết quả thực nghiệm này.
5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?
Lộ trình 10 năm được xác lập qua 3 giai đoạn chiến lược:
- Giai đoạn 1 (Năm 1 - 3): Tích hợp công nghệ học máy (Machine Learning) và thị giác máy tính vào hệ thống phun plasma tự động để giám sát trực tiếp dòng hạt bay (in-flight particle diagnostics: nhiệt độ, vận tốc hạt) và điều khiển vòng lặp kín thời gian thực.
- Giai đoạn 2 (Năm 4 - 6): Hoàn thiện công nghệ biến tính bề mặt sau phun bằng xử lý nhiệt laser (Laser Remelting) và tẩm compozit hạt nano nhằm chế tạo các lớp phủ siêu đặc sít có độ xốp tiệm cận $0%$, chịu ăn mòn axit siêu việt.
- Giai đoạn 3 (Năm 7 - 10): Mở rộng sang các hệ vật liệu gốm đa nguyên tử cấu trúc entropy cao (High-Entropy Oxides - HEOs) và gốm nhiệt động học phục vụ ngành hàng không vũ trụ và động cơ tuabin khí thế hệ mới.
Kết luận
- Nghiên cứu đã làm chủ toàn diện quy trình công nghệ phun phủ plasma trong không khí (APS) để chế tạo lớp phủ gốm composite $Al_2O_3 – 40%TiO_2$ chất lượng cao trên nền thép kết cấu C45 thông qua lớp lót trung gian $NiCr20$.
- Xây dựng thành công hệ phương trình hồi quy đa biến bậc hai mô tả chính xác tương quan toán học giữa 3 thông số công nghệ đầu vào ($I, L, M$) và 4 chỉ tiêu cơ tính phức hợp ($HV, \mu, \tau, P$) với hệ số tương quan $R^2 > 0,90$.
- Lượng hóa chính xác mức độ ảnh hưởng của từng thông số thông qua ANOVA, chứng minh Cường độ dòng điện phun ($I$) chi phối trên $50%$ độ cứng và độ xốp, trong khi Khoảng cách phun ($L$) và Lưu lượng cấp bột ($M$) quyết định độ bền bám dính và hệ số ma sát.
- Ứng dụng thành công chỉ số đánh giá tổng thể (OEC), tìm ra bộ thông số công nghệ tối ưu toàn cục ($I = 550\text{ }A, L = 100\text{ }mm, M = 35\text{ }g/min$), tạo ra lớp phủ có độ cứng $872\text{ }HV$, độ bám dính $38,6\text{ }MPa$, độ xốp $2,1%$ và hệ số ma sát $\mu = 0,29$.
- Chế tạo thành công hệ thống đồ gá phun định vị và khuôn kiểm tra độ bền bám dính trượt chuyên dụng theo tiêu chuẩn JIS H 8666-1980 cho mẫu tấm phẳng, khắc phục hoàn toàn hạn chế đo lường cơ học trước đây.
- Kết quả nghiên cứu là nền tảng khoa học và công nghệ vững chắc, mở ra hướng ứng dụng rộng rãi trong việc chế tạo mới và phục hồi các chi tiết máy chịu mài mòn, ăn mòn khốc liệt, mang lại hiệu quả kinh tế - kỹ thuật to lớn và đóng góp tích cực vào sự phát triển của ngành kỹ thuật cơ khí Việt Nam.