位错的应力场

位错的应力场

假设

位错周围的弹性应力场弹性体假设模型：

1. 晶体是完全弹性体；
2. 晶体是各向同性的；
3. 晶体中没有空隙，由连续介质组成。

柱坐标系下的线弹性力学

平衡方程

$$\begin{array}{r}\begin{array}{rl}& \frac{\partial {\sigma }_{rr}}{\partial r}+\frac{1}{r}\frac{\partial {\sigma }_{\theta r}}{\partial \theta }+\frac{\partial {\sigma }_{zr}}{\partial z}+\frac{{\sigma }_{rr}-{\sigma }_{\theta \theta }}{r}=0\\ & \frac{\partial {\sigma }_{r\theta }}{\partial r}+\frac{1}{r}\frac{\partial {\sigma }_{\theta \theta }}{\partial \theta }+\frac{\partial {\sigma }_{z\theta }}{\partial z}+\frac{2{\sigma }_{r\theta }}{r}=0\\ & \frac{\partial {\sigma }_{rz}}{\partial r}+\frac{1}{r}\frac{\partial {\sigma }_{\theta z}}{\partial \theta }+\frac{\partial {\sigma }_{zz}}{\partial z}+\frac{{\sigma }_{rz}}{r}=0\end{array}\}\end{array}$$

几何方程

$$\begin{array}{r}\begin{array}{rl}& {\epsilon }_{rr}=\frac{\partial u_r}{\partial r}\\ & {\epsilon }_{\theta \theta }=\frac{u_r}{r}+\frac{1}{r}\frac{\partial u_{\theta }}{\partial \theta }\\ & {\epsilon }_{zz}=\frac{\partial u_z}{\partial z}\\ & {\epsilon }_{r\theta }={\epsilon }_{\theta r}=\frac{1}{2}(\frac{1}{r}\frac{\partial u_r}{\partial \theta }+\frac{\partial u_{\theta }}{\partial r}-\frac{u_{\theta }}{r})\\ & {\epsilon }_{rz}={\epsilon }_{zr}=\frac{1}{2}(\frac{\partial u_r}{\partial z}+\frac{\partial u_z}{\partial r})\\ & {\epsilon }_{z\theta }={\epsilon }_{\theta z}=\frac{1}{2}(\frac{\partial u_{\theta }}{\partial z}+\frac{1}{r}\frac{\partial u_z}{\partial \theta })\end{array}\}\end{array}$$

物理方程

$$\begin{array}{r}\left(\begin{array}{c}{\epsilon }_{rr}\\ {\epsilon }_{\theta \theta }\\ {\epsilon }_{zz}\\ {\epsilon }_{r\theta }\\ {\epsilon }_{rz}\\ {\epsilon }_{z\theta }\end{array}\right)=\left(\begin{array}{cccccc}\frac{1}{E}& -\frac{\nu }{E}& -\frac{\nu }{E}& 0& 0& 0\\ -\frac{\nu }{E}& \frac{1}{E}& -\frac{\nu }{E}& 0& 0& 0\\ -\frac{\nu }{E}& -\frac{\nu }{E}& \frac{1}{E}& 0& 0& 0\\ 0& 0& 0& \frac{1}{2G}& 0& 0\\ 0& 0& 0& 0& \frac{1}{2G}& 0\\ 0& 0& 0& 0& 0& \frac{1}{2G}\end{array}\right)\left(\begin{array}{c}{\sigma }_{rr}\\ {\sigma }_{\theta \theta }\\ {\sigma }_{zz}\\ {\sigma }_{r\theta }\\ {\sigma }_{rz}\\ {\sigma }_{z\theta }\end{array}\right)\end{array}$$

刃型位错的应力场

建立如图1所示的刃型位错力学模型。该模型中圆筒的轴线对应刃位错的位错线，圆筒的空心部分相当于位错的中心区。该模型的位移场非常的复杂，目前我看网上也没有相关的推导，甚至连应变场也没有。查阅文献后才能给出其位移场和应变场，具体的文献，可以查阅这篇：Physical Metallurgy.http://dx.doi.org/10.1016/B978-0-444-53770-6.00016-2。

