The paper presents an approach to nonlinear buckling fiber angle optimization of laminated composite shell structures. The approach accounts for the geometrically nonlinear behaviour of the structure by utilizing response analysis up until the critical point. Sensitivity information is obtained efficiently by an estimated critical load factor at a precritical state. In the optimization formulation, which is formulated as a mathematical programming problem and solved using gradient-based techniques, a number of the lowest buckling factors is included such that the risk of "mode switching" during optimization is avoided. The presented optimization formulation is compared to the traditional linear buckling formulation and two numerical examples, including a large laminated composite wind turbine main spar, clearly illustrate the pitfalls of the traditional formulation and the advantage and potential of the presented approach.