【小呆的力学笔记】弹塑性力学的初步认知五：初始屈服条件（3）

上文中，提到了三种屈服条件：1. Tresca屈服条件（即第三强度理论），2. Mises屈服条件（即第四强度理论），3. 最大偏应力屈服条件（即第一强度理论的修正版），那么这三种屈服条件哪个更能贴近物理实际呢？历史上，有两个著名的验证屈服条件的试验，具体如下。

3.7 屈服条件的实验验证一：薄圆管受拉力和内压的联合作用（Lode，1926年）

试验的示意图如下图7所示，薄壁圆管的半径为$R$，壁厚为$h$，并且$h\ll R$(至多1/20)，其中薄壁圆管受内压$p$和端面拉力T的联合作用，那么，取薄壁微元做应力分析如图8。

$$图7试验示意图$$

$$图8薄壁微元示意图$$

由于薄壁圆管内外表面是自由表面，所以${\tau }_{rz}={\tau }_{r\theta }=0$。
极坐标系下的圆筒平衡方程如下式
$$\begin{array}{rl}& \frac{\partial {\sigma }_r}{\partial r}+\frac{1}{r}\frac{\partial {\tau }_{r\theta }}{\partial \theta }+\frac{{\sigma }_r-{\sigma }_{\theta }}{r}+f_r=0\\ & \frac{1}{r}\frac{\partial {\sigma }_{\theta }}{\partial \theta }+\frac{\partial {\tau }_{r\theta }}{\partial r}+\frac{2{\tau }_{r\theta }}{r}+f_{\theta }=0\end{array}\text{\hspace{1em}(40)}$$
当体力为零时，用应力函数$\varphi$表示的应力
$$\begin{array}{rl}& {\sigma }_r=\frac{1}{r}\frac{\partial \varphi }{\partial r}+\frac{1}{r^2}\frac{{\partial }^2\varphi }{\partial {\theta }^2}\\ & {\sigma }_{\theta }=\frac{{\partial }^2\varphi }{\partial r^2}\\ & {\tau }_{r\theta }=-\frac{\partial }{\partial r}(\frac{1}{r}\frac{\partial \varphi }{\partial \theta })=\frac{1}{r^2}\frac{\partial \varphi }{\partial \theta }-\frac{1}{r}\frac{{\partial }^2\varphi }{\partial r\theta }\end{array}\text{\hspace{1em}(41)}$$
应力相容方程：
$$(\frac{{\partial }^2}{\partial r^2}+\frac{1}{r}\frac{\partial }{\partial r}+\frac{1}{r^2}\frac{{\partial }^2}{\partial {\theta }^2}{)}^2\varphi =0\text{\hspace{1em}(42)}$$
在本试验中，结构和载荷都是轴对称的，因此应力函数$\varphi$只是径向坐标$r$的函数，因此式（41）变为
$$\begin{array}{rl}& {\sigma }_r=\frac{1}{r}\frac{\partial \varphi }{\partial r}\\ & {\sigma }_{\theta }=\frac{{\partial }^2\varphi }{\partial r^2}\\ & {\tau }_{r\theta }=0\end{array}\text{\hspace{1em}(43)}$$
同时式（42）变为
$$(\frac{{\partial }^2}{\partial r^2}+\frac{1}{r}\frac{\partial }{\partial r}{)}^2\varphi =0\text{\hspace{1em}(44)}$$
上式为四节常微分方程，通解为
$$\varphi =A\ln r+B{r}^2\ln r+C{r}^2+D\text{\hspace{1em}(45)}$$
应力分量即为
$$\begin{array}{rl}& {\sigma }_r=\frac{A}{r^2}+B(1+2\ln r)+2C\\ & {\sigma }_{\theta }=-\frac{A}{r^2}+B(3+2\ln r)+2C\\ & {\tau }_{r\theta }=0\end{array}\text{\hspace{1em}(46)}$$
极坐标下的几何方程和物理方程为
$$几何方程：\begin{array}{rl}& {\epsilon }_r=\frac{\partial u_r}{\partial r}\\ & {\epsilon }_{\theta }=\frac{u_r}{r}+\frac{1}{r}\frac{\partial u_{\theta }}{\partial \theta }\\ & {\gamma }_{r\theta }=\frac{1}{r}\frac{\partial u_r}{\partial \theta }+\frac{\partial u_{\theta }}{\partial r}-\frac{u_{\theta }}{r}\end{array}\text{\hspace{1em}(47)}$$

