Tổng quan nghiên cứu

Trong ngành cơ khí chế tạo và gia công áp lực, sự hình thành và phát triển của các lỗ hổng vi mô là nguyên nhân cốt lõi dẫn đến hiện tượng phá hủy dẻo kim loại. Theo ước tính từ các nghiên cứu cơ học vật liệu, khoảng 70% đến 85% các vết nứt tế vi trong các chi tiết kết cấu thép xuất phát từ sự tập trung ứng suất quanh các hạt tạp chất hoặc pha thứ hai. Quá trình biến dạng dẻo lớn làm các lỗ hổng ban đầu tăng trưởng, kết hợp lại và làm suy giảm nghiêm trọng khả năng chịu tải của vật thể. Khi các vi lỗ hổng xuất hiện, vật liệu kim loại đặc ban đầu chuyển hóa thành trạng thái vật liệu xốp dẻo, đòi hỏi các mô hình toán học và cơ học chuyên biệt để mô tả chính xác ứng xử cơ tính.

Luận văn thạc sĩ chuyên ngành Cơ học Kỹ thuật với đề tài "Ảnh hưởng độ xốp đến tính chất cơ học của vật liệu xốp dẻo" được thực hiện bởi học viên Dương Thị Hồng Hạnh dưới sự hướng dẫn khoa học của PGS. Trương Tích Thiện tại Trường Đại học Bách Khoa, Đại học Quốc gia Thành phố Hồ Chí Minh. Công trình được triển khai trong thời gian 11 tháng, từ ngày 16/07/2015 đến ngày 17/06/2016 và chính thức hoàn thành vào tháng 07/2016.

Mục tiêu trọng tâm của nghiên cứu là xây dựng hàm chảy giải tích cho vật liệu xốp dẻo thông qua việc mô hình hóa khối vật liệu von Mises dẻo lý tưởng chứa các lỗ hổng hình cầu phân bố tuần hoàn đều ba chiều. Thay vì tìm kiếm hàm chảy trực tiếp theo độ xốp và áp suất thủy tĩnh như các cách tiếp cận truyền thống, tác giả đã phân tích mặt chảy dựa trên hai trạng thái ứng suất đặc biệt: trạng thái tải thuần túy ứng suất lệch và trạng thái tải thuần túy ứng suất thủy tĩnh. Kết quả của luận văn mang ý nghĩa khoa học và thực tiễn sâu sắc, cung cấp công cụ dự báo chính xác độ mềm hóa và giới hạn phá hủy của vật liệu với ngưỡng độ xốp tới hạn đạt khoảng 80%, hỗ trợ tối ưu hóa quy trình gia công kim loại và giảm thiểu tỷ lệ nứt hỏng kết cấu trong sản xuất 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 trên nền tảng của ba lý thuyết cơ học biến dạng dẻo kinh điển: lý thuyết dẻo von Mises cho vật liệu kim loại nền đẳng hướng dẻo lý tưởng, lý thuyết dẻo cho vật liệu xốp dựa trên mô hình Gurson mở rộng và nguyên lý công ảo trong cơ học môi trường liên tục. Vật liệu nền kim loại được giả định có ứng suất chảy ban đầu là 300 MPa, tuân theo tiêu chuẩn chảy dẻo von Mises với bề mặt chảy trong không gian ứng suất chính là một hình trụ tròn có bán kính bằng căn bậc hai của hai phần ba nhân với ứng suất chảy.

Mô hình nghiên cứu vi mô đại diện được thiết lập là một khối vật liệu hình lập phương kích thước cạnh a, chứa các lỗ hổng hình cầu bán kính r phân bố đều theo ba phương không gian. Trong khung lý thuyết này, năm khái niệm và thông số cơ học then chốt được áp dụng xuyên suốt gồm:

  1. Ứng suất tương đương von Mises biểu diễn cường độ của tensor ứng suất lệch.
  2. Ứng suất thủy tĩnh đặc trưng cho trạng thái nén hoặc kéo ba trục đồng đều.
  3. Bất biến thứ hai của tensor ứng suất lệch.
  4. Độ xốp thể tích được xác định bằng tỷ số giữa thể tích lỗ hổng và thể tích toàn phần của phân tố đại diện.
  5. Bán kính lỗ hổng tới hạn xấp xỉ 0.7a, thời điểm các lỗ xốp mở rộng chạm vào nhau và cắt rời hoàn toàn liên kết của phân tố vật liệu.

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

