为更真实地反映负载工况下谐波齿轮的啮合特性,提出了一种基于空载侧隙和周向线性啮合刚度的理论迭代算法,用于计算随负载变化的负载侧隙和啮合力分布。利用精确算法计算了装配状态下的空载侧隙;建立了三维实体单元的柔轮有限元模型,在柔轮齿上逐一施加啮合力,利用接触分析求解了双圆盘波发生器作用下反映柔轮齿廓啮合点周向刚度特性的啮合刚度矩阵;根据装配状态下的空载侧隙,迭代计算了在逐步增大的刚轮转动位移下的负载侧隙、啮合力及负载扭矩。由于理论算法无法考虑啮合刚度的非线性影响,为验证本文算法,利用有限元非线性接触分析方法,通过定义齿廓间的接触关系,数值求解了负载侧隙和啮合力。对比研究表明:在额定工况下,理论算法给出的啮合力分布与有限元分析结果基本吻合,但当负载较大时,理论算法的啮合刚度偏大,造成啮合力幅值偏高。