Abstract: The whole residual stress field of structural components can be reconstructed based on measured results from a limited number of points by the stress function method, but stress function is usually difficult to determine and it is not easy to reconstruct the distribution of several stress components through one stress function. To overcome these shortcomings, a method was developed to simultaneously reconstruct several residual stress component fields by using a general stress function. In this method, the binary Fourier function was used as the Airy stress function. The order of the model was determined by the determination coefficient, and the unknown parameters in the stress function were optimized according to measured results on some points. Finally, the accuracy of this method was verified by simulated residual stress fields of a three-point bending beam, and it was applied to the reconstruction of residual stress field of a real aluminum alloy plate. The results show that this method can simultaneously reconstruct several stress components with high accuracy.