Đồ án tốt nghiệp mô phỏng quá trình hàn ma sát khuấy cho hợp kim nhôm bằng phần mềm abaqus

ĐẠI HỌC BÁCH KHOA HÀ NỘI TRƯỜNG CƠ KHÍ ĐỒ ÁN TỐT NGHIỆP Mơ q trình hàn ma sát khuấy cho hợp kim nhôm phần mềm ABAQUS HÀ NỘI, 3/2022 ĐẠI HỌC BÁCH KHOA HÀ NỘI TRƯỜNG CƠ KHÍ KHOA CƠ KHÍ CHẾ TẠO MÁY CỘNG HỒ XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập – Tự – Hạnh phúc ĐỒ ÁN TỐT NGHIỆP Đầu đề thiết kế Mơ q trình hàn ma sát khuấy cho hợp kim nhôm phần mềm ABAQUS Các số liệu ban đầu:  Thông số vật liệu chi tiết tấm: Hợp kim nhơm AA6061  Kích thước tấm: 30mm×7.5mm×3mm  Thơng số vật liệu dụng cụ: Thép dụng cụ H13  Kích thước dụng cụ: chiều dài chốt 1.6mm, đường kính chốt 2mm, đường kính vai 5.22mm  Tốc độ quay dụng cụ: 1500 rpm  Tốc độ chạy ngang dụng cụ: 150 mm/p – 200 mm/p  Hệ số ma sát: 0.58 – 0.8 Nội dung thuyết minh tính tốn:  Giới thiệu tổng quan phương pháp gia công hàn ma sát khuấy  Giới thiệu tổng quan phương pháp mơ số q trình hàn ma sát khuấy  Xây dựng mơ hình mơ q trình hàn ma sát khuấy phần mềm ABAQUS dựa phương pháp CEL kết hợp FEM  Kết mô thu bao gồm: Miền phân bố nhiệt độ mối hàn Biểu đồ biến thiên nhiệt độ theo miền thời gian Khảo sát biên dạng bề mặt mối hàn Khảo sát hình hành biên dạng khuyết tật Khảo sát tốc độ chuyển động vùng vật liệu tiếp xúc dụng cụ  Nghiên cứu ảnh hưởng thơng số đến q trình hàn ma sát khuấy: thơng số hình học dụng cụ, tốc độ chạy ngang dụng cụ hệ số ma sát Các vẽ  Bản vẽ mơ hình hàn ma sát khuấy hợp kim nhơm AA6061: Ao  Bản vẽ ảnh hưởng thông số đến kết mô số : Ao ĐỒ ÁN TỐT NGHIỆP Hà Nội, ngày tháng năm 2022 Giảng viên hướng dẫn Lời cảm ơn Trong q trình thực đồ án tơi giúp đỡ nhiệt tình thầy Phan Huy Lê Trong thực không khỏi mắc phải sai sót, lời nhận xét, góp ý bổ sung nhằm hoàn thiện đề tài thầy cô điều vô quý giá với Xin chân thành cảm ơn! Tóm tắt nội dung đồ án ĐỒ ÁN TỐT NGHIỆP Hàn ma sát khuấy công nghệ hàn trạng thái rắn Nó trở thành kỹ thuật độc đáo, thích hợp để nối nhiều kim loại khó hàn, đặc biệt hợp kim nhơm có độ bền cao số kim loại khác Do ưu điểm vượt trội so với kỹ thuật hàn nhiệt truyền thống, chẳng hạn khơng nóng chảy biến dạng thấp, hàn ma sát khuấy ứng dụng rộng rãi ngành công nghiệp hàng không vũ trụ, đóng tàu, tơ đường sắt Nghiên cứu kết hợp phân tích Lagrangian-Eulerian phần mềm thương mại ABAQUS để mơ tồn q trình hàn ma sát khuấy cho hợp kim nhơm AA6061 Mơ hình mơ số xây dựng để đánh giá trình truyền nhiệt, hình thành khuyết tật chi tiết hàn Ngoài miền phân bố nhiệt độ, khuyết tật biên dạng mối hàn có mối liên hệ chặt chẽ với thông số trình hàn Nghiên cứu mối quan hệ giúp xác định thông số tối ưu cho trình hàn, làm sở để nâng cao suất, chất lượng thiết kế chế tạo sản phẩm Sinh viên thực Ký ghi rõ họ tên MỤC LỤC DANH MỤC CÁC HÌNH VẼ -ĐỒ THỊ DANH MỤC CÁC BẢNG BIỂU CHƯƠNG GIỚI THIỆU TỔNG QUAN VỀ ĐỀ TÀI 10 1.1 Hàn trạng thái rắn 10 1.2 Một số phương pháp hàn ma sát điển hình 11 1.3 1.2.1 Hàn ma sát quay 11 1.2.2 Hàn ma sát tuyến tính 14 Hàn ma sát khuấy 16 ĐỒ ÁN TỐT NGHIỆP 1.3.1 Nguyên lý hàn ma sát khuấy 16 1.3.2 Vật liệu chiều dày chi tiết hàn 21 1.3.3 Thuật ngữ hàn ma sát khuấy 21 1.3.4 Vùng phân bố nhiệt độ trình hình thành mối hàn 22 1.3.5 Ưu điểm nhược điểm hàn ma sát khuấy 26 1.3.6 Ứng dụng hàn ma sát khuấy công nghiệp 27 CHƯƠNG GIỚI THIỆU TỔNG QUAN PHƯƠNG PHÁP MƠ PHỎNG SỐ CHO Q TRÌNH HÀN MA SÁT KHUẤY 34 2.1 Các bước chung cho q trình mơ số 34 2.2 Mức độ phân tích phương pháp tiếp cận 34 2.3 Mơ hình hình học lắp ghép chi tiết 34 2.4 Thuộc tính định luật cấu thành vật liệu 35 2.4.1 Mơ hình vật liệu Jonson-Cook 35 2.4.2 Mơ hình vật liệu Sheppard-Wright 36 2.5 Các mối quan hệ tương tác 37 2.6 Điều kiện biên 37 2.7 Chia lưới 37 2.8 Mô số 38 CHƯƠNG MƠ PHỎNG Q TRÌNH HÀN MA SÁT KHUẤY BẰNG PHẦN MỀM ABAQUS DỰA TRÊN PHƯƠNG PHÁP CEL 39 3.1 Mơ hình hình học 39 3.2 Mơ hình vật liệu 42 3.3 Chia lưới 44 3.4 Các điều kiện nhiệt độ 46 3.5 Điều kiện biên 52 3.6 Thiết lập tham số liên quan khác 53 CHƯƠNG NGHIÊN CỨU ẢNH HƯỞNG CỦA CÁC THÔNG SỐ ĐẾN QUÁ TRÌNH HÀN MA SÁT KHUẤY 56 4.1 Các đối tượng khảo sát trình hàn ma sát khuấy 56 4.1 Ảnh hưởng tốc độ dụng cụ 60 4.2 Ảnh hưởng biên dạng dụng cụ 64 4.3 Ảnh hưởng hệ số ma sát 71 CHƯƠNG KẾT LUẬN 75 TÀI LIỆU THAM KHẢO 76 ĐỒ ÁN TỐT NGHIỆP DANH MỤC CÁC HÌNH VẼ -ĐỒ THỊ Hình 1.1 Một số phương pháp hàn trạng thái rắn: a) Hàn ma sát b) Hàn siêu âm c) Hàn nổ 11 Hình 1.2 Nguyên lý phương pháp hàn ma sát quay [1] 12 Hình 1.3 Mối ghép chi tiết phương pháp hàn ma sát quay 13 Hình 1.4 Mơ hình hóa hàn ma sát quay phần mềm mơ số [2] 14 Hình 1.5 Nguyên lý phương pháp hàn ma sát tuyến tính [3] 14 Hình 1.6 Mối hàn chi tiết phương pháp hàn ma sát tuyến tính 15 Hình 1.7 Mơ hình hóa q trình hàn ma sát tuyến tính phần mềm ABAQUS [4] 16 ĐỒ ÁN TỐT NGHIỆP Hình 1.8 Sơ đồ minh họa trạng thái trình hàn ma sát khuấy [5] 17 Hình 1.9 Một số biên dạng dụng cụ hàn ma sát khuấy thực tế [6-7] 18 Hình 1.10 Hàn ma sát khuấy thực tế 19 Hình 1.11 Mơ hàn ma sát khuấy phần mềm ABAQUS [8-9] 20 Hình 1.12 Mối hàn hỗn hợp hai độ dày (1mm + mm) mối hàn hồng lên nhiều độ dày 21 Hình 1.13 Các thuật ngữ trình hàn ma sát khuấy [10] 22 Hình 1.14 Sơ đồ phân tán nhiệt bên mối ghép hàn ma sát khuấy [11-12] 23 Hình 1.15 Phân vùng vật liệu vùng hàn qua mơ số [13] 25 Hình 1.16 Ứng dụng hàn ma sát khuấy sản xuất khung xe tơ [14] 27 Hình 1.17 Ứng dụng hàn ma sát khuấy sản xuất vành xe tơ[14] 28 Hình 1.18 Ứng dụng hàn ma sát khuấy phương tiện hàng không 28 Hình 1.19 Ứng dụng hàn ma sát khuấy bể chứa bên Tàu thoi 29 Hình 1.20 Một hệ máy bay phản lực Tập đồn Hàng khơng Eclipse 29 Hình 1.21 Ứng dụng hàn ma sát khuấy cơng nghiệp đường sắt 30 Hình 1.22 Cấu trúc khung nhơm hai lớp đồn Class-395 30 Hình 1.23 Đồn tàu mơ Class-395 giai đoạn sản xuất trình hoàn thiện 31 Hình 1.24 Hàn ma sát khuấy thích hợp cho ván sàn, vách ngăn 31 Hình 1.25 Ứng dụng hàn ma sát khuấy ngành cơng nghiệp hàng hải 32 Hình 1.26 Ứng dụng hàn ma sát khuấy ngành cơng nghiệp hàng hải 32 Hình 1.27 Hàn ma sát khuấy sử dụng cho sản xuất tản nhiệt 33 Hình 2.1 Mơ hình lắp ghép phơi, dao, 35 Hình 3.1 Kích thước miền làm việc mơ số 39 Hình 3.2 Định nghĩa miền phôi dụng cụ Part Manager 40 Hình 3.3 Phân chia vùng trống (vùng đỏ) vùng gán vật liệu (vùng xanh) 40 Hình 3.4 Thiết lập miền gán vật liệu miền trống 41 Hình 3.5 Thơng số hình học dụng cụ 41 Hình 3.6 Thiết lập phần tử rắn đồng cho dụng cụ 42 Hình 3.7 Thiết lập thơng số vật liệu phôi AA6061 thép công cụ H13 44 Hình 3.8 Miền Eulerian chia thành 80640 phần tử EC3D8RT 44 Hình 3.9 Thiết lập phần tử EC3D8RT cho miền Eulerian 45 Hình 3.10 Miền dụng cụ chia thành 5109 phần tử C3D4T 45 Hình 3.11 Thiết lập phần tử C3D4T cho dụng cụ 46 Hình 3.12 Lắp ghép phơi dụng cụ chia lưới hoàn thành 46 Hình 3.13 Thiết lập hệ số tỏa nhiệt mặt chiều rộng phôi 48 ĐỒ ÁN TỐT NGHIỆP Hình 3.14 Thiết lập hệ số tỏa nhiệt mặt chiều dài phôi 48 Hình 3.15 Thiết lập hệ số tỏa nhiệt mặt phôi 49 Hình 3.16 Thiết lập hệ số tỏa nhiệt mặt đáy phôi 50 Hình 3.17 Thiết lập hệ số tỏa nhiệt mặt dụng cụ tiếp xúc 51 Hình 3.18 Thiết lập nhiệt độ ban đầu cho phơi 51 Hình 3.19 Thiết lập nhiệt độ ban đầu cho dụng cụ 52 Hình 3.20 Đặt điều kiện biên cho mơ hình 53 Hình 3.21 Thiết lập ba giai đoạn trình hàn 53 Hình 3.22 Thiết lập hệ số ma sát cho trình hàn ma sát khuấy 54 Hình 3.23 Thiết lập số luồng CPU cho trình mơ 55 Hình 4.1 Mơ hình khảo sát hàn ma sát khuấy 56 Hình 4.2 Ba giai đoạn trình hàn ma sát khuấy 58 Hình 4.3 Ảnh hưởng tốc độ dụng cụ đến phân bố nhiệt độ 59 Hình 4.4 Khảo sát nhiệt độ điểm xung quanh đường hàn theo miền thời gian với thay đổi tốc độ dụng cụ 61 Hình 4.5 So sánh biến dạng vùng hàn miền AS RS 62 Hình 4.6 Khuyết tật sinh bên vùng AS từ nghiên cứu nhóm tác giả P Chauhan [5] 62 Hình 4.7 Phát khuyết tật thông qua khảo sát vùng trống vật liệu 63 Hình 4.8 Khảo sát biên dạng khuyết tật trình FSW 64 Hình 4.9 Biên dạng dụng cụ số 65 Hình 4.10 Ảnh hưởng biên dạng dụng cụ đến phân bố nhiệt độ 67 Hình 4.11 Khảo sát nhiệt độ điểm xung quanh đường hàn theo miền thời gian với hai biên dạng dụng cụ khác 67 Hình 4.12 So sánh biên dạng mối hàn hai vùng AS RS sử dụng dụng cụ số 68 Hình 4.13 Phát khuyết tật thông qua khảo sát vùng trống vật liệu 68 Hình 4.14 Khảo sát biên dạng khuyết tật trình FSW 69 Hình 4.15 Khảo sát vận tốc chuyển động vật liệu xung quanh dụng cụ 71 Hình 4.16 Ảnh hưởng hệ số ma sát đến phân bố nhiệt độ 72 Hình 4.17 Khảo sát nhiệt độ điểm xung quanh đường hàn theo miền thời gian với thay đổi hệ số ma sát 73 Hình 4.18 Phát khuyết tật thơng qua khảo sát vùng trống vật liệu 74 ĐỒ ÁN TỐT NGHIỆP DANH MỤC CÁC BẢNG BIỂU Bảng Các số vật liệu Johnson-Cook cho loại vật liệu khác 36 Bảng Các số vật liệu Sheppard-Wright cho loại vật liệu khác 37 Bảng Giá trị hệ số ma sát cho loại vật liệu khác 37 Bảng Thơng số đặc tính vật liệu hợp kim nhôm AA6061 42 Bảng Hằng số vật liệu cho mơ hình JC vật liệu hợp kim nhôm AA6061 43 Bảng Thông số vật liệu thép dụng cụ H13 43 Bảng Hệ số tỏa nhiệt áp dụng cho cạnh chiều rộng chiều dài phôi 47 Bảng Hệ số tỏa nhiệt áp dụng cho mặt phôi 49 Bảng Hệ số tỏa nhiệt áp dụng cho mặt tiếp xúc dụng cụ 50 ĐỒ ÁN TỐT NGHIỆP CHƯƠNG GIỚI THIỆU TỔNG QUAN VỀ ĐỀ TÀI 1.1 Hàn trạng thái rắn Nhờ phát triển vượt bậc khoa học công nghệ, nhiều phương pháp hàn tiên tiến đời nhằm đáp ứng yêu cầu gia cơng kim loại Nhìn chung, có hai nhóm hàn phổ biến hàn nóng chảy hàn trạng thái rắn Mỗi loại có ưu nhược điểm riêng, dựa vào yêu cầu gia công chi tiết để lựa chọn phương pháp hàn phù hợp Trong quy trình hàn trạng thái rắn, mối hàn tạo nhiệt độ thấp nhiệt độ nóng chảy kim loại Một số phương pháp hàn trạng thái rắn điển hình hàn ma sát, hàn nổ, hàn rèn, hàn áp lực nóng hàn siêu âm Trong phương pháp này, có khuyết tật khơng xảy tượng nóng chảy đông đặc kim loại, đồng thời hạn chế khuyết tật Do đó, kim loại hàn nối giữ đặc tính ban đầu 10 ĐỒ ÁN TỐT NGHIỆP a) v=150mm/p b) v=200mm/p Hình 4.7 Phát khuyết tật thông qua khảo sát vùng trống vật liệu Ngồi nhờ cơng cụ chuyển đổi chế độ suốt cục (Toggle Global Translucency) ta quan sát hình thành khuyết tật cách trực quan thể Hình 4.8 Khuyết tật chạy dọc theo đường hàn, hình thành tồn giai đoạn trình hàn 63 ĐỒ ÁN TỐT NGHIỆP a) b) c) Hình 4.8 Khảo sát biên dạng khuyết tật trình FSW 4.2 Ảnh hưởng biên dạng dụng cụ Trong mục này, ta khảo sát ảnh hưởng biên dạng dụng cụ đến kết thu hình hàn ma sát khuấy Hai mơ hình dụng cụ sử dụng khảo sát Dụng cụ số có biên dạng giả thiết theo Mục 3.1 Hình 64 ĐỒ ÁN TỐT NGHIỆP 3.5 Dụng cụ số có biên dạng sau (Hình 4.9) với đường kính vai 5.22 mm, chốt hình trụ có cạnh chiều cao chốt 1.6mm Trong hai trường hợp khảo sát, thông số hàn giữ nguyên, thay đổi dụng cụ hàn Hình 4.9 Biên dạng dụng cụ số Hình 4.10 mơ tả phân bố nhiệt độ trình hàn ma sát khuấy với hai dụng cụ khác Ta quan sát hai trường hợp, vùng bề mặt tiếp xúc phôi vai dụng cụ nơi tập trung sinh nhiệt cao Do vùng tiếp xúc dụng cụ số rộng so với dụng cụ số nên vùng tập trung nhiệt độ cao rộng Ta quan sát từ mặt cắt mối hàn dụng cụ số 1, miền phân bố nhiệt độ không đối xứng Tuy nhiên mặt cắt mối hàn dụng cụ số 2, ta nhận thấy đối xứng tương đối miền phân bố nhiệt độ Điều khiến cho mối hàn sau trình hàn kết thúc đối xứng thể Hình 4.10b Kết khảo sát chứng tỏ dụng cụ số có đặc tính ưu việt so với dụng cụ số Kết thu từ Hình 4.11 cho thấy, nhiệt độ cực đại điểm A, B C sử dụng dụng cụ số cho giá trị thấp so với sử dụng dụng cụ số Bên cạnh thứ tự nhiệt độ cực đại giảm dần từ A đến C Bên cạnh Hình 4.12 thể biên dạng mối hàn hai miền AS RS gần đối xứng lỗi mối hàn không hình thành khuyết tật sử dụng dụng cụ số Thật vậy, theo kết thu từ mặt cắt ngang mối hàn (như mơ tả Hình 4.13), ta không phát vùng trống vật liệu (vùng đỏ) xuất bên đường hàn Điều chứng tỏ với biên dạng dụng cụ số ta thu mối hàn có chất lượng tốt so với mối hàn số Để nhấn mạnh thêm khẳng định này, ta quan sát cách trực quan thể Hình 4.14 nhờ cơng cụ chuyển đổi chế độ suốt cục (Toggle Global Translucency) 65 ĐỒ ÁN TỐT NGHIỆP a) Dụng cụ số 66 ĐỒ ÁN TỐT NGHIỆP b) Dụng cụ số Hình 4.10 Ảnh hưởng biên dạng dụng cụ đến phân bố nhiệt độ a) Dụng cụ số b) Dụng cụ số Hình 4.11 Khảo sát nhiệt độ điểm xung quanh đường hàn theo miền thời gian với hai biên dạng dụng cụ khác 67 ĐỒ ÁN TỐT NGHIỆP Hình 4.12 So sánh biên dạng mối hàn hai vùng AS RS sử dụng dụng cụ số a) Dụng cụ số b) Dụng cụ số Hình 4.13 Phát khuyết tật thơng qua khảo sát vùng trống vật liệu 68 ĐỒ ÁN TỐT NGHIỆP a) b) c) Hình 4.14 Khảo sát biên dạng khuyết tật q trình FSW Sự khơng hình thành khuyết tật sử dụng dụng cụ số giải thích thơng qua khảo sát véc tơ vận tốc vật liệu (như mơ tả Hình 4.15) Ta dễ dàng quan sát hai giai đoạn (gia cố nhiệt độ) giai đoạn (thực đường hàn), trường hợp sử dụng dụng cụ số 2, vật liệu vùng tiếp xúc dụng cụ chuyển động dễ dàng so với sử dụng dụng cụ số Điều góp phần 69 ĐỒ ÁN TỐT NGHIỆP hạn chế khuyết tật xảy thực trình hàn ma sát khuấy Nguyên nhân chuyển động ma sát phôi bề mặt dụng cụ Với biên dạng dụng cụ khác gây ảnh hưởng đến vận tốc chuyển động vật liệu Kết thu chứng tỏ với biên dạng dụng cụ số đem lại ưu điểm so với dụng cụ có biên dạng số a) Dụng cụ số 1: Giai đoạn b) Dụng cụ số 2: Giai đoạn 70 ĐỒ ÁN TỐT NGHIỆP c) Dụng cụ số 1: Giai đoạn d) Dụng cụ số 2: Giai đoạn Hình 4.15 Khảo sát vận tốc chuyển động vật liệu xung quanh dụng cụ 4.3 Ảnh hưởng hệ số ma sát Hình 4.16 thể phân bố nhiệt độ trình hàn ma sát khuấy hai chế độ hệ số ma sát khác m=0.58 m=0.8 Trong hai trường hợp, vùng tiếp xúc phôi bề mặt vai dụng cụ tập trung nhiệt cao hai miền phân bố nhiệt độ không đối xứng Điều gây biến dạng không sau trình hàn kết thúc 71 ĐỒ ÁN TỐT NGHIỆP a) m=0.8 b) m=0.58 Hình 4.16 Ảnh hưởng hệ số ma sát đến phân bố nhiệt độ 72 ĐỒ ÁN TỐT NGHIỆP Để khảo sát chi tiết thêm ảnh hưởng hệ số ma sát đến truyển nhiệt trình hàn ma sát khuấy, ta kiểm tra biến thiên nhiệt độ theo miền thời gian Kết thu từ Hình 4.17 cho thấy hệ số ma sát giảm từ m=0.8 xuống m=0.58, nhiệt độ cực đại điểm A, B C giảm theo Thứ tự giá trị cực đại điểm A, B C khơng đổi Bên cạnh Hình 4.18 khảo sát khuyết tật thơng qua phát miền trống vật liệu Khi hệ số ma sát giảm cịn m=0.58 khuyết tật hình thành dọc theo đường hàn a) m=0.8 b) m=0.58 Hình 4.17 Khảo sát nhiệt độ điểm xung quanh đường hàn theo miền thời gian với thay đổi hệ số ma sát 73 ĐỒ ÁN TỐT NGHIỆP a) m=0.8 b) m=0.58 Hình 4.18 Phát khuyết tật thông qua khảo sát vùng trống vật liệu 74 ĐỒ ÁN TỐT NGHIỆP CHƯƠNG KẾT LUẬN Mơ hình phần tử hữu hạn nhiệt phát triển dựa phương pháp Coupled Eulerian Lagrangian để mơ quy trình hàn ma sát khuấy hai hợp kim nhơm AA6061 Mơ hình thực Abaqus/Explicit Kết đầu thu nghiên cứu bao gồm nội dung sau:      Khảo sát miền phân bố nhiệt độ qua mặt cắt ngang mối hàn Biểu đồ nhiệt độ biến thiên theo miền thời gian Khảo sát hình thành khuyết tật Khảo sát bề mặt mối hàn sau trình hàn kết thúc Vận tốc chuyển động phần tử tiếp xúc dụng cụ mối hàn Ảnh hưởng tốc độ chạy dụng cụ, biên dạng dụng cụ hệ số ma sát khuấy nghiên cứu rút kết luận sau:  Nhiệt độ tập trung cao lõi hàn tỏa dần miền xung quanh theo thời gian  Vận tốc hàn tăng từ 150mm/p lên 200mm/p nhiệt độ đạt cực đại điểm khảo sát A giảm từ 475oC xuống 410oC  Khi xử dụng dụng cụ số ta nhận thấy đối xứng tương đối miền phân bố nhiệt độ, điều trái ngược so với sử dụng dụng cụ số Khi cho mối hàn đẹp  Khi sử dụng dụng cụ số 2, nhiệt độ cực đại điểm khảo sát giảm khơng hình thành khuyết tật  Hệ số ma sát giảm từ 0.8 xuống 0.58 nhiệt độ cực đại điểm khảo sát giảm, hình thành khuyết tật dọc đường hàn 75 ĐỒ ÁN TỐT NGHIỆP TÀI LIỆU THAM KHẢO [1] García, A M M (2011) BLISK fabrication by linear friction welding Advances in gas turbine technology, 411-434 [2] Seli, H., Awang, M., Ismail, A I M., Rachman, E., & Ahmad, Z A (2013) Evaluation of properties and FEM model of the friction welded mild steel-Al6061alumina Materials Research, 16(2), 453-467 [3] Jedrasiak, P., Shercliff, H R., McAndrew, A R., & Colegrove, P A (2018) Thermal modelling of linear friction welding Materials & Design, 156, 362-369 [4] Okeke, S I., Harrison, N., & Tong, M (2020) Thermomechanical modelling for the linear friction welding process of Ni-based superalloy and verification Proceedings of the Institution of Mechanical Engineers, Part L: Journal of Materials: Design and Applications, 234(5), 796-815 [5] Chauhan, P., Jain, R., Pal, S K., & Singh, S B (2018) Modeling of defects in friction stir welding using coupled Eulerian and Lagrangian method Journal of Manufacturing Processes, 34, 158-166 [6] Rao, C V., Reddy, G M., & Rao, K S (2015) Influence of tool pin profile on microstructure and corrosion behaviour of AA2219 Al–Cu alloy friction stir weld nuggets Defence technology, 11(3), 197-208 [7] Shamsudeen, S., & Dhas, J E R (2018) Optimization of multiple performance characteristics of friction stir welded joint with grey relational analysis Materials Research, 21 [8] Al-Badour, F., Merah, N., Shuaib, A., & Bazoune, A (2013) Coupled Eulerian Lagrangian finite element modeling of friction stir welding processes Journal of Materials Processing Technology, 213(8), 1433-1439 [9] Al-Badour, F., Merah, N., Shuaib, A., & Bazoune, A (2014) Thermomechanical finite element model of friction stir welding of dissimilar alloys The International Journal of Advanced Manufacturing Technology, 72(5), 607-617 [10] Jain, R., Pal, S K., & Singh, S B (2017) Numerical modeling methodologies for friction stir welding process In Computational methods and production engineering (pp 125-169) Woodhead Publishing [11] Grujicic, M., Arakere, G., Yalavarthy, H V., He, T., Yen, C F., & Cheeseman, B A (2010) Modeling of AA5083 material-microstructure evolution during butt friction-stir welding Journal of Materials Engineering and Performance, 19(5), 672-684 [12] Aldanondo, E., Vivas, J., Álvarez, P., & Hurtado, I (2020) Effect of tool geometry and welding parameters on friction stir welded lap joint formation with AA2099-T83 and AA2060-T8E30 aluminium alloys Metals, 10(7), 872 [13] Assidi, M., Fourment, L., Guerdoux, S., & Nelson, T (2010) Friction model for friction stir welding process simulation: Calibrations from welding 76 ĐỒ ÁN TỐT NGHIỆP experiments International Journal of Machine Tools and Manufacture, 50(2), 143155 [14] Besharati-Givi, M K., & Asadi, P (2014) Advances in friction-stir welding and processing Elsevier [15] Kuykendall, K., Nelson, T., & Sorensen, C (2013) On the selection of constitutive laws used in modeling friction stir welding International Journal of Machine Tools and Manufacture, 74, 74-85 [16] Sheppard, T., & Jackson, A (1997) Constitutive equations for use in prediction of flow stress during extrusion of aluminium alloys Materials science and Technology, 13(3), 203-209 [17] Uyyuru, R K., & Kailas, S V (2006) Numerical analysis of friction stir welding process Journal of Materials Engineering and Performance, 15(5), 505518 [18] Tang, J., & Shen, Y (2016) Numerical simulation and experimental investigation of friction stir lap welding between aluminum alloys AA2024 and AA7075 Journal of Alloys and Compounds, 666, 493-500 [19] Grujicic, M., Arakere, G., Pandurangan, B., Ochterbeck, J M., Yen, C F., Cheeseman, B A., & Sutton, M A (2012) Computational analysis of material flow during friction stir welding of AA5059 aluminum alloys Journal of materials engineering and performance, 21(9), 1824-1840 [20] Salloomi, K N., & Al-Sumaidae, S (2021) Coupled Eulerian–Lagrangian prediction of thermal and residual stress environments in dissimilar friction stir welding of aluminum alloys Journal of Advanced Joining Processes, 3, 100052 [21] Meyghani, B., Awang, M B., Momeni, M., & Rynkovskaya, M (2019, April) Development of a finite element model for thermal analysis of friction stir welding (FSW) In IOP conference series: materials science and engineering (Vol 495, No 1, p 012101) IOP Publishing [22] Schmidt, H., & Hattel, J (2004) A local model for the thermomechanical conditions in friction stir welding Modelling and simulation in materials science and engineering, 13(1), 77 [23] Al-Badour, F A K (2012) Numerical and experimental investigations of friction stir welding of tube-tubesheet joints (Doctoral dissertation, King Fahd University of Petroleum and Minerals (Saudi Arabia)) [24] Javadi, M., & Tajdari, M (2006) Experimental investigation of the friction coefficient between aluminium and steel MATERIALS SCIENCE-WROCLAW, 24(2/1), 305 [25] Hou, Z., Sheikh-Ahmad, J., Jarrar, F., & Ozturk, F (2018) Residual stresses in dissimilar friction stir welding of AA2024 and AZ31: experimental and numerical study Journal of Manufacturing Science and Engineering, 140(5) 77