Tổng quan nghiên cứu

Hợp kim nhôm A5083 thuộc hệ nhôm - magiê đóng vai trò then chốt trong ngành công nghiệp đóng tàu cao tốc, hàng không và kết cấu quân sự nhờ trọng lượng riêng chỉ khoảng 2,7 g/cm³ cùng độ bền kéo đạt 275 đến 345 MPa. Tuy nhiên, quá trình gia công hàn nhôm gặp phải rào cản kỹ thuật lớn do lớp màng oxit Al2O3 có nhiệt độ nóng chảy lên tới 2050°C, vượt xa nhiệt độ nóng chảy của nhôm nguyên chất là 660°C. Đồng thời, hệ số dẫn nhiệt cao cùng hệ số giãn nở nhiệt lớn khiến liên kết hàn góc chữ T rất dễ bị cong vênh nghiêm trọng và tích tụ ứng suất dư nhiệt cục bộ.

Vấn đề nghiên cứu trọng tâm đặt ra là việc xác định chính xác trường nhiệt độ biến thiên, mức độ tập trung ứng suất và biên độ biến dạng dư trong liên kết hàn góc chữ T hợp kim nhôm A5083 có chiều dày 8 mm khi thực hiện bằng phương pháp hàn MIG (GMAW). Mục tiêu cụ thể của luận văn là xây dựng mô hình mô phỏng số 3D bằng phương pháp phần tử hữu hạn, giải quyết bài toán cơ nhiệt phi tuyến, từ đó xác lập bộ thông số công nghệ hàn tối ưu nhằm kiểm soát chất lượng mối hàn.

Nghiên cứu được triển khai thực nghiệm và tính toán mô phỏng trong năm 2021 tại Đại học Sư phạm Kỹ thuật Nam Định, ứng dụng trực tiếp kinh nghiệm thực tế từ các dự án chế tạo cabin tàu tuần tra tại các nhà máy đóng tàu ven biển. Kết quả của đề tài mang ý nghĩa thực tiễn to lớn khi giúp các đơn vị sản xuất cắt giảm 15% đến 25% chi phí gia công nắn chỉnh sau hàn, hạ tỷ lệ khuyết tật nứt nhiệt và rỗ khí xuống dưới 5%, đồng thời nâng cao tuổi thọ làm việc của kết cấu vỏ tàu biển.

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 dựa trên nền tảng lý thuyết truyền nhiệt phi tuyến trong vật rắn của Joseph Fourier và phương trình vi phân dẫn nhiệt ba chiều ở trạng thái giả ổn định của Rosenthal - Rykalin. Mô hình nguồn nhiệt thể tích phân bố Gauss bán ellipsoid kép của Goldak được lựa chọn để mô tả chính xác đặc tính cấp nhiệt của hồ quang hàn MIG với 4 thông số hình học: nửa chiều rộng b = 5 mm, chiều sâu ngấu c = 4,5 mm, chiều dài phần trước af = 5 mm và chiều dài phần sau ar = 7 mm.

Bên cạnh đó, khung phân tích nhiệt - đàn hồi - dẻo (Thermo-Elasto-Plastic) được thiết lập dựa trên tiêu chuẩn chảy dẻo Von Mises kết hợp với quy tắc hóa bền đẳng hướng để đánh giá ứng xử cơ học của vật liệu khi thay đổi nhiệt độ. Động học chuyển biến pha và biến dạng nhiệt luyện kim được tích hợp theo mô hình Leblond - Johnson - Mehl - Avrami. Các khái niệm cốt lõi bao gồm: năng lượng đường hàn, vùng ảnh hưởng nhiệt (HAZ), ứng suất dư nhiệt, biến dạng co ngót thể tích và kỹ thuật tiêu diệt - tái sinh phần tử (Element Death and Birth) để mô phỏng sự bổ sung kim loại đắp theo thời gian thực.

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

