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.     

Jurassic igneous rocks of the central Sanandaj–Sirjan zone (Iran) mark a propagating continental rift,not a magmatic arc

Jurassic igneous bodies of the Sanandaj–Sirjan zone (SaSZ) in SW Iran are generally considered as a magmatic arc but critical evaluation of modern geochronology, geochemistry and radiogenic isotopes challenges this conclusion. There is no evidence for sustained igneous activity along the ~1,200 km long SaSZ, as expected for a convergent plate margin; instead activity was brief at most sites and propagated NW at ~20 mm/a. Jurassic igneous rocks define a bimodal suite of gabbro‐diorite and granite. Chemical and isotopic compositions of mafic rocks indicate subcontinental lithospheric mantle sources that mostly lacked subduction‐related modifications. The arc‐like features of S‐type granites reflect massive involvement of Cadomian crust and younger sediments to generate felsic melts in response to mafic intrusions. We conclude that Jurassic SaSZ igneous activity occurred in a continental rift, not an arc. SaSZ igneous rocks do not indicate that subduction along the SW margin of Eurasia began in Jurassic time.     

Dynamic modelling of retrogressive landslides with emphasis on the role of clay sensitivity

This paper presents a detailed numerical study of the retrogressive failure of landslides in sensitive clays. The dynamic modelling of the landslides is carried out using a novel continuum approach, the particle finite element method, complemented with an elastoviscoplastic constitutive model. The multiwedge failure mode in the collapse is captured successfully, and the multiple retrogressive failures that have been widely observed in landslides in sensitive clays are reproduced with the failure mechanism, the kinematics, and the deposition being discussed in detail. Special attention has been paid to the role of the clay sensitivity on each retrogressive failure as well as on the final retrogression distance and the final run‐out distance via parametric studies. Moreover, the effects of the viscosity of sensitive clays on the failure are also investigated for different clay sensitivities.     

Planar deformation features and impact glass in inclusions from the Vredefort Granophyre,South Africa

Abstract— The Vredefort Granophyre represents impact melt that was injected downward into fractures in the floor of the Vredefort impact structure, South Africa. This unit contains inclusions of country rock that were derived from different locations within the impact structure and are predominantly composed of quartzite, feldspathic quartzite, arkose, and granitic material with minor proportions of shale and epidiorite. Two of the least recrystallized inclusions contain quartz with single or multiple sets of planar deformation features. Quartz grains in other inclusions display a vermicular texture, which is reminiscent of checkerboard feldspar. Feldspars range from large, twinned crystals in some inclusions to fine‐grained aggregates that apparently are the product of decomposition of larger primary crystals. In rare inclusions, a mafic mineral, probably biotite or amphibole, has been transformed to very fine‐grained aggregates of secondary phases that include small euhedral crystals of Fe‐rich spinel. These data indicate that inclusions within the Vredefort Granophyre were exposed to shock pressures ranging from <5 to 8–30 GPa. Many of these inclusions contain small, rounded melt pockets composed of a groundmass of devitrified or metamorphosed glass containing microlites of a variety of minerals, including K‐feldspar, quartz, augite, low‐Ca pyroxene, and magnetite. The composition of this devitrified glass varies from inclusion to inclusion, but is generally consistent with a mixture of quartz and feldspar with minor proportions of mafic minerals. In the case of granitoid inclusions, melt pockets commonly occur at the boundaries between feldspar and quartz grains. In metasedimentary inclusions, some of these melt pockets contain remnants of partially melted feldspar grains. These melt pockets may have formed by eutectic melting caused by inclusion of these fragments in the hot (650 to 1610 °C) impact melt that crystallized to form the Vredefort Granophyre.     

Pulsar Wind Nebulae in Egret Error Boxes

A remarkable number of pulsar wind nebulae (PWN) are coincident with EGRET γ-ray sources. X-ray and radio imaging studies of unidentified EGRET sources have resulted in the discovery of at least six new pulsar wind nebulae (PWN). Stationary PWN (SPWN) appear to be associated with steady EGRET sources with hard spectra, typical for γ-ray pulsars. Their toroidal morphologies can help determine the geometry of the pulsar which is useful for constraining models of pulsed γ-ray emission. Rapidly moving PWN (RPWN) with more cometary morphologies seem to be associated with variable EGRET sources in regions where the ambient medium is dense compared to what is typical for the ISM.     

What the Sunspot Record Tells Us About Space Climate

The records concerning the number, sizes, and positions of sunspots provide a direct means of characterizing solar activity over nearly 400 years. Sunspot numbers are strongly correlated with modern measures of solar activity including: 10.7-cm radio flux, total irradiance, X-ray flares, sunspot area, the baseline level of geomagnetic activity, and the flux of galactic cosmic rays. The Group Sunspot Number provides information on 27 sunspot cycles, far more than any of the modern measures of solar activity, and enough to provide important details about long-term variations in solar activity or “Space Climate.” The sunspot record shows: 1) sunspot cycles have periods of 131± 14 months with a normal distribution; 2) sunspot cycles are asymmetric with a fast rise and slow decline; 3) the rise time from minimum to maximum decreases with cycle amplitude; 4) large amplitude cycles are preceded by short period cycles; 5) large amplitude cycles are preceded by high minima; 6) although the two hemispheres remain linked in phase, there are significant asymmetries in the activity in each hemisphere; 7) the rate at which the active latitudes drift toward the equator is anti-correlated with the cycle period; 8) the rate at which the active latitudes drift toward the equator is positively correlated with the amplitude of the cycle after the next; 9) there has been a significant secular increase in the amplitudes of the sunspot cycles since the end of the Maunder Minimum (1715); and 10) there is weak evidence for a quasi-periodic variation in the sunspot cycle amplitudes with a period of about 90 years. These characteristics indicate that the next solar cycle should have a maximum smoothed sunspot number of about 145 ± 30 in 2010 while the following cycle should have a maximum of about 70 ± 30 in 2023.