【小呆的力学笔记】弹塑性力学的初步认知六：后继屈服条件

4. 后继屈服条件

4.1 后继屈服条件

上一章节的屈服条件是在当材料未经受任何塑性变形时且在载荷作用下材料第一次进入屈服应该满足的条件（也叫初始屈服条件），那么对于材料进入屈服后继续加载，材料性能会遵循什么样的规则呢？

回顾在材料的单向拉伸力学试验的讨论中，由下图1应力应变曲线，其中下图1左图为多晶材料的应力应变曲线，存在包晶格效应，表现为一个方向强化对应另一个方向的弱化，右图则没有。不管时哪种应力应变曲线，当材料受到载荷作用时，应力应变先线性增加，到达${\sigma }_s$时，材料进入屈服，随着载荷的进一步增加，应力应变曲线从A点到达M点，此时，如果卸去全部载荷，应力应变遵循MO’卸载到应力为零，再重新加载时，材料应力应变就会沿O’MD变化。那么就相当于M点的应力${\sigma }_s^{{}^{\prime }}$变成了新的屈服极限了，称此为后继屈服点。

同时，材料从OAM到M和从O’M到M，M的应力应变状态是一致的，因此材料载荷连续增加从A增加到M，就可以理解为考虑了塑性变形后材料再从O’M增加到M。那么后继屈服条件方程形式与初始屈服条件应为一致，只需把${\sigma }_s$换成${\sigma }_s^{{}^{\prime }}$即可。

$$图1加载/卸载应力应变曲线$$

对于一般应力状态（${\sigma }_{ij}$），在初始屈服之后遵循规则同单向拉伸一致，即在一般应力状态在初始屈服之后同样满足后继屈服条件，与单向应力状态不同的是，此时后继屈服不是一个点，而是一个面。

我们可以这么想，AM阶段是弹性，那么M就是AM的初始屈服条件了，那么显然可以参考初始屈服条件的形式，此时后继屈服条件也可以写成如下形式
$$f({\sigma }_{ij},{\xi }_{\beta })=0\text{\hspace{1em}(1)}$$
称${\xi }_{\beta }$为内变量，与材料的塑性变形历史相关。当${\xi }_{\beta }=0$时，上式就会化为上章节的屈服条件。随着${\xi }_{\beta }$的变化，后继屈服条件的大小、形状和位置都可能会产生变化。因此后继屈服面是一系列的曲面。

当应力位于屈服曲面内时，应力和应变呈现弹性特点，不产生新的塑性变形，因此${\xi }_{\beta }$不变化，当应力位于屈服曲面上
，当载荷继续增加时，应力的变化会增加新的塑性变形，因此会引起内变量的改变，此时会形成新的屈服曲面
$$f({\sigma }_{ij}+d{\sigma }_{ij},{\xi }_{\beta }+d{\xi }_{\beta })=f({\sigma }_{ij},{\xi }_{\beta })+df({\sigma }_{ij}+d{\sigma }_{ij},{\xi }_{\beta }+d{\xi }_{\beta })=0\text{\hspace{1em}(2)}$$

那么，有
$$df({\sigma }_{ij}+d{\sigma }_{ij},{\xi }_{\beta }+d{\xi }_{\beta })=\frac{\partial f}{\partial {\sigma }_{ij}}d{\sigma }_{ij}+\frac{\partial f}{\partial {\xi }_{\beta }}d{\xi }_{\beta }=0\text{\hspace{1em}(3)}$$
这就是屈服条件的一致性条件。

4.2 强化模型

观察应力应变曲线，材料在超过屈服后有强化阶段，因此后继屈服条件总是与初始屈服条件不一样，为了描述这种材料行为，人们建立了多种强化模型，其中最典型的有等向强化模型和随动强化模型。

4.2.1 等向强化模型

在等向强化模型中，材料进入塑性后在一个方向的应力强化，在相反方向同样进行了强化，同时认为两者的同量强化，即M点和M’点应力绝对值相等。这相当于后继屈服曲面在应力空间中是初始屈服曲面的相似扩大，中心位置和形状都与初始屈服曲面一致。

