Phân tích localized failure - Ứng dụng vào thép, bê tông và bê tông cốt thép

Phân tích hư hỏng cục bộ là kỹ thuật đánh giá sự suy giảm khả năng chịu tải tại vùng cụ thể. Bản chất của hiện tượng này liên quan trực tiếp đến tập trung ứng suất. Khi tải trọng tác dụng, ứng suất không phân bố đều trong cấu kiện. Các vùng có khuyết tật hoặc hình học bất thường chịu ứng suất lớn hơn. Tập trung ứng suất là nguyên nhân chính gây ra localized failure. Trong vật liệu dẻo như thép, biến dạng cục bộ tạo thành vùng co thắt cổ (necking). Ở vật liệu giòn như bê tông, hư hỏng cục bộ biểu hiện qua vết nứt tập trung. Phương pháp phân tích cần xem xét cả yếu tố nhiệt và cơ. Nhiệt độ thay đổi làm biến đổi tính chất vật liệu. Điều này ảnh hưởng trực tiếp đến quá trình strain localization. Mô hình hóa chính xác yêu cầu lý thuyết cơ học vết nứt tiên tiến.

Cơ học vết nứt cung cấp nền tảng lý thuyết cho phân tích hư hỏng cục bộ. Lý thuyết này giải thích cơ chế lan truyền vết nứt trong vật liệu. Fracture mechanics liên hệ giữa năng lượng giải phóng và sức bền vật liệu. Trong bài toán nhiệt-cơ, crack tip chịu cả ứng suất nhiệt và cơ. Cơ học vết nứt mở rộng cho phép đánh giá tương tác này. Griffith và Irwin là hai nhà khoa học tiên phong trong lĩnh vực này. Nghiên cứu của Ngo Van Minh kế thừa và phát triển truyền thống đó. Mô hình cohesive zone thường được sử dụng cho localized failure. Mô hình này giả định hư hỏng xảy ra trong vùng liên kết mỏng. Phương pháp có ưu điểm giảm phụ thuộc vào lưới phần tử. Điều này rất quan trọng cho bài toán strain localization thực tế.

Ba loại vật liệu chính được nghiên cứu: thép, bê tông và bê tông cốt thép. Thép là vật liệu dẻo, có khả năng chịu biến dạng cục bộ lớn trước khi phá hủy. Shear bands trong thép thường xuất hiện ở nhiệt độ cao hoặc tốc độ biến dạng lớn. Bê tông là vật liệu giòn, hư hỏng cục bộ biểu hiện qua vết nứt phức tạp. Vết nứt trong bê tông phát triển theo hướng ứng suất chính lớn nhất. Bê tông cốt thép kết hợp ưu điểm của cả hai vật liệu. Cốt thép ngăn cản sự lan truyền vết nứt trong bê tông. Tuy nhiên, localized failure vẫn xảy ra tại vùng tiếp xúc giữa hai vật liệu. Bài toán coupled thermo-mechanics làm phức tạp thêm quá trình hư hỏng. Mỗi loại vật liệu đòi hỏi mô hình phân tích khác nhau. Luận án đã xây dựng mô hình riêng cho từng loại vật liệu.

Mô hình cohesive zone là công cụ hiệu quả cho phân tích hư hỏng cục bộ. Mô hình giả định vùng liên kết mỏng tồn tại phía trước vết nứt. Quan hệ constitutive liên hệ giữa ứng suất bề mặt và độ mở vết nứt. Đường cong cohesive thường gồm hai phần: tăng tải và giảm tải. Phần tăng tải tương ứng với quá trình tích lũy hư hỏng vi mô. Phần giảm tải thể hiện sự phát triển hư hỏng macro. Diện tích dưới đường cong cohesive bằng năng lượng giải phóng G. Đối với bài toán nhiệt-cơ, tham số cohesive phụ thuộc nhiệt độ. Nhiệt độ cao thường giảm cường độ và năng lượng cohesive. Mô hình này phù hợp cho thép và bê tông cốt thép. Ưu điểm lớn nhất là không cần vết nứt trước tồn tại. Localized failure tự phát triển từ vùng ứng suất tập trung lớn nhất.

