应变度量及其一般变换关系

Strain Measures and Their General Transformation Relations

摘要: 应变是固体力学中最重要的基本概念之一，直接影响几何方程和物理方程的建立，本文从基于变形梯度分解的基本运动学关系和变形张量入手，阐述讨论了Seth-Hill广义应变度量函数的本质和含义，推导了相同构型不同应变度量之间以及不同构型应变度量之间的一般变换关系，为应变度量的变换和应用提供了算法依据。

Abstract: Strain, one of the most important fundamental concepts in mechanics, directly affects the estab-lishment of geometric and constitutive equations. This paper reviews the essential kinematic rela-tions and the different stretch tensors within continuum based on the decomposition of the de-formation gradient. According to the definition of strain, the essence of Seth-Hill general strain measures is discussed. Furthermore, the general transformation relations between different strain measures in the same configuration are derived as well as those in different configurations to provide an algorithm basis for the transformation and application of strain measures.

1. 引言

应变是固体力学中最重要的基本概念之一，它直接影响几何方程和物理方程的建立。基于变形梯度可在初始构型和现时构型上定义不同的应变度量 [1] [2] [3]。

目前，在非线性固体力学中，基于不同变形和应变度量的计算是一个传统但仍然非常活跃的研究领域。Xiao [4]、Ogden [5]、Javier [6] 等学者，分别采用对数应变或考虑应变高阶项讨论材料的非线性特征。在本构建模、数值计算中，有时需要对不同构型、不同应变度量进行的等值变换 [7] [8]。Seth [9] 和Hill [10] 在保证变形张量为充分光滑的单调张量函数基础上，利用伸长张量开创性地定义了通类的应变度量函数，为应变张量的设计提供了指导原则。常见的实用应变，比如线性工程应变、对数应变、采用Lagrange描述的Green应变、采用Euler描述的Almansi应变等，都是Seth-Hill广义应变度量的特定退化形式 [11] [12]。在工程问题中，如何选择合适的应变度量描述实际变形，既是理论上易混淆的问题，也是工程实践中经常引起技术人员困惑的问题。本文从连续介质力学变形度量和应变张量的定义出发，讨论了应变度量的类型及其本质，推导了相同构型上和不同构型上应变度量之间的一般关系，为正确理解应变度量提供理论参考。

2. 变形度量

物体 $\mathcal{B}$ 在 $t_0$ 时刻填满区域 ${\mathcal{B}}_0\subset {ℝ}^3$，运动后t时刻填满区域 ${\mathcal{B}}_t\subset {ℝ}^3$，其运动可用如下映射表示：

$�\left(X,t\right):{\mathcal{B}}_0\times ℝ\to {\mathcal{B}}_t\text{or}\left(X,t\right)\mapsto x$ (1)

位移场的Lagrange描述和Euler描述分别为

$u\left(X,t\right)={�}_t\left(X\right)-X$ (2.1)

$u\left(x,t\right)=x-{�}_t^{-1}\left(x\right)$ (2.2)

相应的变形梯度为初始构型切线空间 $T{\mathcal{B}}_0$ 到现时构型切线空间 $T{\mathcal{B}}_t$ 的切线映射：

$F=\frac{\partial �}{\partial X}:T{\mathcal{B}}_0\to T{\mathcal{B}}_t\text{with}\text{det}\left(F\right)>0$ (3)

引入 $T{\mathcal{B}}_0$ 和 $T{\mathcal{B}}_t$ 的对偶空间 $T^{\text{*}}{\mathcal{B}}_0$ 和 $T^{\text{*}}{\mathcal{B}}_t$ 后，存在以下协变和逆变度规张量：

$G:T{\mathcal{B}}_0\to T^{\text{*}}{\mathcal{B}}_0,G^{-1}:T^{\text{*}}{\mathcal{B}}_0\to T{\mathcal{B}}_0$ (4.1)

$g:T{\mathcal{B}}_t\to T^{\text{*}}{\mathcal{B}}_t,g^{-1}:T^{\text{*}}{\mathcal{B}}_t\to T{\mathcal{B}}_t$ (4.2)

考虑完备性，变形梯度和度规张量满足以下逆协变二阶恒等式：

