Abstract: The constitutive matrix of an anisotropic membrane material applied to triangular membrane elements of Vector Form Intrinsic Finite Element (VFIFE) is derived. The modified constitutive matrix of a wrinkled anisotropic membrane is also derived, and a more reasonable method is proposed for the calculation of the angle between principle stress and local coordinate system. On this basis, a computer program for loading analysis and wrinkling analysis of anisotropic membrane structures is developed, and numerical examples are provided. It is shown from numerical results that both loading analysis and wrinkling analysis of anisotropic membrane structures can be performed well by the developed VFIFE program, verifying the validity of the theoretical derivation and the reliability of the computer program. It is also indicated that anisotropic properties and the wrinkling effect in membrane materials have a significant influence on the structural behavior of membrane structures.