Phương pháp gradient bổ sung chiều dài nội tại vào mô hình constitutive. Chiều dài nội tại kiểm soát độ rộng của vùng strain localization. Điều này ngăn ngừa sự phụ thuộc vào lưới phần tử. Mô hình gradient bậc hai yêu cầu điều kiện biên bổ sung cho biến nội tại. Phase-field approach mô tả hư hỏng qua biến liên tục phi từ 0 đến 1. Giá trị 0 tương ứng vật liệu nguyên vẹn, giá trị 1 tương ứng hư hỏng hoàn toàn. Phương trình lan truyền phase-field tương tự phương trình khuếch tán. Dải trượt được mô tả tự nhiên trong khung phase-field. Phương pháp này đặc biệt hiệu quả cho bài toán nucleation vết nứt mới. Shear bands phát triển tự nhiên mà không cần điều kiện ban đầu đặc biệt. Trong bài toán nhiệt-cơ, phase-field coupled mô tả tương tác nhiệt-hư hỏng.

Tập trung ứng suất là vấn đề số học quan trọng trong phân tích FEM. Phần tử tiêu chuẩn không thể hiện chính xác gradient ứng suất lớn. Phương pháp enrichment bổ sung hàm hình dạng đặc biệt tại vùng singularity. XFEM là phương pháp phổ biến nhất cho bài toán crack propagation. XFEM không yêu cầu lưới phù hợp với hình học vết nứt. Phương pháp này kết hợp tốt với mô hình cohesive zone. Đối với shear bands, cần sử dụng phần tử đặc biệt có chiều dài nội tại. Phần tử F-bar cải thiện hiệu suất tính toán cho biến dạng cục bộ lớn. Tích hợp phần tử giảm bậc giúp tránh hiện tượng locking. Bài toán coupled thermo-mechanics đòi hỏi phần tử phù hợp cho cả trường nhiệt. Luận án đã kiểm tra và so sánh nhiều loại phần tử khác nhau.

Mô hình constitutive của thép phải phản ánh đúng tính chất theo nhiệt độ. Ở nhiệt độ thường, thép tuân theo luật elastoplastic với hardening. Khi nhiệt độ tăng, quá trình creep trở nên quan trọng. Mô hình tổng hợp kết hợp plasticity và creep cho localized failure. Tham số yield stress giảm theo hàm mũ của nhiệt độ. Module đàn hồi cũng giảm đáng kể ở nhiệt độ cao. Hiệu ứng creep dẫn đến biến dạng cục bộ tăng theo thời gian. Shear bands hình thành khi biến dạng creep vượt ngưỡng critical. Mô hình cần tham số nhiệt độ để tính đúng strain localization. Dải trượt trong thép nóng thường rộng hơn thép ở nhiệt độ thường. Điều này ảnh hưởng trực tiếp đến cơ chế phá hủy cuối cùng. Mô hình đã được kiểm chứng bằng thực nghiệm tại ENS Cachan.

Shear bands trong thép là kết quả của quá trình strain localization. Nhiệt độ cao thúc đẩy quá trình này diễn ra nhanh hơn. Cơ chế hình thành gồm ba giai đoạn: nucleation, propagation và coalescence. Giai đoạn nucleation bắt đầu khi tập trung ứng suất vượt ngưỡng local strength. Shear bands lan truyền dọc theo mặt phẳng ứng suất cắt lớn nhất. Khi nhiều shear bands gặp nhau, quá trình coalescence xảy ra. Vết nứt cuối cùng hình thành tại vùng giao nhau. Nhiệt độ gradient làm phức tạp thêm quá trình phát triển. Vùng nóng hơn phá hủy trước, tạo điều kiện cho strain localization. Mô hình coupled thermo-mechanics mô tả chính xác chuỗi sự kiện này. Kết quả mô phỏng phù hợp với hình ảnh thực nghiệm. Phương pháp có ý nghĩa lớn cho đánh giá an toàn kết cấu thép.

