A finite elements-based contact model is developed to predict load distribution along the spline joint interfaces; effects of spline misalignment are investigated along with intentional lead crowning of the contacting surfaces. The effects of manufacturing tooth indexing error on spline load distributions are demonstrated by using the proposed model.