1. Đại lượng cần cộng là tổn thương
Một nút jacket có thể trải qua nhiều tháng biển êm và một số đợt tải lớn. Lấy trung bình các khoảng ứng suất rồi đưa vào đường S–N thường làm mất tác động phi tuyến của tải. Quy tắc Palmgren–Miner gán cho nhóm i phần tổn thương nᵢ/Nᵢ, trong đó nᵢ là số chu kỳ tác dụng và Nᵢ là số chu kỳ theo đường cong ở cùng mức ứng suất. Các phần này được cộng thành D không thứ nguyên. [1]
D = Σ(nᵢ/Nᵢ)
Giả thuyết tuyến tính nằm ở phép cộng tổn thương, không nằm ở quan hệ giữa ứng suất và tuổi thọ. Nếu N = C/Sᵐ thì đóng góp của một chu kỳ tỷ lệ với Sᵐ. Cần giữ định nghĩa S nhất quán: khoảng ứng suất, biên độ hoặc đại lượng khác phải đúng với đường cong đang sử dụng. Ví dụ dưới đây dùng khoảng ứng suất, chưa có hiệu chỉnh ứng suất trung bình và không dùng hệ số từ một phân loại hàn thực tế.
2. Dữ liệu giả định cho một năm hoạt động
Chọn đường một nhánh giả định N = C/S³, với C = 2 × 10¹³ khi S nhập bằng MPa. Một năm quy ước có ba trạng thái: A ở 40 MPa, chu kỳ 8 s, tổng 4.000 giờ; B ở 70 MPa, chu kỳ 10 s, tổng 2.000 giờ; C ở 110 MPa, chu kỳ 12 s, tổng 400 giờ. Trong từng trạng thái, ta thay lịch sử bằng tải tuần hoàn đều với đúng khoảng và chu kỳ đã cho.
Tổng thời gian có tải là 6.400 giờ. Phần 2.360 giờ còn lại được giả định không có chu kỳ ứng suất tại điểm xét. Đây là giả thiết của bài toán minh họa, không có nghĩa kết cấu ngừng sản xuất thì tự động hết tải sóng. Nếu trong thời gian dừng vẫn có dao động, phải bổ sung các nhóm đó vào ngân sách năm.
nᵢ = 3.600hᵢ/Tᵢ
hᵢ tính bằng giờ, Tᵢ bằng giây. Ta nhận nA = 1.800.000; nB = 720.000; nC = 120.000 chu kỳ. Tổng cộng 2.640.000 chu kỳ. Không lấy tổng giờ chia một chu kỳ trung bình đơn giản vì các nhóm có tần suất xuất hiện khác nhau.
3. Tính và truy vết từng phần
Tại 40 MPa, NA = 312.500.000 chu kỳ, DA = 0,005760. Tại 70 MPa, NB = 58.309.038 chu kỳ, DB = 0,012348. Tại 110 MPa, NC = 15.026.296 chu kỳ, DC = 0,007986. Tổng tổn thương năm bằng:
D năm = 0,005760 + 0,012348 + 0,007986
D năm = 0,026094
Nhóm C chỉ chiếm 6,25% thời gian có tải nhưng tạo 30,60% tổn thương. Tuy nhiên nhóm B vẫn đóng góp lớn nhất, khoảng 47,32%. Chỉ tìm trạng thái có ứng suất lớn nhất sẽ bỏ qua việc một mức tải vừa phải xảy ra thường xuyên có thể chi phối tổng tổn thương.
OrcaFlex mô tả quy trình tính tổn thương cho một lần xuất hiện của trường hợp tải, nhân theo thời gian phơi nhiễm, rồi cộng qua các trường hợp tại cùng điểm mỏi. [2] Áp dụng cách kiểm này, một mô phỏng B dài 3 giờ có 1.080 chu kỳ và D mẫu = 0,000018522. Nhân với 2.000/3 cho DB = 0,012348. Không vừa nhân số chu kỳ thành cả năm vừa nhân thêm hệ số phơi nhiễm, vì như vậy thời gian được tính hai lần.
4. Ứng suất tương đương phải bảo toàn tổn thương
Với đường một nhánh và cùng điều kiện hiệu chỉnh, có thể tìm khoảng tương đương Sₑ gây cùng D trong tổng n chu kỳ:
Sₑ = [Σ(nᵢSᵢ³)/Σnᵢ]^(1/3)
Kết quả Sₑ = 58,2535 MPa. Trung bình khoảng ứng suất theo số chu kỳ chỉ bằng 51,3636 MPa. Nếu dùng giá trị trung bình này với toàn bộ 2.640.000 chu kỳ, D còn 0,0178872, thấp hơn kết quả đúng khoảng 31,45%. Sai lệch là hậu quả trực tiếp của hàm lũy thừa, không phải do làm tròn số chu kỳ.
Khái niệm tương đương còn phụ thuộc số chu kỳ tham chiếu. Nếu chọn một triệu chu kỳ thay vì 2,64 triệu, Sₑ phải tăng để giữ D. Do đó báo cáo chỉ ghi “ứng suất tương đương” mà không ghi số chu kỳ và đường cong là chưa đủ. Với đường hai nhánh, không được dùng một số m chung cho mọi nhóm vượt qua điểm gãy; khi cần, giải phương trình n/N(Sₑ) = D trên đúng nhánh của Sₑ.
5. Độ nhạy với tải và thời gian
Tăng riêng khoảng của nhóm C thêm 10%, từ 110 lên 121 MPa, giữ nguyên số chu kỳ, làm DC tăng theo 1,1³ = 1,331. Tổng D năm tăng lên 0,028737366, tức tăng 10,13%. Mức thay đổi toàn bộ nhỏ hơn 33,1% vì chỉ một phần phổ tải bị thay đổi.
Một kịch bản khác giữ tổng 6.400 giờ: chuyển 100 giờ từ nhóm A sang nhóm C. Phần bỏ khỏi A là 100 × 3.600/8 chu kỳ tại 40 MPa; phần thêm vào C là 100 × 3.600/12 chu kỳ tại 110 MPa. D năm mới bằng 0,0279465, tăng khoảng 7,10%. Không cộng thêm 100 giờ C mà giữ nguyên A nếu mục tiêu là khảo sát thay đổi tỷ lệ thời gian trong cùng một năm.
So sánh đóng góp mỗi giờ cũng có ích: A tạo 0,00000144 mỗi giờ, B tạo 0,000006174, C tạo 0,000019965. Một giờ C trong mô hình này gây tổn thương gấp 13,8646 lần một giờ A. Tỷ số khác (110/40)³ vì chu kỳ 12 s và 8 s khác nhau. Dùng mỗi tỷ số ứng suất sẽ bỏ mất tốc độ tích lũy chu kỳ.
6. Từ tổn thương năm đến thời gian và giới hạn
Nếu phổ năm lặp y hệt, mô hình cho D sau 20 năm bằng 0,52188. Thời gian toán học để D đạt 1 là 38,323 năm. Đây là hệ quả của dữ liệu giả định và phép cộng tuyến tính, chưa phải tuổi thọ thiết kế hay thời gian được phép tiếp tục khai thác. Mức tổn thương cho phép, hệ số thiết kế và điều kiện kiểm tra phải lấy từ cơ sở thiết kế áp dụng cho công trình.
Nếu một cấu kiện đã có D lịch sử = 0,18 theo cùng mô hình và định nghĩa, phần còn lại đến mốc 1 là 0,82; chia cho D năm mới được khoảng 31,42 năm. Phép cộng này chỉ có ý nghĩa khi các giai đoạn dùng đường cong, điểm xét và quy tắc xử lý ứng suất tương thích. Sau thay thế hoặc sửa chữa, không tự đặt lịch sử bằng không khi chưa đánh giá phần chi tiết còn lại.
COMSOL nêu rõ Palmgren–Miner không xét thứ tự tải: đưa nhóm tải lớn lên đầu hay cuối vẫn cho cùng đóng góp. [3] Trong ví dụ, đổi thứ tự A–B–C thành C–A–B giữ nguyên D. Kết quả bằng nhau phản ánh giả thuyết mô hình, không chứng minh hai lịch sử thực tế gây trạng thái vật liệu giống nhau. Khi biến dạng dẻo, phát triển vết nứt hoặc tương tác các mức tải trở nên quan trọng, cần phương pháp phù hợp với cơ chế đó.
Bảng tổng hợp nên giữ cả nᵢ, Sᵢ, Nᵢ, Dᵢ, thời gian đại diện và phiên bản đường cong. Nhờ vậy người kiểm tra xác định được sai lệch đến từ phổ biển, đáp ứng kết cấu hay quy đổi thời gian. Chỉ lưu một con số D cuối cùng sẽ làm mất khả năng giải thích vì sao tổng tăng sau khi cập nhật dữ liệu.
Nguồn tham khảo
[1] Calculating Damage with Miner’s Rule — Siemens Digital Industries Software. Truy cập ngày 09/11/2026.
[2] OrcaFlex Help — Fatigue analysis: How damage is calculated — Orcina. Truy cập ngày 09/11/2026.
[3] Fatigue Module 6.4 — Cumulative Damage Model — COMSOL. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Một nút jacket có thể trải qua nhiều tháng biển êm và một số đợt tải lớn. Lấy trung bình các khoảng ứng suất rồi đưa vào đường S–N thường làm mất tác động phi tuyến của tải. Quy tắc Palmgren–Miner gán cho nhóm i phần tổn thương nᵢ/Nᵢ, trong đó nᵢ là số chu kỳ tác dụng và Nᵢ là số chu kỳ theo đường cong ở cùng mức ứng suất. Các phần này được cộng thành D không thứ nguyên. [1]
D = Σ(nᵢ/Nᵢ)
Giả thuyết tuyến tính nằm ở phép cộng tổn thương, không nằm ở quan hệ giữa ứng suất và tuổi thọ. Nếu N = C/Sᵐ thì đóng góp của một chu kỳ tỷ lệ với Sᵐ. Cần giữ định nghĩa S nhất quán: khoảng ứng suất, biên độ hoặc đại lượng khác phải đúng với đường cong đang sử dụng. Ví dụ dưới đây dùng khoảng ứng suất, chưa có hiệu chỉnh ứng suất trung bình và không dùng hệ số từ một phân loại hàn thực tế.
2. Dữ liệu giả định cho một năm hoạt động
Chọn đường một nhánh giả định N = C/S³, với C = 2 × 10¹³ khi S nhập bằng MPa. Một năm quy ước có ba trạng thái: A ở 40 MPa, chu kỳ 8 s, tổng 4.000 giờ; B ở 70 MPa, chu kỳ 10 s, tổng 2.000 giờ; C ở 110 MPa, chu kỳ 12 s, tổng 400 giờ. Trong từng trạng thái, ta thay lịch sử bằng tải tuần hoàn đều với đúng khoảng và chu kỳ đã cho.
Tổng thời gian có tải là 6.400 giờ. Phần 2.360 giờ còn lại được giả định không có chu kỳ ứng suất tại điểm xét. Đây là giả thiết của bài toán minh họa, không có nghĩa kết cấu ngừng sản xuất thì tự động hết tải sóng. Nếu trong thời gian dừng vẫn có dao động, phải bổ sung các nhóm đó vào ngân sách năm.
nᵢ = 3.600hᵢ/Tᵢ
hᵢ tính bằng giờ, Tᵢ bằng giây. Ta nhận nA = 1.800.000; nB = 720.000; nC = 120.000 chu kỳ. Tổng cộng 2.640.000 chu kỳ. Không lấy tổng giờ chia một chu kỳ trung bình đơn giản vì các nhóm có tần suất xuất hiện khác nhau.
3. Tính và truy vết từng phần
Tại 40 MPa, NA = 312.500.000 chu kỳ, DA = 0,005760. Tại 70 MPa, NB = 58.309.038 chu kỳ, DB = 0,012348. Tại 110 MPa, NC = 15.026.296 chu kỳ, DC = 0,007986. Tổng tổn thương năm bằng:
D năm = 0,005760 + 0,012348 + 0,007986
D năm = 0,026094
Nhóm C chỉ chiếm 6,25% thời gian có tải nhưng tạo 30,60% tổn thương. Tuy nhiên nhóm B vẫn đóng góp lớn nhất, khoảng 47,32%. Chỉ tìm trạng thái có ứng suất lớn nhất sẽ bỏ qua việc một mức tải vừa phải xảy ra thường xuyên có thể chi phối tổng tổn thương.
OrcaFlex mô tả quy trình tính tổn thương cho một lần xuất hiện của trường hợp tải, nhân theo thời gian phơi nhiễm, rồi cộng qua các trường hợp tại cùng điểm mỏi. [2] Áp dụng cách kiểm này, một mô phỏng B dài 3 giờ có 1.080 chu kỳ và D mẫu = 0,000018522. Nhân với 2.000/3 cho DB = 0,012348. Không vừa nhân số chu kỳ thành cả năm vừa nhân thêm hệ số phơi nhiễm, vì như vậy thời gian được tính hai lần.
4. Ứng suất tương đương phải bảo toàn tổn thương
Với đường một nhánh và cùng điều kiện hiệu chỉnh, có thể tìm khoảng tương đương Sₑ gây cùng D trong tổng n chu kỳ:
Sₑ = [Σ(nᵢSᵢ³)/Σnᵢ]^(1/3)
Kết quả Sₑ = 58,2535 MPa. Trung bình khoảng ứng suất theo số chu kỳ chỉ bằng 51,3636 MPa. Nếu dùng giá trị trung bình này với toàn bộ 2.640.000 chu kỳ, D còn 0,0178872, thấp hơn kết quả đúng khoảng 31,45%. Sai lệch là hậu quả trực tiếp của hàm lũy thừa, không phải do làm tròn số chu kỳ.
Khái niệm tương đương còn phụ thuộc số chu kỳ tham chiếu. Nếu chọn một triệu chu kỳ thay vì 2,64 triệu, Sₑ phải tăng để giữ D. Do đó báo cáo chỉ ghi “ứng suất tương đương” mà không ghi số chu kỳ và đường cong là chưa đủ. Với đường hai nhánh, không được dùng một số m chung cho mọi nhóm vượt qua điểm gãy; khi cần, giải phương trình n/N(Sₑ) = D trên đúng nhánh của Sₑ.
5. Độ nhạy với tải và thời gian
Tăng riêng khoảng của nhóm C thêm 10%, từ 110 lên 121 MPa, giữ nguyên số chu kỳ, làm DC tăng theo 1,1³ = 1,331. Tổng D năm tăng lên 0,028737366, tức tăng 10,13%. Mức thay đổi toàn bộ nhỏ hơn 33,1% vì chỉ một phần phổ tải bị thay đổi.
Một kịch bản khác giữ tổng 6.400 giờ: chuyển 100 giờ từ nhóm A sang nhóm C. Phần bỏ khỏi A là 100 × 3.600/8 chu kỳ tại 40 MPa; phần thêm vào C là 100 × 3.600/12 chu kỳ tại 110 MPa. D năm mới bằng 0,0279465, tăng khoảng 7,10%. Không cộng thêm 100 giờ C mà giữ nguyên A nếu mục tiêu là khảo sát thay đổi tỷ lệ thời gian trong cùng một năm.
So sánh đóng góp mỗi giờ cũng có ích: A tạo 0,00000144 mỗi giờ, B tạo 0,000006174, C tạo 0,000019965. Một giờ C trong mô hình này gây tổn thương gấp 13,8646 lần một giờ A. Tỷ số khác (110/40)³ vì chu kỳ 12 s và 8 s khác nhau. Dùng mỗi tỷ số ứng suất sẽ bỏ mất tốc độ tích lũy chu kỳ.
6. Từ tổn thương năm đến thời gian và giới hạn
Nếu phổ năm lặp y hệt, mô hình cho D sau 20 năm bằng 0,52188. Thời gian toán học để D đạt 1 là 38,323 năm. Đây là hệ quả của dữ liệu giả định và phép cộng tuyến tính, chưa phải tuổi thọ thiết kế hay thời gian được phép tiếp tục khai thác. Mức tổn thương cho phép, hệ số thiết kế và điều kiện kiểm tra phải lấy từ cơ sở thiết kế áp dụng cho công trình.
Nếu một cấu kiện đã có D lịch sử = 0,18 theo cùng mô hình và định nghĩa, phần còn lại đến mốc 1 là 0,82; chia cho D năm mới được khoảng 31,42 năm. Phép cộng này chỉ có ý nghĩa khi các giai đoạn dùng đường cong, điểm xét và quy tắc xử lý ứng suất tương thích. Sau thay thế hoặc sửa chữa, không tự đặt lịch sử bằng không khi chưa đánh giá phần chi tiết còn lại.
COMSOL nêu rõ Palmgren–Miner không xét thứ tự tải: đưa nhóm tải lớn lên đầu hay cuối vẫn cho cùng đóng góp. [3] Trong ví dụ, đổi thứ tự A–B–C thành C–A–B giữ nguyên D. Kết quả bằng nhau phản ánh giả thuyết mô hình, không chứng minh hai lịch sử thực tế gây trạng thái vật liệu giống nhau. Khi biến dạng dẻo, phát triển vết nứt hoặc tương tác các mức tải trở nên quan trọng, cần phương pháp phù hợp với cơ chế đó.
Bảng tổng hợp nên giữ cả nᵢ, Sᵢ, Nᵢ, Dᵢ, thời gian đại diện và phiên bản đường cong. Nhờ vậy người kiểm tra xác định được sai lệch đến từ phổ biển, đáp ứng kết cấu hay quy đổi thời gian. Chỉ lưu một con số D cuối cùng sẽ làm mất khả năng giải thích vì sao tổng tăng sau khi cập nhật dữ liệu.
Nguồn tham khảo
[1] Calculating Damage with Miner’s Rule — Siemens Digital Industries Software. Truy cập ngày 09/11/2026.
[2] OrcaFlex Help — Fatigue analysis: How damage is calculated — Orcina. Truy cập ngày 09/11/2026.
[3] Fatigue Module 6.4 — Cumulative Damage Model — COMSOL. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.


