Four dimensional (4D) textiles consist of a textile with a printed polymeric grid, where deformation of the grid is induced by introducing residual stresses in the textile, in this case through pre-stretching of the textile substrate prior to printing. The fourth dimension refers to the ability of the structure to change shape over time by changing the residual stress in the textile. In order to design a useful component for a specific application, the material properties of constituents, direction and amount of residual stress, anisotropy of the textile substrate, geometry of the printed polymer, and pattern of the printed grid can all be altered. Due to the large amount of design variables involved, a validated modeling technique that can account for the complex material behavior of the soft, flexible textile under large strains and deformations along with the bifurcation and stability behavior of buckling beams is needed. In this study, an initial model was created to capture the time-independent buckling behavior and compared with experiments for rectangular elements of varying geometry and pre-strain. Once the model was calibrated, it fit experiments well, although additional advancements must be made to predict the nonlinear behavior of a wide variety of architectures with further accuracy. Finally, modeling strategies for large grids are laid out and discussed, and preliminary results are shown. Although this model did not include the transient nature of shape change, it can serve as a precursor for designing 4D Textiles and eventually predicting their time-dependent behavior under loading changes.