Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

Cranny density parameters and porosity measured by elastic wave method in quasi-isotropic cranny rock masses

According to the characteristic of elastic waves propagation in medium and the application of elastic waves method in rock mass engineering, the cranny mass with random crannies was regarded as quasi-isotropic cranny mass. In accordance with the rock rupture mechanics, principle of energy balance and Castiglano＇s theorem, the relationship of effective dynamic parameters of elasticity （E, v, G） and cranny density parameters or porosity was put forward. On this basis, through the theory of elastic waves propagation in isotropic medium, the relationship between the elastic wave velocity and cranny density parameters and porosity was set up. The theoretical research results show that, in this kind of cranny rock masses, there is nonlinear relationships between the effective dynamic parameters of elasticity and wave velocities and the cranny density parameter or porosity; and with the increase of cranny density parameter or porosity of cranny rock masses, the effective dynamic modulus and the elastic wave velocities of cranny rock masses will decrease; and at the same time, when the cranny density parameter or porosity is very small, the effective dynamic modulus of elasticity and the elastic wave velocities change with the cranny density parameter, which can explain the sensitivity of effective elastic parameters and elastic wave velocities to cranny rock masses.

Experimental Study on Mechanical Characteristics of Cracked Rock Mass Reinforced by Bolting and Grouting

The stress hardening characteristics of the reinforced rock mass in uniaxial compression tests were revealed by means of the experimental study on mechanical characteristics of cracked rock mass reinforced by bolting and grouting. And the load-beating mechanism of the reinforced rock mass was perfectly reflected by the experiment. The results can offer some useful advice for support design and stability analysis of deep drifts in unstable strata.

ZONAL DISINTEGRATION MECHANISM OF DEEP CRACK-WEAKENED ROCK MASSES UNDER DYNAMIC UNLOADING

Size and quantity of fractured zone and non-fractured zone are controlled by cracks contained in deep rock masses. Zonal disintegration mechanism is strongly dependent on the interaction among cracks. The strong interaction among cracks is investigated using stress superposition principle and the Chebyshev polynomials expansion of the pseudo-traction. It is found from numerical results that crack nucleation, growth and coalescence lead to failure of deep crack- weakened rock masses. The stress redistribution around the surrounding rock mass induced by unloading excavation is studied. The effect of the excavation time on nucleation, growth, interaction and coalescence of cracks was analyzed. Moreover, the influence of the excavation time on the size and quantity of fractured zone and non-fractured zone was given. When the excavation time is short, zonal disintegration phenomenon may occur in deep rock masses. It is shown from numerical results that the size and quantity of fractured zone increase with decreasing excavation time, and the size and quantity of fractured zone increase with the increasing value of in-situ geostress.

Deformation analysis of transversely isotropic coal-rock mass with porous and cracks

Coal-rock as a typical sedimentary rock has obvious stratification,namely it has transversely isotropic feature.Meanwhile,deformation leads to coal-rock mass having the characteristics of different porous and crack structures as well as local anisotropy.Equivalent axial and circumferential strain' formulas of the pure coal-rock mass specimen with a single crack were derived through the establishment of equivalent mechanical model of standard cylindrical coal-rock specimen,and have been widely used to a variety of media combined different structures containing multiple cracks.The complete stress strain curve of a real coal-rock specimen was obtained by the CTC test.Additionally,according to the comparison with the theoretical value,the theoretical mechanical model could well explain the deformation characteristics of coal-rock mass and verify its validity.Further,following features were analyzed:strain normalized coefficient and elastic modulus(Poisson's ratio) in vertical and parallel direction to the stratification,stratification angle,porosity,pore radius,normal and tangential stiffness of crack,and the relationship of different crack width with different tangential stiffness of crack.Through the analysis above,it substantiate this claim that the theoretical model with better reliability reflects the transversely isotropic nature of the coal-rock and the local anisotropy caused by the porous and cracks.

THE CONSTITUTIVE RELATION OF CRACK-WEAKENED ROCK MASSES UNDER AXIAL-DIMENSIONAL UNLOADING

An accurate and efficient numerical method for solving the crack-crack interaction problem is presented. The method is mainly by means of the dislocation model, stress superposition principle and Chebyshev polynomial expansion of the pseudo-traction. This method can be applied to compute the stress intensity factors of multiple kinked cracks and multiple rows of periodic cracks as well as the overall strains of rock masses containing multiple kinked cracks under complex loads. Many complex computational examples are given. The dependence of the crack-crack interaction on the crack configuration, the geometrical and physical parameters, and loads pattern, is investigated. By comparison with numerical results under confining pressure unloading, it is shown that the crack-crack interaction under axial-dimensional unloading is weaker than those under confining pressure unloading. Numerical results for single faults and crossed faults show that the single faults are more unstable than the crossed faults. It is found from numerical results for different crack lengths and different crack spacing that the interaction among kinked cracks decreases with an increase in length of the kinked cracks and the crack spacing under axial-dimensional unloading.

