晶体塑形自学1—大变形

引言

如下图所示，变形前后构型发生变化，当变化足够下，就可以忽略初始构型和变形后构型的区别，即小变形假设。

而大变形，更准确的说是有限变形，就是在描述应力和应变时考虑构型变化的影响。

小变形假设下，名义应力 $\mathbf{p}$ 约等于 真应力 $\mathbf{\sigma }$
$$\left[\begin{array}{cc}{p}_{11}& 0\\ 0& p_{22}\end{array}\right]\approx \left[\begin{array}{cc}{\sigma }_{11}& 0\\ 0& {\sigma }_{22}\end{array}\right]$$

有限变形下，端面合力保持不变
$$\left[\begin{array}{cc}{p}_{11}& 0\\ 0& p_{22}\end{array}\right]=\left[\begin{array}{cc}1+{ϵ}_{22}& 0\\ 0& 1+{ϵ}_{11}\end{array}\right]\left[\begin{array}{cc}{\sigma }_{11}& 0\\ 0& {\sigma }_{22}\end{array}\right]$$

考虑构型变化后，不能再认为两个应力相同。

应变度量

1、 变形梯度(二维)：$$\mathbf{F}=\left[\begin{array}{cc}\frac{\partial x}{\partial X}& \frac{\partial x}{\partial Y}\\ & \\ \frac{\partial y}{\partial X}& \frac{\partial y}{\partial Y}\end{array}\right]$$
$F$ 的行列式用 $J$ 表示：$J=det(F)$，$J$ 是微元的体积变化（变化后的体积比上初始的体积）。

名义应变：$\mathbf{ϵ}={\mathbf{F}}^{\mathbf{T}}-\mathbf{I}$ (含旋转的话为非对称张量)

名义应变率：
$$\dot{\mathbf{ϵ}}={\dot{\mathbf{F}}}^T=\left[\begin{array}{cc}\frac{\partial v_x}{\partial X}& \frac{\partial v_y}{\partial X}\\ & \\ \frac{\partial v_x}{\partial Y}& \frac{\partial v_y}{\partial Y}\end{array}\right]$$

2、Green应变张量 $E$ 定义：
$$d{s}^2-d{S}^2=2d\mathbf{X}\cdot \mathbf{E}\cdot d\mathbf{X}$$
材料矢量 $d\mathbf{X}$ 长度平方的变化。
\quad
其中：$d{s}^2=d{\mathbf{x}}^T\cdot d\mathbf{x}d{S}^2=d{\mathbf{X}}^T\cdot d\mathbf{X}$
$$\begin{array}{rl}d{\mathbf{x}}^T\cdot d\mathbf{x}-d{\mathbf{X}}^T\cdot d\mathbf{X}& =(\mathbf{F}\cdot d\mathbf{X}{)}^T\cdot (\mathbf{F}\cdot d\mathbf{X})-d{\mathbf{X}}^T\cdot d\mathbf{X}\\ & \\ & =d{\mathbf{X}}^T\cdot {\mathbf{F}}^T\cdot \mathbf{F}\cdot d\mathbf{X}-d{\mathbf{X}}^T\cdot \mathbf{I}\cdot d\mathbf{X}\\ & \\ & =d{\mathbf{X}}^T\cdot ({\mathbf{F}}^T\cdot \mathbf{F}-\mathbf{I})\cdot d\mathbf{X}\end{array}$$
因此 $\mathbf{E}=\frac{1}{2}({\mathbf{F}}^T\cdot \mathbf{F}-\mathbf{I})$ 或 $E_{ij}=\frac{1}{2}(F_{ik}{F}_{kj}-{\delta }_{ij})$

3、 变形率

变形率 $D$ ,也称为速度应变，是变形的率度量。

