杜哈梅尔相似定理

杜哈梅尔（Duhamel）相似定理用于求解不均匀受热引起的热应力问题，通过将其转化为等温弹性力学问题，进而求解。
根据热应力基础概念，正应力方程可整理成以下形式：
$$\{\begin{array}{l}{\sigma }_x=2G{\epsilon }_x+\lambda e-\beta t\\ {\sigma }_y=2G{\epsilon }_y+\lambda e-\beta t\\ {\sigma }_z=2G{\epsilon }_y+\lambda e-\beta t\end{array}\text{\hspace{1em}(1)}$$

式（1）可以分为以下两部分：
$$\{\begin{array}{l}{\sigma }_x={\sigma }_x^{{}^{\prime }}-\beta t\\ {\sigma }_y={\sigma }_y^{{}^{\prime }}-\beta t\\ {\sigma }_z={\sigma }_z^{{}^{\prime }}-\beta t\end{array}\text{\hspace{1em}(2)}$$

其中，${\sigma }_x^{{}^{\prime }}$、${\sigma }_y^{{}^{\prime }}$、${\sigma }_z^{{}^{\prime }}$ 是应变产生的正应力，可整理为：
$$\{\begin{array}{l}{\sigma }_x^{{}^{\prime }}=2G{\epsilon }_x+\lambda e\\ {\sigma }_y^{{}^{\prime }}=2G{\epsilon }_y+\lambda e\\ {\sigma }_z^{{}^{\prime }}=2G{\epsilon }_y+\lambda e\end{array}\text{\hspace{1em}(3)}$$

剪应力则保持不变：

$$\{\begin{array}{l}{\tau }_{xy}={\tau }_{xy}^{{}^{\prime }}=G{\gamma }_{xy}=\frac{E}{2(1+\mu )}{\gamma }_{xy}\\ {\tau }_{yz}={\tau }_{xy}^{{}^{\prime }}=G{\gamma }_{yz}=\frac{E}{2(1+\mu )}{\gamma }_{yz}\\ {\tau }_{zx}={\tau }_{xy}^{{}^{\prime }}=G{\gamma }_{zx}=\frac{E}{2(1+\mu )}{\gamma }_{zx}\end{array}\text{\hspace{1em}(4)}$$

当不计体积力时，即：
$$\nabla \mathbf{p}=\left[\begin{array}{c}\frac{\partial {\sigma }_x}{\partial x}+\frac{\partial {\tau }_{yx}}{\partial y}+\frac{\partial {\tau }_{zx}}{\partial z}\\ \frac{\partial {\tau }_{xy}}{\partial x}+\frac{\sigma {\tau }_y}{\partial y}+\frac{\partial {\tau }_{zy}}{\partial z}\\ \frac{\partial {\tau }_{xz}}{\partial x}+\frac{\partial {\tau }_{yz}}{\partial y}+\frac{\partial {\sigma }_z}{\partial z}\end{array}\right]=0\text{\hspace{1em}(5)}$$

将式（3）、式（4）代入式(5) 中，可得：
$$\nabla \mathbf{p}=\left[\begin{array}{c}\frac{\partial {\sigma }_x^{{}^{\prime }}}{\partial x}+\frac{\partial {\tau }_{yx}^{{}^{\prime }}}{\partial y}+\frac{\partial {\tau }_{zx}^{{}^{\prime }}}{\partial z}-\beta \frac{\partial t}{\partial x}\\ \frac{\partial {\tau }_{xy}^{{}^{\prime }}}{\partial x}+\frac{\sigma {\tau }_y^{{}^{\prime }}}{\partial y}+\frac{\partial {\tau }_{zy}^{{}^{\prime }}}{\partial z}-\beta \frac{\partial t}{\partial x}\\ \frac{\partial {\tau }_{xz}^{{}^{\prime }}}{\partial x}+\frac{\partial {\tau }_{yz}^{{}^{\prime }}}{\partial y}+\frac{\partial {\sigma }_z^{{}^{\prime }}}{\partial z}-\beta \frac{\partial t}{\partial x}\end{array}\right]=0\text{\hspace{1em}(6)}$$

除了必须满足微分平衡方程外，也需要满足边界条件。当忽略表面力（即表面力为零）时：

$$\{\begin{array}{l}\beta tl={\sigma }_x^{{}^{\prime }}l+{\tau }_{yx}^{{}^{\prime }}m+{\tau }_{zx}^{{}^{\prime }}n\\ \beta tm={\sigma }_y^{{}^{\prime }}m+{\tau }_{yx}^{{}^{\prime }}n+{\tau }_{zx}^{{}^{\prime }}l\\ \beta tn={\sigma }_z^{{}^{\prime }}n+{\tau }_{yx}^{{}^{\prime }}l+{\tau }_{zx}^{{}^{\prime }}m\end{array}\text{\hspace{1em}(7)}$$

其中，l、m、n 分别为边界表面在x、y、z方向的法向余弦。通过上式，我们不难看出：

• 在求解热应力问题时，我们可以将其等效为等温弹性问题；
• 假定体积力为：（$-\beta \frac{\partial t}{\partial x}$、$-\beta \frac{\partial t}{\partial x}$、$-\beta \frac{\partial t}{\partial x}$）
• 假定表面力为：（$\beta tl$、$\beta tm$、$\beta tn$）
• 求解该等温弹性问题下的应力${\sigma }_x^{{}^{\prime }}$、${\sigma }_y^{{}^{\prime }}$、${\sigma }_z^{{}^{\prime }}、{\tau }_{xy}^{{}^{\prime }}、{\tau }_{yz}^{{}^{\prime }}、{\tau }_{zx}^{{}^{\prime }}$;
• 对正应力分量${\sigma }_x^{{}^{\prime }}$、${\sigma }_y^{{}^{\prime }}$、${\sigma }_z^{{}^{\prime }}$叠加上该点的流体静压 $-\beta t$；

通过以上步骤即可得出物体不均匀受热引起的热应力。