The dimensional reduction of a hyperelastic coated fabric from 3D to 2D is accomplished asymptotically using the Variational Asymptotic Method (VAM). This method integrates constrained calculus of variations and asymptotics. The VAM bifurcates the analysis into two: a 1D through-the-thickness analysis and a 2D reference surface analysis. The 1D analysis leads to the derivation of an asymptotically correct 3D warping functions and a 2D non-linear constitutive law. The 2D non-linear reference surface analysis utilizes the derived 2D non-linear constitutive law to obtain 2D displacements and strains through the 2D non-linear FEA. The classification of 3D strain energy density into distinct orders is enabled by introducing two intrinsic small parameters: 1) a geometric small parameter denoted by the ratio of thickness to characteristic length (h/l≪1), and 2) a physical small parameter that ensures the largest component of 3D strain is restricted to 20 percent, which is less than 1. The model takes into account both geometric and material nonlinearities. The strain energy function, which describes the anisotropic characteristics of the coated fabric has contributions from the strain energies of fiber, matrix, and fiber-fiber interaction. The findings of the study include analytically derived 3D warping functions, a 2D nonlinear constitutive law, and the prediction of warp and weft stresses and strain for a biaxial loaded tensile specimen. These findings align with the experimental outcomes.