基于数字图像技术的岩土材料有限元-离散元分析

将数字图像技术与Munjiza提出的有限元法-离散元法（FDEM）耦合分析方法结合,研制了可表征岩石真实非均质性的FDEM分析系统,为从细观角度对岩体进行建模以及研究岩体的破裂机制提供了新途径。该系统借助数字图像技术从岩石断面的图像获取岩石材料的真实细观结构,借助成熟的网格剖分技术,将其映射到FDEM计算网格中,从而克服了原有FDEM在考虑材料非均质性所存在的不足。利用该系统,进行了巴西圆盘劈裂试验的数值模拟,再现了花岗岩在荷载作用下的真实破裂过程。数值模拟结果表明,考虑非均质性的岩样应力分布呈现出非对称性,岩石的细观结构对岩石中裂纹的扩展及应力分布有重要影响。     

矿物粒径对花岗岩单轴压缩特性影响的试验与模拟研究

高放废物深部地质处置目前受到世界各国的高度重视。花岗岩是我国高放废物地质处置工程的候选围岩,深入了解处置库花岗岩的强度及破坏特性对于处置系统的设计及性能评价具有十分重要的意义。作为矿物颗粒的集合体,花岗岩是一种由石英、长石和黑云母等矿物组成的非均质岩石,矿物粒径对其宏观力学特性影响明显。以我国高放地质处置库预选区阿拉善花岗岩为例,选取矿物粒径差异明显的似斑状花岗岩和中粒花岗岩两类岩石,采用单轴压缩试验与数值模拟相结合的方式研究了矿物粒径对岩石力学特性的影响。单轴压缩试验在MTS815岩石力学试验系统进行,数值模拟采用基于离散元的颗粒流程序PFC2D完成。数值模拟过程中,以试件表面图像为基础,采用数字图像处理技术获取岩石内部矿物组分的实际空间分布,从而建立了精确反映花岗岩内部矿物种类及其空间位置的数值模型。利用该模型对花岗岩的单轴压缩试验进行了数值模拟,并与试验结果对比,论证了模型的可靠性。试验及模拟结果表明,阿拉善花岗岩破坏形式为脆性张拉破坏,裂纹大多平行于轴压方向,数字图像数值分析方法可真实地反映材料细观结构。矿物粒径对材料力学特性的影响主要表现为:细粒、等粒结构的岩石强度高,粗粒、不等粒结构的岩石强度低。研究成果可为掌握矿物粒径对岩石强度及变形特性的影响提供依据。     

数字图像技术在岩石细观力学研究中的应用

岩石内部的细观组成和结构决定了其在外力作用下的应力-应变状态,进而控制了其宏观力学响应和破坏机制。数字图像处理技术作为一种材料细观尺度上的空间结构精确测量和数字表述手段,已广泛应用于土、岩石及混凝土的细观结构定量分析中。应用数字图像处理进行的岩石细观力学研究是对岩石力学研究方法的革新。现阶段主要研究内容包括：岩石裂隙隙宽的非接触测量,数字表述岩石结构的非均匀性,进行岩石细观力学行为分析,将提取的岩石数字特征值与相应的岩石物性结合以实现岩石流-固耦合研究,建立基于数字图像技术的岩石细观力学数值模拟方法。在研究相关文献的基础上,对数字图像技术在岩石细观力学定量研究中的成果进行客观评述,探讨各种方法的优缺点,分析展望数字图像技术在岩石细观研究领域中的应用前景。     

基于数字图像处理的非均质岩石材料破坏过程模拟

以三峡花岗岩为例,利用有限差分程序FLAC进行非均质岩石材料的破坏过程模拟.首先,利用花岗岩微观结构观测图像确定数值模拟计算网格的合适尺寸;进一步结合花岗岩的矿物鉴定结果,在数字图像处理的基础上,利用模糊C均值聚类方法,实现了数值计算模型与材料真实细观结构的耦合分析,形成了一套较完整的、更加方便快捷的非均质岩石材料数字...     

基于Otsu图像处理的岩石细观模型重构及数值试验

基于大津法(Otsu)多阈值分割方法,利用数字图像处理技术获取了花岗岩细观结构的表征图像,并结合颗粒流程序,重构了反映岩石非均质结构特征的细观模型。通过单轴压缩试验,分析了花岗岩的失稳破坏过程,并研究了细观结构强度对宏观强度及破坏形态的影响。研究发现:花岗岩试件失稳破坏过程中存在4种裂纹形式,加载初期首先产生界面裂纹,其次是云母裂纹,长石和石英在应力接近极限强度时才破坏;当增大胶结界面或矿物强度时,花岗岩的抗压强度随之增加,但增加梯度逐渐减小;胶结界面或3种矿物细观强度的增大会使该结构处裂纹数量减小,而增加其他结构处的裂纹分布;相对于中等强度的长石和高强度的石英,低强度的胶结界面和云母矿物对花岗岩宏观强度和破坏形态的影响较大,是关键细观结构。将Otsu图像处理技术与颗粒流程序相结合,建立反映岩石真实结构的细观模型,为研究岩石非均质性对宏观力学特性的影响提供了一种有效手段。     