Nghiên cứu sử dụng cỡ mẫu thực nghiệm gồm các mẫu liên kết hàn góc chữ T từ phôi nhôm A5083 với kích thước tấm đế và tấm đứng dài 300 mm, rộng 120 mm và dày 8 mm, được chuẩn bị vát mép góc 60° ± 5° với khe hở lắp ghép không quá 2 mm. Phương pháp chọn mẫu là chọn mẫu kỹ thuật định chuẩn theo tiêu chuẩn hàn kết cấu nhôm AWS D1.2M, đại diện cho kết cấu vách ngăn và cabin tàu thủy thực tế.

Lý do lựa chọn phương pháp phân tích số phần tử hữu hạn (FEM) trên phần mềm chuyên dụng SYSWELD (kết hợp VisualMesh và VisualWeld) là khả năng giải quyết bài toán ghép đôi nhiệt - cơ học phi tuyến phức tạp với độ chính xác cao. Phương pháp này cho phép cập nhật liên tục các thông số nhiệt - vật lý phụ thuộc nhiệt độ như nhiệt dung riêng, hệ số dẫn nhiệt và mô đun Young của hợp kim nhôm A5083, điều mà các phương pháp đo đạc thực nghiệm truyền thống không thể theo dõi liên tục trong toàn bộ không gian khối.

Mô hình hình học 3D được phân lưới với 28.000 nút và khoảng 40.000 phần tử khối (Linear Hex 8 và Linear Penta 6). Kích thước lưới vùng mối hàn và vùng HAZ được chia mịn từ 1 mm đến 3 mm, bước thời gian giải được tối ưu hóa ở mức 0,44s đến 0,45s mỗi bước. Quá trình tính toán mô phỏng diễn ra trong tổng thời gian 1000s, bao gồm giai đoạn hàn và làm nguội tự nhiên trong môi trường không khí ở nhiệt độ ban đầu 20°C.

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

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

Mô phỏng trường nhiệt độ cho thấy nhiệt độ đỉnh cực đại tại tâm vũng hàn đạt 1174°C ở nút 12211, đảm bảo kim loại dây hàn ER5183 và kim loại cơ bản nóng chảy hoàn toàn. Sau 1000s làm nguội, nhiệt độ toàn bộ kết cấu hạ xuống mức an toàn là 76°C. Tại các vị trí lân cận mối hàn, nhiệt độ ghi nhận sự phân tầng rõ rệt: nút 4786 (cách chân mối hàn 2 mm) đạt nhiệt độ cực đại 585°C, trong khi nút 18536 (cách chân mối hàn 6 mm) đạt mức nhiệt 475°C.

Về biến dạng hình học, độ cong vênh xuất hiện lớn nhất tại mép tự do không bị gá kẹp, đạt đỉnh ở vị trí giữa chiều dài tấm hàn (tại khoảng cách 150 mm). Giá trị độ võng biến dạng theo mô phỏng là 3,12 mm, rất sát với giá trị đo đạc thực nghiệm là 3,30 mm, mức sai số tương đối chỉ 5,45%.

Phân bố ứng suất tương đương Von Mises tập trung cao nhất tại các khu vực bị ngàm kẹp cố định, đạt 212,2 MPa đối với phương án vát mép một phía và 205,0 MPa đối với phương án vát mép hai phía. Ngược lại, tại các vùng biên không gá kẹp, ứng suất dư giảm mạnh xuống dưới 100 MPa.

Đối với ứng suất pháp theo phương X (vuông góc đường hàn), ứng suất kéo cực đại ở trường hợp vát mép một phía đạt 231,27 MPa, cao hơn trường hợp vát mép hai phía đạt 203,78 MPa (chênh lệch 27,49 MPa, tương đương giảm khoảng 11,9%). Tại hai đầu điểm mồi và kết thúc hồ quang hàn, xuất hiện vùng ứng suất nén cục bộ với biên độ xấp xỉ 130 MPa. Ứng suất pháp kéo dọc trục mối hàn theo phương Y đạt đỉnh 143,8 MPa tại nút 20979, trong khi ứng suất nén đạt 66,12 MPa tại nút 27248.

