A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.     

Weathering process at the natural fission reactor of Bangombé

 The uranium deposits in the basin of Franceville (Gabon) host the only natural fission reactors known in the world. Unique geological conditions favoured a natural fission reaction 2 Ga ago. This was detected by anomalous isotopic compositions of uranium and rare earth elements (REE), which are produced by the fission reaction. In total, 16 reactor zones were found. Most of them are mined out. The reactor zone of Bangombé, is only 10–11 m below the surface. This site has been influenced by surface weathering processes. Six drill cores have been sampled at the site of the reactor zone of Bangombé during the course of the study and only one drill core (BAX 08) hit the core of the reactor. From these data and previous drilling campaigns, the reactor size is estimated to be 10 cm thick, 2–3 m wide and 4–6 m long. The migration of fission products can be traced by the anomalous isotope ratios of REE because of the fission process. The ¹⁴⁹Sm/¹⁴⁷Sm ratio close to the reactor zone is only 0.28 (normal: 0.92) because of the intense neutron capture of ¹⁴⁹Sm and subsequent transmutation, whereas ¹⁴⁷Sm is enriched by the fission reaction. Similar changes in isotopic patterns are detectable on other REE. The isotope ratios of Sm and Nd of whole rock and fracture samples surrounding the reactor indicate that fission-genic REE migrated only a few decimetres above and mainly below the reactor zone. Organic matter (bitumen) seems to act as a trap for fission-genic REE. Additional REE-patterns show less intense weathering with increasing depth in the log profile and support a simple weathering model. Received: 26 November 1999 · Accepted: 2 May 2000     

Effect of element size on the static finite element analysis of steep slopes

The accuracy of the computed stress distribution near the free surface of vertical slopes was evaluated in this study as a function of the element size, including aspect ratio. To accomplish this objective, a parametric study was carried out comparing stresses computed using the finite element method (FEM) to those obtained from a physical model composed of photoelastic material. The results of the study indicate a reasonable agreement between a gelatin model and the FEM model for shear stresses, and an overall good agreement between the two models for the principal stresses. For stresses along the top of the slope, the height of the element tends to be more important than width or aspect ratio, at least for aspect ratios up to 4. In all cases, the greatest difference between the two models occurs in the vicinity of the slope. Specifically, if H is defined as the slope height, an element height of H/10 appears to be adequate for the study of stresses deep within the slope, such as for typical embankment analyses. However, for cases where tensile stresses in the vicinity of the slope face which are critical, such as for the stability analysis of steep slopes, element heights as small as H/32, or higher‐order elements, are necessary. Copyright © 2001 John Wiley & Sons, Ltd.     

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.     

雪宝顶碱性花岗岩岩石地球化学与成矿控制

对雪宝顶碱性花岗岩具有高丰度的W、Sn、Be含量和W、Sn矿化花岗岩的微量元素组合特征进行了分析，结果显示：稀土元素具有∑REE含量与LREE／HREE比值偏低和负Eu异常的特点；元素比值K／Rb，U／Th，Nb／Ta的规律性变化以及矿脉产状特征等表明W-(Sn)-Be脉状矿床形成于岩浆演化晚期，与岩浆期后的热液活动有关。     

辽西中生代火山岩岩浆成因——稀土元素配分及同位素定量模式

本文通过对辽西及邻区中生代盆地火山岩稀土元素配分特征的研究及模拟运算，探讨了其岩浆成因及过程。结果显示南北岩区岩浆成因不同，南岩区岩源为MORB地幔和长期富集型幔熔融混合而成，北岩区为MORB和近期亏损型地幔熔融混合而成。     

安徽铜陵矿集区海底喷流沉积体系的流体包裹体微量元素对比