Nghiên cứu kết hợp chặt chẽ giữa giải tích lý thuyết và mô phỏng số thông qua phương pháp phần tử hữu hạn ba chiều. Để kiểm chứng các biểu thức giải tích được thiết lập, tác giả đã xây dựng ba mô hình phần tử hữu hạn đại diện cho một phần mười sáu phân tố đối xứng của khối vật liệu, bao gồm các cấu hình mạng: 3x3x3 lỗ hổng (tương đương 27 lỗ hổng toàn khối), 4x4x4 lỗ hổng (64 lỗ hổng) và 5x5x5 lỗ hổng (125 lỗ hổng).

Phương pháp chọn mẫu dữ liệu mô phỏng được thực hiện trên 7 mức độ xốp thể tích khác nhau, trải dải từ 0.01 đến 0.50 (tương ứng từ 1% đến 50% độ xốp). Lý do lựa chọn phương pháp phân tích phần tử hữu hạn với phần tử khối biến dạng dẻo suy biến 20 nút là nhằm nắm bắt chính xác sự tập trung ứng suất và biến dạng dẻo cục bộ phi tuyến cao tại vùng biên lỗ hổng, điều mà các công thức giải tích đơn giản hóa không thể phản ánh đầy đủ. Dữ liệu tải và chuyển vị được thu thập dưới hai chế độ đặt tải độc lập: tải đối xứng thuần túy ứng suất lệch (tổng ba ứng suất chính bằng 0) và tải thuần túy thủy tĩnh (ba ứng suất chính bằng nhau). Toàn bộ quá trình tính toán số và đối chuẩn lý thuyết được thực hiện liên tục trong khung thời gian 11 tháng tại phòng thí nghiệm tính toán Cơ học Kỹ thuật.

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

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

Nghiên cứu đã đạt được bốn kết quả mang tính đột phá về mặt lý thuyết và mô phỏng số:

Thứ nhất, luận văn xác định được giá trị độ xốp tới hạn khi vật liệu mất hoàn toàn khả năng chịu tải là 0.8 (tương đương 80%). Khi bán kính lỗ hổng tăng từ 0.5a đến 0.7a, thể tích lỗ hổng tràn ra ngoài biên phân tố lập phương và tại giá trị bán kính r xấp xỉ 0.707a, cấu trúc vật liệu nền bị cắt đứt hoàn toàn. Trong tính toán giải tích đơn giản hóa, tác giả chọn giá trị bán kính tính toán là 0.6a tương ứng với độ xốp tới hạn 0.8.

Thứ hai, thiết lập thành công quy luật giảm tuyến tính của tỷ số ứng suất chảy thuần túy ứng suất lệch theo công thức: tỷ số ứng suất chảy tương đương trên ứng suất chảy nền bằng 1 trừ đi 1.25 lần độ xốp. Kết quả này chứng minh về mặt giải tích nguồn gốc của hệ số hiệu chỉnh q1 bằng 1.25, hoàn toàn trùng khớp với các kết quả mô phỏng số thực nghiệm của Koplik và Needleman (1988).

Thứ ba, xây dựng biểu thức giải tích cải tiến cho ứng suất chảy thuần túy thủy tĩnh kết hợp hàm hiệu chỉnh số mũ bậc 1: tỷ số áp suất thủy tĩnh trên ứng suất chảy nền bằng hai phần ba nhân với hiệu của logarit tự nhiên của nghịch đảo độ xốp trừ đi 1.25 lần độ xốp nhân logarit của 1.25. Biểu thức này khắc phục hoàn toàn nhược điểm của mô hình Gurson gốc, vừa đảm bảo tính chính xác tại vùng độ xốp nhỏ dưới 0.05, vừa thỏa mãn điều kiện triệt tiêu ứng suất tại ngưỡng độ xốp tới hạn 0.8.

