Từ điển học thuật Kỹ thuật và công nghệ

Yếu tố cường độ ứng suất là gì? Bản chất và tính toán

Tiếng Anhstress intensity factor

Tên gọi kháchệ số cường độ ứng suấtstress intensity factorSIF

Yếu tố cường độ ứng suất (stress intensity factor, K) là đại lượng trong cơ học đứt gãy đàn hồi tuyến tính đặc trưng cho độ lớn của trường ứng suất và biến dạng tiệm cận tại đỉnh vết nứt dưới tác dụng của ngoại lực và điều kiện biên hình học.

223 lượt xem Cập nhật 25/8/2026

Yếu tố cường độ ứng suất (stress intensity factor, ký hiệu là K) là đại lượng cơ học then chốt trong lý thuyết cơ học đứt gãy đàn hồi tuyến tính (linear elastic fracture mechanics, LEFM), dùng để định lượng và đặc trưng cho mức độ tập trung của trường ứng suất và trường biến dạng tại vùng lân cận đỉnh vết nứt khi vật thể chịu tác dụng của ngoại lực và các điều kiện biên hình học. Đại lượng này phản ánh trực tiếp cường độ của trường kỳ dị ứng suất tiệm cận, cho phép đánh giá trạng thái nguy hiểm của khuyết tật nứt và dự đoán điều kiện lan truyền không ổn định của vết nứt trong vật liệu rắn đàn hồi.

Bản chất vật lý và cơ sở lý thuyết

Trong vật liệu chứa vết nứt lý tưởng chịu tải trọng ngoài, các phương trình lý thuyết đàn hồi cổ điển cho thấy ứng suất tại mũi vết nứt tiến tới vô cùng theo quy luật kỳ dị. Để giải quyết bài toán này mà không cần xét chi tiết quá trình vi cấu trúc phức tạp tại vùng biến dạng vi mô, George R. Irwin vào năm 1957 đã đề xuất khái niệm hệ số cường độ ứng suất, qua đó cô lập được mức độ gay gắt của tải trọng cơ học và hình học mẫu vào một tham số vô hướng duy nhất là K.

Trong hệ tọa độ cực phẳng đặt tại đỉnh vết nứt, với bán kính r đo từ đỉnh nứt và góc phương vị theta đo từ mặt phẳng chứa vết nứt, phân bố của các thành phần trường ứng suất tiệm cận khi r tiến dần về không được biểu diễn tổng quát theo quan hệ giải tích của Irwin năm 1957:

σij(r,θ)=K2πrfij(θ)\sigma_{ij}(r,\theta) = \frac{K}{\sqrt{2\pi r}} f_{ij}(\theta)

Trong đó:

  • σij\sigma_{ij} là các thành phần của tensor ứng suất tại tọa độ cực tương ứng.
  • KK là hệ số cường độ ứng suất phụ thuộc vào ngoại lực tác dụng và kích thước hình học của mẫu nứt.
  • rr là khoảng cách bán kính từ đỉnh vết nứt tới điểm khảo sát.
  • θ\theta là góc hợp bởi phương khảo sát và trục đối xứng của mặt phẳng vết nứt.
  • fij(θ)f_{ij}(\theta) là các hàm góc không thứ nguyên, hoàn toàn xác định theo chế độ chất tải và góc phương vị theta.

Khi xét thành phần ứng suất kéo pháp tuyến vuông góc với mặt phẳng nứt, biểu thức ứng suất tiệm cận thu gọn về dạng:

σθθ(r,θ)≈K2πrf(θ)\sigma_{\theta\theta}(r,\theta) \approx \frac{K}{\sqrt{2\pi r}} f(\theta)

Đặc tính kỳ dị tỉ lệ nghịch với căn bậc hai của bán kính r cho thấy trường ứng suất tăng vọt khi tiếp cận đỉnh nứt. Tuy nhiên, toàn bộ biên độ của trường kỳ dị này được điều tiết và tỉ lệ thuận tuyến tính với giá trị K. Do đó, chỉ cần xác định được giá trị K cho một kết cấu cụ thể, toàn bộ trường ứng suất và biến dạng cục bộ tại đỉnh vết nứt hoàn toàn được làm rõ.