首先，定义速度梯度 $\mathbf{L}$:
$$\mathbf{L}=\frac{\partial \mathbf{v}}{\partial x}={\left(▽\mathbf{v}\right)}^T={\left(grad\mathbf{v}\right)}^T或L_{ij}=\frac{\partial v_i}{\partial x_j}$$
$$\mathbf{L}=\frac{\partial \mathbf{v}}{\partial x}=\frac{\partial \mathbf{v}}{\partial \mathbf{X}}\frac{\partial \mathbf{X}}{\partial \mathbf{x}}=\dot{\mathbf{F}}\cdot {\mathbf{F}}^{-1}或L_{ij}={\dot{F}}_{ik}{F}_{kj}^{-1}$$
速度梯度张量可分解为对称部分 $\mathbf{D}$ 和反对称部分 $\mathbf{W}$
$$\mathbf{D}=\frac{1}{2}\left(\mathbf{L}+{\mathbf{L}}^T\right)或D_{ij}=\frac{1}{2}\left(\frac{\partial v_i}{\partial x_j}+\frac{\partial v_j}{\partial x_i}\right)$$
$$\mathbf{W}=\frac{1}{2}\left(\mathbf{L}-{\mathbf{L}}^T\right)或W_{ij}=\frac{1}{2}\left(\frac{\partial v_i}{\partial x_j}-\frac{\partial v_j}{\partial x_i}\right)$$
变形率是微小材料线段平方的变化率的度量：
$$\begin{array}{rl}\frac{\partial }{\partial t}(d{s}^2)& =\frac{\partial }{\partial t}(d{\mathbf{x}}^T\cdot d\mathbf{x})=(d{\mathbf{v}}^T\cdot d\mathbf{x})+(d{\mathbf{x}}^T\cdot d\mathbf{v})\\ & \\ & =((\mathbf{L}\cdot d\mathbf{x}{)}^T\cdot d\mathbf{x})+(d{\mathbf{x}}^T\cdot (\mathbf{L}\cdot d\mathbf{x}))\\ & \\ & =(d{\mathbf{x}}^T\cdot {\mathbf{L}}^T\cdot d\mathbf{x})+(d{\mathbf{x}}^T\cdot \mathbf{L}\cdot d\mathbf{x})\\ & \\ & =2(d{\mathbf{x}}^T\cdot \mathbf{D}\cdot d\mathbf{x})\end{array}$$

