A discrete model for an ill-posed nonlinear parabolic PDE

We study a finite-difference discretization of an ill-posed nonlinear parabolic partial differential equation. The PDE is the one-dimensional version of a simplified two-dimensional model for the formation of shear bands via anti-plane shear of a granular medium. For the discretized initial value problem, we derive analytically, and observed numerically, a two-stage evolution leading to a steady-state: (i) an initial growth of grid-scale instabilities, and (ii) coarsening dynamics. Elaborating the second phase, at any fixed time the solution has a piecewise linear profile with a finite number of shear bands. In this coarsening phase, one shear band after another collapses until a steady-state with just one jump discontinuity is achieved. The amplitude of this steady-state shear band is derived analytically, but due to the ill-posedness of the underlying problem, its position exhibits sensitive dependence. Analyzing data from the simulations, we observe that the number of shear bands at time t decays like t^−1/3. From this scaling law, we show that the time-scale of the coarsening phase in the evolution of this model for granular media critically depends on the discreteness of the model. Our analysis also has implications to related ill-posed nonlinear PDEs for the one-dimensional Perona–Malik equation in image processing and to models for clustering instabilities in granular materials.     

Self-organization, localization of shear bands, and aging in loose granular materials

We introduce a mesoscopic model for the formation and evolution of shear bands in loose granular media. Numerical simulations reveal that the system undergoes a nontrivial self-organization process which is governed by the motion of the shear band and the consequent restructuring of the material along it. High density regions are built up, progressively confining the shear bands in localized regions. This results in an inhomogeneous aging of the material with a very slow increase in the mean density, displaying an unusual glassylike system-size dependence.     

Effective temperature and jamming transition in dense, gently sheared granular assemblies

We present extensive computational results for the effective temperature, defined by the fluctuation-dissipation relation between the mean square displacement and the average displacement of grains, under the action of a weak, external perturbation, of a sheared, bi-disperse granular packing of compressible spheres. We study the dependence of this parameter on the shear rate and volume fractions, the type of particle and the observable in the fluctuation-dissipation relation. We find the same temperature for different tracer particles in the system. The temperature becomes independent on the shear rate for slow enough shear suggesting that it is the effective temperature of the jammed packing. However, we also show that the agreement of the effective temperature for different observables is only approximate, for very long times, suggesting that this defintion may not capture the full thermodynamics of the system. On the other hand, we find good agreement between the dynamical effective temperature and a compactivity calculated assuming that all jammed states are equiprobable. Therefore, this definition of temperature may capture an instance of the ergodic hypothesis for granular materials as proposed by theoretical formalisms for jamming. Finally, our simulations indicate that the average shear stress and apparent shear viscosity follow the usual relation with the shear rate for complex fluids. Our results show that the application of shear induces jamming in packings whose particles interact by tangential forces.     

Shear stress fluctuations in the granular liquid and solid phases

We report on experimentally observed shear stress fluctuations in both granular solid and fluid states, showing that they are non-Gaussian at low shear rates, reflecting the predominance of correlated structures (force chains) in the solidlike phase, which also exhibit finite rigidity to shear. Peaks in the rigidity and the stress distribution's skewness indicate that a change to the force-bearing mechanism occurs at the transition to fluid behavior, which, it is shown, can be predicted from the behavior of the stress at lower shear rates. In the fluid state stress is Gaussian distributed, suggesting that the central limit theorem holds. The fiber bundle model with random load sharing effectively reproduces the stress distribution at the yield point and also exhibits the exponential stress distribution anticipated from extant work on stress propagation in granular materials.     

Pressure screening and fluctuations at the bottom of a granular column

We report sets of precise and reproducible measurements on the static pressure at the bottom of a granular column. We make a quantitative analysis of the pressure saturation when the column height is increased. We evidence a great sensitivity of the measurements with the global packing fraction and the eventual presence of shear bands at the boundaries. We also show the limit of the classical Janssen model and discuss these experimental results under the scope of recently proposed theoretical frameworks. Received 7 September 1998 and Received in final form 28 January 1999     

Bifurcations to diversify geometrical patterns of shear bands on granular material

The mechanism to diversify geometrical patterns on granular material was elucidated using a group-theoretic image analysis of patterned shear bands, with associated numerical bifurcation analysis. Pattern formation of granular materials took the course of the evolution of a diamondlike diffuse bifurcation breaking uniformity, followed by further bifurcation, mode jumping, and the formation and disappearance of shear bands through localization. A chaotic explosive increase of possible postbifurcation states was emphasized as a mechanism to diversify geometrical patterns.     

Numerical study of the stress response of two-dimensional dense granular packings