位移场公式如下：
$$\begin{array}{r}\begin{array}{rl}& u_r=-\frac{b}{2\pi }\left(\frac{1-2\nu }{2\left(1-\nu \right)}\sin \theta \ln r-\frac{\sin \theta }{4\left(1-\nu \right)}-\theta \cos \theta \right),\\ & u_{\theta }=-\frac{b}{2\pi }\left(\frac{1-2\nu }{2\left(1-\nu \right)}\cos \theta \ln r+\frac{\cos \theta }{4\left(1-\nu \right)}+\theta \sin \theta \right).\end{array}\}\end{array}$$

根据公式2我们可以容易地得到应变场：
$$\begin{array}{r}\begin{array}{rl}& {\epsilon }_{rr}={\epsilon }_{\theta \theta }=-\frac{b(1-2\nu )}{4\pi (1-\nu )}\frac{\sin \theta }{r},\\ & {\epsilon }_{r\theta }={\epsilon }_{\theta r}=\frac{b}{4\pi (1-\nu )}\frac{\cos \theta }{r},\\ & {\epsilon }_{zz}={\epsilon }_{z\theta }={\epsilon }_{\theta z}={\epsilon }_{rz}={\epsilon }_{zr}=0.\end{array}\}\end{array}$$

最后由公式3给出刃位错在柱坐标系下的应力场公式：
$$\begin{array}{r}\begin{array}{rl}& {\sigma }_{rr}={\sigma }_{\theta \theta }=-\frac{Gb}{2\pi (1-\nu )}\frac{\sin \theta }{r}\\ & {\sigma }_{zz}=-\frac{Gb\nu }{\pi (1-\nu )}\frac{\sin \theta }{r}\\ & {\sigma }_{r\theta }={\sigma }_{\theta r}=\frac{Gb}{2\pi (1-\nu )}\frac{\cos \theta }{r}\\ & {\sigma }_{z\theta }={\sigma }_{\theta z}={\sigma }_{rz}={\sigma }_{zr}=0\end{array}\}\end{array}$$

通过以上公式，我们可以得出以下结论：

1. 径向应力$({\sigma }_{rr})$与切向应力$({\sigma }_{\theta \theta })$大小相同；
2. 当$\theta =0$时有，力${\sigma }_{rr}={\sigma }_{\theta \theta }={\sigma }_{zz}=0$，${\sigma }_{r\theta }={\sigma }_{\theta r}=\frac{Gb}{2\pi (1-\nu )r}$。

螺型位错的应力场

考虑一个长的直线状螺型位错，如图2所示。令圆柱坐标系统的$z$轴沿着位错线的方向。因为在离开位错芯足够远的区域里胡克定律是适用的，所以它的位移场可以用弹性理论来描写。

微元体在径向和切向是没有位移的，即$u_r=u_{\theta }=0$，轴向位移是方位角$\theta$的函数。每经过一个完全的回路$(\theta =2\pi )$就得到一个大小等于柏氏矢量$b$的轴向位移。螺型位错的轴向位移分量由下式给出：
$$\begin{array}{r}{u}_z=\frac{b\theta }{2\pi }\end{array}$$

因此，根据公式2，我们可以看到在螺型位错中只有变应分量${ϵ}_{z\theta }$不为零， 并由下式给出：
$$\begin{array}{r}{\epsilon }_{z\theta }={\epsilon }_{\theta z}=\frac{1}{2}\frac{1}{r}\frac{\partial u_x}{\partial \theta }=\frac{b}{4\pi r}\end{array}$$

下面根据公式3给出螺型位错在柱坐标系下的应力场公式：
$$\begin{array}{r}\begin{array}{rl}& {\sigma }_{z\theta }={\sigma }_{\theta z}=\frac{Gb}{2\pi r}\\ & {\sigma }_{rr}={\sigma }_{\theta \theta }={\sigma }_{zz}={\sigma }_{r\theta }={\sigma }_{\theta r}={\sigma }_{rz}={\sigma }_{zr}=0\end{array}\}\end{array}$$

