smart_hand

第2章:欠驱动与柔顺机构

本章深入探讨灵巧手设计中的两个核心概念:欠驱动机构和柔顺机构。欠驱动通过减少驱动器数量实现低成本、高鲁棒性的自适应抓取;柔顺机构则利用材料变形存储和释放能量,实现被动顺应和安全交互。这两种设计理念的结合,为现代灵巧手带来了革命性的突破,特别是在非结构化环境下的物体操作任务中展现出独特优势。本章将从基础原理出发,逐步深入到工程实现和优化设计,最终通过工业案例和前沿研究展示这些技术的实际应用价值。

2.1 欠驱动原理:自适应抓取与被动顺应

2.1.1 欠驱动的数学定义

欠驱动系统是指独立控制输入数量少于系统自由度数量的机械系统。对于一个具有 $n$ 个自由度的灵巧手,如果只有 $m$ 个独立驱动器($m < n$),则该系统为欠驱动系统。欠驱动度定义为:

\[U = n - m\]

其中 $U > 0$ 表示系统欠驱动。这种设计带来的直接好处是减少了驱动器数量、控制复杂度和系统成本,但代价是失去了对所有自由度的独立控制能力。

2.1.2 自适应抓取机理

欠驱动手指通过机械耦合实现自适应包络抓取。当手指接触物体时,接触力会自动分配到各个关节,实现形状自适应。这个过程可以用虚功原理描述:

\[\tau_{input} = J^T F_{contact}\]

其中 $\tau_{input}$ 是输入扭矩,$J$ 是传动雅可比矩阵,$F_{contact}$ 是接触力向量。传动雅可比矩阵的设计决定了力的传递特性。

2.1.3 传动比设计

欠驱动手指的关键在于合理的传动比设计。对于串联的三关节手指,传动比通常遵循递减原则:

\[r_1 : r_2 : r_3 = 1 : \alpha : \alpha^2\]

其中 $\alpha \in (0.5, 0.8)$ 是传动比系数。这种设计确保指尖具有更大的运动范围和更灵活的适应能力。

     Base                 Tip
      │                    │
    ┌─┴─┐   ┌───┐   ┌───┐ ▼
    │ J1├───┤ J2├───┤ J3│
    └───┘   └───┘   └───┘
     r1      r2      r3

2.1.4 接触顺序与力分配

欠驱动手指的抓取过程可分为三个阶段:

  1. 自由运动阶段:所有关节按预定传动比同步运动
  2. 部分约束阶段:近端关节接触物体后停止,远端关节继续运动
  3. 完全约束阶段:所有关节都接触物体,形成稳定抓取

力分配遵循最小能量原理,可通过拉格朗日乘数法求解:

\(\min \frac{1}{2} F^T W F\) \(s.t. \quad GF = \tau_{ext}\)

其中 $W$ 是权重矩阵,$G$ 是抓取矩阵,$\tau_{ext}$ 是外部负载。

2.2 差动机构设计:行星齿轮与滑轮组

2.2.1 差动传动原理

差动机构允许单一输入驱动多个输出,并根据负载自动分配运动。最典型的是行星齿轮差动机构:

        Sun Gear (Input)
             ○
           ╱ │ ╲
         ○   │   ○  Planet Gears
           ╲ │ ╱
         ┌───┴───┐
         │Carrier│
         └───┬───┘
        Output 1,2

差动方程为: \(\omega_{sun} = \omega_{carrier} + k \cdot \omega_{ring}\)

其中 $k$ 是齿轮比,通过控制太阳轮和齿圈的相对运动实现差动输出。

2.2.2 滑轮-腱绳系统

滑轮组提供了另一种实现欠驱动的方法。N+1腱绳系统使用 $N+1$ 根腱绳控制 $N$ 个关节:

    Motor ──┬── Pulley 1 ── Joint 1
            ├── Pulley 2 ── Joint 2  
            └── Pulley 3 ── Joint 3

腱绳张力与关节扭矩的关系: \(\tau = R \cdot T\)

其中 $R$ 是滑轮半径矩阵,$T$ 是腱绳张力向量。通过设计不同的滑轮半径实现期望的传动比。

2.2.3 耦合矩阵设计

耦合矩阵 $C$ 描述了输入与输出之间的运动关系:

