This work is intended to model the origami metamaterial to predict the static/dynamic response of the origami metamaterials. As part of the modeling, a representative unit cell (RUC) was developed to predict the mechanical properties of origami structures. The actual skeleton of the origami structure should be representative in order to capture the mechanics of deformation, to predict the overall mechanical properties and behavior. Shell theories are suitable for modeling the origami structure since all individual segments are thin or moderately thick compared to each cell size. Thick shell elements adopted for the analysis. This element is applicable for thick and thin shell applications for permitting large strains. Its present formulation includes a stabilization parameter that effectively eliminates artificial deformation modes. The periodic boundary conditions are implemented in a way that the relative displacement between equivalent nodes on two opposite edges is zero. Opposite edges are separated by the characteristic length (L) that is the length of the RUC sides. Fused deposition modeling was used to fabricate the samples using polymeric material, and then the samples tested for mechanical properties to validate the numerical results.