查看“︁应力强度因子”︁的源代码

[[File:Cracktipcoords.svg|thumb|300x300px|裂纹尖端的极坐标]]
'''应力强度因子'''（Stress Intensity Factors），符号为''{{mvar|K}}''，是表征弹性材料的裂纹尖端应力、应变状态控制失稳扩展的参量，是[[断裂力学]]、破壞力学中极其重要的一个参量。在断裂力学中，应力强度因子用于预测由远程载荷或[[残余应力]]引起的裂纹或凹口尖端附近的应力强度，<ref name="ander">{{cite book|last=Anderson|first=T. L.|year=2005|title=Fracture mechanics: fundamentals and applications|publisher=CRC Press}}</ref>是一个应用于均匀、线弹性材料的理论，对脆性材料的破坏非常适用，也可应用于在裂纹尖端表现出小尺度屈服的材料。

''{{mvar|K}}''的大小取决于试样的几何形状、裂纹或凹口的大小和位置，以及材料上载荷的大小和分布。可以写成：<ref>{{Cite book|last=Soboyejo|first=W. O.|url=http://worldcat.org/oclc/300921090|title=Mechanical properties of engineered materials|date=2003|publisher=Marcel Dekker|isbn=0-8247-8900-8|chapter=11.6.2 Crack Driving Force and Concept of Similitude|oclc=300921090|access-date=2024-03-28|archive-date=2021-05-31|archive-url=https://web.archive.org/web/20210531151500/https://www.worldcat.org/title/mechanical-properties-of-engineered-materials/oclc/300921090|dead-url=no}}</ref><ref>{{Cite book|last=Janssen|first=M. (Michael)|url=https://www.worldcat.org/oclc/57491375|title=Fracture mechanics|date=2004|publisher=Spon Press|others=Zuidema, J. (Jan), Wanhill, R. J. H.|isbn=0-203-59686-2|edition=2nd|location=London|pages=41|oclc=57491375}}</ref>

<math>K = \sigma \sqrt{\pi a} \, f(a/W)</math>

其中<math>f(a/W)</math>是试样几何依赖函数，{{mvar|a}}为长度，{{mvar|W}}宽度，{{mvar|&sigma;}}是施加的压力。

[[弹性力学]]理论预测，在[[极坐标系|极坐标]]（<math>r,\theta</math>）中，裂纹尖端附近的应力分布（<math>\sigma_{ij}</math>）具有以下形式：<ref>{{cite book|title=The Stress Analysis of Cracks Handbook|author1=Hiroshi Tada|author-link=Hiroshi Tada (engineer)|author2=P. C. Paris|author2-link=Paul C. Paris|author3=George R. Irwin|author3-link=George Rankine Irwin|publisher=American Society of Mechanical Engineers|edition=3rd|date=February 2000}}</ref>

<math>
   \sigma_{ij}(r, \theta) = \frac {K} {\sqrt{2 \pi r}}\,f_{ij} ( \theta) + \,\,\rm{higher\, order\, terms} </math>

其中{{mvar|K}}是应力强度因子（单位为应力×长度<sup>1/2</sup>），<math>f_{ij}</math>是一个随载荷和几何形状变化的无量纲量。理论上，当{{mvar|r}}趋近于0时，应力<math>\sigma_{ij}</math>趋近于<math>\infty</math>，导致应力奇异性。<ref name="notch3" />然而，实际上，这种关系在靠近裂尖处（{{mvar|r}}很小时）会失效，因为塑性通常发生在超过材料[[屈服强度]]的应力下，此时的线弹性解则不再适用。但如果裂尖的塑性区域相对裂纹长度很小，那么裂尖附近的渐近应力分布则仍然适用。

== 各种模式的应力强度因子 ==
[[File:Fracture_modes_v2.svg|thumb|330x330px|三种模式的裂纹荷载]]
1957年，George Rankine Irwin发现，裂纹周围的应力可以用一个称为应力强度因子的比例因子来表达。他发现，裂纹在任意荷载下可以分解为三种线性独立的裂纹模式。<ref name="suresh04">{{cite book|last1=Suresh|first1=S.|year=2004|title=Fatigue of Materials|publisher=Cambridge University Press|isbn=978-0-521-57046-6}}</ref>这些荷载类型被归类为I型、II型或III型，如图所示。I型是一种拉伸模式，裂纹表面直接分开。II型是一种滑动（平面剪切）模式，裂纹表面相互滑动，方向与裂纹前沿垂直。III型是一种撕裂（反平面剪切）模式，裂纹表面相对移动，并且与裂纹前沿平行。I型是工程设计中最常见的加载类型。

对应于三种不同的[[断裂模式]]，有三个不同的物理量：<math>K_I</math>、<math>K_{II}</math>、<math>K_{III}</math>，不同的下标用于指定三种不同模式的应力强度因子，这些因子的形式定义如下：<ref name="rooke7" />

<math>
\begin{align}
   K_{\rm I} & = \lim_{r\rightarrow 0} \sqrt{2\pi r}\,\sigma_{yy}(r,0) \\
   K_{\rm II} & = \lim_{r\rightarrow 0} \sqrt{2\pi r}\,\sigma_{yx}(r,0) \\
   K_{\rm III} & = \lim_{r\rightarrow 0} \sqrt{2\pi r}\,\sigma_{yz}(r,0) \,.
\end{align}
</math>
{| class="toccolours collapsible collapsed" style="text-align:left"
!Equations for stress and displacement fields
|}

== 能量释放速率与J积分 ==
在[[平面應力|平面应力]]条件下，裂纹在纯I型或纯II型荷载下的应变能释放速率（<math>G</math>）与应力强度因子相关：

<math>
    G_{\rm I} = K_{\rm I}^2\left(\frac{1}{E}\right)
 </math>

<math>
    G_{\rm II} = K_{\rm II}^2\left(\frac{1}{E}\right)
 </math>

其中<math>E</math>是[[杨氏模量]]，<math>\nu</math>是材料[[泊松比]]。假设材料各向同性、分布均匀且为线弹性，则认为裂纹会沿着初始裂纹的方向延伸。

对于[[无穷小应变理论|平面应变]]条件，等效关系为：

<math>
    G_{\rm I} = K_{\rm I}^2\left(\frac{1-\nu^2}{E}\right)\,
 </math>

<math>
    G_{\rm II} = K_{\rm II}^2\left(\frac{1-\nu^2}{E}\right)\,.
 </math>

对于纯III型荷载为：

<math>
    G_{\rm III} = K_{\rm III}^2\left(\frac{1}{2\mu}\right) = K_{\rm III}^2\left(\frac{1+\nu}{E}\right)
 </math>

其中<math>\mu</math>是[[剪切模量]]。对于平面应变中的一般荷载，则成立以下线性组合：

<math>
   G = J = \int_\Gamma \left(W~dx_2 - \mathbf{t}\cdot\cfrac{\partial\mathbf{u}}{\partial x_1}~ds\right) \,.
 </math>

与应力强度因子相对应的表征弹塑性材料的裂纹尖端应力应变场的参量是[[断裂韧性]]J积分，J积分也按照断裂模式分为<math>J_I</math>、<math>J_{II}</math>、<math>J_{III}</math>。

== 临界应力强度因子 ==
{{main|断裂韧性}}应力强度因子<math>K</math>是一个放大应力幅值的参数，其中包括几何参数<math>Y</math>（荷载类型）。任何模式情况下的应力强度与材料上的加载成正比。如果在材料中制作一个非常尖锐的裂纹或V型缺口，可以经验性地确定<math>K_\mathrm{I}</math>的最小值，这是裂纹传播所需的临界应力强度值。在平面应变中确定的I型加载的临界断裂韧性（<math>K_\mathrm{Ic}</math>），称为材料的临界断裂韧性。<math>K_\mathrm{Ic}</math>的单位是应力乘以距离的平方根（例如MN/m<sup>3/2</sup>）。<math>K_\mathrm{Ic}</math>的单位表明，必须在一定的临界距离内达到材料的断裂应力，才能达到<math>K_\mathrm{Ic}</math>，使裂纹扩展。I型临界应力强度因子<math>K_\mathrm{Ic}</math>是断裂力学中最常用的工程设计参数，对桥梁、建筑、飞机甚至钟铃的耐断材料设计十分重要。

=== G-准则 ===
G准则是一种断裂准则，它将临界应力强度因子（或断裂韧性）与三种模式的应力强度因子相关联。这个失效准则可以表述为：<ref name="sih">{{citation|title=Fracture mechanics applied to engineering problems-strain energy density fracture criterion|last1=Sih|first1=G. C.|last2=Macdonald|first2=B.|journal=Engineering Fracture Mechanics|volume=6|number=2|pages=361–386|year=1974|doi=10.1016/0013-7944(74)90033-2}}</ref>

<math>
   K_{\rm c}^2 = K_{\rm I}^2 + K_{\rm II}^2 + \frac{E'}{2\mu}\,K_{\rm III}^2
 </math>

其中，是断裂韧性，<math>E' = E/(1-\nu^2)</math>用于平面应变，<math>E' = E</math>用于平面应力。平面应力的临界应力强度因子通常被写成<math>K_{\rm c}</math>。

== 例子 ==
[[File:CrackInfinitePlate.svg|thumb|294x294px|在I型加载下的无限大平板中的裂纹]]

=== 无限平面：均匀单轴应力 ===
假设长度为<math>2a</math>的直裂纹，在一个具有均匀应力场<math>\sigma</math>的无限平面中，其应力强度因子：<ref name="notch3" /><ref name="rooke7" />

<math>
   K_\mathrm{I}=\sigma \sqrt{\pi a}
</math>
[[File:PennyShapedCrack.svg|thumb|161x161像素|无限域内受单轴张力作用的硬币形裂纹]]

=== 无限域：硬币形裂纹 ===
在无限域内一个半径为<math>a</math>的小裂缝尖端在单张应力σ下的应力强度因子：<ref name="ander2">{{cite book|last=Anderson|first=T. L.|year=2005|title=Fracture mechanics: fundamentals and applications|publisher=CRC Press}}</ref>

<math>
   K_{\rm I} = \frac{2}{\pi}\sigma\sqrt{\pi a} \,.
 </math>
[[File:CrackFinitePlate.svg|thumb|242x242像素|有限平面在I型加载下的裂缝]]

=== 有限平面：均匀单轴应力 ===
如果裂缝位于宽度为<math>2b</math>、高度为<math>2h</math>的有限平面的中心位置，则应力强度因子的近似关系为：<ref name="rooke7" />

<math>
   K_{\rm I} = \sigma \sqrt{\pi a}\left[\cfrac{1 - \frac{a}{2b} + 0.326\left(\frac{a}{b}\right)^2}{\sqrt{1 - \frac{a}{b}}}\right] \,.
 </math>

如果裂纹不位于宽度的中心位置，即<math>d \ne b</math>，则位置'''A'''处的应力强度因子可以通过级数展开来近似得到：<ref name="rooke7" /><ref name="isida66">Isida, M., 1966, ''Stress intensity factors for the tension of an eccentrically cracked strip'', Transactions of the ASME Applied Mechanics Section, v. 88, p.94.</ref>

<math>
   K_{\rm IA} =  \sigma \sqrt{\pi a}\left[1 + \sum_{n=2}^{M} C_n\left(\frac{a}{b}\right)^n\right]
 </math>

从应力强度曲线的拟合中可以基于不同<math>d</math>值，得到因子<math>C_n</math>。<ref name="rooke7" />裂纹'''B'''点可以得到一个类似（但不完全相同）的表达式。在A点和B点的应力强度因子的替代表达式是：<ref name="kathir">{{cite book|last1=Kathiresan|first1=K.|last2=Brussat|first2=T. R.|last3=Hsu|first3=T. M.|year=1984|title=Advanced life analysis methods. Crack Growth Analysis Methods for Attachment Lugs|publisher=Flight Dynamics Laboratory, Air Force Wright Aeronautical Laboratories, AFSC W-P Air Force Base, Ohio|number=AFWAL-TR-84-3080}}</ref>

<math>
   K_{\rm IA} = \sigma\sqrt{\pi a}\,\Phi_A \,\, , K_{\rm IB} = \sigma\sqrt{\pi a}\,\Phi_B
 </math>

其中

<math>
  \begin{align}
   \Phi_A &:= \left[\beta + \left(\frac{1-\beta}{4}\right)\left(1 + \frac{1}{4\sqrt{\sec\alpha_A}}\right)^2\right]\sqrt{\sec\alpha_A} \\
   \Phi_B &:= 1 + \left[\frac{\sqrt{\sec\alpha_{AB}} - 1}{1 + 0.21\sin\left\{8\,\tan^{-1}\left[\left(\frac{\alpha_A - \alpha_B}{\alpha_A + \alpha_B}\right)^{0.9}\right]\right\}}\right]
   \end{align}
 </math>

其中

<math>
  \beta := \sin\left(\frac{\pi\alpha_B}{\alpha_A+\alpha_B}\right) ~,~~ \alpha_A := \frac{\pi a}{2 d} 
   ~,~~ \alpha_B := \frac{\pi a}{4b - 2d} ~;~~ \alpha_{AB} := \frac{4}{7}\,\alpha_A + \frac{3}{7}\,\alpha_B \,.
 </math>

在上述表达试中，<math>d</math>是从裂纹中心到距离点A最近边界的距离。但是当<math>d=b</math>时，上述表达式不会简化为关于中心裂纹的近似表达式。
[[File:CrackEdgeFinitePlate.svg|thumb|261x261像素|单轴应力下有限平面中的边缘裂纹]]

=== 有限平面：边缘裂纹 ===
对于一个尺寸为<math>2h \times b</math>2的平面，包含长度<math>a</math>的无约束边缘裂纹，如果板的尺寸满足<math>h/b \ge 0.5</math>和<math>a/b \le 0.6</math>，则在单轴应力<math>\sigma</math>下，裂纹尖端的应力强度因子为：<ref name="notch3">{{cite journal |last1=Liu |first1=M. |display-authors=etal |year=2015 |title=An improved semi-analytical solution for stress at round-tip notches |url=http://drgan.org/wp-content/uploads/2014/07/032_EFM_2015.pdf |journal=Engineering Fracture Mechanics |volume=149 |pages=134–143 |doi=10.1016/j.engfracmech.2015.10.004 |s2cid=51902898 |access-date=2024-03-28 |archive-date=2018-07-13 |archive-url=https://web.archive.org/web/20180713213957/http://drgan.org/wp-content/uploads/2014/07/032_EFM_2015.pdf |dead-url=no }}</ref>

<math>
   K_{\rm I} = \sigma\sqrt{\pi a}\left[1.122 - 0.231\left(\frac{a}{b}\right) + 10.55\left(\frac{a}{b}\right)^2
    - 21.71\left(\frac{a}{b}\right)^3 + 30.382\left(\frac{a}{b}\right)^4\right] \,.
</math>

对于<math>h/b \ge 1</math>和 <math>a/b \ge 0.3</math>的情况，应力强度因子可近似得：

<math>
  K_{\rm I} = \sigma\sqrt{\pi a}\left[\frac{1 + 3\frac{a}{b}}{2\sqrt{\pi\frac{a}{b}}\left(1-\frac{a}{b}\right)^{3/2}}\right] \,.
</math>
[[File:CrackSlantedPlateBiaxialLoad.svg|thumb|177x177px|双轴加载平面中的倾斜裂缝。]]

=== 无限平面：双轴应力场中的倾斜裂纹 ===
对于一个长度为<math>2a</math>的倾斜裂纹，在双轴向应力场中，应力在<math>y</math>方向为<math>\sigma</math>，在<math>x</math>方向为<math>\alpha\sigma</math>，应力强度因子为：<ref name="rooke7" /><ref name="sih" />

<math>
  \begin{align}
  K_{\rm I} & = \sigma\sqrt{\pi a}\left(\cos^2\beta + \alpha \sin^2\beta\right) \\
  K_{\rm II} & = \sigma\sqrt{\pi a}\left(1- \alpha\right)\sin\beta\cos\beta 
  \end{align}
</math>
其中<math>\beta</math>是裂缝与<math>x</math>轴所成的角度
[[File:CrackFinitePlatePointForce.svg|thumb|203x203px|在局部作用力<math>F_x</math>和<math>F_y</math>下的平面裂纹]]

=== 平面内定点应力的裂缝 ===
一个尺寸为<math>2h \times 2b</math>的平面，含有长度为<math>2a</math>的裂缝。在板的点（<math>x,y</math>）处施加具有分量<math>F_x</math>和<math>F_y</math>的力。

在平面比裂纹尺寸大，且力的位置相对靠近裂纹的情况下，即<math>h \gg a</math>，<math>b \gg a</math>，<math>x \ll b</math>，<math>y \ll h</math>，可将平面视为无限大，这种情况下，裂纹尖端B（<math>x = a</math>）处的<math>F_x</math>的应力强度因子为：<ref name="sih" /><ref name="erdo62">{{citation|title=On the stress distribution in plates with collinear cuts under arbitrary loads|last=Erdogan|first=F.|journal=Proceedings of the Fourth US National Congress of Applied Mechanics|volume=1|pages=547–574|year=1962}}</ref>

<math>
  \begin{align}
  K_{\rm I} & = \frac{F_x}{2\sqrt{\pi a}}\left(\frac{\kappa -1}{\kappa+1}\right)
            \left[G_1 + \frac{1}{\kappa-1} H_1\right] \\
  K_{\rm II} & = \frac{F_x}{2\sqrt{\pi a}}
            \left[G_2 + \frac{1}{\kappa+1} H_2\right] 
  \end{align}
</math>

其中

<math>
   \begin{align}
   G_1 & = 1 - \text{Re}\left[\frac{a+z}{\sqrt{z^2-a^2}}\right] \,,\,\,
   G_2 = - \text{Im}\left[\frac{a+z}{\sqrt{z^2-a^2}}\right] \\
   H_1 & = \text{Re}\left[\frac{a(\bar{z}-z)}{(\bar{z}-a)\sqrt{{\bar{z}}^2-a^2}}\right] \,,\,\,
   H_2 = -\text{Im}\left[\frac{a(\bar{z}-z)}{(\bar{z}-a)\sqrt{{\bar{z}}^2-a^2}}\right] 
   \end{align}
</math>

<math>z = x + iy</math>，<math>\bar{z} = x - iy</math>，<math> \kappa = 3-4\nu</math>用于平面应变，<math>\kappa=    (3-\nu)/(1+\nu)</math>用于平面应力，其中<math>\nu</math>是[[泊松比]]。

在B端点<math>F_y</math>的应力强度因子为：

<math>
  \begin{align}
  K_{\rm I} & = \frac{F_y}{2\sqrt{\pi a}}
            \left[G_2 - \frac{1}{\kappa+1} H_2\right] \\
  K_{\rm II} & = -\frac{F_y}{2\sqrt{\pi a}}\left(\frac{\kappa -1}{\kappa+1}\right)
            \left[G_1 - \frac{1}{\kappa-1} H_1\right] \,.
  \end{align}
</math>

点A（<math>x = -a</math>）处的应力强度因子可根据上述关系确定。对于在位置<math>(x,y)</math> 处施加<math>F_x</math>荷载的应力强度因子为：

<math>
   K_{\rm I}(-a; x,y) = -K_{\rm I}(a; -x,y) \,,\,\, 
   K_{\rm II}(-a; x,y) = K_{\rm II}(a; -x,y) \,.
</math>

<math>F_y</math>的应力强度因子相似：

<math>
   K_{\rm I}(-a; x,y) = K_{\rm I}(a; -x,y) \,,\,\, 
   K_{\rm II}(-a; x,y) = -K_{\rm II}(a; -x,y) \,.
</math>
[[File:LoadedCrackPlate.svg|thumb|203x203px|平面裂纹荷载]]

=== 平面裂纹荷载 ===
如果裂纹受到位于<math>y=0</math>且<math>-a < x < a</math>的点作用力<math>F_y</math>加载，B处的应力强度因子为：<ref name="rooke7">{{cite book|title=Compendium of stress intensity factors|last1=Rooke|first1=D. P.|last2=Cartwright|first2=D. J.|publisher=HMSO Ministry of Defence. Procurement Executive|year=1976}}</ref>

<math>
   K_{\rm I} = \frac{F_y}{2\sqrt{\pi a}}\sqrt{\frac{a+x}{a-x}}\,,\,\,
   K_{\rm II} = -\frac{F_x}{2\sqrt{\pi a}}\left(\frac{\kappa -1}{\kappa+1}\right) \,.
 </math>

如果力在<math>-a < x < a</math>之间均匀分布，则B点的应力强度因子为：

<math>
   K_{\rm I} = \frac{1}{2\sqrt{\pi a}}\int_{-a}^a F_y(x)\,\sqrt{\frac{a+x}{a-x}}\,{\rm d}x\,,\,\,
   K_{\rm II} = -\frac{1}{2\sqrt{\pi a}}\left(\frac{\kappa -1}{\kappa+1}\right)\int_{-a}^a F_y(x)\,{\rm d}x, \,.
 </math>

=== 无限平面：平行裂缝堆叠<ref>{{Cite book|last=Kundu|first=Tribikram|date=2008-01-30|title=Fundamentals of Fracture Mechanics|url=http://dx.doi.org/10.1201/9781439878057|doi=10.1201/9781439878057|isbn=9780429110054}}</ref> ===
如果裂纹间距远远大于裂纹长度（<math>h \gg a</math>），可忽略邻近裂纹之间的相互作用，应力强度因子等于长度为2a的单一裂纹的应力强度因子，裂纹尖端的应力强度因子为：

<math>
  \begin{align}
  K_{\rm I} & = \sigma\sqrt{\pi a}
  \end{align}
</math>

如果裂纹长度远大于间距（<math>a \gg h</math>），则裂纹可以被视为半无限裂纹的堆叠，裂纹尖端应力强度因子为：

<math>
  \begin{align}
  K_{\rm I} & = \sigma\sqrt{h}
  \end{align}
</math>
[[File:CompactTensionSpecimen.svg|thumb|257x257px|用于断裂韧性测试的紧凑拉伸样品]]

=== 紧凑拉伸试样<ref>{{Cite journal |last=范 |first=天佑 |date=1976-07-01 |title=紧凑拉伸试样的K_I |url=http://dx.doi.org/10.1360/csb1976-21-z1-494 |journal=Chinese Science Bulletin |volume=21 |issue=Z1 |doi=10.1360/csb1976-21-z1-494 |issn=0023-074X}}</ref> ===
紧凑拉伸试样裂纹尖端的应力强度因子为：<ref name="bower2" />

<math>
   \begin{align}
   K_{\rm I} & = \frac{P}{B}\sqrt{\frac{\pi}{W}}\left[16.7\left(\frac{a}{W}\right)^{1/2} - 104.7\left(\frac{a}{W}\right)^{3/2}
      + 369.9\left(\frac{a}{W}\right)^{5/2} \right.\\
       & \qquad \left.- 573.8\left(\frac{a}{W}\right)^{7/2} + 360.5\left(\frac{a}{W}\right)^{9/2} \right]
   \end{align}
 </math>

其中<math>P</math>是施加的载荷，<math>B</math>是试样厚度，<math>a</math>是裂纹长度，<math>W</math>是试样的宽度。
[[File:SingleEdgeNotchBending.svg|thumb|176x176px|用于断裂韧性测试的单边缺口弯曲试样（也称为三点弯曲试样）]]

=== 單邊缺口彎曲试样 ===
单边缺口弯曲试样裂纹尖端的应力强度因子为：<ref name="bower2">{{cite book|title=Applied mechanics of solids|author=Bower, A. F.|year=2009|publisher=CRC Press.}}</ref>

<math>
   \begin{align}
   K_{\rm I} & = \frac{4P}{B}\sqrt{\frac{\pi}{W}}\left[1.6\left(\frac{a}{W}\right)^{1/2} - 2.6\left(\frac{a}{W}\right)^{3/2}
      + 12.3\left(\frac{a}{W}\right)^{5/2} \right.\\
       & \qquad \left.- 21.2\left(\frac{a}{W}\right)^{7/2} + 21.8\left(\frac{a}{W}\right)^{9/2} \right]
   \end{align}
 </math>

其中<math>P</math>是施加的载荷，<math>B</math>是试样厚度，<math>a</math>是裂纹长度，<math>W</math>是试样的宽度。

== 相关 ==

* [[断裂力学]]
* [[断裂韧性]]
* [[能量释放率（断裂力学）|能量释放率]]（Strain energy release rate）
* [[J-积分]]（J-integral）
* [[材料失效理论]]（Material failure theory）
* [[帕里斯法则]]

== 参考 ==
[[Category:断裂力学]]