  We consider the eigenvalues of the biharmonic operator subject to several
homogeneous boundary conditions (Dirichlet, Neumann, Navier, Steklov). We show
that simple eigenvalues and elementary symmetric functions of multiple
eigenvalues are real analytic, and provide Hadamard-type formulas for the
corresponding shape derivatives. After recalling the known results in shape
optimization, we prove that balls are always critical domains under volume
constraint.
