A design platform for airship shape optimization with propeller is established with computational fluid dynamics and optimization techniques. Firstly, Hicks-Henne shape functions were used to parameterize hull of the airship. Then,a hybrid dynamic mesh generation method was presented based on inverse distance weighting (IDW) interpolation and transfinite interpolation (TFI). An actuator disk model of propeller was employed to simulate viscous flow around the airship. LOTTE airship was chosen as the initial shape, and drag coefficient was chosen as the objective function. Airship volume was selected as a constraint. The optimized shape based on non-linear programming by quadratic lagrangian (NLPQL) algorithm was obtained. It shows that afterbody shape of the optimized airship is relatively flat and location of the maximum diameter moves forward. Aerodynamic performance of the new shape is significantly improved. Compared with that of the LOTTE hull,drag coefficient of the optimized shape is reduced by 26.88%.