The human knee joint is the intermediate joint of the lower limb that is the largest and most complex joint in the body. Understanding of joint-articulating surface motion is essential for the joint wear, stability, mobility, degeneration, determination of proper diagnosis and so on. However, many studies analyzed the passive motion of the lower limb because of the skin marker artefact and some studies described medial and lateral condyle of a femur as a simple sphere due to the complexity of geometry. Thus, in this paper, we constructed a three-dimensional geometric model of the human knee from the geometry of its anatomical structures using non-uniform B-spline surface fitting as a study for the kinematic analysis of more realistic human knee model. In addition, we developed and verified 6-DOF contact model of the human knee joint using $C^2$ continuous surface of the inferior region of a femur, considering the relative motion of shank to thigh during locomotion.