Tổng quan nghiên cứu

Trong ngành cơ khí chế tạo và xây dựng kết cấu hiện đại, hơn 70% các kết cấu kim loại tấm mỏng được ghép nối bằng các phương pháp hàn nóng chảy. Tuy nhiên, hiện tượng tập trung nhiệt cục bộ trong quá trình hàn thường gây ra các vấn đề nghiêm trọng như biến dạng cong vênh, ứng suất dư và nứt mối hàn, chiếm tới khoảng 35% tổng số nguyên nhân gây phế phẩm liên kết hàn trong sản xuất thực tế. Xuất phát từ thực tiễn đó, luận văn thạc sĩ khoa học chuyên ngành Công nghệ hàn tại Trường Đại học Bách Khoa Hà Nội do tác giả Ngô Thị Thảo thực hiện dưới sự hướng dẫn của TS. Nguyễn Tiến Dương vào năm 2011 đã tập trung giải quyết bài toán: "Nghiên cứu, mô phỏng và tính toán trường nhiệt độ khi hàn tấm mỏng".

Mục tiêu cốt lõi của nghiên cứu là xây dựng mô hình toán học và ứng dụng công nghệ mô phỏng số để xác định chính xác phân bố trường nhiệt độ, kích thước vũng hàn nóng chảy, chiều rộng vùng ảnh hưởng nhiệt cũng như chu trình nung nóng và làm nguội của liên kết hàn giáp mối và hàn đắp. Phạm vi nghiên cứu tập trung vào vật liệu thép cacbon thấp CT38 với phương pháp hàn hồ quang kim loại trong môi trường khí bảo vệ ở vị trí hàn 1G, được tiến hành thực nghiệm tại xưởng cơ khí Trường Đại học Sư phạm Kỹ thuật Hưng Yên. Kết quả nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc, giúp các kỹ sư cơ khí kiểm soát chính xác tốc độ nguội trong dải nhiệt độ tới hạn từ 800°C xuống 500°C, từ đó tối ưu hóa chế độ hàn công nghệ, giảm thiểu hơn 40% biến dạng dư và nâng cao độ bền mỏi cho các kết cấu tấm mỏng trong công nghiệp.

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 được xây dựng dựa trên nền tảng lý thuyết truyền nhiệt của Fourier kết hợp với mô hình nguồn nhiệt đường di chuyển trên tấm phẳng mỏng theo lời giải giải tích kinh điển của Rosenthal. Để xác định phân bố nhiệt độ ở trạng thái giả ổn định xung quanh nguồn nhiệt hàn chuyển động với vận tốc không đổi, lý thuyết ứng dụng hệ thống hàm Bessel loại một và loại hai biến đổi bậc 0. Bên cạnh đó, mô hình nhiệt - luyện kim được áp dụng để giải thích sự biến đổi tổ chức tế vi pha rắn khi kim loại bị nung nóng vượt quá các điểm tới hạn và làm nguội nhanh.

Khung lý thuyết của luận văn xoay quanh 5 khái niệm chuyên ngành then chốt:

  • Trường nhiệt độ: Tập hợp các giá trị nhiệt độ biến thiên theo không gian ba chiều và thời gian trong suốt quá trình nung nóng và làm nguội của vật thể hàn.
  • Mặt đẳng nhiệt và đường đẳng nhiệt: Tập hợp các điểm trong không gian hoặc trên bề mặt vật hàn có cùng một giá trị nhiệt độ tại một thời điểm xác định.
  • Gradien nhiệt độ: Vectơ biểu thị tốc độ và phương hướng biến thiên nhiệt độ lớn nhất, luôn vuông góc với tiếp tuyến của mặt đẳng nhiệt.
  • Vùng ảnh hưởng nhiệt: Vùng kim loại cơ bản tiếp giáp với vũng hàn nóng chảy bị biến đổi tổ chức tế vi và cơ tính dưới tác động của chu trình nhiệt hàn mà không bị nóng chảy hoàn toàn.
  • Chu trình nhiệt và thời gian nguội: Quy luật biến thiên nhiệt độ theo thời gian tại một điểm xác định, trong đó thời gian nguội từ 800°C xuống 500°C là thông số quyết định trực tiếp đến tổ chức pha và độ dẻo dai của mối hàn.

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

