Provenance and Magmatic-Metamorphic Evolution of a Variscan Island-Arc Complex: Constraints from U-Pb Dating, Petrology, and Geospeedometry of the Kyffhauser Crystalline Complex, Central Germany     

ZEH  ARMIN; GERDES  AXEL; WILL  THOMAS M.; MILLAR  IAN L. 《Journal of Petrology》2005,46(7):1393-1420

The Kyffhäuser Crystalline Complex, Central Germany, formspart of the Mid-German Crystalline Rise, which is assumed torepresent the Variscan collision zone between the East Avalonianterrane and the Armorican terrane assemblage. High-precisionU–Pb zircon and monazite dating indicates that sedimentaryrocks of the Kyffhäuser Crystalline Complex are youngerthan c. 470 Ma and were intruded by gabbros and diorites between345 ± 4 and 340 ± 1 Ma. These intrusions had magmatictemperatures between 850 and 900°C, and caused a contactmetamorphic overprint of the sediments at P–T conditionsof 690–750°C and 5–7 kbar, corresponding toan intrusion depth of 19–25 km. At 337 ± 1 Ma themagmatic–metamorphic suite was intruded by granites, syenitesand diorites at a shallow crustal level of some 7–11 km.This is inferred from a diorite, and conforms to P–T pathsobtained from the metasediments, indicating a nearly isothermaldecompression from 5–7 to 2–4 kbar at 690–750°C.Subsequently, the metamorphic–magmatic sequence underwentaccelerated cooling to below 400°C, as constrained by garnetgeospeedometry and a previously published K–Ar muscoviteage of 333 ± 7 Ma. With respect to P–T–D–tdata from surrounding units, rapid exhumation of the KCC canbe interpreted to result from NW-directed crustal shorteningduring the Viséan. KEY WORDS: contact metamorphism; U–Pb dating; hornblende; garnet; Mid-German Crystalline Rise; P–T pseudosection  相似文献   

Discrete element method modeling of inherently anisotropic rocks under uniaxial compression loading          下载免费PDF全文

K. Duan  C.Y. Kwok  M. Pierce 《国际地质力学数值与分析法杂志》2016,40(8):1150-1183

A new numerical approach is proposed in this study to model the mechanical behaviors of inherently anisotropic rocks in which the rock matrix is represented as bonded particle model, and the intrinsic anisotropy is imposed by replacing any parallel bonds dipping within a certain angle range with smooth‐joint contacts. A series of numerical models with β = 0°, 15°, 30°, 45°, 60°, 75°, and 90° are constructed and tested (β is defined as the angle between the normal of weak layers and the maximum principal stress direction). The effect of smooth‐joint parameters on the uniaxial compression strength and Young's modulus is investigated systematically. The simulation results reveal that the normal strength of smooth‐joint mainly affects the behaviors at high anisotropy angles (β > 45°), while the shear strength plays an important role at medium anisotropy angles (30°–75°). The normal stiffness controls the mechanical behaviors at low anisotropy angles. The angle range of parallel bonds being replaced plays an important role on defining the degree of anisotropy. Step‐by‐step procedures for the calibration of micro parameters are recommended. The numerical model is calibrated to reproduce the behaviors of different anisotropic rocks. Detailed analyses are conducted to investigate the brittle failure process by looking at stress‐strain behaviors, increment of micro cracks, initiation and propagation of fractures. Most of these responses agree well with previous experimental findings and can provide new insights into the micro mechanisms related to the anisotropic deformation and failure behaviors. The numerical approach is then applied to simulate the stress‐induced borehole breakouts in anisotropic rock formations at reduced scale. The effect of rock anisotropy and stress anisotropy can be captured. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Experimental characterization of the impact of temperature and humidity on the breakdown of soil water repellency in sandy soils and composts          下载免费PDF全文

Amii Whelan  Cedric Kechavarzi  Frederic Coulon  Stefan H. Doerr 《水文研究》2015,29(8):2065-2073