\[q = C \cdot \theta\]

其中 $q$ 是关节角度向量,$\theta$ 是驱动器角度向量。对于欠驱动系统,$C$ 是非方阵(行数大于列数)。

设计原则:

2.2.4 弹簧-离合器机制

结合弹簧和离合器可实现更复杂的欠驱动行为:

\[F_{joint} = \begin{cases} k(x - x_0) & \text{if } x < x_{limit} \\ F_{max} & \text{if } x \geq x_{limit} \end{cases}\]

这种机制允许关节在接触后继续传递有限的力,提高抓取稳定性。

2.3 柔顺机构基础:伪刚体模型与能量方法

2.3.1 柔顺机构的定义与分类

柔顺机构通过材料的弹性变形实现运动和力的传递,分为:

柔顺度矩阵 $C$ 描述了力与变形的关系: \(\delta = C \cdot F\)

2.3.2 伪刚体模型

伪刚体模型将连续体柔顺机构等效为由刚性杆件和扭转弹簧组成的机构:

  Flexible Beam          Pseudo-Rigid Model
  ═══════════      →     ───○───○───
                         Spring Joints

等效弹簧刚度: \(K_{eq} = \frac{EI}{l \cdot \gamma}\)

其中 $E$ 是弹性模量,$I$ 是截面惯性矩,$l$ 是长度,$\gamma$ 是特征半径系数(通常取0.85)。

2.3.3 能量方法分析

使用卡氏定理计算柔顺机构的变形:

\[\delta_i = \frac{\partial U}{\partial F_i}\]

其中 $U$ 是弹性势能: \(U = \frac{1}{2} \int_V \sigma^T \epsilon \, dV\)

对于细长梁结构,可简化为: \(U = \frac{1}{2} \int_0^L \frac{M^2(x)}{EI} \, dx\)

2.3.4 大变形分析

当变形超过线性范围时,需要考虑几何非线性:

\[\kappa(s) = \frac{d\theta}{ds}\]

其中 $\kappa$ 是曲率,$s$ 是弧长参数。椭圆积分法可用于求解大变形问题:

\(x = \int_0^s \cos\theta(\xi) \, d\xi\) \(y = \int_0^s \sin\theta(\xi) \, d\xi\)

2.3.5 双稳态与多稳态设计

利用屈曲失稳实现双稳态机构:

\[F_{critical} = \frac{\pi^2 EI}{(KL)^2}\]

其中 $K$ 是有效长度系数。通过预应力设计可调节稳态位置:

\[U(\theta) = \frac{1}{2}k\theta^2 - F\cdot l\sin\theta\]

极值点满足: \(\frac{dU}{d\theta} = k\theta - F\cdot l\cos\theta = 0\)

2.4 材料选择:硅胶、TPU与形状记忆合金

2.4.1 弹性体材料特性

柔顺机构的性能很大程度上取决于材料选择。主要考虑的材料参数包括:

超弹性本构模型(Neo-Hookean): \(W = \frac{\mu}{2}(I_1 - 3) + \frac{K}{2}(J - 1)^2\)

其中 $\mu$ 是剪切模量,$K$ 是体积模量,$I_1$ 是第一不变量,$J$ 是体积比。

2.4.2 硅胶材料设计

硅胶(如Dragon Skin、Ecoflex)适用于需要大变形的应用:

材料选择准则:

混合配方优化: \(E_{composite} = \phi_1 E_1 + \phi_2 E_2 + \phi_{12} E_{interface}\)

其中 $\phi_i$ 是体积分数,$E_{interface}$ 考虑界面效应。

2.4.3 热塑性聚氨酯(TPU)

TPU提供了刚度和韧性的平衡:

打印参数对性能的影响:

各向异性模型: \(E_{effective} = E_{\parallel} \cos^4\theta + E_{\perp} \sin^4\theta + G_{12} \sin^2(2\theta)\)

2.4.4 形状记忆合金(SMA)

SMA提供主动驱动能力:

相变温度与应力关系(Clausius-Clapeyron): \(\frac{d\sigma}{dT} = -\rho \frac{\Delta H}{T \epsilon_{tr}}\)

其中 $\Delta H$ 是相变潜热,$\epsilon_{tr}$ 是变形应变。