$$图2等向强化的加载/卸载应力应变曲线$$
回顾初始屈服条件，其方程具备以下形式
$$f({\sigma }_{ij})=f^{\ast }({\sigma }_{ij})-k=0\text{\hspace{1em}(4)}$$
那么等向强化模型的后继屈服条件方程具备以下形式
$$f({\sigma }_{ij},{\xi }_{\beta })=f^{\ast }({\sigma }_{ij})-k({\xi }_{\beta })=0\text{\hspace{1em}(5)}$$
一般常将内变量${\xi }_{\beta }$取为等效塑性应变增量
$$d{\xi }_{\beta }=\sqrt{\frac{2}{3}d{\epsilon }_{kl}^{p}d{\epsilon }_{kl}^p}\text{\hspace{1em}(6)}$$
或者塑性功增量
$$d{\xi }_{\beta }={\sigma }_{kl}d{\epsilon }_{kl}^p=d{W}^p\text{\hspace{1em}(7)}$$
并且，当${\xi }_{\beta }=0$时，式（4）就退化为式（3），即$k({\xi }_{\beta }=0)=k$。

以Mise屈服条件为例，回顾Mise屈服条件，有$\overline{\sigma }-{\sigma }_s=0$，那么对应的后继屈服条件方程为
$$\begin{array}{r}f({\sigma }_{ij},{\xi }_{\beta })=\overline{\sigma }-k({\xi }_{\beta })=0\\ k({\xi }_{\beta }=0)={\sigma }_s\end{array}\text{\hspace{1em}(8)}$$

当内变量${\xi }_{\beta }$取为等效塑性应变增量，$k({\xi }_{\beta })$可以有单向拉伸试验确定。事实上当
$$d{\xi }_{\beta }=d{\epsilon }_1^p$$
在单向拉伸试验中，主应力方向1为拉伸方向，在试验过程中始终不变，因此
$${\xi }_{\beta }=\int d{\epsilon }_1^p={\epsilon }_1^p$$
$k({\xi }_{\beta })=k({\epsilon }_1^p)$确定过程如下图3，下图3左图为单向拉伸应力应变曲线，下图3右图为单向拉伸应力应变曲线去除弹性应变分量，其中在图3中左图切线斜率就是模量E’，图3中右图切线斜率就是前文中提到过的h，不难知道其实$k({\epsilon }_1^p)=h$。

$$图3单向拉伸试验应力应变曲线$$

4.2.2 随动强化模型

等向强化模型不符合很多材料的力学特征，很多材料存在一个方向上强化，其相反方向弱化的特点（即包氏效应），因此对此提出随动强化模型。随动强化模型假设材料一个方向上强化，其相反方向弱化同等数量，即如下图所示，其后继屈服点M与反方向后继屈服点M’的差值保持不变，等于2${\sigma }_s$，那么不难发现其实在应力空间中，就是初始屈服曲面的位置平动形成一系列随动强化模型下的后继屈服曲面。

$$图4随动强化的加载/卸载应力应变曲线$$

那么后继屈服条件的方程为
$$f({\sigma }_{ij},{\xi }_{\beta })=f^{\ast }({\sigma }_{ij}-{\alpha }_{ij})-k=0\text{\hspace{1em}(9)}$$
其中${\alpha }_{ij}$为内变量${\xi }_{\beta }$的分量形式，当${\alpha }_{ij}=0$时，上式退化为初始屈服条件方程，因此其意义就是表示屈服曲面中心平移的度量。不同的${\alpha }_{ij}$演化规律又对应不同的随动强化模型，典型如下
（1）线性随动强化模型
$${\dot{\alpha }}_{ij}=c{\dot{\epsilon }}_{ij}^p\text{\hspace{1em}(10)}$$
式中c是一个材料相关常数。式（9）其实也可以写成增量的形式
$$d{\alpha }_{ij}=cd{\epsilon }_{ij}^p\text{\hspace{1em}(10’)}$$
观察式（10），此式中为各分量形式，因此能够直接积分（塑性等效应变不能直接积分，因为在变形过程中方向会有变化），得到
$${\alpha }_{ij}=c{\epsilon }_{ij}^p\text{\hspace{1em}(11)}$$
代入式（9），那么有
$$f({\sigma }_{ij},{\xi }_{\beta })=f^{\ast }({\sigma }_{ij}-c{\epsilon }_{ij}^p)-k=0\text{\hspace{1em}(12)}$$
以Mise屈服条件为例，其方程为
$$f({\sigma }_{ij},{\xi }_{\beta })=f^{\ast }({\sigma }_{ij}-c{\epsilon }_{ij}^p)-k=\overline{\sigma }({\sigma }_{ij}-c{\epsilon }_{ij}^p)-k=0\text{\hspace{1em}(13)}$$
其中等效应力$\overline{\sigma }$如下所示