4、Green应变率与变形率
$$\dot{\mathbf{E}}=\frac{1}{2}({\dot{\mathbf{F}}}^T\cdot \mathbf{F}+{\mathbf{F}}^T\cdot \dot{\mathbf{F}})$$
$$\mathbf{D}=\frac{1}{2}({\dot{\mathbf{F}}}^T\cdot {\mathbf{F}}^{-1}+{\mathbf{F}}^{-T}\cdot \dot{\mathbf{F}})$$
因此：
$$\dot{\mathbf{E}}=({\mathbf{F}}^T\cdot \mathbf{D}\cdot \mathbf{F}）或\mathbf{D}=(\mathbf{F}\cdot \mathbf{D}\cdot {\mathbf{F}}^T）$$
由此可知：$\dot{\mathbf{E}}$和$\mathbf{D}$是同一个张量分别在初始构型上和当前构型上的分量。

计算不同应变率度量下外力对微元做功的功率

即，应力和应变率的共轭关系

外力对微元做功的功率表示为：$\dot{W}$

已知名义应力乘以名义应变率在初始构型上积分等于功率：
$$\dot{W}={\mathbf{p}:\dot{\mathbf{ϵ}}}^{T}dV=\mathbf{p}:{\dot{\mathbf{F}}}^{\mathbf{T}}dV=\dot{\mathbf{F}}:{\mathbf{p}}^{\mathbf{T}}dV$$

把应变率张量换成变形率 $D$:
$$\dot{W}={\mathbf{F}}^{-T}\cdot \dot{W}\cdot {\mathbf{F}}^{T}dV={\mathbf{F}}^{-T}\cdot \dot{\mathbf{F}}:{\mathbf{p}}^{\mathbf{T}}\cdot {\mathbf{F}}^{T}dV=\mathbf{F}\cdot \mathbf{p}:{\mathbf{F}}^{-T}\cdot \dot{\mathbf{F}}dV$$
$$\dot{W}=\mathbf{F}\cdot \dot{W}\cdot {\mathbf{F}}^{-1}dV=\mathbf{F}\cdot \mathbf{p}:{\dot{\mathbf{F}}}^{\mathbf{T}}\cdot {\mathbf{F}}^{-1}dV$$
因此：
$$\begin{gathered} \dot{W}=\frac{1}{2}\mathbf{F}\cdot \mathbf{p}:\left({\dot{\mathbf{F}}}^{\mathbf{T}}\cdot {\mathbf{F}}^{-1}+{\mathbf{F}}^{-T}\cdot \dot{\mathbf{F}}\right)dV \\ =\mathbf{F}\cdot \mathbf{p}:\mathbf{D}dV=\mathbf{F}\cdot \mathbf{p}:\mathbf{D}{J}^{-1}dv \end{gathered}$$

定义:$\mathbf{\sigma }=J^{-1}\mathbf{F}\cdot \mathbf{p}$(Cauchy应力，真应力); $\mathbf{\tau }=J\mathbf{\sigma }=\mathbf{F}\cdot \mathbf{p}$(Kirchhoff应力)

注：1、$\dot{W}$ 是标量，计算时可以交换位置；2、$dv=JdV$

同理，应变率张量换位Green应变率$\dot{E}$:
$$\dot{W}=\mathbf{p}\cdot {\mathbf{F}}^{-T}:\dot{\mathbf{E}}dV=\mathbf{S}:\dot{\mathbf{E}}dV$$
其中：$S$ 为第二Piola-Kirchhoff应力

客观应力率

最朴素的本构关系：应力=模量矩阵×应变，所以，应力增量=模量矩阵×应变增量。由于应变中不包含旋转，所以通过模量矩阵×应变增量求出的应力增量也不能包含旋转。

旋转造成的应力增量：

纯旋转 ${\hat{\mathbf{\sigma }}}_{\mathbf{X}}$的分量${\sigma }_{X_{ij}}$ 和 ${\mathbf{\sigma }}_{\mathbf{x}}$ 的分量${\sigma }_{ij}$ 相同

$$\begin{array}{rl}\frac{{\mathbf{\sigma }}_{\mathbf{X}}-{\hat{\mathbf{\sigma }}}_{\mathbf{X}}}{\mathbf{\Delta }\mathbf{t}}& =\frac{{\mathbf{\sigma }}_{\mathbf{X}}-\mathbf{R}\cdot {\mathbf{\sigma }}_{\mathbf{X}}\cdot {\mathbf{R}}^{\mathbf{T}}}{\mathbf{\Delta }\mathbf{t}}\\ & =\frac{\partial (\mathbf{R}\cdot {\mathbf{\sigma }}_{\mathbf{X}}\cdot {\mathbf{R}}^{\mathbf{T}})}{\partial \mathbf{R}}\frac{\partial \mathbf{R}}{\partial \mathbf{t}}\\ & =\dot{\mathbf{R}}\cdot {\mathbf{\sigma }}_{\mathbf{X}}\cdot {\mathbf{R}}^{\mathbf{T}}+\mathbf{R}\cdot {\mathbf{\sigma }}_{\mathbf{X}}\cdot {\dot{\mathbf{R}}}^{\mathbf{T}}\\ & =\dot{\mathbf{R}}\cdot {\mathbf{R}}^{\mathbf{T}}\cdot {\mathbf{\sigma }}_{\mathbf{x}}\cdot \mathbf{R}\cdot {\mathbf{R}}^{\mathbf{T}}+\mathbf{R}\cdot {\mathbf{R}}^{\mathbf{T}}\cdot {\mathbf{\sigma }}_{\mathbf{x}}\cdot \mathbf{R}{\dot{\mathbf{R}}}^{\mathbf{T}}\\ & ={\mathbf{\sigma }}_{\mathbf{x}}\cdot {\mathbf{\Omega }}^{\mathbf{T}}+\mathbf{\Omega }\cdot {\mathbf{\sigma }}_{\mathbf{x}}\end{array}$$
其中:$\mathbf{\Omega }=\dot{\mathbf{R}}\cdot {\mathbf{R}}^{\mathbf{T}}$ 角速度张量

由此推出Green-Naghdi客观应力率：
$${\mathbf{\sigma }}^{\nabla \mathbf{G}}=\frac{D\mathbf{\sigma }}{Dt}-\mathbf{\sigma }\cdot {\mathbf{\Omega }}^{\mathbf{T}}-\mathbf{\Omega }\cdot \mathbf{\sigma }或{\sigma }_{ij}^{\nabla G}=\frac{D{\sigma }_{ij}}{Dt}+{\sigma }_{ik}{\Omega }_{kj}-{\Omega }_{ik}{\sigma }_{kj}$$
其中：$\mathbf{\Omega }$为反对称阵，${\mathbf{\Omega }}^T=-\mathbf{\Omega }$
纯旋转下，速度梯度矩阵=转动矩阵(速度梯度矩阵的反对称部分)=角速度矩阵 $\mathbf{L}=\mathbf{W}=\mathbf{\Omega }$，用转动矩阵代替角速度矩阵得到Jaunmann客观应力率：
$${\mathbf{\sigma }}^{\nabla \mathbf{J}}=\frac{D\mathbf{\sigma }}{Dt}-\mathbf{\sigma }\cdot {\mathbf{W}}^{\mathbf{T}}-\mathbf{W}\cdot \mathbf{\sigma }或{\sigma }_{ij}^{\nabla G}=\frac{D{\sigma }_{ij}}{Dt}+{\sigma }_{ik}{W}_{kj}-W_{jk}{\sigma }_{ki}$$
用速度梯度矩阵代替角速度矩阵，由于速度梯度矩阵中含有变形量需要扣除，Truesdell客观应力率：
$$\begin{gathered} {\mathbf{\sigma }}^{\nabla \mathbf{T}}=\frac{D\mathbf{\sigma }}{Dt}-\mathbf{d}\mathbf{i}\mathbf{v}(\mathbf{v})\mathbf{\sigma }-\mathbf{\sigma }\cdot {\mathbf{L}}^{\mathbf{T}}-\mathbf{L}\cdot \mathbf{\sigma } \\ 或{\sigma }_{ij}^{\nabla T}=\frac{D{\sigma }_{ij}}{Dt}-\frac{\partial v_k}{\partial x_k}{\sigma }_{ij}+{\sigma }_{ik}{L}_{kj}-L_{jk}{\sigma }_{ki} \end{gathered}$$

参考书

1. 《Nonlinear Finite Elements for Continua and Structures》 点击下载
2. 《塑性细观力学》 王自强