Coordination and evolution mechanism of fractal crack network in mining rock mass destruction

Based on the natural characters of stratum, complicated geological mining conditions and the essence of mining rock mass destruction, the complexity of rock mass destruction caused by miningw as analyzed. The inner link between rock mass destruction phenomena caused by mining and nonlinear science was revealed. There are numerous cracks in natural rock mass. The cracks’ distribution is irregular and is of statistical fractal structure. Self-organizational nonlinear evolution of the inner structure flaws leads to the rock mass destruction with external force. The evolution includes single fault’s fractal development, formation and evolution of fractal crack network and coordination of fractal crack network, etc. The law of fractal crack network’s evolution was introduced, at the same time, the coordination of fractal crack network was analyzed. Finally, based on coordination the principal equation of mining-caused subsidence of structural rock mass was established and its steady-state solution and unsteady-state solution were found.

Self-catalyzed Effect and Cracking Risk in Mass Concrete Containing Micro-slag

The main results obtained from the experimental and engineering investigation on the heat evolution and cracking risk of a furnace concrete block were presented. The heat evolution of experimental mortars containing micro-slag under different environmental temperatures was instrumented in order to investigate the self-catalyzed effect, which was discovered in engineering. More-over,the thermal stress of the furnace concrete due to heat temperature rise was calculated to evaluate the cracking risk of mass concrete containing micro-slag due to self-catalyzed effect. The experimental results illustrate that with the development of hydration and initial temperature of mixture, the hydra-tion can be also accelerated and temperature of concrete will be continued to rise, which was the self-catalyzed effect. And the thermal stress due to self-catalyzed effect could not result in the cracking of furnace concrete.

Equivalent elastic compliance tensor for rock mass with multiple persistent joint sets:Exact derivation via modified crack tensor

Discontinuities constitute an integral part of rock mass and inherently affect its anisotropic deformation behavior.This work focuses on the equivalent elastic deformation of rock mass with multiple persistent joint sets.A new method based on the space geometric and mechanical properties of the modified crack tensor is proposed,providing an analytical solution for the equivalent elastic compliance tensor of rock mass.A series of experiments validate the capability of the compliance tensor to accurately represent the deformation of rock mass with multiple persistent joint sets,based on conditions set by the basic hypothesis.The spatially varying rules of the equivalent elastic parameters of rock mass with a single joint set are analyzed to reveal the universal law of the stratified rock mass.

Application and prospects of 3D printing in physical experiments of rock mass mechanics and engineering:materials,methodologies and models

The existence of joints or other kinds of discontinuities has a dramatic efect on the stability of rock excavations and engineering.As a result,a great challenge in rock mass mechanics testing is to prepare rock or rock-like samples with defects.In recent years,3D printing technology has become a promising tool in the feld of rock mass mechanics and engineering.This study frst reviews and discusses the research status of traditional test methods in rock mass mechanics tests of making rock samples with defects.Then,based on the comprehensive analysis of previous research,the application of 3D printing technology in rock mass mechanics is expounded from the following three aspects.The frst is the printing material.Although there are many materials for 3D printing,it has been found that 3D printing materials that can be used for rock mass mechanics research are very limited.After research,we summarize and evaluate printing material that can be used for rock mass mechanics studies.The second is the printing methodology,which mainly introduces the current application forms of 3D printing technology in rock mass mechanics.This includes printed precise casting molds and one-time printed samples.The last one is the printing model,which includes small-scale samples for mechanical tests and large-scale physical models.Then,the benefts and drawbacks of using 3D printing samples in mechanical tests and the validity of their simulation of real rock are discussed.Compared with traditional rock samples collected in nature or synthetic rock-like samples,the samples made by 3D printing technology have unique advantages,such as higher test repeatability,visualization of rock internal structure and stress distribution.There is thus great potential for the use of 3D printing in the feld of rock mass mechanics.However,3D printing materials also have shortcomings,such as insufcient material strength and accuracy at this stage.Finally,the application prospect of 3D printing technology in rock mass mechanics research is proposed.

Fracture behaviors of columnar jointed rock mass using interface mechanics theorem

