The layup optimization by genetic algorithm (GA) for the interlaminar strength of laminated composites with free edge is presented. For the calculation of interlaminar stresses of composite laminates with free edges, extended Kantorovich method is applied. In the formulation of GA, repair strategy is adopted for the satisfaction of given constraints. In order to consider the bounded uncertainty of material properties, convex modeling is used. Results of GA optimization with scattered properties are compared with those of optimization with nominal properties. The GA combined with convex modeling can work as a practical tool for maximum interlaminar strength design of laminated composite structures, since uncertainties are always encountered in composite materials and the optimal results can be changed.