Thảo luận kết quả

Hiện tượng chênh lệch nhiệt độ lớn giữa vùng tâm mối hàn (1174°C) và vùng cách chân mối 6 mm (475°C) tạo ra một gradient nhiệt độ cực dốc, nguyên nhân trực tiếp bắt nguồn từ hệ số dẫn nhiệt rất cao của nhôm. Sự co ngót không đồng đều giữa kim loại nóng chảy trong vũng hàn và kim loại cơ bản nguội nhanh hơn đã sinh ra nội ứng suất kéo mạnh tại vùng đường hàn và ứng suất nén ở vùng lân cận.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua đồ thị đường cong biến dạng dọc theo 7 mốc khoảng cách từ 0 mm đến 300 mm (0 mm: 1,9 mm thực nghiệm so với 2,03 mm mô phỏng; 150 mm: 3,3 mm so với 3,12 mm; 300 mm: 1,5 mm so với 2,02 mm). Bảng đối chiếu thực nghiệm và mô phỏng khẳng định độ tin cậy vượt trội của mô hình FEM khi toàn bộ sai số nằm trong giới hạn cho phép của kỹ thuật hàn kết cấu. Đồng thời, biểu đồ phân bố ứng suất Von Mises và ứng suất pháp 3 chiều được trực quan hóa qua thang màu 3D từ xanh dương (ứng suất thấp) đến đỏ đậm (ứng suất cao), giúp kỹ sư nhận diện chính xác các vị trí tập trung ứng suất nguy hiểm.

Khi so sánh với các nghiên cứu của Preston về mô phỏng hàn TIG hợp kim nhôm A2024 và nghiên cứu của Dudzik trên hợp kim A5083, phương pháp hàn MIG trong nghiên cứu này thể hiện ưu thế vượt trội về tốc độ cấp nhiệt và khả năng khống chế chiều sâu ngấu đồng đều. Kết quả này cũng hoàn toàn tương thích với các quy luật lý thuyết về cơ chế hình thành ứng suất dư do co ngót nhiệt của Macherauch và Wohlfahrt.

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

Thứ nhất, chuẩn hóa chế độ thông số công nghệ hàn MIG cho kết cấu hợp kim nhôm A5083 dày 8 mm trong sản xuất thực tế. Đặt năng lượng đường hàn ở mức 750 J/mm cho lớp hàn lót và 800 J/mm cho lớp hàn phủ, duy trì vận tốc hàn ổn định 5 mm/s (tương đương 300 mm/phút) với khí bảo vệ Argon nguyên chất 99,98%. Mục tiêu là kiểm soát độ võng biến dạng dưới 3,0 mm và đảm bảo chiều sâu ngấu đạt 100% tiết diện chân mối ghép. Thời gian chuẩn hóa quy trình dự kiến hoàn thành trong vòng 3 tháng đầu tiên, do đội ngũ kỹ sư công nghệ hàn tại các xưởng gia công cơ khí chủ trì thực hiện.

Thứ hai, thiết kế và tối ưu hóa hệ thống đồ gá kẹp phôi đối xứng. Ứng dụng giải pháp gá kẹp có tấm lót dẫn nhiệt kiểm soát tại vùng đáy liên kết chữ T nhằm giảm tốc độ nung nóng không đều, hướng tới mục tiêu giảm 15% đến 20% ứng suất kéo cục bộ tại các vị trí ngàm. Thời gian triển khai chế tạo và thử nghiệm đồ gá là 6 tháng, do phòng thiết kế công nghệ và xưởng chế tạo đồ gá phối hợp thực hiện.

