A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.     

基于细观试验和离散元法的盐岩力学特性

当前盐岩的宏观力学模型通常是唯象模型,不能很好地解释盐岩受力变形破坏的真正物理基础。盐岩是由于化学沉积而形成的矿物集合体,是一种主要由NaCl和少量杂质组成的多晶体,其变形机制主要由晶粒与晶界的力学特性控制。通过扫描电镜（SEM）,获得盐岩晶粒的微细观结构特征,采用分子动力学方法和纳米压痕技术,确定盐岩晶粒和晶界的微细观力学参数;将盐岩晶粒作为块体,基于Voronoi多边形技术,建立盐岩的微细观数值模型;利用离散元方法,对盐岩试件在单轴压缩和直剪条件下的宏观力学行为进行了数值模拟。数值模拟结果与宏观力学试验结果吻合度高,表明基于盐岩微细观晶粒结构特征并结合离散元数值模拟的方法能够较好地研究盐岩的宏观力学性能及其材料物理基础。     

On the strength of transversely isotropic rocks

Accurate prediction of strength in rocks with distinct bedding planes requires knowledge of the bedding plane orientation relative to the load direction. Thermal softening adds complexity to the problem since it is known to have significant influence on the strength and strain localization properties of rocks. In this paper, we use a recently proposed thermoplastic constitutive model appropriate for rocks exhibiting transverse isotropy in both the elastic and plastic responses to predict their strength and strain localization properties. Recognizing that laboratory‐derived strengths can be influenced by material and geometric inhomogeneities of the rock samples, we consider both stress‐point and boundary‐value problem simulations of rock strength behavior. Both plane strain and 3D loading conditions are considered. Results of the simulations of the strength of a natural Tournemire shale and a synthetic transversely isotropic rock suggest that the mechanical model can reproduce the general U‐shaped variation of rock strength with bedding plane orientation quite well. We show that this variation could depend on many factors, including the stress loading condition (plane strain versus 3D), degree of anisotropy, temperature, shear‐induced dilation versus shear‐induced compaction, specimen imperfections, and boundary restraints.     

3D numerical modelling of bit–rock fracture mechanisms in percussive drilling with a multiple‐button bit

This paper considers numerical modelling of rock fracture induced by dynamic bit–rock interaction in percussive drilling. The work presented here extends the author's earlier research on the topic from the axisymmetric case to 3D case. The numerical method for modelling rock fracture includes a constitutive model for rock and a contact mechanics‐based technique to simulate the bit–rock interaction. The constitutive model is based on a combination of the recent viscoplastic consistency model, the isotropic damage concept and a parabolic compression cap. This model is improved here from its earlier state by calibrating the softening laws using fracture energies G_Ic and G_IIc in tension and compression, respectively. Moreover, the viscosity modulus in tension is calibrated based on the dynamic Brazilian disc test. With these enhancements, the developed method is applied to 3D case of the bit–rock interaction problem assuming one symmetry plane. Single impact with single and multiple‐button bits is simulated. In the latter case, an initial borehole is modelled in order to simulate the usual in‐situ drilling conditions. The different failure types observed in the experiments as well as the interaction between the buttons resulting in chipping are realistically captured in the simulations. Copyright © 2012 John Wiley & Sons, Ltd.     

适用于岩石力学细观实验研究的加载仪

根据岩石工程力学性质研究的需要, 研制了岩石细观力学加载仪( YXJY-5T)。 该仪器设备配装在光学体视显微镜下 , 可以观察岩样在加载过程中, 四个平面( 岩石试样为长方形柱体)变形破坏的全过程, 并获得岩样的应力-应变曲线和相对应的细观结构变化的图像。     

含石量和颗粒破碎对土石混合料强度的影响研究

The characteristics of particle breakage and shear strength of soil-rock aggregate with six rock contents under six normal pressures were studied from macro and mecro perspectives by large-scale direct shear test, particle observation test and particle sieving test. The relationship between macroscopic shear strength properties and mecroscopic particle breakage characteristics was established, thus further revealing the influence mechanism of rock content and particle breakage on the shear strength characteristics of soil-rock aggregate. The results showed that particle breakage mainly occurred near the shear plane. The breakage morphology can be divided into surface grinding, local fracture, complete fracture and complete breakage, resulting from the stress concentration caused by uneven contact forces between particles. Due to particle breakage, the content of fine particles increased, coarse grains decreased, and intermediate grains fluctuated. The relative particle breakage Br increased with the increase of normal pressure ?n or rock content P5, which accorded with the function of two variables. With the increase of normal pressure ?n, the shear strength τ increased nonlinearly and met the modified M-C strength criterion. When the rock content P5 increased, the cohesive force c0 of soil-rock aggregate decreased, the internal friction angle ?0 of soil-rock aggregate increased, and the non-linear parameter Δ? increased. Particle breakage was the direct cause of non-linear strength characteristics of soil-rock aggregate.     

