A superimposing iteration harmonic balance method (SIHB) is presented for the steady-state analysis of strongly nonlinear oscillators.

According to the basic theories of harmony superposition method and linear filtering method,the numerical simulation of wind velocity time history for long-span spatial lattice structure was realized by programming with MATLAB language in such methods as harmony superposition method,FFT-updated harmony superposition method,interpolation method,and auto-regressive method of linear filtering method.

The method of HSM-BTI (harmony superposition method including blade-tower interaction) was proposed in this paper.

According to the numerical simulation methods of wind load in the latest years,the development process of the method is expounded,and the harmony superposition method is introduced briefly,as well as the linear filtering method,the wavelet analysis\'s basic principles in the numerical simulation of wind load and some problems existed.

The wavelet coefficients of the measured wind velocity histories and the simulated data obtained by the conventional harmonic wave superposition method are analyzed by statistic methods such as probability density、kurtosis and skewness analysis.

Procedure of stochastic sequence of wind velocity is generated by wave superposition method, and time domain analysis is performed on the double-layer latticed shell roof structures of Beijing Laoshan Cycling Velodrome under random wind load.