At present, the optimal shape design of structures is of great importance. Problems of this type mostly lead to very complicated nonlinear systems, and therefore new optimization solutions are being sought. One very efficient tool offers inverse variational principles. It starts with formulation of a variational principle under the assumption that the volume (in the 2D area) of the domain variables is constant. It can easily be shown that then the boundary density of potential energy in the optimal state should be constant too. As the differences of the density of the energy at the boundary nodal points may be very large, the new positions of the boundary nodal points are determined in a special scale. Since inverse variational principles, in connection with the boundary element method, seem to be very prospective for solutions of geometric behavior of homogeneous, and also of partly homogeneous (e.g. phase-wise homogeneous), media, this complex study is presented. To prove the ability of the procedure proposed, two examples are solved at the end of this paper. They are selected in such a way that a comparison with statistically determined and statistically undetermined slender beams, obeying Bernoulli-Euler assumptions, is possible. The examples themselves solve an optimization of stretched plates.