Abstract: It is an important means to simulate the deformation process of metal foam with low density numerically for investigating their mechanical behavior and failure. Based on the geometric model of Kelvin containing multiple cells, the finite element method was used to obtain the compressive yield strengths and to simulate the deformation curves of closed-cell metal foams. The deformation processes of the Kelvin foam models were simulated at the same time. During the calculation, the material properties of nonlinearity and geometric nonlinearity were taken into account,also, the influences of the number of cells on the foam properties were discussed. In order to investigate their influences on the mechanical behaviors of foams, the two constitute relations, perfect elastic-plastic relation and linearly hardening elastic-plastic relation, were adopted. The results show that the calculated values agree well with the theoretical results ever reported by other researchers, but the values are higher than the former experimental results.