Nguồn dữ liệu thực nghiệm được thu thập từ bộ mẫu gồm 15 mẫu thử thép cacbon thấp CT38 dạng tấm mỏng có kích thước hình học đồng nhất. Phương pháp chọn mẫu có chủ đích được áp dụng nhằm đại diện điển hình cho nhóm kết cấu tấm mỏng chịu tải trọng phổ biến trong công nghiệp chế tạo cơ khí, đảm bảo khả năng kiểm soát chặt chẽ các thông số đầu vào.

Quy trình nghiên cứu kết hợp giữa tính toán giải tích, mô phỏng số và thực nghiệm kiểm chứng:

  • Thực nghiệm hàn: Sử dụng nguồn hàn GMAW với khí bảo vệ chuyên dụng, thực hiện các đường hàn đắp và hàn giáp mối ở vị trí bằng 1G trên các tấm thép CT38. Hệ thống đầu đo nhiệt độ được bố trí tại các vị trí cách trục mối hàn các khoảng cách xác định để ghi nhận chu trình nhiệt thực tế.
  • Phân tích tế vi và đo kích thước hình học: Mẫu hàn sau khi hoàn thiện được cắt ngang, mài bóng cơ học và tẩm thực bằng dung dịch hóa chất Nital nồng độ 4% để làm hiện rõ ranh giới vũng nóng chảy và vùng ảnh hưởng nhiệt. Kích thước chiều rộng mối hàn và chiều rộng vùng ảnh hưởng nhiệt được đo đạc chính xác qua thiết bị quang học chuyên dụng.
  • Mô phỏng số và tính toán tự động: Thiết lập thuật toán giải tích trên bảng tính Microsoft Excel để tính toán nhanh các phương trình Bessel và nghiệm Rosenthal. Đồng thời, xây dựng mô hình phần tử hữu hạn 2D và 3D trên phần mềm chuyên dụng Ansys để mô phỏng sự phân bố đường đẳng nhiệt và động thái truyền nhiệt theo thời gian thực. Việc lựa chọn kết hợp Ansys và Excel giúp rút ngắn thời gian tính toán công nghệ từ nhiều ngày xuống còn vài phút, đồng thời cung cấp hình ảnh trực quan về trường nhiệt động học. Toàn bộ timeline nghiên cứu được tiến hành tuần tự trong vòng 12 tháng từ khâu mô hình hóa toán học đến kiểm nghiệm thực tế.

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

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

Quá trình tính toán lý thuyết và đối chứng thực nghiệm đã đưa ra 4 phát hiện quan trọng có giá trị thực tiễn cao:

  • Thứ nhất, trường nhiệt độ trên tấm mỏng đạt trạng thái giả ổn định rất nhanh chỉ sau khoảng 3 đến 5 giây kể từ thời điểm mồi hồ quang hàn. Khi đó, hình dạng các đường đẳng nhiệt trở nên cố định tương đối so với nguồn nhiệt di động, tạo thành các đường elip kéo dài về phía sau vũng hàn.
  • Thứ hai, kích thước chiều rộng vũng hàn nóng chảy đo được trên mẫu thực nghiệm đạt giá trị trung bình 6,8 mm, trong khi kết quả mô phỏng bằng phần mềm Ansys cho giá trị 7,1 mm và tính toán giải tích trên Excel cho giá trị 7,2 mm. Mức sai lệch giữa mô phỏng số và thực nghiệm chỉ là 4,41%, khẳng định độ chính xác cao của mô hình phần tử hữu hạn.
  • Thứ ba, chiều rộng vùng ảnh hưởng nhiệt đo đạc từ mẫu tẩm thực hóa chất đạt 2,3 mm, tương thích chặt chẽ với kết quả tính toán lý thuyết là 2,45 mm, sai số tương đối chỉ ở mức 6,52%.
  • Thứ tư, thời gian nguội từ 800°C xuống 500°C ghi nhận được đạt khoảng 4,8 giây. Tốc độ nguội này hoàn toàn nằm trong dải an toàn cho phép của thép CT38, ngăn ngừa triệt để sự xuất hiện của pha mactenxit tôi cứng giòn. Đồng thời, gradien nhiệt độ theo phương ngang vuông góc với trục mối hàn giảm cực nhanh; ở vị trí cách tâm mối hàn 12 mm, nhiệt độ đỉnh cực đại giảm từ trên 1500°C xuống còn 320°C, tương đương mức giảm nhiệt lượng tác động lên tới hơn 78%.

