This contribution focuses on the combined analysis of mistuning and unilateral blade-tip/casing contacts. A two-dimensional (2D) phenomenological finite element (FE) model of an aircraft engine fan stage is considered. It is reduced by means of the Craig–Bampton component mode synthesis method and contact treatment relies on a Lagrange multiplier algorithm within an explicit time-integration scheme. Blade-tip/casing contacts are initiated through the deformed shape of a perfectly rigid casing. Mistuning is accounted for on the blades only. Monte Carlo simulations are carried out in both linear and nonlinear configurations, which allows to compare amplifications predicted in both context due to mistuning. Following a thorough convergence analysis of the proposed numerical strategy, the influence of mistuning level as well as the configuration of the external forcing are investigated. Presented results underline the detrimental consequences of mistuning in a nonlinear structural context, yielding even higher vibration amplifications than in a linear context. A cross-analysis between linear and nonlinear computations reveals that no correlation is found between linear and nonlinear amplifications which suggests that the effect of existing strategies to mitigate vibration amplifications within a linear context may not be suitable within a nonlinear context.