Kết quả mô phỏng đã được kiểm chứng bằng thí nghiệm tiêu chuẩn. Thí nghiệm kéo đơn giản ở nhiệt độ tăng dần được thực hiện đầu tiên. Kết quả cho thấy localized failure xảy ra ở nhiệt độ critical khoảng 500-700°C. Thí nghiệm bending beams xác nhận sự phát triển shear bands trong thép nóng. Hình ảnh DIC (Digital Image Correlation) ghi lại trường biến dạng thực tế. So sánh cho thấy mô hình dự đoán đúng vị trí strain localization. Sai số về tải trọng phá hủy nhỏ hơn 10% so với thực nghiệm. Thí nghiệm trên kết cấu thép thực tế cũng được thực hiện. Kết quả xác nhận tính khả thi của phương pháp phân tích hư hỏng cục bộ. Mô hình coupled thermo-mechanics cải thiện đáng kể độ chính xác. Nghiên cứu đã được công bố trên tạp chí quốc tế uy tín.

Mô hình smeared crack phân bố hư hỏng trên vùng rộng của phần tử. Phương pháp này dễ triển khai trong code phần tử hữu hạn. Tuy nhiên, kết quả phụ thuộc nhiều vào kích thước lưới. Mô hình discrete crack mô tả vết nứt như bề mặt gián đoạn rõ ràng. Phương pháp này chính xác hơn về mặt vật lý. XFEM kết hợp ưu điểm của cả hai phương pháp. Vết nứt lan truyền tự do trong lưới phần tử không cần remeshing. Đối với bài toán nhiệt-cơ, cần mô hình kết hợp cả hai cách tiếp cận. Vùng gần vết nứt sử dụng discrete crack, vùng xa sử dụng smeared crack. Chiều rộng dải hư hỏng được kiểm soát bởi fracture energy. Mô hình đã được áp dụng thành công cho bê tông thường và bê tông cường độ cao. Kết quả phù hợp với thực nghiệm bẻ cong ba điểm và kéo trực tiếp.

Nhiệt độ gây ra nhiều hiệu ứng có hại cho bê tông. Đầu tiên, mất nước dẫn đến co ngót nhiệt và tạo vết nứt vi mô. Thứ hai, sự khác biệt hệ số giãn nở nhiệt giữa đá xi măng và xương liệu. Điều này tạo tập trung ứng suất tại giao diện giữa hai thành phần. Ứng suất nhiệt nội tại thúc đẩy quá trình strain localization. Ở nhiệt độ trên 300°C, bê tông bắt đầu suy giảm cường độ rõ rệt. Localized failure xảy ra sớm hơn và ở tải trọng thấp hơn. Mô hình coupled thermo-mechanics tính đến tất cả hiệu ứng này. Phương trình nhiệt truyền kết hợp với mô hình constitutive hư hỏng. Kết quả cho thấy nhiệt độ gradient ảnh hưởng lớn đến hướng vết nứt. Vết nứt có xu hướng phát triển từ vùng nóng sang vùng lạnh. Phân tích này quan trọng cho thiết kế kết cấu bê tông chịu nhiệt.

Bê tông cốt thép (BTCT) là vật liệu composite phức tạp. Fracture mechanics của BTCT phải xem xét tương tác giữa bê tông và thép. Cốt thép ngăn cản sự lan truyền vết nứt trong bê tông. Tuy nhiên, nhiệt độ cao làm suy giảm cả hai thành phần. Ở nhiệt độ trên 400°C, liên kết giữa bê tông và cốt thép bị phá hủy. Hiện tượng này gọi là bond degradation và thúc đẩy localized failure. Mô hình interface element mô tả giao diện bê tông-cốt thép. Tham số cohesive của giao diện phụ thuộc nhiệt độ. Shear bands xuất hiện tại vùng liên kết bị suy yếu. Mô hình coupled thermo-mechanics cho BTCT rất phức tạp. Cần mô hình hóa đồng thời: hư hỏng bê tông, plasticity thép, phá hủy liên kết. Luận án đã xây dựng framework tổng hợp cho tất cả hiện tượng này.

Thuật toán Newton-Raphson là phương pháp chính cho bài toán phi tuyến. Mỗi bước lặp yêu cầu giải hệ tuyến tính kích thước lớn. Ma trận tangent phải được cập nhật sau mỗi bước lặp. Chiến lược line search cải thiện hội tụ cho bài toán localized failure. Quasi-Newton methods giảm chi phí tính toán cho bài toán lớn. Tích hợp thời gian sử dụng phương pháp Newmark cho bài toán động. Đối với bài toán tĩnh, phương pháp arc-length kiểm soát quá trình tải. Arc-length đặc biệt quan trọng khi strain localization xảy ra. Phần mềm MATLAB và ABAQUS được sử dụng cho implementation. Parallel computing giảm thời gian tính toán cho mô hình 3D. Kiểm tra hội tụ dựa trên residue của cả trường nhiệt và cơ. Chiến lược adaptive time stepping cải thiện hiệu quả tính toán.