Thảo luận kết quả

Sự tương quan chặt chẽ giữa số liệu thực nghiệm và kết quả mô phỏng số (sai số tổng thể luôn dưới ngưỡng 7%) khẳng định tính đúng đắn của việc áp dụng mô hình nguồn nhiệt đường Rosenthal cho liên kết hàn tấm mỏng. Nguyên nhân của mức chênh lệch nhỏ khoảng 4% đến 6% chủ yếu xuất phát từ việc mô hình toán học giải tích đã đơn giản hóa một số điều kiện biên, xem các thông số nhiệt vật lý của vật liệu là hằng số, trong khi trên thực tế nhiệt dung riêng thể tích và độ dẫn nhiệt của thép CT38 biến thiên phi tuyến theo nhiệt độ (nhiệt dung riêng tăng từ 0,11 đến 0,22 cal/g.°C khi nhiệt độ tăng từ nhiệt độ phòng lên điểm nóng chảy).

Các dữ liệu nghiên cứu được biểu diễn trực quan qua biểu đồ đường cong phân bố nhiệt độ theo không gian và bảng đối chiếu thông số hình học. Biểu đồ trường nhiệt độ 2D xuất từ Ansys thể hiện rõ rệt dạng elip bất đối xứng của các đường đẳng nhiệt nóng chảy 1500°C, đường đẳng nhiệt giới hạn vùng ảnh hưởng nhiệt 727°C và đường đẳng nhiệt 500°C. Bảng so sánh đa thông số giữa tính toán giải tích, mô phỏng số và đo đạc thực tế phản ánh toàn diện tác động của vận tốc hàn và công suất hồ quang tới sự thu hẹp hay mở rộng của vùng ảnh hưởng nhiệt. Kết quả này bổ sung cơ sở dữ liệu quan trọng mà các nghiên cứu trước đây về hàn tấm dày chưa giải quyết trọn vẹn.

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

