Apologies if the title's a bit vague; I'll try to explain everything below, and please let me know if I should clarify anything.

I've been looking for a while at variational problems on polytopes. I'm particularly interested in the isoperimetric problem - maximizing volume relative to surface "area" under certain constraints, e.g. over polytopes of a fixed combinatorial type - and more generally, finding some sort of equivalent to surfaces of constant mean curvature in terms of the behaviour of the area functional.

My main goal has been to come up with a "nice" Riemannian metric on the space of polytopes, so that I can look at things like the gradient of area. Eventually I'd like to generalize objects like the shape operator in a way which relates to all this, so that the aforementioned gradient of area works out as being equivalent to mean curvature (I've asked some questions about this problem before, but all of them were vague and based on less understanding than I have currently, so I won't link to anything for now.)

I know that for a sufficiently smooth (hyper)surface in ℝ^{n}, since a normal exists at every point, variations of the surface can be thought of as normal vector fields; if we restrict ourselves to considering variations which preserve smoothness, then the vector field should also be smooth; and we can define an inner product on the space of variations by integrating the pointwise product of any two variations.

It feels like all this should extend to polytopes, but there's a problem in places where the surface bends sharply and no well-defined normal direction exists. The set of such points has measure zero so naively we might expect that we can ignore it when integrating; but I'm not so sure that we can, because we might have a variation where the instantaneous velocity of such points is infinite (even though it is finite for smooth points) and they therefore contribute a nonzero amount to the integral.

What I'd really like is some sort of general theory which encompasses both smooth surfaces and polytopes, as well as other surfaces "in between". It occurred to me to consider general surfaces as limits of smooth ones, but then how can you ensure that all the limits will be well-behaved? e.g. a sequence of smooth surfaces might converge to some other surface, but their normals might not converge at all! Or the normals might converge, but the inner product might not.

I also looked at the boundaries of convex bodies specifically. I know these to be fairly well-behaved, and it feels like it should be possible to extend the smooth theory to these surfaces at the very least. But even convex polytopes are not strictly convex, so they are arbitrarily close to non-convex bodies, so if we're looking at variations the non-convex theory still seems to be important.

My gut feeling is that there should be some space of surfaces which are not quite smooth (in general), but which are still sufficiently close to being smooth (as the boundary of a convex body always is) that it makes sense to define an inner product on variations in a way which generalizes the situation with smooth surfaces, and indeed to do other bits of differential geometry - define the shape operator, equate mean curvature with the gradient of the area functional, etc.

I'm sorely lacking in background knowledge, and also currently without university access to journals. Can anyone point me in the direction of some relevant literature or, better yet, help me gain a bit of insight into the theory I'm stumbling towards here?