A three-dimensional dynamic model of vehicle-road-unsaturated subgrade coupling is established. The semi-analytical and numerical solutions of three-dimensional dynamic response of highway subgrade under vehicle load are presented for the first time. Then, the accuracy of the theoretical model is verified by the field measurement data. Finally, the spatial distribution of dynamic stress and depth of working zone under different influencing factors are studied, and the calculation model of depth of working zone of highway subgrade is proposed. The results show that the dynamic stresses, σx, σy, σz, and τxz increase significantly with the increase of vehicle axle load and decrease significantly with the increase of road thickness under the same subgrade depth. In addition, with the increase of axle load and the decrease of road thickness, the peak values of dynamic stress σx, σy, σz, and τxz decrease more obviously along the depth of roadbed. In the process of vehicle movement, the traditional subgrade working area can be divided into a sensitive area and an influence area of dynamic load. The depth of these two areas is positively correlated with the axle load of the vehicle and negatively correlated with the thickness of the road surface. The advantage of the model in this paper is that the semi-analytical solution of the total stress component of the soil element in the course of traffic load is given.