Eigenfrequency optimization for 3D continua is formulated and exemplified by the geometry and boundary conditions of a thick plate. Numerical finite element models are based on four node tetrahedra and results from subspace iterations give directly the basis for the continuum redesign. The 3D modeling with a large number of elements has the possibility in optimal design to obtain (as found) not only holes but also cavities in the continuum. Sensitivity analysis is presented on the element level with simple physical interpretation of the involved terms. This general result has general value for control of eigenfrequencies. It is found that in the combination of partial differentiation with the chain rule of differentiation, a specific notation is needed and a suggestion is presented. The optimization method is based on a derived optimality criterion, and as such the maximization problem change to a problem of determining a design with uniform values of this criterion. Nonlinear stiffness interpolation may be a physical reality. A two parameter interpolation function is incorporated analytical, also in the sensitivity analysis and the optimality criterion, but without focusing on 1–0 optimal solutions. Two cases of boundary conditions, two cases of total amount of material, and cases of linear and nonlinear stiffness interpolation are studied.