Tổng quan nghiên cứu

Ứng suất dư hình thành trong quá trình gia công cơ, xử lý nhiệt và luyện kim là nguyên nhân trực tiếp dẫn đến biến dạng hình học hoặc phá hủy cấu trúc chi tiết máy. Trong thực tế chế tạo cơ khí chính xác, việc kiểm soát ứng suất dư đóng vai trò then chốt nhằm nâng cao tuổi thọ làm việc của sản phẩm. Phương pháp nhiễu xạ tia X (XRD) từ lâu đã khẳng định ưu thế vượt trội nhờ khả năng kiểm tra không phá hủy, cho kết quả định lượng chính xác và dễ dàng tự động hóa trên dây chuyền công nghiệp. Vị trí đỉnh nhiễu xạ trên phổ thực nghiệm là đại lượng cơ bản để xác định mức độ biến dạng mạng tinh thể, chịu ảnh hưởng đồng thời bởi tổ hợp hệ số LPA bao gồm hiệu ứng Lorentz, hiệu ứng phân cực và hàm hấp thụ bức xạ. Trong đó, hàm hấp thụ có tác động mạnh nhất làm dịch chuyển vị trí đỉnh phổ và làm sai lệch giá trị ứng suất tính toán.

Phần lớn các thiết bị nhiễu xạ tia X thương mại hiện nay vẫn mặc định sử dụng công thức hàm hấp thu trên bề mặt phẳng cho mọi dạng hình học mẫu đo. Khi kích thước chiếu xạ tia X dao động từ 1 milimét vuông đến 100 milimét vuông tác động lên các mẫu có bán kính cong nhỏ, giả thiết mặt phẳng không còn phù hợp, dẫn đến sai số đo đạc vượt ngưỡng 15% đến 20%. Đề tài tập trung giải quyết bài toán cấp thiết: xác định hàm hấp thu tổng quát cho bề mặt cong dạng Ellipsoid sử dụng phương pháp đo psi (side-inclination), cố định góc eta và eta không. Công trình được thực hiện tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh, hoàn thành vào tháng 10 năm 2014, khảo sát trên nền vật liệu sắt với bức xạ đặc trưng Cr-Kalpha có hệ số hấp thụ tuyến tính mu bằng 873. Kết quả mở ra công thức tổng quát có thể suy biến linh hoạt cho mọi bề mặt hình học từ mặt phẳng, mặt trụ đến mặt cầu, nâng cao độ chuẩn xác cho kỹ thuật đo ứng suất dư quang tuyến X.

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 của phương pháp dựa trên lý thuyết tán xạ bức xạ điện từ và hiện tượng giao thoa sóng tia X trong môi trường tinh thể tuân theo định luật Bragg: bước sóng bức xạ bằng hai lần khoảng cách giữa hai mặt phẳng nguyên tử nhân với sin của góc nhiễu xạ theta. Khi chùm tia X chiếu vào vật rắn, các điện tử dao động và phát ra sóng tán xạ thứ cấp. Nếu hiệu quang trình giữa các sóng tán xạ là một số nguyên lần bước sóng, hiện tượng giao thoa tăng cường sẽ xuất hiện, tạo nên các đỉnh nhiễu xạ sắc nét tương ứng với các chỉ số Miller h, k, l của ô cơ sở.

Cường độ đỉnh phổ thu nhận trên thực tế bị suy giảm bởi tổ hợp ba yếu tố hình học và vật lý LPA. Hệ số Lorentz L biểu thị mối quan hệ phụ thuộc vào góc nhiễu xạ theta theo biểu thức nghịch đảo của tích số sin bình phương theta và cosin theta. Hệ số phân cực P mô tả sự phân cực của chùm tia tán xạ, phụ thuộc vào một cộng với bình phương cosin của hai lần góc theta chia cho hai. Yếu tố quan trọng nhất là hàm hấp thụ A, được xây dựng từ quy luật suy giảm cường độ theo chiều dài quãng đường truyền của tia tới và tia phản xạ bên trong lòng vật liệu. Bên cạnh đó, mô hình hình học Ellipsoid bậc hai được thiết lập trong hệ tọa độ ba chiều với hai bán kính xích đạo bằng nhau và một bán kính trục cực, cho phép mô phỏng toàn diện sự tương tác thể tích giữa chùm tia X và vật mẫu cong dưới phương pháp đo psi.

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

