A new tooth pitting modeling method based on matrix equation for evaluating time-varying mesh stiffness. (December 2022)