Meso‐scale particle modeling of concrete deterioration caused by alkali‐aggregate reaction

A meso‐scale particle model is presented to simulate the expansion of concrete subjected to alkali‐aggregate reaction (AAR) and to analyze the AAR‐induced degradation of the mechanical properties. It is the first attempt to evaluate the deterioration mechanism due to AAR using the discrete‐element method. A three‐phase meso‐scale model for concrete composed of aggregates, mortar and the interface is established with the combination of a pre‐processing approach and the particle flow code, PFC2D. A homogeneous aggregate expansion approach is applied to model the AAR expansion. Uniaxial compression tests are conducted for the AAR‐affected concrete to examine the effects on the mechanical properties. Two specimens with different aggregate sizes are analyzed to consider the effects of aggregate size on AAR. The results show that the meso‐scale particle model is valid to predict the expansion and the internal micro‐cracking patterns caused by AAR. The two different specimens exhibit similar behavior. The Young's modulus and compressive strength are significantly reduced with the increase of AAR expansion. The shape of the stress–strain curves obtained from the compression tests clearly reflects the influence of internal micro‐cracks: an increased nonlinearity before the peak loading and a more gradual softening for more severely affected specimens. Similar macroscopic failure patterns of the specimens under compression are observed in terms of diagonal macroscopic cracks splitting the specimen into several triangular pieces, whereas localized micro‐cracks forming in slightly affected specimens are different from branching and diffusing cracks in severely affected ones, demonstrating different failure mechanisms. Copyright © 2012 John Wiley & Sons, Ltd.     

岩石单颗粒压缩破碎的数值模拟分析

岩石颗粒破碎是影响粒状材料剪切强度和变形的最主要因素, 岩石颗粒破碎并不是想象的那么难, 像花岗岩颗粒有时在很小的压应力作用下就可以破碎。岩石单颗粒破碎的物理试验结果常常很离散, 完成大量单颗粒破碎的物理试验费时费力不现实, 采用离散单元法(Discrete element method, DEM)PFC软件模拟单颗粒压缩破碎试验, 既能克服单颗粒破碎物理试验的缺陷, 又能解决单颗粒破碎物理试验工作量大的难题, 是研究单颗粒破碎的理想选择。基于DEM的软件PFC2D, 将粒径为0.075～0.1245mm的基本粒子捆绑成不同粒径的单颗粒, 模拟岩石单颗粒压缩破碎试验, 观察颗粒破碎演化过程, 统计单颗粒破碎强度。计算单颗粒压缩破碎后颗粒分布的分维, 验证单颗粒破碎强度的分形模型和单颗粒破碎强度的尺寸效应。文中引用玄武岩单颗粒破碎试验结果, 与单颗粒破碎的离散单元模拟结果进行比较, 验证单颗粒破碎强度的尺寸效应和修正的Weibull理论的离散单元模拟结果。     

晶体及矿物颗粒大小对岩土材料力学性质的影响

岩石常见较大的晶体或者矿物颗粒,混凝土中是骨料,通过团簇模拟大颗粒的力学行为、团簇可以破裂。根据设计的相同数量、相同位置、不同半径的大颗粒数值单轴压缩试验,在虚拟试验条件下,考察了颗粒大小对材料力学响应的影响。通过分析颗粒材料的破裂形态、裂纹扩展过程、应力-应变曲线和破裂能量演化规律发现：大颗粒具有明显的增强特性,有阻止裂纹扩展的作用,破裂多绕大颗粒发展;增强幅度随颗粒半径的增加呈单增趋势,半径较小时,增强效果不明显。     

干湿循环条件下红砂岩软弱夹层微结构与剪切强度的关联性

极旱极涝频繁交替的极端气候事件会加剧岩土风化,降低强度,影响岩土工程结构稳定。采用直剪试验和扫描电镜技术,研究了对顺层红砂岩边坡稳定性起控制作用的软弱夹层在常温状态、干湿循环条件下（60℃风干）的力学特性及微观结构变化规律,通过相关分析、逐步回归方法得到了影响强度的微结构显著变量,构建了反映二者关联性的回归方程。研究表明：红砂岩软弱夹层剪切强度及黏聚力在1～5次循环过程中急剧降低,5次后变化趋缓,循环次数对内摩擦角变化影响微小;干湿循环次数增加,平均面积、平均直径、平均周长、定向概率熵等4个微结构参数减小,并与黏聚力呈现正相关关系;颗粒个数、平均形状系数、形态分布分形维数3个微观参数与黏聚力呈现负相关关系;分析认为,干湿循环过程中宏观力学特性蜕化及微观结构参数变化之间显著相关,干湿循环造成的岩土颗粒胀缩、碎裂、结构破坏是导致红砂岩软弱夹层剪切强度蜕化的主要原因。     

