In order to obtain the dynamic response of the tubular tower of the horizontal axis wind turbine subjected to the time-varying load, the natural dynamic characteristics of the bending vibration of the tower was analyzed by considering the tower as a cantilever and applying two-node beam element for discretization modeling. Through establishing the structure dynamic equation and calculating the time-varying load on the tower, the dynamic response of the tower was calculated by applying linear acceleration method and modal superposition principle. A correspongding FEM program was developed to carry out the numerical analysis. The dynamic response of the tower of a 1.0 MW wind turbine operating in turbulent wind was simulated by applying the developed program. The obtained displacement, software velocity and acceleration in the simulation time were compared with the results calculated by the "GH Bladed" software, and good agreements were achieved. This verifies the feasibility of the proposed model.