Validation là bước quan trọng để kiểm chứng mô hình toán học. Benchmark đầu tiên là bài toán nhiệt truyền đơn giản trong thanh chịu kéo. Kết quả phân tích và giải tích khớp nhau với sai số <1%. Benchmark thứ hai là bài toán bê tông chịu tải nhiệt tiêu chuẩn. So sánh với kết quả thí nghiệm cho thấy độ chính xác cao. Benchmark thứ ba là bài toán thép ở nhiệt độ cao chịu tải trọng. Mô hình dự đoán đúng thời điểm và vị trí localized failure. Bài toán BTCT phức tạp hơn và cần nhiều benchmark hơn. Thí nghiệm fire resistance được sử dụng cho validation cuối cùng. Kết quả cho thấy mô hình coupled thermo-mechanics cải thiện 30% độ chính xác. So với mô hình uncoupled, mô hình coupled dự đoán tốt hơn. Nghi究 đã được xuất bản trên nhiều tạp chí ISI uy tín.

Vật liệu composite ngày càng phổ biến trong kỹ thuật. CFRP, GFRP và vật liệu nano có tính chất phức tạp hơn. Localized failure trong composite diễn ra theo nhiều chế độ: delamination, fiber breakage, matrix cracking. Mô hình coupled thermo-mechanics cần mở rộng cho vật liệu đa pha. Strain localization trong composite ảnh hưởng bởi hướng sợi và stacking sequence. Shear bands trong composite khác biệt so với vật liệu isotropic. Nghiên cứu hiện tại tập trung vào mô hình multiscale. Mô hình liên kết quy mô nano, micro, meso và macro. Phương pháp phase-field rất phù hợp cho bài toán composite. Fracture mechanics mở rộng cho vật liệu orthotropic đang được phát triển. Kết quả ban đầu rất hứa hẹn cho ứng dụng thực tế.

Ứng dụng chính của phân tích hư hỏng cục bộ là thiết kế công trình an toàn. Kết cấu chịu nhiệt như lò phản ứng hạt nhân đòi hỏi phân tích chính xác. Đường hầm, cầu và tòa nhà cần đánh giá khả năng chịu lửa. Tiêu chuẩn thiết kế hiện tại dựa trên mô hình đơn giản hóa. Mô hình localized failure cung cấp kết quả chi tiết và chính xác hơn. Phần mềm commercial như ABAQUS đã tích hợp nhiều mô hình hư hỏng. Tuy nhiên, cần chuyên gia hiểu rõ lý thuyết để sử dụng hiệu quả. Đánh giá an toàn kết cấu hiện hữu là ứng dụng quan trọng khác. Nhiều cầu và tòa nhà cũ cần kiểm tra khả năng chịu tải hiện tại. Phân tích hư hỏng cục bộ giúp xác định vùng yếu và biện pháp gia cố. Ứng dụng trong ngành dầu khí cho kết cấu ngoài khơi cũng rất lớn.

Nhiều thách thức lớn vẫn tồn tại trong lĩnh vực localized failure. Thứ nhất là vấn đề quy mô: liên kết mô hình vi mô và vĩ mô. Thứ hai là tính toán hiệu quả cho mô hình 3D kích thước lớn. Thứ ba là dữ liệu thực nghiệm đầy đủ cho validation mô hình. Machine learning hứa hẹn cải thiện tốc độ tính toán đáng kể. Deep learning có thể học mối quan hệ phi tuyến phức tạp. Digital twin kết hợp mô hình localized failure với giám sát thời gian thực. Internet of Things cung cấp dữ liệu nhiệt và biến dạng liên tục. Tích hợp tất cả tạo thành hệ thống đánh giá thông minh. Tuy nhiên, cần nhiều nghiên cứu hơn để hiện thực hóa tầm nhìn này. Cộng đồng nghiên cứu quốc tế đang hợp tác chặt chẽ. Các hội nghị quốc tế về fracture mechanics tổ chức thường niên. Tương lai rất hứa hẹn cho lĩnh vực phân tích hư hỏng cục bộ.