Tổng quan nghiên cứu

Phương pháp nhiễu xạ tia X là một trong những kỹ thuật kiểm tra không phá hủy tiên tiến hàng đầu thế giới, sở hữu khả năng phân tích cấu trúc vật chất ở kích thước siêu vi đạt mức 10^-7 mm, vượt trội hoàn toàn so với giới hạn 10^-3 mm của các phương pháp chụp X-quang thông thường. Tại Việt Nam, việc ứng dụng tia X trong phân tích cơ tính và tổ chức tế vi của vật liệu cơ khí như xác định ứng suất dư, tỷ lệ pha, kích thước tinh thể và chiều dày màng phủ vẫn gặp nhiều rào cản do các công đoạn tính toán còn phụ thuộc nặng nề vào thao tác thủ công hoặc sử dụng các công cụ bảng tính rời rạc. Để giải quyết triệt để bài toán này, luận văn thạc sĩ đã tập trung nghiên cứu và phát triển thành công phần mềm chuyên dụng X-Pro 1.0 trên nền tảng ngôn ngữ lập trình C#.

Công trình được thực hiện trong thời gian 06 tháng tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh dưới sự hướng dẫn khoa học của Phó Giáo sư Lê Chí Cương, có sự phối hợp dữ liệu thực nghiệm từ Trung tâm Hạt nhân Thành phố Hồ Chí Minh trực thuộc Viện Năng lượng Nguyên tử Việt Nam. Mục tiêu trọng tâm của đề tài là xây dựng một công cụ thuần Việt có khả năng tự động hóa 100% quy trình xử lý dữ liệu nhiễu xạ, đồng thời phát triển 3 tính năng nâng cao gồm: định lượng tỷ lệ pha của vật liệu 3 pha, đo lường kích thước tinh thể nano và tính toán chiều dày lớp mạ kim loại. Kết quả nghiên cứu mang lại ý nghĩa kinh tế và kỹ thuật to lớn, giúp cắt giảm hơn 80% thời gian xử lý dữ liệu thực nghiệm, đồng thời kiểm soát sai số tính toán ở mức dưới 5%, đặt nền móng quan trọng cho việc tự động hóa kiểm định vật liệu trong nền công nghiệp hiện đại.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nền tảng lý thuyết của nghiên cứu được xây dựng trên hệ thống các định luật vật lý và mô hình toán học giải tích tinh thể học chuẩn mực. Trọng tâm đầu tiên là Định luật Bragg, được thiết lập năm 1933, mô tả mối quan hệ giữa bước sóng tia X và khoảng cách giữa các mặt phẳng nguyên tử thông qua phương trình n nhân lambda bằng 2 nhân d nhân sin theta, làm cơ sở xác định họ mặt phẳng nhiễu xạ hkl. Trọng tâm thứ hai là phương trình Scherrer dùng để định lượng kích thước tinh thể trung bình thông qua độ rộng bán cực đại FWHM của vạch phổ tán xạ. Bên cạnh đó, tác giả áp dụng lý thuyết năng lượng nhiễu xạ toàn phần để xác định tỷ lệ thành phần pha trong vật liệu đa pha, trong đó năng lượng mỗi pha tỷ lệ thuận với tích phân diện tích cường độ trên các mặt phản xạ.

Nhằm chuẩn hóa tín hiệu đầu vào, nghiên cứu ứng dụng mô hình hiệu chỉnh đa tầng bao gồm: thuật toán lọc nhiễu cửa sổ trượt trung bình động với số điểm đối xứng hai phía nL bằng nR bằng 9 điểm, giải thuật hiệu chỉnh đường nền tuyến tính nhằm triệt tiêu tín hiệu nền do tạp chất, và mô hình hiệu chỉnh hệ số hấp thụ phân cực Lorentz-Polarization-Absorption (LPA) theo các phương pháp đo góc nghiêng đồng trục và góc nghiêng bên. Ngoài ra, tọa độ đỉnh phổ được định vị chính xác thông qua 4 phương pháp toán học: phương pháp xấp xỉ Parabola bằng bình phương bé nhất, mô hình hàm Gaussian logarit, phương pháp trọng tâm hình học và phương pháp bề rộng trung bình.

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

Nghiên cứu sử dụng nguồn dữ liệu thực nghiệm thu thập từ các hệ thống máy đo nhiễu xạ tia X chuyên dụng, bao gồm giản đồ nhiễu xạ của mẫu chuẩn Ceri Oxit (CeO2), vật liệu rây phân tử Zeolite ZSM-5 có kích thước hạt tinh thể chuẩn 3 micromet từ Viện Ruder Boskovic, hệ vật liệu đa pha công nghiệp gồm Vonfram Cacbit (WC), Titan Cacbit (TiC), Coban (Co), và chi tiết thép SS400 mạ Niken. Cỡ mẫu nghiên cứu gồm 4 nhóm vật liệu đại diện cho các họ mạng tinh thể lập phương tâm mặt, vật liệu composit đa pha và màng mỏng bề mặt.

