Automatic joint motion planning is very important in robotic wheel hub polishing systems. Higher flexibility is achieved based on the joint configuration with multiple solutions, which means that the robot has kinematic redundancy for machining tasks. Redundant joints can be used to optimize the motion of the robot, but less research has been done on multi-dimensional redundant optimization. In this paper, a 6-axis robot with a 3-axis actuator is designed for wheel hub polishing. We propose an automatic joint motion planning method for a nine-axis industrial robot to achieve the shortest processing time. Firstly, offline programming is designed to generate paths for the complex surface of the hub. In order to reduce the machining path points on the surface of the hub, a improved Douglas-Peucker (DP) algorithm is proposed, which can take into account the change of the path point posture. Secondly, the Greedy Best First Search (GBFS) and Sine cosine algorithm (SCA) are combined to find the optimal joint motion efficiently. Moreover, we use nested SCA for comparison to test whether the combined algorithm can avoid local optima. Finally, the performance and computational efficiency of the method are validated in both simulation and real environments based on the hub surface.