Large wind turbine system is a periodic time-varying system with rigid-flexible coupling multi-bodies. The traditional finite element method cannot solve the singular stiffness matrix produced by the rotation of blades. However, the vector form intrinsic finite element method can effectively solve the geometric deformation of elastic continuum, the nonlinear or discrete constitutive model, the coupling motion of continuum and rigid body, and so on. In this study, a solver program of space beam elements was developed by the vector form intrinsic finite element method with MATLAB code, and verified by two typical examples where one is a cantilever method with a dynamic force acting at the end, and the other is an Euler beam rotating around a fixed axis. The integrated simulation of wind turbine system that consists of tower, rotor blades, and nacelle was then established, and its dynamic response of free vibration under parking was analyzed. The natural frequencies of the turbine system were obtained by modal parameter identification, and they agree well with the results obtained by traditional finite element method. Moreover, the weighted amplitude wave superposition method and proper orthogonal decomposition method as well as B-spline surface interpolation were employed to obtain the wind time series of wind turbine under the normal operation condition. The wind-induced dynamic response of wind turbine system was also calculated by vector form intrinsic finite element method. The results reflect the periodic influence of gravity on the internal forces of blades and the interaction between blades and the tower.