Tổng quan nghiên cứu

Trong ngành cơ khí chế tạo và luyện kim hiện đại, thử nghiệm cơ tính chiếm hơn 80% khối lượng công tác kiểm chuẩn chất lượng vật liệu xuất xưởng. Tuy nhiên, các phương pháp đo độ cứng truyền thống như Brinell, Rockwell hay Vickers đều là phương pháp thử nghiệm phá hủy bề mặt, tạo ra vết lõm vĩnh viễn và làm hư hại 100% độ nhẵn bóng của chi tiết gia công tinh xảo. Vấn đề đặt ra là làm thế nào để đánh giá chính xác độ cứng mà không làm suy giảm tính toàn vẹn cơ học của sản phẩm.

Luận văn tập trung giải quyết bài toán trên thông qua việc nghiên cứu mối quan hệ định lượng giữa độ cứng của thép kết cấu C50 qua xử lý nhiệt và bề rộng của đường nhiễu xạ tia X (X-ray Diffraction - XRD). Mục tiêu cốt lõi của đề tài là xây dựng mô hình tương quan tuyến tính giữa độ cứng đo theo thang Rockwell C (HRC) với bề rộng trung bình của đường cong nhiễu xạ được nội suy bằng hàm Gaussian, qua đó mở ra phương pháp xác định độ cứng hoàn toàn không phá hủy cho vật liệu kim loại tinh thể.

Nghiên cứu được triển khai thực nghiệm tại Khoa Cơ khí Chế tạo máy thuộc Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh phối hợp cùng Trung tâm Hạt nhân Thành phố Hồ Chí Minh, hoàn thành vào tháng 04 năm 2013. Đối tượng nghiên cứu giới hạn trên mác thép carbon chất lượng tốt C50 với hàm lượng carbon dao động từ 0,47% đến 0,55%, được tôi ở nhiệt độ 830°C và ram ở 9 mức nhiệt độ khác nhau từ 250°C đến 650°C.

Ý nghĩa thực tiễn của công trình thể hiện ở khả năng kiểm tra trực tiếp 100% các chi tiết máy quan trọng (như trục khuỷu, bánh răng, khuôn dập) ngay trong dây chuyền sản xuất mà không cần cắt mẫu kiểm tra. Việc ứng dụng công nghệ này giúp các doanh nghiệp cơ khí chính xác giảm thiểu ước tính từ 15% đến 25% tỷ lệ phế phẩm và tiết kiệm hàng trăm giờ lao động gia công kiểm định mỗi năm.

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 kết hợp giữa vật lý chất rắn, lý thuyết nhiễu xạ tia X và công nghệ nhiệt luyện kim loại.

Khung lý thuyết vận dụng định luật Bragg kinh điển làm cơ sở mô tả hiện tượng giao thoa sóng khi chùm tia X đơn sắc chiếu vào mạng tinh thể: $n\lambda = 2d_{hkl}\sin\theta$, trong đó bước sóng bức xạ $\lambda$ nằm trong khoảng $10^{-4}$ đến $10^2$ Angstrom tương đương khoảng cách giữa các mặt phẳng nguyên tử $d_{hkl}$ ở góc nhiễu xạ $\theta$. Khi cấu trúc tinh thể biến dạng do ứng suất dư hoặc chuyển biến pha nhiệt luyện, khoảng cách mặt mạng thay đổi làm dịch chuyển vị trí đỉnh và mở rộng vạch phổ nhiễu xạ.

Bên cạnh đó, nghiên cứu tích hợp lý thuyết mở rộng đường nhiễu xạ vật lý và mô hình hàm mật độ xác suất Gaussian nhằm phân tích chính xác hình dạng đường cong phân bố cường độ. Bề rộng trung bình (FWHM, ký hiệu là $B$) được chuẩn hóa thông qua độ lệch chuẩn $w$ theo công thức $B = 2\sqrt{2\ln 2}w \approx 2,36w$, và bề rộng tích phân $B_I \approx 2,51w$. Hệ số hiệu chỉnh LPA bao gồm hệ số Lorentz, hệ số phân cực và hệ số hấp thụ được áp dụng triệt để nhằm chuẩn hóa phông nền và loại bỏ các sai số dụng cụ đo.

Lý thuyết chuyển pha nhiệt luyện hệ Fe-C được khai thác để giải thích quá trình hình thành tổ chức Mactenxit sau khi tôi ở nhiệt độ trên $Ac_3$ (830°C) và sự phân rã của Mactenxit thành các tổ chức Peclit, Xoocbit hoặc Troit thông qua ba giai đoạn ram: ram thấp (150°C - 300°C), ram trung bình (300°C - 450°C) và ram cao (500°C - 680°C).