铜陵矿集区与铜-金(多金属)矿床有关的热液活动主要有两大体系。即与海西期海底喷流沉积有关的热液体系和与矽卡岩矿化有关的燕山期岩浆热液体系。查明这两类热液体系的流体包裹体特征对区域找矿和矿床成因研究都有实际意义。在包裹体岩相学研究基础上，应用ICP-MS技术和热爆提取方法，研究了新桥、冬瓜山、峙门口、铜官山、朝山等矿床具代表性的热液石英中流体包裹体的微量元素、稀土元素特征。结果表明，两类热液体系在流体包裹体特征上有较大的区别，在流体的微量元素和稀土元素特征方面也很不相同。海底喷流沉积体系的热液石英中流体包裹体与岩浆热液体系的相比。稀土总量较高，LREE／HREE比值较大，δEu不明显。且Mo／(W Sn)比值较高，反映流体中成矿物质的深源特征；Ga／T1、Rb／Cs大。Zr／Hf低，也不同于岩浆热液体系。     

火山岩与地形联合约束下海岛土壤稀土元素地球化学空间分布异质性——以广西涠洲岛为例

海岛土壤地球化学特征的研究对于认识海岛表生地质环境有重要意义。本研究以我国最大的第四纪火山海岛——涠洲岛为案例,重点研究亚热带海岛环境下火山岩风化土壤的稀土元素（REE）地球化学特征,并探讨影响其空间分布的关键因素。结果表明,涠洲岛土壤REE总量（ΣREE）自然变化范围较大,介于59.07~575.30 mg/kg,平均值为285.53 mg/kg,高于中国大陆土壤背景值（163.86 mg/kg）,指示火山岩风化海岛具有较高的土壤REE丰度。反距离权重插值法（IDW）分析指示,岛内土壤的REE含量高低变化呈现出以带状和斑块状为主体形态的空间分布不均一性,推断其成因与母质岩石类型和地形位置等因素密切相关。3种不同岩性的成土母岩的ΣREE含量变化序列为：玄武岩（236.66 mg/kg）>火山碎屑岩（215.58 mg/kg）>含生物碎屑海滩沉积物（18.99 mg/kg）,对应的土壤ΣREE含量具有类似变化规律,即玄武岩土壤（307.20 mg/kg）>火山碎屑岩土壤（297.46 mg/kg）>含生物碎屑海滩沉积物土壤（89.73 mg/kg）。该证据指示母岩是控制岛内土壤REE含量变化的首位因素,玄武岩整体上易形成具有较高REE丰度的土壤,火山碎屑岩次之,而含生物碎屑海滩沉积物反之。此外,本研究还发现火山碎屑岩单一岩性区内,发育的土壤REE含量也会随地形地貌差异而出现显著变化,其原因是地形地貌制约了土壤矿物和化学风化程度,海蚀丘陵比海积平原更利于海岛土壤中REE次生富集作用的发生。

     

鄂尔多斯盆地南部延长组富铀凝灰岩地球化学特征及形成机制

本文系统地研究了鄂尔多斯盆地南部延长组凝灰岩元素和同位素地球化学特征。结果表明：凝灰岩主量元素具有高K、高Al、Si中等至高、低Na的特征; REE总体特征为轻稀土富集、重稀土亏损,∑REE在4857～40212 μg/g 之间变化,Eu呈现出负异常,在034～082内变化,δCe范围在080～208之间; 微量元素中U含量在312～144 μg/g 之间,Ba、U、Th、Hf和Ce呈明显正异常,Nb、P、Ti和Rb呈现负异常。详细的分析表明：U的富集和一些微量元素如Ta、Dy、Lu等相关,Th与∑REE存在明显的正相关性,这些变化可能和盆地中铀矿的富集有关; 凝灰岩的硫同位素总体呈现深源特征。根据以上的研究对凝灰岩的形成原因及形成环境进行了探讨。鄂尔多斯盆地南部延长组凝灰岩来源于秦岭造山带火山作用,其中的铀元素富集是由于砂岩型铀矿中铀向凝灰岩迁移而被还原所致。本研究对盆地砂岩型铀矿的铀运移和富集机理认识具有一定指导意义。     

A New Protocol for the Astrometric Follow-up of Near Earth Asteroids

A new protocol was devised to improve the efficiency of astrometric follow-up observations of Near Earth Asteroids for the accurate determination of their orbits. It was implemented in the activities of the Spaceguard Central Node (SCN, a facility of the Spaceguard Foundation, established with the support of the European Space Agency) in the form of a Priority List. Here we describe this protocol and results obtained during five years of activity (2000–2004).