Thứ tư, phân tích định lượng quá trình mềm hóa của vật liệu trong điều kiện kéo ba trục cho thấy khi độ xốp tăng từ 0.05 lên 0.20, ứng suất chảy tương đương của vật liệu xốp dẻo bị sụt giảm hơn 35% so với vật liệu đặc; khi độ xốp đạt mức 0.40, khả năng chịu tải tương đương giảm đến hơn 65%.

Thảo luận kết quả

Cơ chế vật lý của hiện tượng mềm hóa bắt nguồn từ sự suy giảm diện tích chịu tải hiệu dụng của mạng kim loại và sự gia tăng mức độ tập trung ứng suất lệch tại các mặt cầu lân cận. Trong không gian ứng suất chính, mặt chảy của vật liệu xốp dẻo được biểu diễn dưới dạng mặt tròn xoay ellipsoid đối xứng qua trục thủy tĩnh, với bán trục trên mặt phẳng lệch tỷ lệ với ứng suất chảy lệch và đỉnh trên trục thủy tĩnh bằng căn bậc hai của ba lần ứng suất chảy thủy tĩnh.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua đồ thị so sánh hai chiều giữa tỷ số ứng suất tương đương và tỷ số ứng suất thủy tĩnh, kết hợp bảng đối chiếu sai số tương đối giữa mô hình giải tích đề xuất với ba cấu hình phần tử hữu hạn 3x3x3, 4x4x4 và 5x5x5. So sánh với các công trình quốc tế, mô hình Gurson truyền thống (q1 = 1) chỉ phù hợp với lỗ xốp đơn độc trong phân tố cầu vô hạn và đánh giá thấp độ mềm hóa từ 15% đến 20% ở dải xốp trung bình. Ngược lại, mô hình Tvergaard với q1 = 1.5 lại có xu hướng dự báo độ suy giảm cơ tính quá mức ở vùng độ xốp nhỏ dưới 0.05. Mô hình ellipsoid đề xuất của luận văn với q1 = 1.25 thể hiện sự tương thích vượt trội với dữ liệu phần tử hữu hạn, duy trì sai số tương đối luôn dưới mức 4.5% trên toàn dải biến dạng.

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

Dựa trên các kết quả giải tích và mô phỏng số đã đạt được, luận văn đưa ra bốn nhóm giải pháp và khuyến nghị ứng dụng thực tiễn:

  1. Tích hợp mô hình hàm chảy ellipsoid cải tiến vào các phần mềm mô phỏng phần tử hữu hạn thương mại như Abaqus, ANSYS và LS-DYNA dưới dạng chương trình con người dùng UMAT/VUMAT. Mục tiêu là nâng cao độ chính xác dự báo rạn nứt dẻo trong các bài toán va chạm và gia công tạo hình lên trên 92%. Thời gian triển khai dự kiến từ 6 đến 12 tháng, do các kỹ sư tính toán mô phỏng CAE và các chuyên gia cơ học đảm nhiệm.

  2. Áp dụng quy chuẩn kiểm soát giới hạn độ xốp tới hạn dưới ngưỡng an toàn 0.12 trong các công nghệ gia công biến dạng dẻo nóng và cán kéo kim loại. Biện pháp này giúp các doanh nghiệp luyện kim và cơ khí chế tạo cắt giảm từ 25% đến 30% tỷ lệ phế phẩm do nứt dẻo tế vi. Thời gian áp dụng trong vòng 3 đến 6 tháng tại các phân xưởng dập cán tấm.

  3. Mở rộng nghiên cứu giải tích cho các cấu trúc lỗ xốp có dạng hình học bất đối xứng như ellipsoid dẹt hoặc kéo dài, đồng thời xét đến tương tác biến cứng của vật liệu nền và cơ chế sát nhập lỗ hổng. Mục tiêu là thử nghiệm thành công 3 biến thể hình học lỗ xốp phức tạp trong thời gian từ 12 đến 24 tháng, do các nhóm nghiên cứu chuyên sâu tại các trường đại học kỹ thuật thực hiện.

  4. Xây dựng cơ sở dữ liệu số thực nghiệm chuẩn hóa về thông số cơ tính của ít nhất 10 mác thép kết cấu công nghiệp phổ biến có giới hạn chảy từ 250 MPa đến 650 MPa dưới các trạng thái tải kéo ba trục khác nhau. Đề án này do các viện nghiên cứu vật liệu và trung tâm kiểm định chất lượng triển khai trong lộ trình 18 tháng.

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

