Machine learning-based aerodynamic reduced-order models (ROMs) combine high accuracy with extremely low computational costs, making them highly effective in predicting nonlinear and unsteady bridge aerodynamic forces. Although several machine learning-based nonlinear aerodynamic models have been developed, the majority are built on a single wind speed parameter. However, in nonlinear aerodynamic prediction and aeroelastic analysis of bridges, the variability in incoming wind speed significantly influences the computed results. A ROM relying solely on a single wind speed lacks the ability to accurately forecast the intricate dynamic behaviors arising from changes in wind speed. When the incoming wind speed changes, the model's prediction accuracy significantly decreases. Usually, it is necessary to establish a new database and train a new model, which not only increases time and cost but also greatly reduces the convenience of the ROM. Addressing this challenge, this study proposes a multiple-wind-speed (MWS) nonlinear unsteady bridge aerodynamic model based on the LSTM deep neural network. Taking the Taohuayu Yellow River Bridge in the Henan Province of China as an example, the modeling process of the proposed MWS-ROM is demonstrated, along with non-linear aerodynamic predictions of the deck under various conditions and aerodynamic-elastic analysis of the deck under different wind speeds. The research results show that the trained LSTM network can accurately predict the nonlinear aerodynamic forces of bridges under single and double degrees of freedom vibration conditions. The MWS-ROM performed well in predicting convergent vibrations at low wind speeds and limits cycle oscillations at high wind speeds, aligning closely with results from the CFD full-order model. Compared to CFD, the aerodynamic ROM based on the LSTM network significantly enhances computational efficiency, consequently boosting the convenience and efficiency of bridge flutter analysis. Additionally, the methodology proposed herein can be extended for wind-induced vibration control and response prediction in other types of deck sections.