A novel homotopy-based decoupled wavelet numerical strategy is developed to investigate the large-deflection bending of variable-thickness anisotropic thin plate constituted by heterogeneous material resting on orthotropic foundation. The anisotropic plates are described undergoing sinusoidally initial imperfections with inhomogeneous thickness and exponentially varying elastic modules, while the subgrade reaction involves the nonlinear support effects by Winkler-type foundation along with Pasternak-type orthotropy. Geometric nonlinearity by uneven thickness and large deformation, material nonlinearity by heterogeneity and contact nonlinearity of orthotropic foundation are simultaneously considered with the highly coupled and strongly nonlinear variable-coefficient Von Kármán equations derived. The accuracy and convergence of the present formulation are validated in excellent agreement with published results, while the new results of extreme bending solutions are provided, which are few given by other methods and further verify the effectiveness and superiority of proposed wavelet procedures. Parametric investigations of influences of material orthotropy and inhomogeneity, unevenness of thickness, amplitude of initial imperfection along with orthotropy of nonlinear Winkler–Pasternak foundation have been studied, which are of great significance on the extremely large bending properties.