基于单位分解的扩展径向点插值无网格法

针对线弹性断裂力学问题,提出扩展径向点插值无网格法（X-RPIM）。该方法基于单位分解思想,在传统径向点插值无网格法的位移模式中加入扩展项来描述裂纹两侧的不连续位移场和裂尖奇异场。由于其形函数具有Kronecker ? 函数性质,易于施加本质边界条件。详细描述了X-RPIM不连续位移模式的建立,支配方程的离散形式以及J积分计算混合模式裂纹的应力强度因子的实现过程,讨论了不同积分区域对应力强度因子的影响。数值算例分析证明了该方法在求解断裂问题时的可行性和有效性,同时说明扩展径向点插值无网格法在模拟裂纹扩展问题时具有良好的前景。     

爆炸气体驱动下岩体破裂的有限元-离散元模拟

在原有有限元/离散元（FEM/DEM）耦合分析方法中,实现了一种新的爆破计算模型。该模型考虑了在爆生气体的作用下,随着裂隙的扩展,气体占据的体积不断增大,气体压力逐渐减小这一问题。同时考虑了气体嵌入与爆腔联通的裂隙对裂隙的作用力。克服了原有FEM/DEM方法中的爆破模型仅仅将压力施加于爆腔四周的岩壁上,无法考虑爆生气体嵌入生成的裂隙对裂隙的作用。提出了一种新颖的贯通裂隙网络形成的递归搜索算法,只需通过编写一个简单的递归函数,即可实现复杂裂隙网络的搜索,采用一种简洁的方法完成了对复杂问题的处理。最后通过一个爆破算例,结果表明FEM/DEM方法可以对爆炸过程中应力波的传播及岩体中裂纹的萌生、扩展进行全程捕捉,展现了该方法用于爆破模拟的潜力。     

Cracking risk of partially saturated porous media—Part II: Application to drying shrinkage

Drying of deformable porous media results in their shrinkage, and it may cause cracking provided that shrinkage deformations are hindered by kinematic constraints. Herein, we focus on slow drying of an initially water‐saturated sample of a microheterogeneous poroelastic material damaged by parallel mesocracks. The cracking risk is analyzed by means of the thermodynamics‐based microporoelasticity model described in the companion paper (Part I), which is extended toward consideration of the hierarchical organization of cracked argillite. Drying of a material sample is studied in a framework where macrodisplacements in direction of the crack normal are blocked, while elsewise macrostress‐free boundary conditions prevail. The model implies that the opening/closure behavior of the cracks is governed by an effective pressure, in which the average crack (under)pressure, making the crack opening smaller, competes with the average micropore (under)pressure that makes the crack opening larger. The driving force for crack propagation is a quadratic function of this effective pressure. The model proposes that if drying shrinkage deformations are hindered by kinematic constraints, onset of cracking becomes possible once air entry into the cracks is observed. Copyright © 2009 John Wiley & Sons, Ltd.     

A practical method to estimate crack openings in concrete structures

Numerical modelling of concrete cracking requires robust models able to describe opening and propagation of cracks. Structural concrete codes provide practical relations to describe crack openings. However, these empirical methods were developed for specific structures and cannot be used for general applications. Here, a continuous modelling approach based on damage mechanics is used to compute crack openings in a tie‐beam concrete structure. We propose a post‐processing method to extract crack openings from a continuum damage finite element computation. This method can be applied to all continuum damage/plasticity models. The tie‐beam concrete is characterized by a weak stress gradient; this aspect complicates predictions of crack positions and number. A stochastic method is used to take into account the spatial variability in concrete properties and create a spatially correlated random property field. Copyright © 2010 John Wiley & Sons, Ltd.     

爆破开挖对岩体含水裂纹扩展的扰动机制

