Reliability optimization design of the gear box is proposed. It includes an objective function, 30 design variables and 52 constraints. It is important to note that material properties, geometry parameters and applied loads of the structure are assumed to be normal random variables. Reliability calculation adopts the HL-RF method. The comparison of design parameters demonstrates the proposed method.