Một tổng hư hỏng cần giữ được lịch sử phía sau
Khi một báo cáo ghi chỉ số hư hỏng bằng 0,291, người đọc chưa biết mối hàn đã chịu bao nhiêu chu kỳ, chúng nằm ở những mức ứng suất nào và dự báo tiếp theo giả định vận hành ra sao. Nếu chỉ lưu tổng cuối cùng, rất khó cập nhật khi đổi chế độ sản xuất hoặc sửa đường cong tính toán. Bài này xây dựng một sổ tải nhỏ để làm rõ việc cộng, chuyển phổ và tính lại.
Tất cả số liệu là giả định. Mốc so sánh D* = 0,80 được đặt riêng cho bài học, không phải giới hạn nghiệm thu hay tuổi thọ được phê duyệt. Các kết quả là số khối tải tương đương trong mô hình; không dùng trực tiếp để quyết định tiếp tục vận hành một thiết bị thực tế.
Định nghĩa nhất quán trước khi cộng
Gọi S là khoảng ứng suất của một chu kỳ, bằng ứng suất lớn nhất trừ nhỏ nhất, đơn vị MPa. Nó bằng hai lần biên độ khi dùng định nghĩa biên độ thông thường. Quy trình tính mỏi của Ansys tách bước đếm chu kỳ, lập phổ và áp dụng đường cong S-N; các đại lượng khoảng, biên độ và trung bình phải được phân biệt.[1]
Bài chọn đường cong minh họa một nhánh, áp dụng riêng cho các mức từ 40 đến 160 MPa:
N(S) = 2.000.000(80/S)³
N là số chu kỳ trên đường cong tại khoảng S. Đây không phải đường cong phân loại mối hàn trong tiêu chuẩn. Không đặt ngưỡng bỏ qua mức ứng suất nhỏ trong khoảng đang xét. Với n_i chu kỳ tại S_i, quy tắc cộng tuyến tính được viết:
D = Σ n_i/N(S_i)
D không có đơn vị. Nó cũng không phải phần diện tích vết nứt, xác suất hỏng hay tỷ lệ sức bền tĩnh đã mất. Những đại lượng đó cần mô hình khác. TWI trình bày quy tắc này cùng giả thiết rằng đóng góp của một chu kỳ trong phổ được liên hệ với tuổi thọ ở tải biên độ không đổi.[2]
Khối A: cùng một sổ có hai loại tỷ trọng
Giả sử một khối A gồm 100.000 chu kỳ tại 40 MPa, 20.000 tại 80 MPa và 2.000 tại 160 MPa. Các chu kỳ đã được đếm phù hợp với lịch sử ứng suất; bài không thay việc đếm bằng số mẫu đo. Từ đường cong, N lần lượt là 16.000.000; 2.000.000 và 250.000 chu kỳ.
Đóng góp hư hỏng tương ứng bằng 0,00625; 0,01000 và 0,00800. Tổng mỗi khối là D_A = 0,02425. Tổng số chu kỳ là 122.000, nhưng ba nhóm không đóng góp bằng tỷ lệ số lượng của chúng.
Nhóm 160 MPa chỉ chiếm khoảng 1,64% số chu kỳ, lại tạo 32,989691% tổng D. Nhóm 40 MPa chiếm khoảng 81,97% số chu kỳ và tạo 25,773196% tổng D. Vì vậy cả lập luận “chỉ quan tâm mức cao nhất” lẫn “mức thấp nhiều nên chắc chắn chi phối” đều chưa đủ. Cần nhân tác động trên mỗi chu kỳ với số lần xuất hiện.
Trong sổ tính, nên giữ ít nhất S_i, n_i, N_i, n_i/N_i và tỷ trọng đóng góp. Khi đọc một mức tải bất thường, có thể kiểm tra riêng hàng đó thay vì sửa thủ công tổng cuối bảng. Tổng số chu kỳ là một phép đối chiếu dữ liệu, không thay cho tổng hư hỏng.
Mười hai khối quá khứ và một phổ tương lai mới
Sau 12 khối A, D_quá_khứ = 12 × 0,02425 = 0,291. Nay chế độ giả định chuyển sang khối B: 50.000 chu kỳ tại 40 MPa, 30.000 tại 80 MPa và 5.000 tại 160 MPa. Khối B chỉ có 85.000 chu kỳ nhưng:
D_B = 0,003125 + 0,015 + 0,020
Kết quả 0,038125, lớn hơn A. Nếu sử dụng cùng đường cong và cộng tuyến tính, số khối B tương đương để chạm mốc minh họa được tính:
b = (D* - D_quá_khứ)/D_B
Ta được b = 13,350820. Nếu tiếp tục dùng A, kết quả lại là 20,989691 khối. Hai dự báo không mâu thuẫn; chúng trả lời hai kịch bản vận hành khác nhau. Một số “thời gian còn lại” không ghi phổ tương lai sẽ che khuất khác biệt này.
Sau 13 khối B đầy đủ, tổng D = 0,786625; sau 14 khối là 0,824750. Phần 0,350820 khối chỉ là phần tương đương theo tổng hư hỏng. Nếu khối B bắt đầu bằng nhóm 160 MPa, thời điểm vượt mốc trong khối có thể khác trường hợp nhóm đó nằm cuối. Không được tự đổi phần khối thành tỷ lệ thời gian nếu chưa biết trình tự và tốc độ phát sinh chu kỳ.
Độ nhạy với mức ứng suất của phổ mới
Khảo sát thêm trường hợp mọi mức trong B tăng 10%, giả sử đường cong được mở rộng phù hợp đến mức mới để làm phép so sánh toán học. Với số mũ 3, từng đóng góp nhân 1,1³ = 1,331. Khi đó D_B tăng thành 0,050744375, còn số khối tương đương giảm xuống 10,030668.
Số chu kỳ không đổi, ứng suất chỉ tăng 10%, nhưng mức tiêu hao mỗi khối tăng 33,1%. Đây là hệ quả của độ dốc đường cong đã giả định. Nếu thay độ dốc, phải tính lại hệ số; không dùng 33,1% như một quy tắc chung cho mọi chi tiết hàn.
Một lỗi nhập biên độ vào cột yêu cầu khoảng sẽ giảm S một nửa. Trong nhánh m = 3 này, D bị giảm tám lần. Tổng vẫn có thể trông hợp lý vì không có đơn vị để gây cảnh báo. Việc ghi định nghĩa ngay trên cột dữ liệu có giá trị hơn việc chỉ kiểm định dạng số.
Khi đổi cơ sở S-N, tính lại cả phần đã qua
Giả sử sau một đợt rà soát, đường cong minh họa thay hệ số 2.000.000 bằng 1.600.000, giữ mốc ứng suất và số mũ. Mọi N giảm 20%, nên mọi D tăng 25%. Phần quá khứ trở thành 0,36375; mỗi khối B thành 0,04765625.
Số khối còn lại đến cùng mốc 0,80 nay là 9,154098. Nếu chỉ cập nhật B nhưng giữ D_quá_khứ = 0,291, báo cáo sẽ trộn hai thước đo hư hỏng. Phép cộng khi ấy không còn nhất quán, dù từng phép chia được thực hiện đúng.
Lưu histogram hoặc lịch sử chu kỳ gốc cho phép đánh giá lại toàn bộ. Nếu chỉ lưu D, việc đổi số mũ hoặc dùng đường cong nhiều nhánh có thể không khôi phục được kết quả vì các mức ứng suất đã bị gộp mất.
Rời rạc hóa và những ảnh hưởng mô hình không chứa
COMSOL minh họa rằng tính tại tâm các khoảng phân nhóm thay cho ứng suất thực có thể thay đổi tổng hư hỏng, đôi khi làm giảm đóng góp của một chu kỳ.[3] Trong ví dụ này, ba mức là giá trị chính xác được đặt ra, không phải tâm các khoảng rộng. Với dữ liệu thực, nên khảo sát độ nhạy khi thu hẹp khoảng và giữ cùng quy tắc đếm.
Cộng Miner không nhận biết thứ tự tải. TWI đã so sánh phổ tuần tự và ngẫu nhiên trên mẫu hàn; kết quả và giới hạn của nghiên cứu cho thấy không nên xem giả thiết tuyến tính là định luật đúng cho mọi phổ.[2] Bài không lấy hệ số hiệu chỉnh của thí nghiệm đó áp sang thiết bị khác.
Các ảnh hưởng môi trường, trạng thái ứng suất nhiều trục, sai lệch hình học, vết nứt có sẵn và tương tác quá tải chưa được mô tả. Một sổ tính tốt cần ghi rõ phạm vi đó cùng phiên bản đường cong, định nghĩa ứng suất và kỳ dữ liệu. Khi điều kiện thay đổi, phải cập nhật giả thiết vật lý trước rồi mới cập nhật tổng số.
Nguồn tham khảo
[1] 5.4.1. Formulation of Fatigue — Ansys Motion, Release 2025 R1. Truy cập ngày 09/11/2026.
[2] Fatigue Testing of Full-Scale Girth Welded Pipes under Variable Amplitude Loading, OMAE 2012 — Yan-Hui Zhang và Stephen J. Maddox, TWI. Truy cập ngày 09/11/2026.
[3] Cycle Counting in Fatigue Analysis — Benchmark — COMSOL 6.4. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Khi một báo cáo ghi chỉ số hư hỏng bằng 0,291, người đọc chưa biết mối hàn đã chịu bao nhiêu chu kỳ, chúng nằm ở những mức ứng suất nào và dự báo tiếp theo giả định vận hành ra sao. Nếu chỉ lưu tổng cuối cùng, rất khó cập nhật khi đổi chế độ sản xuất hoặc sửa đường cong tính toán. Bài này xây dựng một sổ tải nhỏ để làm rõ việc cộng, chuyển phổ và tính lại.
Tất cả số liệu là giả định. Mốc so sánh D* = 0,80 được đặt riêng cho bài học, không phải giới hạn nghiệm thu hay tuổi thọ được phê duyệt. Các kết quả là số khối tải tương đương trong mô hình; không dùng trực tiếp để quyết định tiếp tục vận hành một thiết bị thực tế.
Định nghĩa nhất quán trước khi cộng
Gọi S là khoảng ứng suất của một chu kỳ, bằng ứng suất lớn nhất trừ nhỏ nhất, đơn vị MPa. Nó bằng hai lần biên độ khi dùng định nghĩa biên độ thông thường. Quy trình tính mỏi của Ansys tách bước đếm chu kỳ, lập phổ và áp dụng đường cong S-N; các đại lượng khoảng, biên độ và trung bình phải được phân biệt.[1]
Bài chọn đường cong minh họa một nhánh, áp dụng riêng cho các mức từ 40 đến 160 MPa:
N(S) = 2.000.000(80/S)³
N là số chu kỳ trên đường cong tại khoảng S. Đây không phải đường cong phân loại mối hàn trong tiêu chuẩn. Không đặt ngưỡng bỏ qua mức ứng suất nhỏ trong khoảng đang xét. Với n_i chu kỳ tại S_i, quy tắc cộng tuyến tính được viết:
D = Σ n_i/N(S_i)
D không có đơn vị. Nó cũng không phải phần diện tích vết nứt, xác suất hỏng hay tỷ lệ sức bền tĩnh đã mất. Những đại lượng đó cần mô hình khác. TWI trình bày quy tắc này cùng giả thiết rằng đóng góp của một chu kỳ trong phổ được liên hệ với tuổi thọ ở tải biên độ không đổi.[2]
Khối A: cùng một sổ có hai loại tỷ trọng
Giả sử một khối A gồm 100.000 chu kỳ tại 40 MPa, 20.000 tại 80 MPa và 2.000 tại 160 MPa. Các chu kỳ đã được đếm phù hợp với lịch sử ứng suất; bài không thay việc đếm bằng số mẫu đo. Từ đường cong, N lần lượt là 16.000.000; 2.000.000 và 250.000 chu kỳ.
Đóng góp hư hỏng tương ứng bằng 0,00625; 0,01000 và 0,00800. Tổng mỗi khối là D_A = 0,02425. Tổng số chu kỳ là 122.000, nhưng ba nhóm không đóng góp bằng tỷ lệ số lượng của chúng.
Nhóm 160 MPa chỉ chiếm khoảng 1,64% số chu kỳ, lại tạo 32,989691% tổng D. Nhóm 40 MPa chiếm khoảng 81,97% số chu kỳ và tạo 25,773196% tổng D. Vì vậy cả lập luận “chỉ quan tâm mức cao nhất” lẫn “mức thấp nhiều nên chắc chắn chi phối” đều chưa đủ. Cần nhân tác động trên mỗi chu kỳ với số lần xuất hiện.
Trong sổ tính, nên giữ ít nhất S_i, n_i, N_i, n_i/N_i và tỷ trọng đóng góp. Khi đọc một mức tải bất thường, có thể kiểm tra riêng hàng đó thay vì sửa thủ công tổng cuối bảng. Tổng số chu kỳ là một phép đối chiếu dữ liệu, không thay cho tổng hư hỏng.
Mười hai khối quá khứ và một phổ tương lai mới
Sau 12 khối A, D_quá_khứ = 12 × 0,02425 = 0,291. Nay chế độ giả định chuyển sang khối B: 50.000 chu kỳ tại 40 MPa, 30.000 tại 80 MPa và 5.000 tại 160 MPa. Khối B chỉ có 85.000 chu kỳ nhưng:
D_B = 0,003125 + 0,015 + 0,020
Kết quả 0,038125, lớn hơn A. Nếu sử dụng cùng đường cong và cộng tuyến tính, số khối B tương đương để chạm mốc minh họa được tính:
b = (D* - D_quá_khứ)/D_B
Ta được b = 13,350820. Nếu tiếp tục dùng A, kết quả lại là 20,989691 khối. Hai dự báo không mâu thuẫn; chúng trả lời hai kịch bản vận hành khác nhau. Một số “thời gian còn lại” không ghi phổ tương lai sẽ che khuất khác biệt này.
Sau 13 khối B đầy đủ, tổng D = 0,786625; sau 14 khối là 0,824750. Phần 0,350820 khối chỉ là phần tương đương theo tổng hư hỏng. Nếu khối B bắt đầu bằng nhóm 160 MPa, thời điểm vượt mốc trong khối có thể khác trường hợp nhóm đó nằm cuối. Không được tự đổi phần khối thành tỷ lệ thời gian nếu chưa biết trình tự và tốc độ phát sinh chu kỳ.
Độ nhạy với mức ứng suất của phổ mới
Khảo sát thêm trường hợp mọi mức trong B tăng 10%, giả sử đường cong được mở rộng phù hợp đến mức mới để làm phép so sánh toán học. Với số mũ 3, từng đóng góp nhân 1,1³ = 1,331. Khi đó D_B tăng thành 0,050744375, còn số khối tương đương giảm xuống 10,030668.
Số chu kỳ không đổi, ứng suất chỉ tăng 10%, nhưng mức tiêu hao mỗi khối tăng 33,1%. Đây là hệ quả của độ dốc đường cong đã giả định. Nếu thay độ dốc, phải tính lại hệ số; không dùng 33,1% như một quy tắc chung cho mọi chi tiết hàn.
Một lỗi nhập biên độ vào cột yêu cầu khoảng sẽ giảm S một nửa. Trong nhánh m = 3 này, D bị giảm tám lần. Tổng vẫn có thể trông hợp lý vì không có đơn vị để gây cảnh báo. Việc ghi định nghĩa ngay trên cột dữ liệu có giá trị hơn việc chỉ kiểm định dạng số.
Khi đổi cơ sở S-N, tính lại cả phần đã qua
Giả sử sau một đợt rà soát, đường cong minh họa thay hệ số 2.000.000 bằng 1.600.000, giữ mốc ứng suất và số mũ. Mọi N giảm 20%, nên mọi D tăng 25%. Phần quá khứ trở thành 0,36375; mỗi khối B thành 0,04765625.
Số khối còn lại đến cùng mốc 0,80 nay là 9,154098. Nếu chỉ cập nhật B nhưng giữ D_quá_khứ = 0,291, báo cáo sẽ trộn hai thước đo hư hỏng. Phép cộng khi ấy không còn nhất quán, dù từng phép chia được thực hiện đúng.
Lưu histogram hoặc lịch sử chu kỳ gốc cho phép đánh giá lại toàn bộ. Nếu chỉ lưu D, việc đổi số mũ hoặc dùng đường cong nhiều nhánh có thể không khôi phục được kết quả vì các mức ứng suất đã bị gộp mất.
Rời rạc hóa và những ảnh hưởng mô hình không chứa
COMSOL minh họa rằng tính tại tâm các khoảng phân nhóm thay cho ứng suất thực có thể thay đổi tổng hư hỏng, đôi khi làm giảm đóng góp của một chu kỳ.[3] Trong ví dụ này, ba mức là giá trị chính xác được đặt ra, không phải tâm các khoảng rộng. Với dữ liệu thực, nên khảo sát độ nhạy khi thu hẹp khoảng và giữ cùng quy tắc đếm.
Cộng Miner không nhận biết thứ tự tải. TWI đã so sánh phổ tuần tự và ngẫu nhiên trên mẫu hàn; kết quả và giới hạn của nghiên cứu cho thấy không nên xem giả thiết tuyến tính là định luật đúng cho mọi phổ.[2] Bài không lấy hệ số hiệu chỉnh của thí nghiệm đó áp sang thiết bị khác.
Các ảnh hưởng môi trường, trạng thái ứng suất nhiều trục, sai lệch hình học, vết nứt có sẵn và tương tác quá tải chưa được mô tả. Một sổ tính tốt cần ghi rõ phạm vi đó cùng phiên bản đường cong, định nghĩa ứng suất và kỳ dữ liệu. Khi điều kiện thay đổi, phải cập nhật giả thiết vật lý trước rồi mới cập nhật tổng số.
Nguồn tham khảo
[1] 5.4.1. Formulation of Fatigue — Ansys Motion, Release 2025 R1. Truy cập ngày 09/11/2026.
[2] Fatigue Testing of Full-Scale Girth Welded Pipes under Variable Amplitude Loading, OMAE 2012 — Yan-Hui Zhang và Stephen J. Maddox, TWI. Truy cập ngày 09/11/2026.
[3] Cycle Counting in Fatigue Analysis — Benchmark — COMSOL 6.4. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
