Bài toán bắt đầu từ một vết nứt đã có
Phân tích tăng trưởng vết nứt trả lời câu hỏi khác với đọc một đường cong S-N cho mẫu chưa được mô tả bằng kích thước nứt ban đầu. Ở đây, a đã tồn tại tại thời điểm bắt đầu. Ta tính số chu kỳ để a đi từ a₀ đến a_f theo một mô hình tốc độ. Thời gian tạo thành vết nứt trước a₀ không nằm trong tích phân.
Bài dựng dữ liệu để kiểm tra công thức, không cung cấp đặc tính vật liệu của một mối hàn thật. Kích thước cuối 6 mm chỉ là điểm kết thúc tính toán, không phải kích thước nứt cho phép hoặc kích thước phá hủy. Muốn xác định các giới hạn đó cần đánh giá riêng hình học, tải và độ dai thích hợp.
Định luật Paris có miền áp dụng
Dạng mô hình được dùng là:
da/dN = C(ΔK)ᵐ
a có đơn vị m; N là số chu kỳ; da/dN có đơn vị m/chu kỳ. ΔK là khoảng hệ số cường độ ứng suất mode I, đơn vị MPa√m. C và m phải được xác định cùng đơn vị, điều kiện thử và phạm vi dữ liệu.
AFGROW trình bày quan hệ Paris cho vùng trung gian của dữ liệu tăng trưởng nứt tại các tỷ số ứng suất xác định, đồng thời cảnh báo việc ngoại suy khỏi phạm vi phù hợp.[1] Không kéo một đường thẳng trên đồ thị logarit đến vùng ngưỡng hoặc gần mất ổn định rồi coi kết quả tự động đúng.
Giả sử riêng bài học rằng mẫu hình học có hệ số không đổi Y = 1,12 trong khoảng xét, chịu tải kéo với σ_min = 20 MPa và σ_max = 100 MPa. Khoảng ứng suất S = 80 MPa, tỷ số R = 0,20. Bỏ qua đóng nứt, tương tác quá tải, biến đổi môi trường và biến đổi hình dạng mặt nứt. Các hệ số giả định C = 2 × 10⁻¹² và m = 3 được coi là phù hợp riêng mô hình này.
Từ ứng suất đến tốc độ biến thiên theo a
Trong lý tưởng hóa đang xét:
ΔK(a) = YS√(πa)
Thế vào quan hệ Paris và đặt q = C(YS√π)³:
da/dN = qa^(3/2)
Ở bộ đơn vị đã nêu, q = 8,010854314 × 10⁻⁶ m^(-1/2)/chu kỳ. Giữ rõ đơn vị giúp kiểm tra vế phải có m/chu kỳ.
Với a₀ = 0,001 m và a_f = 0,006 m, ΔK tăng từ 5,022072 lên 12,301514 MPa√m dù tải ngoài không đổi. Tốc độ tăng từ 2,533255 × 10⁻¹⁰ lên 3,723109 × 10⁻⁹ m/chu kỳ. Vì vậy không thể lấy tốc độ đầu rồi nhân toàn bộ số chu kỳ để ước lượng tuyến tính suốt quá trình.
Ví dụ của COMSOL cũng sử dụng ΔK trong định luật Paris và lưu ý hệ số C phụ thuộc hệ đơn vị của K.[2] Giá trị hệ số được dùng trong bài này được đặt riêng, không lấy từ ví dụ vật liệu của phần mềm.
Tích phân giải tích và kiểm tra đoạn đường
Tách biến cho mô hình m = 3, Y và S không đổi:
N = 2(a₀^(-1/2) - a_f^(-1/2))/q
Kết quả là 4.671.869,287 chu kỳ từ 1 đến 6 mm. Tính riêng đoạn 1 đến 3 mm được 3.336.812,154 chu kỳ; đoạn 3 đến 6 mm được 1.335.057,133 chu kỳ. Tổng hai đoạn bằng kết quả trực tiếp.
Đoạn tăng 2 mm đầu chiếm khoảng 71,42% số chu kỳ, còn đoạn tăng 3 mm sau chỉ chiếm phần còn lại. Tốc độ tăng theo kích thước làm quá trình không tuyến tính. Chia tổng số chu kỳ theo tỷ lệ chiều dài nứt sẽ sai trong chính mô hình này.
Có thể đảo biểu thức để kiểm a tại một số chu kỳ cho trước:
a(N) = [a₀^(-1/2) - qN/2]⁻²
Biểu thức chỉ dùng khi phần trong ngoặc dương và a còn thuộc miền mô hình. Sau 2.000.000 chu kỳ, a khoảng 1,793648 mm. Không diễn giải điểm kỳ dị khi ngoặc bằng không là dự báo vật lý nứt dài vô hạn; các giới hạn hình học và cơ học sẽ xuất hiện trước đó.
Nhập mm có thể làm C đổi theo cách khó đoán
Nếu chuyển đồng thời a sang mm và ΔK sang MPa√mm, thì a_mm = 1.000a_m, còn ΔK_mm = √1.000 ΔK_m. Để tốc độ vật lý không đổi:
C_mm = 1.000C_m/(√1.000)ᵐ
Với m = 3, C_mm = 6,32455532 × 10⁻¹⁴, khi đầu ra là mm/chu kỳ và K nhập MPa√mm. C giảm về trị số dù đơn vị chiều dài đầu ra nhỏ hơn, vì biến K cũng đã đổi và được lũy thừa.
Nếu chỉ đổi a sang mm nhưng vẫn giữ K bằng MPa√m, hệ số đầu ra lại phải nhân 1.000. Hai cách đều có thể nhất quán, nhưng phải ghi bộ đơn vị. Chỉ sao chép C từ một bảng rồi đổi nhãn m thành mm không bảo toàn mô hình.
Kích thước ban đầu ảnh hưởng lớn
Giữ các thông số và đích 6 mm, nếu a₀ = 2 mm, số chu kỳ chỉ còn 2.359.482,501. Nếu a₀ = 0,5 mm, tích phân cho 7.942.078,040 chu kỳ, với điều kiện giả định mô hình Paris còn phù hợp ở vùng nhỏ hơn đó. Không thể lấy kết quả dài hơn làm bằng chứng rằng khả năng phát hiện 0,5 mm đã được chứng minh.
NASA trình bày cách xem kích thước, vị trí nứt và độ phân tán dữ liệu tăng trưởng là các nguồn bất định; phân tích độ nhạy giúp xác định dữ liệu nào cần được làm rõ.[3] Ở đây, các kịch bản chỉ tách một yếu tố, chưa xây dựng phân phối xác suất hoặc độ tin cậy.
Nếu tăng S từ 80 lên 90 MPa, giữ a₀, a_f và C, số chu kỳ nhân (80/90)³, còn 3.281.203,120. Nếu C tăng gấp đôi, số chu kỳ giảm một nửa. Thay m cần xét lại cả ý nghĩa và đơn vị của C, không giữ nguyên trị số C tùy ý rồi gọi đó là khảo sát vật liệu tương đương.
Tích phân số phải được kiểm hội tụ
Có thể tính N bằng diện tích dưới hàm g(a) = 1/(qa^(3/2)). Dùng quy tắc hình thang với các nút 1, 2, 3, 4, 5, 6 mm cho 5.109.896,609 chu kỳ, cao hơn nghiệm đúng khoảng 9,375847%. Sáu nút trông có vẻ đủ trên bảng nhưng hàm cong mạnh gần điểm đầu.
Bài có nghiệm giải tích nên dùng nó làm phép kiểm độc lập cho chương trình. Khi Y phụ thuộc a hoặc dữ liệu tốc độ là bảng, cần giảm bước và so kết quả cho đến khi sai khác số đủ nhỏ đối với mục đích tính. Hội tụ số không khắc phục được dữ liệu vật liệu không phù hợp.
Chẳng hạn, nếu dựng Y(a) = 1,12(1 + 20a), với a tính bằng m, tích phân số cho khoảng 4.062.008,493 chu kỳ. Hệ số hình học tăng theo a làm thời gian lan nứt giảm. Đây chỉ là hàm thử để minh họa việc phải đưa Y(a) vào trong tích phân, không đại diện một hình học mối hàn đã được xác minh.
Kết quả còn thiếu gì trước khi dùng cho công trình?
Cần kiểm hình dạng nứt, kích thước được đo là chiều sâu hay nửa chiều dài, tải thực, điều kiện tuyến tính đàn hồi và dữ liệu tăng trưởng ở môi trường tương ứng. K_max liên quan đến kiểm mất ổn định; ΔK điều khiển mô hình tăng trưởng đang dùng. Hai đại lượng không thay nhau chỉ vì có cùng đơn vị.
Một khoảng tăng trưởng tính được không tự trở thành chu kỳ kiểm tra. Khả năng phát hiện, bất định kích thước, hậu quả hư hỏng và yêu cầu áp dụng còn phải được xét. Giá trị của bài tính là làm rõ đường đi từ dữ liệu đến số chu kỳ, để các giả thiết yếu có thể được nhận ra và thay bằng bằng chứng phù hợp.
Nguồn tham khảo
[1] Section 5.1.2: Fatigue Crack-Growth Rate (FCGR) Descriptions — AFGROW, Handbook for Damage Tolerant Design. Truy cập ngày 09/11/2026.
[2] Single Edge Crack — COMSOL 6.3. Truy cập ngày 09/11/2026.
[3] Fracture Mechanics Reliability, Guideline GD-AP-2304, October 1995 — NASA Johnson Space Center. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Phân tích tăng trưởng vết nứt trả lời câu hỏi khác với đọc một đường cong S-N cho mẫu chưa được mô tả bằng kích thước nứt ban đầu. Ở đây, a đã tồn tại tại thời điểm bắt đầu. Ta tính số chu kỳ để a đi từ a₀ đến a_f theo một mô hình tốc độ. Thời gian tạo thành vết nứt trước a₀ không nằm trong tích phân.
Bài dựng dữ liệu để kiểm tra công thức, không cung cấp đặc tính vật liệu của một mối hàn thật. Kích thước cuối 6 mm chỉ là điểm kết thúc tính toán, không phải kích thước nứt cho phép hoặc kích thước phá hủy. Muốn xác định các giới hạn đó cần đánh giá riêng hình học, tải và độ dai thích hợp.
Định luật Paris có miền áp dụng
Dạng mô hình được dùng là:
da/dN = C(ΔK)ᵐ
a có đơn vị m; N là số chu kỳ; da/dN có đơn vị m/chu kỳ. ΔK là khoảng hệ số cường độ ứng suất mode I, đơn vị MPa√m. C và m phải được xác định cùng đơn vị, điều kiện thử và phạm vi dữ liệu.
AFGROW trình bày quan hệ Paris cho vùng trung gian của dữ liệu tăng trưởng nứt tại các tỷ số ứng suất xác định, đồng thời cảnh báo việc ngoại suy khỏi phạm vi phù hợp.[1] Không kéo một đường thẳng trên đồ thị logarit đến vùng ngưỡng hoặc gần mất ổn định rồi coi kết quả tự động đúng.
Giả sử riêng bài học rằng mẫu hình học có hệ số không đổi Y = 1,12 trong khoảng xét, chịu tải kéo với σ_min = 20 MPa và σ_max = 100 MPa. Khoảng ứng suất S = 80 MPa, tỷ số R = 0,20. Bỏ qua đóng nứt, tương tác quá tải, biến đổi môi trường và biến đổi hình dạng mặt nứt. Các hệ số giả định C = 2 × 10⁻¹² và m = 3 được coi là phù hợp riêng mô hình này.
Từ ứng suất đến tốc độ biến thiên theo a
Trong lý tưởng hóa đang xét:
ΔK(a) = YS√(πa)
Thế vào quan hệ Paris và đặt q = C(YS√π)³:
da/dN = qa^(3/2)
Ở bộ đơn vị đã nêu, q = 8,010854314 × 10⁻⁶ m^(-1/2)/chu kỳ. Giữ rõ đơn vị giúp kiểm tra vế phải có m/chu kỳ.
Với a₀ = 0,001 m và a_f = 0,006 m, ΔK tăng từ 5,022072 lên 12,301514 MPa√m dù tải ngoài không đổi. Tốc độ tăng từ 2,533255 × 10⁻¹⁰ lên 3,723109 × 10⁻⁹ m/chu kỳ. Vì vậy không thể lấy tốc độ đầu rồi nhân toàn bộ số chu kỳ để ước lượng tuyến tính suốt quá trình.
Ví dụ của COMSOL cũng sử dụng ΔK trong định luật Paris và lưu ý hệ số C phụ thuộc hệ đơn vị của K.[2] Giá trị hệ số được dùng trong bài này được đặt riêng, không lấy từ ví dụ vật liệu của phần mềm.
Tích phân giải tích và kiểm tra đoạn đường
Tách biến cho mô hình m = 3, Y và S không đổi:
N = 2(a₀^(-1/2) - a_f^(-1/2))/q
Kết quả là 4.671.869,287 chu kỳ từ 1 đến 6 mm. Tính riêng đoạn 1 đến 3 mm được 3.336.812,154 chu kỳ; đoạn 3 đến 6 mm được 1.335.057,133 chu kỳ. Tổng hai đoạn bằng kết quả trực tiếp.
Đoạn tăng 2 mm đầu chiếm khoảng 71,42% số chu kỳ, còn đoạn tăng 3 mm sau chỉ chiếm phần còn lại. Tốc độ tăng theo kích thước làm quá trình không tuyến tính. Chia tổng số chu kỳ theo tỷ lệ chiều dài nứt sẽ sai trong chính mô hình này.
Có thể đảo biểu thức để kiểm a tại một số chu kỳ cho trước:
a(N) = [a₀^(-1/2) - qN/2]⁻²
Biểu thức chỉ dùng khi phần trong ngoặc dương và a còn thuộc miền mô hình. Sau 2.000.000 chu kỳ, a khoảng 1,793648 mm. Không diễn giải điểm kỳ dị khi ngoặc bằng không là dự báo vật lý nứt dài vô hạn; các giới hạn hình học và cơ học sẽ xuất hiện trước đó.
Nhập mm có thể làm C đổi theo cách khó đoán
Nếu chuyển đồng thời a sang mm và ΔK sang MPa√mm, thì a_mm = 1.000a_m, còn ΔK_mm = √1.000 ΔK_m. Để tốc độ vật lý không đổi:
C_mm = 1.000C_m/(√1.000)ᵐ
Với m = 3, C_mm = 6,32455532 × 10⁻¹⁴, khi đầu ra là mm/chu kỳ và K nhập MPa√mm. C giảm về trị số dù đơn vị chiều dài đầu ra nhỏ hơn, vì biến K cũng đã đổi và được lũy thừa.
Nếu chỉ đổi a sang mm nhưng vẫn giữ K bằng MPa√m, hệ số đầu ra lại phải nhân 1.000. Hai cách đều có thể nhất quán, nhưng phải ghi bộ đơn vị. Chỉ sao chép C từ một bảng rồi đổi nhãn m thành mm không bảo toàn mô hình.
Kích thước ban đầu ảnh hưởng lớn
Giữ các thông số và đích 6 mm, nếu a₀ = 2 mm, số chu kỳ chỉ còn 2.359.482,501. Nếu a₀ = 0,5 mm, tích phân cho 7.942.078,040 chu kỳ, với điều kiện giả định mô hình Paris còn phù hợp ở vùng nhỏ hơn đó. Không thể lấy kết quả dài hơn làm bằng chứng rằng khả năng phát hiện 0,5 mm đã được chứng minh.
NASA trình bày cách xem kích thước, vị trí nứt và độ phân tán dữ liệu tăng trưởng là các nguồn bất định; phân tích độ nhạy giúp xác định dữ liệu nào cần được làm rõ.[3] Ở đây, các kịch bản chỉ tách một yếu tố, chưa xây dựng phân phối xác suất hoặc độ tin cậy.
Nếu tăng S từ 80 lên 90 MPa, giữ a₀, a_f và C, số chu kỳ nhân (80/90)³, còn 3.281.203,120. Nếu C tăng gấp đôi, số chu kỳ giảm một nửa. Thay m cần xét lại cả ý nghĩa và đơn vị của C, không giữ nguyên trị số C tùy ý rồi gọi đó là khảo sát vật liệu tương đương.
Tích phân số phải được kiểm hội tụ
Có thể tính N bằng diện tích dưới hàm g(a) = 1/(qa^(3/2)). Dùng quy tắc hình thang với các nút 1, 2, 3, 4, 5, 6 mm cho 5.109.896,609 chu kỳ, cao hơn nghiệm đúng khoảng 9,375847%. Sáu nút trông có vẻ đủ trên bảng nhưng hàm cong mạnh gần điểm đầu.
Bài có nghiệm giải tích nên dùng nó làm phép kiểm độc lập cho chương trình. Khi Y phụ thuộc a hoặc dữ liệu tốc độ là bảng, cần giảm bước và so kết quả cho đến khi sai khác số đủ nhỏ đối với mục đích tính. Hội tụ số không khắc phục được dữ liệu vật liệu không phù hợp.
Chẳng hạn, nếu dựng Y(a) = 1,12(1 + 20a), với a tính bằng m, tích phân số cho khoảng 4.062.008,493 chu kỳ. Hệ số hình học tăng theo a làm thời gian lan nứt giảm. Đây chỉ là hàm thử để minh họa việc phải đưa Y(a) vào trong tích phân, không đại diện một hình học mối hàn đã được xác minh.
Kết quả còn thiếu gì trước khi dùng cho công trình?
Cần kiểm hình dạng nứt, kích thước được đo là chiều sâu hay nửa chiều dài, tải thực, điều kiện tuyến tính đàn hồi và dữ liệu tăng trưởng ở môi trường tương ứng. K_max liên quan đến kiểm mất ổn định; ΔK điều khiển mô hình tăng trưởng đang dùng. Hai đại lượng không thay nhau chỉ vì có cùng đơn vị.
Một khoảng tăng trưởng tính được không tự trở thành chu kỳ kiểm tra. Khả năng phát hiện, bất định kích thước, hậu quả hư hỏng và yêu cầu áp dụng còn phải được xét. Giá trị của bài tính là làm rõ đường đi từ dữ liệu đến số chu kỳ, để các giả thiết yếu có thể được nhận ra và thay bằng bằng chứng phù hợp.
Nguồn tham khảo
[1] Section 5.1.2: Fatigue Crack-Growth Rate (FCGR) Descriptions — AFGROW, Handbook for Damage Tolerant Design. Truy cập ngày 09/11/2026.
[2] Single Edge Crack — COMSOL 6.3. Truy cập ngày 09/11/2026.
[3] Fracture Mechanics Reliability, Guideline GD-AP-2304, October 1995 — NASA Johnson Space Center. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