We investigate the Green function of two-dimensional dense random packings of grains in order to discriminate between the different theories of stress transmission in granular materials. Our computer simulations allow for a detailed quantitative investigation of the dynamics which is difficult to obtain experimentally. We show that both hyperbolic and parabolic models of stress transmission fail to predict the correct stress distribution in the studied region of the parameters space. We demonstrate that the compressional and shear components of the stress compare very well with the predictions of isotropic elasticity for a wide range of pressures and porosities and for both frictional and frictionless packings. However, the states used in this study do not include the critical isostatic point for frictional particles, so that our results do not preclude the fact that corrections to elasticity may appear at the critical point of jamming, or for other sample preparation protocols, as discussed in the main text. We show that the agreement holds in the bulk of the packings as well as at the boundaries and we validate the linear dependence of the stress profile width with depth.     

Large-scale flow as an indicator of superplasticity

A model is proposed for superplastic deformation of materials, based on the concept of cooperative grain-boundary slip. The conditions for superplastic deformation are obtained as conditions for coherent shear bands. Analysis of the temperature dependence of the limits of the stress interval for superplastic flow is used as a basis for the introduction of two types of threshold stress that elucidate the cause of the ambiguity in the interpretation of exisiting experimental results. Fiz. Tverd. Tela (St. Petersburg) 39, 2179–2185 (December 1997)     

Strain localisation in granular media

This paper discusses strain localisation in granular media by presenting experimental, full-field analysis of mechanical tests on sand, both at a continuum level, as well as at the grain scale. At the continuum level, the development of structures of localised strain can be studied. Even at this scale, the characteristic size of the phenomena observed is in the order of a few grains. In the second part of this paper, therefore, the development of shear bands within specimen of different sands is studied at the level of the individual grains, measuring grains kinematics with x-ray tomography. The link between grain angularity and grain rotation within shear bands is shown, allowing a grain-scale explanation of the difference in macroscopic residual stresses for materials with different grain shapes. Finally, rarely described precursors of localisation, emerging well before the stress peak are observed and commented.     

Age-dependent transient shear banding in soft glasses

We study numerically the formation of long-lived transient shear bands during shear startup within two models of soft glasses (a simple fluidity model and an adapted "soft glassy rheology" model). The degree and duration of banding depends strongly on the applied shear rate, and on sample age before shearing. In both models the ultimate steady flow state is homogeneous at all shear rates, consistent with the underlying constitutive curve being monotonic. However, particularly in the soft glassy rheology case, the transient bands can be extremely long lived. The banding instability is neither "purely viscous" nor "purely elastic" in origin, but is closely associated with stress overshoot in startup flow.     

Dilatancy and friction in sheared granular media

We introduce a simple model to describe the frictional properties of granular media under shear. We model the friction force in terms of the horizontal velocity and the vertical position z of the slider, interpreting z as a constitutive variable characterizing the contact. Dilatancy is shown to play an essential role in the dynamics, inducing a stick-slip instability at low velocity. We compute the phase diagram, analyze numerically the model for a wide range of parameters and compare our results with experiments on dry and wet granular media, obtaining a good agreement. In particular, we reproduce the hysteretic velocity dependence of the frictional force. Received 16 November 1999     

Effect of shear on an onion texture

Lamellar systems are self-assemblies of surfactant molecules forming planar bilayers separated by layers of solvent. At sufficiently high shear rates, they are known to form spherical objects often referred to as onions. In this paper, we are concerned with the effect of shear on those multi-lamellar vesicles. We measure solvent diffusion by nuclear magnetic resonance (NMR) using a method which is sensitive to the time dependence of mean-squared displacements. This method, combined with NMR velocimetry, allows us to infer onion structure as a function of shear rate, identifying different regimes in which local viscosity is related to the onion size. The role of slip is examined and the stress dependence of wall slip velocities is determined.     

Temporal oscillations of the shear stress and scattered light in a shear-banding-shear-thickening micellar solution

The results of optical and rheological experiments performed on a viscoelastic solution (cetyltrimethylammonium bromide + sodium salicylate in water) are reported. The flow curve has a horizontal plateau extending between two critical shear rates characteristic of heterogeneous flows formed by two layers of fluid with different viscosities. These two bands which also have different optical anisotropy are clearly seen by direct observation in polarized light. At the end of the plateau, apparent shear thickening is observed in a narrow range of shear rates; in phase oscillations of the shear stress and of the first normal stress difference are recorded in a shearing device operating under controlled strain. The direct observation of the annular gap of a Couette cell in a direction perpendicular to a plane containing the vorticity shows that the turbidity of the whole sample also undergoes time dependent variations with the same period as the shear stress. However no banding is observed during the oscillations and the flow remains homogeneous.     

颗粒介质尺度效应的抗剪试验及物理机理分析

