格林应变和柯西应变

Cauchy Strain

在小变形理论中,柯西应变定义为: \[ \varepsilon =\frac{1}{2} \left( U^T+ U\right)=\frac{1}{2} (u_{i,j}+u_{j,i}) \] 这是在学弹性力学时大家所熟知的应变写法。

弹性力学只关心与应力应变有关的形变(体积改变),而不考虑刚体转动。

小应变理论中的应变定义在连续介质力学中可以用变形梯度张量 \(\bf F\) 来表示

\[ \varepsilon_{ij} = \frac{1}{2} (u_{i,j} + u_{j,i})=\frac{1}{2} (F_{ij} + F_{ji}) - \delta_{ij} \]

应变在面对转动时会遇到问题,因为

\[ \pmb \varepsilon = \frac{1}{2}(\pmb F+\pmb F^T)-\pmb I= \frac{1}{2}(\pmb R \cdot \pmb U+\pmb U^T \cdot \pmb R^T)-\pmb I \]

  • 不存在转动时,即 \(\pmb R=\pmb I\) 时,有 \(\pmb{\varepsilon} = {\mathbf U} - {\mathbf I}\) ,此时应变 \(\pmb \varepsilon\) 仅仅与形变有关。
  • 存在刚体转动时\(\pmb R \ne \pmb I\),上式中的 \(\pmb R\) 就会对应变产生影响,即刚体转动产生应变,这显然不对。

实际上小应变理论里独用了一个张量来处理刚体转动问题

\[ {\begin{bmatrix} du \\ dv \\ dw \end{bmatrix}} \, = \, {\begin{bmatrix} 0 & -\omega_z & \omega_y \\ \omega_x & 0 & -\omega_x \\ -\omega_y & \omega_x & 0 \end{bmatrix}} {\begin{bmatrix} dx \\ dy \\ dz \end{bmatrix}} \, + \, {\begin{bmatrix} \varepsilon _x & \frac{1}{2} \gamma _{xy} & \frac{1}{2}\gamma _{xz}\\ \frac{1}{2}\gamma _{yx} & \varepsilon _y & \frac{1}{2}\gamma _{yz}\\ \frac{1}{2}\gamma _{zx} & \frac{1}{2} \gamma _{zy} & \varepsilon _z \end{bmatrix}} {\begin{bmatrix} dx \\ dy \\ dz \end{bmatrix}} \]

Green Strain

而格林应变不受转动矩阵的影响。

由位移公式 \(u=x-X\),有

\[ {\bf F} = \frac{\partial}{\partial \bf X} \left({\bf X} + {\bf u}\right) = { \partial {\bf X} \over \partial {\bf X} } + { \partial {\bf u} \over \partial {\bf X} } = {\bf I} + \frac{\partial{\bf u}}{\partial{\bf X}} \]

这里的 \(\bf U\)(体积改变、旋转)并不是 \(\bf R \cdot \bf U\)\(\bf U\)(体积改变)

有限应变理论中,格林应变

\[ {\bf E} = {1 \over 2} \left( {\bf F}^T \cdot {\bf F} - {\bf I} \right) \]

其中,

\[ {\bf F}^T \cdot {\bf F}=\left( {\bf R} \cdot {\bf U} \right)^T \cdot \left( {\bf R} \cdot {\bf U} \right) = {\bf U}^T \cdot {\bf U} \]

因此格林应变只与形变有关。其中 \(\bf R\) 就是坐标变换矩阵

\[ \bf R^T \cdot \bf R = \bf I \]

\[ {\bf E} = \dfrac{1}{2} \left({\bf F}^T \cdot {\bf F} - {\bf I} \right) \]

写成分量形式,有

\[ E_{ij} = {1 \over 2} \left( {\partial \, u_i \over \partial X_j} + {\partial \, u_j \over \partial X_i} + {\partial \, u_k \over \partial X_i} {\partial \, u_k \over \partial X_j} \right) \]

和柯西应变 \(\varepsilon_{i,j}\) 比较可以发现,差别其实仅仅在于格林应变多了一个第三项,该项是二次项。该项保证了格林应变不受转动的影响。(不过在大变形下就会导致格林应变偏离工程应变 \(e = \dfrac{\Delta L}{L}\))实际上,把格林应变放到小应变假设下,就可以省去二阶以上的高阶项,就有 \(E \approx \varepsilon\)

所以,格林应变的优势就是在于将转动与形变分离开来,只关注形变。