《Actuators》：Structural Optimization of Harmonic Drive Flexspline to Improve Fatigue Lifetime
编辑推荐：
为克服现有谐波驱动（HD）系统研究的局限性，例如缺乏柔轮（FS）结构参数对使用寿命影响的系统探索以及参数优化策略不足，本研究建立了一个综合研究框架，集成了有限元（FE）模拟、响应面优化和实验验证。选取了柔轮的三个关键结构参数，包括杯长（CL）、杯壁厚度（CT）
为克服现有谐波驱动（HD）系统研究的局限性，例如缺乏柔轮（FS）结构参数对使用寿命影响的系统探索以及参数优化策略不足，本研究建立了一个综合研究框架，集成了有限元（FE）模拟、响应面优化和实验验证。选取了柔轮的三个关键结构参数，包括杯长（CL）、杯壁厚度（CT）和齿宽（TW）作为设计变量，研究它们对HD系统位移、von Mises应力和疲劳寿命的影响。构建了可靠的HD系统有限元模型，验证了网格独立性，并使用伺服电机驱动测试平台进行了实验验证。结果表明，模拟与实测位移（4.7%）和应力（4.0%）的相对误差均在可接受的工程范围内（<10%），证实了模型的可靠性。系统分析揭示，柔轮寿命由全局应力分布和结构刚度共同决定，而非仅仅由局部齿根应力决定。增加杯长会降低齿根应力，但会增加杯壁弯曲应力，后者最终主导疲劳失效过程。基于Box-Behnken设计（BBD）响应面分析，确定了最优参数组合为CL=48 mm、TW=6 mm、CT=0.64 mm，实现了低系统应力（77.7 MPa）和长柔轮寿命（2950 h）的双重优化。相对于基线对照组（CL=48 mm、CT=0.46 mm、TW=6 mm），优化配置将系统应力降低了10.1%，疲劳寿命延长了1.7%。本研究阐明了结构参数调节HD系统性能的多因素耦合机制，为HD系统的结构设计和性能改进提供了坚实的理论基础和工程参考。
谐波驱动（HD）系统因其紧凑结构、大传动比和高精度等优势，广泛应用于工业机器人、航空航天及医疗设备等领域。柔轮（FS）作为核心输出部件，在波发生器（WG）作用下经历连续弹性变形，导致齿根和杯壁处产生循环应力集中，疲劳失效是主要失效模式。现有研究多聚焦于降低局部齿根应力，忽视了杯壁应力与全局结构刚度的耦合作用，难以明确结构参数如何协同调节整体寿命。此外，缺乏系统性的多参数优化策略，且有限元（FE）模型多未考虑完整组件接触与非线性行为。为填补这些空白，研究人员在《Actuators》上发表了该研究，旨在通过集成FE模拟、响应面优化与实验验证，阐明杯长（CL）、杯壁厚度（CT）和齿宽（TW）对HD系统位移、von Mises应力及疲劳寿命的影响机制，并寻求最优参数组合以实现低系统应力与长寿命的双重目标。研究结论表明，FS寿命由全局应力分布与结构刚度共同决定，最优组合（CL=48 mm, TW=6 mm, CT=0.64 mm）使系统应力降低10.1%，寿命延长1.7%，为HD系统结构设计提供了理论依据与工程参考。
主要关键技术方法包括：基于Solidworks v2024构建包含柔轮、刚轮、波发生器（含椭圆凸轮和球轴承）的完整HD系统几何模型；采用有限元法（FE）进行网格划分（四面体单元，平均质量0.82），并对柔轮齿和杯壁区域加密（0.3 mm），经网格独立性验证（粗、中、细三套网格，最大von Mises应力相对偏差3.2%）；设置边界条件（刚轮全约束，波发生器输入600 N·m全反向循环扭矩，柔轮自由旋转），并考虑非线性摩擦接触（摩擦系数0.15）和大变形效应；采用Box-Behnken设计（BBD）进行三因素三水平响应面分析，以柔轮疲劳寿命为响应指标，拟合二阶多项式回归模型（R2=0.962）；最后使用伺服电机驱动测试平台（含扭矩传感器、振动速度传感器HG-ZD-20A）进行位移与应力验证，重复测试5次取平均值。
**3.1 HD系统模拟位移分布**：通过FE模拟分析四组（对照组、CL=72 mm组、TW=18 mm组、CT=0.64 mm组）柔轮位移分布。结果表明，对照组柔轮齿位移约0.09 mm，CL=72 mm组最小（0.074 mm），但整体系统位移略高于对照组，说明系统机械响应是组件相互作用的综合结果。增加CL会降低杯壁弯曲刚度，导致杯体变形增大。
**3.2 HD系统模拟应力分布**：FE模拟显示，柔轮齿最大von Mises应力在TW=18 mm组最高（59.2 MPa），而CL=72 mm组最低（25.1 MPa）。但整体系统平均应力在CL=72 mm组最大（104.1 MPa），CT=0.64 mm组最小（77.7 MPa），比CL=72 mm组降低25.3%。所有组别柔轮齿应力均显著低于系统平均应力，表明柔轮齿并非主要承力部件。
**3.3 HD系统模拟寿命**：FE预测的柔轮齿寿命在TW=18 mm组最长（2950 h），CL=72 mm组最短（2510 h）。尽管CL=72 mm组齿根应力最低，但因其杯壁弯曲应力最大（104.1 MPa）且结构刚度过大，导致寿命最短。TW=18 mm组虽齿根应力最高，但寿命最长，证实FS寿命并非由局部齿根应力单独主导，而是由全局应力分布与结构刚度共同决定。
**3.4 HD系统响应面分析及参数优化结果**：基于15组BBD实验，建立FS寿命二阶多项式回归模型（R2=0.962，调整R2=0.947）。方差分析（ANOVA）显示模型高度显著（F=28.15，p<0.0001），失拟项不显著（p=0.183）。CL和CT对寿命影响最显著（p<0.0001），TW边缘显著（p=0.052），交互项CL×CT显著（p=0.012）。优化得到最优参数组合：CL=48 mm、TW=6 mm、CT=0.64 mm，此时系统应力77.7 MPa，FS寿命2950 h，同时实现了重量、惯量及传动精度的平衡。
**3.5 HD系统应力与位移分布验证**：对对照组进行实验验证，测量位移为0.085±0.002 mm，应力为90±5 MPa，与模拟值（0.089 mm、86.4 MPa）的相对误差分别为4.7%和4.0%。所有四组参数配置的位移相对误差均低于10%，满足工程接受标准，确认了FE模型的可靠性。
讨论部分指出，研究结果揭示了CL、CT、TW对HD系统机械性能的多因素耦合机制：增加CL降低齿根应力但升高杯壁弯曲应力，疲劳寿命由全局应力与结构刚度共同决定，而非局部齿根应力。与现有研究一致，但进一步表明了杯壁弯曲疲劳在长杯长时成为主导失效模式。响应面优化提供了定量交互模型，实现了局部应力降低与全局结构完整性的平衡。局限性包括：摩擦系数假设为常数、刚轮全约束、采用单轴疲劳准则、优化设计尚未完全实验验证。未来工作将聚焦瞬态接触动力学分析、多轴疲劳准则应用、摩擦系数实验表征及机器学习加速优化。
结论部分翻译：本研究建立并实验验证了可靠的HD系统有限元模型，融合了核心组件几何特征与界面非线性接触行为。网格收敛性分析确认基线网格（约70,000单元）平衡了效率与精度，基线与加密网格的von Mises应力相对误差低于5%。模拟位移和应力与实测值的相对误差分别为4.7%和4.0%，满足工程验证标准。研究阐明了CL、CT、TW如何调节HD系统的机械响应与疲劳寿命，揭示系统性能取决于组件相互作用而非单一部件行为，且过大的结构刚度会损害疲劳性能即使局部应力较低。通过响应面分析的多目标优化策略获得了最优参数组合（CL=48 mm、TW=6 mm、CT=0.64 mm），实现了低系统应力（77.7 MPa）和长柔轮寿命（2950 h）。与原始对照组设计（CL=48 mm、CT=0.46 mm、TW=6 mm）相比，系统应力降低了10.1%，预测寿命延长了1.7%。最后分析了模拟误差来源并提出了未来优化方向。本研究丰富了HD性能分析理论，为高精度传动应用提供了工程指导。