针对颗粒介质力学特性的颗粒尺度效应研究,选用土矿物颗粒制备不同颗粒尺度的抗剪试样,进行一系列直剪快剪和三轴抗剪试验,测得了不同颗粒粒径和体分比试样的变形曲线及剪应力强度;基于颗粒间微观作用力与重力比值和胞元体模型,首次从微观和细观角度解释颗粒尺度效应的物理机理.结果表明,随着介质中粗颗粒的比例增加和粒径减小,介质变形特性增强,剪应力强度也随之提高;体分比对变形和强度特性的影响比粒径的影响更加显著.基于介质特性尺度效应物理机理分析,提出衡量介质颗粒聚集和摩擦效应的微重比判别参数以及应变梯度和变形协调微裂纹引起颗粒尺度效应的细观机理解释;文中提出的胞元体模型大大减少了颗粒物质体系的计算自由度,为工业和工程设计的计算建模提供一种可行途径.     

Study of the shear-rate dependence of granular friction based on community detection

In this study, computer simulations are performed on three-dimensional granular systems under shear conditions. The system comprises granular particles that are confined between two rigid plates. The top plate is subjected to a normal force and driven by a shearing velocity. A positive shear-rate dependence of granular friction, known as velocity-strengthening, exists between the granular and shearing plate. To understand the origin of the dependence of frictional sliding, we treat the granular system as a complex network, where granular particles are nodes and normal contact forces are weighted edges used to obtain insight into the interiors of granular matter. Community structures within granular property networks are detected under different shearing velocities in the steady state. Community parameters, such as the size of the largest cluster and average size of clusters, show significant monotonous trends in shearing velocity associated with the shear-rate dependence of granular friction. Then, we apply an instantaneous change in shearing velocity. A dramatic increase in friction is observed with a change in shearing velocity in the non-steady state. The community structures in the non-steady state are different from those in the steady state. Results indicate that the largest cluster is a key factor affecting the friction between the granular and shearing plate.     

Interrupted shear of granular media

The response of a granular material during a stop-and-go shear experiment is investigated using an annular shear cell and silicagel powders of different particle sizes. The experimental results are examined on the basis of the Dieterich-Rice-Ruina model for solid friction. In addition to making this analogy with solid friction, we describe a new instability that is observed when restarting shear, where the powder bed is found to slip and compact for short hold times but only dilates for long hold times. The minimum hold time to restore a non-slip behaviour has been investigated for different size particles and normal loadings. The observed dependencies show analogies between this behaviour and the sliding rearrangements seen above the stick-slip threshold.     

A numerical examination of shear banding and simple shear non-coaxial flow rules

Strain localisation and shear band formation is frequently observed during the handling and flow of dense phase particulate materials. However, a complete understanding of how shear bands form and what happens inside shear bands is still lacking. In order to address this problem, discrete particle simulations have been carried out to examine the detailed processes that occur at the grain scale associated with the initiation and development of shear bands. To reliably identify the continuum model applicable within a shear band is difficult due to the small number of particles/contacts involved. However, it is normally accepted that the mode of deformation within a shear band is one of simple shear. Consequently, simple shear simulations have been performed in order to determine the evolution of the stress tensor, dilation rate, and the principal directions of stress and strain-rate. It is demonstrated that the corresponding non-coaxial flow rule is equivalent to that suggested by Tatsuoka et al. (Géotechnique 38 148 (1988)). Furthermore, at fully developed flow when there is no further change in volume, the stress and strain-rate directions are coaxial and the flow rule is that proposed by Hill (The Mathematical Theory of Plasticity (Oxford University Press, 1950) p. 294).     

A non-viscous-featured fractograph in metallic glasses

A fractograph of non-viscous feature but pure shear-offsets was found in three-point bending samples of a ductile Pd–Cu–Si metallic glass. A sustainable shear band multiplication with large plasticity during notch propagation was observed. Such non-viscous-featured fractograph was formed by a crack propagation manner of continual multiple shear bands formation in front of the crack-tip, instead of the conventional rapid fracture along shear bands. With a 2D model of crack propagation by multiple shear bands, we showed that such fracture process was achieved by a faster stress relaxation than shear-softening effect in the sample. This study confirmed that the viscous fracture along shear bands could be not a necessary process in ductile metallic glasses fracture, and could provide new ways to understand the plasticity in the shear-softened metallic glasses.     

Probing shear-band initiation in metallic glasses

In?situ acoustic emission monitoring is shown to capture the initiation of shear bands in metallic glasses. A model picture is inferred from stick-slip flow in granular media such that the origin of acoustic emission is attributed to a mechanism of structural dilatation. By employing a quantitative approach, the critical volume change associated with shear-band initiation in a metallic glass is estimated to be a few percent only. This result agrees with typical values of excess free volume found in the supercooled liquid regime near the glass transition temperature.