SMA弹簧设计: \(F = \frac{Gd^4}{8D^3n} \cdot \delta\)

其中 $G$ 是剪切模量(马氏体:28 GPa,奥氏体:75 GPa),$d$ 是线径,$D$ 是弹簧直径,$n$ 是有效圈数。

2.4.5 复合材料与梯度材料

功能梯度材料(FGM)实现刚度渐变:

\[E(x) = E_{min} + (E_{max} - E_{min}) \cdot f(x)\]

其中 $f(x)$ 是分布函数(线性、指数或幂函数)。

纤维增强设计: \(E_{composite} = V_f E_f + (1-V_f) E_m\)

其中 $V_f$ 是纤维体积分数,下标 $f$ 和 $m$ 分别表示纤维和基体。

2.5 拓扑优化与参数化设计

2.5.1 拓扑优化基础

拓扑优化通过材料分布最优化实现性能目标:

目标函数(最小化柔顺度): \(\min_{\rho} \quad C = U^T K U\) \(s.t. \quad V(\rho) \leq V_{max}\) \(\quad \quad 0 < \rho_{min} \leq \rho \leq 1\)

其中 $\rho$ 是密度分布,$K$ 是刚度矩阵,$U$ 是位移向量。

SIMP方法(Solid Isotropic Material with Penalization): \(E(\rho) = \rho^p E_0\)

其中 $p$ 是惩罚因子(通常取3)。

2.5.2 柔顺机构拓扑优化

对于柔顺机构,需要同时考虑刚度和柔性:

多目标优化: \(\min \quad J = -\frac{u_{out}}{u_{in}} + w \cdot C\)

其中第一项最大化运动放大率,第二项控制结构刚度。

敏感度分析: \(\frac{\partial J}{\partial \rho_e} = -p\rho_e^{p-1} U_e^T K_0 U_e\)

2.5.3 参数化建模方法

贝塞尔曲线描述柔性铰链轮廓:

\[P(t) = \sum_{i=0}^n \binom{n}{i} (1-t)^{n-i} t^i P_i\]

NURBS曲面用于复杂3D结构: \(S(u,v) = \frac{\sum_{i,j} w_{ij} P_{ij} N_i(u) N_j(v)}{\sum_{i,j} w_{ij} N_i(u) N_j(v)}\)

2.5.4 尺寸优化

梯度优化方法:

\[x^{k+1} = x^k - \alpha \nabla f(x^k)\]

约束处理(增广拉格朗日法): \(L_{\rho}(x,\lambda) = f(x) + \lambda^T g(x) + \frac{\rho}{2} ||g(x)||^2\)

2.5.5 制造约束考虑

最小特征尺寸约束: \(\nabla \cdot (\rho \nabla \rho) \geq \beta\)

对称性约束: \(\rho(x,y) = \rho(-x,y)\)

连通性保证(水平集方法): \(\frac{\partial \phi}{\partial t} + V_n |\nabla \phi| = 0\)

2.6 工业案例:Yale OpenHand项目与Soft Robotics夹爪

2.6.1 Yale OpenHand设计哲学

Yale OpenHand项目pioneered了开源欠驱动手设计,其核心创新包括:

SDM手指(Shape Deposition Manufacturing)

关键设计参数:

Model T系列: 采用四连杆差动机构: \(\theta_{distal} = \alpha \cdot \theta_{proximal} + \beta \cdot \theta_{differential}\)

其中 $\alpha = 1.5$,$\beta = 0.8$ 实现指尖优先接触。

2.6.2 Soft Robotics商业化案例

Soft Robotics的mGrip系统展示了柔顺机构的工业应用:

技术特点

控制策略: 压力-变形模型: \(F_{grip} = A \cdot P \cdot (1 - e^{-\kappa \cdot \delta})\)

其中 $A$ 是有效面积,$P$ 是压力,$\kappa$ 是材料常数,$\delta$ 是变形量。

2.6.3 性能对比分析

指标 Yale OpenHand Soft Robotics 传统刚性夹爪
自适应能力 极高
抓取力 15-20N 5-40N 50-200N
位置精度 ±5mm ±10mm ±0.1mm
成本
维护需求 极低

2.6.4 应用场景分析

食品行业

物流分拣

协作机器人

