Yếu tố cường độ ứng suất động (dynamic stress intensity factor - DSIF) là thông số cơ học định lượng mức độ tập trung của trường ứng suất và biến dạng đàn hồi tại vùng lân cận đỉnh vết nứt khi kết cấu chịu tải trọng thay đổi nhanh theo thời gian hoặc khi vết nứt lan truyền với vận tốc cao. Trong cơ học rạn nứt đàn hồi tuyến tính động, đại lượng này mở rộng khái niệm hệ số cường độ ứng suất tĩnh truyền thống bằng cách tính toán trực tiếp ảnh hưởng của quán tính khối lượng và sự lan truyền của các sóng ứng suất đàn hồi trong vật liệu. Bài viết trình bày chi tiết cơ sở lý thuyết, các phương trình trường tiệm cận, phương pháp đo thực nghiệm quang học và ý nghĩa của yếu tố cường độ ứng suất động trong phân tích an toàn kết cấu chịu va đập.
Bản chất vật lý và trường ứng suất động tiệm cận
Trong cơ học rạn nứt đàn hồi tuyến tính cổ điển, trường ứng suất tĩnh tại lân cận đỉnh vết nứt phẳng được biểu diễn qua nghiệm tiệm cận của Williams với đặc trưng kỳ dị tỷ lệ nghịch với căn bậc hai của khoảng cách tính từ đỉnh nứt. Khi vật thể chịu tải trọng động biến đổi nhanh theo thời gian như xung sóng nổ, va đập cơ học hoặc tải trọng địa chấn, gia tốc của các phần tử vật chất không còn có thể bỏ qua. Phương trình cân bằng tĩnh được thay thế bằng phương trình chuyển động động lực học liên tục dạng Navier-Cauchy:
\sigma_{ij,j} = ho rac{\partial^2 u_i}{\partial t^2}
Trong đó là tensor ứng suất, là vector dịch chuyển và là khối lượng riêng của vật liệu. Sự xuất hiện của số hạng quán tính ở vế phải làm xuất hiện hiện tượng truyền sóng đàn hồi trong môi trường rắn, bao gồm sóng dãn dài dọc, sóng biến dạng trượt và sóng bề mặt Rayleigh.
Đối với một vết nứt đứng yên chịu tải trọng động cục bộ, trường ứng suất tiệm cận tại lân cận đỉnh nứt vẫn giữ nguyên tính chất kỳ dị nhưng độ lớn của trường ứng suất lúc này được khuếch đại bởi yếu tố cường độ ứng suất động phụ thuộc thời gian :
\sigma_{ij}(r, heta, t) = rac{K_I(t)}{\sqrt{2\pi r}} f_{ij}( heta) + \dots
Theo chuyên khảo nền tảng của Freund (1990), yếu tố cường độ ứng suất động không đồng nhất với giá trị tính toán từ các phương trình tĩnh. Khi sóng ứng suất va đập truyền đến đỉnh vết nứt, sự phản xạ và giao thoa của các sóng đàn hồi kết hợp với quán tính của khối vật liệu xung quanh tạo ra hiện tượng quá độ mạnh mẽ. Phân tích giải tích của Freund (1990) chỉ ra rằng giá trị đỉnh tức thời của yếu tố cường độ ứng suất động có thể vượt giá trị tĩnh tương ứng từ 20 đến 40% trước khi dao động tắt dần và hội tụ về nghiệm tĩnh bền vững khi năng lượng tiêu tán.
Vết nứt chuyển động nhanh và hàm vận tốc phổ quát
Trường hợp phức tạp và thách thức hơn trong cơ học phá hủy động xuất hiện khi chính vết nứt đang chuyển động với vận tốc cao . Khi đó, trường ứng suất cục bộ phụ thuộc đồng thời vào lịch sử tải trọng, thời gian lan truyền và vận tốc tức thời của đỉnh nứt. Nghiên cứu giải tích kinh điển của Baker (1962) về trường ứng suất động sinh ra bởi một vết nứt chuyển động đã đặt nền móng cho việc phân tích trường tiệm cận chuyển động.
Dựa trên lý thuyết đàn hồi động học, Freund (1990) đã chứng minh rằng yếu tố cường độ ứng suất động của một vết nứt đang lan truyền trong trường tải trọng bất kỳ có thể phân tách thành tích của một hàm vận tốc phổ quát và yếu tố cường độ ứng suất tương đương của vết nứt đứng yên tại cùng chiều dài:
Hàm vận tốc phổ quát có biểu thức giải tích xấp xỉ phụ thuộc vào vận tốc nứt và vận tốc sóng bề mặt Rayleigh trong vật liệu:
k(v) pprox rac{1 - v/c_R}{\sqrt{1 - v/c_d}}
Nghiệm giải tích của Baker (1962) và Freund (1990) cho thấy khi vận tốc lan truyền của vết nứt tiến dần đến vận tốc sóng Rayleigh , hàm vận tốc phổ quát sẽ tiệm cận về 0. Điều này hàm ý rằng dòng năng lượng cơ học truyền vào đỉnh vết nứt bị triệt tiêu hoàn toàn tại vận tốc sóng Rayleigh. Do đó, trong môi trường đàn hồi tuyến tính đồng nhất, vận tốc sóng Rayleigh đóng vai trò như một giới hạn vận tốc vật lý tối đa mà một vết nứt kéo mở mode I không bao giờ có thể vượt qua.
Tiêu chuẩn phá hủy động và độ dai phá hủy
Để đánh giá khả năng chống lại sự rạn nứt của vật liệu dưới điều kiện tải trọng va đập, các nhà kỹ thuật phân biệt ba đại lượng độ dai phá hủy cơ bản tương ứng với các trạng thái động lực học khác nhau:
| Đại lượng | Ký hiệu | Bản chất vật lý | Đặc điểm phụ thuộc |
|---|---|---|---|
| Độ dai phá hủy tĩnh | Khả năng chống nứt dưới tải trọng tĩnh hoặc tăng chậm | Phụ thuộc vào nhiệt độ và cấu trúc vi mô vật liệu | |
| Độ dai khởi nứt động | Giá trị ngưỡng để vết nứt đứng yên bắt đầu phát triển dưới tải động | Phụ thuộc mạnh vào tốc độ tăng ứng suất dK/dt và nhiệt độ | |
| Độ dai nứt lan truyền | Mức cản trở năng lượng khi vết nứt đang di chuyển ở vận tốc v | Hàm phi tuyến tăng nhanh theo vận tốc chuyển động của đỉnh nứt | |
| Độ dai dừng nứt | Giá trị yếu tố cường độ ứng suất tại thời điểm vết nứt dừng lại | Đặc trưng cho giới hạn an toàn ngăn chặn nứt thảm khốc |
Tiêu chuẩn khởi nứt động phát biểu rằng vết nứt đứng yên sẽ bắt đầu lan truyền tại thời điểm khi yếu tố cường độ ứng suất động tức thời đạt đến giá trị độ dai phá hủy khởi nứt động của vật liệu:
Khi vết nứt đã chuyển động, sự cân bằng động học đòi hỏi yếu tố cường độ ứng suất động tác dụng phải liên tục bằng độ dai nứt lan truyền: . Nếu tốc độ cấp năng lượng đàn hồi từ trường biến dạng xung quanh vượt quá năng lượng tiêu tán cực đại của quá trình rạn nứt, vết nứt sẽ đạt đến vận tốc giới hạn (thường bằng 50 đến 60% vận tốc sóng Rayleigh theo phân tích thực nghiệm của Freund, 1990) và xảy ra hiện tượng phân nhánh vết nứt (crack branching) để phân tán năng lượng dư thừa.
Kỹ thuật thực nghiệm đo trường quang học tốc độ cao
Việc đo đạc thực nghiệm yếu tố cường độ ứng suất động gặp nhiều thách thức to lớn do các quá trình biến dạng và lan truyền nứt diễn ra trong khoảng thời gian cực ngắn tính bằng micro-giây. Trước đây, các nhà nghiên cứu chủ yếu dựa vào phương pháp quang đàn hồi động (dynamic photoelasticity) hoặc phương pháp đường tụ quang (method of caustics), vốn đòi hỏi vật liệu mẫu phải trong suốt hoặc có bề mặt phản xạ gương quang học hoàn hảo.
Sự ra đời của phương pháp tương quan ảnh số (Digital Image Correlation - DIC) kết hợp với máy ảnh kỹ thuật số siêu cao tốc đã tạo ra bước đột phá trong thực nghiệm cơ học phá hủy. Trong công trình tiên phong của Kirugulige và cs. (2007), hệ thống chụp ảnh tốc độ hàng trăm nghìn khung hình mỗi giây được sử dụng để theo dõi sự dịch chuyển của các vân đốm ngẫu nhiên trên bề mặt mẫu vật liệu chịu tải va đập. Bằng cách so sánh sự thay đổi hình học giữa các khung hình liên tiếp, thuật toán DIC tính toán chính xác trường dịch chuyển toàn phần hai chiều xung quanh đỉnh vết nứt theo thời gian.
Từ dữ liệu trường dịch chuyển đo được trong thực nghiệm, yếu tố cường độ ứng suất động được trích xuất thông qua việc làm khớp trường dịch chuyển thực nghiệm với các khai triển tiệm cận động học bậc cao. Nghiên cứu thực nghiệm của Kojima và cs. (2021) trên tạp chí Engineering Fracture Mechanics đã ứng dụng thành công kỹ thuật DIC tốc độ cao để đánh giá độ dai phá hủy động và xác định chính xác lịch sử thay đổi của yếu tố cường độ ứng suất động trên các liên kết mặt phân giới vật liệu chịu tải biến dạng cắt tốc độ cao.
Ý nghĩa kỹ thuật và ứng dụng trong thiết kế công trình
Phân tích yếu tố cường độ ứng suất động đóng vai trò sống còn trong việc đảm bảo độ tin cậy và an toàn của các kết cấu kỹ thuật quan trọng chịu tác động của tải trọng xung động. Trong ngành công nghiệp năng lượng hạt nhân, các thùng lò phản ứng bằng thép phải được đánh giá độ dai dừng nứt để ngăn ngừa nguy cơ nứt vỡ thảm khốc khi hệ thống làm mát khẩn cấp kích hoạt gây sốc nhiệt động lực học.
Trong kỹ thuật hàng không vũ trụ và quốc phòng, việc hiểu rõ hành vi của yếu tố cường độ ứng suất động giúp các kỹ sư thiết kế tối ưu hóa cấu trúc tấm chịu va chạm mảnh đạn hoặc va quẹt chim bay ở vận tốc siêu âm. Bằng cách chèn các lớp vật liệu cản sóng hoặc bố trí các vùng biến tính vi mô để kích hoạt cơ chế dừng nứt, kết cấu có thể hấp thụ năng lượng va chạm mà không dẫn đến sự sụp đổ dây chuyền của toàn bộ hệ thống.