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.     

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

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.     

Cokriging nonstationary data

Universal cokriging is used to obtain predictions when dealing with multivariate random functions. An important type of nonstationarity is defined in terms of multivariate random functions with increments which are stationary of orderk. The covariance between increments of different variables is modeled by means of the pseudo-cross-covariance function. Criteria are formulated to which the parameters of pseudo-cross-covariance functions must comply so as to ensure positive-definiteness. Cokriging equations and the induced cokriging equations are given. The study is illustrated by an example from soil science.     

Discrimination of sedimentary facies by association analysis

The ability of association analysis to discriminate sedimentary facies was tested on Purdy's modal analyses of modern sediments of the Great Bahama Bank. Purdy's data set has served in the past as a standard reference for evaluating various multivariate classification algorithms. In order to adapt Purdy's data to association analysis, the percent abundance of the 12 constituents was converted to binary form by dichotomizing each variable on its mean value. The results obtained by association analysis are virtually identical to those obtained by Purdy and other authors. The same four main sedimentary facies were discriminated; 86% of the samples were identically classified (97% when misclassified borderline cases are counted as matches); the total partition variance of the classification is only negligibly greater (4%); and the grouping of the variables yielded the same four groups. The rank order of the three division-attributes responsible for the sample classification is fines, oolites, and corals. Association analysis has been employed by other authors to differentiate meaningful facies groups in studies of ancient reef carbonates, modern reef sediments, and heavy minerals in stream sediments. In all these studies, the results were found to be compatible with those obtained by using the continuous quantitative measurements, indicating that qualitative binary data may often be sufficient for the purpose of facies discrimination in many branches of geology and that association analysis is an effective method for this purpose.     

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

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

REVIEWS OF GEOGRAPHIC SOFTWARE

Reviews of geographic software in this article: DEMO-GRAPHICS: WORLD POPULATIONS AND PROJECTIONS. ESP GAUSS. CEMODEL S. Damus LIMDEP. William H. Greene MICROSTAT 4.1 OTIS PCIPS. (Personal Computer Image Processing System) . H.J. Meyers and R. Bernstein. REGRESSION ANALYSIS OF TIME SERIES (RATS) SPSS/PC+ URBAN DATA MANAGEMENT SOFTWARE (UDMS)     

Distribution and fluxes of metals in the St. Lawrence River from the outflow of Lake Ontario to Québec City

Six anchor stations in the St. Lawrence River from the outflow of Lake Ontario to Quebec City, were occupied for ca. 24 hours in June 1987 during low flow conditions. Samples of water and suspended particulate matter were separated by continuous-flow centrifugation, and were collected every two hours.During this sampling period, fluxes of dissolved forms of zinc, copper and nickel increased in Lac St. François and downstream relative to the sum of the fluxes for the two upstream stations at the outflow of Lake Ontario. Increases in the flux of dissolved zinc and copper were pronounced below Montreal and above Lac St. Pierre. For particulate forms of metals, all five metals show that there are significant inputs in the section of the St. Lawrence River between Lac St. François and the station just above the entrance to Lac St. Pierre.The average concentrations of dissolved cadmium, lead, zinc, copper and nickel ranged from 7–23 ng/l; 9–35 ng/l; 0.434–0.939 g/l; 0.15–0.89 g/l and 0.58–1.12 g/l respectively.Regression analysis of the dissolved and particulate metal concentrations suggests that the concentration of dissolved cadmium, lead, zinc and nickel can be predicted from the regression equation and the determination of particulate metal concentration. This prediction appears to be independent of the suspended particulate matter concentration which varies from ca. 1 to 10 mg/l from the outflow of Lake Ontario to Québec City.     

Sediment transport in analogue flume models compared with real‐world sedimentary systems: a new look at scaling evolution of sedimentary systems in a flume

Evolution of sedimentary systems at large temporal and spatial scales cannot be scaled down to laboratory dimensions by conventional hydraulic Froude scaling. Therefore, many researchers question the validity of experiments aiming to simulate this evolution. Yet, it has been shown that laboratory experiments yield stratigraphic responses to allocyclic forcing that are remarkably similar to those in real‐world prototypes, hinting at scale independency with strong dependence on boundary conditions but weak dependence on the actual sediment transport dynamics. This paper addresses the dilemma by contrasting sediment transport rules that apply in the laboratory with those that apply in real‐world geological systems. It is demonstrated that the generation of two‐dimensional stratigraphy in a flume can be simulated numerically by the non‐linear diffusion equation. Sediment transport theory is used to demonstrate that only suspension‐dominated meandering rivers should be simulated with linear diffusion. With increasing grain‐size (coarse sand to gravel) and shallowness of river systems, the prediction of long‐term transport must be simulated by non‐linear, slope‐dependent diffusion to allow for increasing transport rates and thus change in stratigraphic style. To point out these differences in stratigraphic style, three stages in infill of accommodation have been defined here: (i) a start‐up stage, when the system is prograding to base level (e.g. the shelf edge) with no sediment flux beyond the base‐level point; (ii) a fill‐up stage, when the system is further aggrading while progressively more sediment is bypassing base level with the progression of the infill; and (iii) a keep‐up stage, when more than 90% of the input is bypassing the base level and less than 10% is used for filling the accommodation. By plotting the rate of change in flux for various degrees of non‐linearity (varying the exponent in the diffusion equation) it was found that the error between model and real‐world prototype is largest for the suspension‐dominated prototypes, although never more than 30% and only at the beginning of the fill‐up stage. The error reduces to only 10% for the non‐linear sandy‐gravelly and gravelly systems. These results are very encouraging and open up ways to calibrate numerical models of sedimentary system evolution by such experiments.