Thứ ba, áp dụng công nghệ xử lý khử ứng suất dư sau hàn bằng phương pháp rung siêu âm hoặc gõ búa giảm ứng suất cục bộ tại các vị trí có mức ứng suất vượt trên 200 MPa (đặc biệt là vùng lân cận các điểm mồi hồ quang và nút 20979). Mục tiêu là triệt tiêu từ 30% đến 40% đỉnh ứng suất dư kéo, ngăn ngừa nguy cơ nứt tế vi khi tàu vận hành trong môi trường sóng biển tải trọng động. Thời gian áp dụng quy trình kiểm soát này là 9 tháng, do bộ phận kiểm định chất lượng (QA/QC) và kỹ thuật viên xưởng đóng tàu đảm nhận.

Thứ tư, tích hợp phần mềm mô phỏng SYSWELD vào quy trình thiết kế tiền sản xuất cho các kết cấu tấm dày từ 4 mm đến 12 mm. Mục tiêu nhằm cắt giảm 40% thời gian thử nghiệm hàn mẫu thực tế và tiết kiệm 30% chi phí phôi liệu nhôm phục vụ kiểm tra phá hủy. Kế hoạch đào tạo nhân lực và số hóa quy trình công nghệ được đề xuất thực hiện trong lộ trình 12 tháng, do ban giám đốc kỹ thuật của doanh nghiệp chỉ đạo phối hợp cùng các viện nghiên cứu chuyên ngành.

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

Nhóm thứ nhất là các kỹ sư công nghệ hàn và kỹ sư kết cấu tại các nhà máy đóng tàu vỏ nhôm. Nhóm này có thể ứng dụng trực tiếp bộ thông số năng lượng đường (750 đến 800 J/mm), vận tốc hàn (5 mm/s) và sơ đồ vát mép hai phía để tối ưu hóa quy trình chế tạo cabin, vách ngăn và khung vỏ tàu tuần tra, giúp tiết kiệm chi phí thử nghiệm thực tế.

Nhóm thứ hai là các nhà nghiên cứu, nghiên cứu sinh và giảng viên chuyên ngành cơ khí chế tạo máy, kỹ thuật vật liệu và cơ học tính toán. Luận văn cung cấp hệ thống dữ liệu nhiệt vật lý chi tiết của hợp kim nhôm A5083 cùng mô hình nguồn nhiệt Goldak kép, làm tài liệu tham khảo giá trị để phát triển các bài toán mô phỏng cơ nhiệt phi tuyến phức tạp.

Nhóm thứ ba là cán bộ kiểm định chất lượng (QA/QC), giám định viên đăng kiểm hàng hải và các chuyên viên đánh giá an toàn kết cấu. Các phân tích định lượng về trường ứng suất Von Mises và ứng suất ba phương (X, Y, Z) giúp định vị chính xác những vùng tập trung ứng suất nguy hiểm tiềm ẩn, từ đó xây dựng quy trình kiểm tra không phá hủy (NDT) hiệu quả hơn.

Nhóm thứ tư là học viên cao học và sinh viên ngành kỹ thuật cơ khí, công nghệ hàn. Đây là tài liệu mẫu mực về phương pháp kết hợp nhuần nhuyễn giữa mô phỏng số bằng phần mềm công nghiệp hiện đại và kiểm chứng thực nghiệm bằng mô hình vật lý, hỗ trợ đắc lực cho quá trình học tập và làm đồ án tốt nghiệp.

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

Hợp kim nhôm A5083 có những đặc tính nhiệt nào gây khó khăn lớn nhất cho quá trình hàn? Hợp kim A5083 có lớp màng oxit Al2O3 nóng chảy ở 2050°C trong khi nhôm nóng chảy ở 660°C, đòi hỏi phải dùng dòng điện cực nghịch để phá màng oxit. Đồng thời, hệ số dẫn nhiệt và hệ số giãn nở nhiệt cao làm vũng hàn nguội nhanh, dễ bẫy khí hydro gây rỗ xốp và làm phát sinh biến dạng cong vênh lớn.