$$物理方程：\begin{array}{rl}& {\epsilon }_r=\frac{1}{E}({\sigma }_r-\mu {\sigma }_{\theta })\\ & {\epsilon }_{\theta }=\frac{1}{E}({\sigma }_{\theta }-\mu {\sigma }_r)\\ & {\gamma }_{r\theta }=\frac{2(1+\mu )}{E}{\tau }_{r\theta }\end{array}\text{\hspace{1em}(48)}$$

并且联立几何方程和物理方程，可得位移分量的形式解：
$$\begin{array}{rl}{u}_r& =\frac{1}{E}[-\frac{A}{r}(1+\mu )+2(1-mu)Br(\ln r-1)+(1-3\mu )Br+2Cr(1-\mu )]+I\cos \theta +K\sin \theta \\ u_{\theta }& =\frac{4Br\theta }{E}+Hr-I\sin \theta +K\cos \theta \end{array}\text{\hspace{1em}(49)}$$
由位移的单值条件，有
$$\begin{array}{r}{u}_0=u_{2\pi }=0\to \frac{4Br(0-2\pi )}{E}=0\end{array}\text{\hspace{1em}(50)}$$
那么$B=0$。
应力边界条件为
$$\begin{array}{r}({\sigma }_r{)}_{r=R}=-p,({\sigma }_r{)}_{r=R+h}=0\end{array}\text{\hspace{1em}(51)}$$

不难求得
$${\sigma }_r=-\frac{\frac{R^2}{r^2}-1}{\frac{(R+h{)}^2}{R^2}-1}p,{\sigma }_{\theta }=\frac{\frac{R^2}{r^2}+1}{\frac{(R+h{)}^2}{R^2}-1}p\text{\hspace{1em}(52)}$$
我们来看${\sigma }_r、{\sigma }_{\theta }$两者的大小对比
$$\begin{array}{rl}|\frac{{\sigma }_{\theta }}{{\sigma }_r}|& =\frac{\frac{R^2}{r^2}+1}{\frac{(R+h{)}^2}{R^2}-1}\cdot \frac{\frac{(R+h{)}^2}{R^2}-1}{\frac{R^2}{r^2}-1}\\ & =\frac{\frac{R^2}{r^2}+1}{1-\frac{R^2}{r^2}}\\ & =\frac{R^2+r^2}{r^2-R^2}=\frac{1+(1+\alpha {)}^2}{(1+\alpha {)}^2-1}=\frac{2+2\alpha +{\alpha }^2}{2\alpha +{\alpha }^2}\approx \frac{2+2\alpha }{2\alpha }=1+\frac{1}{\alpha }\end{array}\text{\hspace{1em}(53)}$$
$\alpha \ll 1$，因此可以认为${\sigma }_r\approx 0$

$$\begin{array}{rl}{\sigma }_{\theta }& =\frac{\frac{R^2}{r^2}+1}{\frac{(R+h{)}^2}{R^2}-1}p\\ & =\frac{R^2(\frac{R^2}{r^2}+1)}{2Rh+h^2}p\\ & =\frac{R^2(\frac{R^2}{(1+\alpha {)}^2{R}^2}+1)}{2R\cdot \beta R+(\beta R{)}^2}p\\ & =\frac{\frac{1}{(1+\alpha {)}^2}+1}{2\beta +{\beta }^2}p\\ & \approx \frac{1}{\beta }p=\frac{R}{h}p\end{array}\text{\hspace{1em}(54)}$$

同时，薄壁圆管端面受拉力T，那么根据力平衡可得
$${\sigma }_z\cdot [\pi (R+h{)}^2-\pi R^2]={\sigma }_z\cdot (2\pi Rh+\pi h^2)\approx {\sigma }_z\cdot 2\pi Rh=T\text{\hspace{1em}(55)}$$
因此
$$\begin{array}{rl}& {\sigma }_{\theta }\approx \frac{PR}{h}\\ & {\sigma }_z\approx \frac{T}{2\pi Rh}\end{array}\text{\hspace{1em}(56)}$$

假设载荷$p$和T满足${\sigma }_{\theta }\ge {\sigma }_z\ge {\sigma }_r$，则主应力为
$${\sigma }_1={\sigma }_{\theta },{\sigma }_2={\sigma }_z,{\sigma }_3={\sigma }_r=0\text{\hspace{1em}(57)}$$

