连续介质力学中含张量的应变能积分求解方法咨询

别担心，这个张量积分的逻辑其实和你熟悉的标量积分一脉相承，核心是抓住爱因斯坦求和约定和线性弹性中刚度张量$C_{pqkl}$为常数的特点，一步步拆解就很清晰了：

步骤1：明确求和规则与积分本质 #

首先回忆爱因斯坦求和约定：公式中重复的指标（比如$pq$、$kl$）代表对所有可能的指标值求和（三维问题里$p,q,k,l$取1,2,3）。这里的积分是对整个应变张量分量的多重积分，从应变全为0的初始状态，积分到最终的$\varepsilon_{ij}$状态。

步骤2：代入线性弹性本构关系 #

线性弹性中应力与应变满足$\sigma_{pq} = C_{pqkl}\varepsilon_{kl}$，把它代入应变能积分：
$$\int_0^{\varepsilon_{ij}} \sigma_{pq}d\varepsilon_{pq} = \int_0^{\varepsilon_{ij}} C_{pqkl}\varepsilon_{kl}d\varepsilon_{pq}$$
因为$C_{pqkl}$是常数四阶张量，不随应变变化，所以可以直接把它从积分里提出来：
$$= C_{pqkl} \int_0^{\varepsilon_{ij}} \varepsilon_{kl}d\varepsilon_{pq}$$

步骤3：拆解积分（线性加载的关键） #

小变形下，应变张量的各个分量是独立变量，但这里要注意：我们是从0线性加载到最终应变状态的——加载过程中任意时刻的应变$\varepsilon_{pq}(t) = t\varepsilon_{pq}^{final}$（$t\in[0,1]$），对应的微分$d\varepsilon_{pq} = \varepsilon_{pq}^{final}dt$，同时$\varepsilon_{kl}(t) = t\varepsilon_{kl}^{final}$。

把这个加载过程代入积分项：
$$\int_0^{{\varepsilon_{pq}}{final}} \varepsilon_{kl}(t)d\varepsilon_{pq}(t) = \int_0^1 t\varepsilon_{kl}^{final} \cdot \varepsilon_{pq}^{final}dt = \varepsilon_{kl}^{{final}\varepsilon_{pq}}{final} \int_0^1 tdt = \frac{1}{2}\varepsilon_{kl}\varepsilon_{pq}$$
（这里省略了final上标，默认指最终应变值）

不管指标$(p,q)$和$(k,l)$是否相等，这个积分结果都是$\frac{1}{2}\varepsilon_{kl}\varepsilon_{pq}$，因为线性加载下所有应变分量同步随t线性变化，不存在“独立常数”的情况。

步骤4：合并结果并替换哑标 #

把积分结果代回原式：
$$C_{pqkl} \cdot \frac{1}{2}\varepsilon_{kl}\varepsilon_{pq} = \frac{1}{2}C_{pqkl}\varepsilon_{pq}\varepsilon_{kl}$$
由于爱因斯坦求和中的哑标可以任意替换（把$pq$换成$ij$不影响求和结果），最终就得到了你看到的形式：
$$\frac{1}{2}C_{ijkl}\varepsilon_{kl}\varepsilon_{ij}$$

额外补充：利用张量对称性简化 #

实际工程中，刚度张量$C_{pqkl}$通常自带对称性（比如$C_{pqkl}=C_{qpkl}$对应应力张量对称性，$C_{pqkl}=C_{pqlk}$对应应变张量对称性），这个表达式还可以进一步简化，但核心的积分推导逻辑就是上面的过程。

内容的提问来源于stack exchange，提问作者Élio Pereira