With the increasing demand in ocean engineering field for the ultimate limit state (ULS) analysis of thick plates, the solid (3D) element has been more frequently used in the ULS analysis. Compared to 2D element, 3D element can better consider the transverse shearing effect during the analysis which usually occurs on thick plates, due to its reliable reflection of geometrical characteristics and kinematic equations of the structures. Besides, 3D element can also facilitate the introduction of influencing factors into the ULS analysis, such as the erosions, cracks and residual stress. This may improve the simulation precision of ULS analysis of thick plate models. Until now, there is a lack of 3D element modelling technique in the ULS analysis of flat/stiffened plates. The present study aims to provide helpful information on 3D elements in the ULS analysis of flat plates under axial compression load. A total of 350 plate scenarios were adopted in the parametric study to consider the effects of 3D element shape (αyz). It is found that the element shape significantly influences the ULS analysis of flat plates, where planar-like 3D element is not recommended. An empirical formula for determining the optimal 3D element shape of the finite element (FE) model is proposed based on the parametric results. A guidance for implementing the 3D elements is then documented, which may help engineers further understand the ultimate strength characteristic of very thick flat plates.