Gusts have a significant impact on aircraft and need to be analyzed through flight simulations. The solution for time-domain gust aerodynamic forces stands as a pivotal stage in this process. With the increasing demand for flight simulations within gusty environments, traditional methods related to gust aerodynamics cannot fail to balance computational accuracy and efficiency. A method that can be used to quickly and accurately calculate the time-domain gust aerodynamic force is needed. This study proposes the fitting strip method, a gust aerodynamic force solution method that is suitable for real-time flight simulations. It only requires the current and previous gust information to calculate the aerodynamic force and is suitable for different configurations of aircraft and different kinds of gusts. Firstly, the fitting strip method requires the division of fitting strips and the calculation of the aerodynamic force under calibration conditions. In this study, the double-lattice method and computational fluid dynamics are used to calculate the aerodynamic force of the strips. Then, the amplitude coefficients and time-delay coefficients are obtained through a fitting calculation. Finally, the coefficients and gust information are put into the formula to calculate the gust aerodynamic force. An example of a swept wing is used for validation, demonstrating congruence between the computational results and experimental data across subsonic and transonic speeds, which proves the accuracy of the fitting strip method in both discrete gusts and continuous gusts. Compared with other methods, the fitting strip method uses the shortest time. Furthermore, the results of a calculation for normal-layout aircraft show that this method avoids the shortcomings of the rational function approximation method and is more accurate than the gust grouping method. Concurrently, gust aerodynamic force calculations were performed on aircraft with large aspect ratios and used in a real-time flight simulation.