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.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

岩体三维节点偶对分析及其应用

本文讨论了节点偶对分析的三维扩展。节点偶对分析是针对摩擦滑动节理单元的约束边界进行的序列矢量判定。在三维节理单元的分析中,几何约束和相应的力学判据极为复杂,采用节点偶对等效结合力方法才能迎刃而解。文中对三维转向节理模型的分析得到合理的结果,它表明本文所提出的方法可成功地应用于节理岩体的有限元分析。     

The geochemistry of gem opals as evidence of their origin

Seventy-seven gem opals from ten countries were analyzed by inductively coupled plasma-mass spectrometry (ICP-MS) through a dilution process, in order to establish the nature of the impurities. The results are correlated to the mode of formation and physical properties and are instrumental in establishing the geographical origin of a gem opal. The geochemistry of an opal is shown to be dependant mostly on the host rock, at least for examples from Mexico and Brazil, even if modified by weathering processes. In order of decreasing concentration, the main impurities present are Al, Ca, Fe, K, Na, and Mg (more than 500 ppm). Other noticeable elements in lesser amounts are Ba, followed by Zr, Sr, Rb, U, and Pb. For the first time, geochemistry helps to discriminate some varieties of opals. The Ba content, as well as the chondrite-normalized REE pattern, are the keys to separating sedimentary opals (Ba > 110 ppm, Eu and Ce anomalies) from volcanic opals (Ba < 110 ppm, no Eu or Ce anomaly). The Ca content, and to a lesser extent that of Mg, Al, K and Nb, helps to distinguish gem opals from different volcanic environments. The limited range of concentrations for all elements in precious (play-of-color) compared to common opals, indicates that this variety must have very specific, or more restricted, conditions of formation. We tentatively interpreted the presence of impurities in terms of crystallochemistry, even if opal is a poorly crystallized or amorphous material. The main replacement is the substitution of Si⁴⁺ by Al³⁺ and Fe³⁺. The induced charge imbalance is compensated chiefly by Ca²⁺, Mg²⁺, Mn²⁺, Ba²⁺, K⁺, and Na⁺. In terms of origin of color, greater concentrations of iron induce darker colors (from yellow to “chocolate brown”). This element inhibits luminescence for concentrations above 1000 ppm, whereas already a low content in U (≤ 1 ppm) induces a green luminescence.     

等离子体质谱法测定玄武岩中微量元素三种样品预处理方法的比较

在测定玄武岩的微量元素时,为了保证得到的微量元素结果反映源区特征,需经预处理去除玄武岩的孔洞中充填的碳酸盐等物质。选取雷州半岛新生代玄武岩三个样品,每一个样品用三种预处理方法进行预处理,等离子体质谱法测定其中的微量元素。分析结果表明,用稀硝酸浸泡样品,会导致玄武岩中大多数微量元素严重丢失,得到的数据失真;而用稀盐酸浸泡去除碳酸盐的方法基本不改变玄武岩中微量元素的含量。     

内蒙古敖包吐萤石矿床的Sr、Nd、Pb同位素地球化学特征

敖包吐萤石矿床是内蒙古北部苏莫查干地区单一萤石矿集区中的一个代表性矿床,产于早二叠世大石寨组火山-沉积岩与早白垩世敖包吐花岗岩的接触带上。文章通过分析该矿床岩、矿石的微量元素和稀土元素,揭示出萤石的成矿作用可分为2个阶段,即交代作用和充填作用。交代作用过程中大石寨组的结晶灰岩可能为萤石的形成提供了部分Ca来源,萤石矿石的稀土元素配分模式与海水基本类似,具有Ce负异常;成矿作用后期主要表现为充填作用,形成颗粒粗大的萤石,表现为重稀土元素富集的特征,并随着萤石的沉淀析出,稀土元素总量逐渐下降,反映出成矿流体经历了较长期的演化过程。各地层单元、花岗岩体和萤石矿石的Sr、Nd、Pb同位素研究表明,萤石的放射性同位素组成具有壳、幔源混合的特点,成矿物质来源具有多源性。早白垩世敖包吐花岗岩可能是萤石中F的主要来源,而大石寨组的结晶灰岩则可能提供了Ca。另外,Pb、Nd同位素的极大不均一性,有可能是成矿流体在运移过程中对艾力格庙群放射性组分的选择性吸收的结果。萤石成矿作用与钾玄岩的时空关系暗示了萤石的成矿过程可能是中国东部岩石圈减薄和下地壳的置换地质事件的结果。在构造转型的过程中,燕山中期富碱的酸性花岗岩浆的活动分异出富含F的成矿流体,与幔源流体混合,沿区域重新活化的深大断裂和大石寨组的层间破碎带上升,交代其间的灰岩透镜体,从而形成敖包吐中型萤石矿床。     