$$\overline{\sigma }({\sigma }_{ij}-c{\epsilon }_{ij}^p)=\sqrt{\frac{3}{2}(s_{ij}-c{\epsilon }_{ij}^p)(s_{ij}-c{\epsilon }_{ij}^p)}\text{\hspace{1em}(14)}$$
其中$k$就是${\sigma }_s$。

在简单拉伸时，${\sigma }_{ij}、s_{ij}、{\epsilon }_{ij}^p$如下非零项
$$\begin{gathered} {\sigma }_{11}=\sigma ,s_{11}={\sigma }_{11}-{\sigma }_m=\frac{2}{3}\sigma ,s_{22}=s_{33}=-\frac{1}{3}\sigma ; \\ {\epsilon }_{11}^p={\epsilon }^p,{\epsilon }_{22}^p={\epsilon }_{33}^p=-\frac{1}{2}{\epsilon }^p(塑性时泊松比为\frac{1}{2}) \end{gathered}$$
代入式（14），有
$$\begin{array}{rl}\overline{\sigma }& =\sqrt{\frac{3}{2}(s_{ij}-c{\epsilon }_{ij}^p)(s_{ij}-c{\epsilon }_{ij}^p)}\\ & =\sqrt{\frac{3}{2}[(s_{11}-c{\epsilon }_{11}^p{)}^2+(s_{22}-c{\epsilon }_{22}^p{)}^2+(s_{33}-c{\epsilon }_{33}^p{)}^2]}\\ & =\sqrt{\frac{3}{2}[(\frac{2}{3}\sigma -c{\epsilon }^p{)}^2+(-\frac{1}{3}\sigma +c\cdot \frac{1}{2}{\epsilon }^p{)}^2+(-\frac{1}{3}\sigma +c\cdot \frac{1}{2}{\epsilon }^p{)}^2]}\\ & =\sqrt{\frac{3}{2}(\frac{2}{3}{\sigma }^2-2\sigma c{\epsilon }^p+\frac{3}{2}{c}^2\cdot ({\epsilon }^p{)}^2)}\\ & =\sqrt{\frac{3}{2}\cdot \frac{2}{3}({\sigma }^2-3\sigma c{\epsilon }^p+\frac{9}{4}{c}^2\cdot ({\epsilon }^p{)}^2)}\\ & =\sigma -\frac{3}{2}c\cdot {\epsilon }^p\end{array}\text{\hspace{1em}(15)}$$
将上式代入式（13），有
$$f({\sigma }_{ij},{\xi }_{\beta })=\overline{\sigma }({\sigma }_{ij}-c{\epsilon }_{ij}^p)-k=\sigma -\frac{3}{2}c\cdot {\epsilon }^p-{\sigma }_s=0\text{\hspace{1em}(16)}$$
即
$$\sigma ={\sigma }_s+\frac{3}{2}c\cdot {\epsilon }^p\text{\hspace{1em}(17)}$$

前面章节，由简单拉伸的应力应变曲线
$$\sigma ={\sigma }_s+h{\epsilon }^p\text{\hspace{1em}(18)}$$
得出$c=\frac{2}{3}h$。

（2）Ziegler模型

