Abstract. Two mathematical models are developed within the theoretical framework of large strain elasticity for the determination of upper and lower bounds on the total strain energy of a finitely deformed hyperelastic body in unilateral contact with a rigid surface or with an elastic substrate. The model problems take the form of two continuous optimization problems with inequality constraints, the solutions of which are used to provide an enclosure on the uniform external load acting on the body's surface away from the contact zone. 1. Introduction. In materials analysis and design, contact problems in elasticity are central to the modeling and investigation of many structural systems. In particular, the analysis of soft tissue biomechanics and the design of bioinspired synthetic structures involve nonlinear hyperelastic models for which the mathematical and numerical treatment poses many physical, theoretical, and computational challenges [2,17,29,31,32,33,39].Hyperelastic materials are the class of material models described by a strain energy density function [1,14,15,35,40]. For these materials, boundary value problems are often equivalent to variational problems, which provide powerful methods for obtaining approximate solutions. They can also be used to generate finite element methods for which the numerical analysis stands on the shoulders of the mathematical analysis for the elastic model [6,23,30,34].The complexity of contact problems modeling structural systems is generally associated with the detection of contacts and openings and the resolution of nonlinear equations for contact. For example, in systems formed from linear elastic bodies (e.g., ceramics, metals), surface roughness can impede active contact when a surface is pressed against another, and contact forces cannot be transmitted where surface separation occurs [16,18,19]. By contrast, nonlinear elastic bodies (e.g., rubber, soft tissue) are more pliable and thus capable of attaining more active support through which contact forces can be transmitted effectively [3,8,9,20]. These problems can be formulated and solved in the framework of variational inequalities [27], which originated in the paper by Fichera [11] on the existence and uniqueness of a solution to the celebrated Signorini problem with ambiguous boundary conditions [38]. While the Signorini problem consisted of finding the equilibrium state of an elastic body resting on a rigid frictionless surface, the influence of variational inequalities went beyond the