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.     

Solution to the rotation of the elastic earth by method of rigid dynamics

Hamiltonian mechanics is applied to the problem of the rotation of the elastic Earth. We first show the process for the formulation of the Hamiltonian for rotation of a deformable body and the derivation of the equations of motion from it. Then, based on a simple model of deformation, the solution is given for the period of Euler motion, UT1 and the nutation of the elastic Earth. In particular it is shown that the elasticity of the Earth acts on the nutation so as to decrease the Oppolzer terms of the nutation of the rigid Earth by about 30 per cent. The solution is in good agreement with results which have been obtained by other, different approaches.     

Magnetotelluric investigations for imaging electrical structure of Garhwal Himalayan corridor,Uttarakhand, India

Magnetotelluric investigations have been carried out in the Garhwal Himalayan corridor to delineate the electrical structure of the crust along a profile extending from Indo-Gangetic Plain to Higher Himalayan region in Uttarakhand, India. The profile passing through major Himalayan thrusts: Himalayan Frontal Thrust (HFF), Main Boundary Thrust (MBT) and Main Central Thrust (MCT), is nearly perpendicular to the regional geological strike. Data processing and impedance analysis indicate that out of 44 stations MT data recorded, only 27 stations data show in general, the validity of 2D assumption. The average geoelectric strike, N70°W, was estimated for the profile using tensor decomposition. 2D smooth geoelectrical model has been presented, which provides the electrical image of the shallow and deeper crustal structure. The major features of the model are (i) a low resistivity (<50Ωm), shallow feature interpreted as sediments of Siwalik and Indo-Gangetic Plain, (ii) highly resistive (> 1000Ωm) zone below the sediments at a depth of 6 km, interpreted as the top surface of the Indian plate, (iii) a low resistivity (< 10Ωm) below the depth of 6 km near MCT zone coincides with the intense micro-seismic activity in the region. The zone is interpreted as the partial melting or fluid phase at mid crustal depth. Sensitivity test indicates that the major features of the geoelectrical model are relevant and desired by the MT data.     

基于Poulos弹性理论的被动桩改进算法

传统的Poulos弹性理论仅适合于均质土中土体侧向位移时桩的性状分析,无法考虑土的层状特性。通过引入层状地基中作用一水平集中力的广义Mindlin解和地面作用有竖向荷载时的应力和位移通解,对Poulos方法进行了改进,使之扩展到多层土中,还用于研究堆载条件下的被动桩变形和受力响应。算例分析表明,改进弹性理论要比Poulos方法更为严密、合理,提高了计算精度,应用范围也更广。     

Reduction of hydroelastic responses of a very-long floating structure by a floating oscillating-water-column breakwater system

The hydroelastic responses of a very-long floating structure (VLFS) placed behind a reverse T-shape freely floating breakwater with a built-in oscillating water column (OWC) chamber are analyzed in two dimensions. The Bernoulli–Euler beam equation is coupled with the equations of rigid and elastic motions of the breakwater and the VLFS. The interaction of waves between the floating rigid breakwater and the elastic VLFS is formulated in a consistent manner. It has been shown numerically that the structural deflections of the VLFS can be reduced significantly by a suitably designed reverse T-shape floating breakwater.     

海底剖面声阻抗分布声遥测模拟实验研究

本文介绍一种声遥测海底平面层系声阻抗分布的模拟实验结果。实验在油漆厂冷却水池中进行,采用取样平均迭代的反演方法,对三介质层系的遥测结果与实测值较好地接近。     

小波边缘分析建模波阻抗反演方法在煤层解释中的应用

受限于地震波主波长的1／4,常规地震处理难以分辨较薄煤层。基于模型的波阻抗反演虽然克服了厚、薄煤层计算问题．但因难以得到准确的高频部分,很多情况下反演结果的精度和分辨率往往不能满足地质预测的要求。小波边缘分析建模波阻抗反演方法是从地震数据中提取地震特征参数,用于建模并参与迭代反演。该方法即可以弥补反演过程中高频成分的误差,减少对钻井数据的依赖程度,得到较为合理的初始模型,又可提高反演的精度和分辨率．使反演结果更好的反映实际地下地质情况。利用该方法,对大井-将军庙及新疆红沙泉2个勘探区的侏罗系中统西山窑组含煤地层进行波阻抗反演,反演结果与钻井资料进行比较,其中大井一将军庙0802及0809号孔单煤层厚度误差分别为3％、8％．新疆红沙泉多煤层厚度及煤层之间的夹层厚度预测结果除B6煤层误差较大外,其它误差均小于7％．可见该反演方法的计算结果比较准确。     

Reflection of waves from the boundary of a random elastic semi-infinite medium

In this paper the smooth perturbation technique is employed to investigate the problem of reflection of waves incident on the plane boundary of a semi-infinite elastic medium with randomly varying inhomogeneities. Amplitude ratios have been obtained for various types of incident and reflected waves. It has been shown that an incidentSH orSV type of wave gives rise to reflectedSH, P andSV waves, the main components beingSH andP, SV in the respective cases. The reflected amplitudes have been calculated depending upon the randomness of the medium to the square of the small quantity , where measures the deviation of the medium from homogeneity. An incidentP-type wave produces mainly aP component and also a weakSH component to the order of ². The reflected amplitudes obtainable for elastic media are also altered by terms of the same order. The direction of the reflected wave is influenced by randomness in some cases.     

The elastic properties of single-crystal fayalite as determined by dynamical measurement techniques

We present new elasticity measurements on single-crystal fayalite and combine our results with other data from resonance, pulse superposition interferometry, and Brillouin scattering to provide a set of recommended values for the adiabatic elastic moduliC _ij and their temperature variations. We use a resonance method (RPR) with specimens that were previously investigated by pulse superposition experiments. The nineC _ij of fayalite are determined from three new sets of measurements. One set of our newC _ij data is over the range 300–500 K. We believe that the relatively large discrep ncies found in someC _ij are due in large part to specimen inhomogeneities (chemical and microstructural) coupled with differences in the way various techniques sample, rather than only systematic errors associated with experimental procedures or in the preparations of the specimens.Our recommendeaC _ij's (GPa) and (C _ij/T) _p (GPa/K) are: The resulting values for the isotropic bulk and shear moduli,K _s and , and their temperature derivatives are:K _s=134(4) GPa; =50.7(0.3) GPa; (K _s/T) _p =–0.024(0.005) GPa/K; and (/T) _p =–0.013(0.001) GPa/K. An important conclusion is thatK _s increases as the Fe/(Fe+Mg) ratio in olivine is increased.     

注水煤体应力能 V-f 转化与释放机理研究

文章分析了在煤层注水湿润后煤体弹性能的转化与消耗过程，给出在湿润状态下能量由体变弹性能向形变弹性能转化进而自耗释放的理论模型。指出在坚硬干燥煤体中，注水后在水压力和毛细力的切割作用下，裂隙大量扩展，丰富了自由裂隙，产生更多的裂隙启发(裂隙切口)，使原来近似的刚性煤体变成由裂隙分割的块体组成的组合体或半组合体。裂隙中水的粘滞作用及水分子吸附力使块体问粘性结合、刚度降低、形变增大。煤体形变弹性能自身消耗释放的形式有3种：产生剪应力导致自由裂隙问错动摩擦；强度降低，使裂隙启发被劈裂成自由裂隙；湿润煤体蠕变性增强，产生整体性移动作功。初步给出能量消耗的量化表达式。通过这一理论，能建立煤层注水湿润度与冲击地压灾害发生的更直接的关系，使进一步量化该关系成为可能。从更深层次解释煤层注水防治冲击地压的机理。