Phương pháp chọn mẫu có chủ đích được lựa chọn nhằm mục tiêu kiểm chứng độ chính xác của các thuật toán trong từng điều kiện hình học và cấu trúc tinh thể riêng biệt. Phương pháp phân tích cốt lõi là số hóa thuật toán trên môi trường Visual Studio bằng ngôn ngữ C#, tích hợp giao diện tương tác đồ họa và tự động hóa tính toán giải tích. Lý do lựa chọn phương pháp này xuất phát từ nhu cầu xử lý hàng ngàn điểm dữ liệu đo đạc liên tục chỉ trong vài giây, loại bỏ hoàn toàn các sai số chủ quan do người phân tích khi dựng tiếp tuyến hoặc tích phân thủ công, đồng thời nâng cao độ tin cậy của kết quả phân tích trong khoảng tin cậy 95% với hệ số 1,96 độ lệch chuẩn. Toàn bộ tiến trình nghiên cứu từ xây dựng thuật toán, viết mã nguồn, đo đạc đối chứng đến hoàn thiện phần mềm được tiến hành đồng bộ trong lộ trình 06 tháng từ giữa năm 2015 đến tháng 01 năm 2016.

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

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

Quá trình kiểm nghiệm phần mềm X-Pro 1.0 trên các tập dữ liệu thực tế đã mang lại 4 phát hiện kỹ thuật cốt lõi:

Thứ nhất, thuật toán tự động nhận diện đỉnh phổ dựa trên đạo hàm tiếp tuyến với góc ngưỡng alpha lớn hơn 20 độ kết hợp bộ lọc làm mịn 9 điểm đã loại bỏ hoàn toàn các xung nhiễu nền giả định, nhận diện chính xác 100% các mặt phản xạ Miller hkl của vật liệu mẫu CeO2 tại các góc 2 theta đặc trưng gồm 28,55 độ, 33,08 độ, 47,48 độ và 56,34 độ, độ mịn của tín hiệu tăng hơn 85% so với dữ liệu thô ban đầu.

Thứ hai, phân hệ định lượng vật liệu 3 pha (WC - TiC - Co) khi áp dụng vùng đáy nhiễu xạ bằng 5 lần bề rộng trung bình FWHM đã tính toán thành phần thể tích từng pha với sai số tuyệt đối dưới 2,8% so với phương pháp phân tích hóa học ướt truyền thống, rút ngắn thời gian phân tích từ 48 giờ thử nghiệm hóa học xuống còn dưới 15 giây xử lý trên máy vi tính.

Thứ ba, công thức Scherrer cải tiến sử dụng tỷ số tương đối với mẫu chuẩn Zeolite ZSM-5 kích thước 3 micromet đã định lượng chính xác kích thước hạt tinh thể của các mẫu Zeolite A2 và A1-2, giúp loại bỏ hoàn toàn sự phụ thuộc vào thông số đóng góp sai số cơ học của thiết bị đo.

Thứ tư, phương pháp phân tích chiều dày lớp mạ Niken trên nền thép hợp kim SS400 bằng tia X đạt kết quả tương thích cao, sai lệch chỉ khoảng 4,2% so với phương pháp kiểm tra phá hủy bằng kính hiển vi quang học kim tương ở độ phóng đại 1000 đến 2000 lần và phương pháp đo cảm ứng từ trường trong phạm vi 20% đến 80% thang đo cực đại.

Thảo luận kết quả

Độ chính xác vượt trội của phần mềm X-Pro 1.0 bắt nguồn từ việc tích hợp chặt chẽ khâu hiệu chỉnh nền và hệ số hấp thụ phân cực LPA ngay trong quá trình xử lý toán học. Khi so sánh với các phần mềm quốc tế như Maud, MDI Jade, Axio Vision hay OriginPro, phần mềm X-Pro 1.0 thể hiện ưu thế vượt bậc về tính tiện dụng nhờ giao diện thuần Việt, quy trình phân tích khép kín từ nhập tập tin dữ liệu thô, tự động phân tách đỉnh đến kết xuất báo cáo, thay vì yêu cầu người dùng phải chuyển đổi định dạng và thao tác qua nhiều công cụ trung gian phức tạp.

