Exploring the influence of interparticle friction on critical state behaviour using DEM

Understanding the extent to which discrete element method (DEM) simulations can capture the critical state characteristics of granular materials is important to legitimize the use of DEM in geomechanics. This paper documents a DEM study that considered the sensitivity of the critical state response characteristics to the coefficient of interparticle friction (μ) using samples with gradings that are representative of a real soil. Most of the features that are typically associated with sand behaviour at the critical state were seen to emerge from the DEM simulation data. An important deviation occurs when high μ values (μ ≥ 0.5) are used, as has been the case in a number of prior DEM studies. While there is a systematic variation in the critical state behaviour with μ for μ < 0.5, when μ ≥ 0.5, the behaviour at the critical state seems to be insensitive to further increases in μ. In contrast to observations of conventional soil response, when μ ≥ 0.5, the void ratio at the critical state initially increases with increasing mean effective stress (p′). Analysis of the DEM data and use of simple models of isolated force chains enabled some key observations. When ‘floating’ particles that do not transmit stress are eliminated from the void ratio calculation, the void ratio at the critical state decreases consistently with increasing p′. There is a transition from sliding to rolling behaviour at the contact points as μ increases. Beyond a limiting value of μ, further increases in μ do not increase the buckling resistance of individual strong force chains. Copyright © 2014 John Wiley & Sons, Ltd.     

Intersecting dilated convex polyhedra method for modeling complex particles in discrete element method

This paper describes a new method for representing concave polyhedral particles in a discrete element method as unions of convex dilated polyhedra. This method offers an efficient way to simulate systems with a large number of (generally concave) polyhedral particles. The method also allows spheres, capsules, and dilated triangles to be combined with polyhedra using the same approach. The computational efficiency of the method is tested in two different simulation setups using different efficiency metrics for seven particle types: spheres, clusters of three spheres, clusters of four spheres, tetrahedra, cubes, unions of two octahedra (concave), and a model of a computer tomography scan of a lunar simulant GRC‐3 particle. It is shown that the computational efficiency of the simulations degrades much slower than the increase in complexity of the particles in the system. The efficiency of the method is based on the time coherence of the system, and an efficient and robust distance computation method between polyhedra as particles never intersect for dilated particles. © 2014 The Authors. International Journal for Numerical and Analytical Methods in Geomechanics published by John Wiley & Sons Ltd.     

基于多层次模糊数学的中国低煤阶煤层气选区评价标准——以吐哈盆地为例

中国低煤阶含煤盆地煤层厚度大,储层物性好,煤层气资源丰富,而对于低煤阶煤层气有利区优选评价目前还没有成型的模型,本文以占全国低煤阶煤层气总资源量50%的吐哈盆地为例,运用多层次模糊数学的思路,对资源因素、煤储层因素、保存因素以及次一级影响因素共计19项参数,进行关键要素定量排序,针对吐哈盆地特有的成藏地质条件建立了煤层气选区评价标准,通过"赋值加权求和定量排序"的方法,探讨了中国低煤阶煤层气评价标准及主控因素对煤层气赋存的影响。研究表明:吐哈盆地二级构造单元沙尔湖浅凹陷、大南湖浅凹陷和哈密凹陷是适合煤层气勘探开发的有利目标区。     

Incrementally non‐linear plasticity applied to rock joint modelling

Rock joint constitutive modelling is discussed through two new rock joint constitutive relations and a discrete numerical model. Regarding the constitutive relations, we emphasise the number of ‘tensorial zones’, that is, domains of constitutive incremental linearity; they involve four zones for the first (called ‘quadrilinear’) and an infinite number for the second one (called ‘incrementally nonlinear’). Using these formulations, a large class of loading paths can be considered. Hardening through shearing and relations between the normal and tangential directions of the joint (e.g., dilatancy) can be described. Their predictive abilities are checked. Plastic features are included even if the relations are defined outside the elasto‐plastic formalism. These relations obey, hence, the physical evidence as the plastic limit criterion and flow rule. The flow rule is nonassociated, and the corresponding link with the nonsymmetry of the constitutive matrix is examined. Comparisons between the two relations and the discrete numerical model, that is, a direct numerical simulation, which is fundamentally different, also are discussed within the context of infilled rock joints. Copyright © 2011 John Wiley & Sons, Ltd.     

A general family of counting distributions suitable for modeling cluster phenomena

