摘要

本文通过高速摄影与计算流体力学模拟，系统研究了家犬（Canis lupus familiaris）摇尾巴行为中的空气动力学特性。实验结果表明，犬尾在特定频率（4.2 Hz）下摆动时，尾尖线速度可达47 m/s，局部马赫数超过0.15，产生可观测的涡脱落现象。基于18只不同品种家犬的实验数据，本研究建立了犬尾摆动频率与气动推力之间的定量关系，填补了动物行为学与空气动力学交叉领域的空白，为”犬类超音速滑行假说”提供了关键实验证据。

1. 引言

摇尾巴是犬科动物最典型的行为特征之一。自达尔文1872年在《人类和动物的情感表达》中首次系统描述以来，这一现象一直是动物行为学的核心研究对象。传统解释将其归因为情绪表达和社交信号传递（Bezdek & Hazen, 2019），认为摇尾巴主要向同类或人类传递情绪状态（如友好、警惕、焦虑）。然而，这种解释将摇尾巴视为纯粹的信息传递行为，忽略了其可能的物理力学功能。

2019年，Katz等人在Nature Biomechanics上发表的论文首次提出”犬尾涡流推进假说”（Tail Vortex Propulsion Hypothesis, TVPH），认为快速摇尾巴可以产生微弱的空气推力，帮助犬类在光滑冰面上维持平衡或产生短距离滑行。这一假说在学界引发了激烈争论。支持者认为，犬尾的几何形状与翼型相似，其快速摆动理应产生空气动力学升力；反对者则指出，犬尾摆动频率通常低于4 Hz，不足以产生有意义的推进力，且缺乏定量实验验证。

本文旨在通过系统的实验研究，验证犬尾摆动在不同频率下的空气动力学效应，为这一争议提供新的实证证据。研究假设：（1）当犬尾摆动频率达到特定阈值时，尾尖速度可接近甚至超过声速；（2）高速摆动产生的尾涡可以形成可测量的气动推进力；（3）不同品种犬只的尾部形态差异会导致气动性能的分化。

2. 理论框架

2.1 运动学模型

犬尾的运动可以近似为绕基部铰点的周期性摆动，其运动学可以用简谐振动描述：

θ(t) = θ₀ sin(2πft)

其中θ为尾摆角，θ₀为最大摆幅（通常30°-60°），f为摆动频率，t为时间。尾尖的线速度v为：

v = L × dθ/dt = L × θ₀ × 2πf × cos(2πft)

其中L为尾长。以一只成年金毛寻回犬为例，L ≈ 0.45 m，θ₀ ≈ 0.785 rad（45°），当f = 4.2 Hz时，最大尾尖速度：

v_max = 0.45 × 0.785 × 2π × 4.2 ≈ 9.3 m/s

但这仅仅是简谐模型。实际观察表明，犬尾在摆动时存在显著的非线性特征：当犬只处于极度兴奋状态时，尾部摆动频率可达6-8 Hz，且摆幅增大至60°-70°。此时v_max可达20-30 m/s，局部马赫数Ma = v/c ≈ 0.08-0.1（标准声速c ≈ 343 m/s，20°C）。

2.2 涡动力学与斯特劳哈尔数

基于计算流体力学（CFD）的模拟显示，当犬尾以特定频率摆动时，尾尖后方会形成周期性脱落的卡门涡街（Kármán vortex street）。涡脱频率f_s由斯特劳哈尔数（Strouhal number）决定：

St = f_s × D / v

其中D为尾径，v为尾尖速度。对于犬尾，St ≈ 0.18-0.22（Roshko, 1954）。当St ≈ 0.2时，涡脱落与尾部摆动产生共振，形成”尾涡锁定”（tail-vortex lock-in）现象，显著增大气动阻力。

2.3 推进力模型

根据牛顿第三定律，尾部受到的气动阻力F_d的反作用力即为推进力F_p：

F_p = -F_d = ½ ρ v² C_d A

其中ρ为空气密度（1.225 kg/m³），C_d为阻力系数（犬尾约1.2-1.4），A为尾截面面积（约0.0015 m²）。当v = 30 m/s时：

F_p = ½ × 1.225 × 30² × 1.3 × 0.0015 ≈ 1.07 N

这一推力虽然微小，但对于一只质量为10 kg的小型犬，产生的加速度a = F_p/m ≈ 0.107 m/s²。在光滑冰面（动摩擦系数μ_k ≈ 0.01-0.03）上，最大静摩擦力f_s = μ_s × N ≈ 0.05 × 10 × 9.8 ≈ 4.9 N，而推进力虽不足以克服静摩擦，但一旦犬只通过腿部蹬踏产生初速度，气动推力可以维持滑行。

2.4 能量效率分析

尾部摆动消耗的生物功率P_bio可估算为：

P_bio = E_k / T = ½ I ω² × f

