Tổng quan nghiên cứu

Trong ngành cơ khí chế tạo và kết cấu thép hiện đại, công nghệ hàn nóng chảy chiếm hơn 65% khối lượng công việc liên kết kim loại. Tuy nhiên, sự tập trung nguồn nhiệt cục bộ với nhiệt độ hồ quang vượt quá 6000°C tạo ra chu trình nhiệt khắc nghiệt, gây biến thiên nhiệt độ phức tạp và dẫn đến hiện tượng ứng suất dư cùng biến dạng nhiệt nghiêm trọng. Luận văn thạc sĩ khoa học chuyên ngành Hàn và Công nghệ Kim loại thực hiện tại Trường Đại học Bách khoa Hà Nội năm 2009 tập trung giải quyết bài toán mô phỏng hình học vũng hàn và sự phân bố trường nhiệt độ trong quá trình hàn hồ quang kim loại trong môi trường khí bảo vệ (GMAW/MAG).

Mục tiêu cụ thể của nghiên cứu là xây dựng công cụ tính toán số học tự động nhằm xác định chính xác kích thước hình học vũng hàn (chiều sâu ngấu, chiều rộng bề mặt), thể tích vũng hàn nóng chảy, quỹ đạo các đường bao đẳng nhiệt, chiều rộng vùng ảnh hưởng nhiệt (HAZ) cũng như thông số tốc độ nguội của liên kết hàn. Đề tài khảo sát trên đối tượng vật liệu thép carbon thấp CT38 với hai trường hợp hình học điển hình gồm tấm mỏng chiều dày 5 mm hàn 1 lớp và tấm dày 10 mm vát mép chữ V hàn 3 lớp tại vị trí giáp mối 1G. Nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm thiểu trên 40% chi phí thử nghiệm hàn mẫu thực tế, rút ngắn 60% thời gian tối ưu hóa chế độ công nghệ và khống chế sai lệch kích thước vũng hàn thực nghiệm xuống dưới 6%.

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 ứng dụng nền tảng lý thuyết truyền nhiệt trong hàn nóng chảy, trọng tâm là mô hình nguồn nhiệt di động của Rosenthal và lý thuyết phân bố nhiệt của Grong trong môi trường vật thể bán vô hạn và tấm phẳng hữu hạn. Quá trình truyền nhiệt dẫn xuất từ định luật Fourier được giải ở trạng thái cân bằng nhiệt giả ổn định (quasi-steady state), mô tả mối quan hệ vi phân giữa gradient nhiệt độ, vận tốc dịch chuyển nguồn nhiệt và các thông số nhiệt vật lý của vật liệu.

Bốn khái niệm cốt lõi được định hình và phân tích xuyên suốt đề tài bao gồm:

  1. Trường nhiệt độ hàn: Sự phân bố giá trị nhiệt độ tức thời trong không gian 2D và 3D của kết cấu hàn theo hàm vị trí và thời gian.
  2. Đường bao đẳng nhiệt: Quỹ đạo tập hợp tất cả các điểm có cùng nhiệt độ tại một thời điểm khảo sát, đặc biệt là đường bao nhiệt độ nóng chảy ở ngưỡng khoảng 1500°C của thép CT38.
  3. Vùng ảnh hưởng nhiệt (HAZ): Khu vực kim loại cơ bản tiếp giáp vũng hàn bị biến đổi cấu trúc tế vi và cơ tính do tác động của nhiệt độ vượt quá 727°C mà không bị nóng chảy.
  4. Thời gian và tốc độ nguội t8/5: Khoảng thời gian làm nguội của mối hàn trong dải nhiệt độ nguy hiểm từ 800°C xuống 500°C, đóng vai trò quyết định đến tổ chức pha kim loại và nguy cơ nứt nguội. Các thông số nhiệt lý như hệ số dẫn nhiệt khoảng 45 W/m.K và nhiệt dung riêng 480 J/kg.K được thiết lập chuẩn xác trong thuật toán.

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

Nghiên cứu áp dụng phương pháp kết hợp mô phỏng số học giải tích trên máy tính và kiểm chứng thực nghiệm đối chứng đa thông số. Nguồn dữ liệu đầu vào bao gồm các thông số chế độ hàn tự động GMAW: cường độ dòng điện hàn trong khoảng 140A đến 260A, điện áp hồ quang từ 18V đến 28V, vận tốc hàn từ 20 cm/phút đến 45 cm/phút, kết hợp hệ số hiệu dụng nhiệt hồ quang đạt 75% đến 85%.

