Frictionless bead packs have macroscopic friction, but no dilatancy

The statement of the title is shown by numerical simulation of homogeneously sheared packings of frictionless, nearly rigid beads in the quasistatic limit. Results coincide for steady flows at constant shear rate γ in the limit of small γ and static approaches, in which packings are equilibrated under growing deviator stresses. The internal friction angle ϕ, equal to 5.76 $\pm$ 0.22 degrees in simple shear, is independent on the average pressure P in the rigid limit. It is shown to stem from the ability of stable frictionless contact networks to form stress-induced anisotropic fabrics. No enduring strain localization is observed. Dissipation at the macroscopic level results from repeated network rearrangements, like the effective friction of a frictionless slider on a bumpy surface. Solid fraction Φ remains equal to the random close packing value ≃ 0.64 in slowly or statically sheared systems. Fluctuations of stresses and volume are observed to regress in the large system limit, and we conclude that the same friction law for simple shear applies in the large psystem limit if normal stress or density is externally controlled. Defining the inertia number as I = γ m/(aP), with m the grain mass and a its diameter, both internal friction coefficient $\mu$∗ = tan ϕ and volume 1/Φ increase as powers of I in the quasistatic limit of vanishing I, in which all mechanical properties are determined by contact network geometry. The microstructure of the sheared material is characterized with a suitable parametrization of the fabric tensor and measurements of connectivity and coordination numbers associated with contacts and near neighbors.Comment: 19 pages. Additional technical details may be found in v

Crack channelling mechanisms in brittle coating systems under moisture or temperature gradients

Abstract: Crack channelling is predicted in a brittle coating-substrate system that is subjected to a moisture or temperature gradient in the thickness direction. Competing failure scenarios are identified, and are distinguished by the degree to which the coating-substrate interface delaminates, and whether this delamination is finite or unlimited in nature. Failure mechanism maps are constructed, and illustrate the sensitivity of the active crack channelling mechanism and associated channelling stress to the ratio of coating toughness to interfacial toughness, to the mismatch in elastic modulus and to the mismatch in coefficient of hygral or thermal expansion. The effect of the ratio of coating to substrate thickness upon the failure mechanism and channelling stress is also explored. Closed-form expressions for the steady-state delamination stress are derived, and are used to determine the transition value of moisture state that leads to unlimited delamination. Although the results are applicable to coating-substrate systems in a wide range of applications, the study focusses on the prediction of cracking in historical paintings due to indoor climate fluctuations, with the objective of helping museums developing strategies for the preservation of art objects. For this specific application, crack channelling with delamination needs to be avoided under all circumstances, as it may induce flaking of paint material. In historical paintings, the substrate thickness is typically more than ten times larger than the thickness of the paint layer; for such a system, the failure maps constructed from the numerical simulations indicate that paint delamination is absent if the delamination toughness is larger than approximately half of the mode I toughness of the paint layer. Further, the transition between crack channelling with and without delamination appears to be relatively insensitive to the mismatch in the elastic modulus of the substrate and paint layer. The failure maps developed in this work may provide a useful tool for museum conservators to identify the allowable indoor humidity and temperature fluctuations for which crack channelling with delamination is prevented in historical paintings

Crack channelling mechanisms in brittle coating systems under moisture or temperature gradients

Abstract: Crack channelling is predicted in a brittle coating-substrate system that is subjected to a moisture or temperature gradient in the thickness direction. Competing failure scenarios are identified, and are distinguished by the degree to which the coating-substrate interface delaminates, and whether this delamination is finite or unlimited in nature. Failure mechanism maps are constructed, and illustrate the sensitivity of the active crack channelling mechanism and associated channelling stress to the ratio of coating toughness to interfacial toughness, to the mismatch in elastic modulus and to the mismatch in coefficient of hygral or thermal expansion. The effect of the ratio of coating to substrate thickness upon the failure mechanism and channelling stress is also explored. Closed-form expressions for the steady-state delamination stress are derived, and are used to determine the transition value of moisture state that leads to unlimited delamination. Although the results are applicable to coating-substrate systems in a wide range of applications, the study focusses on the prediction of cracking in historical paintings due to indoor climate fluctuations, with the objective of helping museums developing strategies for the preservation of art objects. For this specific application, crack channelling with delamination needs to be avoided under all circumstances, as it may induce flaking of paint material. In historical paintings, the substrate thickness is typically more than ten times larger than the thickness of the paint layer; for such a system, the failure maps constructed from the numerical simulations indicate that paint delamination is absent if the delamination toughness is larger than approximately half of the mode I toughness of the paint layer. Further, the transition between crack channelling with and without delamination appears to be relatively insensitive to the mismatch in the elastic modulus of the substrate and paint layer. The failure maps developed in this work may provide a useful tool for museum conservators to identify the allowable indoor humidity and temperature fluctuations for which crack channelling with delamination is prevented in historical paintings

Solid behavior of anisotropic rigid frictionless bead assemblies

We investigate the structure and mechanical behavior of assemblies of frictionless, nearly rigid equal-sized beads, in the quasistatic limit, by numerical simulation. Three different loading paths are explored: triaxial compression, triaxial extension and simple shear. Generalizing recent results [1], we show that the material, despite rather strong finite sample size effects, is able to sustain a finite deviator stress in the macroscopic limit, along all three paths, without dilatancy. The shape of the yield surface is adequately described by a Lade-Duncan (rather than Mohr-Coulomb) criterion. While scalar state variables keep the same values as in isotropic systems, fabric and force anisotropies are each characterized by one parameter and are in one-to-one correspondence with principal stress ratio along all three loading paths.The anisotropy of the pair correlation function extends to a distance between bead surfaces on the order of 10% of the diameter. The tensor of elastic moduli is shown to possess a nearly singular, uniaxial structure related to stress anisotropy. Possible stress-strain relations in monotonic loading paths are also discussed