Ba chế độ tải trọng cơ bản trong cơ học đứt gãy

Theo phân loại chuẩn mực của cơ học đứt gãy được tổng hợp bởi Anderson năm 2017, biến dạng cục bộ tại đỉnh vết nứt dưới mọi điều kiện tải trọng phức tạp đều có thể phân tích thành tổ hợp tuyến tính của ba chế độ tải trọng cơ bản, ký hiệu tương ứng bằng các chỉ số La Mã I, II và III.

Chế độ I: Mở vết nứt

Chế độ I (Mode I, opening mode) xảy ra khi ngoại lực tác dụng vuông góc với mặt phẳng của vết nứt, làm cho hai bề mặt của bờ vết nứt tách xa nhau theo phương đối xứng. Đây là chế độ phổ biến nhất trong thực tế kỹ thuật và cũng là chế độ nguy hiểm nhất dẫn đến hiện tượng phá hủy giòn thảm khốc của các kết cấu chịu kéo hoặc uốn. Hệ số cường độ ứng suất tương ứng được ký hiệu là KIK_I.

Chế độ II: Trượt trong mặt phẳng

Chế độ II (Mode II, in-plane shear mode hay sliding mode) xuất hiện khi tải trọng tác dụng gây ra ứng suất cắt trong mặt phẳng nứt, làm hai bờ vết nứt trượt tương đối với nhau theo hướng vuông góc với tiền duyên vết nứt. Hệ số cường độ ứng suất tương ứng được ký hiệu là KIIK_{II}.

Chế độ III: Xé ngoài mặt phẳng

Chế độ III (Mode III, out-of-plane shear mode hay tearing mode) đặc trưng cho hiện tượng ứng suất cắt tác dụng làm hai bờ vết nứt trượt tương đối theo phương song song với tiền duyên vết nứt (ngoài mặt phẳng). Hiện tượng này thường gặp trong các chi tiết trục quay chịu mômen xoắn lớn hoặc trong các tấm chịu cắt ngang. Hệ số cường độ ứng suất tương ứng được ký hiệu là KIIIK_{III}.

Chế độ chất tải Ký hiệu cường độ ứng suất Đặc điểm chuyển vị tại bờ nứt Thành phần ứng suất chủ đạo
Chế độ I (Mode I - Mở) KIK_I Hai bờ nứt tách rời đối xứng theo phương vuông góc mặt nứt Ứng suất kéo pháp tuyến σyy\sigma_{yy}
Chế độ II (Mode II - Trượt) KIIK_{II} Hai bờ nứt trượt tương đối trong mặt phẳng nứt Ứng suất cắt trong mặt phẳng σxy\sigma_{xy}
Chế độ III (Mode III - Xé) KIIIK_{III} Hai bờ nứt trượt tương đối ngoài mặt phẳng song song tiền duyên Ứng suất cắt ngoài mặt phẳng σxz\sigma_{xz}

Công thức tính toán và hệ số hiệu chỉnh hình học

Đối với một vết nứt chiều dài xác định trong một tấm vô hạn chịu ứng suất kéo đều danh nghĩa ở vô cực, hệ số cường độ ứng suất Chế độ I có dạng giải tích kinh điển của Irwin:

KI=σπaK_I = \sigma \sqrt{\pi a}

Trong đó σ\sigma là ứng suất danh nghĩa tác dụng từ xa và aa là nửa chiều dài vết nứt xuyên tâm (hoặc toàn bộ chiều dài đối với vết nứt mép bề mặt).

Đối với các kết cấu có kích thước hữu hạn và hình học phức tạp, hiệu ứng tương tác của các bề mặt tự do và điều kiện biên làm sai lệch phân bố ứng suất so với trường hợp tấm vô hạn. Khi đó, một hệ số hiệu chỉnh hình học không thứ nguyên, thường ký hiệu là Y hoặc f, được đưa vào công thức:

