The problem of topological optimization of a symmetrical double-shear adhesive joint has been solved. The suggested mathematical model of a joint with variable thickness generalizes the classic Holland–Reissner model. The shape of the doubler is described by means of the Bezier curve. Seeking parameters in the optimization problem are coordinates of reference points of the Bezier curve. Both joint length and doubler cross-section area can be considered an objective function. The restriction is applied on stress in adhesive film and doubler. The direct problem of finding the joint stress state at given geometric parameters was solved using the finite difference method. The genetic algorithm was used to solve the optimization problem. In order to improve the convergence of the genetic algorithm, the island model of evolution is suggested, which ensures quick evolution selection and stability of obtained results. The model problem is solved.