Time-varying mesh stiffness calculation by a novel loaded tooth contact analysis model considering EHL contact for face-gear drive