KI=σπa YK_I = \sigma \sqrt{\pi a}\,Y

Theo sổ tay tra cứu của Tada, Paris và Irwin năm 2000, đối với mẫu tấm chịu kéo đơn trục có vết nứt mép với chiều dài nứt a và bề rộng tấm W, hàm hiệu chỉnh hình học Y được biểu diễn bằng đa thức phụ thuộc vào tỉ số a/W:

Y=1.12−0.23aW+10.55(aW)2−21.72(aW)3+30.39(aW)4Y = 1.12 - 0.23\frac{a}{W} + 10.55\left(\frac{a}{W}\right)^2 - 21.72\left(\frac{a}{W}\right)^3 + 30.39\left(\frac{a}{W}\right)^4

Trong các thí nghiệm đo độ dai phá hủy chuẩn sử dụng mẫu kéo nhỏ gọn (compact tension, CT), công thức chuẩn xác định hệ số cường độ ứng suất có dạng:

KI=FBW f ⁣(aW)K_I = \frac{F}{B \sqrt{W}}\,f\!\bigl(\tfrac{a}{W}\bigr)

Trong đó:

  • FF là lực kéo tác dụng lên mẫu thử.
  • BB là chiều dày của mẫu thử.
  • WW là chiều rộng hiệu dụng của mẫu thử đo từ tâm lỗ chất tải đến mép sau.
  • f(a/W)f(a/W) là hàm hình học chuẩn hóa của mẫu kéo nhỏ gọn.

Mối quan hệ nhiệt động và tích phân J

Năm 1968, James R. Rice đã xây dựng nền tảng giải tích của tích phân không phụ thuộc đường bao J (J-integral), đại diện cho tốc độ giải phóng năng lượng biến dạng đàn hồi trên một đơn vị diện tích bề mặt vết nứt mới tạo thành. Trong phạm vi vật liệu đàn hồi tuyến tính, Rice năm 1968 đã chứng minh quan hệ tương đương trực tiếp giữa tích phân J và hệ số cường độ ứng suất:

