Abstract. The chapter deals with FE simulations of incipient granular silo flow within enhanced elasto-plasticity. The calculations were performed with a micropolar elasto-plastic constitutive model. The quasi-static, dynamic and rapid flow was studied. A shear zone formation was taken into account in analyses.
Onset of Quasi-static Flow in Model SiloThe FE analysis of the onset of quasi-static mass flow of granular materials in a plane strain model silo with parallel walls with controlled outlet velocity (h = 0.5 m, b = 0.2 m) (Chapter 5.1) was performed using a micro-polar elastoplastic constitutive law (Chapter 4.4). During FE-calculations, the effect of the initial density, modulus of elasticity and mean grain diameter of the solid, initial stress state, wall roughness, wall stiffness, wall flexibility and wall imperfection, and pressure level on the material behaviour was investigated. In addition, the effect of both large deformations and curvatures, and micro-polar constants was analysed.For plane strain calculations, a FE-mesh with quadrilateral finite elements composed of four diagonally crossed triangles was applied to avoid volumetric locking and spurious element behaviour. Totally, 3000 triangular elements (Fig. 7.1) with linear shape functions for the displacements and the Cosserat rotation were used. Symmetry with respect to the centreline was taken into account. In order to describe realistically the interface behaviour along the wall, the FE-mesh was strongly refined at the wall. The width of quadrilateral elements close to the wall was equal to 0.5 mm, 1.5 mm and 3 mm, respectively. The height of all quadrilateral elements was 10 mm.The calculations were performed with both a changing and a constant elastic modulus E calculated by Eq. 4.86. In the first case, one assumed E ≅ 500p (with e 0 = 0.60-83 and C s = 0.003-0.004, Wu 1992)) and in the second case E = 500×4.0=2000 kPa (with p ≅ 4.0 kPa on the basis of initial calculations).As an initial stress state before the bottom lowering, the stresses σ ij after filling according to a slice method by Janssen (1895) for plane strain were assumed (with couple stresses being zero and shear stresses linearly decreasing towards the axis of symmetry). One chose on the basis of model tests: K = 0.29 and ϕ w = 43 o (very 308 7 FE Simulations Based on Enhanced Elasto-Plasticity rough walls, dense sand), K = 0.39 and ϕ w = 39° (very rough walls and medium dense sand), K= 0.39 and ϕ w =32° (very rough walls and initially loose sand), K = 0.27 and ϕ w = 36° (rough walls and dense sand), K = 0.24 and ϕ w = 15º (smooth walls and initially dense sand). The FE-calculations were begun with