Mô hình nguồn nhiệt thể tích Goldak kép mang lại ưu thế gì so với nguồn nhiệt điểm truyền thống? Mô hình Goldak kép mô phỏng nguồn nhiệt dạng khối elip với hai nửa bán cầu độc lập (trước và sau tâm hồ quang), phản ánh chính xác sự phân bố mật độ năng lượng thực tế của hồ quang hàn MIG. Mô hình này khắc phục hoàn toàn nhược điểm cho ra nhiệt độ vô hạn tại tâm của nguồn nhiệt điểm Rosenthal.

Sai số biến dạng giữa mô phỏng số trên SYSWELD và đo đạc thực nghiệm là bao nhiêu? Độ lệch biến dạng uốn lớn nhất tại điểm giữa mép tấm (khoảng cách 150 mm) giữa mô phỏng đạt 3,12 mm và thực nghiệm đạt 3,30 mm chỉ là 0,18 mm, tương đương sai số tương đối 5,45%. Mức sai số rất nhỏ này khẳng định độ tin cậy cao của mô hình phần tử hữu hạn.

Năng lượng đường hàn tối ưu được xác lập trong nghiên cứu là bao nhiêu? Năng lượng đường tối ưu cho liên kết hàn góc chữ T nhôm A5083 dày 8 mm được xác lập là 750 J/mm cho lớp hàn lót và 800 J/mm cho lớp hàn phủ, kết hợp với vận tốc hàn ổn định 5 mm/s dưới sự bảo vệ của khí Argon có độ tinh khiết 99,98%.

Quy cách vát mép hai phía mang lại lợi ích gì vượt trội so với vát mép một phía? Phương án vát mép hai phía giúp phân bố dòng nhiệt và kim loại đắp cân bằng hơn, làm giảm ứng suất tương đương Von Mises cực đại từ 212,2 MPa xuống 205,0 MPa (giảm 3,4%) và hạ thấp ứng suất kéo pháp phương X từ 231,27 MPa xuống 203,78 MPa, qua đó hạn chế tối đa nguy cơ nứt chân mối hàn.

Kết luận

  • Luận văn đã xây dựng thành công mô hình phần tử hữu hạn 3D hoàn chỉnh trên phần mềm SYSWELD, mô phỏng chính xác trường nhiệt độ và ứng xử cơ nhiệt của liên kết hàn góc chữ T hợp kim nhôm A5083 dày 8 mm.
  • Xác định trường nhiệt độ đỉnh đạt 1174°C tại tâm vũng hàn và làm rõ chu trình nhiệt tại các vị trí vùng ảnh hưởng nhiệt cách chân mối hàn từ 2 mm đến 6 mm với dải nhiệt từ 475°C đến 585°C.
  • Định lượng chính xác độ võng biến dạng cực đại đạt 3,12 mm tại vị trí giữa mép tự do của tấm hàn, tương thích chặt chẽ với kết quả đo thực nghiệm (3,30 mm) với sai số chỉ 5,45%.
  • Đánh giá toàn diện bức tranh phân bố ứng suất dư, xác định mức ứng suất Von Mises đỉnh đạt 205,0 MPa (vát mép hai phía) tập trung chủ yếu tại vị trí ngàm kẹp và nằm trong ngưỡng an toàn chịu lực của vật liệu.
  • Đề xuất bộ thông số công nghệ hàn MIG chuẩn mực (năng lượng đường 750 đến 800 J/mm, vận tốc hàn 5 mm/s) có giá trị ứng dụng thực tiễn cao cho ngành cơ khí đóng tàu vỏ nhôm tại Việt Nam.

Trong lộ trình từ 6 đến 12 tháng tới, các doanh nghiệp cơ khí đóng tàu và viện nghiên cứu nên tiếp tục mở rộng mô hình cho các kết cấu tấm cong và liên kết nhiều đường hàn phức tạp. Hãy khai thác ngay giải pháp công nghệ và mô hình mô phỏng số này để tối ưu hóa chi phí sản xuất và nâng tầm chất lượng sản phẩm công nghiệp!