This study focuses on the aero-elastic-wake coupling behavior of wind turbines. A newly developed code called FALM (Flexible Actuator Line Model), which combines the wake simulation ability of the actuator line mothed (ALM) and the structural simulation ability of the flexible multibody dynamics method, was employed to achieve these simulations. The power output and thrust, the natural frequency and deformation and the wake characteristics of a single NREL 5MW wind turbine were studied to validate this code. It shows that nonlinear effects such as spin softening and stress stiffening were fully considered by FALM and it can also guarantee a reliable prediction of power and thrust of wind turbines. Furthermore, a case of two-wind-turbines with the inlet wind speed of 14m/s (exceeding the rated wind speed) were carried out to study the aero-elastic-wake coupling behavior in a wind farm. It shows that FALM is able to simulate multiple wind turbines with power control system involved. The pitch control process of the upstream wind turbine was predicted and the dynamic loads of the downstream wind turbine caused by the wake effect were studied. These results will contribute to the study of reducing the fatigue load caused by wake effect.