则有
$${\mu }_{\sigma }=\frac{2{\sigma }_2-{\sigma }_1-{\sigma }_3}{{\sigma }_1-{\sigma }_3}=\frac{T-\pi R^{2}p}{\pi R^{2}p}\text{\hspace{1em}(58)}$$
当T=0时，${\mu }_{\sigma }=-1$，此时${\theta }_{\sigma }=-30^{。}$，此时即为下图红圈，对应的就是材料的单向拉伸试验。
当$T=\pi R^{2}p$时，${\mu }_{\sigma }=0$，此时${\theta }_{\sigma }=0^{。}$，此时即为下图x轴与屈服截面曲线的交点（绿色圈代表此时三种屈服条件）。
此时，${\sigma }_{\theta }=\frac{pR}{h},{\sigma }_z=\frac{pR}{2h},{\sigma }_r=0$，如果减去静水应力$\frac{pR}{2h}$，那么${\sigma }_{\theta }=\frac{pR}{2h},{\sigma }_z=0,{\sigma }_r=-\frac{pR}{2h}$，此时对应的就是材料的纯剪切试验。

当$T=2\pi R^{2}p$时，${\mu }_{\sigma }=1$，此时${\theta }_{\sigma }=30^{。}$此时即为下图黄圈，对应的就是材料的单向压缩试验。

因此只要载荷满足$0\le T\le 2\pi R^{2}p$的范围内改变T和p，就可以得到$-30^{。}\sim 30^{。}$的屈服曲线。

$$图9三种屈服条件$$

Lode曾对铁、铜等材料进行了一系列试验，现在来对比三个屈服条件和试验的符合程度。

不失一般性，假设${\sigma }_1\ge {\sigma }_2\ge {\sigma }_3$，那么对于Tresca屈服条件，有
$$\frac{{\sigma }_1-{\sigma }_3}{{\sigma }_s}=1\text{\hspace{1em}(59)}$$
对于Mises屈服条件$f({\sigma }_{ij})=J_2-k_2^2=\frac{1}{3}{\sigma }_s^2-k_2^2=0$

$${\sigma }_s=\sqrt{3}{k}_2=\sqrt{3}\sqrt{J_2}=\sqrt{3}\cdot \frac{r_{\sigma }}{\sqrt{2}}=\sqrt{\frac{3}{2}}{r}_{\sigma }=\sqrt{\frac{3}{2}}\sqrt{x^2+y^2}\text{\hspace{1em}(60)}$$
其中x、y详见上章节
$${\sigma }_s=\sqrt{\frac{3}{2}}\sqrt{x^2+y^2}=\sqrt{\frac{3}{2}}x\sqrt{1+(\frac{y}{x}{)}^2}=\sqrt{\frac{3}{2}}\cdot \frac{{\sigma }_1-{\sigma }_3}{\sqrt{2}}\cdot \sqrt{1+(\frac{{\mu }_{\sigma }}{\sqrt{3}}{)}^2}$$
那么
$$\frac{{\sigma }_1-{\sigma }_3}{{\sigma }_s}=\frac{2}{\sqrt{3}}\cdot \sqrt{\frac{1}{1+(\frac{{\mu }_{\sigma }}{\sqrt{3}})}}=\frac{2}{\sqrt{3+{\mu }_{\sigma }^2}}\text{\hspace{1em}(61)}$$

最后来看最大偏应力屈服条件，由假设${\sigma }_1\ge {\sigma }_2\ge {\sigma }_3$，可得$s_1\ge s_2\ge s_3$，所以
$${\mu }_{\sigma }=\frac{2{\sigma }_2-{\sigma }_1-{\sigma }_3}{{\sigma }_1-{\sigma }_3}=\frac{3s_2}{s_1-s_3}=\frac{-3(s_1+s_3)}{s_1-s_3}\text{\hspace{1em}(62)}$$
由上式可以求出
$$\frac{s_1}{s_3}=\frac{{\mu }_{\sigma }-3}{{\mu }_{\sigma }+3}\text{\hspace{1em}(63)}$$
上一节内容中以及得到以下结论
$$\begin{gathered} s_1=\frac{2}{3}{\sigma }_s，-1\le {\mu }_{\sigma }\le 0 \\ {s}_3=-\frac{2}{3}{\sigma }_s，0\le {\mu }_{\sigma }\le 1 \end{gathered}$$