Về mặt trình diễn thông tin, toàn bộ dữ liệu phân tích có thể được mô hình hóa trực quan thông qua các giản đồ phổ tán xạ trong dải góc 2 theta từ 10 đến 90 độ, tích hợp các bảng tra cứu tổng bình phương chỉ số Miller h^2 + k^2 + l^2 và đường hồi quy tuyến tính sin bình phương Psi trong phân tích ứng suất. Kết quả kiểm định thống kê cho thấy các giá trị tính toán đều nằm trong dải dung sai cho phép với khoảng tin cậy 95%, khẳng định tính khả thi và độ tin cậy tuyệt đối của phần mềm trong việc thay thế các công cụ thương mại đắt tiền.

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

Để phát huy tối đa giá trị thực tiễn của phần mềm X-Pro 1.0 trong nghiên cứu khoa học và sản xuất công nghiệp, 4 giải pháp cụ thể được đề xuất triển khai:

Thứ nhất, mở rộng thư viện cơ sở dữ liệu tinh thể học tự động: Nhóm phát triển phần mềm cần chủ trì bổ sung kho dữ liệu chuẩn từ 50 mẫu cơ bản lên hơn 5.000 cấu trúc tinh thể của các hợp kim và khoáng vật phổ biến, tích hợp định dạng cơ sở dữ liệu quốc tế ICDD trong thời gian 12 tháng tới, nhằm tăng tỷ lệ tự động nhận diện pha lên trên 98%.

Thứ hai, tối ưu hóa thuật toán tách đỉnh phổ chập: Các kỹ sư thuật toán cần nghiên cứu tích hợp mô hình hàm Pseudo-Voigt hoặc hàm Pearson VII vào mô-đun phân tích đỉnh trong vòng 06 tháng, nhằm nâng cao độ chính xác phân tách các dải phổ chồng lấn phức tạp thêm 15% đối với các vật liệu kích thước nanomet.

Thứ ba, xúc tiến thương mại hóa và chuyển giao công nghệ: Ban Giám hiệu nhà trường phối hợp cùng Trung tâm Hạt nhân Thành phố Hồ Chí Minh triển khai ứng dụng thử nghiệm phần mềm tại ít nhất 10 phòng thí nghiệm vật liệu, viện nghiên cứu và doanh nghiệp cơ khí chế tạo trong lộ trình 18 tháng, hướng tới mục tiêu nội địa hóa công cụ phần mềm kỹ thuật.

Thứ tư, chuẩn hóa quy trình đào tạo và kiểm định NDT: Bộ môn Khoa học Vật liệu cần biên soạn bộ tài liệu chuẩn hóa quy trình phân tích nhiễu xạ tia X tích hợp phần mềm X-Pro 1.0 trong vòng 09 tháng, giúp cắt giảm thời gian đào tạo kỹ thuật viên phòng thí nghiệm từ 120 giờ xuống dưới 35 giờ thao tác chuẩn mực.

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

Luận văn thạc sĩ là tài liệu khoa học chuyên sâu, cung cấp giá trị ứng dụng thiết thực cho 4 nhóm đối tượng trọng điểm:

Thứ nhất, các kỹ sư cơ khí và chuyên gia luyện kim tại các nhà máy chế tạo: Ứng dụng quy trình đo chiều dày lớp mạ Niken trên thép SS400 và phương pháp định lượng pha hợp kim cứng WC-TiC-Co để giám sát chất lượng nhiệt luyện, cơ tính và độ bền bỉ của lớp phủ bảo vệ chi tiết máy.

Thứ hai, các nhà khoa học và nghiên cứu sinh ngành vật liệu nano và hóa học: Khai thác phương trình Scherrer cải tiến kết hợp mẫu chuẩn Zeolite 3 micromet để đo lường chính xác kích thước hạt nano của vật liệu xúc tác, vật liệu hấp phụ mà không cần đầu tư hệ thống đo quang sai phức tạp.

Thứ ba, các kỹ sư phần mềm khoa học và lập trình viên công nghiệp: Tham khảo kiến trúc mã nguồn C#, thuật toán xử lý tín hiệu số, giải thuật lọc nhiễu cửa sổ trượt và các phương pháp xấp xỉ toán học Parabola, Gaussian phục vụ phát triển các phần mềm kỹ thuật chuyên sâu.

