Contact pressure distribution over the active area of proton exchange membrane fuel cell (PEMFC) has significant effects on the performance of PEMFCs. Even clamping pressure over the membrane electrode assembly (MEA) affects contact resistance, characteristics of porous media and sealing task. This paper develops a PEM fuel cell model to study the contact pressure distribution over the membrane electrode assembly using finite element model. At first, the three-dimensional model of a single cell was reduced to a two-imensional model to decrease the calculation time. After validation of the obtained results via the pressure sensitive films, the effect of some parameters such as thickness and material of the end plates, sealant hardness, number of the stack's cells and position of the cell on the contact pressure distribution over the MEA were investigated. The results reveal that optimizing mentioned parameters leads to design PEM fuel cells with proper contact pressure distribution over the MEA. Copyright (C) 2015, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.