For a special geological structure of columnar jointed rock mass(CJRM),its mechanical properties are strongly affected by the columnar joints.To describe the fracture behaviors of CJRM using the basic theories of interface mechanics for composite materials,the interface stresses of the vertical and horizontal joints,which are the two primary joints in the CJRM under triaxial compression,are studied,and their mathematical expressions are derived based on the superposition principle.Based on the obtained interface stresses of the vertical and horizontal joints in the CJRM,the crack initiation of the joint interface in the CJRM is studied using the maximum circumferential stress theory in fracture mechanics.Moreover,based on this investigation,the fracture behaviors of CJRM are analyzed.According to the results of similar material physical model tests for the CJRM,the theoretical study is verified.Finally,the influence of the mechanical parameters of the CJRM on the joint interface stress is discussed comprehensively.

Spatial Thermal Crack Control in Mass Concrete

The finite element software,MIDAS is used to predict the distribution of temperatures and,analyzes the cracking control methods within a hydrating mass concrete.The temperature control of mass concrete has great significance in assuring the project quality.Adiabatic or semi adiabatic temperature measurement is mostly used for measuring and controlling the temperature fluctuation during construction.The temperature distribution produced by the finite element thermal analysis of the model is used to quantify the maximum allowable internal temperature difference before crack initiation on concrete.This study analyzes the data from one high-rise structure project in Shanghai are used to verify the finite element model developed.Results suggest that reliance on a limiting maximum temperature differential to control cracking in massive concrete applications should be supplemented with a requirement for analysis showing the calculated spatial temperature and stress response to the predicted temperature distribution within the concrete,to ensure that the induced tensile stresses will not exceed the tensile strength of the concrete and so minimize the risk of having thermal cracks at early age.

冰充填裂隙对冻结岩体压缩破坏特征的影响研究

冻结岩体的力学性质是决定冻结地层工程施工与运营安全的核心因素,而冻结岩体的力学性质由裂隙结构和冻结温度等条件决定。本文在常温和冻结条件下开展了含预制裂隙砂岩的力学实验,通过高速摄影技术对试样表面的裂纹起裂和扩展过程进行了监测,基于实验结果着重讨论了冰充填裂隙对冻结岩体压缩破坏特征的影响机制。结果表明:(1)冻结没有改变试样强度随裂隙倾角的变化趋势,随裂隙倾角的增大,试样的强度均呈先减后增的趋势。强度最小值出现在裂隙倾角α=60°时,最大值出现在α=90°时;(2)在冻结状态下,试样的起裂模式随裂隙倾角的增加整体呈现“拉伸起裂→剪切起裂→拉伸起裂”的变化趋势,此外,相同裂隙倾角的试样在冻结状态下起裂的应力水平均高于在常温状态下;(3)冻结裂隙砂岩的力学性质受冻结作用和裂隙倾角双重控制。并进一步分析了冰填充裂隙对冻结岩体压缩破坏特征的影响,认为冰均起到了支撑、填充和黏结作用。在裂隙倾角较小时(0°~30°),裂隙冰的法向压力较高,强化效应主要来源于冰的支撑作用;当裂隙倾角较大(45°~75°)时,裂隙冰上的法向压力较低,强化效应主要来源于冰-岩界面黏结作用;裂隙倾角(近90°)与加载方向平行时,裂隙冰的法向应力变为拉应力,强化效应主要来源于冰-岩界面黏结作用。

基于数字图像相关方法的煤岩表面应变场及裂纹演化特征

基于数字图像相关方法,开展了煤岩体试样单轴压缩破坏的裂纹演化特征及破坏过程研究,分析了试样破坏过程的表面应变场、变形局部化带的拉伸位移演化和错动位移演化特征,精细化分析煤岩体变形破坏过程。研究结果表明:(1)根据试样的应力-应变曲线特征,将试样分为缓慢破坏型和瞬间破坏型。缓慢破坏型试样的破坏过程较为平稳,裂纹宽度较小。瞬间破坏型试样破坏过程较为剧烈,表面裂纹宽度增加,表现出明显的岩爆现象。(2)加载过程中,缓慢破坏型试样先产生应变集中后,裂纹沿变形局部化带扩展延伸并逐渐贯通;瞬间破坏型试样在应力峰值后试样表面瞬间产生宏观裂纹。(3)试样表面应变演化特征与加载过程中应力变化特征对应,反映了裂纹扩展规律。缓慢破坏型试样在加载过程中表面应变变化过程较为缓慢,并在后期呈现波动状态;瞬间破坏型试样在应力峰值点处,拉伸应变瞬间突增达到峰值,导致试样破坏。

破碎岩体隧道洞口开挖诱导滑坡与古滑坡耦合分析

