Wind loads are key considerations in the structural design of steel roof structures, especially for large span ones. The analysis of wind loads on large span steel roof structure (LSSRS) requires large amounts of calculations. Due to combined effects of horizontal and vertical winds, the wind induced vibrations of LSSRS are analyzed with the frequency domain method as the first application of the method for the analysis of wind responses of LSSRS. A program is developed to analyze the wind-induced vibrations due to a combination of wind vibration modes. The program, which predicts the wind vibration coefficient and wind pressure acting on the LSSRS, is designed with input and output interfaces to other finite element software, resulting in preferably solving the wind load analytical problem in the design of LSSRS. The effectiveness and accuracy of the proposed method and the program are verified by numerical analyses of practical projects.