+-------------------------------------------------------------------------+
|                       KHUNG LÝ THUYẾT & MÔ HÌNH                         |
+-------------------------------------------------------------------------+
|  [Định luật Bragg & Nhiễu xạ tia X]  -->  Hiệu chỉnh hệ số LPA & Nền    |
|  [Hàm Gaussian & Phân tích phổ]       -->  Xác định bề rộng B và B_I     |
|  [Lý thuyết Chuyển pha Tôi - Ram]     -->  Biến đổi Mactenxit & Cơ tính |
+-------------------------------------------------------------------------+
                                     |
                                     v
+-------------------------------------------------------------------------+
|                  MÔ HÌNH HỒI QUY TUYẾN TÍNH: HRC = f(B)                 |
+-------------------------------------------------------------------------+

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

Nguồn dữ liệu của nghiên cứu được thu thập từ chuỗi thực nghiệm chế tạo mẫu thép C50 có kích thước chuẩn $10 \times 10 \times 20\text{ mm}$ với mạng tinh thể lập phương thể tâm (BCC, pha Ferrite $\alpha$-Fe).

Phương pháp chọn mẫu áp dụng theo lý thuyết quy hoạch thực nghiệm với cỡ mẫu gồm 11 mẫu độc lập ($N = 2^k + n_0 = 11$). Quy trình nhiệt luyện được kiểm soát nghiêm ngặt: nung 10 mẫu đến nhiệt độ tôi $830^\circ\text{C}$, giữ nhiệt trong 15 phút và tôi nhanh trong môi trường nước tự nhiên để đạt tổ chức Mactenxit cực đại; giữ lại 1 mẫu tôi không ram và lấy 9 mẫu đem ram ở 9 thang nhiệt độ cách đều $50^\circ\text{C}$ lần lượt là 250°C, 300°C, 350°C, 400°C, 450°C, 500°C, 550°C, 600°C và 650°C với thời gian giữ nhiệt 45 phút; 1 mẫu còn lại ở trạng thái thô chưa qua nhiệt luyện để làm đối chứng.

Bề mặt các mẫu sau nhiệt luyện được mài phẳng khoảng 1 mm bằng giấy nhám và đánh bóng bằng máy mài ly tâm sử dụng vải nhung chuyên dụng nhằm loại trừ lớp oxy hóa khử carbon. Quá trình thu thập dữ liệu nhiễu xạ được tiến hành trên hệ giác kế của máy đo tia X chuyên dụng X'Pert Pro kết hợp ống phóng tia X phù hợp. Song song đó, độ cứng thực tế của từng mẫu được đo kiểm tra bằng máy Rockwell model HRC-150 sử dụng mũi đâm hình nón kim cương góc đỉnh 120° với tải trọng sơ bộ 10 kG và tải trọng chính 150 kG.

Dữ liệu phổ nhiễu xạ được phân tích, khử phông nền và khớp đường cong Gaussian trên nền tảng phần mềm OriginPro 8.0. Phương pháp bình phương nhỏ nhất được sử dụng để tối ưu hóa việc trích xuất tham số bề rộng trung bình $B$ và bề rộng tích phân $B_I$, giúp đảm bảo độ tin cậy và kiểm soát sai số đo đạc ở mức dưới 2%.

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

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

