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.     

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.     

Morphology and Luminosity Segregations in Galaxy Groups

We study morphology and luminosity segregation of galaxies in groups. We analyze the two catalogs of (∼2×400) groups which have been identified in the Nearby Optical Galaxy sample, by means of hierarchical and percolation `friends-of-friends' methods. We find that earlier-type (brighter) galaxies are more clustered and lie closer to the group centers, both in position and in velocity, than later-type (fainter) galaxies. Spatial segregations are stronger than kinematical segregations. These effects are generally detected at the ≳3-sigma level, with the exception of morphological segregation in velocity, which is the weakest effect. Our main results are confirmed by the analysis of statistically more reliable groups (with at least five members), and are strengthened by the detection of segregation in both hierarchical and percolation catalogs. Luminosity segregation is shown to be independent of morphology segregation. Our conclusions agree with a continuum of segregation properties of galaxies in systems, from low-mass groups to massive clusters. This revised version was published online in July 2006 with corrections to the Cover Date.     

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

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

The carbohydrates in relation to mineralogic and granulometric composition of surface sediments in the karst estuary (River Krka estuary,Yugoslavia)

The investigation of mineral, granulometric and chemical composition of sediments of the River Krka estuary (Yugoslavia) were performed in order to elucidate the origin of the sediments and the pattern of sedimentation. Estuarine surface sediments were found to be fine-grained with a bimodal distribution. Environmental conditions in estuarine sediments favour conservation of the organic matter (anoxic conditions). The carbohydrates in the sediments were investigated to determine whether they are of terrigenous or authigenous origin. Glucose, galactose, mannose, xylose, rhamnose, glucosamine and glucuronic acid were detected in the sediments. Their mutual relationship indicates a preferentially terrigenous source of sedimented organic material in estuarine sediments.     

Experimental solidification of anhydrous latitic and trachytic melts at different cooling rates: The role of nucleation kinetics

Two sets of cooling experiments were run at atmospheric conditions for two anhydrous starting latitic and trachytic melts: 1) five cooling rates (25, 12.5, 3, 0.5, and 0.125 °C/min) between 1300° and 800 °C, and 2) a 0.5 °C/min cooling rate from 1300 °C with quench temperatures at 1200°, 1100°, 1000° and 900 °C. Trachytic run-products are invariably glassy. Nucleation is also suppressed in the latitic run-products at the three highest cooling rates. Conversely, in the 0.5 and 0.125 °C/min runs, latites have a crystal content of  90 vol.%. The phases are: plagioclase, clinopyroxene, glass and iron-bearing oxide (in order of abundance). The variable quench temperatures, investigated by coupling experiments with Pt wire and Pt capsule sample containers in set 2, again did not produce crystallization of trachyte, whereas latitic samples are characterized by  10 vol.% of oxides, pyroxenes and plagioclase (in order of appearance), at temperature < 1000 °C. Effects of (preferential) heterogeneous nucleation on sample holders, of superheating degree, and chemical species loss during cooling are absent for both melt compositions. The difference of solidification paths between these two silicate melts can be ascribed only to their small chemical differences. In comparison with calculated equilibrium conditions all the experimental latitic and trachytic run-products revealed strong kinetic effects, interpretable in the light of the nucleation theory. The glass-forming ability (GFA) of trachyte is higher, whereas their critical cooling rate (Rc) is lower (< 0.125 °C/min), in comparison to latitic melts (Rc > 0.5 °C/min). The experimental results carried out in this study can be applied to lava flows and domes; trachytic lavas are able to flow for longer period with respect to latitic ones in a metastable condition. Glass-rich terrestrial lavas, i.e. obsidians, can be the result of sluggish nucleation kinetics due to the relative high polymerisation of evolved silicate melts.     

Plant distribution at the mobile dune scale and its relevance to soil properties and topographic features

A survey was conducted in an 11-year recovery mobile dune (RMD11) and a 20-year recovery mobile dune (RMD20), in Horqin Sandy Land, Northern China, to determine plant distribution at the mobile dune scale and its relevance to soil properties and topographic features. The results showed that (1) vegetation cover and species number increased from dune top to bottom in the restoration process of mobile dune; (2) the average value of soil organic C, total N, pH, relative height of sampling site, very fine sand content and soil water contents (40−60 and 60−80 cm) of RMD11 were less than that of RMD20, respectively, and there were significant differences (P < 0.05) between the two dunes; (3) soil resources were redistributed by shrub restoration and relative height of sampling site on dune. The distribution of sand pioneer plant, Agriophyllum squarrosum, was positively related to the relative height of sampling site and soil water content, while that of other herbaceous plants was positively related to soil nutrients in the restoration process of mobile dune. These results suggest that at mobile dune scale, plant distributions are determined by a combination of soil properties and topographic feature. Much effort should be made to preserve the interdune lowland and to improve the level of soil nutrients on mobile dune.     

The relationship between hydrodynamic properties and weathering of soils derived from volcanic rocks,Galapagos Islands (Ecuador)

The aim of this interdisciplinary study is to examine a component of the hydrological cycle in Galapagos by characterizing soil properties. Nine soil profiles were sampled on two islands. Their physical and hydrodynamic properties were analyzed, along with their mineralogical composition. Two groups of soils were identified, with major differences between them. The first group consists of soils located in the highlands (>350 m a.s.l.), characterized by low hydraulic conductivity (<10⁻⁵ m s⁻¹) and low porosity (<25%). These soils are thick (several meters) and homogeneous without coarse components. Their clay fraction is considerable and dominated by gibbsite. The second group includes soils located in the low parts of the islands (<300 m a.s.l.). These soils are characterized by high hydraulic conductivity (>10⁻³ m s⁻¹) and high porosity (>35%). The structure of these soils is heterogeneous and includes coarse materials. The physical properties of the soils are in good agreement with the variations of the rainfall according to the elevation, which appears as the main factor controlling the soil development. The clayey alteration products constrain soils physical and hydrodynamic properties by reducing the porosity and consequently the permeability and also by increasing water retention.     

基于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.