Timing of UHP metamorphism in the Hong’an area,western Dabie Mountains,China: evidence from zircon U–Pb age,trace element and Hf isotope composition

The Hong’an area (western Dabie Mountains) is the westernmost terrane in the Qinling-Dabie-Sulu orogen that preserves UHP eclogites. The ages of the UHP metamorphism have not been well constrained, and thus hinder our understanding of the tectonic evolution of this area. LA-ICPMS U–Pb age, trace element and Hf isotope compositions of zircons of a granitic gneiss and an eclogite from the Xinxian UHP unit in the Hong’an area were analyzed to constrain the age of the UHP metamorphism. Most zircons are unzoned or show sector zoning. They have low trace element concentrations, without significant negative Eu anomalies. These metamorphic zircons can be further subdivided into two groups according to their U–Pb ages, and trace element and Lu–Hf isotope compositions. One group with an average age of 239 ± 2 Ma show relatively high and variable HREE contents (527 ≥ Lu_N ≥ 14) and ¹⁷⁶Lu/¹⁷⁷Hf ratios (0.00008–0.000931), indicating their growth prior to a great deal of garnet growth in the late stage of continental subduction. The other group yields an average age of 227 ± 2 Ma, and shows consistent low HREE contents and ¹⁷⁶Lu/¹⁷⁷Hf ratios, suggesting their growth with concurrent garnet crystallization and/or recrystallization. These two groups of age are taken as recording the time of prograde HP to UHP and retrograde UHP–HP stages, respectively. A few cores have high Th/U ratios, high trace element contents, and a clear negative Eu anomaly. These features support a magmatic origin of these zircon cores. The upper intercept ages of 771 ± 86 and 752 ± 70 Ma for the granitic gneiss and eclogite, respectively, indicate that their protoliths probably formed as a bimodal suite in rifting zones in the northern margin of the Yangtze Block. Young Hf model ages (T _DM1) of magmatic cores indicate juvenile (mantle-derived) materials were involved in their protolith formation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

An experimental study of the grain-scale processes of peridotite melting: implications for major and trace element distribution during equilibrium and disequilibrium melting

The grain-scale processes of peridotite melting were examined at 1,340°C and 1.5 GPa using reaction couples formed by juxtaposing pre-synthesized clinopyroxenite against pre-synthesized orthopyroxenite or harzburgite in graphite and platinum-lined molybdenum capsules. Reaction between the clinopyroxene and orthopyroxene-rich aggregates produces a melt-enriched, orthopyroxene-free, olivine + clinopyroxene reactive boundary layer. Major and trace element abundance in clinopyroxene vary systematically across the reactive boundary layer with compositional trends similar to the published clinopyroxene core-to-rim compositional variations in the bulk lherzolite partial melting studies conducted at similar P–T conditions. The growth of the reactive boundary layer takes place at the expense of the orthopyroxenite or harzburgite and is consistent with grain-scale processes that involve dissolution, precipitation, reprecipitation, and diffusive exchange between the interstitial melt and surrounding crystals. An important consequence of dissolution–reprecipitation during crystal-melt interaction is the dramatic decrease in diffusive reequilibration time between coexisting minerals and melt. This effect is especially important for high charged, slow diffusing cations during peridotite melting and melt-rock reaction. Apparent clinopyroxene-melt partition coefficients for REE, Sr, Y, Ti, and Zr, measured from reprecipitated clinopyroxene and coexisting melt in the reactive boundary layer, approach their equilibrium values reported in the literature. Disequilibrium melting models based on volume diffusion in solid limited mechanism are likely to significantly underestimate the rates at which major and trace elements in residual minerals reequilibrate with their surrounding melt. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.