Currently developed underwater gliders can be roughly divided into the two types:traditional configuration and unconventional configuration. As a type of underwater gliders with unconventional configuration, a blended-wing-body (BWB) underwater glider has better fluid dynamic performances because of its unique shape. However, it is difficult to design the shape of the BWB underwater glider that has excellent hydrodynamic performances. Therefore, it is of great significance to optimize its shape, which this paper carries out by using the free-form deformation (FFD). The complete and automatic shape optimization framework is established by jointly using FFD parameterization method, CFD solver, optimization algorithm and mesh deformation method. The framework is used to optimize the shape of a BWB underwater glider. The average drag coefficient of the BWB underwater glider during its sinking and floating in one working period is used as the objective function to optimize its shape, with the volume constraints considered. The optimization results show that the gliding performance of the BWB underwater glider is remarkably enhanced.