岩体裂纹的水力劈裂是岩体开裂渗漏甚至施工涌水的重要影响因素之一,也是岩土工程界的研究热点。从断裂力学角度分析了爆破开挖对岩体含水裂纹扩展的扰动作用,结果表明,爆破开挖扰动下,岩体含水裂纹的扩展,与爆炸应力波强度及其入射角、地应力的大小与方向、孔隙水压大小、裂纹的倾角及断裂韧度等因素相关;爆炸应力波的作用,相当于增大了岩体裂纹中的孔隙水压力,每1 cm/s的峰值振动速度相当于增大100 kPa的孔隙水压力,爆破振动速度越大,所产生的爆破扰动荷载越大;岩体开挖引起的岩体裂纹近区地应力及其孔隙水压力的变化,对裂纹的失稳与扩展具有较复杂的影响,可改变裂纹的失稳扩展模式。     

Modeling the connectivity and intersection of hydraulically loaded cracks with in situ fractures in rock

The interaction of an advancing hydraulically loaded crack and in situ fracture network can yield highly complex patterns. We model the connectivity of cells in a finite element domain and in a fracture network by a simplicial complex data structure. The complete adjacency information between cells is determined by one level down facet and one level up cofacet neighborhood information. Combined with a disjoint set data structure, explicit algorithms are derived to efficiently track network connectivity and load transfer between independent fracture sets. We also propose an approach to regularize the application of hydraulic load to newly intersected in situ cracks to smoothen the transition of pressure on intersected cracks from ambient to hydraulic pressure and to avoid the sudden loading of the entire length of these cracks. Numerical results demonstrate the performance of crack connectivity and load transfer models and the effect of regularization model. The results show that as the angle between an incoming hydraulically loaded crack and an in situ crack increases, the effect of in situ crack shifts from slight realignment to diversion/offsetting of the loaded crack. As the angle difference approaches the normal angle, the loaded crack tends to penetrate through the in situ crack. The proposed schemes are also used for transient simulation of 2D reservoirs with multiple perforations surrounded by in situ cracks with and without a bias in the distribution of their orientation. It was shown that from 2 perforations with angles closer to in situ cracks at low loading rates to all perforations at higher loading rates can result in active hydraulic crack propagation. The h‐adaptive method of asynchronous space‐time discontinuous Galerkin method is used to exactly track complex fracture patterns in these dynamic fracture simulations.     

钢衬钢筋混凝土压力管道外包混凝土的裂缝控制研究

基于钢衬钢筋混凝土压力管道在设计荷载下外包混凝土将带裂缝工作的基本承载特性,以缓解开裂与限裂间的矛盾为指导思想,将减小坝后背管外包混凝土裂缝宽度和提高开裂管道的耐久性作为研究目的,结合模型试验和有限元法,研究了一些裂缝控制措施的应用效果,以及有效的计算内水压力下钢衬外包混凝土裂缝宽度的公式。结果表明：总用钢量相同时,减薄钢衬,加大钢筋用量,减小裂缝宽度效果明显;减薄钢衬外包混凝土厚度,也可减小裂缝宽度;性能可靠的防水涂料具备黏接性强和变形能力高的基本性质,可进一步研究其对提高管道耐久性的长期作用;采用П-780-83规范中的裂缝宽度公式,计算结果与模型试验结果较为符合。     

Numerical simulation on the cracking and failure law of compacted clay lining in landfill closure cover system