Nội dung và phương pháp luận của luận văn mang lại giá trị thiết thực cho bốn nhóm đối tượng chính:

  1. Học viên cao học và nghiên cứu sinh chuyên ngành Cơ học Kỹ thuật, Cơ học Vật rắn và Kỹ thuật Vật liệu. Luận văn cung cấp phương pháp tiếp cận mẫu mực trong việc thiết lập mặt chảy dẻo từ các trạng thái tải đặc biệt và kỹ thuật mô phỏng phần tử hữu hạn 3D cho các phân tố vi mô chứa khuyết tật.

  2. Kỹ sư R&D và chuyên viên phân tích mô phỏng CAE tại các tập đoàn sản xuất ô tô, hàng không và cơ khí chế tạo máy. Luận văn cung cấp công thức hàm chảy tin cậy để thiết lập các mô hình phá hủy vật liệu trong phân tích dập sâu kim loại tấm và mô phỏng va đập kết cấu với độ tin cậy đạt trên 95%.

  3. Giảng viên và các nhà khoa học tại các trường đại học kỹ thuật và viện nghiên cứu chuyên ngành cơ học tính toán. Đây là tài liệu tham khảo học thuật giá trị phục vụ công tác giảng dạy các học phần Cơ học biến dạng dẻo, Cơ học phá hủy và Phương pháp phần tử hữu hạn nâng cao.

  4. Kỹ sư quản lý chất lượng và luyện kim tại các nhà máy sản xuất phôi thép và hợp kim. Nghiên cứu giúp kỹ sư nắm rõ cơ chế suy giảm cơ tính do mật độ lỗ xốp vi mô sinh ra từ tạp chất sunfua và cacbit, từ đó tối ưu hóa quy trình khử khí và xử lý nhiệt để duy trì độ xốp vật liệu luôn dưới ngưỡng 0.05.

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

  1. Tại sao luận văn lại chọn cách tìm hàm chảy thông qua hai ứng suất chảy đặc biệt thay vì tìm trực tiếp? Việc phân tách bài toán thành hai trạng thái tải thuần túy ứng suất lệch và thuần túy thủy tĩnh giúp xác định chính xác hai điểm mút hình học giới hạn của mặt chảy trên mặt phẳng lệch và trục thủy tĩnh. Cách tiếp cận này biến bài toán tìm hàm hai biến phức tạp thành việc xác định các tham số hình học rõ ràng của mặt ellipsoid, giúp công thức giải tích đơn giản và dễ ứng dụng hơn.

  2. Ý nghĩa vật lý thực sự của hệ số hiệu chỉnh q1 bằng 1.25 trong mô hình là gì? Hệ số q1 bằng 1.25 tương ứng trực tiếp với giá trị độ xốp tới hạn ff bằng 0.8. Đây là thời điểm mà các vi lỗ hổng hình cầu phân bố trong mạng lập phương phát triển đạt bán kính từ 0.6a đến 0.7a, khiến các phân tố vật liệu kim loại nền bị cắt rời hoàn toàn và mất 100% khả năng chịu tải, khắc phục giả định không thực tế của mô hình cũ khi cho rằng vật liệu chỉ mất tải khi độ xốp bằng 1.

  3. Mô hình đề xuất trong luận văn có ưu điểm gì vượt trội so với mô hình Gurson năm 1977? Mô hình Gurson gốc được xây dựng cho một lỗ hổng đơn độc trong phân tố hình cầu vô hạn nên đánh giá thấp mức độ tương tác giữa các lỗ xốp liền kề, dẫn đến sai số lớn khi độ xốp vượt quá 0.05. Mô hình của luận văn mô tả chính xác sự phân bố tuần hoàn 3 chiều của mạng lỗ xốp, mang lại độ chính xác cao hơn với sai số đối chuẩn phần tử hữu hạn dưới 4.5%.

  4. Kỹ thuật phần tử hữu hạn được thiết lập như thế nào để đảm bảo tính chính xác của dữ liệu đối chứng? Tác giả đã mô hình hóa 1/16 phân tố lập phương đối xứng sử dụng phần tử khối 20 nút suy biến với 3 cấu hình phân bố mạng lỗ hổng gồm 3x3x3, 4x4x4 và 5x5x5. Vật liệu nền được gán ứng suất chảy 300 MPa và khảo sát qua 7 mức độ xốp từ 0.01 đến 0.50 dưới các điều kiện biên đặt tải đối xứng nghiêm ngặt để xác định chính xác điểm chảy dẻo.

  5. Kết quả của luận văn có thể được áp dụng ngay vào quy trình dập vuốt kim loại tấm không? Có thể áp dụng trực tiếp thông qua việc cài đặt hàm chảy ellipsoid và tiêu chuẩn độ xốp tới hạn vào các thuật toán mô phỏng quá trình tạo hình kim loại. Nhờ việc dự báo chính xác vùng tập trung biến dạng dẻo và thời điểm bắt đầu nứt tế vi, quy trình công nghệ dập vuốt có thể giảm thiểu từ 20% đến 30% nguy cơ rách mép hoặc co thắt cục bộ.

