In the stability analysis of footing, the determination of the ultimate load also known as collapse load is important. In this study, the collapse load has been calculated for a strip footing of width B surfaced over anisotropic undrained clay underlain by a circular void. The study follows lower bound finite element limit analysis incorporating an adaptive meshing technique, and results obtained using a finite element model have been compared with the available solutions from conventional bearing capacity theory and also with those from previous research. The shear strength of undrained clay was assumed to be directionally dependent and varies linearly with depth. The influence of various parameters such as anisotropic and heterogeneous coefficients, diameter of void to width of footing ratio, horizontal and vertical distance of void from the centerline of the footing, and ground surface, respectively, has been evaluated in this study. The numerical solutions of collapse load have been presented in non-dimensional form, known as bearing capacity factor Nc, which can be directly used by the site/practitioner engineers. The magnitude of Nc increases significantly with the increase in anisotropic and heterogeneous parameters. The heterogeneous and anisotropic behavior in clay has been taken into account by varying the undrained strength linearly with depth and changing the shear strength in relation to the major principal stress direction, respectively, and the effects due to voids are taken into account by varying horizontal distance from centerline of model, vertical depth, and diameter of voids. Typical failure patterns have also been discussed.