基于优化的蜗杆型齿轮加工中通过刀具几何与进给协调实现扭转控制的方法

《Machines》：An Optimization-Based Approach to Twist Control Through Tool Geometry and Feed Coordination in Worm-Type Gear Generation

　　

在精密齿轮制造中，齿面上的纵向鼓形修形通常通过蜗杆型展成过程中应用对角进给来实现，使用诸如变齿厚滚刀（variable-tooth-thickness hobs）和可修整磨削蜗杆（dressable grinding worms）等刀具。然而，由于展成刀具的几何参数与机床进给设置在底层展成运动学中强烈耦合，精确的扭转控制仍然困难。此外，在标准刀具状态附近，直接数值优化变得不可靠，因为对角进给系数的灵敏度退化，常规线性化求解器可能失效。为解决这些问题，本研究提出了一种针对扭转约束的蜗杆型齿轮展成的多变量优化框架。研究人员开发了一种迭代奇异值分解（SVD）方案来构建和更新灵敏度矩阵，同时引入热启动连续策略以克服局部奇异性并提高数值稳健性。还提出了两种对角进给系数的闭环表达式，分别对应最小SVD地形残差和最小齿面扭转，作为实际有用的初始估计。对60个案例的参数扫描进行的数值验证显示，在所测试范围内最大相对误差低于1.6%。该框架协调了刀具几何设计与对角进给选择，以生成具有指定鼓形修形特性、同时满足特定扭转要求并限制所需对角位移的齿面。数值示例表明，相对于Hsu和Fong的分析设置，迭代框架将均方根（RMS）地形误差从1.14 μm降低到0.027 μm。这些结果表明，所提方法为先进CNC齿轮展成中的扭转控制和工艺参数设计提供了可靠的计算基础。从制造角度看，由于三个设计标准仅通过调整机床上的对角进给比即可实现，单个展成刀具设计可满足一系列鼓形修形和扭转要求而无需换刀，从而减少生产中的装夹和刀具成本。

**论文解读：基于优化的蜗杆型齿轮加工中通过刀具几何与进给协调实现扭转控制的方法**

**研究背景与问题**
在精密齿轮制造领域，蜗杆型展成工艺（如连续展成磨削、滚齿等）是生产精密圆柱齿轮的高效方法。然而，实际传动中轴的挠曲、装配误差和负载变形常导致接触斑偏向齿边。为局域化接触区并缓解边缘接触，通常在齿面上引入纵向修形（如鼓形修形）。但传统的纵向修形方法（如中心距变化修形）会引发齿面扭转问题，即修形量沿齿深方向分布不均，这源于展成运动学中刀具几何与机床进给设置的强耦合。早期研究通过变齿厚刀具结合对角进给可减少扭转，但现有方法多为分析性参数选择或局部校正策略，缺乏同时优化刀具几何与进给参数的数值稳健框架。特别是，当刀具处于标准均匀厚度状态时，对角进给对生成齿面的一阶灵敏度严格为零（敏感性退化），导致直接线性化优化失效。Hsu和Fong虽推导了对角进给系数的分析表达式，但未处理这种退化作为数值优化问题，且其系数在低螺旋角工件齿轮等极限情况下发散。针对这些不足，本研究提出了一个多变量优化框架，以在指定修形要求下可靠控制扭转。

**研究概述与结论**
研究人员开发了一套结合闭环分析指导与迭代数值精化的优化框架，用于蜗杆型齿轮展成中的扭转控制。该框架包含三个核心组件：首先识别了标准刀具状态下灵敏度退化的结构性质；其次引入迭代奇异值分解（SVD）方案与热启动连续策略以克服退化并提高数值稳健性；最后提出了两种对角进给系数的闭环表达式，分别对应最小SVD地形残差和最小齿面扭转。通过60个案例参数扫描验证，所提表达式（公式28和30）的最大相对误差分别低于0.6%和1.6%。数值示例显示，对于典型乘用车螺旋齿轮，迭代SVD方案将均方根（RMS）地形误差从使用Hsu-Fong系数的1.14 μm降至0.027 μm，同时最小扭转系数进一步将扭转指数从6.52 μm降至5.40 μm（代价是RMS误差增至0.31 μm）。该框架还考虑了扭转约束与磨削蜗杆可用位移长度的几何约束，形成了设计可行区间。论文发表于《Machines》。

**主要关键技术方法**
本研究基于展成运动学模型，将蜗杆型展成工具与工件齿轮建模为交叉轴螺旋齿轮副。关键方法包括：（1）灵敏度矩阵构建：通过前向差分计算7个设计变量（6个刀具几何修形系数和1个对角进给系数）对90个齿面偏差点的雅可比矩阵；（2）迭代奇异值分解（SVD）方案：在每次迭代中计算灵敏度矩阵的SVD，对低于相对阈值（10^-6倍最大奇异值）的奇异值截断，利用Moore-Penrose伪逆求解最小范数最小二乘校正，并加入阻尼因子防止过冲；（3）热启动连续策略：在工件齿面评估网格上，从中心纵向截面开始，将前一网格点的收敛解作为下一网格点的初始猜测，以稳定包络求解器；（4）闭环表达式推导：基于模态分解，通过参数扫描拟合得到最小残差和最小扭转的对角进给系数公式（依赖于工件齿轮工作螺旋角^β_w、刀具导程角^γ_t和法向压力角^α_n）；（5）扭转约束工作流：先由闭环表达式确定初始对角进给系数，再通过迭代SVD精化完整设计向量，并检查扭转指数是否满足允许上限。本研究无样本队列来源，全部为数值仿真。