平台中心角对岩石抗拉强度测定影响的数值分析

岩石抗拉强度是岩石的一个重要力学参数,在测试岩石抗拉强度的平台巴西盘劈裂试验中,平台中心角对测试结果会产生重要影响,但测试中对此角度尚无严格规定,从而造成抗拉强度测试结果有很大的离散性。从细观的角度着手,应用数字图像技术对岩石材料内部不同细观介质的空间分布进行精确测量和数值表达,建立与有限元网格之间的相互映射关系,并引入到原有的岩石破坏过程分析系统RFPA中,建立了能反映岩石非均匀性的平台巴西盘数值模型。通过对花岗岩断面图像进行处理,运用数值模型研究了平台中心角对平台巴西盘轴线应力分布、劈裂破坏模式和抗拉强度测定的影响,确定平台中心角的合理值应在20～30?的范围内,这与试验结果相一致,从机制上进一步分析验证了平台中心角的取值范围,为测试中平台中心角的选取提供了理论依据。     

基于CT图像的花岗岩矿物组分与细观结构分析

获取准确的岩石矿物组分与细观结构特征,对从细观层面认识岩石宏观力学性质具有重要意义。CT扫描技术是研究岩石矿物细观结构的有力工具。本文采用高分辨率三维CT扫描系统,获得了花岗岩CT扫描图像;采用阈值分割的方法,综合X射线衍射试验结果确定合理的分割阈值,实现了对花岗岩CT灰度图像的三值化分析;最终获得了花岗岩的主要矿物含量与细观结构特征。研究结果表明:选择合理的分割阈值,可实现花岗岩主要矿物石英、长石、云母的自动识别和定量分析;基于三维重构模型,可获得不同矿物的形状、粒径和空间分布特征。本文方法对定量测试岩石矿物组分和认识其细观结构具有一定的参考意义和价值。     

基于变分水平集方法获取岩石的细观结构

岩土工程材料的细观结构特征会在很大程度上影响其宏观力学行为。随着数字图像处理技术的逐渐成熟,基于变分水平集方法的偏微分方程（PDE）方法在图像处理中越来越受重视。在岩石数字图像的基础上,运用基于变分水平集方法的C-V模型,结合图像的背景填充技术实现花岗岩图片的3相分割（长石、石英和云母）,为进一步力学分析提供条件。上述方法可以量化岩石细观结构,最终获得能够准确反映岩石内部特征的图片,为将来的有限元分析提供真实几何模型。处理结果表明,通过选取适当的迭代次数和数字图像处理技术,该法能够获得较高质量的岩石细观结构量化图。     

岩石细观统计渗流模型研究(Ⅱ)：实例分析

进一步分析了第Ⅰ部分中提出的岩石细观统计渗流模型[1]。讨论了由细观统计渗流模型得到岩石宏观渗透系数的方法。为了利用细观统计渗流模型对岩石损伤破坏过程中的渗透性演化进行分析,引入了各向异性损伤模型,用以对岩石的细观损伤和裂纹的扩展过程进行力学分析,并讨论了细观统计渗流模型与各向异性损伤模型的联合分析方法。对Lac du Bonnet花岗岩损伤过程中的渗透性演化进行了计算分析和验证。结果表明,细观统计渗流模型可以较精确地描述裂纹充分贯通前岩石渗透性的演化。     

岩石细观统计渗流模型研究(II):实例分析

进一步分析了第I部分中提出的岩石细观统计渗流模型[1]。讨论了由细观统计渗流模型得到岩石宏观渗透系数的方法。为了利用细观统计渗流模型对岩石损伤破坏过程中的渗透性演化进行分析,引入了各向异性损伤模型,用以对岩石的细观损伤和裂纹的扩展过程进行力学分析,并讨论了细观统计渗流模型与各向异性损伤模型的联合分析方法。对Lac du Bonnet花岗岩损伤过程中的渗透性演化进行了计算分析和验证。结果表明,细观统计渗流模型可以较精确地描述裂纹充分贯通前岩石渗透性的演化。     

Digital Image Based Approach for Three-Dimensional Mechanical Analysis of Heterogeneous Rocks

