火星车车轮与控制系统协同设计优化方法
1. 火星车车轮与控制系统协同设计概述
在火星探测任务中，火星车的移动性能直接决定了科学探测的范围和效率。传统设计方法通常将车轮机械设计与控制系统开发分离进行，这种割裂的设计流程往往导致系统整体性能受限。我们团队开发的这套基于贝叶斯优化的协同设计框架，创新性地将车轮几何参数与转向控制参数纳入统一优化流程，通过高保真度闭环仿真实现了系统级性能优化。
这个设计框架的核心突破在于解决了三个关键难题：首先，采用连续介质地形力学模型(CRM)替代传统的离散元方法(DEM)，将单次仿真的计算时间从20小时缩短到10分钟量级；其次，建立了包含8个设计变量的参数化模型，涵盖车轮半径、宽度、轮齿特征等机械参数以及转向PID控制参数；最后，开发了多目标优化函数，平衡了行驶速度、轨迹跟踪精度和能耗三个关键性能指标。
2. 技术方案设计与实现
2.1 仿真平台架构
我们基于Chrono物理引擎构建了完整的仿真系统，主要包含三个核心模块：
车辆动力学模型：精确复现了1/6比例火星车(ART)的双横臂悬架、转向机构和传动系统。模型通过实验数据校准，最大程度保证了仿真真实性。
连续介质地形模型(CRM)：采用光滑粒子流体动力学(SPH)方法模拟可变形地形，将土壤视为连续介质而非离散颗粒。这种方法的计算效率比DEM高2-3个数量级，同时保持了足够的精度。
协同仿真框架：通过Chrono::FSI模块实现车辆与地形的双向耦合。系统采用非阻塞并行计算，车辆状态和接触力每0.5ms同步一次，确保实时交互的准确性。
关键技术细节：CRM模型使用Jaumann客观应力率描述土壤本构关系，结合μ(I)流变学模型捕捉颗粒材料的率相关摩擦特性。空间离散采用RK2格式，时间步长2×10⁻⁴秒，每个SPH粒子对应一个GPU线程实现并行计算。
2.2 参数化建模
2.2.1 车轮几何参数化
我们采用5个关键参数描述车轮几何特征：
- 外半径(ro)：0.05-0.14m
- 宽度比(wr)：宽度与外半径比值(0.7-1.4)
- 轮齿比(gr)：轮齿高度与外半径比值(0-0.3)
- 轮齿数量(ng)：2-16个
- 轮齿角度(αg)：45-135°
这种参数化设计在保证设计自由度的同时，考虑了3D打印制造约束。特别值得注意的是轮齿角度参数，通过允许轮齿扇形展开，在不增加基座宽度的情况下提高了有效接触面积。
2.2.2 控制系统参数化
转向控制系统采用PID控制器，优化以下参数：
- 比例增益(Kp,s)：0-15
- 积分增益(Ki,s)：0-5
- 微分增益(Kd,s)：0-5
油门控制参数保持固定，因为在松软地形上系统通常工作在扭矩饱和状态，油门控制对整体性能影响有限。
2.3 贝叶斯优化框架
2.3.1 目标函数设计
复合目标函数包含三个加权项：
J = 10×(ws·telapsed/tideal + wt·erms/eideal + wp·Pavg/Pideal)其中：
- 时间项(telapsed)：实际行驶时间与理想时间(5m/s)的比值
- 跟踪项(erms)：轨迹跟踪的RMS误差，参考值0.01m
- 功率项(Pavg)：平均机械功率，参考值20W
权重配置根据优化阶段动态调整：直线牵引测试(ws=1)，联合优化(ws=0.4, wt=0.2, wp=0.4)，纯控制优化(wt=1)。
2.3.2 优化流程
初始化阶段：使用Sobol序列生成1200-1800个初始样本，均匀覆盖设计空间。
代理模型训练：采用具有Matérn-5/2核的高斯过程建模目标函数，整数参数ng在优化中视为连续变量，评估时取整。
样本采集：使用LogNoisyEI采集函数平衡探索与开发，每次迭代选择最具潜力的候选设计进行仿真评估。
迭代优化：重复更新代理模型和采集新样本，直到达到3000次仿真预算(约5-9天计算时间)。
3. 优化结果与分析
3.1 车轮单独优化
在直线牵引测试中，我们对车轮几何参数进行了独立优化，获得了以下重要发现：
参数敏感性分析：
- 外半径(ro)贡献了77%的主要效应和87%的总效应方差，是影响牵引性能的决定性因素
- 宽度比(wr)贡献约7%主要效应，表现出明显的尺寸效应
- 轮齿参数(gr, ng, αg)合计贡献约13%效应，其中轮齿高度比数量更重要
最优参数分布：
- 外半径最优区间0.09-0.12m，过大反而降低性能
- 宽度比趋向上限1.4，说明宽轮在松软地形优势明显
- 轮齿数量14个、高度比0.3时获得最佳牵引力
实验验证：
排名 仿真时间(s) 实测时间(s) 1 (最优) 2.31 8.0 100 2.63 10.3 1000 3.20 15.0
虽然绝对时间存在差异，但性能排序完全一致，验证了仿真优化的有效性。
3.2 联合优化对比
我们比较了两种优化策略在正弦轨迹测试中的表现：
联合优化：
- 计算预算：3000次仿真(1800 Sobol + 1200 BO)
- 优势：发现非直观的设计组合
- 局限：控制器参数可辨识性低(Kd,s仅2.4%贡献)
分阶段优化：
- 阶段1：2000次仿真优化车轮(ws=0.5, wp=0.5)
- 阶段2：1000次仿真优化控制(wt=1)
- 优势：控制器参数贡献度提升至91%(Kd,s)
- 结果：跟踪误差降低37%，能耗减少22%
实测表明，分阶段策略更适合工程实践，因为它更符合传统的"先机械后控制"开发流程，且能保证各子系统参数得到充分优化。
4. 工程实践建议
基于大量仿真和实验，我们总结了以下实用经验：
车轮设计准则：
- 松软地形优先考虑外半径(0.1m左右最佳)
- 增加宽度比可显著改善稳定性但会增加质量
- 轮齿高度比数量更重要，建议0.25-0.3
- 轮齿角度45°时牵引性能最佳
控制参数调试：
- 微分增益对跟踪精度影响最大
- 比例增益过高会导致振荡
- 积分项在低速移动中作用有限
仿真设置技巧：
- SPH粒子间距建议0.5cm
- 时间步长采用CFL条件自动调整
- 边界条件使用Adami插值法
硬件实现注意事项：
- 3D打印车轮需考虑层间结合强度
- 实际土壤参数与仿真存在差异，建议保留15%余量
- 在轨运行时需考虑重力差异(火星重力约地球1/3)
这套框架已经开源发布，包含完整的仿真模型、优化脚本和CAD模板。对于想尝试的研究者，建议从简化模型入手，先验证单轮性能再扩展到整车。我们也提供了参数敏感性分析工具，可快速识别关键设计变量，大幅提高优化效率。