A three-dimensional multibody frictional contact problem in a CANDU nuclear fuel rod is studied in this thesis. To account for the effect of friction on fuel rod deformations, a computational model is developed based on an existing model. For the contact algorithm, a set of complementary equations involving complementarity variables inequality is introduced to deal with the three-dimensional frictional contact. The algorithm is implemented into a computer program. Three cases are studied to show the feasibility and efficiency of the numerical scheme for handling frictional contact among multiple bodies. The application of the research work offers the assessment of deflection as to safety evaluation for the nuclear industry under different scenarios.