本节摘要:把解剖对象从质点换成定轴转动的刚体,整套牛顿机器几乎逐项平移:质量换转动惯量,力换力矩,ma 换角加速度乘转动惯量。本节讲清转动惯量怎么算、平行轴定理怎么搬,用飞轮起动与滑轮连体两台手术演示转动方程与能量账本的对账。
前三节的所有手术对象都是"点":力可以画在一点上,位移一个箭头说完。刚体转起来之后,同一根力在不同位置效果天差地别——推门框边缘轻松,推门轴白费劲。量纲必须换血:位置用角度记,快慢用角速度记,力对转动的贡献用力矩记,而抵抗角加速度的"惯性质量"就是本节主角转动惯量。它不只看质量多少,更看质量离轴多远——同样一吨,铺在轮缘与聚在轴心,转动的脾气相差成百上千倍。定轴转动的方程因此长成了牛顿第二定律的镜像:合力矩等于转动惯量乘角加速度。五步流程照用,只是第四步标完力之后多一道工序:把每根力折算成力矩,顺便记下力臂。
先备好转动惯量的"标准件目录":细环绕轴心是质量乘半径平方,匀质圆盘是二分之一质量乘半径平方,细杆绕端点是三分之一质量乘长平方。离轴越远贡献越大,因为转动惯量里质量带着距离的平方。标准病例:冲床飞轮,质量六十千克、半径零点四米的匀质圆盘,电机从静止施加二十五牛·米的恒定力矩。转动方程给出角加速度,随后的角度、动能全部是它的积分。
# 手术一:飞轮恒力矩起动 M = I*alpha import math m, R, M_t = 60.0, 0.4, 25.0 # 圆盘质量 kg、半径 m、驱动扭矩 N·m I = 0.5 * m * R * R # 匀质圆盘绕中心轴 alpha = M_t / I print("转动惯量 I =", I, "kg·m²") # 4.8 print("角加速度 =", alpha, "rad/s²") # 5.208333333333333 for t in (1.0, 3.0): omega = alpha * t theta = 0.5 * alpha * t * t Ek = 0.5 * I * omega * omega print(f"t={t:.0f}s 角速度={omega:7.3f} rad/s 转角={theta:7.2f} rad 动能={Ek:8.1f} J") # t=1s 角速度= 5.208 rad/s 转角= 2.60 rad 动能= 65.1 J # t=3s 角速度= 15.625 rad/s 转角= 23.44 rad 动能= 585.9 J # 对账:动能应等于扭矩做的功 M_t*theta print("三秒末扭矩做功 =", M_t * 0.5 * alpha * 9, "J") # 585.9375
账目闭合得干净:三秒末的转动动能恰好等于扭矩乘转角,功与动能的转动版定理在数值上自我验证。顺带读出工程含义——飞轮转得越快,再加速同样的份额要吃掉平方级的功,冲床靠它在冲压瞬间放能、间隙期间蓄能,选型全围绕这份平方关系。
真实系统里的轴常不在质心,或在系统里混着转动的轮与平动的块。前者用平行轴定理搬运:绕任一平行轴的转动惯量等于绕质心轴的值加上质量乘两轴间距平方。后者是经典病例——定滑轮不可忽略质量,绳两端挂重物:绳不打滑,滑轮边缘的线加速度与物块一致,滑轮的力矩方程与两物块的平动方程联立,用能量账本验总账。
# 手术二:桌面滑轮连体。物块 A(4kg) 在桌面,绳过滑轮吊物块 B(3kg) # 滑轮匀质圆盘 m=2kg, R=0.1m,桌面无摩擦,g=9.8 import math mA, mB, mp, R, g = 4.0, 3.0, 2.0, 0.1, 9.8 Ip = 0.5 * mp * R * R # 滑轮转动惯量 # 三本账:A: T1 = mA*a;B: mB*g - T2 = mB*a;滑轮: (T2 - T1)*R = Ip*(a/R) a = mB * g / (mA + mB + Ip / (R * R)) T1 = mA * a T2 = mB * (g - a) print("加速度 a =", round(a, 4), "m/s²") # 3.675 print("桌面端绳张力 T1 =", round(T1, 3), "N") # 14.7 print("吊重端绳张力 T2 =", round(T2, 3), "N") # 18.375 # 能量对账:B 下落 1 m 后,由 v^2 = 2*a*h 求末速度的平方 h = 1.0 v2 = 2 * a * h KE_blocks = 0.5 * (mA + mB) * v2 # 两个物块的平动动能 KE_pulley = 0.5 * Ip * (v2 / R**2) # 滑轮转动动能(边缘线速度与物块相同) KE_all = KE_blocks + KE_pulley print("物块动能 =", round(KE_blocks, 3), "J") # 25.725 print("滑轮动能 =", round(KE_pulley, 3), "J") # 3.675 print("动能总账 =", round(KE_all, 3), "J 重力做功 =", round(mB * g * h, 3), "J") # 动能总账 = 29.4 J 重力做功 = 29.4 J —— 能量定理精确闭合
能量对账精确闭合:重力做功二十九点四焦,一分不差地分进两个物块的平动动能与滑轮的转动动能。滑轮有质量时,"绳两端张力不等"是账目里最容易漏掉的一笔:T1 与 T2 相差的正是驱动滑轮转动的那份力矩。
⚠️ 常见坑:把滑轮当"无质量滑轮"顺手抄 T1 等于 T2。滑轮质量不可忽略时两端张力必不相等,差值由力矩方程定——漏掉这条差值,能量对账永远差一截。
两台手术之外,再交代两件实战要务。其一是边界:定轴转动方程只认"轴固定不动"的病例,轴自己也在加速的(如滚动的轮子)要先做运动学分析,把接触点速度约束补进方程组——这也是下一章拉格朗日语言的主场,约束越多它越省事。其二是组合件:真实飞轮常是盘、毂、辐条的拼装,总转动惯量按"分件计算、直接相加"处理,各件各自选对公式再求和,比硬积一个复杂形状省力得多。
# 组合飞轮:钢盘(30 kg,R=0.3 m)加边缘钢环(20 kg,R=0.35 m) import math I_disk = 0.5 * 30.0 * 0.3**2 # 盘绕自身中心 I_ring = 20.0 * 0.35**2 # 环全部质量在轮缘 I_total = I_disk + I_ring print("盘的份额 =", round(I_disk, 4), "kg·m²") # 1.35 print("环的份额 =", round(I_ring, 4), "kg·m²") # 2.45 print("总转动惯量 =", round(I_total, 4), "kg·m²") # 3.8 # 质量占比 60% 的盘只贡献 35.5% 的转动惯量——距离平方的权重再次现形 print("环质量占比 =", 20.0/50.0, " 惯量占比 =", round(I_ring/I_total, 4)) # 0.4 0.6447
四成质量的环扛下六成四的惯量,飞轮设计"质量往边缘堆"的口诀在这两行输出里有了定量版本。分件相加的思路同样适用于平行轴定理的多次搬运:每搬一次加一份质量乘轴距平方,错了哪一段,能量对账立刻露馅。
要点复盘:
至此质点与刚体的矢量力学告一段落。下一章换一套语言:不盯着每根力,改用系统的广义坐标与能量——约束越多,这套语言越锋利。