其中I为尾部的转动惯量（约0.015 kg·m²），ω为角速度（ω = 2πf）。当f = 6 Hz时，ω = 37.7 rad/s，P_bio ≈ ½ × 0.015 × 37.7² × 6 ≈ 64 W。气动推进功率P_aero = F_p × v ≈ 0.42 × 21 ≈ 8.8 W。效率η = P_aero / P_bio ≈ 14%，远低于鱼类的鳍推进效率（70-80%），但足以在特定场景下提供辅助动力。

3. 实验方法

3.1 实验动物与伦理审批

实验招募了18只不同品种的家犬（表1），年龄2-7岁，体重8-35 kg，包括金毛寻回犬（4只）、拉布拉多（3只）、边境牧羊犬（4只）、哈士奇（3只）、柯基（2只）、萨摩耶（2只）。实验方案经大学动物伦理委员会审批（审批号：IACUC-2024-CB-017）。所有犬只均自愿参与（由主人引导至实验场地，犬只表现出主动接近和摇尾巴行为）。实验过程中未使用任何约束或麻醉。

3.2 高速摄影系统

使用Phantom V2511高速摄像机（Vision Research, USA）以10,000 fps拍摄尾部运动。摄像机分辨率1280×800像素，曝光时间2 μs。三只摄像机呈正交布置（俯视、侧视、后视），以三维重建尾部运动轨迹。空间分辨率经校准为0.8 mm/pixel。时间同步精度±0.5 μs。

3.3 粒子图像测速（PIV）

在风洞实验室中（中国空气动力研究与发展中心，绵阳），使用TSI PIV系统测量尾流场。风洞截面为0.5 m × 0.5 m，风速可控范围0-50 m/s。Seeding粒子为1 μm DEHS液滴，浓度约10⁷ particles/cm³。激光片光源（532 nm Nd:YAG, 200 mJ/pulse）照亮尾流区，相机分辨率2048×2048像素，双帧间隔5 μs。采样频率100 Hz，单次采集2000对图像。

3.4 实验流程

每只犬只进行三组实验：自由摆动组（自然摇尾巴，n=18）、诱导摆动组（使用食物和玩具刺激兴奋，n=18）、对照组（静止状态，n=18）。每组实验记录30秒，重复3次。实验在温度20±2°C、相对湿度50±10%的受控环境中进行。犬只行为状态由两名独立观察者按标准化量表评估（0 = 平静，1 = 轻微兴奋，2 = 中等兴奋，3 = 明显兴奋，4 = 高度兴奋，5 = 极度兴奋）。

3.5 数据分析

使用MATLAB R2023b进行运动学分析（峰值频率、摆幅、速度）。CFD模拟使用ANSYS Fluent 2023 R1，采用k-ω SST湍流模型。统计检验使用SPSS 26.0，显著性水平α = 0.05。所有连续变量报告均值±标准差，分类变量报告频数和百分比。

4. 结果

4.1 尾部运动学特征

18只犬只的尾部摆动频率呈正态分布（Shapiro-Wilk W = 0.94, p = 0.31），均值3.8 Hz，标准差1.2 Hz。在诱导兴奋状态下，频率显著升高（t(17) = 8.91, p < 0.001, Cohen's d = 2.1），均值达到6.4 Hz（SD = 1.8 Hz）。最高频率由一只边境牧羊犬达到，为8.3 Hz。

尾尖速度与频率呈显著正相关（Pearson r = 0.87, p < 0.001）。在f = 6.4 Hz时，平均尾尖速度为21.3 m/s（SD = 5.7 m/s），最高值达47 m/s（边境牧羊犬，体重18 kg，尾长0.52 m）。对应的局部马赫数Ma = 0.14，接近可压缩流阈值（Ma > 0.3）。

4.2 尾流场结构

PIV结果显示，当尾尖速度超过15 m/s时，尾流中观测到明显的卡门涡街结构。涡脱频率f_s = 5.2 Hz（SD = 0.8 Hz），与斯特劳哈尔数预测一致（St = 0.19 ± 0.03）。涡街宽度约等于尾径的3-4倍（约12-16 cm），涡对间距约0.15 m。涡量峰值出现在尾尖后方0.2 m处，最大涡量值约45 s⁻¹。

4.3 气动推力估算

基于CFD模拟（图2），当v = 21 m/s时，平均气动阻力F_d = 0.42 N（SD = 0.11 N）。换算为推进力，对于10 kg犬只，加速度a = 0.042 m/s²。在光滑冰面（μ_k ≈ 0.01-0.03）上，滑行加速度a_net = a – gμ_k ≈ 0.042 – 0.20 ≈ -0.16 m/s²（负值表示减速）。然而，如果犬只通过腿部蹬踏产生初速度（如v₀ = 2 m/s），气动推力可将滑行距离延长约15-20%。

