含单裂隙岩体蠕变破裂裂纹扩展演化规律研究

岩体内部存在大量不同尺寸裂隙，形成不同结构形式的裂隙岩体，裂隙的存在会降低岩体的力学性能，对于岩体工程的整体稳定性具有很大的影响。本研究基于滑动翼裂纹模型对恒定载荷下含单裂隙岩体蠕变变形进行了理论推导和分析，基于MAPLE软件直观给出恒定载荷下裂隙角度、长度对蠕变变形的影响，采用连续-非连续数值分析软件GDEM,对预制单裂隙在不同倾角、不同长度情况下岩体的蠕变破裂裂纹扩展进行数值模拟。通过理论分析和数值模拟可以看出，岩体蠕变速率随着裂隙长度的增大而增大；随着裂隙角度的增加出现先增大后减小的趋势，其中裂隙倾角为30°时蠕变速率最大。含单裂隙岩体蠕变破裂裂纹扩展规律的研究可为岩体工程蠕变破裂规律研究及蠕变治理提供一定的理论依据和参考。     

非贯通裂隙岩体损伤演化率相关性及变形特征

含非贯通裂隙岩体是自然界中岩体的主要赋存形式，其裂隙几何特征对岩体的强度及变形均产生显著影响。应变率对岩体的损伤演化及黏滞效应也具有显著的率相关性。首先，运用模型元件的方法，将非贯通裂隙岩体动态破坏过程视为具复合损伤、静态弹性特性、动态黏滞特性的非均质点组成，对黏弹性响应的Maxwell体进行改进，将细观损伤体与裂隙损伤演化的宏观损伤体根据等效应变假设并联组成宏细观复合损伤体，构建综合考虑岩体宏细观缺陷的动态损伤模型；其次，基于断裂力学及应变能理论，对岩体宏观裂隙动态扩展的能量机制进行分析，综合考虑初始裂隙应变能、裂隙动态损伤演化过程应变能、裂隙闭合应变能，得到裂隙岩体宏观动态损伤变量计算公式；最后，将模型计算结果与实验结果进行比较，模型计算结果与实验结果吻合较好，证明了模型的合理性，同时利用模型讨论了裂隙倾角、应变率、岩石性质对岩体变形特征的影响规律。     

考虑裂隙间相互作用情况下围压卸荷过程应力应变关系

岩体的稳定性和变形特性主要决定于裂隙，同时裂隙间的相互作用对岩体的稳定和变形产生显著的影响。裂隙岩体在加载和卸荷条件下的力学特性有显著的区别。为此本文首次利用位错模型法结合叠加原理研究在围压卸荷条件下裂隙间的相互作用对岩体的变形的影响问题。文中推导了考虑裂隙间的相互作用情况下裂隙岩体围压卸荷过程的应力应变关系及应力强度因子表达式，且进行了数值计算。     

结构化岩体采动裂隙分布规律与分形性实验研究

采动岩体裂隙形成、扩展和分布均十分复杂，而初始结构面的存在使其更加复杂。本文借助相似材料模型实验，模拟了结构化岩体采动裂隙的分布规律，应用分形几何理论研究采动岩体裂隙分布的自相似性。研究表明：具有初始结构面的岩体采动后形成的裂隙网络具有很好的统计自相似分形性质，分形维数可以综合描述采动岩体裂隙化程度，且地表沉陷特征值与分形维数具有明确的数量对应关系。     

裂隙岩体强度试验单点法及其前景

本文在回顾裂隙岩体强度试验的基础上，用大量原位测试数据对裂隙岩体直剪破坏机理和裂隙岩体直剪破坏的基本规律及其成因进行了深入探讨；总结出一套行之有效的试验技术；给出了裂隙岩体直剪破坏时的充分必要条件并加以量化。     

综放顶煤爆破能量的分形研究

基于综放开采预爆破弱化坚硬顶工艺,利用分形理论建立岩石断耗能量与块度分维值之间的理论关系,推导出岩石爆破块度分数维与炸药单耗之间关系,进行了室内相似材料模型爆破实验,验证了上述理论关系,利用爆破岩石表面裂隙分布的分形维数值作了衡量指标,优化爆破孔网参数。     

不同应力差下非线性蠕变研究

深部工程处于复杂的环境场(应力场、温度场和渗流场)之中,引起的地下工程围岩流变破坏事故屡见不鲜。本文通过三轴蠕变实验,根据岩体流变力学和蠕变原理对影响深部岩体蠕变的因素进行了研究。用实验数据进行数值计算,得出了深部软岩巷道围岩在高地应力作用下的非线性蠕变模型,建立了应力差-时间-蠕变的定量关系,利用MATLAB软件得到了不同情况下的蠕变-时间曲线;应力差-时间-蠕变曲面。从而验证了非线性流变理论,找到了引起蠕变的主要因素应力差。并利用应力差理论分析了深部巷道围岩的蠕变破坏机理。     

