本节摘要:参考系一转,牛顿第二定律就"缺账":旋转系里补上离心力,物体相对参考系再运动就还得补科里奥利力。两笔惯性力没有施力者,是坐标选择的记账成本,但账面效应真实可测——离心机、台风旋向、傅科摆全在替它们签字。本节给出两笔账的公式、量级与地球现场的兑现。
第一章审过参考系的资格:牛顿第二定律只在惯性系里精确成立。可地球自转、车厢转弯,舞台本身就在动,站上去的人怎么办?办法是补账:把参考系的加速度"折算成假想的力"加进方程,牛顿定律照用不误。折算出两笔:参考系旋转带来的离心力,大小为质量乘角速度平方乘半径,指向离轴方向;物体在旋转系里还有相对速度时,再补一笔科里奥利力,大小为两倍质量乘角速度乘速度乘纬度正弦,方向垂直于速度(北半球指向运动方向的右侧)。
两笔力的共同点是"没有施力者"——你找不到哪根弹簧或哪块肌肉在推。它们是惯性在非惯性账本上的体现:惯性系里看,物体只是想沿直线走;旋转系里看,这条直线被记录成弯曲的轨迹,于是账本上必须虚构一笔力来圆账。认清这一点,公式就不会用错场合。
先看一小一大两台病例。小病例:转盘上放半公斤的滑块,离轴零点八米,转盘每秒转零点四八圈——滑块随盘做圆周运动需要向心力,站在盘上的人却看到滑块静止,这笔向心力在旋转账本上以离心力名义出现,由摩擦支付。大病例:地球自转带来的离心力在地表各处打了折扣,赤道最足、两极为零。
# 手术一:离心力 F = m*omega^2*r;转盘与地球两笔账 import math m, r, f = 0.5, 0.8, 0.48 # 滑块 kg、半径 m、转速 圈/s omega = 2 * math.pi * f F_cent = m * omega**2 * r print("角速度 =", round(omega, 4), "rad/s") # 3.0159 print("离心力 =", round(F_cent, 4), "N") # 3.6383 mu_need = F_cent / (m * 9.8) print("维持不动需求摩擦系数 =", round(mu_need, 4)) # 0.7425 # 地球自转:omega = 7.292e-5 rad/s,赤道离心加速度 omega^2*R w_earth, R_earth = 7.292e-5, 6.371e6 a_cent = w_earth**2 * R_earth print("赤道离心加速度 =", round(a_cent, 5), "m/s²") # 0.03388 print("与重力之比 =", round(a_cent / 9.8, 5)) # 0.00346 # 纬度 45 度处再打 cos(45) 的折扣 print("45 度纬度离心/重力 =", round(a_cent * math.cos(math.radians(45)) / 9.8, 5)) # 0.00244
转盘上的滑块离临界只差一步:需求摩擦零点七四,换成橡胶面才稳得住,转盘再快半圈滑块立刻甩出——离心分离机的转速设计、游乐设施的刺激预算,算的都是这一笔。地球那笔则人人头顶都在享受:赤道的离心折扣约千分之三,这就是同一物体在赤道称得最轻的物理原因——称出来的"重量"已经扣掉了离心账。
物体在旋转系里一走动,第二笔账就开始计息。标准病例:北纬四十五度,一枚炮弹以每秒三百米平飞,科里奥利加速度多大、往哪边偏。再放大到天气尺度:大气里任何长距离流动都被科里奥利力持续掰弯,北半球一律右偏——台风逆时针旋转、信风的风向偏转、傅科摆的进动,全是同一笔账的宏观兑现。
# 手术二:科氏加速度 a = 2*omega*v*sin(纬度),方向北半球右偏 import math w_earth = 7.292e-5 v, lat = 300.0, 45.0 a_cor = 2 * w_earth * v * math.sin(math.radians(lat)) print("科氏加速度 =", round(a_cor, 5), "m/s²") # 0.03094 # 与重力之比:千分之三量级,远程弹道必须修正 print("与重力之比 =", round(a_cor / 9.8, 5)) # 0.00316 # 飞行 60 s 的横向偏移(近似匀加速):0.5*a*t^2 t = 60.0 drift = 0.5 * a_cor * t * t print("一分钟横向偏移 =", round(drift, 1), "m") # 55.7 # 傅科摆进动周期 = 23.93 h / sin(纬度) T_foucault = 23.93 / math.sin(math.radians(45.0)) print("45 度傅科摆进动周期 =", round(T_foucault, 2), "小时") # 33.84 # 赤道上 sin(0)=0:傅科摆不进动,科氏力水平分量为零 print("赤道进动周期 =", 23.93 / math.sin(math.radians(0.0001)), "≈ 无穷大")
一分钟偏移五十六米——远程炮击不修正科氏账就永远打不中,这正是第一次世界大战海军射击表的经典修正项。傅科摆的进动周期在巴黎纬度约三十二小时,摆在替站在地上的人作证:转的不是摆,是脚下的地球。
⚠️ 常见坑:给"离心力"找施力者。它没有出处——在惯性系的账本上根本不存在这一笔,只有向心力;两套账本各画各的受力图,混着画必然多出一笔幽灵力。
离心力最成功的工程兑现是离心分离:混合液在转鼓里旋转,重颗粒需要的向心力大,"甩"向外壁;轻的留在内侧——分离系数(离心加速度与重力之比)决定分离效率。家用洗衣机的脱水、实验室的离心机、牛奶的奶油分离,同一套物理,参数差三个量级。
# 离心分离:分离系数 = omega^2*r/g import math g = 9.8 for rpm, r, tag in ((1200, 0.12, "洗衣机脱水"), (15000, 0.08, "实验室离心机")): omega = 2 * math.pi * rpm / 60 ratio = omega**2 * r / g print(f"{tag}: {rpm} rpm 半径 {r} m 分离系数 = {ratio:,.0f} g") # 洗衣机脱水: 1200 rpm 半径 0.12 m 分离系数 = 193 g # 实验室离心机: 15000 rpm 半径 0.08 m 分离系数 = 20,142 g # 两万倍重力把细胞器都压到管底——半径小就靠转速补 print("同样 0.08 m,想要 40000 g 需要 rpm =", round(60*math.sqrt(40000*g/0.08)/(2*math.pi))) # 需要 rpm = 21138
洗衣机上百倍重力已经够把衣物甩干,实验室两万倍重力才能压沉细胞器——分离系数对转速平方依赖,转速翻倍、效率翻两番,这就是离心机越造越快的物理理由。
要点复盘:
三界病例到此全部清点。下一章进入收尾工序:把全套流程交给代码,用数值方法与极限检验给每台手术上最后一道保险。