Nghiên cứu sử dụng phương pháp giải tích hình học kết hợp mô phỏng số trên môi trường phần mềm Matlab 2010a và Autocad 2007. Khung tính toán phân chia khối Ellipsoid thành các phân tố thể tích vô cùng bé r nhân dr nhân d omega nhân dz cắt bởi các mặt phẳng song song với mặt phẳng đáy. Quãng đường truyền tia tới AB và tia phản xạ BC tới phân tố nhiễu xạ sâu bên trong mẫu được xác định chính xác thông qua các hệ thức lượng trong tam giác không gian, khắc phục triệt để giả thiết dừng chùm tia ở độ sâu giới hạn mang 96% năng lượng như các nghiên cứu trước đây.

Về cỡ mẫu và dữ liệu khảo sát, nghiên cứu thiết lập mô hình tính toán số trên 2 bộ kích thước hình học chuẩn bao gồm mẫu có bán kính Ra bằng 4 centimet, Rb bằng 7 centimet và mẫu có Ra bằng 6 centimet, Rb bằng 9 centimet. Dữ liệu quét phổ thực hiện trên dải góc nhiễu xạ hai theta từ 80 độ đến 180 độ với 3 mức góc nghiêng psi lần lượt là 10 độ, 20 độ và 30 độ. Toàn bộ quá trình chọn mẫu tham số được tối ưu hóa nhằm bao quát đầy đủ các biên độ cong từ cong thoải đến cong gắt trong chi tiết máy công nghiệp. Thuật toán tích phân 3 lớp được giải bằng phương pháp số trên phần mềm chuyên dụng trong khoảng thời gian từ cuối năm 2013 đến tháng 10 năm 2014, đảm bảo độ hội tụ cao và triệt tiêu sai số làm tròn số học.

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

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

Thứ nhất, nghiên cứu đã thiết lập thành công phương trình vi phân và công thức tích phân 3 lớp tổng quát xác định hàm hấp thu A cho bề mặt Ellipsoid dưới phương pháp đo psi với diện tích chiếu xạ không bị giới hạn. Biểu thức này tích hợp đầy đủ các tham số quang học góc tới gamma, góc phản xạ beta, tọa độ hình học phân tố và hàm bán kính cong biến thiên dọc theo trục cực.

Thứ hai, kết quả mô phỏng số trên nền vật liệu sắt với hệ số hấp thụ tuyến tính mu bằng 873 cho thấy quy luật biến thiên rõ nét: khi góc nhiễu xạ hai theta tăng dần từ 80 độ đến 180 độ, hệ số hấp thụ A tăng đều theo đường phi tuyến. Khi tăng góc nghiêng mẫu psi từ 10 độ lên 30 độ, giá trị hấp thu tổng thể bị suy giảm từ 25% đến 40% trên toàn dải quét.

Thứ ba, sự thay đổi kích thước biên dạng hình học tác động trực tiếp đến hàm hấp thu. Khi mở rộng bán kính cong từ Ra bằng 4 centimet, Rb bằng 7 centimet lên Ra bằng 6 centimet, Rb bằng 9 centimet, đường cong hàm hấp thụ có xu hướng tịnh tiến và tiệm cận dần về đường đặc tính của bề mặt phẳng. Độ lệch giá trị hấp thu giữa bề mặt cong Ellipsoid và bề mặt phẳng đạt mức chênh lệch cao nhất khoảng 18.5% tại góc nghiêng 30 độ.

Thứ tư, mô hình toán học đạt độ tương thích hoàn hảo khi thực hiện suy biến giới hạn. Khi cho hai bán kính Ra bằng Rb, công thức tự động thu về dạng hàm hấp thu trên mặt cầu. Khi cho Rb tiến tới vô cùng, phương trình hoàn toàn trùng khớp với hàm hấp thụ mặt trụ. Khi cả Ra và Rb đều tiến tới vô cùng, kết quả suy biến chính xác về công thức mặt phẳng kinh điển của Cullity và Lê Chí Cương.

Thảo luận kết quả

Sự sai khác đáng kể giữa đặc tuyến hấp thụ trên mặt phẳng và mặt cong bắt nguồn từ cơ chế hình học quang tuyến. Trên mặt cong Ellipsoid, khoảng cách từ bề mặt mẫu đến phân tố tán xạ thay đổi phi tuyến theo vị trí điểm chiếu, khiến tổng quãng đường truyền tia AB cộng BC bị co ngắn hoặc kéo dài tùy theo góc quay quang học. Điều này làm cho thể tích tán xạ hiệu dụng thực tế tham gia vào quá trình tạo đỉnh phổ không còn đồng nhất như trên mặt phẳng chuẩn.

