热装应力模拟器 返回
结构力学模拟器

热装应力模拟器 — 厚壁圆筒的 Lamé 解

基于拉梅厚壁圆筒解,可视化实心轴与外套筒热装时产生的接触压力、套筒周向应力、轴内压缩应力与摩擦传递扭矩。改变过盈与尺寸即可直观判断应力是否在安全范围内。

参数设置
接触界面直径 D_b
mm
套筒外径 D_o
mm
过盈量 δ(直径基)
mm
弹性模量 E
GPa

假设实心轴(r_i = 0)、配合长度 L = 50 mm、摩擦系数 μ = 0.15(钢-钢 干燥)。轴与套筒为同种材料(同 E 与 ν)。

计算结果
界面接触压力 p
套筒内壁最大 σ_θ(拉)
轴内 σ_r = σ_θ(压缩)
最大传递扭矩 T(μ=0.15)
热装断面与径向应力分布

左:热装断面(蓝=轴、橙=套筒)/右:σ_r(r)(橙)与 σ_θ(r)(蓝)沿径向的分布,界面 r=r_b 处 σ_θ 跳变

理论与主要公式

同种材料(同 E 与 ν)的实心轴(半径 r_i=0 至 r_b)与外套筒(r_b 至 r_o)的热装由拉梅(Lamé)厚壁圆筒解描述。δ 为直径基过盈量。

界面接触压力 p(实心轴情形):

$$p = \frac{E\,\delta}{2\,r_b}\cdot\frac{r_o^2 - r_b^2}{2\,r_o^2}$$

套筒内(r_b ≤ r ≤ r_o)的径向应力与周向应力:

$$\sigma_r(r) = \frac{p\,r_b^2}{r_o^2 - r_b^2}\left(1 - \frac{r_o^2}{r^2}\right),\quad \sigma_\theta(r) = \frac{p\,r_b^2}{r_o^2 - r_b^2}\left(1 + \frac{r_o^2}{r^2}\right)$$

套筒内壁(r=r_b)最大拉应力与实心轴内的一致应力:

$$\sigma_{\theta,\max} = p\,\frac{r_o^2 + r_b^2}{r_o^2 - r_b^2},\qquad \sigma_r = \sigma_\theta = -p\ (\text{shaft})$$

摩擦传递的最大扭矩(μ 为摩擦系数,L 为配合长度):

$$T = 2\pi\,\mu\,p\,r_b^2\,L$$

默认参数(D_b=100, D_o=150, δ=0.1 mm, E=210 GPa, L=50 mm, μ=0.15)下,p ≈ 58.3 MPa,σ_θ,max ≈ 151.7 MPa(拉),轴内 −58.3 MPa(压缩),T ≈ 6.87 kN·m。

什么是热装应力模拟器

🙋
「热装」不就是把套筒加热膨胀套到轴上吗?为什么冷却后能固定得那么牢?
🎓
大致来说,套筒冷却时想要收缩,但被略大的轴挡住了。这部分「想缩却缩不回去」的量就是过盈 δ,全部转化为界面处的弹性应变,从而产生接触压力 p。把上面模拟器的「过盈量 δ」从 0.05 拉到 0.20 mm,p 的卡片几乎线性上升——拉梅厚壁圆筒解里 p 关于 δ 就是一次的。
🙋
那是不是过盈越大,能传的扭矩越大?
🎓
原理上是。传递扭矩 $T = 2\pi\mu p r_b^2 L$ 与 p 成正比。但请看右边的蓝色曲线——套筒内壁的周向应力 σ_θ。过盈翻倍,p 翻倍,σ_θ 也按相同比例上升。默认参数下 σ_θ ≈ 152 MPa,若把 δ 加到 0.20 mm,σ_θ 会超过 300 MPa,远高于普通碳素钢的许用拉应力(约 200〜250 MPa),套筒会屈服甚至胀裂。「扭矩能力」与「套筒破坏」始终是对拉关系。
🙋
那把套筒做厚一点呢?把 D_o 拉大应力应该会降下来吧?
🎓
观察很准。把 D_o 从 150 拉到 300 mm,σ_θ 确实下降了,但下降速度越来越慢。由 $\sigma_{\theta,\max} = p\frac{r_o^2+r_b^2}{r_o^2-r_b^2}$ 可知,r_o → ∞ 时 σ_θ,max 收敛到 p——无论套筒多厚,内壁周向应力都不会小于接触压力。工程上 D_o ≈ 1.5~2.0×D_b 是「重量与应力的甜点」,默认 D_o/D_b = 1.5 正是这个经济区间。
🙋
弹性模量 E 改变会怎么样?比如黄铜(E≈100 GPa)和钢(E≈210 GPa)的差别?
🎓
p 与 E 成正比。同样过盈下,黄铜套筒只能产生一半的接触压力,传递扭矩也减半。有趣的是 σ_θ/p 只取决于几何(r_o、r_b),与 E 无关,所以「相对许用应力的安全率」并不随 E 变化。结论是:弹性模量较小的材料要靠加大过盈才能产生同样的接触压力。

