Python温度应力 温度应力计算公式

文章目录

• 1.热应力概述
• 2.广义胡克定律
• 3.位移平衡方程

1.热应力概述

物体的变形不仅仅由外力作用引起，温度的变化也能够引起变形，称为热变形：如热胀冷缩。

但是单有温度的变化，不一定就在物体内部产生应力，只有当温度变化所引起的膨胀或收缩受到约束时，才会在物体内产生应力。这种无外力作用而是由于温度变化引起的热变形受到约束而产生的应力，称为热应力或温度应力。

• 外部变形约束：物体受到外部条件影响，不能自由伸缩，如两端固定等；
• 内部变形约束：物体内部温度分布不均匀，如物体内外、上下存在温差，导致零件翘曲；
• 相互变形约束：不同零件组合在一起，热膨胀系数不同，相互影响；

2.广义胡克定律

对于各向同性物体，在自由膨胀情况下，假设温度变化为 t，热膨胀率为 $\alpha$，则：
$\begin{aligned} &\varepsilon_x=\varepsilon_y=\varepsilon_z=\alpha t\\ &\gamma_{xy}=\gamma_{yz}=\gamma_{zx}=0 \end{aligned}\tag1$

应变一般由两部分组成：一部分为由温度变化引起，另一部分为应力变化引起。现将热应力引入弹性力学中的胡克定律中( 这一部分内容可参考应力应变基础理论分析 )，可得：
$\begin{cases}\varepsilon_x=\frac{1}{E}[\sigma_x-\mu(\sigma_y+\sigma_z)]+\alpha t\\[1.5ex] \varepsilon_y=\frac{1}{E}[\sigma_y-\mu(\sigma_x+\sigma_z)]+\alpha t\\[1.5ex] \varepsilon_z=\frac{1}{E}[\sigma_z-\mu(\sigma_x+\sigma_y)]+\alpha t\\[1.5ex] \gamma_{xy}=\frac{\tau_{xy}}{G}\\[1.5ex] \gamma_{yz}=\frac{\tau_{yz}}{G}\\[1.5ex] \gamma_{zx}=\frac{\tau_{zx}}{G} \end{cases}\tag2$

已知体积应力 $\Theta =\sigma_x+\sigma_y+\sigma_z\tag3$，
剪切弹性模量$G=\frac{E}{2(1+\mu)}\tag4$
式（2）中正应变可整理为：

$\begin{cases} \varepsilon_x=\frac{1}{2G}(\sigma_x-\frac\mu{1+\mu}\Theta)+\alpha t\\[1.5ex] \varepsilon_y=\frac{1}{2G}(\sigma_y-\frac\mu{1+\mu}\Theta)+\alpha t\\[1.5ex] \varepsilon_z=\frac{1}{2G}(\sigma_z-\frac\mu{1+\mu}\Theta)+\alpha t \end{cases}\tag5$

将应力转变为应变的函数：
$\begin{cases}\sigma_x=2G\varepsilon_x+\frac{\mu}{1+\mu}\Theta-2G\alpha t\\[1.5ex] \sigma_y=2G\varepsilon_y+\frac{\mu}{1+\mu}\Theta-2G\alpha t\\[1.5ex] \sigma_z=2G\varepsilon_y+\frac{\mu}{1+\mu}\Theta-2G\alpha t\\[1.5ex] \tau_{xy}=G\gamma_{xy}=\frac E{2(1+\mu)}\gamma_{xy}\\[1.5ex] \tau_{yz}=G\gamma_{yz}=\frac E{2(1+\mu)}\gamma_{yz}\\[1.5ex] \tau_{zx}=G\gamma_{zx}=\frac E{2(1+\mu)}\gamma_{zx} \end{cases}\tag6$

由式（5)，对于体积应变 $e$，
$e=\varepsilon_x+\varepsilon_y+\varepsilon_z=\frac{1-2\mu}{1+\mu}\frac{\Theta}{2G}+3\alpha t\tag7$

即
$\Theta=\frac{E}{1-2\mu}(e-3\alpha t)\tag8$

根据式（6)、式 (7),将正应力公式整理为以下形式：

$\begin{cases}\sigma_x=2G\varepsilon_x+\lambda e -\beta t\\[1.5ex] \sigma_y=2G\varepsilon_y+\lambda e -\beta t\\[1.5ex] \sigma_z=2G\varepsilon_y+\lambda e -\beta t \end{cases}\tag9$

其中，$\lambda$

$\lambda=\frac{E\mu}{(1+\mu)(1-2\mu)}\tag{10}$

$\beta$ 为热应力系数：
$\beta=\alpha(3\lambda+2G)\tag {11}$

3.位移平衡方程

我们假设热应力为 $p$，则：
$\nabla p = \left[\begin{matrix} \frac{\partial\sigma_{x}}{\partial x}+\frac{\partial\tau_{yx}}{\partial y}+\frac{\partial\tau_{zx}}{\partial z} \\[1.5ex] \frac{\partial\tau_{xy}}{\partial x}+\frac{\sigma\tau_{y}}{\partial y}+\frac{\partial\tau_{zy}}{\partial z} \\[1.5ex] \frac{\partial\tau_{xz}}{\partial x}+\frac{\partial\tau_{yz}}{\partial y}+\frac{\partial\sigma_{z}}{\partial z} \end{matrix}\right] \tag{12}$

仅受应力影响的应变微分方程位移形式表述如下：

$\varepsilon=\begin{bmatrix}\varepsilon_x \\\varepsilon_y \\\varepsilon_z\\ \gamma_{xy}\\\gamma_{yz}\\\gamma_{zx} \end{bmatrix}=\begin{bmatrix} \frac{\partial u}{\partial x} \\\frac{\partial v}{\partial y}\\ \frac{\partial w}{\partial z} \\ \frac{\partial v}{\partial x}+\frac{\partial u}{\partial y} \\\frac{\partial w}{\partial y}+\frac{\partial v}{\partial z} \\\frac{\partial u}{\partial z}+\frac{\partial w}{\partial x} \end{bmatrix} \tag{13}$

结合式（6）、（9）、（12）、（13）可得：

$\nabla p = \left[\begin{matrix} (\lambda+G)\frac{\partial e}{\partial x}+G\nabla^2u-\beta\frac{\partial t}{\partial x} \\[1.5ex] (\lambda+G)\frac{\partial e}{\partial y}+G\nabla^2v-\beta\frac{\partial t}{\partial y} \\[1.5ex] (\lambda+G)\frac{\partial e}{\partial z}+G\nabla^2w-\beta\frac{\partial t}{\partial z} \end{matrix}\right] \tag{14}$

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

$\begin{cases} \sigma_{x}l+\tau_{yx}m+\tau_{zx}n=0\\[1.5ex] \sigma_{y}m+\tau_{yx}n+\tau_{zx}l =0\\[1.5ex] \sigma_{z}n+\tau_{yx}l+\tau_{zx}m=0 \end{cases} \tag{15}$

其中，l、m、n 分别为边界表面在x、y、z方向的法向余弦。