Phương pháp chọn mẫu thực nghiệm được thực hiện trên 10 mẫu phôi thép tấm CT38 chuẩn quy cách công nghiệp. Mẫu tấm mỏng có kích thước chiều dày 5 mm hàn giáp mối không vát mép; mẫu tấm dày 10 mm được gia công vát mép góc 60 độ, khe hở đáy 1.5 mm và thực hiện hàn nhiều lớp. Thiết bị thực nghiệm sử dụng máy hàn Inverter 19M 350 tích hợp trên Robot hàn tự động 6 bậc tự do AX-MV6, sử dụng dây hàn ER70S-6 đường kính 1.2 mm và khí bảo vệ CO2 tinh khiết với lưu lượng định mức 14 lít/phút.

Thuật toán mô phỏng được lập trình toàn diện trên môi trường MATLAB, xây dựng giao diện tương tác tính toán trực quan, cho phép giải nhanh phương trình truyền nhiệt vi phân và xuất biểu đồ đẳng nhiệt tự động. Nghiên cứu lựa chọn giải pháp giải tích mô phỏng trên MATLAB thay vì phương pháp phần tử hữu hạn thương mại nặng nề nhằm rút ngắn 70% thời gian xử lý dữ liệu mà vẫn đảm bảo độ tin cậy kết quả đạt trên 94%.

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

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

Quá trình mô phỏng và đối soát thực nghiệm đã mang lại những kết quả định lượng mang tính chuẩn hóa cao cho công nghệ hàn GMAW thép carbon:

Thứ nhất, đối với liên kết hàn tấm mỏng chiều dày 5 mm (hàn 1 lớp không vát mép), mô hình số học xác định bề rộng vũng hàn đạt 7.2 mm, chiều sâu ngấu 5.1 mm và bề rộng vùng ảnh hưởng nhiệt HAZ là 2.4 mm. Kết quả đo thực nghiệm tế vi cho thấy kích thước bề rộng thực tế đạt 7.5 mm và chiều sâu ngấu đạt 5.3 mm, độ lệch giữa mô phỏng và thực nghiệm chỉ ở mức 4.0% đến 4.3%.

Thứ hai, đối với liên kết hàn tấm dày 10 mm (vát mép chữ V, hàn 3 lớp), mô phỏng trường nhiệt 3D xác định biên dạng vũng hàn có bề rộng phủ ngoài đạt 12.8 mm và chiều sâu nóng chảy tích lũy đạt 10.2 mm. Gradient nhiệt độ tại tâm vũng hàn vượt mức 2200°C và giảm dốc xuống ngưỡng chuyển pha 727°C tại biên HAZ, với sai số kích thước so với mẫu cắt thực tế chỉ 5.6%.

Thứ ba, sự phân kỳ mạnh mẽ về tốc độ nguội giữa hai chế độ truyền nhiệt được xác lập rõ ràng. Thời gian nguội t8/5 ở tấm mỏng 5 mm đạt 8.5 giây (do cơ chế truyền nhiệt phẳng 2D thoát nhiệt chậm hơn), trong khi ở tấm dày 10 mm thời gian t8/5 chỉ còn 4.2 giây (do khối kim loại lớn gây tản nhiệt nhanh 3D), làm tốc độ nguội tăng lên hơn 102%.

Thứ tư, nghiên cứu xác định quy luật tương quan phi tuyến giữa năng lượng đường và thể tích kim loại nóng chảy. Khi tăng năng lượng đường từ 8.5 kJ/cm lên 14.2 kJ/cm, thể tích vũng hàn tăng thêm 42.5%, đồng thời độ rộng dải biến đổi nhiệt HAZ mở rộng thêm 31.8%, làm tăng nguy cơ thô đại hạt tổ chức kim loại.

Thảo luận kết quả

Nguyên nhân của sự sai khác nhỏ giữa mô phỏng và thực tế (dưới ngưỡng sai số 6%) xuất phát từ việc mô hình giả định các hệ số truyền nhiệt và nhiệt dung riêng là hằng số đẳng hướng, đồng thời chưa xét đến hiện tượng đối lưu chất lỏng mạnh mẽ sinh ra bởi lực điện từ (lực Lorentz) trong lòng bể chảy kim loại lỏng.

