Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.Nghiên cứu biến dạng đàn hồi – phi tuyến của kim loại, hợp kim xen kẽ hai và ba thành phần.
BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM HÀ NỘI NGUYỄN ĐỨC HIỀN NGHIÊN CỨU BIẾN DẠNG ĐÀN HỒI – PHI TUYẾN CỦA KIM LOẠI, HỢP KIM XEN KẼ HAI VÀ BA THÀNH PHẦN Chuyên ngành: Vật lí lý thuyết Vật lí tốn Mã số: 9.44.01.03 TĨM TẮT LUẬN ÁN TIẾN SĨ VẬT LÍ HÀ NỘI - 2023 LUẬN ÁN ĐƯỢC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC SƯ PHẠM HÀ NỘI Người hướng dẫn khoa học: PGS TS Nguyễn Quang Học PGS.TS Hồng Văn Tích Phản biện 1: GS TS Bạch Thành Công – Trường Đại học KHTN – ĐHQG Hà Nội Phản biện 2: GS TS Vũ Văn Hùng – Trường Đại học Giáo dục – ĐHQG Hà Nội Phản biện 3: PGS.TS Nguyễn Hồng Quang – Viện Vật lý – Viện Hàn lâm KH&CN Việt Nam Luận án bảo vệ Hội đồng đánh giá luận án cấp Trường họp Trường Đại học Sư phạm Hà Nội vào hồi ngày tháng năm Có thể tìm hiểu luận án tại: - Thư viện Quốc gia Việt Nam - Thư viện Trường Đại học Sư phạm Hà Nội MỞ ĐẦU Lý chọn đề tài Kim loại hợp kim vật liệu có lịch sử phát triển lâu đời sử dụng phổ biến ngành cơng nghiệp sống thực tiễn hàng ngày Có hai loại hợp kim hợp kim thay (HKTT) hợp kim xen kẽ (HKXK) Thép dạng HKXK điển hình có tầm quan trọng đặc biệt xây dựng, giao thông vận tải, chế tạo máy Một HKXK Fe FeC gọi thép cacbon chiếm tỷ trọng lớn ngành công nghiệp thép Fe HKXK FeSi, FeC, FeH chiếm phần lớn lõi Trái Đất thiên thể Các tính chất nhiệt động, đàn hồi, nóng chảy chúng cung cấp cho thông tin thành phần, cấu trúc, tiến hóa…của Trái Đất thiên thể Những thông tin kim loại Au, Cu, Fe, Cr…và số hợp kim chúng có ý nghĩa đặc biệt quan trọng để thiết kế chế tạo thiết bị phục vụ nhu cầu sống người Việc nghiên cứu tính chất học bao gồm trình biến dạng loại vật liệu kim loại, hợp kim có HKXK thu hút quan tâm đặc biệt nhiều nhà nghiên cứu lý thuyết thực nghiệm Các kết nghiên cứu đóng vai trị quan trọng khoa học, công nghệ sống người Nhiều cơng trình dùng PPTKMM làm cơng cụ nghiên cứu cơng bố tạp chí khoa học có uy tín ngồi nước điều khẳng định thành công PPTKMM Hiệu bật PPTKMM thể hai điểm Thứ PPTKMM giải tốt tốn phi điều hịa tinh thể Thứ hai sử dụng PPTKMM, đại lượng nhiệt động đàn hồi mô tả biểu thức giải tích thuận tiện cho tính số dễ dàng so sánh kết tính số theo PPTKMM với số liệu thực nghiệm kết tính tốn theo phương pháp khác Trước đây, tác giả Hung Hoa sử dụng PPTKMM để xây dựng lý thuyết biến dạng đàn hồi phi tuyến kim loại HKTT hai thành phần [4,59-61] Trong cơng trình [51], Hoc cộng sử dụng PPTKMM để nghiên cứu ảnh hưởng nhiệt độ nồng độ nguyên tử xen kẽ đến tính chất nhiệt động HKXK hai thành phần có cấu trúc LPTK Tuy nhiên, việc nghiên cứu biến dạng HKXK PPTKMM tập trung vào biến dạng đàn hồi, vận tốc sóng đàn hồi, tính chất biến dạng phi tuyến ảnh hưởng biến dạng đến khuếch tán chưa nghiên cứu Với tất lý kể trên, lựa chọn đề tài luận án “Nghiên cứu biến dạng đàn hồi-phi tuyến kim loại, hợp kim xen kẽ hai ba thành phần” Đây lĩnh vực nhiều vấn đề chưa nghiên cứu hoàn thiện PPTKMM đồng thời kết nghiên cứu phương pháp lý thuyết khác cịn sai lệch nhiều so với thực nghiệm Chính vậy, việc nghiên cứu biến dạng kim loại HKXK nhiều thành phần khơng có ý nghĩa mặt lý thuyết mà cịn có tính thời khoa học giá trị thực tiễn Mục đich, đối tượng phạm vi nghiên cứu Nghiên cứu biến dạng đàn hồi – phi tuyến, ảnh hưởng biến dạng lên khuếch tán HKXK hai ba thành phần với cấu trúc LPTD LPTK tính đến ảnh hưởng nhiệt độ, áp suất nồng độ nguyên tử xen kẽ Chúng giới hạn phạm vi nghiên cứu khoảng nhiệt độ từ 0K đến nhiệt độ nóng chảy kim loại hợp kim, khoảng nồng độ nguyên tử thay từ đến 10%, khoảng nồng độ nguyên tử xen kẽ từ đến 5% khoảng áp suất lên tới 100 GPa Phạm vi nghiên cứu thuận lợi cho việc tính số dễ so sánh với kết nghiên cứu thực nghiệm kết tính tốn phương pháp lý thuyết khác Phương pháp nghiên cứu Phương pháp nghiên cứu phương pháp thống kê mômen số phương pháp tính gần phương pháp cầu phối vị, phương pháp gần phân bố nồng độ nguyên tử, Sử dụng phần mềm Maple, Origin Graph Grabber để lập trình, tính số vẽ đồ thị Các kết tính tốn lí thuyết lập bảng vẽ đồ thị nhằm đánh giá kết nhận được, so sánh với thực nghiệm kết tính tốn phương pháp lý thuyết khác Ý nghĩa khoa học thực tiễn luận án Luận án cung cấp nhiều thông tin biến dạng đàn-dẻo ảnh hưởng biến dạng lên khuếch tán HKXK Kết giải tích cho biết phụ thuộc nhiệt độ, áp suất, nồng độ thành phần hợp kim môđun đàn hồi, số đàn hồi, vận tốc sóng đàn hồi, mật độ lượng biến dạng, ứng suất thực cực đại, giới hạn biến dạng đàn hồi hệ số khuếch tán Bên cạnh đó, luận án cịn phụ thuộc độ biến dạng hệ số khuếch tán HKXK Luận án góp phần phát triển PPTKMM nghiên cứu tính chất vật liệu nói chung nghiên cứu biến dạng đàn-dẻo, ảnh hưởng biến dạng lên khuếch tán kim loại HKXK nói riêng Một số kết tính số luận án tư liệu tham khảo nhằm định hướng cho nghiên cứu thực nghiệm đồng thời mở khả nghiên cứu biến dạng đàn-dẻo kim loại HKXK có cấu trúc khác ví dụ LGXC Những đóng góp luận án Luận án xây dựng biểu thức giải tích đại lượng đặc trưng biến dạng đàn hồi phi tuyến môđun đàn hồi, số đàn hồi, vận tốc sóng đàn hồi, mật độ lượng biến dạng, ứng suất thực cực đại giới hạn biến dạng đàn hồi phụ thuộc vào nhiệt độ, áp suất nồng độ thành phần HKXK hai ba thành phần Luận án xây dựng biểu thức giải tích hệ số khuếch tán phụ thuộc vào nhiệt độ, áp suất, nồng độ thành phần độ biến dạng HKXK hai ba thành phần Một số kết tính số đại lượng biến dạng PPTKMM gần thực nghiệm so với kết tính số phương pháp lý thuyết khác Các kết tính số khác PPTKMM dùng để dự đoán, định hướng thực nghiệm tương lai Cấu trúc luận án Ngoài phần Mở đầu, Kết luận, Tài liệu tham khảo phụ lục nội dung luận án trình bày chương Nội dung luận án báo cáo hội nghị khoa học chuyên ngành quốc gia quốc tế Hội nghị Vật lí lý thuyết Quốc gia lần thứ 42 (NCTP42), Cần Thơ, 31/7-3/8/2017 Hội nghị Vật lí chất rắn khoa học vật liệu Quốc gia lần thứ 10 (SPMS2017), Huế, 19-21/10/2017 Hội nghị Quốc tế Hiệp hội Khoa học vật liệu tính tốn châu Á Mơ hình hóa nhiều cấp độ vật liệu phát triển bền vững (ACCMS-TM 2018), Đại học Quốc gia Hà Nội – 9/9/2018 công bố báo tạp chí khoa học nước quốc tế CHƯƠNG 1: TỔNG QUAN VỀ BIẾN DẠNG CỦA KIM LOẠI VÀ HỢP KIM 1.1 Hợp kim xen kẽ Trong HKXK, kim loại thành phần quan trọng nhất, thường chiếm nồng độ từ 90% trở lên gọi kim loại (kim loại mẹ hay kim loại sở) Ngồi kim loại chính, HKXK cịn có thành phần khác thường phi kim Đó tác nhân tạo hợp kim chiếm nồng độ nhỏ Đôi HKXK hợp chất xen kẽ nguyên tố tạo thành hợp chất liên kết hóa học Hợp chất xen kẽ thường dung dịch rắn nguyên tử nguyên tố hỗn hợp với giống chất tan dung môi Các phi kim quan trọng HKXK H, B, C, N, O, Si, Cấu trúc tinh thể HKXK cấu trúc tinh thể kim loại định Các HKXK hai thành phần loại hợp kim mà thành phần kim loại cịn thành phần khác ngun tử phi kim có kích thước nhỏ nằm khoảng trống nút mạng nguyên tử kim loại Trong thực tế, hợp kim thường có nhiều thành phần với nồng độ khác HKXK ba thành phần loại hợp kim mà ngồi thành phần kim loại chính, cịn thành phần thứ hai ngun tử phi kim có kích thước nhỏ xen kẽ nút mạng nguyên từ kim loại thành phần thứ ba nguyên tử phi kim xen kẽ khác nguyên tử kim loại khác thay số vị trí nguyên tử kim loại 1.2 Lý thuyết biến dạng 1.2.1 Biến dạng đàn hồi Dưới tác dụng ngoại lực, vật rắn bị biến dạng, tức bị biến đổi hình dạng kích thước Khi biến dạng, điểm hay nguyên tử vật rắn dịch chuyển Biến dạng đàn hồi biến dạng vật rắn tác dụng ngoại lực mà sau cất tải, vật rắn trở lại hình dạng, kích thước ban đầu 1.2.2 Biến dạng phi tuyến Biến dạng phi tuyến (biến dạng phi đàn hồi, biến dạng dư) biến dạng vật rắn tác dụng ngoại lực mà sau cất tải, biến dạng không bị vật rắn khơng trở lại hình dạng, kích thước ban đầu Điều xảy ngoại lực (tải) phải đủ lớn Biến dạng phi tuyến khơng làm thay đổi thể tích vật biến dạng 1.2.3 Sóng đàn hồi vật rắn Các vận tốc sóng dọc sóng ngang có dạng Vd = 2C44 + C12 C44 , Vn = (1.28) 1.2.4 Ảnh hưởng biến dạng lên khuếch tán Sự phụ thuộc hệ số khuếch tán D vào ứng suất kéo lưỡng trục có dạng 2 r m m V + V − V/ / , Dx ( ) = Dx ( = 0) exp (1.31) kBT Dx ( ) hệ số khuếch tán theo phương x hệ chịu tác dụng ứng suất lưỡng trục , V r thể tích hồi phục V/m/ thành phần thể tích dịch chuyển theo phương song song với phương khuếch tán 1.3 Một số phương pháp nghiên cứu chủ yếu Có nhiều phương pháp lý thuyết khác nghiên cứu biến dạng đàn hồi phi tuyến kim loại hợp kim lý thuyết điều hòa, lý thuyết chuẩn điều hòa, phương pháp động lực học phân tử, phương pháp phần tử hữu hạn, phương pháp tính tốn từ nguyên lí (ab initio), phương pháp Hamiltonien liên kết chặt, phương pháp lý thuyết phiếm hàm mật độ, phương pháp học máy, phương pháp hàm Green mạng, phương pháp tính tốn giản đồ pha, phương pháp ngun tử nhúng biến dạng, 1.4 Phương pháp thống kê mômen Kˆ n +1 a = Kˆ n a Qˆ n +1 a + Kˆ n an +1 a − B2m i 2m )! m =0 ( 2m Kˆ n (2m) an +1 (1.57 ) a Kết luận chương Trong chương 1, chúng tơi trình bày tổng quan lý thuyết biến dạng phương pháp nghiên cứu biến dạng đàn hồi - phi tuyến kim loại hợp kim, đề cập đến lý thuyết biến dạng đàn hồi, biến dạng phi tuyến, vận tốc sóng đàn hồi ảnh hưởng biến dạng lên khuếch tán vật liệu; nội dung ưu nhược điểm phương nghiên cứu biến dạng phổ biến Chương cịn giới thiệu cơng thức truy chứng liên hệ mômen bậc cao qua mômen bậc thấp dùng để xác định lượng tự hệ CHƯƠNG 2: NGHIÊN CỨU BIẾN DẠNG ĐÀN HỒI VÀ PHI TUYẾN CỦA HỢP KIM XEN KẼ AC 2.1 Mơ hình hợp kim lượng tự Helmholtz Hình 2.1 Mơ hình HKXK AC với cấu trúc LPTK (a) LPTD(b) Năng lượng tự Helmholtz HKXK AC có dạng AC = c X X − TScAC , (2.1) X X = U0 X + X + k X + 2 k X4 4 X YX 3 YX 1 + X = 3 x X 2 1X X YX − YX 1 + + YX − 1X + 2 1X X 1 + (1 + YX ) , ( ) + ln (1 − e ) , Y x coth x −2 x X (2.2) X cA = − 7cC , cA1 = 2cC , cA2 = 4cC mạng LPTK, cA = − 15cC , cA1 = 6cC , cA2 = 8cC mạng LPTD X = A, C, A1, A2 , 2.2 Năng lượng liên kết, thông số hợp kim khoảng gần trung bình hai nguyên tử Đối với mạng LPTK, u0 A = 4 AA ( r1A ) + 3 AA ( r2 A ) , r2 A = kA = r1 A , (2.4) 2 d AA ( r1A ) d AA ( r1A ) d AA ( r2 A ) + + + 3r1A dr1A dr12A dr22A + d AA ( r2 A ) , r2 A dr2 A (2.5) d AA ( r1A ) d AA ( r1A ) d AA ( r1 A ) d AA ( r1 A ) + + + 2 54 9r1A dr1A dr1A dr1A 9r1A dr1A 9r1A 1A = + d AA ( r2 A ) d AA ( r2 A ) d AA ( r2 A ) + - , 24 dr2 A dr24A 4r22A dr22A 4r2 A (2.6) d AA ( r1A ) d AA ( r1A ) d AA ( r1A ) d AA ( r2 A ) + + ( 2.7 ) dr1A 2r2 A dr14A 3r12A dr12A 3r13A dr23A 2A = u0C = AC ( r1C ) + 2 AC ( r2C ) , r2C = 2r1C , d AC ( r1C ) d AC ( r2C ) d AC ( r2C ) + + , r1C dr1C r2C dr2C dr2C (2.8) kC = C = (1C + 2C ) , 1C = + (2.10) d AC (r1C ) d AC (r1C ) d AC (r2C ) - + + dr1C 48 8r12C dr12C 8r1C dr24C d 3 AC (r2C ) d 2 AC (r2C ) d AC (r2C ) + , 2 8r2C dr2C dr2C 16r2C dr2C 16r23C 2C = + (2.9) (2.11) d 3 AC (r1C ) d 2 AC (r1C ) d AC (r2C ) - + + 4r1C dr1C dr13C 2r1C dr12C 2r1C d 3 AC (r2C ) d 2 AC (r2C ) d AC (r2C ) − + , 2 4r2C dr2C dr2C 4r2C dr2C 4r2C ( ) u0 A1 = u0 A + 3 A1C r1 A1 , k A1 = k A + ( ) d 2 A1C r1A1 dr12A1 + ( (2.13) ( ) d A1C r1A1 , r1 A1 dr1 A1 ) A1 = 1A1 + A1 , 1A1 = 1A + (2.12) (2.14) (2.15) d A1C d A1C − + 24 dr14A 6r1 A1 dr13A 1 + d A1C d A1C − , 2 4r1A1 dr1A1 4r1A1 dr1A1 (2.16) d A1C (r1A1 ) , 4r1 A1 dr13A1 γ 2A1 = A + (2.17) ( ) u0 A2 = u0 A + 6 A2 C r1 A2 , k A2 = k A + ( )+ d 2 A2 C r1 A2 dr12A2 ( (2.18) ( ) d A2 C r1 A2 , r1 A2 dr1 A2 ) A2 = A2 + A2 , 1A2 = 1A + (2.19) (2.20) ( ) d A2 C r1A2 + 24 dr14A 2 15 d A2 C (r1A2 ) 15 d A2 C (r1A2 ) + − , dr1 A2 dr12A2 4r1A2 4r1A2 A2 = A + d A2 C (r1A2 ) d A2 C (r1A2 ) + + 4r1A2 dr14A2 dr13A2 + d A2 C (r1A2 ) d A2 C (r1A2 ) - , dr1 A2 dr12A2 8r12A2 8r1A2 Đối với mạng LPTD, u0 A = 6 AA ( r1A ) + 3 AA ( r2 A ) , r2 A = kA = 1A = (2.21) d 2 AA ( r1A ) dr12A + 2r1A , (2.22) (2.23) d AA ( r1 A ) d AA ( r2 A ) d AA ( r2 A ) + + , (2.24) r1A dr1A r2 A dr2 A dr22A d AA ( r1A ) d AA ( r1A ) d AA ( r1A ) d AA ( r1A ) + + + 2 24 4r1A dr1A dr1A dr1A 8r1A dr1A 8r1A + 2A = d AA ( r2 A ) d AA ( r2 A ) d AA ( r2 A ) + - , (2.25) 24 dr2 A dr24A 4r22A dr22A 4r2 A d AA ( r1 A ) d AA ( r1 A ) d AA ( r1 A ) d AA ( r1 A ) + + − + 2 4r1A dr1A dr1A dr1A 8r1A dr1A 8r1A + d AA ( r2 A ) d AA ( r2 A ) d AA ( r2 A ) + (2.26) 2 2r2 A dr2 A dr2 A 4r2 A dr2 A 4r2 A 1A2 = 1A + d A2 C (r1A2 ) d A2 C (r1A2 ) + 54 9r1A2 dr14A2 dr13A2 − A2 = A + 2 d A2 C (r1A2 ) d A2 C (r1A2 ) + , 2 dr1 A2 dr1A2 9r1A2 9r1A2 (2.40) d A2 C (r1A2 ) d A2 C (r1A2 ) d A2 C (r1A2 ) + - , ( 2.41) 2 dr1 A2 dr1A2 3r1A2 dr1A2 3r1A u0 k Pv = −r1 + θx coth x , 2k r1 r1 (2.42) u0 ω k Pv = −r1 + 4k r1 r1 r1C ( P,T ) = r1C ( P, 0) + y A1 (P,T ),r1A (P,T ) = r1A (P, 0) + y A (P,T ), r1A1 ( P, T ) = r1C ( P, T ), r1A2 ( P, T ) = r1A2 ( P, 0) + yC (P, T ) (2.43) (2.44) r1A ( P,T ) = r1A ( P, 0) + y( P,T ),r1A ( P, 0) = (1 − cC ) r1A (P, 0) + cC r1A (P, 0), y ( P,T ) = c A y A ( P,T ) + cC yC ( P,T ) + c A1 y A1 ( P,T ) + c A2 y A2 ( P,T ), r1A ( P, 0) = 3r1C ( P, 0) , c A = − 7cC , c A1 = 2cC , c A2 = 4cC (2.45) mạng LPTK, r1A ( P, 0) = r1C ( P, 0) ,c A = − 15cC ,c A1 = 6cC , cA2 = 8cC mạng LPTD 2.3 Biến dạng đàn hồi EYAC = EYA cX X 2 X 2 A , (2.51) Ta tìm mơđun đàn hồi khác, số đàn hồi vận tốc truyền sóng đàn hồi theo công thức (2.52) đến (2.59) 2.4 Biến dạng phi tuyến (2.61) r1FX ( P, 0) = r1X ( P, 0)(1 + ) , (2.62) r1FX ( P, T ) = r1X ( P, T ) + r1X ( P, )( + ) , F AC = c X F X − TScACF , X 10 (2.63) AC = AC f AC ( ) = X + c X X 2 XF F 2v AC 1 F − v v AC AC AC , 1+ (2.64) F 2 r01 + F X v AC F F 2r01F X + FX r r1X T 1X f AC ( ) = CAC AC , ( ) XF F r1X + T 2r01 X , T (2.72) (2.73) f AC ( F ) = f AC max = CAC AC max F f ACmax , C AC F (2.74) AC max f AC max = 1+ F C AC F (1 + F ) (2.75) ACmax = 1AC max = (2.71) F = AC − AC AC 1+ F f AC ( 0,2 ) (2.77) , , (2.78) e AC , + e (2.80) AC 0,2 0,2 ACe = AC (2.76) ( F ) lAC max = AC C AC = , e AC (2.81) + e 2.5 Kết tính số thảo luận kim loại hợp kim EYAC e = AC Bảng 2.7 EY(cC) (GPa) FeC P = T = 300 K tính PPTKMM theo TN Speich cộng (1972)[115] E c cC (%) PPTKMM TN [115] δ (%) 1,4 2,3 208,2 208,2 198,5 204,8 3,1 194,6 201,5 3,4 188,8 197,9 4,6 185,9 193,8 4,1 11 (1010Pa/%) -0,6309 -0,6180 2,1 Các kết tính số Fe, FeSi, FeH, FeC, Au, Cu, CuSi tổng kết Bảng 2.2 đến Bảng 2.10, Bảng 2.13 đến 2.20 minh họa Hình 2.2 đến Hình 2.19 7000 6800 6600 Vd (ms-1) 6400 6200 6000 5800 PPTKMM TN Shibazaki cộng (2016) TN Antonangeli Ohtani (2015) TN Decremps cộng (2014) TN Liu cộng (2014) 5600 5400 5200 5000 10 P (GPa) Hình 2.3 Vd ( P ) Fe T = 300K tính PPTKMM TN Antonangeli [11], Decremps [25], Liu [80] Shibazaki [111] Bảng 2.8 EY(cH) (GPa) FeH P = T = K tính PPTKMM ab initio Psiachos cộng (2011)[104] cH (%) PPTKMM ab initio [104] (%) 222,8 216,9 211,1 205,4 199,8 194,2 229,2 225,9 222,6 219,4 216,1 212,8 2,79 3,98 5,17 6,38 5,01 8,74 180 PPTKMM TN Santra cộng (2014) TN Ledbetter Naimon (1974) 160 EY (GPa) 140 120 100 80 cSi(%) Hình 2.11 EY (cSi ) CuSi T = 300K P = tính PPTKMM theo TN Ledbetter Naimon [75], Santra cộng [109] 12 18 1200 16 cSi = cSi = 2% cSi = 5% 14 800 1(MPa) f(GPa) 12 1000 10 600 400 Fe, PPTKMM Fe, TN tia X Smith cộng (2020) Fe, TN hình kĩ thuật số Smith cộng (2020) FeSi2%, PPTKMM FeSi5%, PPTKMM (b) 200 (a) 0 (%) (%) Hình 2.13 (a) f ( , cSi ) (b) ( , cSi ) FeSi T = 300 K, P = tính PPTKMM từ TN [113] Kết luận chương Bằng PPTKMM, xây dựng lý thuyết biến dạng đàn hồi phi tuyến HKXK hai thành phần với cấu trúc lập phương, rút biểu thức giải tích tổng quát lượng tự Helmholtz, độ dời hạt khỏi nút mạng, khoảng lân cận gần trung bình hai nguyên tử, đại lượng biến dạng đàn hồi môđun đàn hồi đẳng nhiệt đoạn nhiệt, môđun Young, môđun nén khối, môđun trượt, số đàn hồi, vận tốc sóng dọc, vận tốc sóng ngang, đại lượng biến dạng phi tuyến mật độ lượng biến dạng, ứng suất thực cực đại, giới hạn biến dạng đàn hồi để xác định đường cong ứng suất - độ biến dạng Kết giải tích đại lượng đặc trưng cho biến dạng thu phụ thuộc vào nhiệt độ, áp suất nồng độ nguyên tử xen kẽ Kết lý thuyết tính số, phân tích đánh giá kim loại Fe, Au, Cu HKXK FeSi, FeH, FeC, AuSi, CuSi phầm mềm Maple, Origin Từ kết tính tốn thu được, rút hai qui luật chung Thứ môđun đàn hồi EY, K, G, số đàn hồi C11, C12, C44, vận tốc sóng đàn hồi Vn, Vd, ứng suất thực cực đại 1max , giới hạn biến dạng đàn hồi e HKXK hai thành phần với cấu trúc lập phương giảm nồng độ nguyên tử xen kẽ tăng Thứ hai qui luật biến dạng phụ thuộc nhiệt độ, áp suất HKXK hai thành phần với cấu trúc lập phương tương tự 13 kim loại HKXK cấu trúc Các đại lượng biến dạng đàn hồi, vận tốc truyền sóng đàn hồi đại lượng biến dạng phi tuyến kim loại HKXK hai thành phần với cấu trúc lập phương giảm tăng nhiệt độ tăng tăng áp suất Nhiều kết tính số PPTKMM kim loại điều kiện nhiệt độ, áp suất khác có phù hợp tốt với số liệu thực nghiệm kết tính tốn khác Một số kết tính số môđun đàn hồi PPTKMM HKXK FeC, FeH, CuSi mô đun đàn hồi điều kiện nhiệt độ, nồng độ nguyên tử xen kẽ khác có phù hợp tốt với số liệu thực nghiệm kết tính tốn khác Các kết tính số chưa có số liệu so sánh kết nguồn tư liệu tham khảo để dự đoán, định hướng thực nghiệm tương lai CHƯƠNG 3: NGHIÊN CỨU BIẾN DẠNG ĐÀN HỒI VÀ PHI TUYẾN CỦA HỢP KIM XEN KẼ ABC VỚI CẤU TRÚC LẬP PHƯƠNG 3.1 Mơ hình hợp kim lượng tự Helmholtz Hình 3.1 Mơ hình HKXK ba thành phần ABC với cấu trúc (a) LPTK (b) LPTD Đối với mơ hình HKXK ba thành phần ABC cấu trúc LPTK với điều kiện nồng độ nguyên tử xen kẽ C nhỏ nhiều so với nồng độ nguyên tử thay B nồng độ nguyên tử thay B nhỏ nhiều so với nồng độ nguyên tử kim loại A, nguyên tử xen kẽ C nằm tâm mặt, nguyên tử kim loại A gọi A1 nằm tâm khối thay nguyên tử B nguyên tử kim loại A gọi A2 nằm đỉnh sở lập phương (xem Hình 3.1a) Đối với mơ hình HKXK ba thành phần ABC cấu trúc LPTD với 14 điều kiện nồng độ nguyên tử xen kẽ C nhỏ nhiều so với nồng độ nguyên tử thay B nồng độ nguyên tử thay B nhỏ nhiều so với nồng độ nguyên tử kim loại A, nguyên tử xen kẽ C nằm tâm khối, nguyên tử kim loại A gọi A1 nằm tâm mặt thay nguyên tử B nguyên tử kim loại A gọi A2 nằm đỉnh ô sở lập phương (xem Hình 3.2b) Xét hệ HKXK ABC lý tưởng với cấu trúc lập phương Số nguyên tử hợp kim N có NB nguyên tử thay B, NC nguyên tử xen kẽ C, N A1 nguyên tử A1, N A2 nguyên tử A2 NA = N – NB – NC – N A2 – N A1 nguyên tử A Gọi X lượng tự Helmholtz ứng với nguyên tử X (là nguyên tử A, B, A1, A2 C) Năng lượng tự Helmholtz hợp kim ABC hoàn toàn trật tự ABC = AC + cB ( B − A ) + TScAC − TScABC (3.1) AC = N c A A + cC C + c A1 A1 + c A2 A2 − TScAC ( ) lượng tự Helmholtz HKXK AC theo công thức (2.1), NX (X = A, B, C, A1, A2) số nguyên tử X, c X = NX nồng độ N nguyên tử X, T nhiệt độ, ScABC entrơpi cấu hình hợp kim ABC, ScAC tương ứng entrơpi cấu hình HKXK AC, cA = − 7cC , cA1 = 2cC , cA2 = 4cC cấu trúc LPTK cA = − 15cC , cA1 = 6cC , cA2 = 8cC cấu trúc LPTD 3.2 Khoảng lân cận gần trung bình a ABC = c AC a AC BTAC B + cB aB TB , BT = cAC BTAC + cB BTB BT BT 3.3 Biến dạng đàn hồi EYABC = EYAC + cB ( EYA − EYB ) , (3.2) (3.9) Xác định mơđun Young ta xác định môđun đàn hồi khác, số đàn hồi vận tốc truyền sóng đàn hồi theo cơng thức (3.10) đến (3.17) 3.4 Biến dạng phi tuyến (3.19) r1FX ( P, 0) = r1X ( P, 0)(1 + ) , 15 r1FX ( P, T ) = r1X ( P, T ) + r1X ( P, )( + ) , (3.20) F F ABC = AC + cB BF − AF + TScACF − TScABCF , (3.21) ( ) ABC = ABC f ABC ( ) = c X X + 2 F 2v ABC 2 F F X2 r1 X X 1 F − v v ABC ABC F 2r01 X T ( ) −cB A F − v v ABC ABC + 2 F 2v ABC 2 F F A2 r1A ABC , 1+ ( 2 r01F X + F v ABC F + FX r 1X ) XF F r1X AF F r1A F + FA r 1A + T 2r01 A , T f ABC ( ) = CABC ABC , f ABC ( F ) = f ABC max = CABC ABC max F (3.26) (3.27) (3.28) f ABCmax , C ABC F (3.29) ABC max f ABC max = 1+ F C ABC F (1 + F ) (3.30) ABCmax = 1ABC max = + 2r01 X − T F 2 r01 + F A v ABC F 2r01 A T (3.22) F = ABC , − ABC ( F ) lABC max = ABC C ABC = AC 1+ F f ABC ( 0,2 ) , (3.32) , (3.33) e ABC , + e (3.34) ABC 0,2 0,2 ABCe = ABC (3.31) EYABC e = ABC e ABC + e 16 (3.36) 240 280 220 260 200 240 180 K (GPa) 300 220 200 160 140 PPTKMM Ab initio Zhang 2010 Ab initio Olsson 2003 TN Heintze 2009 TN Speich 1972 180 (a) 160 PPTKMM Ab initio Zhang, 2010 Ab initio Olsson, 2003 TN Heintze, 2009 TN Specich, 1972 120 (b) 100 140 80 10 cCr (%) cCr (%) Hình 3.2 (a) EY ( cCr ) (b) K ( cCr ) FeCr T = 298K, P = tính PPTKMM, ab initio Olsson (2003)[100], ab initio Zhang (2010)[131] từ TN Heintze (2009)[40] Specich (1972)[115] 35 cSi = cSi = 1% cSi = 3% cSi = 5% 30 EY (1010 Pa) EY (GPa) 2.6 Kết tính số thảo luận hợp kim Các kết tính số cho FeCrSi, AuCuSi tóm tắt Bảng 3.1 đến Bảng 3.18 minh họa Hình 3.2 đến Hình 3.10 25 20 15 10 10 20 30 40 50 60 70 P (GPa) Hình 3.3 EY(P, cSi) AuCuSi cCu = 10% T = 300K 17 10 18 750 16 cSi=0 cSi=2% cSi=5% 14 700 f(GPa) f(GPa) 12 10 650 600 cSi=0 cSi=2% cSi=5% 550 (a) 2 (b) 500 7 (%) (%) Hình 3.4 (a) f ( , cSi ) (b) ( , cSi ) FeCrSi cCr =10%, T = 300K P = 25 800 P=0 P = 2,55 GPa P = 4,88 GPa P = 9,47 GPa P = 18,78 GPa 20 P=0 P = 2,55 GPa P = 4,88 GPa P = 9,47 GPa P = 18,78 GPa 700 600 500 f(GPa) 1(MPa) 15 10 400 300 200 100 (a) (b) 0 10 12 10 12 (%) (%) Hình 3.9 (a) f ( , P) (b) ( , P ) AuCuSi cCu =10%, cSi = 1% T = 300K Kết luận chương Bằng PPTKMM, xây dựng lý thuyết biến dạng đàn hồi phi tuyến HKXK ba thành phần với cấu trúc lập phương rút biểu thức giải tích tổng quát lượng tự Helmholtz, độ dời hạt khỏi nút mạng, khoảng lân cận gần trung bình hai nguyên tử, đại lượng biến dạng đàn hồi môđun đàn hồi đẳng nhiệt đoạn nhiệt, môđun Young, môđun nén khối, môđun trượt, số đàn hồi, vận tốc sóng dọc, vận tốc sóng ngang, đại lượng biến dạng phi tuyến mật độ lượng biến dạng, ứng suất thực cực đại, giới hạn biến dạng đàn hồi để xác định đường cong ứng suất - độ biến dạng Kết giải tích đại lượng đặc trưng cho 18 biến dạng thu phụ thuộc vào nhiệt độ, áp suất, nồng độ nguyên tử thay nồng độ nguyên tử xen kẽ Kết lý thuyết tính số, phân tích đánh giá HKXK FeCrSi AuCuSi phầm mềm Maple, Origin Từ kết tính tốn thu được, rút hai qui luật chung Thứ môđun đàn hồi EY, K, G, số đàn hồi C11, C12, C44, vận tốc sóng đàn hồi Vn, Vd , ứng suất thực cực đại 1max , giới hạn biến dạng đàn hồi e HKXK ba thành phần với cấu trúc lập phương giảm nồng độ nguyên tử xen kẽ tăng Qui luật tương tự qui luật HKXK hai thành phần với cấu trúc lập phương chương Thứ hai đại lượng biến dạng đàn hồi, vận tốc t sóng đàn hồi đại lượng biến dạng phi tuyến HKXK ba thành phần với cấu trúc lập phương giảm tăng nhiệt độ tăng tăng áp suất Qui luật tương tự qui luật HKXK hai thành phần qui luật kim loại A có cấu trúc Thứ ba đại lượng biến dạng đàn hồi, vận tốc truyền sóng đàn hồi đại lượng biến dạng phi tuyến HKXK ba thành phần với cấu trúc lập phương thay đổi chậm theo nồng độ nguyên tử thay Qui luật tương tự qui luật HKTT hai thành phần với cấu trúc Trong số trường hợp giới hạn nồng độ nguyên tử xen kẻ khơng, kết tính số PPTKMM so sánh với kết tính tốn khác số liệu thực nghiệm Đối với HKTT FeCr nồng độ Cr khác nhau, nhiều kết tính tốn môđun đàn hồi PPTKMM phù hợp tốt với số tính tốn ab initio số liệu thực nghiệm Các kết tính số chưa có số liệu so sánh kết nguồn tư liệu tham khảo để dự đoán, định hướng thực nghiệm tương lai CHƯƠNG 4: NGHIÊN CỨU ẢNH HƯỞNG CỦA BIẾN DẠNG LÊN SỰ KHUẾCH TÁN ĐỐI VỚI KIM LOẠI, HỢP KIM XEN KẼ AC VÀ ABC VỚI CẤU TRÚC LẬP PHƯƠNG 4.1 Ảnh hưởng độ biến dạng lên khuếch tán kim loại HKXK AC Trong mơ hình HKXK AC với cấu trúc lập phương chúng tơi, ngun tử xen kẽ nằm vị trí tâm mặt tâm khối mạng LPTK mạng LPTD Do đó, chúng tham gia tạo nút khuyết nguyên từ kim loại Từ điều kiện cực tiểu nhiệt động 19 Gibbs HKXK AC có khuyết tật (nút khuyết) suy nồng độ nút khuyết cân HKXK AC có khuyết tật c A gvf ( A) + c A gvf ( A1 ) + c A gvf ( A2 ) g f ( AC ) , ( 4.1) nv = exp − v nvA exp − k T k T Bo Bo Q(0, T ) DX ( ) = D0 X ( ) exp − (4.34) , kBoT ( AC V r + V m AC AC AC D0 X ( ) = D0 X (0, T ) exp kBoT Q(0, T ) DX ( P) = D0 X ( P) exp − , kBoT PV a D0 X ( P) = D0 X (0, T ) exp − AC k T Bo ) , (4.35) (4.41) (4.42) 4.2 Ảnh hưởng độ biến dạng lên khuếch tán HKXK ABC Ta có tất kết trình bày mục 4.1 Tuy nhiên với HKXK ABC ta có biểu thức c A gvf ( A) + c A gvf ( A1 ) + c A gvf ( A2 ) c g f ( B) exp − B v nv = exp − , kBoT kBoT RABC = N 1 − n n + n ( B v v X (4.43) − 1) c X X + nv n1c X X(1) − X ( ) − 1 − nv n1 + nv ( BA − 1) cB A + nv n1cB A(1) − T ScABC * + ScABC , (4.44) Đối với mạng LPTK V r (X ) = N ( ) (4.46) ( ) (4.47) R3 a ABC − a3ABC Đối với mạng LPTK V r (X ) = N R3 a ABC − a3ABC 4.3 Kết tính số thảo luận kim loại hợp kim Các kết tính số tổng kết bảng từ Bảng 4.1 đến Bảng 4.10 minh họa hình vẽ từ Hình 4.2 đến Hình 4.5 20 2000 Au PPTKMM Au TN Markin (1957) Au TN Herzig (1978) AuSi0,01 PPTKMM 1750 D(10-11 m2s-1) 1500 1250 1000 750 500 250 950 1000 1050 1100 1150 1200 1250 1300 1350 T (K) Hình 4.3 D(T)(D(cm2/s)) Au, AuSi0,01 P = tính PPTKMM theo TN Markin (1957)[83] Herzig (1978)[41] 0.30 0.050 Au AuSi0,01 AuSi0,03 0.045 0.040 0.28 D0(10-2 cm2/s) 0.035 D0(cm2/s) Cu CuSi0,01 CuSi0,03 0.29 0.030 0.025 0.020 0.015 0.27 0.26 0.25 0.24 0.010 (a) 0.23 0.005 0.000 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 0.22 1.0 (b) 1.5 2.0 2.5 3.0 3.5 4.0 4.5 (%) (%) Hình 4.4 D0(, cSi) (a) AuSi (b) CuSi T = 900K P = tỉnh PPTKMM Kết luận chương Trong chương 4, xây dựng biểu thức giải tích đại lượng đặc trưng cho khuếch tán lượng kích hoạt Q, hệ số trước hàm mũ D0, hệ số khuếch tán D phụ thuộc vào nhiệt độ áp suất Rút biểu thức giải tích cho phép xác định phụ thuộc ứng suất (hay độ biến dạng ) hệ số khuếch tán D hệ số trước hàm mũ D0 Theo kết tính số cho đại lượng đặc trưng cho khuếch tán nói kim loại, HKXK cấu trúc LPTD Au, Cu, AuSi CuSi, rút số qui luật chung Thứ giá trị lượng kích hoạt, hệ số khuếch tán hệ số trước mũ tăng nhiệt độ tăng Giá trị hệ số khuếch tán đáng kể vùng nhiệt độ cao Còn vùng nhiệt độ thấp, hệ số khuếch tán có giá trị nhỏ Thứ 21 5.0 hai hệ số khuếch tán hệ số trước hàm mũ giảm áp suất độ biến dạng (hay ứng suất ) tăng Thứ ba nguyên tử xen kẽ Si làm giảm hệ số khuếch tán, giảm hệ số trước hàm mũ tăng lượng kích hoạt Nhiều kết tính tốn phù hợp tốt với số liệu thực nghiệm, tính tốn CALPHAD ab initio Các kết tính số chưa có số liệu so sánh kết làm tư liệu tham khảo để định hướng, dự đoán thực nghiệm tương lai KẾT LUẬN CHUNG Trong luận án này, sử dụng PPTKMM để nghiên cứu trình biến dạng đàn hồi, biến dạng phi tuyến đồng thời nghiên cứu ảnh hưởng biến dạng lên trình khuếch tán kim loại, HKXK hai ba thành phần với cấu trúc lập phương Các kết luận án bao gồm: Rút biểu thức giải tích mơđun đàn hồi E, K, G, số đàn hồi C11, C12, C44 vận tốc sóng dọc Vd ngang Vn kim loại, HKXK hai ba thành phần với cấu trúc lập phương Thiết lập biểu thức giải tích mật độ lượng biến dạng, phụ thuộc ứng suất thực vào độ biến dạng, ứng suất thực cực đại, giới hạn biến dạng đàn hồi kim loại, HKXK hai ba thành phần với cấu trúc lập phương Xây dựng lý thuyết khuếch tán theo chế nút khuyết kim loại, HKXK hai ba thành phần với cấu trúc lập phương rút biểu thức giải tích đại lượng đặc trưng cho khuếch tán lượng kích hoạt, hệ số trước hàm mũ, hệ số khuếch tán thể tích kích hoạt phụ thuộc vào nhiệt độ áp suất Xây dựng biểu thức giải tích cho phép xác định phụ thuộc độ biến dạng hệ số khuếch tán hệ số trước hàm mũ Kết giải tích thu đại lượng đặc trưng cho biến dạng đàn hồi phi tuyến áp dụng tính số, thảo luận rút qui luật biến thiên theo nhiệt độ, áp suất nồng độ thành phần cho ba kim loại Fe, Au, Cu, ba HKXK hai thành phần FeSi, FeC, FeH, CuSi hai HKXK ba thành phần FeCrSi AuCuSi Áp dụng tính số, thảo luận rút qui luật biến thiên theo nhiệt độ, áp suất đại lượng đặc trưng cho khuếch tán kim loại Au, Cu hai HKXK hai thành phần CuSi, AuSi Chúng rút ba qui luật sau: 22 Thứ đại lượng biến dạng đàn hồi, vận tốc truyền sóng đàn hồi đại lượng biến dạng phi tuyến HKXK hai ba thành phần với cấu trúc lập phương giảm nồng độ nguyên tử xen kẽ tăng Thứ hai đại lượng biến dạng đàn hồi, vận tốc truyền sóng đàn hồi đại lượng biến dạng phi tuyến kim loại, HKXK hai ba thành phần với cấu trúc lập phương giảm tăng nhiệt độ, tăng tăng áp suất Thứ ba đại lượng biến dạng đàn hồi, vận tốc truyền sóng đàn hồi đại lượng biến dạng phi tuyến HKXK ba thành phần với cấu trúc lập phương thay đổi chậm theo nồng độ nguyên tử thay Các biểu thức thu PPTKMM luận án có dạng giải tích thuận lợi cho việc tính số, đánh giá số liệu phần mềm Maple, Origin, Hầu hết kết tính tốn PPTKMM luận án có sai số 10% so với thực nghiệm phù hợp với tính toán theo phương pháp khác ab initio, MD, CALPHAD, phù hợp với số liệu thực nghiệm Phần lớn kết tính số PPTKMM HKXK hai ba thành phần nhiệt độ, áp suất nồng độ thành phần khác mới, có tính dự báo, định hướng thực nghiệm tương lai HƯỚNG NGHIÊN CỨU TIẾP THEO Áp dụng kết lý thuyết luận án để tính số cho đại lượng đặc trưng cho biến dạng đàn hồi phi tuyến HKXK hai ba thành phần khác với hợp kim tính số luận án Xây dựng lý thuyết biến dạng đàn hồi phi tuyến HKXK hai thành phần với cấu trúc LGXC Xây dựng lý thuyết biến dạng đàn hồi phi tuyến HKXK hai ba thành phần xét đến khuyết tật điểm mạng tinh thể hợp kim DANH MỤC CƠNG TRÌNH ĐÃ CƠNG BỐ LIÊN QUAN ĐẾN LUẬN ÁN Hoc N Q and Hien N D (2018), Study on elastic deformation of substitution alloy AB with interstitrial atom C and BCC structure under pressure, 42nd Vietnam National Conference on Theoretical Physics 23 (NCTP-42), IOP Conf Series: Journal of Physics: Conf Series, 1034(2018)012005, IOP Publishing DOI: 10.1088/1742-6596/1034/1/012005 Hoc N Q, Tinh B D and Hien N D (2019), Elastic moduli and elastic constants of alloy AuCuSi with FCC structure under pressure, High Temperature Materials and Processes, 38, 264-272 DOI: 10.1515/htmp-2018-0027 Hoc N Q, Hoa N T and Hien N D (2019), Build the Theory of Nonlinear Deformation for BCC and FCC Substitutional Alloys AB with Interstitial Atom C under Pressure, HNUE Journal of Science, Natural Sciences, 64(6), 45-56 DOI: 10.18173/2354-1059.2019-0030 Nguyễn Quang Học, Nguyễn Đức Hiền, Nguyễn Thị Hòa, Lê Lan Hương Nguyễn Ngọc Lê (2022), Sự khuếch tán kim loại Au, Cu hợp kim xen kẽ AuSi, CuSi, HNUE Journal of Science, Natural Sciences, 67(1), 27-37 DOI: 10.18173/2354-1059.2022-0004 Nguyễn Quang Học, Nguyễn Đức Hiền, Nguyễn Thị Hòa, Phạm Phương Uyên Trịnh Hồng Ngọc (2022), Biến dạng đàn hồi vận tốc sóng đàn hồi hợp kim tam nguyên FeCrSi AuCuSi, HNUE Journal of Science, Natural Sciences, 67(1), 38-53 DOI: 10.18173/2354-1059.2022-0005 Hoc N Q and Hien N D (2022), Nonlinear deformation of Au, AuSi and AuCuSi with FCC structure under pressure, HNUE Journal of Science, Natural Sciences, 67(2), 40-54 DOI: 10.18173/2354-1059.2022-0021 Hoc N Q, Hien N D and Vi T K (2022), Elastic deformation and elastic wave velocity of iron and its binary interstitial alloys, Temperature, pressure and interstitial atom concentration dependences, Physica B: Condensed Matter, 644, 414134 DOI: https://doi.org/10.1016/j.physb.2022.414134 24