Nearly zero energy buildings (nZEBs) are considered as a promising solution to mitigate the energy and environmental problems. A proper sizing of the nZEB systems (e.g. HVAC systems, PV panels, wind turbines and batteries) is essential for achieving the desirable level of thermal comfort, energy balance and grid dependence. Parameter uncertainty, component degradation and maintenance are three crucial factors affecting the nZEB system performances and should be systematically considered in system sizing. Until now, there are some uncertainty-based design methods been developed, but most of the existing studies neglect component degradation and maintenance. Due to the complex impacts of degradation and maintenance, proper sizing of nZEB systems considering multiple criteria (i.e. thermal comfort, energy balance and grid dependence) is still a great challenge. This paper, therefore, proposes a robust design method of nZEB systems using genetic algorithm (GA) which takes into account the parameter uncertainty, component degradation and maintenance. The nZEB life-cycle cost is used as the fitness function, and the user’ performance requirements on thermal comfort, energy balance and grid dependence are defined as three constraints. This study can help improve the designers’ understanding of the impacts of uncertainty, degradation, and maintenance on the nZEB life-cycle performances. The proposed method is effective in minimizing the nZEB life-cycle cost through designing the robust optimal nZEB systems sizes and planning the optimal maintenance scheme, meanwhile satisfying the user specified constraints on thermal comfort, energy balance, and grid dependence during the whole service life.