裂隙岩体损伤局部化破坏分岔模型及其应用

采用概率统计方法分析节理裂隙岩体的几何特征，定义了反映裂隙岩体几何特征的组构张量．基于不可逆热力学理论，通过裂纹扩展的细观分析，得出了损伤的发展机理和演化方程，把损伤演化和裂隙的几何特征的变化联系起来，建立了弹塑性损伤本构关系．为分析含有节理裂隙岩体在发生局部化破坏时的特征，通过对发生局部化时的裂隙岩体的分析。构造了适用于节理裂隙岩体局部化分析的不连续分岔模型．利用非线性规划数值解法，可以得出局部化破坏的方向特征．在有限元方法中，根据该模型给出了节理裂隙岩体相关的算例，分析表明该模型用于分析裂隙岩体的局部化破坏是有效的．     

岩体系统渗流场与应力场耦合的广义双重介质模型的应用研究

本文通过实验资料分析, 提出了裂隙岩体应力与渗流之间的分形几何关系式。以研究区岩体系统结构为基础, 提出了岩体系统应力场与渗流场耦合的广义双重介质模型, 并将该模型应用于某水利水电工程的渗流与控制中, 对坝区岩体进行了应力场和渗流场的耦合分析, 取得了比较满意的结果。     

岩体风化的综合分带研究

岩体的风化不仅会引起岩石物质成份、化学成份变化 ,而且也引起岩体裂隙、岩体完整性的改变。因而风化岩体中裂隙数量、岩体块度的变化与岩体的风化程度具有较好的对应关系。论文研究了用岩体裂隙间距、岩体的完整性指标、岩石质量指标对岩体风化程度进行分带 ,较好地将岩体风化分带与岩体结构、岩体工程特性紧密地结合起来。在岩体风化程度的定量划分上作了新的探索。     

层状介质中双硬币形裂纹对扭转波的散射──远场解

本文在[1]的基础上,研究了多层介质中由于双硬币形裂纹引起的对弹性扭转波散射的远场特性.文中应用Hankel积分变换和Abel变换,将问题最后归结为求解一组第二类Fredholm积分方程,并导出了用积分形式给出的散射位移场表达式.最后运用围道积分技术和渐近分析的方法,对散射位移场在远离裂纹时的主要性态进行了分析讨论.     

The deformation and fracture of balloons

Inflation of balloons provides a straightforward way of achieving large biaxial deformations. Previous studies have shown that when a balloon bursts, crack propagation occurs at very high speed – much higher than would be expected from the low strain modulus and elastic wave velocity of the rubber. The present paper is concerned with studies of the deformation and fracture of cylindrical balloons. On inflation, the deformations of such a balloon pass through an unstable region but subsequently increase monotonically with pressure. In this relatively high pressure region, the ratio of the longitudinal and circumferential extension ratios is broadly in accord with expectations from high-strain elasticity theory when the ratio of the corresponding stresses is taken into account. On bursting, crack speeds up to around 300 m/s in this region. It is shown that these speeds are in accord with large increase in incremental moduli for the highly-strained rubber. Marked changes in crack tip profile observed at very high crack speeds are consistent with control of the rate of growth by inertia rather than by the viscoelastic properties of the rubber (as is believed to be the case at lower speeds). Consistent with this, various elastomers having different glass transition temperatures show similar crack growth behaviour in the very high speed region.     

Ultrasonic sizing of embedded vertical cracks in ferritic steel welds

Embedded vertical cracks are very often encountered in structural components necessitating reliable procedure for determining their size. Specialized ultrasonic techniques are required for estimation of through thickness height of these cracks [Mater. Eval. (1970) 28; J. NDI 26 (1977) 320; NND (1975) 146]. The conventional pulse-echo technique has serious limitations for quantitative estimation of vertical cracks because it is subject to variations due to reflectivity, coupling factors etc. The present study aims to improve the accuracy of embedded crack height measurement based on transit time rather than the reflected amplitude [NDT Int. (1982) 315; J. JSNDT 27 (1978) 118]. A unique method has been developed for generating known sizes of fatigue cracks ranging from 2.5 to 18 mm in height and 15 mm in length. These cracks were generated in rectangular steel specimens and subsequently embedded at predetermined locations using welding. These welded specimens containing the embedded vertical cracks were examined using ultrasonic techniques based on transit time like (a) flaw tip echo method and (b) mode conversion method. Using both the techniques excellent agreement was obtained between the estimated crack height and actual crack height with an accuracy of ±1 mm.     

Crack propagation in glass coatings under expanding spherical contact