**研究结果**
**第4.1节：基于默认案例的五种方法比较**
通过比较常规中心距修形、Hsu-Fong均匀系数^c_h、最小残差系数^c_res、最小扭转系数^c_tw和迭代SVD方案（初始化自^c_res），研究人员发现：（1）常规方法产生13.59 μm扭转（约为30 μm鼓形修形的45%），RMS误差为8.62 μm；（2）Hsu-Fong系数将扭转降至7.59 μm，RMS误差为1.14 μm；（3）最小残差系数^c_res（-2.018）将RMS误差降至0.027 μm（改善42倍），扭转降至6.52 μm；（4）最小扭转系数^c_tw（-2.283）将扭转进一步降至5.40 μm，但RMS误差增至0.31 μm；（5）迭代SVD方案收敛至与^c_res相同的操作点，验证了闭环表达式的准确性。收敛速度迅速，通常在1-3次迭代内完成。

**第4.2节：参数空间验证**
对参数扫描（工件齿轮工作螺旋角^β_w从15°到30°，刀具工作螺旋角^β_t从-5°到-15°，法向模数^m_n从1.5 mm到4.5 mm，法向压力角^α_n从15°到30°，工件齿轮齿数^Z_w从20到80等）的验证表明：闭环表达式^c_res最大相对误差为0.6%，^c_tw最大为1.6%；两种表达式对齿数、模数、面宽和目标鼓形幅度的敏感性几乎为零。手性符号约定（同手性取负号，异手性取正号）通过单个反手性配置验证。

**第4.3节：极限情况与实际边界**
Hsu-Fong系数在^β_w = ^γ_t时发散（如低螺旋角工件齿轮用多头磨削蜗杆），而本研究提出的闭环表达式在所有有效参数范围内有界。在直齿轮极限（^β_w → 0）下，两种表达式发散，但此时对角线进给无法产生纵向修形，扭转概念无意义，框架不适用。

**第4.4节：扭转约束可行区域**
通过扫描区间[^c_tw, ^c_res]并绘制均方根误差与扭转指数的帕累托前沿，研究人员发现两者呈光滑权衡关系。例如，对于默认案例，施加6 μm扭转约束将可行区间限制在约[-2.10, -2.05]内，对应RMS误差约0.15-0.31 μm。

**第4.5节：磨削蜗杆几何约束**
考虑磨削蜗杆可用位移长度^L_w时，刀具施加的对角进给系数为^c_w = ±^L_w/^b（^b为面宽）。对于默认案例，^c_res和^c_tw对应所需蜗杆位移约94 mm和107 mm。若现有蜗杆位移不足，则需接受更高扭转、利用刀具其他修形系数补偿，或采用多道次策略。最终可行区间为扭转约束区间与几何约束区间的交集。

**讨论与结论**
讨论部分指出：灵敏度退化是结构性的而非数值性的，对角进给方向始终是设计空间中的“软方向”，闭环表达式确定了该方向上的唯一极小值，迭代SVD负责其余六个“硬方向”的精化。该框架在运动学层面与工艺无关，适用于滚齿、连续展成磨削、刮齿等工艺，但实际实现依赖于具体工艺（如刀具修形方式、电子齿轮箱编程）。局限性包括：闭环表达式基于经验观察和数值验证，缺乏严格分析推导；仅验证了单一纵向鼓形修形目标；手性符号约定验证不充分；未考虑交叉角调制作为独立变量；全部为计算验证，缺乏实验确认。结论部分翻译如下：本研究开发了一种基于优化的蜗杆型齿轮展成扭转控制框架，结合了闭环分析指导与迭代数值精化。主要贡献包括：证明了标准均匀厚度刀具状态下对角进给系数灵敏度的严格退化是结构性的；提出了迭代奇异值分解方案结合热启动连续策略以克服退化；提出了两种对角进给系数的闭环表达式（最小残差和最小扭转），并在测试范围内数值验证；制定了扭转约束工作流，适用于多种蜗杆型展成工艺；数值示例表明迭代SVD可将RMS误差从1.14 μm降至0.027 μm，最小扭转系数进一步将扭转指数从6.52 μm降至5.40 μm（代价是RMS误差增至0.31 μm）；揭示了最小化地形残差与最小化扭转之间的权衡，两种系数界定了设计空间的单参数族，允许从业者根据主导设计需求选择操作点。未来工作将致力于从啮合方程直接推导闭环表达式，扩展至更一般的目标地形，并纳入交叉角调制作为额外设计变量。