$I=F^{-1}\cdot F=G^{-1}\cdot G:T{\mathcal{B}}_0\to T{\mathcal{B}}_0$ (5.1)

$i=F\cdot F^{-1}=g\cdot g^{-1}:T{\mathcal{B}}_t\to T{\mathcal{B}}_t$ (5.2)

根据极分解定理，有

$F=R\cdot U=v\cdot R$ (6)

其中，旋转张量 $R$ 正交， $\det \left(R\right)=1$ ；右拉伸张量 $U$ 和左拉伸张量 $v$ 均对称正定； $F$ 、 $U$ 、 $v$ 具有相同的特征值 ${\lambda }_i\left(i=1,2,3\right)$。

协变和逆变度规张量存在如下关系：

$R^{\text{T}}\cdot g\cdot R=G$ (7)

应用谱分解定理，式(6)中的各个张量可用特征值 ${\lambda }_i$ 和特征向量 $N_i$ 、 $N^i$ 、 $n_i$ 、 $n^i$ 表示为

$F={\displaystyle \underset{i=1}{\overset{3}{\sum }}{\lambda }_i{n}_i\otimes N^i},R={\displaystyle \underset{i=1}{\overset{3}{\sum }}{n}_i\otimes N^i},U={\displaystyle \underset{i=1}{\overset{3}{\sum }}{\lambda }_i{N}_i\otimes N^i},v={\displaystyle \underset{i=1}{\overset{3}{\sum }}{\lambda }_i{n}_i\otimes n^i}$ (8)

由式(4)的映射可知，特征向量存在关系：

$N_i=G^{-1}\cdot N^i,N^i=G\cdot N_i$ (9.1)

$n_i=g^{-1}\cdot n^i,n^i=g\cdot n_i$ (9.2)

而且满足：

$N^i\cdot G^{-1}\cdot N^i=N_i\cdot G\cdot N_i=1$ (10.1)

$n^i\cdot g^{-1}\cdot n^i=n_i\cdot g\cdot n_i=1$ (10.2)

根据极分解定理，得

$n_i=R\cdot N_i,N_i=R^{\text{T}}\cdot n_i$ (11.1)

$n^i=R\cdot N^i,N^i=R^{\text{T}}\cdot n^i$ (11.2)

在上述关系的基础上，基于线元平方长度的变化，可以导出以下不同的变形张量：

$C=F^{\text{T}}\cdot g\cdot F=U^{\text{T}}\cdot G\cdot U={\displaystyle \underset{i=1}{\overset{3}{\sum }}{\lambda }_i^2{N}^i\otimes N^i}$ (12.1)

$B=F^{-1}\cdot g^{-1}\cdot F^{-\text{T}}=U^{-1}\cdot G^{-1}\cdot U^{-\text{T}}={\displaystyle \underset{i=1}{\overset{3}{\sum }}{\lambda }_i^{-2}{N}_i\otimes N_i}$ (12.2)

$c=F^{-\text{T}}\cdot G\cdot F^{-1}=v^{-\text{T}}\cdot g\cdot v^{-1}={\displaystyle \underset{i=1}{\overset{3}{\sum }}{\lambda }_i^{-2}{n}^i\otimes n^i}$ (12.3)

$b=F\cdot G^{-1}\cdot F^{\text{T}}=v\cdot g^{-1}\cdot v^{\text{T}}={\displaystyle \underset{i=1}{\overset{3}{\sum }}{\lambda }_i^2{n}_i\otimes n_i}$ (12.4)

式(12)定义的变形张量均对称正定，采用主轴表述法，用特征值表征主方向上的主伸长量，而特征向量是与构型对应的变形主轴方向，这种表示法为变形给出了一个很形象的几何解释。显然，

$B\cdot C=I,b\cdot c=i$ (13)

这些变形张量之间的关系如表1。

Table 1. Relations between deformation tensors

表1. 变形张量之间的关系

说明：表中INV表示互逆，EV表示特征值相同，EV⁻¹表示特征值互为倒数。

3. Hill-Seth广义应变度量

与变形度量相对应，应变度量同样也可以表达在不同构型上，相应应变度量的几何解释当然也就不同，由此产生了不同的应变张量。