Discrete element modelling of scaled railway ballast under triaxial conditions

The aim of this study is to demonstrate the use of tetrahedral clumps to model scaled railway ballast using the discrete element method (DEM). In experimental triaxial tests, the peak friction angles for scaled ballast are less sensitive to the confining pressure when compared to full-sized ballast. This is presumed to be due to the size effect on particle strength, whereby smaller particles are statistically stronger and exhibit less abrasion. To investigate this in DEM, the ballast is modelled using clumps with breakable asperities to produce the correct volumetric deformation. The effects of the quantity and properties of these asperities are investigated, and it is shown that the strength affects the macroscopic shear strength at both high and low confining pressures, while the effects of the number of asperities diminishes with increasing confining pressure due to asperity breakage. It is also shown that changing the number of asperities only affects the peak friction angle but not the ultimate friction angle by comparing the angles of repose of samples with different numbers of asperities

Combined Discrete-Continuum Analysis for Ballasted Rail Tracks

A study on the load-deformation behaviour of railway ballast aggregates subjected to cyclic loadings using a combined discrete-continuum modelling approach is presented. Discrete ballast particles are simulated in the DEM and the continuum-based subgrade is simulated by the FDM. Interface elements are generated to transmit contact forces and displacements between the two domains (i.e. discrete and continuum) whereby the DEM exchanges contact forces to the FDM, and then the FDM transfers the displacement back to the DEM. Distributions of contact forces, coordination number, stress contours on the subgrade and corresponding number of broken bonds (representing ballast breakage) are analysed

Railway-induced ground vibrations – a review of vehicle effects

This paper is a review of the effect of vehicle characteristics on ground- and track borne-vibrations from railways. It combines traditional theory with modern thinking and uses a range of numerical analysis and experimental results to provide a broad analysis of the subject area. First, the effect of different train types on vibration propagation is investigated. Then, despite not being the focus of this work, numerical approaches to vibration propagation modelling within the track and soil are briefly touched upon. Next an in-depth discussion is presented related to the evolution of numerical models, with analysis of the suitability of various modelling approaches for analysing vehicle effects. The differences between quasi-static and dynamic characteristics are also discussed with insights into defects such as wheel/rail irregularities. Additionally, as an appendix, a modest database of train types are presented along with detailed information related to their physical attributes. It is hoped that this information may provide assistance to future researchers attempting to simulate railway vehicle vibrations. It is concluded that train type and the contact conditions at the wheel/rail interface can be influential in the generation of vibration. Therefore, where possible, when using numerical approach, the vehicle should be modelled in detail. Additionally, it was found that there are a wide variety of modelling approaches capable of simulating train types effects. If non-linear behaviour needs to be included in the model, then time domain simulations are preferable, however if the system can be assumed linear then frequency domain simulations are suitable due to their reduced computational demand

Interface damage model for fatigue-driven fracture

A fatigue interface damage model is presented that is based on combining a fatigue evolution law with a static interface damage model. The fatigue model elegantly enables the simulation of crack initiation and propagation in a computationally efficient and accurate way, accounting for mixed-mode loading conditions and fatigue loading of variable amplitude. The main features of the fatigue model are explained and demonstrated by means of an illustrative numerical example

Formulation and numerical implementation of micro-scale boundary conditions for particle aggregates

Novel numerical algorithms are presented for the implementation of micro-scale boundary conditions of particle aggregates modelled with the discrete element method. The algorithms are based on a servo-control methodology, using a feedback principle comparable to that of algorithms commonly applied within control theory of dynamic systems. The boundary conditions are defined in accordance with the large deformation theory, and are imposed on a frame of boundary particles surrounding the interior granular micro-structure. Following the formulation presented in Miehe et al. (Int J Numer Methods Eng 83(8–9): 1206–1236, 2010), first three types of classical boundary conditions are considered, in accordance with (1) a homogeneous deformation and zero particle rotation (D), (2) a periodic particle displacement and rotation (P), and (3) a uniform particle force and free particle rotation (T). The algorithms can be straightforwardly combined with commercially available discrete element codes, thereby enabling the determination of the solution of boundary-value problems at the micro-scale only, or at multiple scales via a micro-to-macro coupling with a finite element model. The performance of the algorithms is tested by means of discrete element method simulations on regular monodisperse packings and irregular polydisperse packings composed of frictional particles, which were subjected to various loading paths. The simulations provide responses with the typical stiff and soft bounds for the (D) and (T) boundary conditions, respectively, and illustrate for the (P) boundary condition a relatively fast convergence of the apparent macroscopic properties under an increasing packing size. Finally, a homogenization framework is derived for the implementation of mixed (D)–(P)–(T) boundary conditions that satisfy the Hill–Mandel micro-heterogeneity condition on energy consistency at the micro- and macro-scales of the granular system. The numerical algorithm for the mixed boundary conditions is developed and tested for the case of an infinite layer subjected to a vertical compressive stress and a horizontal shear deformation, whereby the response computed for a layer of cohesive particles is compared against that for a layer of frictional particles