Abstract: In order to improve the accuracy of burden surface prediction in a blast furnace with bell-less top, a mathematical model of burden distribution considering burden motion was established. Based on the analysis of burden motion in the furnace, the burden motion was pointed as an important factor affecting the burden formation, and a 1:10 dimension scale bell-less top distributor model was used to analyze the effect of different burden components' velocity on heap formation. Then a heap formation prediction model considering burden motion was established, the numerical method was employed to determine the heap position and the heap shape, and at last the model was applied to predict the burden surface. Results show that the burden motion is the main reason causing the difference between the heaping angles in both sides of the heap, the variations of area of the heap cross section and heap profile. The heap profile is constructed by lines and curves in the heap formation prediction model, and the important model parameters, including the heaping angles and the length of curved transition area, all take burden motion into account. Thus the profile constructed is close to the actual heap, and the accurate predicted burden surface is achieved after applying this model to burden distribution prediction.