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.     

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

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.     

The mass ratio of Charon to Pluto from Hubble Space Telescope astrometry with the fine guidance sensors

The mass ratio of Charon to Pluto is a basic parameter describing the binary system and is necessary for determining the individual masses and densities of these two bodies. Previous measurements of the mass ratio have been made, but the solutions differ significantly (Null et al., 1993; Young et al., 1994; Null and Owen, 1996; Foust et al., 1997; Tholen and Buie, 1997). We present the first observations of Pluto and Charon with a well-calibrated astrometric instrument—the fine guidance sensors on the Hubble Space Telescope. We observed the motion of Pluto and Charon about the system barycenter over 4.4 days (69% of an orbital period) and determined the mass ratio to be 0.122±0.008 which implies a density of 1.8 to 2.1 g cm⁻³ for Pluto and 1.6 to 1.8 g cm⁻³ for Charon. The resulting rock-mass fractions for Pluto and Charon are higher than expected for bodies formed in the outer solar nebula, possibly indicating significant postaccretion loss of volatiles.     

Spectroscopic studies of three novae

We present our spectroscopic observations of the novae V1425, V1493, and V1494 Aql carried out with the 125-cm telescope at the Crimean Station of the Sternberg Astronomical Institute in the wavelength range 4000–11000 Å. We measured the emission-line intensities, determined the nova shell expansion velocities from the line profile FWHMs and components, and estimated the interstellar reddening from the first members of the Balmer series. The chemical composition of the nova shells is analyzed. Nitrogen and oxygen were found to be overabundant in V1425 and V1494 Aql; the helium abundance turned out to be normal in the two stars.     

Pulsar magnetic field for an infinitely conductive accretion disk

The magnetic-field distribution outside a flat, infinitely conductive unbounded disk in the field of a point magnetic dipole is determined. A relationship is established between the problem of magnetic-field determination and the problem of the flow of an ideal incompressible fluid around an infinitely thin disk.