通过以上公式，我们可以得出以下结论：

1. 螺型位错周围是简单的纯剪切，应变具有径向对称性；
2. 在包含位错线的任何晶向平面上剪应力都是$Gb/2\pi r$，与$\theta$角无关；
3. 螺型位错的应力场可用位错周围一定尺寸的圆柱体表示。

不过，你以为就这么简单了吗？那你就想错了以上情况是对于无限长圆柱体的。有限长螺型位错${\sigma }_{z\theta }$应力分量在圆柱体的端面产生扭矩项，使得自由表面处的边界条件：
$$\begin{array}{r}{\sigma }_{ij}\mathrm{d}{S}_j=0\end{array}$$
没有实现，因此，应力分量${\sigma }_{z\theta }$和${\sigma }_{\theta z}$必须通过确保圆柱体端面的自由应力状态的项来校正。读者可以自行尝试一下公式推导，在这篇文章接下来的内容中都默认模型是无限长的。

位错的弹性能

位错周围点阵畸变引起弹性应力场导致晶体能量的增加，这部分能量即为位错的应变能。包括两部分：$E_{total}=E_{core}+E_{el}$。
\begin{enumerate}
\item 位错核心能$E_{core}$：在位错核心几个原子间距$r_o=2\left|b\right|=2b$以内的区域，滑移面两侧原子间的错排能即相当于位错核心能。错排能约占位错能的$1/10$，可忽略；
\item 弹性应变能$E_{el}$：在位错核心区以外，长程应力场作用范围所具有的能量，约占位错能的$9/10$。
\end{enumerate}

在掌握了应力和应变场的情况下，可以使用公式11计算与螺型位错和刃型位错相关的弹性能。
$$\begin{array}{r}{E}_{el}=\iiint w\mathrm{d}V\end{array}$$
其中：
$$\begin{array}{r}w=\frac{1}{2}{\sigma }_{ij}{\epsilon }_{ij}\end{array}$$

如果这样做遇到了困难，我们还可以使用做功法计算：
$$\begin{array}{r}\mathrm{d}{E}_{el}={\sigma }_{ij}(2\pi r\mathrm{d}x)\mathrm{d}{b}^{\prime }\end{array}$$
\subsection{刃型位错的应变能}
对于刃型位错（如图3），其应变能是由${\sigma }_{rr},{\sigma }_{\theta \theta },{\sigma }_{r\theta },{\sigma }_{\theta r}$贡献的，虽然有${\sigma }_{zz}$，但是其并没有导致应变，所以对应变能没有贡献。应变能积分表达式如下：

$$\begin{array}{rl}{E}_{Edge}& =\iiint \left(\frac{1}{2}{\sigma }_{rr}{\epsilon }_{rr}+\frac{1}{2}{\sigma }_{\theta \theta }{\epsilon }_{\theta \theta }+\frac{1}{2}{\sigma }_{r\theta }{\epsilon }_{r\theta }+\frac{1}{2}{\sigma }_{\theta r}{\epsilon }_{\theta r}\right)\mathrm{d}V\\ & ={\int }_0^1\mathrm{d}z{\int }_0^{2\pi }r\mathrm{d}\theta {\int }_{r_0}^R\left(\frac{1}{2}{\sigma }_{rr}{\epsilon }_{rr}+\frac{1}{2}{\sigma }_{\theta \theta }{\epsilon }_{\theta \theta }+\frac{1}{2}{\sigma }_{r\theta }{\epsilon }_{r\theta }+\frac{1}{2}{\sigma }_{\theta r}{\epsilon }_{\theta r}\right)\mathrm{d}r\end{array}$$