So với công bố của Taizo Oguri chỉ khảo sát hạn hẹp trong dải góc hai theta từ 152 độ đến 170 độ và giới hạn độ sâu thẩm thấu mang 96% năng lượng, mô hình nghiên cứu này bao quát trọn vẹn toàn bộ thể tích tán xạ trong dải góc rộng từ 80 độ đến 180 độ, loại bỏ hoàn toàn sai số cắt cụt thể tích. Dữ liệu nghiên cứu được trình bày một cách trực quan thông qua các hệ đồ thị biểu diễn tương quan giữa hệ số hấp thụ A và góc hai theta theo từng họ đường cong góc nghiêng psi. Đồng thời, bảng tổng hợp sai số định lượng cho phép các kỹ sư dễ dàng tra cứu và hiệu chỉnh vị trí đỉnh nhiễu xạ cực đại một cách chuẩn xác.

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

Thứ nhất, tích hợp trực tiếp thuật toán hàm hấp thu Ellipsoid vào vi chương trình điều khiển của các thiết bị đo nhiễu xạ tia X chuyên dụng. Mục tiêu giảm sai số tính toán ứng suất dư từ mức 12% hiện hành xuống dưới 3%, hướng tới hoàn thiện triển khai thử nghiệm trên các hệ thống đo lường công nghiệp trong lộ trình 6 đến 12 tháng do các nhóm nghiên cứu phát triển thiết bị chủ trì.

Thứ hai, chuẩn hóa quy trình đo lường không phá hủy cho các chi tiết cơ khí dạng trục cong, bề mặt bi cầu, rãnh then và bánh răng chịu tải trọng cao. Việc áp dụng chuẩn hiệu chỉnh mới sẽ giúp nâng tỷ lệ phát hiện ứng suất dư bất lợi lên mức 99%, bắt đầu áp dụng tại các trung tâm kiểm định chất lượng và nhà máy cơ khí chính xác từ quý 1 năm 2025.

Thứ ba, phát triển module phần mềm tự động trích xuất bán kính cong cục bộ chi tiết thông qua công nghệ quét laser 3D trước khi định vị chùm tia X. Giải pháp này giúp cắt giảm 50% thời gian thiết lập thông số đo đạc ban đầu, do các kỹ sư cơ điện tử và chuyên gia đo lường quang học phối hợp hoàn thiện trong thời gian 18 tháng.

Thứ tư, mở rộng cơ sở dữ liệu hệ số hấp thụ cho hơn 10 nhóm vật liệu hợp kim công nghiệp phổ biến khác như hợp kim titan, nhôm hàng không và thép hợp kim cao cấp trong giai đoạn 2025 đến 2027, tạo tiền đề vững chắc cho việc ứng dụng phương pháp nhiễu xạ tia X vào các ngành công nghiệp hàng không và năng lượng.

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

Nhóm kỹ sư R&D và chuyên gia kiểm tra không phá hủy (NDT) trong ngành công nghiệp ô tô, hàng không và đóng tàu: Tài liệu cung cấp cơ sở toán học chuẩn xác để hiệu chuẩn máy đo XRD, phục vụ kiểm soát chất lượng ứng suất dư sau gia công nhiệt luyện cho các trục khuỷu, vòng bi và cánh tuabin phức tạp.

Nhóm các nhà khoa học, nghiên cứu sinh và học viên cao học chuyên ngành Kỹ thuật Cơ khí, Khoa học Vật liệu: Luận văn là tài liệu tham khảo chuyên sâu về phương pháp giải tích hình học không gian, mô hình tương tác bức xạ tia X với tinh thể và kỹ thuật lập trình mô phỏng số bằng phần mềm Matlab.

Nhóm doanh nghiệp sản xuất và chế tạo thiết bị đo lường quang học: Cung cấp thuật toán cốt lõi để nâng cấp phần mềm phân tích phổ, mở rộng tính năng thương mại của các dòng máy đo nhiễu xạ đơn tinh thể sang phân khúc đo bề mặt cong đa dạng.

