Metamaterial structures that exhibit negative stiffness (MSNS) behavior have shown the potential to be used as energy dissipation systems in different applications. Negative stiffness honeycomb structures consist of multiple pre-buckled beams that snapping through from one pre-buckled mode to another buckled mode, when subjected to transverse loads, displaying negative stiffness, and dissipating a significant portion of the input energy. This paper presents finite element models (FEMs) developed to investigate the performance of negative stiffness pre-buckled beams and honeycombs. The FEMs were validated using experimental results of single pre-buckled beams and honeycomb structures. The experimental testing of the pre-buckled beams was presented in the current study while the results of the honeycomb were available in the literature. The FEMs have the ability to capture the transition of the elements from one mode of buckling to another mode as well as the force thresholds, energy dissipation values, and residual displacements. The FEMs were extended to investigate the effect of different parameters such as boundary condition, apex height-to-beam thickness ratio, and constraining the second mode of buckling on the behavior of negative stiffness honeycombs.