In this work, a numerical study of the gas-solid flow in a gas cyclone is carried out by use of the combined discrete element method (DEM) and computational fluid dynamics (CFD) model where the motion of discrete particles phase is obtained by DEM which applies Newton s equations of motion to every individual particle and the flow of continuum fluid by the traditional CFD which solves the Navier-Stokes equations at a computational cell scale. The model successfully captures the key flow features in a gas cyclone, such as the strands flow pattern of particles, and the decrease of pressure drop and tangential velocity after loading solids. The effect of solid loading ratio is studied and analysed in terms of gas and solid flow structures, and the particle-gas, particle-particle and particle-wall interaction forces. It is found that the gas pressure drop increases first and then decreases when solids are loaded. The reaction force of particles on gas flow is mainly in the tangential direction and directs mainly upward in the axial direction. The reaction force in the tangential direction will decelerate gas phase and the upward axial force will prevent gas phase from flowing downward in the near wall region. The intensive particle-wall collision regions mainly locate in the wall opposite to the cyclone inlet and the cone wall. Moreover, as the solid loading ratio increases, number of turns travelled by solids in a cyclone decreases especially in the apex region of the cyclone while the width of solid strands increases, the pressure drop and tangential velocity decrease, the high axial velocity region moves upwards, and the radial flow of gas phase is significantly dampened.