According to the mass conservation equation of solid phase, water transport, air migration and enthalpy balance equations in a deformable and unsaturated porous medium, and the mechanics model, a thermal‐hydro‐mechanical (THM) coupling mathematical model is developed in order to accurately predict the crack propagation of the compacted clay lining (CCL) in a landfill closure cover system. Numerical simulation analysis of the cracking and failure of the CCL in a landfill closure cover system under complex conditions (landfill gas pressure, variation of soil surface temperature, moisture loss of coverage soil, and uneven settlement of landfilled MSW (municipal solid waste)) is performed. The results show that the cracks are discovered firstly on the soil surface and then extend to the deeper cover system and became bigger. The cracks had broken through the closure cover system at 30 days. Over the depth range of 0.06 m, the change rate of settlement is found to be the highest because of the strong influence of temperature. The rise of temperature around the cracks is fast when the CCL is broken through by the cracks. There is a hollow area near the cracks over the depth range from 0.15 to 0.7 m, and the temperature stabilize at around 11°C without any obvious reduction. Under the condition of evaporation, the main crack in the CCL is similar to a pumping well, which causes more rapid moisture content reduction around the cracks. Vertical deformation increases linearly with depth in the process of deformation and cracking of CCL induced by uneven settlement of landfilled MSW. Copyright © 2014 John Wiley & Sons, Ltd.     

水平集和无网格耦合法在裂纹扩展中的应用

提出了一种模拟裂纹扩展的水平集和无网格耦合方法。由于水平集和无网格方法都是建立在离散节点上,因而可以很自然地实现耦合。在该方法中,两个在裂尖处相互正交的水平集不仅用于描述裂纹的几何形态和裂尖位置,而且用于建立无网格伽辽金法（简称EFGM）不连续近似函数中的Heaviside跳跃项和裂尖处的Westergaard扩展项。当裂纹扩展时,则由水平集更新算法确定新裂纹的位置。水平集和无网格耦合法无需使用可视法、衍射法或透明法,克服了这些方法在裂尖处人为引入的不连续且能很好地再生 奇异场;而且节点影响域不受裂纹线切割的影响,在计算中往往使用较小的影响域,保持了整体刚度矩阵的带状、稀疏性;另外,水平集简化了扩展节点的选取和附加函数的建立,其更新过程无需求解演化方程,实现简单且易于编程。数值算例表明本文方法具有较高的计算精度,其模拟的裂纹扩展路径与试验结果吻合得很好,从而验证了本文方法的正确性和可行性。     

基于高密度电阻率法的土体干缩裂隙动态发育过程精细监测研究

土体在气候作用下发育干缩裂隙是一种常见的自然现象,裂隙的存在会极大弱化土体的工程性质,诱发许多岩土与地质工程问题。为了实时掌握黏性土中干缩裂隙网络的发展状态,提出了一套基于高密度电阻率层析成像技术（ERT）的土体干缩裂隙动态发展过程精细监测方法。分别开展模型试验及原位试验,利用自行研制的测定系统持续采集电流-电位差数据,随后利用自行开发的有限元法电阻率层析成像（FemERT）系统进行数据处理,获取了裂隙网络在不同发育阶段的空间分布特征。结果表明：（1）ERT可以实现土体裂隙发育过程的精细监测,具备监测三维裂隙网络几何形态的能力,裂隙宽度的识别精度达到毫米级,裂隙深度的识别精度达到厘米级;（2）ERT的感度分布特征解释了裂隙发育对于土体电阻率的影响规律,测定电阻值时程曲线因裂隙产生位置的不同而呈现不同的变化规律;（3）反演电阻率及其相对变化率（Rev）可以直观表征裂隙网络在不同阶段的空间几何形态,凸显裂隙动态发育过程对于土体导电性的影响。     

不同围压时含孔洞模型裂缝扩展的 连续—非连续数值模拟

为了有效地模拟连续介质向非连续介质的转化,发展了一种拉格朗日元方法、变形体离散元方法及虚拟裂纹模型耦合且考虑四边形单元沿对角线开裂的连续—非连续方法。利用该方法,模拟了不同围压时含孔洞模型的变形—开裂过程,统计了拉裂缝及剪裂缝区段数目,监测了一些单元的最大主应力。研究结果表明:当围压较小时,初始拉裂缝首先出现在孔洞的顶、底部,然后向模型的上、下端扩展,在初始拉裂缝的左、右两侧的拉应力集中区中产生远场拉裂缝,随后在孔洞的左、右两侧出现剪裂缝,最后,剪裂缝贯穿模型;当围压较大时,远场拉裂缝数量较少,未充分发展,远场拉裂缝与剪裂缝的发展阶段的界限不分明。含孔洞模型的最大承载力的下降是由于孔洞左、右两侧的剪裂缝向外扩展造成的。随着围压的增加,开始出现初始拉裂缝的时步数目增大,初始拉裂缝两侧的远场拉裂缝数目变少、出现变晚。      

