In view of the problem that the output thrust of solenoid actuator is affected by various factors and is difficult to calculate in actual working conditions, this paper proposes a semi-analytical model constructed by magnetic field subdomain method with internal and external boundary conditions in cylindrical coordinate system for calculation, and the general solution equations of magnetic vector potential for each subdomain are derived and solved by MATLAB. Taking a push-pull electromagnet as an example, the finite element simulation and experimental comparative analysis are carried out. The correctness and applicable conditions of the subdomain method are illustrated by comparing the gradient plot of magnetic vector potential, inductance curve and electromagnetic force. By calculating the output thrust after considering the core gravity, spring force and electromagnetic force, it is shown that this method has the advantage of computational flexibility compared with the finite element method, and it is easier to write special algorithms according to various working conditions to calculate the important parameters in engineering applications.