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.     

Micromechanical approach to effective viscoelastic properties of micro‐fractured geomaterials

The aim of this paper is to formulate a micromechanics‐based approach to non‐aging viscoelastic behavior of materials with randomly distributed micro‐fractures. Unlike cracks, fractures are discontinuities that are able to transfer stresses and can therefore be regarded from a mechanical viewpoint as interfaces endowed with a specific behavior under normal and shear loading. Making use of the elastic‐viscoelastic correspondence principle together with a Mori‐Tanka homogenization scheme, the effective viscoelastic behavior is assessed from properties of the material constituents and damage parameters related to density and size of fractures. It is notably shown that the homogenized behavior thus formulated can be described in most cases by means of a generalized Maxwell rheological model. For practical implementation in structural analyses, an approximate model for the isotropic homogenized fractured medium is formulated within the class of Burger models. Although the approximation is basically developed for short‐term and long‐term behaviors, numerical applications indicate that the approximate Burger model accurately reproduce the homogenized viscoelastic behavior also in the transient conditions.     

西秦岭温泉花岗岩体岩石学特征及岩浆混合标志

温泉花岗岩体由酸性端元的寄主岩石和暗色微细粒镁铁质包体群及基性岩墙群组成。无岩浆混合作用或岩浆混合作用较弱区段，寄主岩石以似斑状二长花岗岩为主．显示正常的花岗岩结构构造岩浆混合作用强烈区段。岩石的异常结构构造十分发育．矿物之间自形程度差异显著．常见包晶反应、包含结构、交代边、熔蚀边、交代蚕食的港湾状结构构造及交代缝合线、矿物镶边、斜长石异常环带和矿物残留等，多见指示岩浆混合的标志性矿物针状磷灰石。暗色微粒包体中多见寄主二长花岗岩中的捕掳晶。包体的形态、结构构造以及与寄主岩石强烈地成分交换等均是岩浆混合作用的标志。     

Smeared crack approach: back to the original track

This paper briefly reviews the formulations used over the last 40 years for the solution of problems involving tensile cracking, with both the discrete and the smeared crack approaches. The paper focuses on the smeared approach, identifying as its main drawbacks the observed mesh‐size and mesh‐bias spurious dependence when the method is applied ‘straightly’. A simple isotropic local damage constitutive model is considered, and the (exponential) softening modulus is regularized according to the material fracture energy and the element size. The continuum and discrete mechanical problems corresponding to both the weak discontinuity (smeared cracks) and the strong discontinuity (discrete cracks) approaches are analysed and the question of propagation of the strain localization band (crack) is identified as the main difficulty to be overcome in the numerical procedure. A tracking technique is used to ensure stability of the solution, attaining the necessary convergence properties of the corresponding discrete finite element formulation. Numerical examples show that the formulation derived is stable and remarkably robust. As a consequence, the results obtained do not suffer from spurious mesh‐size or mesh‐bias dependence, comparing very favourably with those obtained with other fracture and continuum mechanics approaches. Copyright © 2006 John Wiley & Sons, Ltd.     

浙西北江山─绍兴断裂带构造演化特征

江山—绍兴断裂带是中国南方最重要的构造带之一。它是江南陆块与华夏陆块神功运动形成的碰撞拼贴带。它表现出明显的多期变形、构造叠加的特征：神功运动形成近东西向平卧褶皱；印支运动沿断裂带有向北西强烈的推覆挤压。断裂带不同时期的活动反映了不同特征的构造环境，形成不同的构造变形特征，并对两侧的地质历史具有明显的控制作用。     

GRADATION AND ASSESSMENT OF DAMAGE DEGREE OF NATURAL LANDSCAPE IN XING, AN COUNTY, GUANGXI

A new technique designed to help quantify the degree of damage to the landscape from one area to another shows a close relationship between population density and the degree of landscape damage. The technique establishes a scale of damage from 0 to 5 (zero = no damage; 5 = severe damage) using data from aerial photographs, land-use maps, and field data. The related formula allows one to compare the relative degree of damage across regions using a combination of an absolute index, a theoretical index, a relative index, and population density. Xing'an County is used to demonstrate the technique.     

甘新蒙北山地区成矿带划分和基本特征对比

在对甘新蒙北山地区成矿地质背景和成矿在时空上的联系性等综合研究分析基础上,运用新理论、新认识,以板块构造单元划分为基础,将本区划分为3个Ⅲ级成矿带、11个Ⅳ级成矿带,并对Ⅲ级成矿带特征进行了阐述,以期为矿产勘查和地质大调查工作部署提供新的思路和依据。进一步就主要矿床(点)、构造单元、矿种、赋矿建造、控矿因素、成矿类型、成矿时代等方面进行了对比,认为各Ⅳ级成矿带在区域上具有相对稳定的展布范围,具有相同的构造环境和成矿条件、类似的矿种和成矿时代,存在一定的相似性、可比性,并且各成矿带具有一定的潜在找矿价值。     

Damage identification method using harmonic excitation force considering both modelling and measurement errors

A method of structural damage identification using harmonic excitation force is presented. It considers the effects of both measurement and modelling errors in the baseline finite element model. Damage that accompanies changes in structural parameters can be estimated for a damaged structure from the change between measured vibration responses and ones calculated from the analytical model of the intact structure. In practice, modelling errors exist in the analytical model due to material and geometric uncertainties and a reduction in the degrees of freedom as well as measurement errors, making identification difficult. To surmount these problems, bootstrap hypothesis testing, which enables statistical judgment without information about these errors, was introduced. The method was validated by numerical simulation using a three‐dimensional frame structure and real vibration data for a three‐storey steel frame structure. Copyright © 2005 John Wiley & Sons, Ltd.     

THE CHARACTERISTIC OF MARINE ENVIRONMENT IN LINGDINGYANG ESTUARY

1 HYDROLOGIC FEATURES Lingdingyang Estuary, located at the middle south of Guangdong Province, is a bell-shaped estuary with a north-south direction. Its area is about 2100km2. The north of Qi′ao Island and Inner-Lingding Island, and the south of Humen are grouped as Neilingdingyang Estuary, having an area of 1041km2. Affected by topography, runoff and tide, its dynamic condition is very complicated. Different water areas have different hydrologic features. The topography under …