The growth of transverse cracks under expanding spherical contact in a model system consisted of soda-lime glass bonded to a polycarbonate substrate is observed in situ from below or from the polished edge of the bilayer. Abrasion or chemical etching is employed on the coating surfaces to control the initial fracture. In the limit case of monoliths, the crack mouth becomes fully engulfed by the expanding contact, which results in a much steeper crack angle compared to the classical Hertzian cone case. As the coating thickness is reduced, flexure stresses are set in the coating which drive the cone crack to well away from the contact circle and initiate semi-elliptical-like radial cracks at the subsurface, right under the contact. Common to all three fracture modes is an initial unstable propagation phase following by a stable growth, with detrimental failure associated with severe damage to the top surface and/or delamination at the coating/substrate interface taking place at loads several times the fracture initiation loads.LEFM in conjunction with a large-strain FEM contact code is used to study the post-initiation fracture, with the crack path controlled by the principal stress trajectory or zero-mode II S.I.F. The analysis exposes the leading geometric and material parameters in each fracture mode, which may be useful in the design of bilayer structures for optimal mechanical performance. The well-known Auerbach law governing the initial fracture of monoliths is found to apply also to the bilayer crack systems within a certain range of the problem parameters. The numerical prediction for the crack profiles and the fracture envelopes generally collaborate well with the tests.     

First fundamental problems of anisotropic elastic plane weakened by periodic collinear cracks

Using the method of complex functions, we discuss the first fundamental problems of an anisotropic infinite elastic plane weakened by periodic collinear cracks and with periodic boundary loads on both sides of the cracks. This problem was considered by Cai [Engineering Fracture Mechanics 46（1）, 133-142 （1993）]. However, the previous method is imperfect. Therefore, the results are incorrect. Here, we revise the method and give a correct solution.     

The non-linear dynamic response of the Euler-Bernoulli beam with an arbitrary number of switching cracks

In this study the non-linear dynamic response of the Euler-Bernoulli beam in presence of multiple concentrated switching cracks (i.e. cracks that are either fully open or fully closed) is addressed. The overall behaviour of such a beam is non-linear due to the opening and closing of the cracks during the dynamic response; however, it can be regarded as a sequence of linear phases each of them characterised by different number and positions of the cracks in open state. In the paper the non-linear response of the beam with switching cracks is evaluated by determining the exact modal properties of the beam in each linear phase and evaluating the corresponding time history linear response through modal superposition analysis. Appropriate initial conditions at the instant of transition between two successive linear phases have been considered and an energy control has been enforced with the aim of establishing the minimum number of linear modes that must be taken into account in order to obtain accurate results. Some numerical applications are presented in order to illustrate the efficiency of the proposed approach for the evaluation of the non-linear dynamic response of beams with multiple switching cracks. In particular, the behaviour under different boundary conditions both for harmonic loading and free vibrations has been investigated.     

Multiple cracking of elastic coatings under transient thermal load

A periodic array of cracks in an elastic coating bonded to a homogeneous substrate is considered. The medium is subjected to mechanical loads and/or thermal loads. The problem is formulated in terms of a singular integral equation with the crack face displacement as the unknown variable. In addition to the time-varying stress intensity factors and stresses for various parameters of the problem, the effect of periodic cracking on the relaxation of the transient stress on the coating surface is discussed. Solution techniques for a single elastic layer and an elastic coating bonded to an infinite substrate are given. It is found that, if the crack density attains a saturation value, the transient thermal stress in the medium could be released significantly, suggesting that further cracking is difficult.     

Parameters of the Edge Effect as a Function of the Mechanical Characteristics of a Composite with an Interfacial Crack

The paper studies the dependence of edge effects in a laminated composite of regular structure with a periodic system of symmetric interfacial cracks on the mechanical characteristics of the composite components. Loading the composite induces constant strain in the reinforcement direction. The problem is solved in exact formulation using the linear elastic equations for piecewise-homogeneous media and a criterion for quantitative evaluation of edge effects. An approximate solution of the problem is found by the mesh approach. A difference scheme for the mixed problem under consideration is derived. The edge-effect zones for normal stresses are constructed. The maximum values of the normal stress concentration factors are determined     

Influence of Young's Moduli of a Laminated Composite with a Periodic System of Cracks on Edge Effects

A problem on edge effects in a laminated composite with a periodic system of symmetric contact cracks is studied. Loading of the composite induces constant deformation in the reinforcement direction beyond the edge effect. The problem is solved in an exact formulation under the linear theory of elasticity for piecewise-homogeneous media. An approximate solution is found by the mesh approach using the concept of a base scheme. For stresses _ij , the shape and size of the edge effect zone are studied depending on the ratio of Young's moduli of the composite components, which varies over a wide range     

用样条边界元与奇异积分方程法联合求解多裂纹柱体的扭转

通过间解的分离，本文将径向多裂纹柱体的导曲函两个调和函数表示，使问题归为解一组混混合型积分方程。针对方程的特点，本文联合使用三次样条边界法与奇异积分方程的数值方法对所得方程建立了数值法，并对裂纹相交情形作了特殊处理。最后对工程中感兴趣的一些典型的多裂纹柱体的扭转作了例题计算，结果表明，本文方法具有收敛快，精度高的特点。