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.     

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

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年通信消费情况。最后，根据预测结果对通信网络建设提出了一些有价值的建议。     

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

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

Investigating the mechanisms leading to the deglaciation of past continental northern hemisphere ice sheets with the CLIMBER–GREMLINS coupled model

A coupling procedure between a climate model of intermediate complexity (CLIMBER-2.3) and a 3-dimensional thermo-mechanical ice-sheet model (GREMLINS) has been elaborated. The resulting coupled model describes the evolution of atmosphere, ocean, biosphere, cryosphere and their mutual interactions. It is used to perform several simulations of the Last Deglaciation period to identify the physical mechanisms at the origin of the deglaciation process. Our baseline experiment, forced by insolation and atmospheric CO₂, produces almost complete deglaciation of past northern hemisphere continental ice sheets, although ice remains over the Cordilleran region at the end of the simulation and also in Alaska and Eastern Siberia. Results clearly demonstrate that, in this study, the melting of the North American ice sheet is critically dependent on the deglaciation of Fennoscandia through processes involving switches of the thermohaline circulation from a glacial mode to a modern one and associated warming of the northern hemisphere. A set of sensitivity experiments has been carried out to test the relative importance of both forcing factors and internal processes in the deglaciation mechanism. It appears that the deglaciation is primarily driven by insolation. However, the atmospheric CO₂ modulates the timing of the melting of the Fennoscandian ice sheet, and results relative to Laurentide illustrate the existence of threshold CO₂ values, that can be translated in terms of critical temperature, below which the deglaciation is impeded. Finally, we show that the beginning of the deglaciation process of the Laurentide ice sheet may be influenced by the time at which the shift of the thermohaline circulation from one mode to the other occurs.     

CATASTROPHIC ECO-ENVIRONMENTAL CHANGE IN THE SONGNEN PLAIN, NORTHEASTERN CHINA SINCE 1900S

Although the Songnen Plain in the northeastern China was developed relatively late in the temperate zone of the world, its eco-environment has changed greatly. This paper analyzes the changes of land cover and the rates and trends ofdesertification during the past 100 years in the Songnen Plain. According to the macroscopic analysis, we find that the eco-environment in the plain has reached to the threshold of catastrophic change since the 1950s. The Thorn Needle Catastrophic Model was used to determine and validate this conclusion. Human activities, including large-scale construction projects, such as huge dams and dikes, and excessive grazing were the primary factors contributing to regional eco-environmental catastrophe. And irrational reclamation of the wilderness also affected the eco-environmental change. The results reveal the complex human-land interactions.     

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.     

One zone modeling of irregular variability of stellar convective envelope

One zone modeling of the irregular variability of red super-giants is intended with regard to the nonlinear coupling of finite amplitude pulsation with convection. The nonlocal mixing length is employed for the evaluation of the convective flux, the turbulent pressure and the turbulent power of temperature fluctuations. The radial pulsation and the Boussinesq convection are assumed for simplicity. The one zone is defined as the layer having the entropy maximum and the minimum at the bottom and at the top, respectively. The quasi-adiabatic approximation is consistent with this definition in fixing the zone to the same mass range. The spatial derivatives are evaluated under the assumption of homologous changes with the equilibrium homologous parameters. Then, a set of 6 simultaneous first order nonlinear ordinary differential equations are obtained as the one zone representation of the irregular variability of the convective envelope.     

A MESO-β SCALE SIMULATION OF THE EFFECTS OF BOREAL FOREST ECOSYSTEM ON THE LOWER ATMOSPHERE

Based on the Intensive Field Campaign(IFC-1)data of Boreal Ecosystem-Atmosphere Study(BOREAS).a three-dimensional meso-β scale model is used to simulate the effect of boreal forests onthe lower atmosphere.A fine horizontal resolution of 2 km×2 km is used in order to distinguish thevegetative heterogeneity in the boreal region.A total of 20×25 grid points cover the entire sub-modeling area in BOREAS' South Study Area(SSA).The ecosystem types and their coverage ineach grid square are extracted from the North American Land Cover Characteristics Data Base(NALCCD)generated by the U.S.Geographical Survey(USGS)and the University of Nebraska-Lincoln(UNL).The topography of the study area is taken from the Digital Elevation Map(DEM)of USGS.The model outputs include the components of the energy balance budget within the canopyand at the ground.the turbulence parameters in the atmospheric boundary layer and the wind.temperature and humidity profiles extending up to a height of 1500 m.In addition to the fine timeand spatial step,the unique feature of the present model is the incorporation of both dynamic andbiological effects of the Boreal forest into the model parameterization scheme.The model resultscompare favorably with BOREAS' IFC-1 data in 1994 when the forest was in the luxuriant growingperiod.