Khi so sánh với các công bố khoa học sử dụng nguồn nhiệt điểm đơn thuần, giải thuật tích hợp nguồn nhiệt diện và nguồn nhiệt khối trong nghiên cứu này cho độ chính xác cao hơn 15% tại khu vực quá độ nhiệt. Dữ liệu trường nhiệt được hiển thị trực quan thông qua các biểu đồ 2D dạng elip đẳng nhiệt kéo dài về phía sau hướng dịch chuyển của mỏ hàn và biểu đồ 3D dạng mặt cong đẳng nhiệt dạng quả lê. Bảng tổng hợp so sánh các thông số hình học (chiều rộng B, chiều sâu H, diện tích tiết diện F) chứng minh tính tương thích cao giữa phần mềm tính toán trên MATLAB và thực tế xưởng sản xuất cơ khí.

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

Dựa trên kết quả mô phỏng và thực nghiệm, luận văn đề xuất bốn nhóm giải pháp kỹ thuật nhằm nâng cao chất lượng liên kết hàn và hiện đại hóa quy trình sản xuất:

  1. Thiết lập quy chuẩn bộ thông số hàn tự động trên Robot hàn AX-MV6 với dòng điện 180A đến 210A, điện áp 22V đến 24V, tốc độ dịch chuyển 30 cm/phút đến 35 cm/phút cho thép CT38 chiều dày 5 mm. Mục tiêu giữ chiều rộng vùng HAZ dưới ngưỡng 2.5 mm, giảm 25% độ biến dạng co ngót ngang của kết cấu. Kỹ sư công nghệ hàn tại các doanh nghiệp kết cấu cơ khí cần triển khai áp dụng trong khung thời gian 3 đến 6 tháng.
  2. Tích hợp trực tiếp thuật toán mô phỏng trường nhiệt MATLAB vào hệ thống điều khiển PLC của dây chuyền hàn tự động nhằm tự động bù trừ nhiệt độ khi hàn nhiều lớp. Đích nhắm giảm 80% tỷ lệ phế phẩm do quá nhiệt hoặc thiếu ngấu, thực hiện bởi các phòng R&D và viện nghiên cứu tự động hóa trong lộ trình 12 tháng.
  3. Áp dụng chế độ gia nhiệt sơ bộ ở nhiệt độ 100°C đến 150°C đối với các liên kết thép dày trên 10 mm hàn nhiều lớp. Giải pháp này giúp kéo dài thời gian nguội t8/5 từ 4.2 giây lên dải tối ưu 10 đến 14 giây, loại bỏ 95% nguy cơ hình thành pha Martensite giòn gây nứt nguội. Bộ phận kiểm soát chất lượng (QC) tại các nhà máy đóng tàu và bồn bể áp lực cần áp dụng ngay trong quý 1 năm tiếp theo.
  4. Chuẩn hóa hình học vát mép góc 60 độ với khe hở đáy từ 1.5 mm đến 2.0 mm cho mối hàn giáp mối 3 lớp, giúp tăng 15% hệ số điền đầy kim loại mối hàn và tiết kiệm 12% khối lượng dây hàn ER70S-6 tiêu hao. Tổ trưởng sản xuất và kỹ thuật viên chế tạo máy chịu trách nhiệm giám sát thi công thường nhật.

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

Nội dung toàn diện của luận văn mang lại giá trị học thuật và ứng dụng sâu rộng cho nhiều nhóm đối tượng:

  1. Sinh viên và học viên cao học ngành Kỹ thuật Cơ khí, Hàn và Công nghệ Kim loại: Tài liệu cung cấp cơ sở lý thuyết chuẩn xác về truyền nhiệt trong hàn, cùng hướng dẫn chi tiết về phương pháp xây dựng thuật toán mô phỏng trường nhiệt bằng phần mềm tính toán kỹ thuật.
  2. Kỹ sư công nghệ hàn và chuyên viên thiết kế kết cấu thép: Tài liệu hỗ trợ trực tiếp việc tính toán dự báo biến dạng, ứng suất dư và lựa chọn chế độ hàn GMAW/MAG tối ưu cho các kết cấu tấm dày từ 5 mm đến 10 mm mà không cần thử nghiệm phá hủy tốn kém.
  3. Giảng viên và cán bộ nghiên cứu tại các trường đại học kỹ thuật: Luận văn là tài liệu tham khảo giá trị phục vụ công tác giảng dạy các học phần Lý thuyết quá trình hàn, Tự động hóa công nghệ hàn và Mô phỏng số trong cơ khí.
  4. Quản lý sản xuất tại các nhà máy gia công cơ khí, đóng tàu và kết cấu công trình: Hỗ trợ chuyển đổi số quy trình lập kế hoạch hàn, tối ưu hóa mức tiêu hao vật liệu và ứng dụng hiệu quả cánh tay Robot hàn công nghiệp trong sản xuất hàng loạt.

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

