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.     

Age and early metamorphic history of the Sanbagawa belt: Lu–Hf and P–T constraints from the Western Iratsu eclogite

Two distinct age estimates for eclogite-facies metamorphism in the Sanbagawa belt have been proposed: (i) c.  120–110 Ma based on a zircon SHRIMP age for the Western Iratsu unit and (ii) c.  88–89 Ma based on a garnet–omphacite Lu–Hf isochron age from the Seba and Kotsu eclogite units. Despite the contrasting estimates of formation ages, petrological studies suggest the formation conditions of the Western Iratsu unit are indistinguishable from those of the other two units—all ∼20 kbar and 600–650 °C. Studies of the associated geological structures suggest the Seba and Western Iratsu units are parts of a larger semi-continuous eclogite unit. A combination of geochronological and petrological studies for the Western Iratsu eclogite offers a resolution to this discrepancy in age estimates. New Lu–Hf dating for the Western Iratsu eclogite yields an age of 115.9 ± 0.5 Ma that is compatible with the zircon SHRIMP age. However, petrological studies show that there was significant garnet growth in the Western Iratsu eclogite before eclogite facies metamorphism, and the early core growth is associated with a strong concentration of Lu. Pre-eclogite facies garnet (Grt1) includes epidote–amphibolite facies parageneses equilibrated at 550–650 °C and ∼10 kbar, and this is overgrown by prograde eclogite facies garnet (Grt2). The Lu–Hf age of c.  116 Ma is strongly skewed to the isotopic composition of Grt1 and is interpreted to reflect the age of the pre-eclogite phase. The considerable time gap ( c.  27 Myr) between the two Lu–Hf ages suggests they may be related to separate tectonic events or distinct phases in the evolution of the Sanbagawa subduction zone.     

Optimal Filtering Algorithm for Stochastic 2-D FMM Ⅱ with Multiplicative Noise

1　Introduction　Inrecentyearstherehasbeengrowingresearchin terestinthetwo dimensional (2 D )systemtheory .The 2 Dsystemsmayfindapplicationsinareassuchasmarineseismicdataprocessingandimageprocessing .Althoughmoreandmorevaluableresultshavebeengained ,mos…     

居民收入与通信网络市场经济效益研究--以河北省通信网络为实例

利用Logisitic曲线模型和扩展支撑模型对居民可支配收入预测和单项消费预测，以河北省居民的收入及其消费支出为实例，研究了城镇居民和农民的可支配收入分配和单项边际消费倾向，合理建立了预测模型，并理论计算了1999～2003年的通信消费方面的边际消费倾向、预测了2004～2008年通信消费情况。最后，根据预测结果对通信网络建设提出了一些有价值的建议。     

Kuiper带小天体动力演化

Kuiper带是指太阳系内位于离太阳 30～ 50AU一个区域。 1 992年该区域陆续发现了一群半径在几十到几百公里的小天体。这些小天体在Kuiper带的分布是极其不均匀的。Kuiper带小天体的发现对人们认识太阳系的形成与演化有重要的意义。本文回顾了近年来国际上在Kuiper带小天体动力演化方面的研究 ,着重分析了目前国际上几种用以解释其非均匀分布的动力学机制 ,并提出目前该领域的一些尚未解决的问题。     

A high resolution,multiproxy Late‐glacial record of climate change and intrasystem responses in northwest England

A lacustrine carbonate sequence from Hawes Water, Lancashire, UK, has been studied using stable isotopic, lithological, pollen and mineral magnetic analysis. The data reveal four abrupt climatic oscillations in the Late‐glacial Interstadial leading up to the onset of the Loch Lomond Stadial. The data also point to climatic warming relatively early within the stadial, ca. 12 500 GRIP yr, prior to the onset of the Holocene. The oxygen isotope record is taken as a signature of climate forcing against which the response of the lake‐system can be monitored. By adopting this approach it is revealed that the response of the biological system to the rapid climatic oscillations is non‐linear and primarily a function of the antecedent conditions. A significant end‐Devensian isotopic excursion (A) is matched by only minor changes in the cold‐adapted floras and faunas. During the warmer interstadial, the response of the biological ecosystem (events B–D) is clearly influenced by thresholds: major changes in the catchment vegetation associated with relatively minor oscillations in the isotopic signature. The stratigraphical patterns reveal significant lag effects between the onset of climate deterioration and resulting changes in vegetation. Copyright © 2002 John Wiley & Sons, Ltd.     

流褶层及韧变带研究

流褶层与韧变带是地壳拉伸变形，顺层固态流变作用下的产物。流褶层是以原始层理为变形面或再经递进变形的褶皱变形岩层或岩石共生组合层位。韧变带具明显的层控性，受岩石成分和应变程度控制，不同环境和不同成分岩石的韧变带具有相异的组合型式和变形机制，井具有一定的递变规律。流褶层和韧变带可分属不同层位，但流褶层可实现向韧变带的转化。     

康古尔塔格缝合带的变形机制

北天山东段康古尔塔格带是晚古生代塔里木板块和准噶尔板块碰撞的结果。它是一条复杂的、强烈的高应变带．并具有独特的变形机制、应变序列以及构造变形。本文运用构造－地层研究方法对该碰撞带的构造特征加以分析和研究。     

Tidal Breakup of Comets

In the context of the survival of periodic comets of different origins, rotational breakup and tidal disruption could be important, especially of the short period comets injected from the Kuiper belt. This is because long-period comets from the distant Oort cloud tend to be subject to thermal stress and volatile 'explosion' far more severely. A simple calculation using the Öpik method of random planetary close encounters was performed to estimate the probability of tidal disruption of comets and scattered Kuiper belt objects (SKBOs) during their orbital migration. It was found that a large fraction of the short period comets and SKBOs might have been internally fragmented by single or multiple close encounters with the outer planets.     

利用高光谱数据进行植被生化成分反演方法研究

高光谱数据包含着丰富光谱信息,能够定量地分析物质成分^[1]。由高光谱数据,可以运用多元逐步线性回归方法反演植被生化组分含量,从而达到监测植被生长状况的目的^[2]。本文具体介绍了多元逐步线性回归方法,及其在氮、木质素和纤维素含量估算中的应用。