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.     

Spatial patch structure and adaptive strategy for desert shrub of Reaumuria soongorica in arid ecosystem of the Heihe River Basin

In many arid ecosystems, vegetation frequently occurs in high-cover patches interspersed in a matrix of low plant cover. However, theoretical explanations for shrub patch pattern dynamics along climate gradients remain unclear on a large scale. This context aimed to assess the variance of the Reaumuria soongorica patch structure along the precipitation gradient and the factors that affect patch structure formation in the middle and lower Heihe River Basin (HRB). Field investigations on vegetation patterns and heterogeneity in soil properties were conducted during 2014 and 2015. The results showed that patch height, size and plant-to-patch distance were smaller in high precipitation habitats than in low precipitation sites. Climate, soil and vegetation explained 82.5% of the variance in patch structure. Spatially, R. soongorica shifted from a clumped to a random pattern on the landscape towards the MAP gradient, and heterogeneity in the surface soil properties (the ratio of biological soil crust (BSC) to bare gravels (BG)) determined the R. soongorica population distribution pattern in the middle and lower HRB. A conceptual model, which integrated water availability and plant facilitation and competition effects, was revealed that R. soongorica changed from a flexible water use strategy in high precipitation regions to a consistent water use strategy in low precipitation areas. Our study provides a comprehensive quantification of the variance in shrub patch structure along a precipitation gradient and may improve our understanding of vegetation pattern dynamics in the Gobi Desert under future climate change.

     

Crystal-Melt Separation and the Development of Isotopic Heterogeneities in Hybrid Magmas

If a magma is a hybrid of two (or more) isotopically distinctend-members, at least one of which is partially crystalline,separation of melt and crystals after hybridization will leadto the development of isotopic heterogeneities in the magmaas long as some of the pre-existing crystalline material (antecrysts)retains any of its original isotopic composition. This holdstrue whether the hybridization event is magma mixing as traditionallyconstrued, bulk assimilation, or melt assimilation. Once a magma-scaleisotopic heterogeneity is formed by crystal–melt separation,it is essentially permanent, persisting regardless of subsequentcrystallization, mixing, or equilibration events. The magnitudeof the isotopic variability resulting from crystal–meltseparation can be as large as that resulting from differentialcontamination, multiple isotopically distinct sources, or insitu isotopic evolution. In one model, a redistribution of one-thirdof the antecryst cargo yielded a crystal-enriched sample with⁸⁷Sr/⁸⁶Sr of 0·7058, whereas the complementary crystal-poorsample has ⁸⁷Sr/⁸⁶Sr of 0·7068. In other models, crystal-richsamples are enriched in radiogenic Sr. Isotopic heterogeneitiescan be either continuous (controlled by the modal distributionof crystals and melt) or discontinuous (when there is completeseparation of crystals and liquid). The first case may be exemplifiedby some isotopically zoned large-volume rhyolites, formed bythe eruptive inversion of a modally zoned magma chamber. Inthe latter case, the isotopic composition of any (for example)interstitial liquid will be distinct from the isotopic compositionof the bulk crystal fraction. The separation of such an interstitialliquid may explain the presence of isotopically distinct late-stageaplites in plutons. Crystal–melt separation provides anadditional option for the interpretation of isotopically zonedor heterogeneous magmas. This option is particularly attractivefor systems whose chemical variation is otherwise explicableby fractionation-dominated processes. Non-isotopic chemicalheterogeneities can also develop in this fashion. KEY WORDS: isotopic heterogeneity; zoning; hybrid magma; crystal separation; Sr isotopes; aplite; rhyolite     

冈底斯新近纪钾质火山作用:消减沉积物折返的地球化学与Sr-Nd-Pb同位素证据

