The fluid seepage in saturated zone of subgrade promotes the migration of fine particles in the filler, resulting in the change of pore structure and morphology of the filler and the deformation of solid skeleton, which affects the fluid seepage characteristics. Repeatedly, the muddy interlayer, mud pumping, and other diseases are finally formed. Based on the theory of two-phase seepage, the theory of porous media seepage, and the principle of effective stress in porous media, a two-phase fluid-solid coupling mathematical model in saturated zone of subgrade considering the effects of fine particles migration is established. The mathematical model is numerically calculated with the software COMSOL Multiphysics^®. The two-phase seepage characteristics and the deformation characteristics of the solid skeleton in saturated zone of the subgrade are studied. The research results show that the volume fraction of fine particles first increases then decreases and finally becomes stable with the increase of time, due to the continuous erosion and migration of fine particles in saturated zone of the subgrade. The volume fraction of fine particles for the upper part of the subgrade is larger than that for the lower part of the subgrade. The porosity, the velocity of fluid, the velocity of fine particles, and the permeability show a trend of increasing first and then stabilizing with time; the pore water pressure has no significant changes with time. The vertical displacement increases first and then decreases slightly with the increase of time, and finally tends to be stable. For the filler with a larger initial volume fraction of fine particles, the maximum value of the volume fraction of fine particles caused by fluid seepage is larger, and the time required to reach the maximum value is shorter. It can be concluded that the volume fraction of fine particles in the subgrade filler should be minimized on the premise that the filler gradation meets the requirements of the specification in actual engineering.

Leakage is a common defect of embankment dam slopes, and is usually accompanied by internal erosion. This study establishes a mathematical two-phase seepage coupling model that explicitly considers the effects of internal erosion. With the help of finite element software COMSOL Multiphysics^®, we use this model to study the characteristics of fluid seepage, the fraction of leakage volume that is made up of migrated fine particles (migrated fine particles volume fraction), porosity, permeability, and displacement in the leakage process of a three-dimensional embankment dam slope, as well as to study the influence of the water level and leakage outlet size on the aforementioned features. Our results show that water seepage velocity gradually increases with time, especially at the downstream leakage outlet. Therefore, the erosion and migration of fine particles occur primarily at the downstream leakage outlet, resulting in a significant increase in the migrated fine particles volume fraction, porosity, and permeability. In addition, the maximum migrated fine particles volume fraction, the maximum porosity, and the maximum permeability in the slope increase nonlinearly with time. Curves for maximum displacement with a strength reduction factor can be divided into three stages: the early stable stage, the midterm nonlinear growth stage (creeping state), and the later rapid growth stage (instability state). We also find that the rise of the water level promotes the erosion and migration of fine particles on the slope and that the maximum migrated fine particles volume fraction and the maximum porosity increase with the water level. The values of the strength reduction factor for the creeping state decrease as the water level rises. Furthermore, the larger the size of the downstream leakage outlet, the larger the maximum migrated fine particles volume fraction and maximum porosity. However, the relationship between maximum displacement and strength reduction factor is nearly unaffected by the size of the leakage outlet.

The fluid seepage in local-saturated zone of subgrade promotes the migration of fine particles in the filler, resulting in the change of pore structure and morphology of the filler and the deformation of solid skeleton, which affects the fluid seepage characteristics. Repeatedly, the muddy interlayer, mud pumping and other diseases are finally formed. Based on the theory of two-phase seepage, the theory of porous media seepage, and the principle of effective stress in porous media, a two-phase fluid-solid coupling mathematical model in local-saturated zone of subgrade considering the effect of fine particles migration is established. The mathematical model is numerically calculated with the software COMSOL Multiphysics○R, the two-phase seepage characteristics and the deformation characteristics of the solid skeleton in local-saturated zone of the subgrade are studied. The research results show that due to the continuous erosion and migration of fine particles in local-saturated zone of the subgrade, the volume fraction of fine particles first increases then decreases and finally becomes stable with the increase of time. And the volume fraction of fine particles for the upper part of the subgrade is larger than that for the lower part of the subgrade. The porosity, the velocity of fluid, the velocity of fine particles, and the permeability show a trend of increasing first and then stabilizing with time; the pore water pressure has no significant changes with time. The vertical displacement increase first and then decrease slightly with the increase of time, and finally tend to be stable. For a filler with a larger initial volume fraction of fine particles, the maximum value of the volume fraction of fine particles caused by fluid seepage is larger, and the time required to reach the maximum value is shorter. It can be concluded that in actual engineering, the volume fraction of fine particles in the subgrade filler should be minimized on the premise that the filler gradation meets the requirements of the specification.