The paper presents a formulation for multidisciplinary design optimization of vessels, subject to uncertain operating conditions. The formulation couples the multidisciplinary design analysis with the Bayesian approach to decision problems affected by uncertainty. In the present context, the design specifications are no longer given in terms of a single operating design point, but in terms of probability density function of the operating scenario. The optimal configuration is that which maximizes the performance expectation over the uncertain parameters variation. In this sense, the optimal solution is "robust" within the stochastic scenario assumed. Theoretical and numerical issues are addressed and numerical results in the hydroelastic optimization of a keel fin of a sailing yacht are presented.