Carbon fiber–reinforced epoxy composites are widely used in the marine and aviation sectors due to their superior mechanical properties. As research into these materials progresses, the demand for materials with even more complex and superior properties has increased. Traditional trial-and-error methods are limited to optimizing single factors and are plagued by heavy workloads and prohibitive costs when applied to multifactor composite designs. Addressing these problems, this study introduces a nonparametric multiobjective optimization process for composite materials. This approach employs a Gaussian process (GP) to establish a multiobjective regression model and uses nondominated sorting genetic algorithm II (NSGA II) to determine the Pareto front. Pareto solutions were ranked in accordance with their similarity to ideal solutions and on the basis of predetermined preference weights. The results demonstrate variations in optimal outcomes contingent on the weight values assigned. The differences between the optimized results and experimental validations reached 12.9%, with the minimum deviation being 0.6%. The effectiveness of the proposed method was demonstrated. The strategy combines GP and NSGA II to extend design methods for small samples and materials with multiple attributes. This method enables the simultaneous optimization of multiple conflicting objectives, offering greater efficiency than addressing them individually. •A nonparametric model was established to simultaneously optimize multiple properties.•Combinatorial weighting was introduced to achieve balanced Pareto-optimal solutions with superior performance.•Pareto-optimal solutions are determined considering various problem cases.•The properties of the carbon fiber-reinforced composite exhibited high consistency between experiments and predictions.