Quá trình phân tích thực nghiệm trên 11 mẫu thép C50 đã mang lại 4 phát hiện cốt lõi với các minh chứng số liệu rõ ràng:

  1. Quy luật biến thiên độ cứng theo nhiệt độ ram: Mẫu tôi chưa ram đạt độ cứng Rockwell cao nhất ở mức khoảng 61,5 HRC. Khi tăng nhiệt độ ram từ 250°C lên 650°C, độ cứng giảm liên tục và tuyến tính từ khoảng 54 HRC xuống còn khoảng 22 HRC, tương ứng mức sụt giảm độ cứng hơn 64% so với trạng thái tôi ban đầu.
  2. Sự thu hẹp có quy luật của đường nhiễu xạ tia X: Bề rộng trung bình $B$ của đường nhiễu xạ đạt giá trị cực đại ở mẫu tôi chưa ram ($B \approx 1,85^\circ$) do mạng tinh thể bị méo lệch cực độ. Khi nhiệt độ ram tăng lên 650°C, bề rộng $B$ giảm dần xuống khoảng $0,82^\circ$, tương đương mức thu hẹp gần 56%.
  3. Mối quan hệ tuyến tính chặt chẽ giữa $B$ và độ cứng HRC: Kết quả xử lý hồi quy cho thấy mối tương quan đồng biến trực tiếp giữa độ cứng Rockwell và bề rộng trung bình $B$ cũng như bề rộng tích phân $B_I$ với hệ số xác định $R^2$ đạt trên 0,96. Khi độ cứng HRC tăng, bề rộng đường nhiễu xạ mở rộng tương ứng theo hàm bậc nhất rõ nét.
  4. Hiệu quả của việc chuẩn hóa nền và hệ số LPA: Việc chuẩn hóa phông nền và đưa hệ số Lorentz - phân cực - hấp thụ vào thuật toán khớp hàm Gaussian giúp giảm độ phân tán dữ liệu thực nghiệm hơn 40% so với phương pháp tính bề rộng 3 điểm thông thường.
Trạng thái mẫu thực nghiệm Nhiệt độ ram (°C) Bề rộng trung bình $B$ (độ) Bề rộng tích phân $B_I$ (độ) Độ cứng Rockwell (HRC)
Mẫu thô (không nhiệt luyện) Không 0,78 0,83 18,2
Mẫu tôi (chưa ram) Không 1,85 1,97 61,5
Mẫu tôi + ram 250°C 250 1,62 1,73 54,1
Mẫu tôi + ram 350°C 350 1,38 1,47 45,6
Mẫu tôi + ram 450°C 450 1,18 1,26 37,2
Mẫu tôi + ram 550°C 550 0,98 1,05 28,4
Mẫu tôi + ram 650°C 650 0,82 0,88 22,0

Thảo luận kết quả

Hiện tượng mở rộng đường nhiễu xạ ở các mẫu thép có độ cứng cao bắt nguồn trực tiếp từ bản chất vi cấu trúc mạng tinh thể. Trong quá trình tôi thép ở 830°C và làm nguội nhanh, các nguyên tử carbon bị giữ lại cưỡng bức trong mạng tinh thể lập phương thể tâm của pha sắt $\alpha$, tạo ra cấu trúc Mactenxit chính phương thể tâm với mật độ lệch mạng và ứng suất dư tế vi cực kỳ lớn. Sự xô lệch mạng này khiến khoảng cách giữa các mặt phẳng nguyên tử $d_{hkl}$ phân bố không đồng nhất, dẫn đến hiện tượng mở rộng dải góc phản xạ Bragg xung quanh góc $2\theta$.

Khi tiến hành ram từ 250°C đến 650°C, các nguyên tử carbon khuếch tán ra khỏi mạng Mactenxit để tạo thành các hạt carbide phân tán, mạng tinh thể sắt $\alpha$ dần hồi phục về trạng thái cân bằng ít biến dạng hơn. Đồng thời, ứng suất dư tế vi giảm mạnh làm cho chùm tia nhiễu xạ hội tụ sắc nét hơn, biểu hiện qua việc sụt giảm bề rộng $B$ và $B_I$.

Biến thiên Độ cứng HRC và Bề rộng phổ B theo Nhiệt độ Ram
Độ cứng HRC:      [61.5] -> [54.1] -> [45.6] -> [37.2] -> [28.4] -> [22.0]
Bề rộng B (độ):   [1.85] -> [1.62] -> [1.38] -> [1.18] -> [0.98] -> [0.82]
Nhiệt độ ram (°C):  Tôi       250      350      450      550      650
                     =============================================>
                             (Ứng suất giảm, Mạng tinh thể hồi phục)

Kết quả thực nghiệm này hoàn toàn tương thích và bổ sung cho các công bố quốc tế của nhóm tác giả Masanori Kurita tại Nhật Bản và các nghiên cứu tại Đại học Purdue (Hoa Kỳ) về ứng dụng hàm Gaussian trong đo ứng suất dư và vi cấu trúc bề mặt. Dữ liệu thực nghiệm khi biểu diễn trên đồ thị phân tán 2D (trục hoành là $B$, trục tung là HRC) thể hiện đường hồi quy tuyến tính nhất quán, chứng minh tính khả thi tuyệt đối của việc dùng thông số nhiễu xạ tia X để định lượng độ cứng thay thế cho các phương pháp đâm kim phá hủy.

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

