In each vibration source of shunt reactor, the winding vibration caused by Lorentz force should not be ignored. In this paper, the electromagnetic-mechanical coupling field simulation model was established based on the structure of disks and blocks of shunt reactor, and the simulation modeling method was verified by an actual reactor manufactured according to the scaled model. The vibration caused by the alternating Lorentz force acting on the winding was simulated and analyzed. The results show that the axial Lorentz force is symmetrical up and down, whose maximum values are distributed at the top and bottom of the winding, and the values approaches zero at about 1/6, 1/2, 5/6 height. The winding deformation and amplitude are mainly axial, and the overall distribution along the winding axis is M-shape. With the increase of the elastic modulus or the number of winding blocks, the axial amplitudes of disks show a downward trend, but the decline rate gradually slows down. In addition to vibration attenuation performance, the design of winding blocks should also consider the investment costs and thermal constraints. The research results could provide a basis for simulation and control of reactor winding vibration.