Soil water repellency is a widespread phenomenon with the capacity to alter hydrological and geomorphological processes. Water repellency decays with time, and the consequences are only of concern during the timescale at which the water repellency persists. This study aimed to characterize the influence of temperature and humidity on the breakdown of water repellency. Apparent contact angle measurements were carried out on samples consisting of sand treated with stearic acid as well as naturally repellent dune sands and composts. Temperature and humidity were controlled using a cooled incubator and a purpose designed enclosed box in which humidity could be raised or lowered. Results showed the contact angle of the stearic‐acid‐treated sands decayed with time and that there was a significant increase with stearic acid concentration. For all samples, the decay in apparent contact angle could be described with a continuous breakdown model. The stearic‐acid‐treated sands showed a significant increase in contact angle with relative humidity at a temperature of 10 and 20 °C. These differences diminished with increasing temperature. Similar results were seen for the dune sands and composts. Despite the influence of temperature and humidity on contact angles, there was no significant change in the rate at which the contact angle decayed in any sample. Absolute humidity was found to provide a more relevant indicator than relative humidity when assessing the influence of humidity on repellency over a range of temperatures. The contact angle initially increased with absolute humidity before plateauing owing to the confounding effect of temperature. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Elastic behavior of weakly cemented contact          下载免费PDF全文

V. Langlois 《国际地质力学数值与分析法杂志》2015,39(8):854-860

We investigate the elastic behavior of weakly cemented contact. We show that the radial distribution of stresses and the stiffnesses of a cemented contact are governed by the ratio a/RΛ, where R, a, and Λ are, respectively, the grain radius, the contact size, and the ratio of the elastic moduli of cement and grains. Moreover, we show that a cemented contact is always less stiff than a Hertzian contact having a similar size. Finally, we propose accurate approximate expressions of the contact stiffnesses. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Low‐P and high‐T metamorphism of basalts: Insights from the Sudbury impact melt sheet aureole and thermodynamic modelling     

Taus R. C. Jrgensen  Douglas K. Tinkham  C. Michael Lesher 《Journal of Metamorphic Geology》2019,37(2):271-313

Low‐pressure and high‐temperature (LP–HT) metamorphism of basaltic rocks, which occurs globally and throughout geological time, is rarely constrained by forward phase equilibrium modelling, yet such calculations provide valuable supplementary thermometric information and constraints on anatexis that are not possible to obtain from conventional thermometry. Metabasalts along the southern margin of the Sudbury Igneous Complex (SIC) record evidence of high‐grade contact metamorphism involving partial melting and melt segregation. Peak metamorphic temperatures reached at least ~925°C at ~1–3 kbar near the SIC contact. Preservation of the peak mineral assemblage indicates that most of the generated melt escaped from these rocks leaving a residuum characterized by a plagioclase–orthopyroxene–clinopyroxene–ilmenite‐magnetite±melt assemblage. Peak temperatures reached ~875°C up to 500 m from the SIC lower contact, which marks the transition to metabasalts that only experienced incipient partial melting without melt loss. Metabasalts ~500 to 750 m from the SIC contact are characterized by a similar two‐pyroxene mineral assemblage, but typically contain abundant hornblende that overgrew clino‐ and orthopyroxene along an isobaric cooling path. Metabasalts ~750 to 1,000 m from the SIC contact are characterized by a hornblende–plagioclase–quartz–ilmenite assemblage indicating temperatures up to ~680°C. Mass balance and phase equilibria calculations indicate that anatexis resulted in 10–20% melt generation in the inner ~500 m of the aureole, with even higher degrees of melting towards the contact. Comparison of multiple models, experiments, and natural samples indicates that modelling in the Na₂O–CaO–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O₂ (NCFMASHTO) system results in the most reliable predictions for the temperature of the solidus. Incorporation of K₂O in the most recent amphibole solution model now successfully predicts dehydration melting by the coexistence of high‐Ca amphibole and silicate melt at relatively low pressures (~1.5 kbar). However, inclusion of K₂O as a system component results in prediction of the solidus at too low a temperature. Although there are discrepancies between modelling predictions and experimental results, this study demonstrates that the pseudosection approach to mafic rocks is an invaluable tool to constrain metamorphic processes at LP–HT conditions.  相似文献   

Tensionless–frictionless interaction of flexible annular foundation with a transversely isotropic multi‐layered half‐space          下载免费PDF全文

Morteza Eskandari‐Ghadi  Ghasem Gorji‐Bandpey  Azizollah Ardeshir‐Behrestaghi  Seyed Masoud Nabizadeh 《国际地质力学数值与分析法杂志》2015,39(2):155-174

