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.     

Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.     

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.     

有限元法与伪谱法混合求解弹性波动方程

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.     

论中国盐溶角砾岩的特征和分类

中国的蒸发岩矿床非常丰富,与蒸发岩在成因上密切相关的盐溶角砾岩极其发育。本文简述了中国盐溶角砾岩的分布,总结了盐溶角砾岩的一般特征,划分出盐溶角砾岩的主要类型,探讨了盐溶角砾岩与盐泉的关系。     

Petrography and textural development of inorganic and biogenic lithotypes in a relict barite tufa deposit at Flybye Springs, NT, Canada

A relict mound of Holocene barite (BaSO₄) tufa underlies the Flybye Springs, a small, barium‐rich, cold sulphur spring system in the Northwest Territories of Canada. The tufa is composed of relatively pure barite with ≤0·34 wt% Ca²⁺ and ≤0·77 wt% Sr²⁺. The mound is made up of coated bubble, raft, undulatory sheet, stromatolitic, coated grain and detrital conglomerate barite tufa. Although previously unreported in barite, these lithotypes are akin to facies found in many carbonate spring deposits. Raft and ooid‐coated grain tufa was formed via ‘inorganic’ barite precipitation in spring water ponds and tributaries where rapid oxidation of sulphide to sulphate established barite supersaturation. Undulatory sheet tufa may have formed by the reaction of dissolved barium with sulphate derived from the oxidation of extracellular polysaccharide‐rich colloidal sulphur films floating in oxygenated, barite‐saturated spring water ponds. Coated bubble, oncoid‐coated grain and stromatolitic tufa with filamentous microfossils was formed in close association with sulphur‐tolerant microbes inhabiting dysoxic and oxygenated spring water tributaries and ponds. Adsorption of dissolved barium to microbial extracellular polysaccharide probably facilitated the development of these ‘biogenic’ lithotypes. Detrital conglomerate tufa was formed by barite cementation of microdetrital tufa, allochthonous lithoclasts and organic detritus, including caribou hair. Biogenic textures, organic artefacts and microfossils in the Flybye barite tufa have survived diagenetic aggradational recrystallization and precipitation of secondary cements, indicating the potential for palaeoecological information to be preserved in barite in the geological record. Similarities between the Flybye barite tufa and carbonate spring deposits demonstrate that analogous textures can develop in chemical sedimentary systems with distinct mineralogy, biology and physiochemistry.     

Occurrence and evolution of the Xiaotangshan hot spring in Beijing, China

Thermal groundwater occurs in bedrock aquifers consisting of the dolomite of the Wumishan Group of the Jixianin System and the Cambrian carbonate in the Xiaotangshan geothermal field near the northern margin of the North China Plain, China. The hot water in the geothermal field of basin-type discharges partly in the form of the Xiaotangshan hot spring under natural conditions. The hot water has TDS of less than 600 mg/L and is of Na·Ca-HCO₃ type. The geothermal water receives recharge from precipitation in the mountain area with elevation of about 500 m above sea level to the north of the spring. Thermal groundwater flows slowly south and southeast through a deep circulation with a residence time of 224 years estimated with the Ra–Rn method. The Xiaotangshan hot spring dried up in the middle of the 1980s owing to the increasing withdrawal of the hot water in the geothermal field in the past decades. The water level of the geothermal system still falls continually at an annual average rate of about 2 m, although water temperature changes very little, indicating that the recharge of such a geothermal system of basin-type is limited. Over-exploitation has a dramatic impact on the geothermal system, and reduction in exploitation and reinjection are required for the sustainable usage of the hot water.     

基于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)初步探讨了黄海暖流路径的季节和年际变异,并指出这种变异与北向风的强弱密切相关;(2)通过分析济州岛西侧海域混合水的去向,进一步确认了部分混合水绕济州岛运行,并进入济州海峡这一事实.同时,初步揭示进入黄海的混合水,即黄海暖流水,含有更多的东海陆架水成分.