在冈底斯岩浆岩带的羊八井地区,林子宗火山岩系上部出现小规模的粗面质火山岩.岩石学与地球化学研究表明粗面岩与时空密切相关的大体积安山质-英安质-流纹质火山岩属于不同的火山岩系:碱性钾质系列与钙碱性系列.详细的地球化学研究证明林子宗晚期小体积钾质岩具有独立的岩浆源区,而非早期钙碱性系统低压岩浆过程演化的产物.粗面质熔岩SiO2含量为62.91%～64.63%,具有高K2O(7.52%～8.05%)、Al2O3(16.64%-17.35%)、低TiO2(0.59%～0.68%)与MsO(0.15%～0.77%).钾质熔岩富集LILE与LREE,亏损HFSE,具有高Rb/Sr(1.1～2.3)、Th/La(0.59～0.80)、Th/U(6.3～10.9)和低Ce/Pb(4.1～6.2)与Ba/Th(4.3～14.5)比值,其同位素组成变化范围分别是87Sr/86Sr(I)=0.7068～0.7075、143Nd/144Nd(I)=0.51241～0.51252、206Pb/204Pb=18.87～18.95、207Pb/204Pb=15.63～15.70和208Pb/204Pb=39.24～39.68.粗面岩的地球化学与同位素特征表明冈底斯新近纪钾质火山岩来源于角闪岩相的消减沉积物熔融.因此,羊八井新近纪钾质火山熔岩提供了印度-亚洲大陆碰撞早期消减沉积物折返的证据.     

青藏高原北缘阿尔金走滑边界的侧向扩展——甘肃北山晚新生代走滑构造与地壳稳定性分析

对北山地区遥感影像和野外地质特征的分析表明,自阿尔金断裂带向NW方向依次出露三危山-双塔断裂、大泉断裂和红柳河断裂。这些断裂近于平行,且同为左行走滑断裂,具有相似的展布特征,空间走向均为NE40-50°,断裂系末端均发育“树枝状”分支断层．在断层夹块之间形成“多米诺”构造,构成了北山地区主要的构造样式。断层谷地沉积物分析和断层泥ESR年代学测试结果表明,三危山-双塔断裂形成于上新世（N2k）,大泉断裂形成于早更新世（1．2-1．5Ma）．而北山地区分支断层和次级断层的活动在400ka之后。对北山地区断裂构造几何学和年代学的研究表明．阿尔金断裂系晚新生代以来向NW方向的侧向扩展．是阿尔金走滑边界重要的生长方式。北山地区特殊的走滑构造组合样式．使该地区的构造变形难于在某条断层上聚集能量,而分散在若干条次级断层上的位移量又微乎其微,该地区成为“最稳定的活动区”。     

软粘土中插入式大直径薄壁圆筒结构的一种极限状态分析模型

根据插入式大直径薄壁圆筒结构的工作特点建立了其合理的破坏机制，以此为基础建立了受横向荷载作用的插入式大直径薄壁圆筒结构的上限法极限分析模型，从中可以得到作用于插入式大直径薄壁圆筒结构上的横向极限荷载近似的上限解；并根据上限法模型计算结果，提出了计算受横向荷载作用的插入式大圆筒结构侧壁极限抗力的简化修正模型，利用此修正模型，基于极限平衡法，可以求解作用于插入式大直径薄壁圆筒结构上的横向极限荷载。     

Behavior of Pile Groups under Lateral Load

Based on investigation and model tests, and in combination with the research work on group effect for pile groups under lateral loads relating to the code of fixed offshore platforms, a series of studies have been performed on the behavior and failure mechanism of laterally loaded pile groups, critical pile spacing inducing group effect, lateral bearing capacity of pile groups and its main influence factors, the stress-strain relationship for single piles and pile groups and so on. Some new laws about non-uniformity of load distribution in the longitudinal direction of pile groups and load-deflection (p - y) curves for pile groups have been discovered, and an empirical formula is presented in order to remedy the defect of current calculating methods at home and abroad. These results can be used for reference in the design of pile foundation under lateral loads.