机械结构设计是轮足机械臂机器人性能的决定性因素。本章深入探讨结构刚度与重量的优化权衡、振动特性对控制系统的影响、万向节锁的机械规避策略、线缆布局的动力学考虑，以及先进材料的工程应用。这些看似纯机械的设计决策，直接影响着算法的实现难度、控制带宽的上限和系统的鲁棒性。对于AI工程师而言，理解这些硬件约束是设计高性能控制算法的前提。
机器人的结构刚度直接决定了其定位精度和动态响应特性。刚度不足会导致：
而过度追求刚度则导致：
关键设计指标是比刚度（specific stiffness）：
\[E/\rho\]其中 $E$ 是材料的弹性模量，$\rho$ 是密度。
拓扑优化
基于有限元分析（FEA）的拓扑优化可以在给定设计域内找到最优材料分布。优化目标函数：
\[\min_{\mathbf{x}} \quad c(\mathbf{x}) = \mathbf{U}^T\mathbf{K}\mathbf{U}\]约束条件： \(V(\mathbf{x}) \leq V_0\) \(\mathbf{K}\mathbf{U} = \mathbf{F}\)
其中 $c$ 是柔度（刚度的倒数），$\mathbf{U}$ 是位移向量，$\mathbf{K}$ 是刚度矩阵，$V$ 是体积，$\mathbf{F}$ 是外力向量。
截面优化
对于梁式结构，惯性矩 $I$ 的优化至关重要：
\[I = \int_{A} y^2 \, dA\]空心矩形截面相比实心截面，在相同重量下可提供更高的抗弯刚度。优化准则：
\[\frac{I}{A} \propto \frac{h^2}{t}\]其中 $h$ 是截面高度，$t$ 是壁厚。
机器人关节的等效刚度包含多个串联环节：
\[\frac{1}{K_{eq}} = \frac{1}{K_{motor}} + \frac{1}{K_{gear}} + \frac{1}{K_{bearing}} + \frac{1}{K_{structure}}\]典型数值（以7自由度机械臂为例）：
最弱环节决定了系统整体刚度，通常是减速器和结构。
实际刚度测量方法：
刚度矩阵标定：
\[\mathbf{F} = \mathbf{K}_c \Delta \mathbf{X}\]其中 $\mathbf{K}_c$ 是笛卡尔刚度矩阵（6×6），可通过雅可比变换从关节刚度得到：
\[\mathbf{K}_c = \mathbf{J}^{-T} \mathbf{K}_j \mathbf{J}^{-1}\]机器人结构的振动特性由其质量矩阵 $\mathbf{M}$ 和刚度矩阵 $\mathbf{K}$ 决定。特征值问题：
\[(\mathbf{K} - \omega_i^2 \mathbf{M})\boldsymbol{\phi}_i = 0\]其中 $\omega_i$ 是第 $i$ 阶固有频率，$\boldsymbol{\phi}_i$ 是对应振型。
对于典型6轴机械臂：
控制系统的带宽必须远低于最低固有频率，经验准则：
\[f_{control} < \frac{f_{natural}}{3}\]这意味着如果第一阶固有频率是30 Hz，控制带宽应限制在10 Hz以下。
相位裕度考虑
在固有频率附近，系统相位快速变化，容易导致不稳定。Nyquist稳定性要求：
\[PM > 30° \quad \text{at} \quad f = f_{natural}\]被动阻尼
添加粘弹性材料或调谐质量阻尼器（TMD）：
\(m_{TMD} = \mu \cdot m_{structure}\) \(f_{TMD} = \frac{f_{target}}{1 + \mu}\)
其中 $\mu$ 是质量比，通常取0.01-0.05。
主动振动控制
输入整形（Input Shaping）技术，通过预先设计的脉冲序列消除残余振动：
\[u(t) = \sum_{i=1}^{n} A_i \delta(t - t_i)\]零振动（ZV）整形器参数： \(A_1 = \frac{1}{1+K}, \quad A_2 = \frac{K}{1+K}\) \(t_1 = 0, \quad t_2 = \frac{T_d}{2}\)
其中 $K = e^{-\zeta\pi/\sqrt{1-\zeta^2}}$，$\zeta$ 是阻尼比，$T_d$ 是阻尼周期。
现代设计方法将结构和控制器同时优化：
\[\min_{p, K} J = \int_0^{\infty} (\mathbf{x}^T\mathbf{Q}\mathbf{x} + \mathbf{u}^T\mathbf{R}\mathbf{u}) dt\]约束：
这种协同设计可以在不增加重量的前提下，通过优化质量分布提高可控性。
当使用欧拉角表示姿态时，某些构型下会失去一个自由度。以ZYX欧拉角为例，当 $\theta_y = \pm 90°$ 时：
\[\mathbf{R} = \mathbf{R}_z(\phi) \mathbf{R}_y(\pm 90°) \mathbf{R}_x(\psi)\]此时 $\phi$ 和 $\psi$ 的效果相同，系统退化为2自由度。
雅可比矩阵的奇异性：
\[\det(\mathbf{J}) = \cos(\theta_y) \rightarrow 0 \quad \text{when} \quad \theta_y \rightarrow \pm 90°\]冗余自由度设计
添加第7个自由度（如7轴机械臂），利用零空间避开奇异构型：
\[\dot{\mathbf{q}} = \mathbf{J}^+ \dot{\mathbf{x}} + (\mathbf{I} - \mathbf{J}^+ \mathbf{J})\mathbf{z}\]其中 $\mathbf{z}$ 是零空间投影向量，可用于优化次要目标。
关节限位设计
机械限位避免进入奇异区域：
\[\theta_{limit} = \theta_{singular} - \Delta\theta_{safety}\]安全裕度 $\Delta\theta_{safety}$ 通常取10-15°。
偏置设计
通过非正交轴设计避免万向节锁：
传统设计： 优化设计：
Z Z'
| /
| /
Y---X Y'---X'
正交系 偏置5-10°
四元数表示避免了万向节锁，但需要正确的插值方法。球面线性插值（SLERP）：
\[q(t) = \frac{\sin((1-t)\Omega)}{\sin(\Omega)}q_0 + \frac{\sin(t\Omega)}{\sin(\Omega)}q_1\]其中 $\Omega = \arccos(q_0 \cdot q_1)$。
四元数误差用于姿态控制：
\[\mathbf{e}_o = \mathbf{q}_d \otimes \mathbf{q}_c^{-1}\]转换为轴角表示用于控制律：
\[\boldsymbol{\omega}_{cmd} = 2 \cdot \text{sign}(e_0) \cdot [e_1, e_2, e_3]^T\]现代协作机器人（如Universal Robots）采用的策略：
奇异性度量指标：
\[\mu = \frac{\sigma_{min}}{\sigma_{max}}\]其中 $\sigma$ 是雅可比矩阵的奇异值。当 $\mu < 0.01$ 时触发奇异性保护。
线缆在关节运动时产生的反作用力矩：
\[\tau_{cable} = k_{cable} \cdot \theta + c_{cable} \cdot \dot{\theta}\]其中 $k_{cable}$ 是线缆刚度，$c_{cable}$ 是阻尼系数。
多芯线缆的等效刚度：
\[k_{total} = \sum_{i=1}^{n} k_i \cdot r_i^2\]其中 $r_i$ 是第 $i$ 根线缆到旋转中心的距离。
中空轴设计
将线缆穿过关节中心，最小化力臂：
┌─────┐
│ │ 外壳
│ ┌─┐ │
│ │o│ │ 中空轴
│ └─┘ │ 线缆通道
└─────┘
优势：
拖链系统设计
线缆长度计算：
\[L_{cable} = L_0 + \theta_{max} \cdot r + L_{safety}\]疲劳寿命估算（基于S-N曲线）：
\[N = A \cdot \sigma^{-b}\]其中 $N$ 是循环次数，$\sigma$ 是应力幅值，$A$、$b$ 是材料常数。
前馈补偿
基于线缆模型的力矩补偿：
\[\tau_{motor} = \tau_{desired} + \hat{\tau}_{cable}(\theta, \dot{\theta})\]自适应估计
使用递归最小二乘（RLS）在线估计线缆参数：
\[\hat{\mathbf{p}}_{k+1} = \hat{\mathbf{p}}_k + \mathbf{K}_k(\tau_{measured} - \boldsymbol{\phi}_k^T \hat{\mathbf{p}}_k)\]其中 $\mathbf{p} = [k_{cable}, c_{cable}]^T$，$\boldsymbol{\phi} = [\theta, \dot{\theta}]^T$。
完全无线化面临的技术挑战：
部分无线化方案：
| 特性 | 碳纤维复合材料 (T700) | 铝合金 (7075-T6) | 钛合金 (Ti-6Al-4V) |
|---|---|---|---|
| 密度 (g/cm³) | 1.6 | 2.8 | 4.4 |
| 抗拉强度 (MPa) | 2550 | 570 | 950 |
| 弹性模量 (GPa) | 135 | 72 | 114 |
| 比强度 | 1594 | 204 | 216 |
| 比刚度 | 84 | 26 | 26 |
| 疲劳性能 | 优秀 | 一般 | 良好 |
| 成本 (相对) | 10-20× | 1× | 30-50× |
层合板设计
根据主应力方向优化纤维铺层角度：
\[[0°/±45°/90°]_s\]准各向同性铺层，平衡各方向性能。
对于扭转载荷为主的部件：
\[[±45°]_n\]失效准则
Tsai-Wu失效准则：
\[F_1\sigma_1 + F_2\sigma_2 + F_{11}\sigma_1^2 + F_{22}\sigma_2^2 + 2F_{12}\sigma_1\sigma_2 + F_{66}\tau_{12}^2 < 1\]预浸料成型
孔隙率控制在<1%以保证强度。
连接设计
碳纤维与金属连接的关键问题：
总拥有成本（TCO）模型：
\[TCO = C_{material} + C_{manufacturing} + C_{energy} \cdot T_{operation}\]碳纤维的高初始成本可通过以下方式回收：
临界点分析：当机器人臂长>1.5m或载荷>20kg时，碳纤维通常更具成本效益。
宇树B2是国产高性能机械臂的代表，其轻量化设计体现了先进的工程理念。
动力学参数优化
通过遗传算法优化连杆参数：
\[\min J = w_1 \cdot m_{total} + w_2 \cdot \frac{1}{f_1} + w_3 \cdot \sigma_{max}\]约束条件：
模块化设计理念
标准接口定义：
这种模块化允许快速配置不同自由度和负载能力的机器人。
密度法（SIMP）优化问题：
\[\min_{\rho} \quad c = \sum_{e=1}^{N} \rho_e^p \mathbf{u}_e^T \mathbf{k}_e \mathbf{u}_e\]约束： \(\sum_{e=1}^{N} v_e \rho_e \leq V^*\) \(\mathbf{K}\mathbf{U} = \mathbf{F}\) \(0 < \rho_{min} \leq \rho_e \leq 1\)
其中 $\rho_e$ 是单元密度，$p$ 是惩罚因子（通常取3）。
金属3D打印（SLM/EBM）允许：
设计约束：
骨小梁启发的多孔结构：
\[E^* = E_s \cdot \left(\frac{\rho^*}{\rho_s}\right)^n\]其中 $E^*$ 是多孔结构模量，$n≈2$ 对于开孔泡沫。
这种结构可实现：
机械结构设计对机器人性能的影响是全方位的：
关键设计准则：
这些硬件特性直接影响算法设计：
练习3.1 一个3自由度机械臂，各关节刚度分别为 $K_1 = 10^5$ Nm/rad、$K_2 = 5×10^4$ Nm/rad、$K_3 = 2×10^4$ Nm/rad。如果末端施加100N的力，臂长均为0.5m，计算最坏情况下的末端位置偏差。
提示：考虑机械臂完全伸展的构型。
练习3.2 某机器人结构的第一阶固有频率为25 Hz，阻尼比 $\zeta = 0.05$。设计一个ZV输入整形器，计算脉冲幅值和时间间隔。
提示：使用ZV整形器公式。
练习3.3 碳纤维管和铝合金管，外径50mm，壁厚2mm，长度1m。计算两者的弯曲刚度比和重量比。 （$E_{CF} = 135$ GPa，$\rho_{CF} = 1.6$ g/cm³；$E_{Al} = 70$ GPa，$\rho_{Al} = 2.8$ g/cm³）
提示：使用薄壁圆管的惯性矩公式。
练习3.4 设计一个机械臂连杆的拓扑优化问题。连杆长500mm，两端承受弯矩±100 Nm，中间承受50N垂直力。目标是最小化柔度，体积约束为原始体积的40%。写出优化问题的数学形式并讨论求解策略。
提示：考虑SIMP方法和敏感度分析。
练习3.5 一个7自由度机械臂在构型 $\mathbf{q} = [0, 45°, 0, 90°, 0, 45°, 0]$ 时，计算其可操作度指标。如果第5关节接近±90°会发生腕部奇异，设计一个基于零空间运动的奇异性规避算法。
提示：使用雅可比矩阵的条件数和零空间投影。
练习3.6 机器人关节使用扁平柔性电缆，宽度30mm，包含20根导线。电缆在±180°范围内反复弯曲，弯曲半径50mm。估算电缆的疲劳寿命，并提出改进设计。
提示：使用Coffin-Manson关系和应变计算。
练习3.7 使用四元数设计一个机械臂的姿态控制器。当前姿态 $\mathbf{q}_c = [0.7071, 0, 0.7071, 0]$，目标姿态 $\mathbf{q}_d = [1, 0, 0, 0]$。计算姿态误差和控制力矩。
提示：使用四元数乘法和误差四元数。
❌ 错误：只关注静态刚度，忽略动态特性
✅ 正确：同时考虑静态刚度和固有频率
静态分析 → 确保变形<0.1mm
模态分析 → 确保f₁>30Hz
谐响应分析 → 识别共振峰
❌ 错误：盲目使用碳纤维
✅ 正确：基于具体需求选择
臂长<1m、载荷<10kg → 铝合金
臂长>1.5m、载荷>20kg → 碳纤维
需要阻尼 → 复合材料夹层
❌ 错误：线缆长度刚好够用
✅ 正确：预留30%余量，考虑疲劳
L_cable = L_min × 1.3 + π × r × θ_max/180
❌ 错误：用高增益控制器抑制振动
✅ 正确：陷波滤波器 + 输入整形
# 陷波滤波器设计
notch_freq = 25 # Hz
Q = 10
# 在控制环路中滤除共振频率
❌ 错误：检测到奇异性时紧急停止
✅ 正确：提前规避 + 平滑过渡
if manipulability < 0.05:
activate_null_space_motion()
reduce_speed_gradually()
这份检查清单可作为设计评审的标准模板，确保机械结构设计的完整性和可靠性。