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.     

渤、黄海生态环境模拟的参数敏感度空间差异分析

随着海洋生态系统模型的发展，生态变量增多，众多生物过程参数量值的确定成为制约生态环境模拟的瓶颈问题，生态系统结构区域性要求模型中的生态参数具有区域差异。为探究不同海区的关键参数及参数敏感度的空间差异，本研究在渤、黄海建立了ROMS-CoSiNE物理–生物耦合的高分辨率生态系统模型，并对13种生态参数的敏感度空间分布进行分析。结果表明:南黄海中部与渤海及近岸海域的敏感度差异较大。渤海敏感度最大的参数为决定光合速率的浮游植物P-I曲线初始斜率，其次为浮游动物捕食半饱和常数和浮游动物最大捕食率。而南黄海中部敏感度最大的参数为浮游动物最大捕食率，其次为浮游植物死亡率和浮游植物P-I曲线初始斜率。结合敏感度分布及浮游植物生物量收支得出，渤海水体透明度较南黄海偏低、浮游植物生长光限制较强，是引起浮游植物P-I曲线初始斜率敏感度在渤海高于黄海的主要原因。浮游动物最大捕食率及浮游植物死亡率的敏感度空间差异，受渤、黄海浮游植物生物量差异的影响，与生态系统中的高度非线性特征有关。     

An evaluation of Guided Regularized Random Forest for classification and regression tasks in remote sensing

New Earth observation missions and technologies are delivering large amounts of data. Processing this data requires developing and evaluating novel dimensionality reduction approaches to identify the most informative features for classification and regression tasks. Here we present an exhaustive evaluation of Guided Regularized Random Forest (GRRF), a feature selection method based on Random Forest. GRRF does not require fixing a priori the number of features to be selected or setting a threshold of the feature importance. Moreover, the use of regularization ensures that features selected by GRRF are non-redundant and representative. Our experiments based on various kinds of remote sensing images, show that GRRF selected features provides similar results to those obtained when using all the available features. However, the comparison between GRRF and standard random forest features shows substantial differences: in classification, the mean overall accuracy increases by almost 6% and, in regression, the decrease in RMSE almost reaches 2%. These results demonstrate the potential of GRRF for remote sensing image classification and regression. Especially in the context of increasingly large geodatabases that challenge the application of traditional methods.     

杭州湾中部实测波浪特性分析

利用杭州湾中部一年实测波浪资料,分析该地区的波参数统计变化特性,采用最小二乘法拟合得出相关参数之间的关系式,并统计分析了频谱特性。结果表明:杭州湾中部以小浪和轻浪为主,在各月分布较为均匀,年平均周期为2.97 s;最大一次波浪过程为冷空气影响所致,影响强度大于台风莫拉克;常浪向分布在东北至东南向,其中又以小浪出现较多,对应的波周期以2~4 s为主;强浪向主要分布在西北和东北方向,对应的波周期主要分布在3~5 s;多数特征波参数之间相关关系较好;由谱分析得知风浪占多数,且以单峰为主。研究结果可为相关工程和理论研究提供基础资料。     

Isotope hydrological study of mean transit time in the granitic Strengbach catchment (Vosges massif,France): application of the FlowPC model with modified input function

Measurements of ¹⁸O concentrations in precipitation, soil solution, spring and runoff are used to determine water transit time in the small granitic Strengbach catchment (0·8 km²; 883–1146 m above sea level) located in the Vosges Mountains of northeastern France. Water transit times were calculated by applying the exponential, exponential piston and dispersion models of the FlowPC program to isotopic input (rainfall) and output (spring and stream water) data sets during the period 1989–95. The input function of the model was modified compared with the former version of the model and estimated by a deterministic approach based on a simplified hydrological balance. The fit between observed and calculated output data showed marked improvements compared with results obtained using the initial version of the model. An exponential piston version of the model applied to spring water indicates a 38·5 month mean transit time, which suggests that the volume in the aquifer, expressed in water depth, is 2·4 m. A considerable thickness (>45 m) of fractured bedrock may be involved for such a volume of water to be stored in the aquifer. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

利用Lageos 1 SLR资料解算1990-2001年的地球定向参数

利用 12年的Lageos 1卫星激光测距资料 (1990 - 2 0 0 1)解算得到了地球定向参数 (EOP) ,将该序列的结果与同期的EOP(IERS)C0 4进行比较 ,其外符精度为 :极移XP— 0 .4 0mas,YP— 0 .4 2mas ,日长变化Dr— 0 .0 35ms。     

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

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

APPLICATION OF GEOGRAPHICAL PARAMETER DATABASE TO ESTABLISHMENT OF UNIT POPULATION DATABASE

Now GIS is turning into a good tool in handling geographical, economical, and population data, so we can obtain more and more information from these data. On the other hand, in some cases, for a calamity, such as hurricane, earthquake, flood, drought etc., or a decision-making, such as setting up a broadcasting transmitter, building a chemical plant etc., we have to evaluate the total population in the region influenced by a calamity or a project. In this paper, a method is put forward to evaluate the population in such special region. Through exploring the correlation of geographical parameters and the distribution of people in the same region by means of quantitative analysis and qualitative analysis, unit population database (1km× 1km) is established. In this way, estimating the number of people in a special region is capable by adding up the population in every grid involved in this region boundary. The geographical parameters are obtained from topographic database and DEM database on the scale of     

Influence of Structural Non-Stationarity of Surface Roughness on Morphological Characterization and Mechanical Deformation of Rock Joints

Summary Structural non-stationarity of surface roughness affects accurate morphological characterization as well as mechanical behaviour of rock joints at the laboratory scale using samples with a size below the stationarity threshold. In this paper, the effect of structural non-stationarity of surface roughness is investigated by studying the scale dependence of surface roughness and mechanical behaviour of rock joints. The results show that the structural non-stationarity mainly affects the accurate characterization of the surface roughness of the fracture samples. It also controls the amount and location of the contact areas during shear tests, which in turn affects the mechanical properties and asperity degradation of the samples. It is concluded that for accurate determination of the morphological and mechanical properties of rock joints at laboratory and field scales, samples with size equal to or larger than the stationarity threshold are required. Author’s address: Nader Fardin, Rock Mechanics Group, Department of Mining Engineering, Faculty of Engineering, University of Tehran, P.O. Box: 11365/4563, Tehran, Iran