The initial geo-stress field is an important initial condition for finite element analysis of structure-soil interaction. The accuracy of initial geo-stress balance results directly affects the accuracy of numerical analysis results of structure. The operation of the automatic geo-stress balance method is the simplest，and the results of it are greatly better than other calculation methods. However，owing to being confined seriously by the foundation and structure type，the soil foundation parameters and contact conditions，the application of it in practice is very limited. We adopt the software of ABAQUS based on following the actual timeline for geo-stress balance in the process of numerical simulation strictly，apply the way of setting life and death elements to implement the replacement of the foundation-structure space element before and after the initial geo-stress balance，establish a general automatic geo-stress balance method and its procedures，and adopt a series of typical complex soil-wharf structure calculation models to verify the proposed method. The results show that the automatic geo-stress balance method and its procedures established in this paper are easy to implement but with the highest accuracy compared with other methods，and it can conveniently and effectively achieve the initial geo-stress balance for various complex calculation models.