土石混合体压缩性的三维颗粒力学研究

土石混合体是由作为骨料的砾石或卵石与作为填充料的黏土或砂土组成的特殊工程地质材料,无论天然形成或人工合成,其力学特性与均质土体或岩体相比均有较大差异。以昆明新机场高填方中大量采用土石混合体填料为工程背景,在三维颗粒流程序中对特定的土石混合体试样的侧限压缩试验进行数值模拟,比较了不同级配条件下土石混合体模型的微观结构及基本力学物理性能。通过比较侧限压缩模量和压缩后孔隙率,研究了含石率对土石混合体骨架效应等结构性的影响,并提出了工程中建议采用的土石比例区间,在颗粒流方法用于实际工程方面进行了有益的探索。     

Relationships between texture and mechanical properties of hybrid rocks from the Jaala–Iitti complex, southeastern Finland

The relationships between petrographical and mechanical properties of rock aggregate raw materials from the hybridised, subvolcanic Jaala–Iitti complex, southeastern Finland, were investigated. Petrography was quantified from polished thin sections with a polarising microscope to determine the modal composition and grain size distribution, and resistance to fragmentation and abrasion were determined. Abundance of fine-grained minerals (especially of hornblende), fine grain size-dominated grain size distribution, uniform spatial dispersion of hornblende crystals, and intense micrographic intergrowth texture with interlocking grain boundaries were found to have the greatest positive influence on the mechanical properties. The results showed the potentiality of hybridised rocks as raw materials for high quality aggregates that can resist fragmentation and abrasion.     

Discrete modelling of time‐dependent rockfill behaviour

This modelling study deals with the time‐dependent behaviour of rockfill media, which is of particular interest during the life of rockfill dams. Breakage of rock blocks and crack propagation are the main processes responsible for rockfill creep and collapse. The modelling procedure presented here is performed on two scales: on the rock block scale, where the grain is taken to be an assembly of rigid particles initially endowed with cohesive bonds, and on the rockfill scale, which is taken to involve a set of breakable grains interacting via contact and friction processes. The grain breakage process is described in term of a thermodynamically consistent damage interface model, where the damage is a gradual delayed process. This model was implemented in a non‐smooth contact dynamics code. The effects of the main parameters involved were analysed by performing numerical studies. The ability of the model to predict the creep behaviour of rockfill media is confirmed by presenting several simulations. Copyright © 2008 John Wiley & Sons, Ltd.     

Microplane constitutive model for porous isotropic rocks

The paper deals with constitutive modelling of contiguous rock located between rock joints. A fully explicit kinematically constrained microplane‐type constitutive model for hardening and softening non‐linear triaxial behaviour of isotropic porous rock is developed. The microplane framework, in which the constitutive relation is expressed in terms of stress and strain vectors rather than tensors, makes it possible to model various microstructural physical mechanisms associated with oriented internal surfaces, such as cracking, slip, friction and splitting of a particular orientation. Formulation of the constitutive relation is facilitated by the fact that it is decoupled from the tensorial invariance restrictions, which are satisfied automatically. In its basic features, the present model is similar to the recently developed microplane model M4 for concrete, but there are significant improvements and modifications. They include a realistic simulation of (1) the effects of pore collapse on the volume changes during triaxial loading and on the reduction of frictional strength, (2) recovery of frictional strength during shearing, and (3) the shear‐enhanced compaction in triaxial tests, manifested by a deviation from the hydrostatic stress–strain curve. The model is calibrated by optimal fitting of extensive triaxial test data for Salem limestone, and good fits are demonstrated. Although these data do not cover the entire range of behaviour, credence in broad capabilities of the model is lend by its similarity to model M4 for concrete—an artificial rock. The model is intended for large explicit finite‐element programs. Copyright © 2002 John Wiley & Sons, Ltd.     

节理岩体卸荷强度特性的试验研究

卸荷状态工程岩体的强度特性与传统的加载岩体有本质的区别。目前对于岩石卸荷强度特性的研究比较多见,但对于含节理的岩体在卸荷应力状态下的强度特性,尤其是节理面对其强度的影响研究实属少见。本文通过岩体三轴模拟试验研究了卸荷应力状态下节理岩体的破坏特征,突出考虑了节理面性质对卸荷强度的影响,在此基础上提出了节理岩体卸荷强度准则的一般表达式并进行了对比验证。     

Multiscale modeling of a sensitive marine clay