根据Hill-Seth通类应变度量函数的定义原则，物质应变定义在初始构型 ${\mathcal{B}}_0$ 上，采用Lagrange描述，与右伸长张量 $U$ 共轴。物质应变既可以采用变形张量 $C$ 定义协变的应变张量 $E_{\left(m\right)}\left(U;C;G\right)$，也可以采用变形张量 $B$ 定义逆变的应变张量 $K_{\left(m\right)}\left(U;B;G^{-1}\right)$，应变度量函数如下：

$E_{\left(m\right)}\left(U;C;G\right)=\{\begin{array}{l}{\displaystyle \underset{i=1}{\overset{3}{\sum }}\frac{1}{m}\left({\lambda }_i^m-1\right)N^i\otimes N^i}=\frac{1}{m}\left[C^{\frac{m}{2}}-G\right]=\frac{1}{m}\left(U^m-I\right)\cdot G\text{with}m>0\\ {\displaystyle \underset{i=1}{\overset{3}{\sum }}\ln {\lambda }_i{N}^i\otimes N^i}=\frac{1}{2}\ln C=\ln U\cdot G\text{if}m=0\end{array}$ (14)

$K_{\left(m\right)}\left(U;B;G^{-1}\right)=\{\begin{array}{l}{\displaystyle \underset{i=1}{\overset{3}{\sum }}\frac{1}{m}\left(1-{\lambda }_i^{-m}\right)N_i\otimes N_i}=\frac{1}{m}\left[G^{-1}-B^{\frac{m}{2}}\right]=\frac{1}{m}\left(I^{-1}-U^{-m}\right)\cdot G^{-1}\text{with}m>0\\ {\displaystyle \underset{i=1}{\overset{3}{\sum }}\ln {\lambda }_i{N}_i\otimes N_i}=-\frac{1}{2}\ln B=-\ln \left(U^{-1}\right)\cdot G^{-1}\text{if}m=0\end{array}$ (15)

空间应变定义在现时构型 ${\mathcal{B}}_t$ 上，采用Euler描述，与左伸长张量 $v$ 共轴。空间应变可以分别采用变形张量 $c$ 或 $b$ 定义协变的应变张量 $e_{\left(m\right)}\left(v;c;g\right)$ 和逆变的应变张量 $k_{\left(m\right)}\left(v;b;g^{-1}\right)$，应变度量函数如下：

$e_{\left(m\right)}\left(v;c;g\right)=\{\begin{array}{l}{\displaystyle \underset{i=1}{\overset{3}{\sum }}\frac{1}{m}\left(1-{\lambda }_i^{-m}\right)n^i\otimes n^i}=\frac{1}{m}\left[g-c^{\frac{m}{2}}\right]=\frac{1}{m}\left(i-v^{-m}\right)\cdot g\text{with}m>0\\ {\displaystyle \underset{i=1}{\overset{3}{\sum }}\ln {\lambda }_i{n}^i\otimes n^i}=-\frac{1}{2}\ln c=-\ln \left(v^{-1}\right)\cdot g\text{if}m=0\end{array}$ (16)

$k_{\left(m\right)}\left(v;b;g^{-1}\right)=\{\begin{array}{l}{\displaystyle \underset{i=1}{\overset{3}{\sum }}\frac{1}{m}\left({\lambda }_i^m-1\right)n_i\otimes n_i}=\frac{1}{m}\left[b^{\frac{m}{2}}-g^{-1}\right]=\frac{1}{m}\left(v^m-i^{-1}\right)\cdot g^{-1}\text{with}m>0\\ {\displaystyle \underset{i=1}{\overset{3}{\sum }}\ln {\lambda }_i{n}_i\otimes n_i}=\frac{1}{2}\ln b=\ln v\cdot g^{-1}\text{if}m=0\end{array}$ (17)

式(14)~(17)中的m是非负整数。

Hill-Seth应变度量函数本质上表示了线元的幂次(广义几何长度)改变相对于不同参考构架的变化率。具体来讲， $E_{\left(m\right)}\left(U;C;G\right)$ 表示线元的幂次改变相对于初始构型的协变度规 $G$ 的变化率； $e_{\left(m\right)}\left(v;c;g\right)$ 表示线元的幂次改变相对于现时构型的协变度规 $g$ 的变化率；而 $K_{\left(m\right)}\left(U;B;G^{-1}\right)$ 则表示线元幂次变化相对于初始构型的逆变度规 $G^{-1}$ 的变化率；同理， $k_{\left(m\right)}\left(v;b;g^{-1}\right)$ 表示线元幂次变化相对于现时构型的逆变度规 $g^{-1}$ 的变化率。Hill-Seth通类应变度量函数是单调函数，显然， ${\lambda }_i=1$ 时对应主轴方向的主应变为零。

4. 相同构型协变与逆变应变度量之间的关系

前述Hill-Seth应变度量函数族中，定义在同一构型上相同m值的应变族存在协变张量与逆变张量之间的转换关系。

因为 $B$ 与 $C$ 互逆，可得

$\begin{array}{c}{C}^{\frac{m}{2}}\cdot K_{\left(m\right)}\cdot G=\frac{1}{m}{C}^{\frac{m}{2}}\cdot \left[G^{-1}-B^{\frac{m}{2}}\right]\cdot G\\ =\frac{1}{m}\left[C^{\frac{m}{2}}\cdot G^{-1}\cdot G-C^{\frac{m}{2}}\cdot B^{\frac{m}{2}}\cdot G\right]\\ =\frac{1}{m}\left[C^{\frac{m}{2}}\cdot I-I\cdot G\right]=\frac{1}{m}\left[C^{\frac{m}{2}}-G\right]\end{array}$

即，

$E_{\left(m\right)}=C^{\frac{m}{2}}\cdot K_{\left(m\right)}\cdot G$ (18.1)

类似地，因为 $b$ 与 $c$ 互逆，可得

$\begin{array}{c}{c}^{\frac{m}{2}}\cdot k_{\left(m\right)}\cdot g=\frac{1}{m}{c}^{\frac{m}{2}}\cdot \left[b^{\frac{m}{2}}-g^{-1}\right]\cdot g\\ =\frac{1}{m}\left[c^{\frac{m}{2}}\cdot b^{\frac{m}{2}}\cdot g-c^{\frac{m}{2}}\cdot g^{-1}\cdot g\right]\\ =\frac{1}{m}\left[i\cdot g-c^{\frac{m}{2}}\cdot i\right]=\frac{1}{m}\left[g-c^{\frac{m}{2}}\right]\end{array}$

即，

$e_{\left(m\right)}=c^{\frac{m}{2}}\cdot k_{\left(m\right)}\cdot g$ (18.2)

关系式(18.1)和(18.2)也可以写成以下等价形式：

$K_{\left(m\right)}=B^{\frac{m}{2}}\cdot E_{\left(m\right)}\cdot G^{-1}$ (19.1)

$k_{\left(m\right)}=b^{\frac{m}{2}}\cdot e_{\left(m\right)}\cdot g^{-1}$ (19.2)

将式(12.1)、(12.3)分别代入式(18.1)、(18.2)， $E_{\left(m\right)}$ 与 $K_{\left(m\right)}$ 、 $e_{\left(m\right)}$ 与 $k_{\left(m\right)}$ 的变换关系可以分别用右拉伸张量 $U$ 和左拉伸张量 $v$ 表示为

$\begin{array}{c}{E}_{\left(m\right)}={\left(U^{\frac{m}{2}}\right)}^{\text{T}}\cdot G\cdot U^{\frac{m}{2}}\cdot K_{\left(m\right)}\cdot G\\ ={\left(U^{\frac{m}{2}}\right)}^{\text{T}}\cdot G\cdot K_{\left(m\right)}\cdot {\left(U^{\frac{m}{2}}\right)}^{\text{T}}\cdot G\\ ={\left(U^{\frac{m}{2}}\right)}^{\text{T}}\cdot G\cdot K_{\left(m\right)}\cdot G\cdot U^{\frac{m}{2}}\end{array}$ (20.1)

$\begin{array}{c}{e}_{\left(m\right)}={\left(v^{-\frac{m}{2}}\right)}^{\text{T}}\cdot g\cdot v^{-\frac{m}{2}}\cdot k_{\left(m\right)}\cdot g\\ ={\left(v^{-\frac{m}{2}}\right)}^{\text{T}}\cdot g\cdot k_{\left(m\right)}\cdot {\left(v^{-\frac{m}{2}}\right)}^{\text{T}}\cdot g\\ ={\left(v^{-\frac{m}{2}}\right)}^{\text{T}}\cdot g\cdot k_{\left(m\right)}\cdot g\cdot v^{-\frac{m}{2}}\end{array}$ (20.2)

同理，式(19.1)和(19.2)也可以用拉伸张量表示为

$K_{\left(m\right)}=U^{-\frac{m}{2}}\cdot G^{-1}\cdot E_{\left(m\right)}\cdot G^{-1}\cdot {\left(U^{-\frac{m}{2}}\right)}^{\text{T}}$ (21.1)

$k_{\left(m\right)}=v^{\frac{m}{2}}\cdot g^{-1}\cdot e_{\left(m\right)}\cdot g^{-1}\cdot {\left(v^{\frac{m}{2}}\right)}^{\text{T}}$ (21.2)

5. 不同构型协变应变张量之间或逆变应变张量之间的关系

初始构型和现时构型上应变张量之间的变换称为“拉回”和“推前”操作。协变张量的“拉回”和“推前”分别记为 ${\varphi }^{\leftarrow }(\cdot )$ 和 ${\varphi }^{\to }(\cdot )$，逆变张量的“拉回”和“推前”分别记为 ${\chi }^{\leftarrow }(\cdot )$ 和 ${\chi }^{\to }(\cdot )$，具体操作运算如下：

${\varphi }_{\left(F\right)}^{\leftarrow }(\cdot )=F^{\text{T}}\cdot (\cdot )\cdot F,{\varphi }_{\left(F\right)}^{\to }(\cdot )=F^{-\text{T}}\cdot (\cdot )\cdot F^{-1}$ (22)

${\chi }_{\left(F\right)}^{\leftarrow }(\cdot )=F^{-1}\cdot (\cdot )\cdot F^{-\text{T}},{\chi }_{\left(F\right)}^{\to }(\cdot )=F\cdot (\cdot )\cdot F^{\text{T}}$ (23)

由此可得，同为协变张量的 $E_{\left(m\right)}$ 与 $e_{\left(m\right)}$ “拉回”和“推前”关系为

$E_{\left(m\right)}={\varphi }_{\left(F\right)}^{\leftarrow }\left(e_{\left(m\right)}\right)=F^{\text{T}}\cdot e_{\left(m\right)}\cdot F$ (24.1)

$e_{\left(m\right)}={\varphi }_{\left(F\right)}^{\to }\left(E_{\left(m\right)}\right)=F^{-\text{T}}\cdot E_{\left(m\right)}\cdot F^{-1}$ (24.2)

同为逆变张量的 $K_{\left(m\right)}$ 与 $k_{\left(m\right)}$ “拉回”和“推前”关系为

$K_{\left(m\right)}={\chi }_{\left(F\right)}^{\leftarrow }\left(k_{\left(m\right)}\right)=F^{-1}\cdot k_{\left(m\right)}\cdot F^{-\text{T}}$ (25.1)

$k_{\left(m\right)}={\chi }_{\left(F\right)}^{\to }\left(K_{\left(m\right)}\right)=F\cdot K_{\left(m\right)}\cdot F^{\text{T}}$ (25.2)

6. 不同构型协变与逆变应变张量之间的关系

比较式(14)~(17)中用主轴法表示的谱分解应变度量函数， $E_{\left(m\right)}$ 与 $k_{\left(m\right)}$ 、 $K_{\left(m\right)}$ 与 $e_{\left(m\right)}$ 具有完全相同的特征值，只是主轴方向不同，前后二者之间的“拉回”和“推前”变换只涉及旋转张量 $R$。

将(7)和(9.1)代入(14)用主轴表示的应变度量函数中，有

$E_{\left(m\right)}=\underset{i=1}{\overset{3}{{\displaystyle \sum }}}\frac{1}{m}\left({\lambda }_i^m-1\right)G\cdot N_i\otimes N_i\cdot G=\underset{i=1}{\overset{3}{{\displaystyle \sum }}}\frac{1}{m}\left({\lambda }_i^m-1\right)R^{\text{T}}\cdot g\cdot R\cdot N_i\otimes N_i\cdot R^{\text{T}}\cdot g\cdot R$

$E_{\left(m\right)}$ 做 $R$ 旋转推前操作，

$\begin{array}{c}{R}^{-\text{T}}\cdot E_{\left(m\right)}\cdot R^{-1}=\underset{i=1}{\overset{3}{{\displaystyle \sum }}}\frac{1}{m}\left({\lambda }_i^m-1\right)R^{-\text{T}}\cdot R^{\text{T}}\cdot g\cdot R\cdot N_i\otimes N_i\cdot R^{\text{T}}\cdot g\cdot R\cdot R^{-1}\\ =\underset{i=1}{\overset{3}{{\displaystyle \sum }}}\frac{1}{m}\left({\lambda }_i^m-1\right)g\cdot R\cdot N_i\otimes N_i\cdot R^{\text{T}}\cdot g\\ =\underset{i=1}{\overset{3}{{\displaystyle \sum }}}\frac{1}{m}\left({\lambda }_i^m-1\right)g\cdot n_i\otimes n_i\cdot g=g\cdot k_{\left(m\right)}\cdot g\end{array}$

从而得到 $E_{\left(m\right)}$ 与 $k_{\left(m\right)}$ 的变换关系为

(26.1)

类似地，对做旋转推前，可得

(26.2)

对和做旋转拉回，可得

(26.3)

(26.4)

式(26.1)~(26.4)所表示的不同构型协变与逆变应变张量的转换关系也可以通过式(22)和(23)所表示的完整“拉回”和“推前”操作，用左右伸长张量写成如下形式：

(27.1)

(27.2)

(27.3)

(27.4)

7. 常用参考架构下的简化

应变度量通常引入参考架构用应变张量的分量组合表示，应变张量与参考架构无关，但应变张量分量与参考架构有关。引入度规张量以后，式(14)~(17)的应变度量对初始构型和现时构型采用何种参考架构没有约束和限制，采用不同的坐标系，实际上就是选择了不同的基矢和度规，由此推导的上述变换关系不受参考构架的影响，既可以采用直线坐标系，也可以选择曲线坐标系，坐标轴可以正交，也可以斜交。

根据度规张量的几何意义，其分量等于两个对应坐标轴的切向基矢的点积，它反映了非正交曲线坐标系的坐标轴夹角和弯曲程度。正交坐标系度规张量的副元分量都为零，斜角坐标系则不全为零；工程实用中常用的参考架构主要是正交坐标系，比如直角坐标系、球坐标系、柱坐标系，此时，度规张量除了主元分量非零以外，其他分量都为零。球坐标系和柱坐标系是曲线坐标系，度规张量的主元分量不全为1，大小与坐标轴的曲率有关，协变度规张量与逆变度规张量对应的主元分量呈倒数关系。

对于笛卡尔直角坐标系，每个坐标轴上的量度一致，对应的构型空间为欧氏空间，均匀且各向同性，协变基矢和逆变基矢完全一致，都是正交单位矢量，度规张量为单位张量，即

此时，考虑到变形张量和特征值相同，虽然式(14)和式(17)的应变度量仍具有不同的物理意义和几何意义，但二者数值是相等的；同理，式(15)和式(16)的应变度量数值也相同。

所以对于笛卡尔应变张量，协变张量和分别与逆变张量和等价，

此时无需再把应变张量分成协变和逆变，通常仅考虑基于初始构型的应变度量和基于现时构型的应变度量。

8. 结语

本文基于Seth-Hill广义应变度量函数，推导了由此衍生的同一构型上的应变度量，即，与、与的协变与逆变变换关系，以及不同构型上应变度量，即，与、与以及与、与的“推前”与 “拉回”变换关系，这些变换关系不受参考构架的约束和限制，有利于拓宽材料本构建模和非线性计算的手段，为复杂工程问题的数值计算提供算法依据。

非线性连续介质力学的本构理论中，应变度量和应力度量存在能量共轭关系，二者成对用于建立能量泛函。本文仅从理论上探索了不同应变度量之间的变换，各类应变度量共轭应力与之配套的相应变换有待于进一步研究。

基金项目

国家自然科学基金项目(51578405)资助。