Thứ tư, giảng viên, học viên cao học và sinh viên các khối ngành kỹ thuật cơ khí, vật lý kỹ thuật: Sử dụng luận văn làm tài liệu học thuật tham khảo mẫu mực về kỹ thuật kiểm tra không phá hủy (NDT), phương pháp mô hình hóa toán học và quy trình kiểm chứng thực nghiệm chuẩn khoa học.

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

  1. Phần mềm X-Pro 1.0 có những ưu thế nổi bật nào so với các phần mềm quốc tế như MDI Jade hay Maud? Phần mềm X-Pro 1.0 được thiết kế với giao diện tiếng Việt thân thiện, tích hợp đồng bộ 4 chức năng cốt lõi gồm phân tích đỉnh phổ, định lượng 3 pha, đo kích thước tinh thể và xác định chiều dày mạ trong cùng một môi trường. Điều này giúp giảm hơn 75% thời gian thao tác so với việc kết hợp rời rạc các phần mềm ngoại nhập.

  2. Thuật toán làm mịn dữ liệu trong luận văn hoạt động như thế nào để loại bỏ xung nhiễu? Luận văn áp dụng giải thuật lọc trung bình động với cửa sổ trượt đối xứng gồm 9 điểm bên trái và 9 điểm bên phải điểm đo. Phương pháp này triệt tiêu hoàn toàn các gai nhiễu tín hiệu ngẫu nhiên nhưng vẫn bảo toàn nguyên vẹn diện tích tích phân và vị trí tâm đỉnh phổ thực tế của mẫu đo.

  3. Vì sao nghiên cứu lại chọn ngưỡng góc tiếp tuyến lớn hơn 20 độ để xác định chân mặt nhiễu xạ? Thực nghiệm trên nhiều mẫu phổ cho thấy khi góc nghiêng tiếp tuyến đạt từ 15 đến 30 độ, đặc biệt là ngưỡng mặc định 20 độ, điểm tính toán vừa chạm vào chân của đỉnh phổ nhiễu xạ hkl. Giá trị này giúp thuật toán nhận diện chính xác chân phổ mà không bị nhầm lẫn với các dao động ngẫu nhiên của đường nền.

  4. Phương pháp xác định kích thước tinh thể Zeolite trong phần mềm có điểm gì đột phá? Điểm mới của phương pháp là việc thiết lập tỷ số tương đối giữa mẫu nghiên cứu với mẫu chuẩn Zeolite ZSM-5 có kích thước 3 micromet. Kỹ thuật này triệt tiêu hoàn toàn đại số sai số đóng góp từ thiết bị đo, cho phép xác định trực tiếp kích thước hạt nano mà không cần đo phổ chuẩn thiết bị phức tạp.

  5. Độ tin cậy của phép đo chiều dày lớp mạ bằng tia X được kiểm chứng như thế nào? Chiều dày lớp mạ Niken đo bằng phần mềm được đối chiếu đồng thời với phương pháp kim tương ở độ phóng đại 1000 đến 2000 lần và phương pháp cảm ứng từ trường. Kết quả cho thấy độ chính xác đạt khoảng tin cậy 95% trong giới hạn dung sai 1,96 độ lệch chuẩn, khẳng định tính khả thi vượt trội trong kiểm định công nghiệp.

Kết luận

  • Phát triển thành công phần mềm chuyên dụng X-Pro 1.0 bằng ngôn ngữ lập trình C# với giao diện tiếng Việt trực quan, làm chủ hoàn toàn công nghệ phân tích vật liệu bằng nhiễu xạ tia X.
  • Hoàn thiện hệ thống thuật toán tự động hóa xử lý tín hiệu nhiễu xạ từ khâu làm mịn dữ liệu với cửa sổ 9 điểm, hiệu chỉnh nền, hiệu chỉnh LPA đến nhận diện đỉnh phổ hkl chính xác 100%.
  • Bổ sung 3 tính năng phân tích nâng cao gồm xác định tỷ lệ vật liệu 3 pha WC-TiC-Co với sai số dưới 3%, đo kích thước tinh thể Zeolite bằng mẫu chuẩn 3 micromet và tính chiều dày mạ Niken đạt độ tin cậy 95%.
  • Đóng góp giải pháp công nghệ quan trọng giúp nâng cao năng suất xử lý dữ liệu lên hơn 80%, tiết kiệm chi phí đầu tư phần mềm ngoại nhập cho các phòng thí nghiệm trong nước.
  • Kế hoạch tiếp theo trong 12 tháng tới sẽ tập trung nâng cấp cơ sở dữ liệu phổ chuẩn và tích hợp giải thuật tách đỉnh phổ chập Pseudo-Voigt. Quý độc giả, các nhà nghiên cứu và doanh nghiệp hãy kết nối ngay với nhóm tác giả để tiếp cận bản cài đặt thử nghiệm và thúc đẩy hợp tác chuyển giao công nghệ.