基于局部单元劈裂的FEM/DEM自适应分析方法

为了模拟岩体中裂纹的萌生、扩展,Munjiza提出了有限元法/离散元法(FEM/DEM)耦合分析方法。因为裂纹是沿单元边界进行扩展的,亦即裂纹扩展具有网格依赖性,为获得较好的裂纹扩展形态,需要划分密集的初始网格。为解决上述难题,基于FEM/DEM耦合分析方法,提出了基于局部单元动态劈裂的FEM/DEM自适应分析方法,以克服裂纹扩展形态对网格的依赖性。该方法在最初建模时无需划分很密的初始网格,随着荷载的施加,对裂纹尖端附近的局部单元进行动态劈裂,为裂纹的后续扩展提供了更多可能的扩展方向,使得裂纹扩展不必沿着初始网格的单元边界扩展,即可以沿着单元内部进行扩展,裂纹扩展形态更为平滑,与实际情况更为接近。同时相对原FEM/DEM耦合分析方法一开始就划分很密的网格而言,新方法可以划分较为稀疏的初始网格,计算成本降低。最后,通过巴西劈裂算例与原FEM/DEM耦合分析方法对比,分析表明,新方法在一定程度上克服了裂纹扩展形态对初始网格的依赖性。     

A three-dimensional heat transfer and thermal cracking model considering the effect of cracks on heat transfer

A simple three-dimensional heat transfer model is developed to consider the hindering effect of cracks on heat transfer. The 3D heat transfer model can also be applied to numerical methods such as the combined finite-discrete element method (FDEM), discrete element method (DEM), discontinuous deformation analysis (DDA), the numerical manifold method (NMM), and the finite element method (FEM) to construct thermo-mechanical coupling models that allow these methods to solve thermal cracking problems and dynamically consider the hindering effect of cracks on heat transfer. In the 3D heat transfer model, the continuous-discontinuous medium is discretized into independent tetrahedral elements, and joint elements are inserted between adjacent tetrahedral elements. Heat transfer calculations for continuous-discontinuous media are converted to heat conduction in tetrahedral elements and the heat exchange between the adjacent tetrahedral elements through the joint element. If the joint element between adjacent tetrahedral elements breaks (ie, a crack generates), the heat exchange coefficient of the joint element is reduced to account for the hindering effect of cracks on heat conduction. Then the model and the FDEM are combined to build a thermo-mechanical coupling model to simulate thermal cracking. The thermally induced deformation, stress, and cracking are investigated by the thermo-mechanical coupling model, and the numerical results are compared with analytical solutions or experimental results. The 3D heat transfer model and thermo-mechanical model can provide a powerful tool for simulating heat transfer and thermal cracking in a continuous-discontinuous medium.     

红层泥质胶结物干缩裂缝深度与间距预测模型

红层泥质胶结物在干燥条件下会因含水率变化产生开裂,裂缝进一步发展贯通最终导致红层泥岩崩解破坏。目前对红层泥岩崩解机制的研究多为定性描述,缺乏定量分析方法。因此,定量研究红层泥质胶结物开裂的产生及扩展对深化红层泥岩崩解理论研究具有重要意义。基于线弹性断裂力学及非饱和土力学,建立了考虑含水率变化的红层泥质胶结物裂缝深度及间距计算公式,并对公式参数敏感性进行了分析。应用CT试验数据对PFC2D数值模型进行了校准,并基于干缩开裂数值试验结果,对裂缝深度及间距计算公式进行了对比验证。结果表明：提出的计算公式能较好地预测红层泥质胶结物的开裂深度及间距,裂缝深度随含水率减小呈先快速增加再缓慢增长的特点。研究成果有助于加深对红层泥岩崩解的量化理论分析,并为预测工程中红层泥岩的开裂程度提供参考。     