A transversely isotropic multi‐layered half‐space, with axis of material symmetry perpendicular to the free surface, supports a flexible either annular or solid circle foundation. The contact area of the foundation and the half‐space is considered to be both frictionless and tensionless. The foundation is assumed to be affected by a vertical static axisymmetric load. Detailed analysis of the interaction of these two systems with different thickness of layers is the target of this paper. With the use of ring load Green's functions for both the foundation and the continuum half‐space, an integral equation accompanied with some inequalities is introduced to model the complex BVP. With the incorporation of ring‐shape FEM, we are capable of capturing both regular and singular solution smoothly. The validity of the combination of the analytical and numerical method is proved with comparing the results of this paper with a number of benchmark cases of both linear and nonlinear interaction of circular and annular foundation with half‐space. Some new illustrations are presented to portray the aspect of the anisotropy and layering of the half‐space. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

Numerical and theoretical research on stress distribution in the loosening zone of the trapdoor problem     

Luju Liang  Changjie Xu 《国际地质力学数值与分析法杂志》2019,43(7):1426-1447

The traditional theory of soil arching effect was developed on the assumption that stress distribution in the loosening zone is uniform. However, because of the deflection of principal stress' direction, the stress distribution in the loosening zone is actually ununiform. For the evaluation of principal stress axis deflection and stress redistribution, a discrete element method numerical model of trapdoor problem is established for the simulation of soil arching effect. Based on the numerical results, an arc shape of major principal stress trajectory and uniform horizontal stress distribution at the same depth of the loosening zone are adopted. An analytical model is raised to estimate the average loosening earth pressure acting on the trapdoor and stress distribution in the loosening zone at a limit state. In addition, comparison studies are carried out between the predictions of the proposed solutions and discrete element method numerical results as well as available model test results, thereby validating the accuracy of the proposed theoretical model. Both numerical and theoretical results indicate that the vertical stress distribution in the loosening zone is obviously ununiform. The load acting in the middle of loosening zone is transferred toward two sides so that the vertical stress distribution in loosening zone is concave.  相似文献   

Parallelized combined finite-discrete element (FDEM) procedure using multi-GPU with CUDA     

Quansheng Liu  Weiqin Wang  Hao Ma 《国际地质力学数值与分析法杂志》2020,44(2):208-238

This paper focuses on the efficiency of finite discrete element method (FDEM) algorithmic procedures in massive computers and analyzes the time-consuming part of contact detection and interaction computations in the numerical solution. A detailed operable GPU parallel procedure was designed for the element node force calculation, contact detection, and contact interaction with thread allocation and data access based on the CUDA computing. The emphasis is on the parallel optimization of time-consuming contact detection based on load balance and GPU architecture. A CUDA FDEM parallel program was developed with the overall speedup ratio over 53 times after the fracture from the efficiency and fidelity performance test of models of in situ stress, UCS, and BD simulations in Intel i7-7700K CPU and the NVIDIA TITAN Z GPU. The CUDA FDEM parallel computing improves the computational efficiency significantly compared with the CPU-based ones with the same reliability, providing conditions for achieving larger-scale simulations of fracture.  相似文献   

基于部分最小二乘岭估计的粗差定值定位     

金丽宏  吕玉婷  汪耀  许运鹏  张小龙 《测绘工程》2016,25(6):61-64

在利用部分最小二乘原理进行粗差定值定位时,模型的法方程矩阵可能存在病态性,使得到的粗差定值定位结果不可靠。文中针对观测数据包含多个粗差且法方程病态问题,利用岭估计处理病态问题,建立部分最小二乘岭估计的粗差定值定位方法,给出粗差搜索步骤,利用迭代算法实现多个粗差的定值和定位。通过模拟算例分析部分最小二乘法、部分最小二乘岭估计在粗差搜索方面的效果,从另一个角度探讨粗差处理方法,推广现有的误差理论,证明文中方法的有效性。  相似文献   

重力场数学拟合模型解的影响因素探究     

鹿荻  陈锴 《测绘科学》2016,41(3):28-33

为了提高地球重力场模型不适定方程求解的精度,该文采用谱分析方法从级数展开阶数、数据采样率及数据缺失量3个方面探索影响数学拟合效果的根本因素:从常用的三角级数及勒让德级数模型出发,引出重力场拟合模型球谐函数模型,观察在改变级数展开阶数、数据采样率及数据缺失量等情况下所对应设计矩阵谱结构的变化,并从微观上研究影响误差分配的有关因素及最小奇异值对误差的决定性作用,为探求重力场模型解不准的原因及实现更高精度的全球重力场模型的建立提供参考。  相似文献