Compared with static shape control, dynamic shape control is more difficult due to its time-varying nature. To improve the control precision and reduce the number of actuators needed, an integrated design optimization model of structure and control is proposed. The following three sets of design variables are considered to minimize the time domain variance between the controlled and desired dynamic shapes: the time-varying actuator voltages, the actuator layout and the structural ply parameters. To satisfy the engineering performance requirements to the greatest extent possible, constraints on the structural mass, the energy, the maximum transient voltage and the number of actuators are considered. Because the control voltages are varying in time, a two-level optimization strategy is adopted. In the inner optimization problem, the Newmark integral method is applied to derive discrete time domain expressions for the shape control equations, and then, the Kuhn-Tucker condition is introduced to calculate the optimal time-varying voltage distribution. To solve the outer optimization problem, because of the coexistence of discrete variables and continuous variables, a simulated annealing algorithm is used. To address the shape control of complicated curved shells, a finite element formulation for an eight-node laminated curved shell element with piezoelectric actuators is derived. Numerical examples show that the proposed integrated design optimization method can significantly improve the control effect and that the optimization of the structural ply parameters plays an important role. Moreover, the control system can be simplified by taking the minimum number of actuators as the objective function when the control accuracy allows.