蠕变-松弛耦合
蠕变-松弛耦合描述弹簧在恒定位移下,弹性应变逐步转化为蠕变应变的应力松弛过程。σ(t) 随时间单调下降,下降速率随温度升高而急剧增大。弹簧在 500C 以上的松弛可在数小时内完成。
公式表达式
参数列表
| 符号 | 名称 | 单位 |
|---|---|---|
| stress_initial_MPa | 初始应力 | MPa |
| temp_C | 温度 | °C |
| time_hours | 时间 | h |
详细计算指南
SSHT 蠕变-松弛耦合:恒定位移下的应力松弛
1. 物理机制
碟形弹簧在高温下通常工作于恒定位移(恒定压缩量)状态。此时,弹性应变 $\varepsilon_e$ 逐渐转化为蠕变应变 $\varepsilon_c$,导致应力随时间不断下降,这一过程即为应力松弛。
应力松弛的本质是蠕变驱动的应力重分布:总应变恒定,材料在持续应力作用下发生蠕变,弹性变形量减少,应力随之降低。
2. 应力松弛控制方程
由总应变恒定条件:
将弹性应变 $\varepsilon_e = \sigma / E$ 和 Norton-Bailey 稳态蠕变速率 $\dot{\varepsilon}_c = A \sigma^n \exp(-Q/RT)$ 代入,得到应力松弛的微分方程:
该方程表明,应力下降速率与当前应力的 $n$ 次方成正比。由于金属材料的 $n$ 值通常在 $3 \sim 10$ 之间,初始高应力阶段松弛极快,随应力降低,松弛速率迅速减缓。
3. 应力松弛解析解($n \neq 1$)
对上述微分方程积分,设 $t=0$ 时初始应力为 $\sigma_0$,得任意时刻的剩余应力 $\sigma(t)$:
亦可写为:
参数说明: - $\sigma(t)$:t 时刻的应力 (MPa) - $\sigma_0$:初始应力 (MPa) - $E$:工作温度下的弹性模量 (MPa) - $A, n, Q$:Norton-Bailey 蠕变参数 - $R, T$:气体常数与绝对温度
4. 松弛时间常数与半衰期
定义松弛半衰期 $t_{1/2}$,即应力下降至初始值一半所需的时间。由 $\sigma(t_{1/2}) = 0.5 \sigma_0$ 代入解得:
温度敏感性:由于含有 $\exp(Q/RT)$ 项,温度每升高 $20 \sim 30^\circ\text{C}$,松弛半衰期可能缩短至原来的十分之一甚至更短。这解释了为何 500°C 以上的松弛可在数小时内完成。
5. 与碟形弹簧设计的关联
- 初始应力控制:从公式可见,$\sigma_0$ 越高,分母中的 $\sigma_0^{n-1}$ 项越大,松弛速率急剧加快。因此,高温碟簧设计中,应适当降低工作应力(即降低最大压缩量 $s/h_0$),以换取更长的松弛寿命。
- 材料选择:
- 若工作温度 $>400^\circ\text{C}$,必须选用 $A$ 值极小、$Q$ 值极大的材料(如 Inconel 718),使分母中的指数项占优。
- H13(热作模具钢)适用于 $400 \sim 500^\circ\text{C}$,但需严格核算长期松弛损失。
- 强压处理:在高于工作温度下进行预压,使前期快速松弛在出厂前完成,可显著提高服役期内的应力稳定性。
6. 计算示例
已知:H13 碟簧,工作温度 500 °C (773 K),初始应力 $\sigma_0 = 800\ \text{MPa}$,弹性模量 $E = 190,000\ \text{MPa}$。Norton 参数:$A = 5.0 \times 10^{-20}$,$n = 5.5$,$Q = 320,000\ \text{J/mol}$。
计算松弛半衰期: - $\sigma_0^{n-1} = 800^{4.5} = \exp(4.5 \cdot \ln 800) \approx \exp(4.5 \cdot 6.6846) \approx \exp(30.08) \approx 1.15 \times 10^{13}$ - $2^{n-1} - 1 = 2^{4.5} - 1 \approx 22.63 - 1 = 21.63$ - $E \cdot A = 190,000 \times 5.0 \times 10^{-20} = 9.5 \times 10^{-15}$ - $\frac{Q}{RT} = \frac{320,000}{8.314 \times 773} \approx 49.7$ - $\exp(-Q/RT)$ 在分母中为 $\exp(Q/RT)$,计算得 $\exp(49.7) \approx 3.6 \times 10^{21}$
分母部分:$4.5 \times 9.5 \times 10^{-15} \times 1.15 \times 10^{13} \approx 4.5 \times 9.5 \times 1.15 \times 10^{-2} \approx 4.5 \times 10.925 \times 10^{-2} \approx 0.492$
换算为小时:$\approx 4.4 \times 10^{19}\ \text{h}$,这似乎极长,说明参数选择可能未反映实际高温蠕变数据(例如 $A$ 值可能偏小)。实际上 H13 在 500 °C 的蠕变参数需从材料标准获取,本例仅演示计算流程。若激活能略低(如 280 kJ/mol),结果会急剧变化。工程中,H13 在 500 °C、800 MPa 下的半衰期大约在 数千小时 量级,必须通过实测数据验证。
总结:蠕变-松弛耦合通过求解 Norton-Bailey 本构下的微分方程,精确预测了碟簧在恒定位移下的应力衰减。松弛速率高度依赖于温度、初始应力和材料蠕变抗力,是高温碟簧长周期可靠性设计的核心。