Mô hình nhiệt của Rosenthal và Grong trong luận văn có ưu điểm gì nổi bật? Mô hình cho phép tính toán giải tích trực tiếp trường nhiệt độ và kích thước vũng hàn ở trạng thái giả ổn định chỉ trong 2 đến 3 giây tính toán. Sai số dự báo hình học so với thực tế duy trì ở mức dưới 5.6%, giúp tiết kiệm hơn 70% thời gian phân tích so với các mô hình phần tử hữu hạn thông thường.

Tại sao tốc độ làm nguội t8/5 ở tấm dày 10 mm lại nhanh hơn đáng kể so với tấm mỏng 5 mm? Tấm dày 10 mm kích hoạt cơ chế truyền nhiệt không gian 3D, trong đó khối kim loại xung quanh đóng vai trò như một bộ tản nhiệt vô hạn. Điều này làm tốc độ tản nhiệt tăng hơn 102%, rút ngắn thời gian t8/5 xuống còn 4.2 giây so với 8.5 giây ở tấm mỏng 5 mm truyền nhiệt phẳng 2D.

Phần mềm mô phỏng viết trên MATLAB có khả năng ứng dụng thực tế ra sao? Chương trình cung cấp giao diện nhập trực tiếp các thông số dòng điện từ 140A đến 260A, điện áp, tốc độ hàn và kích thước vát mép. Hệ thống tự động xuất đồ thị đường bao đẳng nhiệt và kích thước vũng hàn, hỗ trợ kỹ sư lựa chọn chế độ hàn chính xác ngay tại phân xưởng mà không cần cài đặt phần mềm FEM phức tạp.

Làm thế nào để kiểm soát kích thước vùng ảnh hưởng nhiệt HAZ khi hàn thép CT38? Cần kiểm soát năng lượng đường trong khoảng tối ưu từ 8.5 kJ/cm đến 11.0 kJ/cm bằng cách tăng tốc độ hàn lên mức 30 cm/phút đến 35 cm/phút và duy trì điện áp hồ quang 22V đến 24V. Giải pháp này giúp khống chế bề rộng vùng HAZ dưới mức 2.5 mm, bảo toàn cơ tính liên kết.

Việc hàn 3 lớp trên tấm dày 10 mm vát mép chữ V cần lưu ý những thông số nhiệt nào? Cần chú ý kiểm soát nhiệt độ giữa các lớp hàn ở mức 120°C đến 150°C và duy trì năng lượng đường ổn định ở mỗi lượt hàn. Điều này giúp ngăn ngừa hiện tượng tích tụ nhiệt quá mức làm hạt kim loại thô đại hoặc tản nhiệt quá nhanh gây nứt tế vi chân mối hàn.

Kết luận

  • Luận văn đã xây dựng thành công chương trình tính toán số học trên MATLAB mô phỏng chính xác hình học vũng hàn và phân bố trường nhiệt độ khi hàn GMAW.
  • Kiểm chứng thực nghiệm trên thép CT38 khẳng định độ chính xác cao của mô hình với sai số kích thước vũng hàn và vùng HAZ chỉ từ 4.0% đến 5.6%.
  • Nghiên cứu làm sáng tỏ sự khác biệt bản chất về tốc độ nguội t8/5 giữa cơ chế truyền nhiệt 2D trên tấm 5 mm và truyền nhiệt 3D trên tấm 10 mm.
  • Đóng góp giải pháp công nghệ định lượng giúp tối ưu hóa chế độ hàn tự động bằng Robot AX-MV6, nâng cao chất lượng và độ bền mỏi cho kết cấu thép.
  • Lộ trình tiếp theo hướng tới tích hợp trí tuệ nhân tạo và học máy để hiệu chỉnh thông số nguồn nhiệt theo thời gian thực trên dây chuyền sản xuất công nghiệp 4.0.