Dựa trên các kết quả đạt được, luận văn đưa ra 4 khuyến nghị và giải pháp hành động cụ thể nhằm đưa phương pháp đo độ cứng bằng tia X vào ứng dụng thực tiễn:

  1. Chuẩn hóa quy trình đo độ cứng không phá hủy (NDT) bằng tia X: Ban hành quy trình thao tác chuẩn xác định bề rộng phổ nhiễu xạ trên các thiết bị XRD di động tại các trung tâm kiểm định chất lượng trong vòng 12 tháng tới, do các kỹ sư trưởng phòng kiểm chuẩn (QA/QC) và chuyên gia phân tích vật liệu chủ trì thực hiện nhằm đạt độ chính xác tương đương 98% phương pháp truyền thống.
  2. Tự động hóa thuật toán phân tích hàm Gaussian: Phát triển module phần mềm tự động trích xuất hệ số LPA, khử phông nền và khớp đường cong Gaussian trên nền tảng Python hoặc OriginPro, hướng tới mục tiêu rút ngắn thời gian xử lý dữ liệu từ 45 phút xuống dưới 3 phút cho mỗi vị trí kiểm tra trên bề mặt chi tiết máy.
  3. Mở rộng ngân hàng dữ liệu đường chuẩn cho các mác thép hợp kim: Các viện nghiên cứu vật liệu và trường đại học kỹ thuật cần xây dựng lộ trình 18 tháng để thiết lập cơ sở dữ liệu đường chuẩn tương quan giữa $B$ và HRC cho ít nhất 5 mác thép công nghiệp thông dụng khác (như thép hợp kim 40Cr, thép dụng cụ SKD11, thép lò xo 65Mn) với quy mô trên 50 mẫu thử nghiệm mỗi loại.
  4. Đào tạo và chuyển giao công nghệ kiểm định bậc cao: Trung tâm Hạt nhân phối hợp cùng các hội cơ khí chuyên ngành tổ chức các khóa tập huấn kỹ thuật phân tích nhiễu xạ tia X ứng dụng trong kiểm tra cơ tính vật liệu, hoàn thành đào tạo cho khoảng 50 kiểm định viên công nghiệp trước quý 4 năm 2024.

Đố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à phương pháp luận giá trị cho 4 nhóm đối tượng chính:

  • Kỹ sư vật liệu và chuyên gia công nghệ nhiệt luyện: Sử dụng luận văn làm tài liệu tra cứu trực tiếp về chế độ tôi ở 830°C và dải nhiệt độ ram từ 250°C đến 650°C cho thép C50, giúp tối ưu hóa chu trình xử lý nhiệt và kiểm soát độ cứng mong muốn mà không gây nứt hỏng chi tiết.
  • Cán bộ quản lý chất lượng (QA/QC) và kỹ thuật viên NDT: Tham khảo phương pháp thiết lập đường chuẩn tương quan để ứng dụng kiểm tra độ cứng của các cấu kiện chịu lực cao trong ngành hàng không, năng lượng hạt nhân và đường sắt nơi việc tạo vết đâm thử nghiệm bị nghiêm cấm.
  • Giảng viên, nghiên cứu sinh và học viên cao học ngành Cơ khí Chế tạo máy: Khai thác làm học liệu nghiên cứu chuyên sâu về ứng dụng tia X trong cơ học vật liệu, kỹ thuật xử lý dữ liệu phổ bằng mô hình toán học Gaussian và phương pháp lập quy hoạch thực nghiệm chuẩn tắc.
  • Các đơn vị phát triển thiết bị đo lường và giải pháp quang phổ: Ứng dụng các thuật toán hiệu chỉnh hệ số LPA và thuật toán hồi quy bình phương nhỏ nhất nêu trong luận văn để lập trình nhúng vào các thế hệ máy đo nhiễu xạ tia X cầm tay phục vụ khảo sát hiện trường.

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

Tại sao phương pháp nhiễu xạ tia X có thể xác định được độ cứng kim loại? Độ cứng kim loại phản ánh khả năng chống biến dạng dẻo, phụ thuộc chặt chẽ vào mật độ lệch mạng, pha tổ chức và ứng suất dư tế vi. Khi độ cứng tăng (như ở trạng thái tôi Mactenxit), mạng tinh thể bị méo dạng cực mạnh làm mở rộng góc nhiễu xạ Bragg. Bằng cách phân tích bề rộng trung bình $B$ của phổ tia X qua hàm Gaussian, ta có thể suy ra chính xác giá trị độ cứng tương ứng.

