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.     

深水铺管船托管架疲劳分析及动态监测点选取研究

作为S型铺管作业的关键性装备,托管架长期承受着交变载荷的作用,随着海洋开发向超深水发展,结构疲劳破坏问题不容忽视。采用疲劳谱分析的方法并结合线性累积损伤理论,对托管架结构频域下的变形进行了分析,计算了正常海况下和极端海况下托管架疲劳损伤度并对疲劳寿命进行了评估。研究发现托管架在正常海况下作业符合安全要求,在极端海况下局部结构会受到破坏。并从托管架结构安全监测角度,筛选了结构疲劳分析关键点位,为监测点位的选取提供了依据。     

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.     

Source rates and ion recycling rates for Na and K in Mercury's atmosphere

The supply rates of Na and K to the atmosphere of Mercury by processes acting on the extreme surface—thermal vaporization, photon-stimulated desorption (PSD), and ion-sputtering—are limited by the rates at which atoms can be supplied to the extreme surface by diffusion from inside the regolith grains. Supply rates to the atmosphere are further regulated by ion retention and by gardening rates that supply new grains to the surface. We consider the limits on supply of sodium and potassium atoms to the atmosphere, and rates of photoion recycling to the surface. Thermal vaporization rates are severely limited by the ability of atoms to diffuse to the surface of the grain. Therefore, the diffusion-limited thermal vaporization rates on Mercury's surface are comparable to or less than the PSD rates. Ion sputtering is primarily due to highly ionized heavy ions, even though they represent a small fraction of the solar wind. We have shown that up to 60% of the Na photoions are deposited on the surface of Mercury. Ion recycling to the surface can have a long-term effect on the regolith abundance if an average recycling pattern persists such that more ions return to a particular area than are launched there. It is unknown whether the formation of latitude bands of >100% ion retention persist on average despite a rapidly changing magnetosphere. The total exospheric column of sodium observed at Mercury between 1997 to 2003 varied by a factor of 2-3 from perihelion to aphelion.     

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.     

Outflow channel sources, reactivation, and chaos formation, Xanthe Terra, Mars

The undulating, warped, and densely fractured surfaces of highland regions east of Valles Marineris (located north of the eastern Aureum Chaos, east of the Hydraotes Chaos, and south of the Hydaspis Chaos) resulted from extensional surface warping related to ground subsidence, caused when pressurized water confined in subterranean caverns was released to the surface. Water emanations formed crater lakes and resulted in channeling episodes involved in the excavation of Ares, Tiu, and Simud Valles of the eastern part of the circum-Chryse outflow channel system. Progressive surface subsidence and associated reduction of the subsurface cavernous volume, and/or episodes of magmatic-driven activity, led to increases of the hydrostatic pressure, resulting in reactivation of both catastrophic and non-catastrophic outflow activity. Ancient cratered highland and basin materials that underwent large-scale subsidence grade into densely fractured terrains. Collapse of rock materials in these regions resulted in the formation of chaotic terrains, which occur in and near the headwaters of the eastern circum-Chryse outflow channels. The deepest chaotic terrain in the Hydaspis Chaos region resulted from the collapse of pre-existing outflow channel floors. The release of volatiles and related collapse may have included water emanations not necessarily linked to catastrophic outflow. Basal warming related to dike intrusions, thermokarst activity involving wet sediments and/or dissected ice-enriched country rock, permafrost exposed to the atmosphere by extensional tectonism and channel incision, and/or the injection of water into porous floor material, may have enhanced outflow channel floor instability and subsequent collapse. In addition to the possible genetic linkage to outflow channel development dating back to at least the Late Noachian, clear disruption of impact craters with pristine ejecta blankets and rims, as well as preservation of fine tectonic fabrics, suggest that plateau subsidence and chaos formation may have continued well into the Amazonian Period. The geologic and paleohydrologic histories presented here have important implications, as new mechanisms for outflow channel formation and other fluvial activity are described, and new reactivation mechanisms are proposed for the origin of chaotic terrain as contributors to flooding. Detailed geomorphic analysis indicates that subterranean caverns may have been exposed during chaos formation, and thus chaotic terrains mark prime locations for future geologic, hydrologic, and possible astrobiologic exploration.     

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.