A procedure to simplify the generation of basis vectors for shape optimization of complex structures is presented. The approach involves the creation of geometric regions (termed domain elements). Control perturbations are applied to corner and/or midside nodes of the domain. These control perturbations are used together with the domain geometry to automatically generate perturbations for all of the nodes in the region using standard linear and/or quadratic isoparametric interpolation functions. Multiple types of domain elements are presented to facilitate the creation of basis vectors of a wide range of structures.