Dựa trên các kết luận rút ra từ mô hình trường nhiệt độ, 4 khuyến nghị cụ thể được đề xuất nhằm ứng dụng trực tiếp vào sản xuất công nghiệp:

  • Ứng dụng bảng tính tự động Excel và phần mềm Ansys vào quy trình thiết kế công nghệ hàn tại các nhà máy cơ khí nhằm kiểm soát tốc độ hàn trong khoảng tối ưu từ 25 đến 35 cm/phút, giúp giảm 40% biến dạng uốn cong vênh ở các kết cấu tấm mỏng trong vòng 6 tháng tới. Chủ thể thực hiện: Kỹ sư công nghệ và chuyên viên thiết kế đồ gá hàn.
  • Chuẩn hóa thông số chế độ hàn GMAW đối với thép cacbon thấp CT38 với dòng điện hàn duy trì từ 120A đến 160A và điện áp hồ quang từ 18V đến 22V, đảm bảo thời gian nguội từ 800°C xuống 500°C luôn đạt trên 4,5 giây để triệt tiêu hoàn toàn nguy cơ nứt nguội mối hàn. Chủ thể thực hiện: Bộ phận quản lý chất lượng và kỹ thuật viên vận hành thiết bị hàn.
  • Nâng cấp mô hình mô phỏng trường nhiệt độ 2D hiện tại lên mô hình không gian 3D kết hợp phân tích ứng suất và biến dạng nhiệt trong lộ trình 12 tháng, phục vụ việc dự báo trước độ co ngót và ứng suất dư cục bộ cho các kết cấu hàn phức tạp. Chủ thể thực hiện: Nhóm nghiên cứu mô phỏng tại các viện cơ khí và phòng nghiên cứu phát triển của doanh nghiệp.
  • Xây dựng hệ thống tài liệu hướng dẫn và cơ sở dữ liệu tra cứu nhanh trường nhiệt độ cho các nhóm vật liệu tấm mỏng khác như thép không gỉ SUS304 và hợp kim nhôm trong thời gian 18 tháng tới. Chủ thể thực hiện: Hiệp hội Công nghệ hàn phối hợp cùng các trường đại học kỹ thuật chuyên ngành.

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

  • Kỹ sư công nghệ hàn và chuyên viên R&D cơ khí: Khai thác phương pháp tính toán giải tích và mô phỏng Ansys để thiết lập chế độ hàn tối ưu, rút ngắn 50% thời gian thử nghiệm thực tế và giảm thiểu chi phí phôi mẫu khi thiết kế quy trình hàn tấm mỏng.
  • Giảng viên và học viên cao học chuyên ngành Cơ khí, Luyện kim và Công nghệ Vật liệu: Sử dụng toàn bộ khung lý thuyết truyền nhiệt hàn, mô hình giải tích nguồn nhiệt Rosenthal và thuật toán lập trình bảng tính làm giáo trình tham khảo chuyên sâu trong giảng dạy và nghiên cứu học thuật.
  • Doanh nghiệp sản xuất kết cấu kim loại tấm, đóng tàu và chế tạo ô tô: Áp dụng quy trình kiểm soát trường nhiệt và thời gian nguội để nâng cao năng suất dây chuyền hàn tự động, giảm tỷ lệ sản phẩm lỗi hỏng do biến dạng nhiệt xuống dưới mức 3%.
  • Chuyên viên kiểm định chất lượng và giám định mối hàn: Sử dụng dữ liệu dự báo kích thước vùng ảnh hưởng nhiệt để định vị chính xác vị trí kiểm tra khuyết tật vi mô, độ cứng bề mặt và phân bố tổ chức kim loại trên các liên kết hàn quan trọng.

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

Truyền nhiệt khi hàn tấm mỏng có điểm gì khác biệt cốt lõi so với hàn tấm dày? Khi hàn tấm mỏng, sự truyền nhiệt diễn ra chủ yếu theo mô hình hai chiều trên mặt phẳng tấm, dòng nhiệt truyền theo phương chiều dày gần như đồng nhất và tổn thất nhiệt qua bề mặt ra môi trường chiếm tỷ trọng lớn hơn. Do đó, các đường đẳng nhiệt bị kéo dài hơn về phía sau nguồn nhiệt và tốc độ nguội thường chậm hơn khoảng 20% đến 30% so với hàn trên vật thể khối hoặc tấm dày.

Tại sao khoảng nhiệt độ nguội từ 800°C xuống 500°C lại quyết định độ bền của liên kết hàn? Khoảng nhiệt độ từ 800°C xuống 500°C là dải nhiệt diễn ra quá trình chuyển pha từ austenit thành các tổ chức pha ferit, peclit, bainit hoặc mactenxit trong thép cacbon. Nếu thời gian nguội quá ngắn (dưới 2 giây), mối hàn dễ hình thành cấu trúc mactenxit giòn gây nứt nóng hoặc nứt nguội; nếu duy trì được khoảng 4 đến 6 giây như trong nghiên cứu, mối hàn đạt được tổ chức dẻo dai tối ưu.

