热套应力模拟器 返回
结构分析模拟器

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

实心轴与外套筒热套配合中产生的接触压、套筒内周向应力、轴内压缩应力、摩擦可传递扭矩使用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处σ_θ跳跃

理论·主要公式

使用Lamé(拉梅)厚壁圆筒解处理同一材料(相同E和ν)的实心轴(半径r_i=0~r_b)和外套筒(r_b~r_o)的热套。δ为直径基础的过盈量。

界面接触压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{轴})$$

摩擦最大可传递扭矩(μ为摩擦系数,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的数值会怎么变——应该会几乎成比例地跳跃。根据Lamé厚壁圆筒解,p确实与δ呈一次线性关系。
🙋
那过盈量越大能传递的扭矩就越多呗?
🎓
原理上是这样。可传递扭矩公式是$T = 2\pi\mu p r_b^2 L$,与p成正比。但你看右边图表中那条蓝线——套筒内面的周向应力σ_θ。当过盈增大时,p上升,σ_θ也会以相同比例跳升。在默认参数下σ_θ约152 MPa,如果把δ加倍到0.20 mm,σ_θ就会超过300 MPa。钢材的许用应力一般是200~250 MPa,超过了套筒就会塑性变形,甚至断裂。所以设计时存在「可传递扭矩」与「套筒破坏风险」的权衡。
🙋
如果把套筒做厚一点,应力能降下来吗?试着把外径D_o的滑块从150 mm 拉到300 mm 看看会怎样。
🎓
这个观察很好。确实能看到σ_θ下降。但增幅越来越小。看公式$\sigma_{\theta,\max} = p\frac{r_o^2+r_b^2}{r_o^2-r_b^2}$,当r_o→∞时,σ_θ,max→p,有个渐近极限。也就是说,套筒再怎么厚,内面应力也不会低于接触压p本身。在工程实践中,通常D_o/(D_b)取1.5~2.0为「强度与重量的好平衡点」。模拟器的默认1.5就是参考这个原则。
🙋
最后问一个,如果改变杨氏模量E,比如黄铜(E≈100 GPa)和钢(E≈210 GPa)有什么不同?
🎓
E与p成正比,起作用很直接。用同样过盈量,黄铜套筒的接触压会是钢的一半,可传递扭矩也是一半。不过有趣的是,应力比σ_θ/p只取决于几何形状(r_o、r_b),与E无关。所以改变E不影响「材料强度余度」的评估。只是说,用低模量材料时,如果要传递同样的扭矩,过盈量得取得更大。这就是为什么E小的材料在同样负载下需要更大过盈的原因。

常见问题

最终接触压和过盈量的计算都用Lamé厚壁圆筒解,但组装方法不同。压入是用液压机沿轴向把套筒推进去,易在接触面留下划伤,而且过盈大时压入力巨大。热套是把套筒加热到150~300℃膨胀后,以间隙配合的状态套进轴,冷却后收缩——这样组装力接近零,大过盈量也容易实现。大直径、大扭矩传递时优先选热套。
所需膨胀量是「过盈量δ+组装余隙(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℃。但若套筒已经过回火处理,加热温度要低于回火温度,通常不超过200℃为宜。
本模拟器默认实心轴(r_i=0)。对于中空轴,通用Lamé公式为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²阶增。高速机械需要取更大初始过盈,或者设计「离脱转速」(接触压消失的转速)要远高于许可工作转速。本模拟器计算的是静止状态下的值,高速旋转的设计需另外进行遐心力补正。

实际应用案例

铁路车轮与车轴:热套的典型应用就是铁道车辆车轮。车轮直接热套在车轴上,无需键或焊接就能传递巨大的驱动和制动扭矩。长期使用中过盈量会略微减小(「蠕变」现象),需要定期用超声波探伤或压痕测量检查配合状态。高铁级列车对过盈精度要求更高。

大型齿轮、联轴器的轴安装:冶金轧钢机、船舶推进轴的大齿轮多采用无键热套,避免键槽处的应力集中。为便于组装和分解,有时配合SKF油膜联轴器(油压膨胀式继手)等特殊装置。

汽轮机、燃气轮机转子装配:蒸汽汽轮机和燃气轮机的动叶片盘装配到轴上时,运行中遐心力和热膨胀都要考虑。冷态过盈取得过大容易在启动时因热应力导致盘片断裂;过小则高速时易发生分离和振动。这里模拟器计算的静态应力只是设计的起点,还要做详细的热应力和动力学分析。

滚动轴承外圈与轴承座配合:滚动轴承外圈在轴承座内的「过盈配合」会给轴承施加额外的周向预紧。过大的过盈会压缩轴承内部间隙,导致早期磨损。各轴承型号对配合等级都有JIS规定的过盈上下限。

常见误区与注意

最常见的误区是认为「热套像胶水一样均匀粘合」。实际上界面接触压理想情况下均匀分布,但配合长L的端部会出现「端部效应」应力集中。套筒端部的σ_θ可能是中间的1.3~2.0倍,疲劳破坏几乎总是从端部开始。设计中常做法是对套筒端部做倒角、或使端部逐渐变细(锥形)来分散应力。本模拟器显示的是「中间截面」的应力值,需留意这一点。

其次常见的是过盈量用「直径基」还是「半径基」混淆。本模拟器和JIS公差表用的都是直径基:δ = D_轴 − D_孔。但某些教科书用「半径基」δ_r = δ/2来推导公式,两者混用会导致2倍的计算错误。从图纸公差推算时一定要确认是「直径差」,并对照参考书中的公式定义。

第三个危险的简化是「只要接触压不超过屈服应力就安全」。周向应力σ_θ是拉应力,往往是接触压的数倍。在默认几何(D_o/D_b=1.5)下,σ_θ,max≈2.6×p,所以p=100 MPa时σ_θ,max已达260 MPa。需要检验的是套筒内面σ_θ,max在许用拉应力以内。对脆性材料(铸铁等)应改用最大主应力判定而不是Von Mises应力,更加安全。

使用指南

  1. 输入实心轴外径Db(mm)、套筒内径Do(mm),设置过盈量Δ(μm)
  2. 指定材料的杨氏模量E(GPa)。例如调质钢为206 GPa、铝青铜为110 GPa
  3. 模拟器会从Lamé厚壁圆筒解自动计算接触压p、套筒周向应力σθ、轴内压缩应力、可传递扭矩T(μ=0.15)
  4. 检验界面接触压是否超过材料屈服应力、套筒周向应力是否在许用值以内

具体计算例

调质钢轴(Db=40 mm、E=206 GPa)配铝青铜套筒(Do=50 mm、E=110 GPa),过盈Δ=30 μm,计算结果为接触压p≈183 MPa、套筒内面周向应力σθ≈456 MPa、轴内压缩应力≈183 MPa,最大可传递扭矩T≈1382 N·m。若套筒拉伸强度500 MPa,则安全。

工程实务注意事项

  1. 过盈过大会导致轴或套筒塑性变形,Lamé解失效。直径Db≤50 mm的轴通常限制Δ≤50 μm
  2. 热套时温差会产生热膨胀,改变实际过盈。轴加热80℃时,钢的线膨胀系数11.5×10⁻⁶/℃,实效过盈会变化
  3. 若严格区分两种材料的泊松比,可用Do代入重新计算
  4. 可传递扭矩假定摩擦系数μ=0.15。更保险的设计应用μ=0.10~0.12重新评估