在$-1\le {\mu }_{\sigma }\le 0$时，
$$\frac{{\sigma }_1-{\sigma }_3}{{\sigma }_s}=\frac{s_1-s_3}{{\sigma }_s}=\frac{s_1}{{\sigma }_s}(1-\frac{s_3}{s_1})=\frac{s_1}{{\sigma }_s}(1-\frac{{\mu }_{\sigma }+3}{{\mu }_{\sigma }-3})=\frac{4}{3-{\mu }_{\sigma }}\text{\hspace{1em}(64)}$$
在$0\le {\mu }_{\sigma }\le 1$时，
$$\frac{{\sigma }_1-{\sigma }_3}{{\sigma }_s}=\frac{s_1-s_3}{{\sigma }_s}=\frac{s_3}{{\sigma }_s}(\frac{s_1}{s_3}-1)=\frac{s_3}{{\sigma }_s}(\frac{{\mu }_{\sigma }-3}{{\mu }_{\sigma }+3}-1)=\frac{4}{3+{\mu }_{\sigma }}\text{\hspace{1em}(65)}$$
那么，在$-1\le {\mu }_{\sigma }\le 1$时，$$\frac{{\sigma }_1-{\sigma }_3}{{\sigma }_s}=\frac{4}{3+|{\mu }_{\sigma }|}\text{\hspace{1em}(66)}$$

$$图10实验结果以及三种屈服条件$$
将试验结果以及三种屈服条件绘制在一张图，如上图，不难发现实验结果更接近Mises屈服条件。同时可以发现所有结果都在Tresca屈服条件和最大偏应力条件中间，即Tresca屈服条件为下届，最大偏应力条件为上届。

3.8 屈服条件的实验验证二：薄圆管受拉力和扭矩的联合作用（Taylor-Quinney,1931）

试验的示意图如下图11所示，薄壁圆管的半径为$R$，壁厚为$h$，并且$h\ll R$(至多1/20)，其中薄壁圆管受端面拉力T和扭矩M的联合作用。

$$图11试验示意图$$
对薄圆管在扭矩作用下进行的截面受力分析，示意图如下图12，那么不难得出
$${\int }_0^A{\tau }_z\cdot dA\cdot R={\int }_0^{2\pi }{\tau }_{z}R\cdot \frac{(R+h{)}^2-R^2}{2}\cdot d\theta ={\tau }_z\pi R\cdot ((R+h{)}^2-R^2)=M\text{\hspace{1em}(67)}$$
由薄圆管的$h\ll R$(至多1/20)，可得
$${\tau }_z=\frac{M}{\pi R\cdot ((R+h{)}^2-R^2)}=\frac{M}{\pi R\cdot (2Rh+h^2)}\approx \frac{M}{2\pi R^{2}h}\text{\hspace{1em}(68)}$$

$$图12受力分析示意图$$

薄壁圆管端面受拉力T，上文中已得到
$${\sigma }_z\cdot [\pi (R+h{)}^2-\pi R^2]={\sigma }_z\cdot (2\pi Rh+\pi h^2)\approx {\sigma }_z\cdot 2\pi Rh=T$$
因此，在拉力T和扭矩M作用下，应力分量为
$$\begin{array}{rl}& {\sigma }_{\theta z}\approx \frac{M}{2\pi R^{2}h}\\ & {\sigma }_z\approx \frac{T}{2\pi Rh}\end{array}\text{\hspace{1em}(69)}$$

其应力张量近似为
$$\left[\begin{array}{ccc}0& 0& 0\\ 0& 0& {\sigma }_{\theta z}\\ 0& {\sigma }_{\theta z}& {\sigma }_z\end{array}\right]\text{\hspace{1em}(70)}$$
那么其主应力为
$$\begin{array}{rl}& {\sigma }_1=\frac{{\sigma }_z}{2}+\frac{1}{2}\sqrt{{\sigma }_z^2+4{\sigma }_{\theta z}^2}\\ & {\sigma }_2=0\\ & {\sigma }_3=\frac{{\sigma }_z}{2}-\frac{1}{2}\sqrt{{\sigma }_z^2+4{\sigma }_{\theta z}^2}\end{array}\text{\hspace{1em}(71)}$$
那么Lode参数为
$${\mu }_{\sigma }=\frac{2{\sigma }_2-{\sigma }_1-{\sigma }_3}{{\sigma }_1-{\sigma }_3}=\frac{-{\sigma }_z}{\sqrt{{\sigma }_z^2+4{\sigma }_{\theta z}^2}}=\frac{-T}{\sqrt{T^2+4M^2/R^2}}\text{\hspace{1em}(72)}$$