2.7 高级话题:基于机器学习的柔顺机构逆向设计

2.7.1 逆向设计问题定义

传统的柔顺机构设计是正向过程:给定结构参数预测性能。逆向设计则相反:给定期望性能,自动生成满足要求的结构。这是一个典型的逆问题,存在以下挑战:

问题形式化: \(\mathcal{D}^* = \arg\min_{\mathcal{D}} \mathcal{L}(f(\mathcal{D}), \mathcal{P}_{target})\)

其中 $\mathcal{D}$ 是设计参数,$f$ 是正向仿真,$\mathcal{P}_{target}$ 是目标性能,$\mathcal{L}$ 是损失函数。

2.7.2 深度生成模型方法

变分自编码器(VAE)架构

编码器学习设计空间的潜在表示: \(q_\phi(z|x) = \mathcal{N}(\mu_\phi(x), \sigma_\phi^2(x))\)

解码器生成设计参数: \(p_\theta(x|z) = \mathcal{N}(\mu_\theta(z), \sigma_\theta^2(z))\)

损失函数: \(\mathcal{L} = \mathbb{E}_{q_\phi(z|x)}[\log p_\theta(x|z)] - D_{KL}(q_\phi(z|x)||p(z))\)

条件生成对抗网络(cGAN)

生成器损失: \(\mathcal{L}_G = \mathbb{E}_{z,c}[\log(1-D(G(z|c),c))] + \lambda||\mathcal{P}(G(z|c)) - c||^2\)

其中 $c$ 是性能条件,第二项确保生成设计满足性能要求。

2.7.3 强化学习优化策略

将设计过程建模为马尔可夫决策过程(MDP):

策略梯度方法(PPO): \(\mathcal{L}^{CLIP}(\theta) = \mathbb{E}_t[\min(r_t(\theta)\hat{A}_t, \text{clip}(r_t(\theta), 1-\epsilon, 1+\epsilon)\hat{A}_t)]\)

其中 $r_t(\theta) = \frac{\pi_\theta(a_t s_t)}{\pi_{\theta_{old}}(a_t s_t)}$,$\hat{A}_t$ 是优势函数估计。

2.7.4 物理信息神经网络(PINN)

结合物理约束的神经网络训练:

\[\mathcal{L} = \mathcal{L}_{data} + \lambda_1 \mathcal{L}_{PDE} + \lambda_2 \mathcal{L}_{BC}\]

PDE约束(弹性力学): \(\mathcal{L}_{PDE} = ||\nabla \cdot \sigma + f||^2\)

边界条件约束: \(\mathcal{L}_{BC} = ||u|_{\Gamma_D} - u_0||^2 + ||\sigma \cdot n|_{\Gamma_N} - t_0||^2\)

2.7.5 实例:软体手指形状优化

问题设定: 设计软体手指形状以实现特定抓取力分布。

网络架构

Input: Force Distribution (100D)
    ↓
Encoder: FC(256) → FC(128) → FC(64)
    ↓
Latent Space (32D)
    ↓
Decoder: FC(64) → FC(128) → FC(256)
    ↓
Output: Finger Geometry (Bezier Control Points, 20D)

训练数据生成

  1. 参数采样:Latin Hypercube Sampling
  2. FEA仿真:Abaqus/COMSOL
  3. 数据增强:几何变换、噪声注入

结果验证

2.7.6 未来发展方向

多尺度设计: 同时优化宏观结构和微观材料分布: \(\mathcal{D} = \{\mathcal{D}_{macro}, \mathcal{D}_{meso}, \mathcal{D}_{micro}\}\)

自适应材料编程: 学习时变刚度分布: \(E(x,t) = NN_\theta(x, t, \tau)\)

其中 $\tau$ 是任务参数。

仿真-实物迁移: 域适应技术缩小sim-to-real gap: \(\mathcal{L}_{DA} = \mathcal{L}_{task} + \lambda \cdot \mathcal{L}_{MMD}(P_{sim}, P_{real})\)

本章小结