以大栗树高速公路隧道洞口施工为工程背景,研究了在考虑古滑坡因素下洞口工作面边坡失稳及其加固技术。利用钻孔测斜数据(BCD)和Midas GTS折减强度有限元分析模型,分别推测了最危险滑动面,得知两者滑动面位置相似。并采用不平衡推力传递法,基于BCD推断滑动面计算安全系数。以BCD和有限元法的滑动面为基础,将滑动面近似假设为对数螺旋线,考虑地形因素,建立考虑多台阶边坡上限分析模型,研究了不同加固措施对边坡稳定的影响。结果表明,3种方法计算出的边坡安全系数结果相近。

超宽鱼腹式现浇箱梁混凝土裂缝控制

近年来,国内站城综合体TOD开发成为一种新型建筑模式,地铁、市政道路、商业综合体等多种大体量建筑形态交叉同步施工成为工程实践难题。基于此,结合已经实施的成都博览城综合交通枢纽工程,解析福州路桥鱼腹式多幅四箱室闭合型超宽箱梁,采取优化混凝土配合比、U形冷却管布置、后浇带布置、支架沉降控制、箱梁浇筑顺序等5种措施,解决大体积C60梁体混凝土防开裂难题,并取得明显实践成果,获得国家优质工程奖和土木工程詹天佑奖。

地铁车站大体积混凝土裂缝产生原因及解决措施分析

地铁建设近年来流行于各大中城市,地铁车站建设工程大体积混凝土施工比较普遍,在混凝土凝固过程中由于其体积收缩以及各种压力的作用,会造成大量的墙体裂缝;大面积的裂缝会导致建筑渗漏,给地铁车站带来了很大威胁。本文从荷载、位移、温度三个方面分析了地铁车站大体积混凝土结构的开裂成因;根据裂缝的大小、深度、成因等几个角度对裂缝进行了划分;从设计、材料、施工三个角度,对如何防止混凝土开裂进行了探讨。本文研究成果对保障地铁车站大体积混凝土施工质量有一定的参考意义。

BFRP网格对大体积混凝土早期温度裂缝控制研究

黄茅海跨海通道索塔在早期养护中表面开裂现象严重,索塔是一种典型的大体积混凝土,其早期出现的开裂现象大部分是由较高的温度应力导致的,即温度裂缝。因此,为了对早期温度裂缝进行控制,同时考虑到海洋环境的腐蚀情况,工程采用BFRP网格对索塔表面进行裂缝控制,控裂效果显著,开裂现象被有效缓解。BFRP网格虽然成功应用于大体积混凝土的温度裂缝控制,但是,目前BFRP网格控裂机理尚不清晰。基于此,采用ABAQUS有限元软件模拟BFRP网格对大体积混凝土中的温度裂缝控制过程,其中裂缝是通过扩展有限元(XFEM)实现。模拟结果表明:BFRP网格可以有效控制大体积混凝土温度裂缝面积、体积以及开裂单元数量。随着BFRP网格网孔的减小,温度裂缝面积、体积以及开裂单元数量进一步下降,网孔为5 cm的BFRP网格将裂缝控制为仅有1条。BFRP网格可以有效控制裂缝主要是由于BFRP网格是一种双向受力结构并且具有强节点,在与混凝土黏结之后,发挥了协调混凝土应力发展的作用,有效地均匀化混凝土的应力分布,进而均匀化裂缝发展,最终达到减少裂缝的效果。

江阴靖江长江隧道现浇大体积混凝土的抗裂技术与应用

在原材料优选和配合比优化的基础上,通过协同调控水化速率与膨胀历程,设计制备了低温升高抗裂混凝土(M0不掺抗裂剂、M1和M2掺抗裂剂),对比研究了其力学性能、绝热温升和自生体积变形。在此基础上,依托江阴靖江长江隧道现浇大体积C40混凝土工程,采用基于水化-温度-湿度-约束多场耦合机制的收缩模型对该工程不同结构部位的开裂风险进行了定量评估。试验结果表明:掺入适量抗裂剂能有效改善混凝土的绝热温升和自生体积变形,且对抗压强度影响不大。开裂风险评估计算结果表明:底板、侧墙可选择M1配合比,并控制入模温度不大于28℃,同时在侧墙中心布置冷却水管;顶板可选择M2配合比,并控制倒角处和其他部位的入模温度分别不大于28℃和32℃。工程现场持续跟踪监测6个月,未发现混凝土出现开裂、渗漏,达到预期目标。

氧化镁膨胀剂对大体积混凝土裂缝控制与应用

城市副中心综合交通枢纽工程存在底板浇筑方量大、施工气温较高、大体积混凝土浇筑等难点,采用M型氧化镁膨胀剂配制补偿收缩混凝土,通过原材料的选择、优化配合比设计、控制绝热温升及施工养护等措施,形成整体性裂缝控制技术,可较大程度降低混凝土开裂风险,保障工程质量。