Kết luận

  • Luận văn đã giải quyết thành công bài toán xác định hàm chảy giải tích cho vật liệu xốp dẻo dựa trên mô hình mạng lỗ hổng hình cầu phân bố đều 3 chiều trong không gian phân tố lập phương.
  • Chứng minh về mặt giải tích ý nghĩa vật lý của hệ số hiệu chỉnh q1 bằng 1.25 gắn liền với độ xốp tới hạn 80% (f = 0.8), thời điểm vật liệu bị cắt rời hoàn toàn bởi các vi lỗ hổng có bán kính đạt khoảng 0.7a.
  • Đề xuất hàm ứng suất chảy thủy tĩnh cải tiến với hàm mũ bậc 1, giải quyết triệt để sự mâu thuẫn giữa độ chính xác ở dải độ xốp nhỏ và điều kiện biên triệt tiêu ứng suất ở dải xốp lớn.
  • Kiểm chứng toàn diện mô hình lý thuyết với kết quả phần tử hữu hạn trên 3 cấu hình phân tố đối xứng (3x3x3, 4x4x4, 5x5x5) trên 7 dải độ xốp với sai số tương đối luôn được kiểm soát dưới mức 4.5%.
  • Định hình lộ trình nghiên cứu tiếp theo giai đoạn 2026-2028 hướng tới việc lập trình thư viện phần tử UMAT mã nguồn mở và mở rộng mô hình cho vật liệu xốp biến cứng dị hướng.

Công trình luận văn thạc sĩ của tác giả Dương Thị Hồng Hạnh là tài liệu học thuật và kỹ thuật mẫu mực, đóng góp quan trọng vào kho tàng cơ học biến dạng dẻo và tính toán kết cấu kim loại tại Việt Nam. Các kỹ sư, nhà khoa học và doanh nghiệp hoạt động trong lĩnh vực cơ khí chính xác được khuyến khích nghiên cứu, trích dẫn và ứng dụng ngay các biểu thức hàm chảy của đề tài để nâng cao chất lượng thiết kế sản phẩm và tối ưu hóa quy trình công nghệ tạo hình kim loại.