Summary This paper presents a digital image based approach for three-dimensional (3-D) numerical simulation and failure analysis of rocks by taking into account the actual 3-D heterogeneity. Digital image techniques are adopted to extract two-dimensional (2-D) material heterogeneity from material surface images. The 2-D image mesostructures are further extrapolated to 3-D cuboid mesostructures by assuming the material surface as a representation of the inner material heterogeneity within a very small depth. The iterative milling and scanning system is set up to generate the 3-D rock mesostructures. A Hong Kong granite specimen is used as an example to demonstrate the procedure of 3-D mesostructure establishment. The mechanical responses and failure process under the conventional Brazilian tensile test condition are examined through numerical analyses. The stress distribution, crack propagation process and failure model of heterogeneous material cases are simulated with a finite difference software. The numerical results indicate that material heterogeneity plays an important role in determining the failure behavior of rocks under external loading.     

Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

花岗岩低温热力效应参数及强度规律研究

长期处于变低温环境中的岩石热力学效应主要表现为材料的力学特性参数、热力系数和饱和冻结状态孔隙冰胀系数等随温度改变而发生变化,这种变化对岩体的变形及强度特性均有显著影响。根据变物性参数的非线性热弹性理论,建立了考虑热力系数和冰胀力系数的本构方程,给出了单轴压缩条件下热力效应系数随温度变化的分析格式和确定方法;借助于花岗岩不同低温干燥和饱和冻结状态的单轴压缩试验资料,获得了花岗岩低温热力效应与温度的关系,并探讨了其对花岗岩抗压强度的影响特性。分析表明,在变低温环境下花岗岩热力系数和冰胀力系数均随相对温差的增大而逐渐减小,热力系数降低速率小于冰胀系数降低速率,在同一温差下热力系数大于冰胀系数;热应力和冰胀应力与温度改变量呈非线性增长关系;花岗岩抗压强度在两种状态下均呈增大趋势,但主要以热应力为主,热力系数引起的试样轴向应力增量总是大于冰胀力引起的轴向应力增量。所给方法及研究成果可用于长期处于变低温状态下材料的力学性能研究及其工程应用。     

Factors Affecting Crack Initiation in Low Porosity Crystalline Rocks

Crack initiation in uniaxial compressive loading of rocks occurs well before the peak strength is reached. The factors that may influence the onset of cracking and possible initiating mechanisms were explored using a discrete element numerical approach. The numerical approach was based on grain-based model that utilized the Voronoi tessellation scheme to represent low porosity crystalline rocks such as granite. The effect of grain size distribution (sorting coefficient ranging from 1.5 to 1.03), grain size (average grain size ranging from 0.75 to 2.25 mm), and the heterogeneities of different mineral grains (quartz, K-feldspar, plagioclase) on the onset of cracking were examined. The modelling revealed that crack initiation appears to be a tensile mechanism in low porosity rocks, and that shear cracking along grain boundaries is only a prominent mechanism near the peak strength. It was also shown that the heterogeneity introduced by the grain size distribution had the most significant effect on peak strength and crack initiation stress. The peak strength ranges from 140 to 208 MPa as the grain size distribution varies from heterogeneous to uniform, respectively. However, the ratio of crack initiation to peak stress showed only minor variation, as the heterogeneity decreases. The other factors investigated had only minor effects on crack initiation and peak strength, and crack initiation ratio.     

A three-dimensional numerical investigation of the fracture of rock specimens containing a pre-existing surface flaw

Three-dimensional surface crack initiation and propagation in two kinds of heterogeneous rocks were numerically investigated via parallel finite element analysis using a supercomputer. Numerically simulated rock specimens containing a pre-existing flaw were subjected to uniaxial compression until failure. The initiation and propagation of wing cracks, anti-wing cracks, and shell-like cracks were reproduced by numerical simulations. The numerically simulated results demonstrate that the further propagation of wing cracks and shell-like cracks stop due to their wrapping (curving) behavior in three-dimensional spaces, even if the applied loads continue to increase. Furthermore, rock heterogeneity could significantly influence crack propagation patterns and the peak uniaxial compressive strengths of rock specimens. Moreover, anti-wing cracks only appeared in relatively heterogeneous rocks, and the peak uniaxial compressive strengths of the specimens were observed to depend on the inclination of the pre-existing flaw. Finally, the mechanism of surface crack propagation is discussed in the context of numerically simulated anti-plane loading tests, wherein it was identified that Mode III loading (anti-plane loading) does not lead to Mode III fracture in rocks due to their high ratio of uniaxial compressive strength to tensile strength. This finding could explain the lateral growth of an existing flaw in its own plane, which is a phenomenon that has not been observed in laboratory experiments.     