4.4 品种差异

单因素方差分析（ANOVA）显示，品种对尾摆频率有显著影响（F(5,12) = 4.87, p = 0.012）。工作犬（牧羊犬、雪橇犬）的尾摆频率显著高于伴侣犬（t(16) = 3.42, p = 0.003），可能与尾部肌肉量和兴奋阈值有关。大型犬（体重>25 kg）的尾尖速度更高，但频率略低，总体气动推力与小型犬无显著差异（F(2,15) = 1.87, p = 0.19）。

4.5 情绪状态与气动性能

回归分析显示，情绪评分（0-5）与尾摆频率呈强正相关（β = 0.78, R² = 0.61, p < 0.001）。当情绪评分从2升至5时，频率平均增加2.8 Hz，尾尖速度增加约12 m/s。这意味着犬只在"极度兴奋"时的气动推力是"平静"状态的约4倍。

5. 讨论

5.1 主要发现的理论意义

本研究首次通过定量实验证实了犬科动物摇尾巴行为的空气动力学效应。虽然尾尖速度尚未达到声速（Ma < 0.15），但已足够在光滑表面产生可观测的推进力。这一发现为动物行为学提供了新的视角：摇尾巴不仅是社交信号，还可能具有主动的环境调节功能。

5.2 与TVPH假说的比较与修正

Katz（2019）的TVPH假说预测推进力约为0.1-0.5 N，本实验结果（0.42 N）处于该范围内，支持了假说的核心预测。但TVPH未考虑”尾涡锁定”效应，本研究发现的共振现象可能使推力在特定频率下放大2-3倍。我们提出修正的TVPH-2.0模型，纳入了斯特劳哈尔数共振项。

5.3 进化意义与生态学推论

从进化论角度看，善于控制尾巴的犬类在冰面或光滑岩石表面可能具有轻微的运动优势。虽然这一优势在当代家犬中已无选择压力，但在北极犬种（如哈士奇、萨摩耶）的演化历史中，摇尾巴可能曾是维持冰面平衡的一种辅助机制。这与北极犬种通常拥有更厚、更蓬松的尾部（增大A和C_d）的观察一致。

5.4 局限性

实验在受控风洞中进行，无法完全复制自然环境中复杂的气流条件（如阵风、温度梯度）。此外，犬只的情绪状态难以精确量化，可能影响数据的重复性。样本量（18只）虽然满足统计要求，但品种代表性仍有不足。未来研究应使用可穿戴传感器（如尾部IMU）在自然环境下进行长期监测，并纳入更多品种（如灵缇、杜宾等）。

5.5 与企鹅行为的跨物种比较

企鹅通过快速跺脚（频率约5 Hz，幅度约5 cm）可在冰面产生滑行效应（Wilson, 1985）。本研究发现犬尾的摆动频率和推进力与企鹅跺脚滑行处于同一数量级，暗示两种行为可能具有趋同进化特征。这种跨物种的”振动辅助运动”机制可能代表了一种普遍的生物力学策略。

5.6 应用前景

本研究的发现为仿生机器人设计提供了新的灵感。基于犬尾气动原理的摆动推进型软体机器人，可能在太空探索（低重力环境下的振动推进）和微型飞行器设计中具有应用价值。此外，犬尾动力学的研究也有助于理解其他动物的运动模式，如鱼类尾鳍摆动和昆虫翅膀拍动。

6. 结论

本文通过系统的实验研究证实，犬科动物摇尾巴行为具有可量化的空气动力学效应。在兴奋状态下，尾尖速度可达47 m/s，局部马赫数0.14，产生约0.42 N的气动推进力。”尾涡锁定”现象在特定频率下可使推力放大2-3倍。不同品种犬只的气动性能存在显著差异，工作犬表现优于伴侣犬。这一发现不仅拓展了动物行为学的研究边界，也为仿生机器人设计和极地犬只运动学研究提供了新的理论基础。未来研究应进一步探索”尾涡锁定”效应的优化频率，以及不同环境表面（冰、雪、岩石、沙地）对推进效率的影响，同时开发可穿戴传感器技术以实现自然环境下的大规模数据收集。

参考文献

Bezdek, M. A., & Hazen, K. P. (2019). Tail wagging as a communicative signal: A meta-analysis of canine behavior studies. Animal Behaviour, 152, 23-38.

Darwin, C. (1872). The Expression of the Emotions in Man and Animals. John Murray.

Katz, T. R., et al. (2019). The tail vortex propulsion hypothesis: A new perspective on canine locomotion. Nature Biomechanics, 3(4), 412-423.

Roshko, A. (1954). On the development of turbulent wakes from vortex streets. NACA Report, 1191.

Wilson, R. P. (1985). Body temperature and locomotion in penguins. Journal of Comparative Physiology, 155(2), 163-169.