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.     

山东省区域地质构造演化探讨

山东省区域地质构造演化分为5个阶段.①陆核形成阶段形成太古宙高级区,地壳分异成稳定的花岗岩穹窿和活动的绿岩带,第一次克拉通化完成.②陆块发生形成阶段地壳向刚性发展,在华北陆核硅铝壳的基础上先后有3次张开、闭合裂谷作用,第二次克拉通化完成.这一阶段演化在鲁西地区主要表现为挤压作用,形成大量造山花岗岩;鲁东地区则以拉张作用为主,形成海槽,产生沉积.③秦昆洋形成演化阶段四堡期沿鲁东南部地壳拉张,在华北板块与扬子板块间形成秦昆洋.晋宁期秦昆洋关闭,华北板块与扬子板块对接碰撞,沿胶南造山带产生大量同碰撞花岗岩,同时产生超高压变质作用及形成丰富多彩的碰撞构造.晋宁运动最终形成统一的原始中国古陆,第三次克拉通化完成.④陆块发展阶段鲁西地区地壳频繁升降,形成广泛的海相及海陆交互相沉积;鲁东地区则以造山抬升为主,地层沉积较少.⑤滨太平洋发展阶段该阶段的主要特征是断块构造发育,形成盆岭构造格局,产生大陆边缘花岗岩,构造体系由古亚洲构造域转向滨太平洋构造域.     

渭河盆地北缘断裂带活动特征的初步研究

本文从渭河盆地北缘断裂的形成和活动时代,活动特征等资料出发,结合北缘断裂带及整个盆地历史地震活动和新生界地层发育特征的综合分析,对北缘断裂带的活动期次,主要断层的运动幅度和滑动速率及其时空演变规律和机制等问题进行了探讨。文章指出,北缘断裂带的形成是一个由盆地中心向北逐渐扩展的过程,自上新世起,断层活动明显有东强西弱的变化特征,而且扩展方向也发生了偏转。这一转变及活动强度的东西差异与山西剪切带对渭河盆地的影响密切相关     

U–Pb SHRIMP zircon geochronology and T–t–d history of the Kampa Dome, southern Tibet

Structural, petrographic and geochronologic studies of the Kampa Dome provide insights into the tectonothermal evolution of orogenic crust exposed in the North Himalayan gneiss domes of southern Tibet. U–Pb ion microprobe dating of zircons from granite gneiss exposed at the deepest levels within the dome yields concordia ²⁰⁶Pb/²³⁸U age populations of 506 ± 3 Ma and 527 ± 6 Ma, with no evidence of new zircon growth during Himalayan orogenesis. However, the granite contains penetrative deformation fabrics that are also preserved in the overlying Paleozoic strata, implying that the Kampa granite is a Cambrian pluton that was strongly deformed and metamorphosed during Himalayan orogenesis. Zircons from deformed leucogranite sills that cross-cut Paleozoic metasedimentary rocks yield concordant Cambrian ages from oscillatory zoned cores and discordant ages ranging from ca. 491–32 Ma in metamict grains. Since these leucogranites clearly post-date the metasedimentary rocks they intrude, the zircons are interpreted as xenocrysts that are probably derived from the Kampa granite. The Kampa Dome formed via a series of progressive orogenic events including regional ~ N–S contraction and related crustal thickening (D₁), predominately top-to-N ductile shearing and crustal extension (D₂), top-to-N brittle–ductile faulting and related folding on the north limb of the dome, localized top-to-S faulting on the southern limb of the dome, and crustal doming (D₃), and continued N–S contraction, E–W extension and doming (D₄). Structural and geochronologic variability amongst adjacent North Himalayan gneiss domes may reflect changes in the magnitude of crustal exhumation along the North Himalayan antiform, possibly relating to differences in the mid-crustal geometry of the exhuming fault systems.     

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.粗面岩的地球化学与同位素特征表明冈底斯新近纪钾质火山岩来源于角闪岩相的消减沉积物熔融.因此,羊八井新近纪钾质火山熔岩提供了印度-亚洲大陆碰撞早期消减沉积物折返的证据.     

空间刚架结构极限强度分析

本文用理想结构单元法对空间刚架结构进行了极限强度分析，理想结构单元法在于将结构分成最大可能的单元。以简单的形式归并几何非线性与材料非线性。按增量形式加载，直到达到极限强度。本文根据Ｒａｓｈｅｄ理论导出了对称的切向方程组，节省了存贮单元、改善了破坏准则的描述。算例结果与非线性有限单元法结果十分吻合。     

On the structural evolution of the Central Trough in the Norwegian and Danish sectors of the North Sea

The Central Trough of the North Sea is not a simple rift graben. It is an elongated area of regional subsidence which was initiated in mid Cretaceous times and continued to subside through to the late Tertiary. Its form is not representative of pre-mid Cretaceous tectonics.In Late Permian times the North Sea was divided into a northern and southern Zechstein basin by the E-W trending Mid North Sea-Ringkøbing-Fyn High. The latter was dissected by a narrow graben trending NNW through the Tail End Graben and the Søgne Basin. The Feda Graben was a minor basin on the northern flank of the Mid North Sea High at this time. This structural configuration persisted until end Middle Jurassic times when a new WNW trend separated the Tail End Graben from the Søgne Basin. Right lateral wrench movement on this new trend caused excessive subsudence in the Tail End and Feda Grabens while the Søgne Basin became inactive.Upper Jurassic subsidence trends continued during the Early Cretaceous causing the deposition of large thicknesses of sediments in local areas along the trend. From mid Cretaceous times the regional subsidence of the Central Trough was dominant but significant structural inversions occurred in those areas of maximum Early Cretaceous and Late Jurassic subsidence.