This paper examines the mechanical behavior of a sensitive marine clay. Various laboratory tests on intact and reconstituted samples of Guinea Gulf marine clay were performed under isotropic compression and drained triaxial compression at constant confining stresses. Microstructure analysis on intact and reconstituted samples was also carried out under different loading conditions. The effect of inter‐aggregates bonding on mechanical properties is discussed. Based on experimental analysis, a new modeling method is proposed. In this approach, the clay is regarded as an assembly of aggregates of clay particles. An inter‐aggregate contact law is introduced relating contact forces to aggregates relative displacements. The deformation of the assembly can be obtained by integrating the movement of the inter‐aggregate contacts in all orientations. Thus, the effect of inter‐aggregates bonds and debonding is considered in a direct way. The model is evaluated through comparisons between the predicted and measured results on Guinea Gulf marine clay. The evolutions of local stresses, strains, and bonds in inter‐aggregates planes are discussed to explain the anisotropy induced by the applied loading. Copyright © 2010 John Wiley & Sons, Ltd.     

考虑颗粒破碎引起级配演变的粗粒料屈服函数研究

级配作为粗粒料的重要物理特性,显著影响着粗粒料本身的力学性质。准确预测土体在加载过程中的级配演变,是有效分析土体结构全寿命周期强度和变形特性的基础。引入Einav的分形破碎理论,认为土体在加载过程中的颗粒破碎耗能增量正比于颗粒分形破碎率增量,结合考虑颗粒破碎的能量平衡方程,选取已有文献中的试验数据,对粗粒料加载过程中的级配演变进行了研究,提出了一种预测试样加载过程中级配演变的方法。在此基础上,建立了一个可以反映粗粒料级配演变的屈服函数,对该屈服面函数进行了初步探讨。每一个屈服面对应土体的一个级配演化曲线;屈服面是土体剪切应变的等值面,屈服面的轨迹和剪应变大小密切相关。     

Micromechanical inspection of incremental behaviour of crushable soils

In granular soils grain crushing reduces dilatancy and stress obliquity enhances crushability. These are well-supported specimen-scale experimental observations. In principle, those observations should reflect some peculiar micromechanism associated with crushing, but which is it? To answer that question the nature of crushing-induced particle-scale interactions is here investigated using an efficient DEM model of crushable soil. Microstructural measures such as the mechanical coordination number and fabric are examined while performing systematic stress probing on the triaxial plane. Numerical techniques such as parallel and the newly introduced sequential probing enable clear separation of the micromechanical mechanisms associated with crushing. Particle crushing is shown to reduce fabric anisotropy during incremental loading and to slow fabric change during continuous shearing. On the other hand, increased fabric anisotropy does take more particles closer to breakage. Shear-enhanced breakage appears then to be a natural consequence of shear-enhanced fabric anisotropy. The particle crushing model employed here makes crushing dependent only on particle and contact properties, without any pre-established influence of particle connectivity. That influence does not emerge, and it is shown how particle connectivity, per se, is not a good indicator of crushing likelihood.

     

An enhanced constitutive model for crushable granular materials

Studies in the past have tried to reproduce the mechanical behaviour of granular materials by proposing constitutive relations based on a common assumption that model parameters and parameters describing the properties, including gradation of individual grains are inevitably linked. However successful these models have proved to be, they cannot account for the changes in granular assembly behaviour if the grains start to break during mechanical loading. This paper proposes to analyse the relation between grading change and the mechanical behaviour of granular assembly. A way to model the influence of grain breakage is to use a critical state‐based model. The influence of the amount of grain breakage during loading, depending on the individual grain strength and size distribution, can be introduced into constitutive relations by means of a new parameter that controls the evolution of critical state with changes in grain size distribution. Experimental data from a calcareous sand, a quartz sand, and a rockfill material were compared with numerical results and good‐quality simulations were obtained. The main consequences of grain breakage are increased compressibility and a gradual dilatancy disappearance in the granular material. The critical state concept is also enriched by considering its overall relation to the evolution of the granular material. Copyright © 2009 John Wiley & Sons, Ltd.     

Correlation of mineralogical and textural characteristics with engineering properties of selected granitic rocks from Turkey

Granitic rocks show a variety of engineering properties that may affect quarrying operations, tunneling, mining, slope stability and the use of rock as a construction material. The physical and mechanical properties are a function of the mineralogical and textural characteristics of the rock. The purpose of this study is to apply correlation analysis to investigate the relationships between petrographical and engineering properties of granitic rocks. A variety of granitic rock samples from different parts of Turkey were subjected to petrographic studies. The same samples were then tested to determine specific gravity, dry and saturated unit weight, water absorption, effective and total porosity, sonic velocity, Schmidt hardness, point load strength index, uniaxial compressive strength, tensile strength and modulus of elasticity. The relationships between these properties and the petrographical characteristics are described by simple regression analyses. The study revealed that the influence of the textural characteristics on the engineering properties appears to be more important than the mineralogy. It also determined that the types of contacts, grain (mineral) shape and size significantly influence the engineering properties of the granitic rocks.