This paper extends sublaminate-based variable kinematics plate finite elements towards the buckling analysis of composite structures. Robust locking-free 4-node and 8-node interpolations are used to build the finite element matrices in terms of fundamental nuclei, invariant with respect to the employed kinematic model of the composite plate. Geometric nonlinearities are accounted for in the von Kármán sense. The classical linearized stability analysis is mainly conducted for sandwich panels, for which different kinematic assumptions are employed for the skins and the core. Global buckling of the sandwich panel as well as short-wavelength wrinkling of the skins are investigated by referring to a variety of case studies, including homogeneous and laminated skins as well as isotropic and orthotropic cores. Convergence studies are performed to establish the minimum number of elements for the local instabilities to be grasped. The proposed computational framework requires a simple 2D mesh and is capable of providing quasi-3D response patterns. This is as demonstrated by numerous case studies that address the transition from global to local buckling, the onset of wrinkling in flat sandwich panels with anisotropic skins under various loading conditions as well as the local face sheet instability occurring in sandwich panels working in bending. It is concluded that he proposed FEM-based tool is computationally efficient and can be advantageously employed in pre-sizing design phases without resorting to full 3D models.