当$M=0$，$T>0$时，这对应简单拉伸实验，此时${\mu }_{\sigma }=-1$。
当$T=0$，$M\ne 0$时，这对应纯剪切实验，此时${\mu }_{\sigma }=0$。
改变T和M的比值，便可得到$-1\le {\mu }_{\sigma }\le 0$的各种应力状态。

Tresca屈服条件为
$${\tau }_{max}=\frac{{\sigma }_1-{\sigma }_3}{2}=\frac{1}{2}\sqrt{{\sigma }_z^2+4{\sigma }_{\theta z}^2}=\frac{{\sigma }_s}{2}\text{\hspace{1em}(73)}$$
上式可变换成
$$(\frac{{\sigma }_z}{{\sigma }_s}{)}^2+4(\frac{{\sigma }_{\theta z}}{{\sigma }_s}{)}^2=1\text{\hspace{1em}(74)}$$

Mises屈服条件为
$$J_2-\frac{{\sigma }_s^2}{3}=\frac{1}{6}[({\sigma }_1-{\sigma }_2{)}^2+({\sigma }_1-{\sigma }_3{)}^2+({\sigma }_2-{\sigma }_3{)}^2]-\frac{{\sigma }_s^2}{3}=0\text{\hspace{1em}(75)}$$
将式（71）代入上式
$$\begin{array}{rl}& \frac{1}{6}[({\sigma }_1-{\sigma }_2{)}^2+({\sigma }_1-{\sigma }_3{)}^2+({\sigma }_2-{\sigma }_3{)}^2]-\frac{{\sigma }_s^2}{3}\\ =& \frac{1}{6}[{\sigma }_1^2+({\sigma }_1-{\sigma }_3{)}^2+{\sigma }_3^2]-\frac{{\sigma }_s^2}{3}\\ =& \frac{1}{6}[2(\frac{{\sigma }_z}{2}{)}^2+2(\frac{1}{2}\sqrt{{\sigma }_z^2+4{\sigma }_{\theta z}^2}{)}^2+(\sqrt{{\sigma }_z^2+4{\sigma }_{\theta z}^2}{)}^2]-\frac{{\sigma }_s^2}{3}\\ =& \frac{1}{6}(2{\sigma }_z^2+6{\sigma }_{\theta z}^2)-\frac{{\sigma }_s^2}{3}=0\end{array}\text{\hspace{1em}(76)}$$
跟Tresca屈服条件一样，式（76）也可以变化为式（74）形式如下
$$(\frac{{\sigma }_z}{{\sigma }_s}{)}^2+3(\frac{{\sigma }_{\theta z}}{{\sigma }_s}{)}^2=1\text{\hspace{1em}(77)}$$

最后看最大偏应力屈服条件，先看主偏应力应为

$$\begin{array}{rl}& s_1=\frac{{\sigma }_z}{6}+\frac{1}{2}\sqrt{{\sigma }_z^2+4{\sigma }_{\theta z}^2}\\ & s_2=-\frac{{\sigma }_z}{3}\\ & s_3=\frac{{\sigma }_z}{6}-\frac{1}{2}\sqrt{{\sigma }_z^2+4{\sigma }_{\theta z}^2}\end{array}\text{\hspace{1em}(78)}$$

那么最大偏应力屈服条件为
$$3s_1=3(\frac{{\sigma }_z}{6}+\frac{1}{2}\sqrt{{\sigma }_z^2+4{\sigma }_{\theta z}^2})=2{\sigma }_s\text{\hspace{1em}(79)}$$
同样，上式也可以化为如下形式
$$\frac{1}{4}(\frac{{\sigma }_z}{{\sigma }_s})+\frac{3}{4}\sqrt{(\frac{{\sigma }_z}{{\sigma }_s}{)}^2+4(\frac{{\sigma }_{\theta z}}{{\sigma }_s}{)}^2}=1\text{\hspace{1em}(80)}$$

将实验结果画在一张图中，可见实验结果更接近Mises屈服条件，最大偏应力屈服条件其次。