For the efficient reliability analysis, Bi-direction two-point approximation(BTPA) method is developed which solves shortcomings of conventional two-point approximation(TPA) methods that generate an approximate surface with low accuracy or sometimes do an unstable approximate surface. The conventional reliability based design optimization(RBDO) methods require high computational cost compared with the deterministic design optimization(DO) methods. To overcome the computational inefficiency of RBDO, the approximate reliability analysis approaches on the TPA surface are proposed. Using these FORM and SORM analysis strategies, multi-point aerodynamic-structure interacted shape design optimizations with uncertainty are performed very efficiently.