高低温冻融循环条件下花岗岩力学性能试验研究

高低温冻融循环条件下岩石的力学特性是岩石力学与工程中的一个重要课题。通过试验研究了花岗岩在不同温度幅值的高低温冻融循环条件下,单轴抗压强度、弹性模量、峰值变形、应力-应变曲线随循环次数变化的规律,并基于Loland损伤力学模型,推导出岩石类材料损伤变量表达式,揭示了花岗岩力学性能劣化与循环次数、循环温度幅值的关系,为高低温冻融循环条件下岩石力学性能的研究提供一定的试验、理论基础     

对一些角隅模型的认识

基于花岗岩的真三轴强度试验结果和双剪强度理论,考虑岩石在三轴压缩状态下和在三轴伸长状态下不同的强度特性,导出了花岗岩在三轴压缩状态下的 函数关系式和在三轴伸长状态下的 函数关系式,并分别将目前常用的 Gudehus-Arygris、William-Warnke、郑颖人、俞茂宏-刘风羽等角隅模型的屈服强度计算结果与花岗岩的真三轴强度试验结果进行了对比,得出了这些角隅模型能较为准确地预测屈服强度等一些有益的认识。     

不同温度条件下饱水风化花岗岩强度及变形特性分析

为研究温度和围压对风化花岗岩抗压强度、 剪切强度参数及变形特性的影响规律, 以新疆天山某一矿区的风化花岗岩为研究对象, 对不同温度（15、 -5、 -15 ℃）、 不同围压（0、 4、 7、 10 MPa）条件下的饱水风化花岗岩进行单轴和三轴压缩试验。结果表明： 相同温度条件下, 围压在0 ~ 10 MPa变化时, 风化花岗岩三轴抗压强度随围压线性增大。在相同围压下, 抗压强度随温度的降低而明显提高, 风化花岗岩的黏聚力c随温度降低而增大, 内摩擦角φ随温度的降低呈增长趋势。弹性模量随围压的增大不断提高, 但随着温度的降低, 增长幅度逐渐减小。泊松比也随温度降低和围压的增加呈增大趋势。温度和围压对风化花岗岩试样的破坏形态影响机制不同。温度降低使矿物颗粒及内部的微裂纹和间隙收缩, 进而胶结强度增大; 荷载和围压的增大会使岩石内部形成微裂隙并逐渐贯通, 同时孔隙冰破碎裂隙充分接触, 进一步增加摩擦力提高岩石强度。     

Numerical modeling of rock failure under dynamic loading with polygonal elements

Polygonal finite elements are gaining an increasing attention in the computational mechanics literature, but their application in rock mechanics is very rare. This paper deals with numerical modeling of rock failure under dynamic loading based on polygonal finite elements. For this end, a damage-viscoplastic constitutive model for rock based on the Mohr-Coulomb criterion with the Rankine criterion as a tensile cutoff is employed and implemented with the polygonal finite element method. Moreover, the mineral mesostructure or rock is described by randomly mapping groups of polygonal elements representing the constituent minerals into a global mesh and assigning these groups with the corresponding mineral material properties. The performance of the polygonal elements is compared with that of the linear and quadratic triangular and bilinear quadrilateral elements in numerical simulations of controlled shear band formation under uniaxial compression and lateral splitting failure in the dogbone tension test. Numerical simulations of uniaxial tension and compression tests as well as dynamic Brazilian disc test under increasing loading rates demonstrate that the present approach predicts the correct failure modes as well as the dynamic increase in strength of rock.     

非均匀岩石介质单轴压缩强度及变形破裂规律的数值模拟

以有限差分法(FDM)为计算框架,利用Weibull分布描述细观单元弹性模量和单轴抗压强度的分布特征,采用弹塑性应变软化本构模型描述细观单元的力学响应,进而建立一种模拟非均匀岩石介质破裂的数值模型。采用该数值模型探讨了单轴压缩时细观均质度m及细观结构对数值试样宏观特性的影响。结果表明,（1）随着细观均质度提高,数值模型的非线性特征逐渐减弱,脆性逐渐增强;宏观峰值强度及弹性模量逐渐增大,峰值强度与lnm呈线性关系,而弹性模量与1/m为线性关系;数值试样表现出由塑性流动破坏至剪切破坏进而为张拉破坏的破坏模式;（2）当细观均质度一定时,细观结构或细观单元空间排列是决定岩石力学行为波动性的主要因素;应力–应变曲线峰前阶段对细观单元的空间排列不敏感,但峰值强度附近及峰后阶段对细观单元的空间排列比较敏感。