This paper developed a strength constrained topology optimization method for hyperelastic structures with large deformation-induced frictionless contact. The Neo-Hookean hyperelastic constitutive equation is adopted as the material model that incorporates both material and geometric nonlinearity. The large deforming contact is described by the node-to-segment algorithm. The topology optimization model is developed under the SIMP framework. Given the SIMP-related fictitious domain, we develop a nodal variable-based interpolation scheme to build smooth correlation between the contact force and nodal density variables, thus enabling the smooth transition and variation of contact conditions. Moreover, to guarantee the structure strength, local strain energy constraints are formulated and aggregated by the P-norm function. The details of sensitivity analysis are explicitly explained. At last, several benchmark examples are conducted that validate the proposed method.