Thin solar sail membranes of very large span are being envisioned for near-term space missions. One major design issue that is inherent to these very flexible structures is the formation of wrinkling patterns. Structural wrinkles may deteriorate a solar sail's performance and, in certain cases, structural integrity. In this paper, a geometrically nonlinear, updated Lagrangian shell formulation is employed using the ABAQUS finite element code to simulate the formation of wrinkled deformations in thin-film membranes. The restrictive assumptions of true membranes, i.e. Tension Field theory (TF), are not invoked. Two effective modeling strategies are introduced to facilitate convergent solutions of wrinkled equilibrium states. Several numerical studies are carried out, and the results are compared with recent experimental data. Good agreement is observed between the numerical simulations and experimental data. Introduction The exploitation of solar energy for the purpose of near-term space exploration presents a viable and attractive possibility in the minds of NASA scientists and engineers. The specific propulsion technology is called solar sails. Very large, ultra-low-mass, thinpolymer film (gossamer) structures are now being designed and tested for a wide variety of space exploration missions. The difficulties associated with conducting tests in a weightless environment place greater emphasis on the reliance on computational methods. Other gossamer structures that possess similar structural characteristics include inflatable space antennas, sun shields, and radars. A concise overview of gossamer structures and related technologies can be found, for example, in [1].A solar sail can gain momentum from incidence of sunlight photons at its surface. Since the momentum carried by an individual photon is very small, the solar sail must have a large surface area and a low mass, so that sufficient acceleration can be generated. Also, a solar sail requires a highly reflective surface that has minimal wrinkling and billowing under operational conditions. In the presence of significant wrinkling and billowing, the solar sail may lose efficiency, as compared to a flat sail, due to the oblique incidence of photons. The billowing problem may be overcome by applying sufficient tension to the sail membrane. Higher tensile stresses, however, are commonly accompanied by greater amplitude and longer structural wrinkles.A thin-film solar sail is a classical two-dimensional structure, with a thickness that is much smaller than its lateral dimensions. Since the bending stiffness is negligibly small compared to its membrane stiffness, the load-carrying capability using thin, low-modulus films is predominantly due to tensile membrane stresses. One key response phenomenon intrinsic to this class of structures is structural wrinkling. Structural wrinkles are local post-buckling patterns that are manifested by geometrically large transverse deformations whose magnitudes are much larger than the membrane thickness. Their formation is...