将公式(5),(6)带入括号中，可以计算得到刃型位错的弹性能量密度为：
$$\begin{array}{rl}\sum w& =\frac{G{b}^2}{8{\pi }^2(1-\nu {)}^2}\frac{1-2\nu {\sin }^2\theta }{r^2}\\ & =\frac{G{b}^2}{8{\pi }^2(1-\nu {)}^2}\frac{1-\nu +\nu \cos 2\theta }{r^2}\end{array}$$
将其带入公式(15)就可以得到刃型位错弹性能表达式：
$$\begin{array}{rl}{E}_{Edge}& =\frac{G{b}^2}{8{\pi }^2(1-\nu {)}^2}{\int }_0^1\mathrm{d}z{\int }_0^{2\pi }r\mathrm{d}\theta {\int }_{r_0}^R\left(\frac{1-\nu +\nu \cos 2\theta }{r^2}\right)\mathrm{d}r\\ & =\frac{G{b}^2}{8{\pi }^2(1-\nu {)}^2}{\int }_0^{2\pi }(1-\nu +\nu \cos 2\theta )\mathrm{d}\theta {\int }_{r_0}^R\frac{1}{r}\mathrm{d}r\\ & =\frac{G{b}^2}{8{\pi }^2(1-\nu {)}^2}2\pi (1-\nu )\ln \left(\frac{R}{r_0}\right)\\ & =\frac{G{b}^2}{4\pi (1-\nu )}\ln \left(\frac{R}{r_0}\right)\end{array}$$

螺型位错的应变能

螺型位错应变能的推导就更加简单了如图4，读者可以自行尝试，这里给出其结果：

$$\begin{array}{r}{E}_{Spiral}=\frac{\mu b^2}{4\pi }\ln \frac{R}{r_o}\end{array}$$
\subsection{混合位错的弹性能}
值得注意的是，如公式(5)(6)(8)(9)所述，螺型位错和刃型位错的应力和应变张量没有公共项，即两个张量是正交的。因此，在各向同性线性弹性的框架中，叠加原理适用，混合位错的应力和应变场通过独立计算与柏氏矢量的螺型和刃型分量相关联的应力和应变场，并将两者的贡献相加而获得。
$$\begin{array}{r}{W}_{\text{mix }}=\left[\frac{G(b\sin \alpha {)}^2}{4\pi (1-\nu )}+\frac{G(b\cos \alpha {)}^2}{4\pi }\right]\ln \left(\frac{R}{{\mathrm{r}}_0}\right)\end{array}$$

张力

张力是存在于面/线缺陷的表面/界面中的，下面先给出几个概念：

• 表面能$(\mathrm{J}/{\mathrm{m}}^2)$：单位面积表面原子相对于内部原子的能量增加。
• 表面张力$(\mathrm{N}/\mathrm{m})$：表面内使单位长度表面收缩的力。
• 位错应变能$(\mathrm{J}/\mathrm{m})$：单位长度位错(线)相对于完美晶体的能量增加。
• 位错线张力$(\mathrm{N})$：位错线内部使位错线长度收缩的力。

表面能和表面张力、位错应变能和位错线张力是同一事物的不同一侧面，两者的量纲与数值必然相同。所以以上求得的应变能也就是位错线张力。
位错线张力一段曲线状的位错线受一个指向曲率中心的力。图5表示一个长度为$S$，曲率半径为$R$的位错线。作用在这段位错线上的径向力是

$$\begin{array}{r}\frac{\text{力}}{\text{单位长度}}=\frac{2\tau \sin \theta }{S}\simeq \frac{2\tau \theta }{S}=\frac{2\tau }{S}\frac{S}{2R}=\frac{\tau }{R}\end{array}$$

图6绘出了作用在表面张力为$\gamma$的物质的球缺表面上的径向力。作用在球缺单位面积上的力是：

$$\begin{array}{r}\frac{\text{力}}{\text{单位面积}}=\frac{2\pi r\gamma \sin \theta }{A}\simeq \frac{2\pi r\gamma \theta }{\pi r^2}=\frac{2\pi r\gamma }{\pi r^2}\frac{r}{R}=\frac{2\gamma }{R}\end{array}$$

OK！搞定！