Finite element analysis(FEA) is the most popular numerical method to simulate plasticity-induced fatigue crack closure and can predict fatigue crack closure behavior. Finite element analysis under plane stress state using 4-node isoparametric elements is performed to investigate the detailed closure behavior of fatigue cracks and the numerical results are compared with experimental results. The mesh of constant size elements on the crack surface can not correctly predict the opening level for fatigue crack as shown in the previous works. The crack opening behavior for the size mesh with a linear change shows almost flat stress level after a crack tip has passed by the monotonic plastic zone. The prediction of crack opening level presents a good agreement with published experimental data regardless of stress ratios, which are using the mesh of the elements that are in proportion to the reversed plastic zone size considering the opening stress intensity factors. Numerical interpolation results of finite element analysis can precisely predict the crack opening level. This method shows a good agreement with the experimental data regardless of the stress ratios and kinds of materials.