This paper presents an analytical model of the electrical contact resistance between the carbon paper gas diffusion layers (GDLs) and the graphite bipolar plates (BPPs) in a proton exchange membrane (PEM) fuel cell. The model is developed based on the classical statistical contact theory for a PEM fuel cell, using the same probability distributions of the GDL structure and BPP surface profile as previously described in Wu et al. [Z. Wu, Y. Zhou, G. Lin, S. Wang, S.J. Hu.J. Power Sources 182 (2008) 265-269] and Zhou et al. [Y. Zhou, G. Lin, A.J. Shih, S.J. Hu, J. Power Sources 163 (2007) 777-783]. Results show that estimates of the contact resistance compare favorably with experimental data by Zhou et al. [Y. Zhou, G. Lin, A.J. Shih, S.J. Hu. J. Power Sources 163 (2007) 777-783]. Factors affecting the contact behavior are systematically studied using the analytical model, including the material properties of the two contact bodies and factors arising from the manufacturing processes. The transverse Young's modulus of chopped carbon fibers in the GDL and the surface profile of the BPP are found to be significant to the contact resistance. The factor study also sheds light on the manufacturing requirements of carbon fiber GDLs for a better contact performance in PEM fuel cells. (C) 2009 Elsevier B.V. All rights reserved.