Ứng suất danh nghĩa chưa mô tả được đầu nứt
Hai tấm chịu cùng ứng suất kéo có thể có trạng thái rất khác nếu một tấm chứa vết nứt sắc. Cơ học phá hủy đưa kích thước và hình học vết nứt vào bài toán, thay vì chỉ so ứng suất trung bình với giới hạn chảy. Bài này xét một mô hình đàn hồi tuyến tính để đọc đúng hệ số cường độ ứng suất K, không thiết lập điều kiện chấp nhận khuyết tật cho một thiết bị dầu khí.
Toàn bộ trị số vật liệu, tải và kích thước dưới đây là giả định. Vết nứt xuyên chiều dày nằm giữa tấm rất rộng, tổng chiều dài bằng 2a. Ứng suất kéo xa vết nứt vuông góc với mặt nứt; hai mặt được phép mở. Không có tiếp xúc mặt nứt, tải hỗn hợp hoặc biến dạng dẻo lớn trong mô hình ban đầu.
K mô tả biên độ trường gần đầu nứt
Trong chế độ mở I, dạng dẫn đầu của ứng suất gần đầu nứt tỷ lệ K_I/√(2πr), nhân với hàm góc thích hợp. r là khoảng cách tới đầu nứt. AFGROW giải thích rằng tải và hình học đi vào biên độ K, còn mô hình đàn hồi chỉ thích hợp khi vùng biến dạng dẻo cục bộ đủ nhỏ trong miền đàn hồi bao quanh.[2]
K_I có đơn vị ứng suất nhân căn chiều dài, thường dùng MPa√m. Nó không phải một ứng suất theo MPa, cũng không phải hệ số tập trung ứng suất K_t không thứ nguyên. Giá trị ứng suất lý tưởng tăng vô hạn khi r tiến về không là đặc điểm toán học của đầu nứt sắc; không có nghĩa vật liệu thật chịu được ứng suất vô hạn.
Với hình học giữa tấm vô hạn đã chọn, tài liệu MIT cho:[1]
K_I = σ√(πa)
σ là ứng suất xa vết nứt, a là nửa chiều dài toàn vết nứt. Với hình học khác, thường viết K_I = Yσ√(πa), nhưng Y phụ thuộc cấu tạo, cách định nghĩa a và cách xác định ứng suất. Không tự đặt Y = 1 cho một vết nứt mép, nứt bán elip hoặc nứt ở chân mối hàn.
Tính một trường hợp và kiểm hai lỗi đơn vị
Chọn σ = 120 MPa và tổng chiều dài vết nứt 2a = 8 mm. Khi ấy a = 4 mm = 0,004 m:
K_I = 120√(π×0,004)
K_I = 13,451979 MPa√m
Nếu nhập toàn chiều dài 8 mm vào chỗ a, kết quả thành 19,023971 MPa√m, cao hơn √2 lần. Đây là lỗi định nghĩa hình học. Nếu giữ a = 4 nhưng vẫn ghi đơn vị đầu ra MPa√m, kết quả thành 425,388924, cao hơn √1.000 lần. Thực chất phép tính đó đang dùng MPa√mm và chưa đổi đơn vị.
Một hồ sơ kiểm tra cần giữ bản vẽ thể hiện đầu nứt nào tới đầu nứt nào, mặt cắt nào dùng và độ sâu nào được đo. Cột “kích thước nứt” không đủ nếu không biết đó là a, 2a, chiều sâu hay chiều dài chỉ thị trên bề mặt. Cùng con số 8 mm có thể thuộc các bài toán hoàn toàn khác nhau.
Độ nhạy tải và kích thước không ngang nhau
Giữ hình học tương tự và Y không đổi, tăng σ 10% làm K tăng đúng 10%. Tăng a 10% chỉ làm K tăng √1,10 lần, khoảng 4,881%. Đối với thay đổi nhỏ:
dK/K = dσ/σ + 0,5 da/a
Nếu Y cũng thay đổi theo a thì phải thêm dY/Y. Vì vậy, hệ số một phần hai không phải độ nhạy chung cho mọi vết nứt trong một tấm hữu hạn. Khi nứt tiến sát mép, việc tiếp tục dùng Y hằng có thể bỏ qua tác động hình học đang tăng nhanh.
Giả định σ nằm trong khoảng 114–126 MPa, còn a nằm trong khoảng 3,5–4,5 mm. Do K tăng theo cả hai biến trong mô hình này, hai góc biên cho K_min = 11,954015 và K_max = 14,981377 MPa√m. Khoảng này là kết quả của các giới hạn đầu vào đã đặt; nó chưa có ý nghĩa xác suất nếu chưa biết các giới hạn được xác định thế nào.
Một cách mô tả khác dùng độ không đảm bảo chuẩn độc lập u_σ = 6 MPa và u_a = 0,30 mm. Truyền tuyến tính theo đạo hàm cho:
(u_K/K)² = (u_σ/σ)² + (u_a/2a)²
u_K bằng 0,840749 MPa√m, tức 6,25% K gốc. Hai cách mô tả này không thay thế nhau: khoảng biên xét tổ hợp cực trị, còn độ không đảm bảo chuẩn dựa trên một mô hình thống kê. Không gọi 13,45±0,84 là khoảng bảo đảm tuyệt đối.
So với một mức tới hạn giả định
Đặt riêng một mức K_c = 40 MPa√m để khảo sát toán học. Không gắn trị số này với một mác thép hoặc gọi nó là kết quả thử K_Ic. Với Y = 1 và σ = 120 MPa, nghiệm K_I = K_c là:
a_c = (K_c/σ)²/π
a_c = 35,367765 mm
Đó là nửa chiều dài, nên tổng chiều dài tương ứng bằng 70,735530 mm. Nếu giữ a = 4 mm thay vì giữ σ, ứng suất tại cùng mức K_c là 356,824823 MPa. Hai phép tính đảo này dùng cùng phương trình nhưng trả lời hai câu hỏi khác nhau.
Khi a_c tăng tới hàng chục milimét, phải quay lại kiểm bề rộng thật và phần vật liệu còn lại. Một công thức giả định tấm vô hạn không tự chứng minh điều kiện hình học vẫn đúng tại nghiệm cuối. Tương tự, việc K hiện tại nhỏ hơn K_c không chứng minh đã xét chảy dẻo toàn tiết diện, tăng trưởng mỏi, môi trường hoặc tải bất thường.
Ước lượng vùng dẻo để nhìn giới hạn mô hình
Đặt ứng suất đàn hồi dẫn đầu ở phía trước đầu nứt bằng một giới hạn chảy giả định σ_y = 500 MPa. Phép ước lượng bậc đầu theo ứng suất phẳng trong tài liệu MIT cho:[1]
r_y = (K_I/σ_y)²/(2π)
Với K gốc, r_y = 0,1152 mm, bằng 2,88% a. Đây là một thang chiều dài từ trường đàn hồi, không phải hình dạng vùng dẻo đầy đủ và không phải tiêu chí chứng nhận phép thử. Ràng buộc theo chiều dày, tái phân bố ứng suất và hóa bền không nằm trong phép đặt bằng đơn giản này.
Nếu tăng σ lên 240 MPa, K gấp đôi còn r_y tăng bốn lần lên 0,4608 mm, tức 11,52% a. Trị số ứng suất danh nghĩa vẫn nhỏ hơn 500 MPa nhưng vùng phi tuyến đã lớn hơn tương đối. Quan sát này giải thích vì sao chỉ so σ với σ_y chưa đủ kiểm tra điều kiện của cơ học phá hủy đàn hồi.
Phân biệt K_Ic với một số độ dai bất kỳ
Trang công khai ASTM E399-24 mô tả K_Ic trong điều kiện đàn hồi tuyến tính, ràng buộc kéo cao và vùng dẻo nhỏ; tính hợp lệ của phép thử cần được đánh giá theo phương pháp tương ứng. Trang này cũng nêu giới hạn đối với phá hủy cleavage và biến dạng dẻo đáng kể.[3] Không lấy một kết quả Charpy theo joule để thế trực tiếp vào công thức K theo MPa√m.
Khi dùng dữ liệu, cần giữ nhiệt độ, môi trường, tốc độ thử, hướng lấy mẫu và trạng thái vật liệu. Một giá trị trong bảng không tự đại diện cho kim loại mối hàn, vùng ảnh hưởng nhiệt và vật liệu nền cùng lúc. Bài toán ví dụ cố ý giữ K_c là một giả thiết để không biến nó thành khuyến nghị vật liệu.
Chồng tải cần cùng một trường hợp hình học
Nếu có thêm một ứng suất kéo đều 40 MPa, cùng phương và cùng hình học đàn hồi đang xét, tổng K bằng 160√(πa) = 17,935972 MPa√m. Kết quả này chỉ đúng vì hai tải có cùng quan hệ với trường đầu nứt. Một phân bố ứng suất dư tự cân bằng hoặc một trường uốn không thể luôn thay bằng giá trị đỉnh rồi cộng vào σ đều.
Báo cáo có thể kiểm tra lại phải nêu hình học nứt, tải, biểu thức Y, đơn vị, cách xác định dữ liệu độ dai và các cơ chế chưa xét. K là cầu nối hữu ích giữa tải và vết nứt trong phạm vi mô hình của nó. Một phép so duy nhất giữa hai số K không thay toàn bộ đánh giá tính toàn vẹn của cấu kiện.
Nguồn tham khảo
[1] Introduction to Fracture Mechanics — David Roylance, June 14, 2001 — Massachusetts Institute of Technology — OpenCourseWare. Truy cập ngày 09/11/2026.
[2] DTD Handbook 2.2.1 — Stress Intensity Factor: What It Is — AFGROW — Damage Tolerant Design Handbook. Truy cập ngày 09/11/2026.
[3] ASTM E399-24 — Standard Test Method for Linear-Elastic Plane-Strain Fracture Toughness of Metallic Materials, public scope and significance — ASTM International. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Hai tấm chịu cùng ứng suất kéo có thể có trạng thái rất khác nếu một tấm chứa vết nứt sắc. Cơ học phá hủy đưa kích thước và hình học vết nứt vào bài toán, thay vì chỉ so ứng suất trung bình với giới hạn chảy. Bài này xét một mô hình đàn hồi tuyến tính để đọc đúng hệ số cường độ ứng suất K, không thiết lập điều kiện chấp nhận khuyết tật cho một thiết bị dầu khí.
Toàn bộ trị số vật liệu, tải và kích thước dưới đây là giả định. Vết nứt xuyên chiều dày nằm giữa tấm rất rộng, tổng chiều dài bằng 2a. Ứng suất kéo xa vết nứt vuông góc với mặt nứt; hai mặt được phép mở. Không có tiếp xúc mặt nứt, tải hỗn hợp hoặc biến dạng dẻo lớn trong mô hình ban đầu.
K mô tả biên độ trường gần đầu nứt
Trong chế độ mở I, dạng dẫn đầu của ứng suất gần đầu nứt tỷ lệ K_I/√(2πr), nhân với hàm góc thích hợp. r là khoảng cách tới đầu nứt. AFGROW giải thích rằng tải và hình học đi vào biên độ K, còn mô hình đàn hồi chỉ thích hợp khi vùng biến dạng dẻo cục bộ đủ nhỏ trong miền đàn hồi bao quanh.[2]
K_I có đơn vị ứng suất nhân căn chiều dài, thường dùng MPa√m. Nó không phải một ứng suất theo MPa, cũng không phải hệ số tập trung ứng suất K_t không thứ nguyên. Giá trị ứng suất lý tưởng tăng vô hạn khi r tiến về không là đặc điểm toán học của đầu nứt sắc; không có nghĩa vật liệu thật chịu được ứng suất vô hạn.
Với hình học giữa tấm vô hạn đã chọn, tài liệu MIT cho:[1]
K_I = σ√(πa)
σ là ứng suất xa vết nứt, a là nửa chiều dài toàn vết nứt. Với hình học khác, thường viết K_I = Yσ√(πa), nhưng Y phụ thuộc cấu tạo, cách định nghĩa a và cách xác định ứng suất. Không tự đặt Y = 1 cho một vết nứt mép, nứt bán elip hoặc nứt ở chân mối hàn.
Tính một trường hợp và kiểm hai lỗi đơn vị
Chọn σ = 120 MPa và tổng chiều dài vết nứt 2a = 8 mm. Khi ấy a = 4 mm = 0,004 m:
K_I = 120√(π×0,004)
K_I = 13,451979 MPa√m
Nếu nhập toàn chiều dài 8 mm vào chỗ a, kết quả thành 19,023971 MPa√m, cao hơn √2 lần. Đây là lỗi định nghĩa hình học. Nếu giữ a = 4 nhưng vẫn ghi đơn vị đầu ra MPa√m, kết quả thành 425,388924, cao hơn √1.000 lần. Thực chất phép tính đó đang dùng MPa√mm và chưa đổi đơn vị.
Một hồ sơ kiểm tra cần giữ bản vẽ thể hiện đầu nứt nào tới đầu nứt nào, mặt cắt nào dùng và độ sâu nào được đo. Cột “kích thước nứt” không đủ nếu không biết đó là a, 2a, chiều sâu hay chiều dài chỉ thị trên bề mặt. Cùng con số 8 mm có thể thuộc các bài toán hoàn toàn khác nhau.
Độ nhạy tải và kích thước không ngang nhau
Giữ hình học tương tự và Y không đổi, tăng σ 10% làm K tăng đúng 10%. Tăng a 10% chỉ làm K tăng √1,10 lần, khoảng 4,881%. Đối với thay đổi nhỏ:
dK/K = dσ/σ + 0,5 da/a
Nếu Y cũng thay đổi theo a thì phải thêm dY/Y. Vì vậy, hệ số một phần hai không phải độ nhạy chung cho mọi vết nứt trong một tấm hữu hạn. Khi nứt tiến sát mép, việc tiếp tục dùng Y hằng có thể bỏ qua tác động hình học đang tăng nhanh.
Giả định σ nằm trong khoảng 114–126 MPa, còn a nằm trong khoảng 3,5–4,5 mm. Do K tăng theo cả hai biến trong mô hình này, hai góc biên cho K_min = 11,954015 và K_max = 14,981377 MPa√m. Khoảng này là kết quả của các giới hạn đầu vào đã đặt; nó chưa có ý nghĩa xác suất nếu chưa biết các giới hạn được xác định thế nào.
Một cách mô tả khác dùng độ không đảm bảo chuẩn độc lập u_σ = 6 MPa và u_a = 0,30 mm. Truyền tuyến tính theo đạo hàm cho:
(u_K/K)² = (u_σ/σ)² + (u_a/2a)²
u_K bằng 0,840749 MPa√m, tức 6,25% K gốc. Hai cách mô tả này không thay thế nhau: khoảng biên xét tổ hợp cực trị, còn độ không đảm bảo chuẩn dựa trên một mô hình thống kê. Không gọi 13,45±0,84 là khoảng bảo đảm tuyệt đối.
So với một mức tới hạn giả định
Đặt riêng một mức K_c = 40 MPa√m để khảo sát toán học. Không gắn trị số này với một mác thép hoặc gọi nó là kết quả thử K_Ic. Với Y = 1 và σ = 120 MPa, nghiệm K_I = K_c là:
a_c = (K_c/σ)²/π
a_c = 35,367765 mm
Đó là nửa chiều dài, nên tổng chiều dài tương ứng bằng 70,735530 mm. Nếu giữ a = 4 mm thay vì giữ σ, ứng suất tại cùng mức K_c là 356,824823 MPa. Hai phép tính đảo này dùng cùng phương trình nhưng trả lời hai câu hỏi khác nhau.
Khi a_c tăng tới hàng chục milimét, phải quay lại kiểm bề rộng thật và phần vật liệu còn lại. Một công thức giả định tấm vô hạn không tự chứng minh điều kiện hình học vẫn đúng tại nghiệm cuối. Tương tự, việc K hiện tại nhỏ hơn K_c không chứng minh đã xét chảy dẻo toàn tiết diện, tăng trưởng mỏi, môi trường hoặc tải bất thường.
Ước lượng vùng dẻo để nhìn giới hạn mô hình
Đặt ứng suất đàn hồi dẫn đầu ở phía trước đầu nứt bằng một giới hạn chảy giả định σ_y = 500 MPa. Phép ước lượng bậc đầu theo ứng suất phẳng trong tài liệu MIT cho:[1]
r_y = (K_I/σ_y)²/(2π)
Với K gốc, r_y = 0,1152 mm, bằng 2,88% a. Đây là một thang chiều dài từ trường đàn hồi, không phải hình dạng vùng dẻo đầy đủ và không phải tiêu chí chứng nhận phép thử. Ràng buộc theo chiều dày, tái phân bố ứng suất và hóa bền không nằm trong phép đặt bằng đơn giản này.
Nếu tăng σ lên 240 MPa, K gấp đôi còn r_y tăng bốn lần lên 0,4608 mm, tức 11,52% a. Trị số ứng suất danh nghĩa vẫn nhỏ hơn 500 MPa nhưng vùng phi tuyến đã lớn hơn tương đối. Quan sát này giải thích vì sao chỉ so σ với σ_y chưa đủ kiểm tra điều kiện của cơ học phá hủy đàn hồi.
Phân biệt K_Ic với một số độ dai bất kỳ
Trang công khai ASTM E399-24 mô tả K_Ic trong điều kiện đàn hồi tuyến tính, ràng buộc kéo cao và vùng dẻo nhỏ; tính hợp lệ của phép thử cần được đánh giá theo phương pháp tương ứng. Trang này cũng nêu giới hạn đối với phá hủy cleavage và biến dạng dẻo đáng kể.[3] Không lấy một kết quả Charpy theo joule để thế trực tiếp vào công thức K theo MPa√m.
Khi dùng dữ liệu, cần giữ nhiệt độ, môi trường, tốc độ thử, hướng lấy mẫu và trạng thái vật liệu. Một giá trị trong bảng không tự đại diện cho kim loại mối hàn, vùng ảnh hưởng nhiệt và vật liệu nền cùng lúc. Bài toán ví dụ cố ý giữ K_c là một giả thiết để không biến nó thành khuyến nghị vật liệu.
Chồng tải cần cùng một trường hợp hình học
Nếu có thêm một ứng suất kéo đều 40 MPa, cùng phương và cùng hình học đàn hồi đang xét, tổng K bằng 160√(πa) = 17,935972 MPa√m. Kết quả này chỉ đúng vì hai tải có cùng quan hệ với trường đầu nứt. Một phân bố ứng suất dư tự cân bằng hoặc một trường uốn không thể luôn thay bằng giá trị đỉnh rồi cộng vào σ đều.
Báo cáo có thể kiểm tra lại phải nêu hình học nứt, tải, biểu thức Y, đơn vị, cách xác định dữ liệu độ dai và các cơ chế chưa xét. K là cầu nối hữu ích giữa tải và vết nứt trong phạm vi mô hình của nó. Một phép so duy nhất giữa hai số K không thay toàn bộ đánh giá tính toàn vẹn của cấu kiện.
Nguồn tham khảo
[1] Introduction to Fracture Mechanics — David Roylance, June 14, 2001 — Massachusetts Institute of Technology — OpenCourseWare. Truy cập ngày 09/11/2026.
[2] DTD Handbook 2.2.1 — Stress Intensity Factor: What It Is — AFGROW — Damage Tolerant Design Handbook. Truy cập ngày 09/11/2026.
[3] ASTM E399-24 — Standard Test Method for Linear-Elastic Plane-Strain Fracture Toughness of Metallic Materials, public scope and significance — ASTM International. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