岩石剪切细观开裂演化与贯通机制分析

利用自主研发的煤岩剪切细观试验装置,以砂岩为研究对象,开展了不同饱和度数条件下的岩石剪切细观开裂演化与贯通机制试验研究,分析了砂岩剪切细观开裂演化过程和贯通机制,探讨了饱水系数对细观开裂和裂纹形态特征的影响。研究结果表明：剪切过程中,伴随着应力降,裂纹一般首先出现在剪切面中部,几条裂纹间断分布,距离较近的裂纹相互连通形成大裂纹,在应力继续下降的一瞬间,各裂纹连通形成贯通面;水的弱化作用影响了岩石的强度,使得试件底部出现较大的次级裂纹;饱水系数越大,分叉裂纹越多,使得试件的表面裂纹形态更加复杂,断裂面弯折段的数目越多,下端部次级裂纹越发育,整体损伤范围越大;细观上,局部独立裂纹贯通过程伴随着胶结颗粒岩桥的梁式旋转和弯曲折断,水的弱化作用,使得岩桥内部和外部微裂纹增多;随着饱水系数的增加,岩石剪切过程中出现的独立和次级裂纹的数目均呈增加趋势。     

Digital spatial cracking behaviors of fine-grained sandstone with precracks under uniaxial compression

X-ray computed tomography (CT) imaging and digital image correlation techniques are applied to study spatial cracking behaviors of sandstone under uniaxial compression, in which the angle between precracks is 45°, 90°, and 135° and the crack depth is 7.5 mm and 10 mm, respectively. Layered anisotropy damages and spatial cracking evolution are quantitatively analyzed by the defined digital layered anisotropy index and digital damage ratio, respectively. Three cases with different array of precracks evidence the depth effects of precracks on spatial crack propagation. Results show that the failure process of samples is first controlled by the coalescence of surface cracks in 2D space and then the samples are failed by the propagation of coalesced cracks (shear cracks with different shapes). The crack types for samples with precrack depth of 7.5 mm are all shear cracks for Cases 1‑3. Nevertheless, the crack types for samples with precrack depth of 10 mm are, respectively, the half X-shape crack for Case 1, X-shape crack for Case 2, and double shell crack for Case 3. The precrack has a significant promotion effect on the failure process when the angle between the two precracks is β = 90°, and the precrack has little to no effect on the failure process when the angle between the two precracks is β = 135°. As the depth of precrack increases to 10 mm, the crack types are changed in this study. The peak strength of sample subjected to uniaxial compression decreases with increasing depth of precracks, implying the decrease of the rock strength by the discontinuity.     

A two‐dimensional coupled hydromechanical discontinuum model for simulating rock hydraulic fracturing

Within the framework of our discontinuous deformation analysis for rock failure algorithm, this paper presents a two‐dimensional coupled hydromechanical discontinuum model for simulating the rock hydraulic fracturing process. In the proposed approach, based on the generated joint network, the calculation of fluid mechanics is performed first to obtain the seepage pressure near the tips of existing cracks, and then the fluid pressure is treated as linearly distributed loads on corresponding block boundaries. The contribution of the hydraulic pressure to the initiation/propagation of the cracks is considered by adding the components of these blocks into the force matrix of the global equilibrium equation. Finally, failure criteria are applied at the crack tips to determine the occurrence of cracking events. Several verification examples are simulated, and the results show that this newly proposed numerical model can simulate the hydraulic fracturing process correctly and effectively. Although the numerical and experimental verifications focus on one unique preexisting crack, because of the capability of discontinuous deformation analysis in simulating block‐like structures, the proposed approach is capable of modeling rock hydraulic fracturing processes. Copyright © 2014 John Wiley & Sons, Ltd.     

