1. 软体执行器SSA的建模基础与材料特性
软体执行器(Soft Structural Actuator, SSA)作为新一代柔性驱动装置，其核心优势在于能够模拟生物肌肉的柔顺性和适应性。与传统刚性执行器相比，SSA采用超弹性材料构建，通过内部气压变化实现大变形运动，特别适用于人机交互场景。本章将深入解析SSA的建模理论基础和关键材料特性。
1.1 超弹性材料本构关系
热塑性聚氨酯(TPU)作为SSA的核心材料，展现出典型的超弹性力学行为。在单轴拉伸测试中，NinjaTek公司的TPU材料可承受高达600%的应变而不发生塑性变形。这种非线性应力-应变关系需要通过应变能函数W来描述，其一般形式为：
W = W(I₁, I₂, I₃)其中I₁、I₂、I₃为应变不变量，通过主伸长比λ计算得到：
I₁ = λ₁² + λ₂² + λ₃² I₂ = λ₁²λ₂² + λ₂²λ₃² + λ₃²λ₁² I₃ = λ₁²λ₂²λ₃²注意：主伸长比λ定义为变形后长度L与原始长度L₀的比值(λ=L/L₀)，在小应变条件下可近似表示为λ≈1+ε，其中ε为工程应变。
1.2 本构模型对比与选择
通过ASTM D638标准试样的拉伸实验(图1)，我们对比了五种主流超弹性模型对TPU材料的拟合效果：
| 模型名称 | 应变能函数形式 | 适用场景 |
|---|---|---|
| Neo-Hookean | W=C₁₀(I₁-3) | 小变形近似 |
| Mooney-Rivlin | W=C₁₀(I₁-3)+C₀₁(I₂-3) | 中等变形 |
| Yeoh | W=ΣCᵢ(I₁-3)ⁱ (i=1-3) | 大变形压缩 |
| Ogden | W=Σ(μₐ/αₐ)(λ₁^αₐ+λ₂^αₐ+λ₃^αₐ-3) | 多模式变形 |
| Arruda-Boyce | 基于统计力学模型 | 分子链网络变形 |
实验数据表明，五参数Mooney-Rivlin模型对TPU85和TPU95的拟合优度(R²)分别达到0.998和0.997，其完整表达式为：
W = C₁₀(I₁-3) + C₀₁(I₂-3) + C₁₁(I₁-3)(I₂-3) + C₂₀(I₁-3)² + C₀₂(I₂-3)²通过ANSYS Workbench的曲线拟合功能，获得的材料参数如表1所示：
表1 TPU85与TPU95的Mooney-Rivlin参数(单位：MPa)
| 参数 | TPU85值 | TPU95值 |
|---|---|---|
| C₁₀ | -3.1992 | -28.763 |
| C₀₁ | 6.977 | 42.995 |
| C₂₀ | 0.0281 | 0.10499 |
| C₁₁ | -0.074972 | -6.6676 |
| C₀₂ | 0.92155 | 9.138 |
| D₁ | 0 | 0 |
实操技巧：材料测试时应确保试样制备工艺与最终产品一致，包括打印温度(TPU建议210-220℃)、层厚(0.2mm)和填充率(100%)，以保证数据可比性。
2. SSA几何建模与运动学分析
2.1 折叠结构运动原理
SSA的创新性在于其独特的折叠式波纹管设计(图2)，通过几何参数化建模可预测三种基本运动模式：
- 轴向伸展：折叠展开导致整体长度增加
- 轴向收缩：折叠压紧使长度减小
- 弯曲运动：非对称折叠产生曲率变化
关键几何参数包括：
- 折叠长度S
- 折叠角θ
- 波纹管壁厚t
- 约束层厚度tr
- 折叠数量N
2.2 位移建模与验证
2.2.1 伸展位移模型
单个折叠单元的伸展位移δ_ext_single可表示为：
δ_ext_single = 2S(1 - sinθ)对于N个折叠的SSA，总伸展位移为：
δ_ext = N × δ_ext_single = 2NS(1 - sinθ)实验验证显示，当θ=30°、S=15mm时，理论计算与实测值的平均误差<5%。值得注意的是，随着θ增大，伸展位移非线性减小，当θ→90°时，sinθ→1，导致δ_ext→0。
2.2.2 收缩位移模型
收缩模式下，总位移δ_cont计算公式为：
δ_cont = 2N × S × sinθ与伸展模式相反，增大θ会提升收缩位移，但实际测试发现当θ>45°时，折叠结构会出现自锁现象，导致有效位移降低。
2.2.3 弯曲运动建模
弯曲曲率半径ρ由下式确定：
ρ = L(r + tr) / (N × δ_ext_single)其中L为约束侧长度，r为波纹管内径。通过控制两侧折叠的非对称展开，可实现精确的曲率控制。
设计经验：实际应用中建议将θ控制在20°-40°范围内，既可保证足够的位移输出，又避免结构失稳。折叠数量N通常取6-12个，过多会导致响应速度下降。
3. 静态力学模型与参数影响
3.1 静态力平衡方程
在准静态条件下，SSA的输出力F_out由气压作用力与轴向刚度力平衡决定：
F_out = Σ(F₁ + F₂y - F₃y) - F_stiff其中：
- F₁：端盖作用力 (F₁ = P × A₁)
- F₂y：外壁轴向分力
- F₃y：内壁轴向分力
- F_stiff：轴向刚度力
3.2 轴向刚度特性
通过多项式拟合获得的轴向刚度力模型：
F_stiff = 4.1481×10⁻⁴y³ + 1.2865×10⁻²y² + 2.0789y - 0.2246对应刚度为：
K_stiff = dF_stiff/dy = 1.24443×10⁻³y² + 2.5730×10⁻²y + 2.0789关键设计参数的影响规律：
- 壁厚增加20%→ 刚度提升35-50%
- 约束层增加1层→ 刚度提升约25%
- 材料硬度从85A到95A→ 刚度提升3-4倍
3.3 压力-位移特性
静态测试数据表明(图3)：
- 低压区(0-50kPa)：位移随压力线性增加
- 中压区(50-150kPa)：出现明显非线性
- 高压区(>150kPa)：位移趋于饱和
典型性能指标：
- 最大伸展位移：57.5mm (L13型号)
- 阻塞力：85N (TPU85材料)
- 能效比：1.2J/mm
4. 动态响应与控制系统
4.1 动态微分方程
基于牛顿第二定律建立的SSA动力学模型：
Mÿ + bẏ + (1.24443×10⁻³y³ + 2.5730×10⁻²y² + 2.0789y) = (A₁+A₂-A₃)P其中：
- M：负载质量
- b：阻尼系数(TPU85约0.8 N·s/m)
- P：控制压力
4.2 PID控制器设计
采用根轨迹法整定的PID控制律：
u(t) = K_p e(t) + K_i ∫e(τ)dτ + K_d de(t)/dt优化后的参数：
- K_p = 12.5
- K_i = 0.8
- K_d = 0.25
性能指标：
- 上升时间：1.0s
- 调节时间：2.0s
- 超调量：<5%
4.3 动态性能测试
4.3.1 频率响应
带宽测试结果(图4)：
- -3dB带宽：0.65Hz
- 最大响应频率：5Hz(位移衰减至3mm)
- 极限频率：20Hz(阀件限制)
4.3.2 阶跃响应
10mm阶跃输入时：
- 上升时间：0.8s
- 稳定时间：1.0s
- 稳态误差：0μm
- 抗扰恢复时间：0.2s(4.5kg负载)
4.3.3 轨迹跟踪
- 斜坡跟踪(1.5mm/s)：误差<50μm
- 正弦跟踪(0.5Hz)：相位滞后<20°
- 复杂轨迹RMSE：≤51.9μm
调试心得：比例阀相比电磁阀虽然响应稍慢(上升时间增加15%)，但能显著降低压力波动，使位移波动幅度从±100μm降至±25μm。对于微米级定位应用，建议采用伺服阀方案。
5. 工程应用与性能优化
5.1 穿戴辅助设备应用
SSA在穿戴式移动辅助设备(WMAD)中的独特优势：
- 结构兼容性：中空设计便于肢体穿套
- 安全交互：本质柔顺避免机械伤害
- 功率密度：单位重量出力比达15N/g
临床测试表明，在膝关节辅助场景下：
- 步态周期跟踪误差<3°
- 能耗降低22%相比传统电机驱动
- 用户舒适度评分提升35%
5.2 参数优化策略
根据应用需求的设计导向：
| 性能需求 | 优化方向 | 典型参数组合 |
|---|---|---|
| 大位移输出 | 减小θ，增加S，降低硬度 | θ=20°, S=20mm, TPU85 |
| 高力输出 | 增加壁厚，多层约束 | t=2mm, 3约束层, TPU95 |
| 快速响应 | 减少N，优化气路 | N=6, 高速阀, 短管路 |
| 高精度控制 | 提高刚度，精密阀控 | TPU95, 伺服阀, PID优化 |
5.3 故障诊断与维护
常见问题处理经验：
位移异常减小
- 检查气密性(肥皂水检测)
- 验证材料塑性变形(永久变形>5%需更换)
- 清理折叠结构异物
响应速度下降
- 检查阀件响应(阶跃测试)
- 排查管路堵塞(压降测试)
- 评估材料老化(硬度变化)
控制振荡
- 重新整定PID参数(ziegler-nichols法)
- 增加阻尼(调节排气节流)
- 检查传感器安装(消除机械间隙)
在实际部署中发现，定期(每200小时)进行以下维护可延长使用寿命30%以上：
- 折叠结构润滑(硅基润滑剂)
- 气路过滤器更换
- 材料应力松弛处理(24小时无负载恢复)
通过三年期的跟踪测试，优化后的SSA执行器在康复机器人应用中实现了平均无故障工作时间(MTBF)超过5000小时，验证了其工程可靠性。未来将通过材料改性(如纳米复合材料)和结构创新(多自由度折叠)进一步提升性能边界。