本章系统介绍了欠驱动与柔顺机构在灵巧手设计中的应用。核心要点包括:

  1. 欠驱动原理:通过减少驱动器实现自适应抓取,关键在于合理的传动比设计和力分配策略
  2. 差动机构:行星齿轮和滑轮组提供了灵活的运动耦合方案
  3. 柔顺机构:利用材料弹性实现运动传递,伪刚体模型简化了分析过程
  4. 材料选择:从超软硅胶到形状记忆合金,不同材料适合不同应用场景
  5. 优化设计:拓扑优化和参数化方法实现性能最优
  6. 工业实践:Yale OpenHand和Soft Robotics展示了技术的实用性
  7. 智能设计:机器学习方法开启了逆向设计新范式

关键公式回顾:

练习题

基础题

习题 2.1 一个三关节欠驱动手指,每个关节的运动范围为0-90度。如果采用传动比 $r_1:r_2:r_3 = 1:0.7:0.49$,当驱动器转动120度时,计算各关节的理论转角(假设无负载)。

提示:考虑传动比的累积效应

参考答案 根据传动比关系: - 第一关节:$\theta_1 = 120° \times 1 = 120°$,但受限于运动范围,实际 $\theta_1 = 90°$ - 剩余30°传递到下一级 - 第二关节:$\theta_2 = 30° \times 0.7 + 90° \times 0.7 = 84°$ - 第三关节:$\theta_3 = 84° \times 0.7 = 58.8°$ 实际角度:$\theta_1 = 90°$,$\theta_2 = 84°$,$\theta_3 = 58.8°$

习题 2.2 设计一个柔性铰链,要求在5N的力作用下产生15度的转角。如果使用TPU材料(E=50MPa),铰链长度为10mm,计算所需的截面惯性矩。

提示:使用伪刚体模型,$K_{eq} = \frac{EI}{l \cdot \gamma}$,其中$\gamma = 0.85$

参考答案 扭转刚度要求:$K = \frac{M}{\theta} = \frac{5N \times 0.01m}{15° \times \pi/180} = 0.191 N·m/rad$ 根据伪刚体模型: $$I = \frac{K \times l \times \gamma}{E} = \frac{0.191 \times 0.01 \times 0.85}{50 \times 10^6} = 3.25 \times 10^{-11} m^4$$ 对于矩形截面(宽度w=5mm): $$I = \frac{wh^3}{12}$$ $$h = \sqrt[3]{\frac{12I}{w}} = 0.55mm$$

习题 2.3 一个滑轮-腱绳系统有3个滑轮,半径分别为10mm、15mm、20mm。如果腱绳张力为10N,计算各关节产生的扭矩。

提示:扭矩=半径×张力

参考答案 各关节扭矩: - 关节1:$\tau_1 = 0.01m \times 10N = 0.1 N·m$ - 关节2:$\tau_2 = 0.015m \times 10N = 0.15 N·m$ - 关节3:$\tau_3 = 0.02m \times 10N = 0.2 N·m$ 扭矩比:$\tau_1 : \tau_2 : \tau_3 = 1 : 1.5 : 2$

挑战题

习题 2.4 设计一个欠驱动两指夹爪,要求能够自适应抓取直径20-100mm的圆柱形物体。请确定: a) 手指长度和关节数量 b) 传动比设计 c) 驱动器扭矩要求(抓取力目标20N)

提示:考虑几何约束和力传递效率

参考答案 a) 手指设计: - 采用3关节设计,每段长度40mm,总长120mm - 可覆盖20-100mm直径范围 b) 传动比: - 采用递减传动比:1:0.6:0.36 - 确保远端关节有足够灵活性 c) 驱动器扭矩: 最不利情况(100mm直径): - 力臂约50mm - 所需扭矩:$\tau = 20N \times 0.05m = 1 N·m$ - 考虑传动效率(0.7):$\tau_{motor} = 1/0.7 = 1.43 N·m$

习题 2.5 使用拓扑优化设计一个柔顺放大机构,输入位移1mm,输出位移5mm。给定设计域100mm×50mm,体积约束30%。写出优化问题的数学形式并讨论求解策略。

提示:考虑几何增益和结构刚度的平衡

参考答案 优化问题形式: $$\max_{\rho} \quad GA = \frac{u_{out}}{u_{in}}$$ $$s.t. \quad K(\rho)U = F$$ $$\quad \quad \sum \rho_i v_i \leq 0.3 V_{total}$$ $$\quad \quad 0.001 \leq \rho_i \leq 1$$ 求解策略: 1. 使用SIMP方法进行材料插值 2. 敏感度分析采用伴随法 3. 优化算法:MMA(Method of Moving Asymptotes) 4. 过滤技术避免棋盘格现象 5. 投影方法确保0-1设计