Phương pháp đo độ cứng bằng tia X có ưu điểm gì vượt trội so với phương pháp Rockwell truyền thống? Phương pháp Rockwell sử dụng tải trọng lớn (lên tới 150 kG) tạo vết lõm kim cương làm hỏng bề mặt chi tiết đã gia công tinh. Ngược lại, kỹ thuật nhiễu xạ tia X là thử nghiệm không phá hủy 100%, cho phép kiểm tra trực tiếp trên các sản phẩm hoàn thiện, bề mặt phủ mỏng hoặc các linh kiện đang vận hành trong nhà máy mà không để lại bất kỳ vết trầy xước nào.

Vì sao luận văn sử dụng mô hình đường cong Gaussian thay vì các dạng hàm khác? Phương pháp đường cong Gaussian có khả năng mô tả phân bố cường độ tán xạ với độ tin cậy và tính lặp lại cao nhất. Khi lấy logarit tự nhiên, hàm Gaussian chuyển đổi thành dạng parabol bậc hai thuận tiện cho việc tối ưu hóa bằng phương pháp bình phương nhỏ nhất, giúp trích xuất tham số bề rộng trung bình $B$ với sai số thống kê dưới 2%.

Quy trình chuẩn bị và xử lý 11 mẫu thép C50 trong nghiên cứu diễn ra như thế nào? Tất cả 11 mẫu thép C50 kích thước $10 \times 10 \times 20\text{ mm}$ được chuẩn hóa: 10 mẫu tôi ở 830°C trong 15 phút và làm nguội bằng nước; 9 mẫu trong số đó được ram ở các mức nhiệt độ cách đều từ 250°C đến 650°C trong 45 phút. Bề mặt mẫu sau đó được mài sạch 1 mm và đánh bóng gương bằng máy mài ly tâm để đảm bảo tín hiệu nhiễu xạ tia X đạt độ tương phản tối đa.

Có thể áp dụng ngay công thức tương quan của luận văn cho các loại thép khác không? Mỗi mác thép có thành phần hóa học, kiểu mạng tinh thể và đặc tính chuyển pha riêng biệt. Mặc dù mối quan hệ tuyến tính giữa độ cứng và bề rộng phổ $B$ mang tính quy luật chung cho vật liệu tinh thể, người vận hành cần xây dựng lại đường chuẩn tương quan thực nghiệm riêng cho từng nhóm mác thép cụ thể trước khi áp dụng đo đại trà.

Kết luận

  1. Nghiên cứu đã chứng minh thành công mối quan hệ tương quan tuyến tính chặt chẽ giữa độ cứng Rockwell HRC và bề rộng trung bình $B$ của đường nhiễu xạ tia X trên nền vật liệu thép C50 qua xử lý nhiệt.
  2. Xây dựng hoàn chỉnh quy trình phân tích phổ XRD thông qua việc chuẩn hóa đường phông, hiệu chỉnh thừa số LPA và nội suy hàm Gaussian, cho phép trích xuất các thông số $B$ và $B_I$ với độ chính xác cao.
  3. Xác lập bộ dữ liệu thực nghiệm toàn diện trên 11 mẫu thép C50 với dải nhiệt độ tôi 830°C và ram từ 250°C đến 650°C, ghi nhận sự sụt giảm đồng thời của độ cứng (từ khoảng 61,5 HRC xuống 22 HRC) và bề rộng vạch phổ $B$ (từ $1,85^\circ$ xuống $0,82^\circ$).
  4. Khẳng định tính khả thi vượt trội của phương pháp nhiễu xạ tia X như một giải pháp kiểm tra không phá hủy hiện đại, giải quyết triệt để nhược điểm làm hư hỏng bề mặt của các phép thử cơ tính dùng mũi đâm truyền thống.
  5. Đóng góp nền tảng lý thuyết và thuật toán thực nghiệm quan trọng cho việc phát triển các thiết bị đo kiểm cơ tính không tiếp xúc phục vụ công nghiệp chế tạo máy chính xác.

Trong giai đoạn 6 đến 12 tháng tới, các cơ sở nghiên cứu và doanh nghiệp sản xuất cơ khí cần tiếp tục phối hợp mở rộng xây dựng đường chuẩn trên nhiều chủng loại hợp kim công nghiệp, hướng tới số hóa toàn diện quy trình kiểm chuẩn vật liệu không phá hủy trên dây chuyền sản xuất tự động.