J=KI2E′J = \frac{K_I^2}{E'}

Trong đó mô đun đàn hồi tương đương E′E' phụ thuộc chặt chẽ vào trạng thái ứng suất - biến dạng của vật thể:

  • Đối với trạng thái ứng suất phẳng (plane stress, thường gặp ở tấm rất mỏng): E′=EE' = E.
  • Đối với trạng thái biến dạng phẳng (plane strain, áp dụng cho cấu kiện dày): E′=E1−ν2E' = \frac{E}{1-\nu^2}, với EE là mô đun đàn hồi Young và ν\nu là hệ số Poisson của vật liệu.

Khi đó biểu thức hoàn chỉnh cho điều kiện biến dạng phẳng được biểu diễn là:

J=KI2E′(E′=E1−ν2)J = \frac{K_I^2}{E'}\quad \bigl(E' = \tfrac{E}{1-\nu^2}\bigr)

Mối liên hệ này nối liền cách tiếp cận năng lượng toàn cục của Griffith với cách tiếp cận trường ứng suất cục bộ của Irwin, tạo nên sự thống nhất toán học hoàn chỉnh cho cơ học đứt gãy.

Ứng dụng trong đánh giá độ bền mỏi và định luật Paris

Khi kết cấu chịu tải trọng chu kỳ thay đổi theo thời gian, vết nứt sẽ lan truyền dần dần qua từng chu kỳ chất tải dưới giá trị tới hạn tĩnh. Hiện tượng này được gọi là sự phát triển vết nứt mỏi (fatigue crack growth). Năm 1963, Paul Paris và Fazil Erdogan đã phát hiện ra rằng động học phát triển vết nứt mỏi chịu sự chi phối quyết định của biên độ biến thiên hệ số cường độ ứng suất ΔK=Kmax⁡−Kmin⁡\Delta K = K_{\max} - K_{\min}.

Mối quan hệ định lượng giữa tốc độ tăng trưởng chiều dài vết nứt trên mỗi chu kỳ tải da/dNda/dN và biên độ ΔK\Delta K được thể hiện qua phương trình định luật Paris năm 1963:

dadN=C (ΔK)m\frac{da}{dN} = C\,(\Delta K)^m

Trong đó:

  • aa là chiều dài vết nứt tại chu kỳ khảo sát.
  • NN là số chu kỳ chất tải lặp.
  • ΔK\Delta K là khoảng biến thiên của hệ số cường độ ứng suất trong một chu kỳ tải trọng.
  • CC và mm là các thông số thực nghiệm đặc trưng cho từng loại vật liệu và môi trường thử nghiệm.

Định luật Paris cho phép các kỹ sư tính toán tích phân để dự báo chính xác số chu kỳ làm việc an toàn còn lại của kết cấu trước khi vết nứt đạt đến kích thước tới hạn gây phá hủy đột ngột.

Phương pháp xác định hệ số cường độ ứng suất

Để đánh giá an toàn kết cấu trong thực tế, các phương pháp xác định hệ số cường độ ứng suất được chia thành phương pháp giải tích, phương pháp số và phương pháp thực nghiệm:

  • Phương pháp hàm ứng suất và biến đổi giải tích: Sử dụng hàm ứng suất phức Westergaard hoặc phương pháp ánh xạ bảo giác của Muskhelishvili để tìm nghiệm giải tích cho các hình học đối xứng chuẩn.
  • Phương pháp phần tử hữu hạn (FEA): Là công cụ số phổ biến nhất hiện nay trong thiết kế kỹ thuật, tính toán K thông qua việc nội suy trường chuyển vị cục bộ quanh đỉnh nứt, kỹ thuật đóng vết nứt ảo (VCCT), hoặc đánh giá tích phân J trên các miền bao quanh mũi nứt.
  • Phương pháp thực nghiệm quang đàn hồi (Photoelasticity): Sử dụng mẫu thử từ vật liệu quang học đẳng hướng trong suốt, quan sát các dải vân đẳng sắc dưới ánh sáng phân cực để suy ra trực tiếp gradien ứng suất và giá trị K tại mũi nứt.
  • Phương pháp đo biến dạng bằng điện trở (Strain Gauge): Dán các cảm biến đo biến dạng chuyên dụng tại các vị trí xác định gần đỉnh nứt để đo biến dạng đàn hồi và tính toán ngược lại hệ số K.

Ranh giới áp dụng và chuyển tiếp sang cơ học đứt gãy đàn dẻo

Lý thuyết hệ số cường độ ứng suất được xây dựng dựa trên giả thiết vật liệu đàn hồi tuyến tính hoàn hảo. Tuy nhiên, trong thực tế các kim loại và hợp kim kết cấu luôn xảy ra hiện tượng chảy dẻo cục bộ tại đỉnh vết nứt do sự tập trung ứng suất vượt quá giới hạn chảy của vật liệu, tạo nên một vùng biến dạng dẻo (plastic zone) bao quanh đỉnh nứt.

Theo Anderson năm 2017, điều kiện bắt buộc để yếu tố cường độ ứng suất K giữ được ý nghĩa vật lý chính xác là điều kiện chảy dẻo quy mô nhỏ (small-scale yielding, SSY). Điều này đòi hỏi kích thước bán kính của vùng biến dạng dẻo phải nhỏ hơn rất nhiều so với chiều dài vết nứt, chiều dày của mẫu và khoảng cách tới các biên tự do của kết cấu.

Khi vật liệu có độ dẻo cao hoặc kết cấu làm việc ở nhiệt độ cao khiến vùng dẻo phát triển lớn, trường ứng suất không còn tuân theo quy luật kỳ dị đàn hồi tuyến tính. Trong trường hợp này, việc sử dụng đơn thuần hệ số cường độ ứng suất K sẽ dẫn đến sai số nghiêm trọng và không còn bảo toàn tính chất đặc trưng. Kỹ sư phải chuyển sang áp dụng các tham số của cơ học đứt gãy đàn dẻo (elastic-plastic fracture mechanics, EPFM), tiêu biểu là tích phân J mở rộng hoặc độ mở đỉnh vết nứt (crack tip opening displacement, CTOD).

Ý nghĩa thực tiễn trong kỹ thuật và kiểm định an toàn

Trong thực tiễn kỹ thuật công nghiệp và xây dựng, yếu tố cường độ ứng suất là tham số cốt lõi để thực hiện triết lý thiết kế chấp nhận khuyết tật (damage tolerance design):

  • Đánh giá độ dai phá hủy vật liệu: Giá trị tới hạn của hệ số cường độ ứng suất trong điều kiện biến dạng phẳng khi bắt đầu xảy ra sự lan truyền vết nứt không ổn định được định nghĩa là độ dai phá hủy KIcK_{Ic}. Đây là một hằng số cơ tính cơ bản của vật liệu, được tiêu chuẩn hóa trong các phép thử cơ học.
  • Xác định kích thước vết nứt giới hạn: Cho phép thiết lập kích thước vết nứt lớn nhất cho phép tồn tại trong kết cấu mà không gây nguy cơ gãy vỡ tức thời dưới cấp tải trọng thiết kế lớn nhất.
  • Lập kế hoạch kiểm tra không phá hủy (NDT): Kết hợp định luật Paris với hệ số K để xác định chu kỳ kiểm tra định kỳ bằng siêu âm, chụp ảnh bức xạ hoặc dòng điện xoáy cho các công trình trọng điểm như cầu thép, bồn áp lực, giàn khoan dầu khí và đường ống dẫn khí áp suất cao tại Việt Nam.

Câu hỏi thường gặp

Yếu tố cường độ ứng suất khác gì so với hệ số tập trung ứng suất?

Hệ số tập trung ứng suất (Kt) áp dụng cho các góc lượn hoặc lỗ khoét có bán kính góc lượn hữu hạn trong lý thuyết đàn hồi cổ điển, biểu thị tỉ số ứng suất cực đại so với ứng suất danh nghĩa. Ngược lại, yếu tố cường độ ứng suất (K) áp dụng cho vết nứt sắc nhọn lý thuyết nơi ứng suất tiệm cận vô cùng, đặc trưng cho biên độ của trường kỳ dị ứng suất tại đỉnh nứt.

Khi nào thì vết nứt bắt đầu lan truyền không ổn định?

Vết nứt sẽ lan truyền không ổn định dẫn đến sự phá hủy giòn kết cấu khi hệ số cường độ ứng suất tác dụng đạt hoặc vượt quá độ dai phá hủy tới hạn của vật liệu (KI >= KIc) trong điều kiện biến dạng phẳng.

Tại sao yếu tố cường độ ứng suất Chế độ I (Mode I) được quan tâm nhiều nhất?

Chế độ I (mở vết nứt dưới tải trọng kéo pháp tuyến) là chế độ phổ biến nhất trong thực tế kỹ thuật và có xu hướng làm vết nứt lan truyền nhanh nhất với năng lượng giải phóng thấp nhất so với chế độ trượt và xé.

Tài liệu tham khảo

  1. Irwin, G. R. (1957). Analysis of stresses and strains near the end of a crack traversing a plate. Journal of Applied Mechanics, 24(3), 361-364. DOI: 10.1115/1.4011547
  2. Anderson, T. L. (2017). Fracture Mechanics: Fundamentals and Applications (4th ed.). CRC Press. DOI: 10.1201/9781315370293
  3. Paris, P., & Erdogan, F. (1963). A critical analysis of crack propagation laws. Journal of Basic Engineering, 85(4), 528-534. DOI: 10.1115/1.3656900
  4. Rice, J. R. (1968). A path independent integral and the approximate analysis of strain concentration by notches and cracks. Journal of Applied Mechanics, 35(2), 379-386. DOI: 10.1115/1.3601206
  5. Tada, H., Paris, P. C., & Irwin, G. R. (2000). The Stress Analysis of Cracks Handbook (3rd ed.). ASME Press. DOI: 10.1115/1.801535