页岩水力压裂过程的态型近场动力学模拟研究

水力压裂技术广泛使用于页岩气开采工程中。为了分析压裂过程中多裂缝扩展形成复杂裂缝网的机制,尝试将态型近场动力学理论引入页岩水平井水力压裂过程的力学建模与数值仿真,在物质点间相互作用力模型中加入等效水压力项以实现在新生裂缝面上跟踪施加水压力,建立了水力压裂过程的近场动力学分析模型。通过模拟页岩储层的水力压裂过程,可得到复杂的裂缝扩展路径、裂缝网络的形成过程以及裂缝扩展受射孔间距及页岩天然裂缝和层理的影响。研究结果表明:射孔间距过小会造成起裂干扰,使中心射孔的裂缝扩展受到抑制;在压裂压力一定的情况下适当增大射孔间距,可以显著增强页岩压裂形成裂缝网的能力;压裂过程中水平层理面可能张开形成水平裂缝,且天然裂缝会诱导形成更复杂的垂直裂缝。模型和方法可为页岩水力压裂过程和机制研究及工程实践提供参考。     

双向压缩条件下闭合裂隙岩体断裂破坏机制及渗透压环境判定准则

在研究双向压缩条件下压剪复合型裂纹应力分布特征及断裂破坏机制基础上,考虑渗透压对初始裂隙面上有效正应力的影响,提出高低渗透压环境的判定准则,并基于滑动裂纹模型理论及最大周向拉应力破坏准则,得到不同渗透压环境下初始裂隙尖端微裂纹起裂特征与规律。研究结果表明：压剪复合应力条件下,初始裂隙尖端发育微裂纹的最优倾角与裂隙面摩擦系数直接相关,随裂隙面摩擦系数的增大,最优初始裂隙倾角由45°起逐渐增大;低渗透压条件下,渗流场的存在使裂纹面摩擦系数发生弱化,进而使得最优初始裂隙倾角向45°靠近,而渗透压直接降低裂隙面上有效正应力且与裂隙倾角无关,其仅仅影响裂隙体材料的初裂强度;高渗透压条件下,初始裂隙面由压剪复合应力状态转化为拉剪复合应力状态,并在拉剪复合应力场作用下,尖端微裂纹起裂角随KI/KII的不断增大,由70.5°逐渐趋近于0°。     

Polygonal grain‐based distinct element modeling for mechanical behavior of brittle rock

The microstructure of rock was numerically reproduced by a polygonal grain‐based model, and its mechanical behavior was examined by performing the uniaxial compression test and Brazilian tests via the Universal Distinct Element Code. The numerical results of the model demonstrated good agreement with the experimental results obtained with rock specimens in terms of the stress–strain behavior, strength characteristics, and brittle fracture phenomenon. An encouraging result is that the grain‐based model‐Universal Distinct Element Code model can reproduce a low ratio of tensile to compressive strength of 1/20 to 1/10 without the need for an additional process. This finding is ascribed to the fact that the geometrical features of polygons can effectively capture the effects of angularity, finite rotation, and interlocking of grains that exist in reality. A numerical methodology to monitor the evolution of micro‐cracks was developed, which enabled us to examine the progressive process of the failure and distinguish the contribution of tensile cracking to the process from that of shear cracking. From the observations of the micro‐cracking process in reference to the stress–strain relation, crack initiation stress, and crack damage stress, it can be concluded that the failure process of the model closely resembles the microscopic observations of rock. We also carried out a parametric study to examine the relationships between the microscopic properties and the macroscopic behavior of the model. Depending on the micro‐properties, the model exhibited a variety of responses to the external load in terms of the strength and deformation characteristics, the evolution of micro‐cracks, and the post‐peak behavior. Copyright © 2016 John Wiley & Sons, Ltd.