Mức độ tin cậy của phần mềm Ansys khi mô phỏng trường nhiệt độ hàn tấm mỏng là bao nhiêu? Nghiên cứu chỉ ra sai số giữa mô phỏng Ansys và đo đạc thực nghiệm chỉ đạt 4,41% đối với chiều rộng vũng hàn (7,1 mm so với 6,8 mm) và 6,52% đối với vùng ảnh hưởng nhiệt (2,45 mm so with 2,3 mm). Mức sai số dưới 7% này chứng minh mô phỏng số trên Ansys hoàn toàn đủ độ tin cậy để ứng dụng vào thiết kế công nghệ hàn công nghiệp.

Vật liệu thép CT38 có những ưu thế và thách thức gì khi thực hiện hàn hồ quang tấm mỏng? Thép CT38 là loại thép cacbon thấp có hàm lượng cacbon dưới 0,22% nên tính hàn rất tốt và ít nhạy cảm với hiện tượng nứt nguội tự nhiên. Tuy nhiên, khi hàn ở dạng tấm mỏng, độ cứng vững của kết cấu thấp khiến thép CT38 rất dễ bị biến dạng uốn và cong vênh nếu gradien nhiệt độ phân bố không đều.

Làm thế nào để ứng dụng công cụ Excel trong luận văn vào thực tế phân xưởng mà không cần cài đặt phần mềm chuyên dụng? Bảng tính Excel được lập trình tự động hóa các công thức tích phân Rosenthal và hàm Bessel, cho phép người dùng chỉ cần nhập cường độ dòng điện (120A đến 180A), điện áp (18V đến 24V) và tốc độ hàn để nhận ngay kết quả chiều rộng vũng hàn và kích thước vùng ảnh hưởng nhiệt trong vòng 30 giây mà không cần đòi hỏi bản quyền phần mềm mô phỏng phức tạp.

Kết luận

  • Xây dựng thành công mô hình giải tích nguồn nhiệt di động kết hợp hàm Bessel trên bảng tính Excel, cho phép tính toán nhanh và chính xác trường nhiệt độ hàn tấm mỏng.
  • Ứng dụng hiệu quả phương pháp phần tử hữu hạn trên phần mềm Ansys để trực quan hóa động thái truyền nhiệt, hình dạng vũng hàn và các đường đẳng nhiệt 2D/3D.
  • Thực nghiệm kiểm chứng trên thép cacbon thấp CT38 với phương pháp hàn GMAW cho độ sai lệch kích thước hình học so với mô phỏng số ở mức rất thấp, chỉ từ 4,41% đến 6,52%.
  • Xác định chính xác quy luật phân bố vùng ảnh hưởng nhiệt có chiều rộng 2,3 mm và kiểm soát thời gian nguội từ 800°C xuống 500°C đạt mức tối ưu 4,8 giây, đảm bảo cơ tính mối hàn.
  • Đóng góp giải pháp công nghệ toàn diện giúp các doanh nghiệp cơ khí làm chủ bài toán kiểm soát nhiệt độ, hạn chế tối đa khuyết tật và ứng suất dư khi hàn kết cấu mỏng.

Lộ trình nghiên cứu mở rộng trong 6 đến 12 tháng tới sẽ tập trung hoàn thiện mô hình tương tác nhiệt - cơ đa trường trên các loại vật liệu hợp kim đặc biệt. Hãy chủ động ứng dụng ngay các giải pháp mô phỏng số và kiểm soát trường nhiệt độ từ công trình nghiên cứu này để nâng cao năng suất, tối ưu hóa chi phí sản xuất và nâng tầm chất lượng cho các sản phẩm cơ khí tấm mỏng trong doanh nghiệp.