Stress distribution in two-dimensional silos

Simulations of a polydispersed two-dimensional silo were performed using molecular dynamics, with different numbers of grains reaching up to 64 000, verifying numerically the model derived by Janssen and also the main assumption that the walls carry part of the weight due to the static friction between grains with themselves and those with the silo's walls. We vary the friction coefficient, the radii dispersity, the silo width, and the size of grains. We find that the Janssen's model becomes less relevant as the the silo width increases since the behavior of the stresses becomes more hydrostatic. Likewise, we get the normal and tangential stress distribution on the walls evidencing the existence of points of maximum stress. We also obtained the stress matrix with which we observe zones of concentration of load, located always at a height around two thirds of the granular columns. Finally, we observe that the size of the grains affects the distribution of stresses, increasing the weight on the bottom and reducing the normal stress on the walls, as the grains are made smaller (for the same total mass of the granulate), giving again a more hydrostatic and therefore less Janssen-type behavior for the weight of the column.