A numerical method is developed to study the effect of the compression deformation of the gas diffusion layer (GDL) on the performance of the proton exchange membrane fuel cell (PEMFC). The GDL compression deformation, caused by the clamping force, plays an important role in controlling the performance of PEMFC since the compression deformation affects the contact resistance, the GDL porosity distribution, and the cross-section area of the gas channel. In the present paper, finite element method (FEM) is used to first analyze the ohmic contact resistance between the bipolar plate and the GDL, the GDL deformation, and the GDL porosity distribution. Then, finite volume method is used to analyze the transport of the reactants and products. We investigate the effects of the GDL compression deformation, the ohmic contact resistivity, the air relative humidity, and the thickness of the catalyst layer (CL) on the performance of the PEMFC. The numerical results show that the fuel cell performance decreases with increasing compression deformation if the contact resistance is negligible, but an optimal compression deformation exists if the contact resistance is considerable. (C) 2007 Elsevier B.V. All rights reserved.