This subroutine models the dynamic load from four wheels passing over a bridge deck. The load F is applied based on the wheel's position, which moves over time with a fixed increment of 0.01s. The force is applied when the vertical position Z is within specific zones (2.2-2.5 m and 3.5-3.8 m) and the horizontal position X matches the calculated displacement. If the conditions are met, a force of 1 is applied; otherwise, the force is zero.

The corresponding Abaqus .inp file and subroutine file are attached in a ZIP archive for your convenience.