The two-dimensional spatial distribution of precious stones, such as diamonds in alluvial and coastal deposits, shows a high degree of clustering. Usually, stones tend to gather in relatively small clusters or traps, made by potholes, gullies, or small depressions in the rough bedcock. Therefore, when taking samples of such deposits, discrete distributions of the number of stones counted in each sample yield an extreme skewness. Most samples have no stones, whereas samples containing a few hundred stones are not unusual. This paper constructs a model and a method for fitting a new and general family of counting distributions based on the Neyman-Scott cluster model and the mixed Poisson process, which can be used to model a differing degree of clustering. General recursion equations for the discrete probabilities of these distributions are derived. Application of this model to simulated data shows that information such as cluster size, number of point events per cluster, and number of clusters per measurement unit can be extracted easily from this model. Fitting the model to data of two real diamond deposits of a totally different nature—small rich clusters of Namibia versus larger but less rich clusters of Guinea—demonstrates its flexibility.     

论城镇地籍测量中面积量算问题及其处理方法

实践表明 ,现行有关面积量算成果检核办法对于解析坐标法面积量算并不完全适应。为此文章对现行面积量算限差的制定及其不完全适合解析法面积量算检核原因等方面进行分析 ,并针对解析坐标法面积量算检核方法提出了一些处理方法     

Propagation of drought through groundwater—a new approach using linear reservoir theory

The effect of drought on groundwater heads and discharge is often complex and poorly understood. Therefore the propagation of a drought from groundwater recharge to discharge and the influence of aquifer characteristics on the propagation was analysed by tracking a drought in recharge through a linear reservoir. The recharge was defined as a sinusoid function with a period of 1 year. The decrease in recharge owing to drought was simulated by multiplying the recharge during 1 year with a drought fraction between 0 and 1, which represents a decrease in the recharge of 100 to 0%, respectively. The droughts were identified using the threshold level approach, with a threshold that is constant in time. For this case analytical formulations were derived, which express the drought duration and deficit in the groundwater discharge in terms of the decrease in recharge, the reservoir coefficient that characterizes aquifer properties and the height of the threshold level. The results showed that the delay in the groundwater system caused a shift of the main part of the decrease in recharge from the high‐flow to the low‐flow period. This resulted in an increase in drought deficit for discharge compared with the drought deficit for recharge. Also the development of multiyear droughts caused an increase in drought deficit. The attenuation in the groundwater system caused a decrease in drought deficit. In most cases the net effect of these processes was an increase of drought deficit as a result of the propagation through groundwater. Only for small droughts the deficit decreased from recharge to discharge. The amount of increase or decrease depends on the reservoir coefficient and the severity of the drought. Under most conditions a maximum in the drought deficit occurred for a reservoir coefficient of around 200 days. Copyright © 2003 John Wiley & Sons, Ltd.     

An Analytical Solution of the Lagrange Equations Valid also for Very Low Eccentricities: Influence of a Central Potential

This paper presents an analytic solution of the equations of motion of an artificial satellite, obtained using non singular elements for eccentricity. The satellite is under the influence of the gravity field of a central body, expanded in spherical harmonics up to an arbitrary degree and order. We discuss in details the solution we give for the components of the eccentricity vector. For each element, we have divided the Lagrange equations into two parts: the first part is integrated exactly, and the second part is integrated with a perturbation method. The complete solution is the sum of the so-called “main” solution and of the so-called “complementary” solution. To test the accuracy of our method, we compare it to numerical integration and to the method developed in Kaula (Theory of Satellite Geodesy, Blaisdell publ. Co., New York. 1966), expressed in classical orbital elements. For eccentricities which are not very small, the two analytical methods are almost equivalent. For low eccentricities, our method is much more accurate.     

Review of validation of the discontinuous deformation analysis (DDA) method

Over the last decade, researchers in the discontinuous deformation analysis (DDA) community have dedicated a great deal of effort to document the accuracy of the method by performing validation studies. This paper contains a summary of more than 100 published and unpublished validation studies which comprise the body of DDA validation information to which the authors have access. The studies are grouped into three general categories: (a) validation with respect to analytical solutions, (b) validation with respect to results of other numerical techniques, and (c) validation with respect to laboratory and field data. Three general techniques for validation are described: qualitative assessment visually examining runtime behaviour of simulations, semi‐quantitative assessment comparing numerical results of simulations, and quantitative where numerical simulation results are evaluated in detail with respect to similar analytical, laboratory or field results. We find that for many of the problems addressed by the papers in this review, DDA performs more than adequately for engineering analysis. Copyright © 2006 John Wiley & Sons, Ltd.