Nhóm giảng viên tại các trường đại học kỹ thuật công nghệ: Sử dụng làm học liệu giảng dạy cho các học phần chuyên đề ứng suất dư, vật lý chất rắn ứng dụng và kỹ thuật phân tích kết cấu vật liệu kim loại.

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

Tại sao yếu tố hấp thụ lại đóng vai trò quyết định nhất trong bộ ba LPA khi xác định ứng suất dư? Hệ số Lorentz và phân cực chủ yếu phụ thuộc vào góc nhiễu xạ và có dạng giải tích đơn giản ít biến động. Ngược lại, hàm hấp thụ phụ thuộc trực tiếp vào hình học biên dạng mẫu đo và góc nghiêng quang học, gây méo mó và dịch chuyển đỉnh phổ nhiễu xạ mạnh nhất nếu không được hiệu chỉnh đúng.

Phương pháp đo psi (side-inclination) có ưu thế gì so với phương pháp đo omega? Phương pháp đo psi thiết lập mặt phẳng nhiễu xạ vuông góc với hướng đo ứng suất, giúp ổn định chiều dài vết chiếu tia X trên bề mặt mẫu cong khi thay đổi góc nghiêng, từ đó hạn chế hiện tượng che khuất tia ở các góc quét lớn và đơn giản hóa mô hình tích phân hình học.

Sai số sẽ ra sao nếu tiếp tục áp dụng công thức mặt phẳng cho chi tiết cong Ellipsoid? Theo kết quả mô phỏng với hệ số hấp thụ 873, việc bỏ qua độ cong có thể gây sai lệch giá trị hấp thụ lên đến 18.5% ở góc nghiêng 30 độ, dẫn tới việc xác định sai lệch vị trí đỉnh hai theta và tạo ra sai số ứng suất dư từ 10% đến 20%.

Nghiên cứu này đã khắc phục điểm hạn chế cốt lõi nào của công trình Taizo Oguri? Mô hình trước đây chỉ khảo sát dải góc hẹp từ 152 độ đến 170 độ và giới hạn độ sâu dừng ở lớp hấp thụ 96% năng lượng. Nghiên cứu này quét trọn vẹn dải góc từ 80 độ đến 180 độ và tích phân toàn bộ thể tích tán xạ thực tế, loại trừ hoàn toàn sai số cắt cụt biên.

Công thức tổng quát trong đề tài có dùng được cho các loại bức xạ khác ngoài Cr-Kalpha không? Hoàn toàn có thể áp dụng. Người dùng chỉ cần thay đổi hệ số hấp thụ tuyến tính mu tương ứng với từng loại bức xạ như Cu-Kalpha, Co-Kalpha hoặc Mo-Kalpha và đặc tính nguyên tử của vật liệu đo vào hàm tích phân số đã được thiết lập.

Kết luận

  • Thiết lập thành công mô hình giải tích và công thức tích phân ba lớp xác định hàm hấp thu tổng quát A cho bề mặt Ellipsoid sử dụng phương pháp đo psi.
  • Chứng minh tính bao quát của công thức thông qua khả năng suy biến chính xác về các dạng hình học cơ bản gồm mặt cầu, mặt trụ và mặt phẳng.
  • Xác định định lượng quy luật biến thiên của hệ số hấp thụ trên nền vật liệu sắt với hệ số hấp thụ 873 trên dải góc quét từ 80 độ đến 180 độ và góc nghiêng từ 10 độ đến 30 độ.
  • Khẳng định mức độ ảnh hưởng của độ cong hình học với sai lệch hấp thụ cực đại lên đến 18.5%, giải quyết triệt để bài toán sai số tồn tại lâu năm trong các thiết bị đo phổ thương mại.
  • Đóng góp cơ sở lý luận và thuật toán thực tiễn vững chắc cho sự phát triển của lĩnh vực kiểm tra không phá hủy và đo lường ứng suất dư tại Việt Nam.

Công trình tạo tiền đề kỹ thuật quan trọng để nâng cấp hệ thống phần mềm phân tích phổ nhiễu xạ thế hệ mới trong giai đoạn 2025 đến 2026. Các kỹ sư cơ khí, chuyên gia vật liệu và nhà phát triển thiết bị đo lường được khuyến khích khai thác, ứng dụng bộ công thức tổng quát này vào quy trình kiểm chuẩn thực tế nhằm tối ưu hóa độ chính xác và nâng cao chất lượng sản phẩm chế tạo.