Ứng suất dùng để đọc đường cong phải được định nghĩa
Một mô hình phần tử hữu hạn có thể trả về ứng suất danh nghĩa, ứng suất kết cấu gần mối hàn và đỉnh ứng suất tại một góc hình học. Các con số không tự thuộc cùng một đường cong S-N. Khi thay cách trích ứng suất mà giữ nguyên đường cong, kết quả tuổi thọ có thể thay đổi chỉ do phương pháp xử lý dữ liệu.
Nghiên cứu của Zhang và Maddox phân biệt cách đánh giá theo ứng suất danh nghĩa, hot-spot và phương pháp cục bộ. Đường cong theo danh nghĩa hoặc hot-spot từ mẫu hàn đã chứa ảnh hưởng nhất định của hình học khía tại nơi khởi đầu nứt; một phương pháp ứng suất cục bộ phải đi cùng cơ sở đường cong phù hợp.[1] Không coi ứng suất lớn nhất mà phần mềm hiển thị là đầu vào mặc định cho mọi phương pháp.
Bài sử dụng trường ứng suất được dựng và các phép ngoại suy để kiểm tra độ nhạy. Khoảng cách, hệ số và đường cong đều là giả định phục vụ phân tích, không phải hướng dẫn chọn điểm theo một tiêu chuẩn cụ thể.
Ngoại suy từ hai điểm: giữ rõ tọa độ và dấu
Xét một đường trên bề mặt tấm, vuông góc với chân hàn. Đặt x = 0 tại vị trí cần ngoại suy; x dương đi ra xa chân. Với tấm dày t = 20 mm, chọn hai điểm minh họa x₁ = 8 mm và x₂ = 20 mm. Ở một trạng thái tải đàn hồi, thành phần ứng suất cùng phương tại hai điểm là s₁ = 112 MPa và s₂ = 94 MPa.
Giả định ứng suất kết cấu trên đoạn được biểu diễn bằng đường thẳng s(x) = a + bx. Giải hai phương trình cho giá trị tại x = 0:
s_h = (x₂s₁ - x₁s₂)/(x₂ - x₁)
Kết quả s_h = 124 MPa, lớn hơn cả hai giá trị đo. Điều này phù hợp với gradient giảm khi đi ra xa. Phép tính không phải lấy trung bình 103 MPa. Hệ số của hai số đo lần lượt là 5/3 và -2/3; hệ số âm là đặc trưng của ngoại suy ngoài khoảng quan sát.
Maddox mô tả cách xác định hot-spot từ phân bố ứng suất và lưu ý loại vị trí hàn, vị trí chân trên mặt tấm hay mép tấm có thể cần cách xử lý khác nhau.[2] Vì vậy hai điểm trong ví dụ chỉ chứng minh số học. Chưa thể dùng chúng cho một liên kết khác nếu chưa chọn phương pháp đánh giá thích hợp.
Ngoại suy khuếch đại thành phần bất định nào?
Giả sử mỗi s_i có độ không đảm bảo chuẩn u = 2 MPa, hai sai số độc lập, tọa độ được xem là chính xác. Do biểu thức tuyến tính, độ không đảm bảo của kết quả là:
u_h² = (5u₁/3)² + (2u₂/3)²
Ta được u_h = 3,590110 MPa. Tăng số chữ số báo cáo không làm giảm giá trị này. Nếu nhân hệ số phủ 2 chỉ để khảo sát, khoảng mở rộng tương ứng là ±7,180220 MPa; chưa gán mức bao phủ xác suất nếu không có giả thiết phân phối phù hợp.
Nếu hai giá trị cùng chịu một độ lệch hiệu chuẩn b giống nhau, sai số đầu ra bằng (5/3 - 2/3)b = b. Thành phần chung không bị khuếch đại như hai thành phần độc lập. Ngược lại, hai sai số bằng nhau về độ lớn nhưng trái dấu sẽ cộng bất lợi. Vì thế phải giữ cấu trúc tương quan, không gộp mọi nguồn sai số thành hai cột độc lập.
Sai số vị trí đo là vấn đề riêng: đổi x làm đổi hệ số ngoại suy và ứng suất lấy tại vị trí đó. Không thể xử lý đầy đủ chỉ bằng việc cộng thêm sai số áp suất hay lực vào s_i.
Nếu trường thực có độ cong
Đặt z = x/t không có đơn vị. Dựng một trường kết cấu giả định:
s(z) = 120 - 30z + 10z²
Các hệ số có đơn vị MPa. Giá trị tại z = 0 đúng bằng 120 MPa. Tại z = 0,4 và 1, giá trị lần lượt là 109,6 và 100 MPa. Ngoại suy tuyến tính từ hai điểm này lại cho 116 MPa, thấp hơn đích 4 MPa.
Có thể suy tổng quát: với s = a + bz + cz², giao điểm của đường qua z₁ và z₂ tại gốc bằng a - cz₁z₂. Sai lệch do độ cong là -cz₁z₂. Nó không biến mất khi máy đo hoàn toàn chính xác. Đây là sai khác mô hình, cần phân biệt với độ không đảm bảo đo.
Thêm điểm z = 0,7 cho s = 103,9 MPa và khớp đúng đa thức bậc hai sẽ khôi phục 120 MPa trong bài dựng. Điều này không có nghĩa luôn nên dùng bậc cao nhất; với số đo nhiễu hoặc vùng chịu ảnh hưởng khía rất gần chân, mô hình phù hợp cần được kiểm tra độc lập.
Hệ quả trên một đường cong minh họa
Giả sử riêng trường hợp tải tỷ lệ đã xác định được khoảng hot-spot S_h = 124 MPa và khoảng danh nghĩa tương ứng 80 MPa. Hệ số chuyển cho trường hợp này là K = 1,55. Đặt đường cong giả định phù hợp với định nghĩa hot-spot:
N = 2.000.000(100/S_h)³
Kết quả N = 1.048.974,523 chu kỳ. Nếu xét trường độ cong ở phần trước và nhập 116 thay vì 120 MPa vào cùng dạng m = 3, tuổi thọ được tính cao hơn theo hệ số (120/116)³ = 1,107056, tức khoảng 10,71%. Sai khác 3,33% trong ứng suất đã thành sai khác lớn hơn trong số chu kỳ.
Không nhân K thêm lần nữa nếu S_h đã là giá trị sau chuyển. Ngược lại, cũng không bỏ ảnh hưởng kết cấu cần thiết khi đầu vào còn là danh nghĩa. Vấn đề nằm ở định nghĩa của từng cột, không nằm ở việc luôn chọn con số lớn hơn.
Không có một hệ số duy nhất cho mọi tổ hợp tải
SACS định nghĩa SCF khi chuyển từ ứng suất danh nghĩa sang hot-spot và đánh giá các trạng thái tải trước khi hình thành khoảng ứng suất.[3] Một hệ số của tải kéo không nhất thiết dùng được cho tải uốn tại cùng vị trí.
Dựng hai thành phần cùng phương ứng suất, cùng tần số. Thành phần danh nghĩa thứ nhất có biên độ 40 MPa và hệ số 1,8; thành phần thứ hai có biên độ 25 MPa và hệ số 2,6. Tại điểm xét, biên độ riêng sau chuyển là 72 và 65 MPa. Nếu lệch pha 90 độ:
s_h(t) = 72 sin(ωt) + 65 cos(ωt)
Biên độ tổng bằng √(72² + 65²) = 97 MPa; khoảng bằng 194 MPa. Nếu cùng pha, khoảng là 274 MPa; nếu ngược pha, còn 14 MPa. Ba trường hợp có cùng biên độ riêng và SCF nhưng lịch sử tổng khác nhau.
Cộng hai khoảng riêng luôn cho 274 MPa, tương đương giả định cùng pha trong bài này. So với trường hợp lệch 90 độ, hư hỏng trên nhánh m = 3 bị nhân khoảng 2,817387. Đây là phép so cho hai thành phần cùng phương; tải nhiều trục cần tiêu chí phù hợp hơn việc cộng tùy ý các đại lượng vô hướng.
Hồ sơ cần lưu để kết quả có thể kiểm lại
Nên lưu đường trích ứng suất, hệ trục, tọa độ các điểm, loại phần tử, trạng thái tải và công thức chuyển. Khi tinh lưới, kiểm sự ổn định của đại lượng đã định nghĩa, đồng thời giữ phương pháp trích nhất quán. Thay đổi đỉnh ứng suất ngay tại góc sắc không tự chứng minh hot-spot đã thay đổi cùng tỷ lệ.
Một kết luận mỏi cần đồng thời xác định vị trí có khả năng nứt, đại lượng ứng suất, đường cong tương ứng và lịch tải. Ngoại suy và SCF là các bước nối giữa mô hình kết cấu với dữ liệu mỏi; chúng chỉ có ý nghĩa khi đường nối ấy được mô tả đầy đủ.
Nguồn tham khảo
[1] Estimating Long-Endurance Fatigue Strength of Girth-Welded Pipes using Local Stress Approach, OMAE 2014 — Yanhui Zhang và Steve J. Maddox, TWI. Truy cập ngày 09/11/2026.
[2] Recommended Hot-Spot Stress Design S-N Curves for Fatigue Assessment of FPSOs, 2001 — S. J. Maddox, TWI. Truy cập ngày 09/11/2026.
[3] Fatigue — Definitions and Overview — Bentley SACS 25.0.2. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Một mô hình phần tử hữu hạn có thể trả về ứng suất danh nghĩa, ứng suất kết cấu gần mối hàn và đỉnh ứng suất tại một góc hình học. Các con số không tự thuộc cùng một đường cong S-N. Khi thay cách trích ứng suất mà giữ nguyên đường cong, kết quả tuổi thọ có thể thay đổi chỉ do phương pháp xử lý dữ liệu.
Nghiên cứu của Zhang và Maddox phân biệt cách đánh giá theo ứng suất danh nghĩa, hot-spot và phương pháp cục bộ. Đường cong theo danh nghĩa hoặc hot-spot từ mẫu hàn đã chứa ảnh hưởng nhất định của hình học khía tại nơi khởi đầu nứt; một phương pháp ứng suất cục bộ phải đi cùng cơ sở đường cong phù hợp.[1] Không coi ứng suất lớn nhất mà phần mềm hiển thị là đầu vào mặc định cho mọi phương pháp.
Bài sử dụng trường ứng suất được dựng và các phép ngoại suy để kiểm tra độ nhạy. Khoảng cách, hệ số và đường cong đều là giả định phục vụ phân tích, không phải hướng dẫn chọn điểm theo một tiêu chuẩn cụ thể.
Ngoại suy từ hai điểm: giữ rõ tọa độ và dấu
Xét một đường trên bề mặt tấm, vuông góc với chân hàn. Đặt x = 0 tại vị trí cần ngoại suy; x dương đi ra xa chân. Với tấm dày t = 20 mm, chọn hai điểm minh họa x₁ = 8 mm và x₂ = 20 mm. Ở một trạng thái tải đàn hồi, thành phần ứng suất cùng phương tại hai điểm là s₁ = 112 MPa và s₂ = 94 MPa.
Giả định ứng suất kết cấu trên đoạn được biểu diễn bằng đường thẳng s(x) = a + bx. Giải hai phương trình cho giá trị tại x = 0:
s_h = (x₂s₁ - x₁s₂)/(x₂ - x₁)
Kết quả s_h = 124 MPa, lớn hơn cả hai giá trị đo. Điều này phù hợp với gradient giảm khi đi ra xa. Phép tính không phải lấy trung bình 103 MPa. Hệ số của hai số đo lần lượt là 5/3 và -2/3; hệ số âm là đặc trưng của ngoại suy ngoài khoảng quan sát.
Maddox mô tả cách xác định hot-spot từ phân bố ứng suất và lưu ý loại vị trí hàn, vị trí chân trên mặt tấm hay mép tấm có thể cần cách xử lý khác nhau.[2] Vì vậy hai điểm trong ví dụ chỉ chứng minh số học. Chưa thể dùng chúng cho một liên kết khác nếu chưa chọn phương pháp đánh giá thích hợp.
Ngoại suy khuếch đại thành phần bất định nào?
Giả sử mỗi s_i có độ không đảm bảo chuẩn u = 2 MPa, hai sai số độc lập, tọa độ được xem là chính xác. Do biểu thức tuyến tính, độ không đảm bảo của kết quả là:
u_h² = (5u₁/3)² + (2u₂/3)²
Ta được u_h = 3,590110 MPa. Tăng số chữ số báo cáo không làm giảm giá trị này. Nếu nhân hệ số phủ 2 chỉ để khảo sát, khoảng mở rộng tương ứng là ±7,180220 MPa; chưa gán mức bao phủ xác suất nếu không có giả thiết phân phối phù hợp.
Nếu hai giá trị cùng chịu một độ lệch hiệu chuẩn b giống nhau, sai số đầu ra bằng (5/3 - 2/3)b = b. Thành phần chung không bị khuếch đại như hai thành phần độc lập. Ngược lại, hai sai số bằng nhau về độ lớn nhưng trái dấu sẽ cộng bất lợi. Vì thế phải giữ cấu trúc tương quan, không gộp mọi nguồn sai số thành hai cột độc lập.
Sai số vị trí đo là vấn đề riêng: đổi x làm đổi hệ số ngoại suy và ứng suất lấy tại vị trí đó. Không thể xử lý đầy đủ chỉ bằng việc cộng thêm sai số áp suất hay lực vào s_i.
Nếu trường thực có độ cong
Đặt z = x/t không có đơn vị. Dựng một trường kết cấu giả định:
s(z) = 120 - 30z + 10z²
Các hệ số có đơn vị MPa. Giá trị tại z = 0 đúng bằng 120 MPa. Tại z = 0,4 và 1, giá trị lần lượt là 109,6 và 100 MPa. Ngoại suy tuyến tính từ hai điểm này lại cho 116 MPa, thấp hơn đích 4 MPa.
Có thể suy tổng quát: với s = a + bz + cz², giao điểm của đường qua z₁ và z₂ tại gốc bằng a - cz₁z₂. Sai lệch do độ cong là -cz₁z₂. Nó không biến mất khi máy đo hoàn toàn chính xác. Đây là sai khác mô hình, cần phân biệt với độ không đảm bảo đo.
Thêm điểm z = 0,7 cho s = 103,9 MPa và khớp đúng đa thức bậc hai sẽ khôi phục 120 MPa trong bài dựng. Điều này không có nghĩa luôn nên dùng bậc cao nhất; với số đo nhiễu hoặc vùng chịu ảnh hưởng khía rất gần chân, mô hình phù hợp cần được kiểm tra độc lập.
Hệ quả trên một đường cong minh họa
Giả sử riêng trường hợp tải tỷ lệ đã xác định được khoảng hot-spot S_h = 124 MPa và khoảng danh nghĩa tương ứng 80 MPa. Hệ số chuyển cho trường hợp này là K = 1,55. Đặt đường cong giả định phù hợp với định nghĩa hot-spot:
N = 2.000.000(100/S_h)³
Kết quả N = 1.048.974,523 chu kỳ. Nếu xét trường độ cong ở phần trước và nhập 116 thay vì 120 MPa vào cùng dạng m = 3, tuổi thọ được tính cao hơn theo hệ số (120/116)³ = 1,107056, tức khoảng 10,71%. Sai khác 3,33% trong ứng suất đã thành sai khác lớn hơn trong số chu kỳ.
Không nhân K thêm lần nữa nếu S_h đã là giá trị sau chuyển. Ngược lại, cũng không bỏ ảnh hưởng kết cấu cần thiết khi đầu vào còn là danh nghĩa. Vấn đề nằm ở định nghĩa của từng cột, không nằm ở việc luôn chọn con số lớn hơn.
Không có một hệ số duy nhất cho mọi tổ hợp tải
SACS định nghĩa SCF khi chuyển từ ứng suất danh nghĩa sang hot-spot và đánh giá các trạng thái tải trước khi hình thành khoảng ứng suất.[3] Một hệ số của tải kéo không nhất thiết dùng được cho tải uốn tại cùng vị trí.
Dựng hai thành phần cùng phương ứng suất, cùng tần số. Thành phần danh nghĩa thứ nhất có biên độ 40 MPa và hệ số 1,8; thành phần thứ hai có biên độ 25 MPa và hệ số 2,6. Tại điểm xét, biên độ riêng sau chuyển là 72 và 65 MPa. Nếu lệch pha 90 độ:
s_h(t) = 72 sin(ωt) + 65 cos(ωt)
Biên độ tổng bằng √(72² + 65²) = 97 MPa; khoảng bằng 194 MPa. Nếu cùng pha, khoảng là 274 MPa; nếu ngược pha, còn 14 MPa. Ba trường hợp có cùng biên độ riêng và SCF nhưng lịch sử tổng khác nhau.
Cộng hai khoảng riêng luôn cho 274 MPa, tương đương giả định cùng pha trong bài này. So với trường hợp lệch 90 độ, hư hỏng trên nhánh m = 3 bị nhân khoảng 2,817387. Đây là phép so cho hai thành phần cùng phương; tải nhiều trục cần tiêu chí phù hợp hơn việc cộng tùy ý các đại lượng vô hướng.
Hồ sơ cần lưu để kết quả có thể kiểm lại
Nên lưu đường trích ứng suất, hệ trục, tọa độ các điểm, loại phần tử, trạng thái tải và công thức chuyển. Khi tinh lưới, kiểm sự ổn định của đại lượng đã định nghĩa, đồng thời giữ phương pháp trích nhất quán. Thay đổi đỉnh ứng suất ngay tại góc sắc không tự chứng minh hot-spot đã thay đổi cùng tỷ lệ.
Một kết luận mỏi cần đồng thời xác định vị trí có khả năng nứt, đại lượng ứng suất, đường cong tương ứng và lịch tải. Ngoại suy và SCF là các bước nối giữa mô hình kết cấu với dữ liệu mỏi; chúng chỉ có ý nghĩa khi đường nối ấy được mô tả đầy đủ.
Nguồn tham khảo
[1] Estimating Long-Endurance Fatigue Strength of Girth-Welded Pipes using Local Stress Approach, OMAE 2014 — Yanhui Zhang và Steve J. Maddox, TWI. Truy cập ngày 09/11/2026.
[2] Recommended Hot-Spot Stress Design S-N Curves for Fatigue Assessment of FPSOs, 2001 — S. J. Maddox, TWI. Truy cập ngày 09/11/2026.
[3] Fatigue — Definitions and Overview — Bentley SACS 25.0.2. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