习题 2.6 分析一个硅胶软体手指在大变形下的行为。手指长度80mm,截面10mm×5mm,材料Neo-Hookean模型(μ=0.1MPa)。当指尖受到5N横向力时,估算指尖位移。

提示:需要考虑几何非线性

参考答案 大变形悬臂梁分析: 1. 线性估算(小变形): $$\delta = \frac{FL^3}{3EI} = \frac{5 \times 0.08^3}{3 \times 0.3 \times 10^6 \times 4.17 \times 10^{-10}} = 6.8mm$$ 2. 考虑大变形修正: 使用椭圆积分方法,实际位移约为线性预测的1.5-2倍 $$\delta_{actual} \approx 10-14mm$$ 3. 几何非线性效应: - 轴向缩短 - 应力刚化 - 接触自锁可能

习题 2.7(开放性思考题) 如何结合欠驱动和柔顺机构的优势,设计一个能够操作易碎物品(如鸡蛋)的机械手?讨论设计要点、材料选择和控制策略。

参考思路 设计要点: 1. **结构设计**: - 欠驱动实现形状自适应 - 柔顺材料提供被动顺应 - 多指设计增加稳定性 2. **材料选择**: - 接触面:超软硅胶(Shore 00-20) - 支撑结构:TPU或半刚性材料 - 可考虑充气结构 3. **控制策略**: - 力控制而非位置控制 - 阻抗控制实现柔顺交互 - 触觉反馈检测接触状态 4. **安全机制**: - 机械限位防止过度挤压 - 压力释放阀(气动) - 被动柔顺吸收冲击

习题 2.8 设计一个基于形状记忆合金的自适应刚度手指。要求在常温下柔软(便于被动适应),加热后变硬(提供抓取力)。给出设计方案和性能估算。

参考答案 设计方案: 1. **结构设计**: - SMA丝嵌入硅胶基体 - 螺旋布置增加行程 - 分段控制实现局部刚度调节 2. **材料参数**: - SMA丝径:0.5mm - 相变温度:45-55°C - 硅胶基体:Shore A 30 3. **性能估算**: - 刚度变化比:10-20倍 - 响应时间:加热2-3s,冷却5-10s - 最大收缩应变:4-6% - 驱动应力:200-300 MPa 4. **控制方法**: - PWM控制加热电流 - 温度反馈维持相变 - 预测控制补偿滞后

常见陷阱与错误

设计陷阱

  1. 过度依赖仿真
    • 错误:完全相信FEA结果,忽略材料非线性和制造偏差
    • 正确:仿真指导设计,实验验证性能
  2. 忽视疲劳失效
    • 错误:只考虑静态强度
    • 正确:进行疲劳分析,特别是柔性铰链
  3. 传动比设计不当
    • 错误:传动比过大导致远端关节无法充分运动
    • 正确:递减传动比,保证运动协调

材料选择错误

  1. 硬度选择不当
    • 错误:一味追求柔软,导致承载能力不足
    • 正确:根据应用场景平衡柔顺性和刚度
  2. 忽略环境因素
    • 错误:不考虑温度、湿度对材料性能的影响
    • 正确:选择环境稳定性好的材料或添加保护
  3. 材料相容性
    • 错误:不同材料直接粘接,界面失效
    • 正确:使用合适的粘接剂或机械连接

制造问题

  1. 公差累积
    • 错误:忽略制造公差对机构性能的影响
    • 正确:鲁棒性设计,关键尺寸严格控制
  2. 支撑结构设计
    • 错误:3D打印时支撑不当导致变形
    • 正确:优化打印方向,设计易去除支撑

控制误区

  1. 刚度补偿不足
    • 错误:不考虑柔顺变形对精度的影响
    • 正确:建立变形模型,实时补偿
  2. 响应速度期望过高
    • 错误:期望软体机构有刚性机构的响应速度
    • 正确:合理设定性能指标,优化动态特性

最佳实践检查清单

设计阶段

材料选择

仿真验证

原型制作

性能测试

优化迭代

系统集成