The reliability based design optimization of a long-span suspension bridge is performed considering probabilistic flutter constraint accounting for the stochastic nature of parameters such as the extreme wind velocity and experimentally obtained flutter derivatives. The objective of this study is to seek the minimum volumes of the main cables and bridge girder by varying the main cable area and girder plate thicknesses while satisfying the required safety level under flutter. Three methods of Reliability Index Approach, Performance Measure Approach and Sequential Optimization and Reliability Assessment are utilized to solve this problem. The main code to perform the entire RBDO process is written in Matlab (R), which is linked to Abaqus finite element models to obtain structural responses and FLAS code developed by our research group to compute aeroelastic responses of the bridge. The three RBDO methods were applied to the Messina Bridge project to examine their numerical performance and computational efficiency. (C) 2015 Elsevier Ltd. All rights reserved.