常见问题

两者最终的接触压力都由同一个拉梅厚壁圆筒解给出,但装配方式不同。压入用液压机轴向推入,对过盈和表面光洁度有限制,容易擦伤配合面。热装将套筒加热到 150~300 °C 使其膨胀,几乎无装配力即可套上,冷却后即获得大过盈牢固夹紧。大直径、大扭矩传动多用热装,因为装配力不随过盈量增加。
所需热膨胀 = 「过盈 δ + 装配余量 0.05~0.10 mm」。钢的线膨胀系数 α ≈ 12×10⁻⁶ /K,所以温升 ΔT = (δ + 余量)/(α·D_b)。例如 D_b = 100 mm、δ = 0.10 mm、余量 0.05 mm,ΔT ≈ 125 K,套筒加热到约 150 °C 即可。对调质件,加热温度应远低于回火温度,一般以 200 °C 为上限。
本模拟器假设实心轴(r_i=0)。空心轴需使用一般化拉梅式 p = (E·δ)/(2·r_b) × (r_o²−r_b²)(r_b²−r_i²) / (2·r_b²·(r_o²−r_i²))。空心轴增加了内侧的弹性顺度,所以相同过盈下接触压力较低。同时空心轴内壁(r=r_i)会出现附加的压缩周向应力,应力校核需要同时检查轴和套筒。
高速旋转件(涡轮转子、飞轮等)的套筒因离心力向外胀大,会减小有效过盈。松脱量按 ρω²r_b² 的量级增长,因此高速机械须在静态过盈上预留较大裕量,或将「离脱速度」(接触压力消失的转速)设计得明显高于工作转速。本模拟器仅给出静态值,高速场合需另作离心补偿。

实际工程应用

铁路轮对:热装的经典案例。车轮直接热装到车轴上,不用键也不用焊接,就能传递牵引和制动扭矩。长期使用过程中过盈会出现缓慢「蠕变」,定期用超声波探伤和压痕测量来评估配合状态。新干线级别的高速车辆对过盈和公差要求更严格。

大型齿轮与联轴器的轴连接:轧钢机大齿轮、船舶推进轴等多采用无键热装,避免键槽处的应力集中。模拟器给出的最大传递扭矩通常除以 1.5~2.0 的安全系数。对超大直径,常采用液压扩张联轴器(如 SKF 油膜联轴器)以方便拆装。

涡轮转子组装:蒸汽轮机、燃气轮机的轮盘热装到轴上时,需考虑运行时的离心力与热膨胀。常温过盈过大会在启动热应力下导致轮盘裂纹,过小则在高转速下脱离引起振动。本模拟器的静态应力评估正是这种设计研究的起点。

滚动轴承外圈与轴承座:滚动轴承外圈以过盈配合压入轴承座,由接触压力产生的周向应力会令外圈向内收缩,使内部游隙减小。若过盈过大会缩短轴承寿命,ISO/JIS 按配合等级规定了各种允许过盈范围。

常见误解与注意事项

最常见的误解是把热装当成一致黏接的均匀联接。理想分析中接触压力是均匀的,但实际配合在配合长度 L 的两端会出现「端部效应」。套筒边缘的周向应力可能比中部高 1.3〜2.0 倍,疲劳裂纹几乎总是从端部开始。设计上常在套筒端做倒角,或对配合长度做局部锥形过渡以分散端部应力。本模拟器只给出「中部断面」的应力,需要意识到这一点。

第二个常见错误是混淆直径基与半径基过盈。本模拟器及 JIS/ISO 公差表均使用「直径基 δ = D_shaft − D_hole」。某些教科书使用「半径基 δ_r = δ/2」,两者混用会带来 2 倍的接触压力估算误差。由图纸公差计算时务必确认使用的是「直径差」,再核对所参考公式。

最后,仅用「接触压力不超过屈服」做单一判定是危险的。周向应力 σ_θ 是拉应力,最大值往往是接触压力的若干倍。默认几何(D_o/D_b = 1.5)下 σ_θ,max ≈ 2.6×p,p = 100 MPa 就能让 σ_θ,max 达到 260 MPa。需要校核的是套筒内壁的 σ_θ,max 是否小于材料的许用拉应力。对于铸铁等脆性材料,应使用最大主应力而非 von Mises 应力来判断,更偏安全。