$${\dot{\alpha }}_{ij}=\dot{\eta }({\sigma }_{ij}-{\alpha }_{ij})(\dot{\eta }>0)\text{\hspace{1em}(19)}$$
其增量型方程为
$$d{\alpha }_{ij}=d\eta \cdot ({\sigma }_{ij}-{\alpha }_{ij})(d\eta >0)\text{\hspace{1em}(19’)}$$
上式代入一致性条件，可得
$$\begin{array}{rl}& \frac{\partial f}{\partial ({\sigma }_{ij}-{\alpha }_{ij})}\cdot \frac{\partial ({\sigma }_{ij}-{\alpha }_{ij})}{\partial {\sigma }_{ij}}d{\sigma }_{ij}+\frac{\partial f}{\partial ({\sigma }_{ij}-{\alpha }_{ij})}\cdot \frac{\partial ({\sigma }_{ij}-{\alpha }_{ij})}{\partial {\alpha }_{ij}}d{\alpha }_{ij}\\ =& \frac{\partial f}{\partial ({\sigma }_{ij}-{\alpha }_{ij})}\cdot 1\cdot d{\sigma }_{ij}+\frac{\partial f}{\partial ({\sigma }_{ij}-{\alpha }_{ij})}\cdot -1\cdot d{\alpha }_{ij}=0\end{array}\text{\hspace{1em}(20)}$$
随动强化模型后继屈服条件的方程有
$$f({\sigma }_{ij},{\xi }_{\beta })=f^{\ast }({\sigma }_{ij}-{\alpha }_{ij})-k=0$$
因此，有
$$\frac{\partial f}{\partial ({\sigma }_{ij}-{\alpha }_{ij})}=\frac{\partial f}{\partial {\sigma }_{ij}}=\frac{\partial f^{\ast }}{\partial {\sigma }_{ij}}\text{\hspace{1em}(21)}$$

将上式代入式（20），有
$$\frac{\partial f^{\ast }}{\partial {\sigma }_{ij}}d{\sigma }_{ij}-\frac{\partial f^{\ast }}{\partial {\sigma }_{ij}}d{\alpha }_{ij}=\frac{\partial f^{\ast }}{\partial {\sigma }_{ij}}(d{\sigma }_{ij}-d{\alpha }_{ij})=0\text{\hspace{1em}(22)}$$
将式（19’）代入上式，则
$$\frac{\partial f^{\ast }}{\partial {\sigma }_{ij}}(d{\sigma }_{ij}-d{\alpha }_{ij})=\frac{\partial f^{\ast }}{\partial {\sigma }_{ij}}[d{\sigma }_{ij}-d\eta \cdot ({\sigma }_{ij}-{\alpha }_{ij})]=0\text{\hspace{1em}(23)}$$

$$d\eta =\frac{\partial f^{\ast }}{\partial {\sigma }_{ij}}d{\sigma }_{ij}/[\frac{\partial f^{\ast }}{\partial {\sigma }_{kl}}({\sigma }_{kl}-{\alpha }_{kl})]\text{\hspace{1em}(24)}$$
那么，${\alpha }_{ij}=\int d\eta \cdot ({\sigma }_{ij}-{\alpha }_{ij})$
同样将上式代入式（13），就可得到相应的后继屈服条件方程
$$f({\sigma }_{ij},{\xi }_{\beta })=\overline{\sigma }({\sigma }_{ij}-{\alpha }_{ij})-k=0$$
具体形式先不展开了，因为在暂时用不上。

4.2.3 两种强化模型的讨论

在$\pi$平面上，两种模型的形象比较，下图中C1为初始Mise屈服条件的屈服曲线，C2为等向强化模型的初始Mise屈服条件的屈服曲线，C3为随动强化模型的初始Mise屈服条件的屈服曲线，其中$\alpha$为${\alpha }_{ij}$。

图中可以看出等向强化模型对应屈服面等比例扩大，但是中心位置不变，随动强化对应的是中心位置的移动，但是半径不变。

$$图5两种强化模型的屈服曲线（\pi 平面）$$

Ivey曾对铝合金薄圆管进行了拉扭实验，但随着剪应力的增加，整个后继屈服面都向剪应力增加的方向移动，接近于随动强化模型。

$$图6实验的后继屈服曲线$$