考虑Hansbo渗流的二维Biot固结有限元分析

为进一步深入研究弹性饱和黏性土地基的二维固结机制,引入Hansbo渗流方程描述固结过程中的非达西渗流,修正Biot二维固结方程。基于加权残数法,给出相应的有限元数值求解格式。通过和饱和黏性土一维非达西渗流固结理论有限体积法数值结果的对比,证明数值计算方法的有效性。在此基础上,探讨Hansbo渗流参数对二维地基固结进程的影响。计算结果表明,在固结初期,Hansbo渗流将增强Mandel-Cryer效应,增大孔压的峰值,并延长孔隙水压力达到峰值的时间;在固结中后期,整个土层存在孔隙水压力滞后现象。同时,Hansbo渗流将阻碍地基沉降的发展。而且,上述影响会随着Hansbo渗流参数的增大而更加明显。     

赣南北部黄陂河流域离子型稀土矿地区水质与健康风险评价

水质安全和健康是保障中国赣南老区乡村振兴发展的重要因素。赣南离子型稀土矿长期开发利用,导致浸矿剂和矿体中重金属元素等危害人体健康的物质进入水循环系统,给周边乡村饮用水卫生安全带来了潜在风险。目前,针对当地复垦后稀土矿及周边地区水质和健康风险缺乏系统调查评价,本文以《生活饮用水卫生标准》(GB 5749—2006)为评价依据,选择赣南北部黄陂河流域典型离子型稀土矿及周边的水体开展调查研究,采用电感耦合等离子体发射光谱/质谱等技术测定锰、镉等元素含量,采用水质指数(WQI)、危害商(HQ)、致癌风险(CR)评价了锰和铅等9种指标及其健康风险。结果表明：地表水中的异常指标有氨氮(平均值750μg/L)、锰(平均值207μg/L),地下水中的异常指标有氨氮(平均值4533μg/L)、锰(平均值4009μg/L);世界卫生组织(WHO)公布的Ⅰ类致癌物砷在地表水及地下水均未见异常。WQI显示研究区内85.7%的地表水适宜饮用。地表水及地下水中氨氮的HQ平均值<1,对人类健康没有不良影响;地下水中锰的HQ平均值>1,可能会对人类健康产生不良影响。地表水及地下水中致癌元素砷的CR值分布区...     

Experimental determination of REE partition coefficients between zircon,garnet and melt: a key to understanding high‐T crustal processes

The partitioning of rare earth elements (REE) between zircon, garnet and silicate melt was determined using synthetic compositions designed to represent partial melts formed in the lower crust during anatexis. The experiments, performed using internally heated gas pressure vessels at 7 kbar and 900–1000 °C, represent equilibrium partitioning of the middle to heavy REE between zircon and garnet during high‐grade metamorphism in the mid to lower crust. The D_REE (zircon/garnet) values show a clear partitioning signature close to unity from Gd to Lu. Because the light REE have low concentrations in both minerals, values are calculated from strain modelling of the middle to heavy REE experimental data; these results show that zircon is favoured over garnet by up to two orders of magnitude. The resulting general concave‐up shape to the partitioning pattern across the REE reflects the preferential incorporation of middle REE into garnet, with D_Gd (zircon/garnet) ranging from 0.7 to 1.1, D_Ho (zircon/garnet) from 0.4 to 0.7 and D_Lu (zircon/garnet) from 0.6 to 1.3. There is no significant temperature dependence in the zircon–garnet REE partitioning at 7 kbar and 900–1000 °C, suggesting that these values can be applied to the interpretation of zircon–garnet equilibrium and timing relationships in the ultrahigh‐T metamorphism of low‐Ca pelitic and aluminous granulites.     

A high strain‐rate constitutive model for sand and its application in finite‐element analysis of tunnels subjected to blast

The paper presents a constitutive model for simulating the high strain‐rate behavior of sands. Based on the concepts of critical‐state soil mechanics, the bounding surface plasticity theory and the overstress theory of viscoplasticity, the constitutive model simulates the high strain‐rate behavior of sands under uniaxial, triaxial and multi‐axial loading conditions. The model parameters are determined for Ottawa and Fontainebleau sands, and the performance of the model under extreme transient loading conditions is demonstrated through simulations of split Hopkinson pressure bar tests up to a strain rate of 2000/s. The constitutive model is implemented in a finite‐element analysis software Abaqus to analyze underground tunnels in sandy soil subjected to internal blast loads. Parametric studies are conducted to examine the effect of relative density and type of sand and of the depth of tunnel on the variation of stresses and deformations in the soil adjacent to the tunnels. Copyright © 2012 John Wiley & Sons, Ltd.     

Geochemical Aspects of Hypogene Hydrothermal Alteration Zones in the Chol‐qeshlaghi Area,NW Iran: Constrains on REEs

Chol-qeshlaghi altered area lies in the northwestern part of the post-collisional Urumieh-Dokhtar magmatic arc, NW Iran. Pervasive silicic, argillic, phyllic and propylitic altered zones appears to be intimately affiliated to the fluids derivative of upper Oligocene Khankandi granodiorite. This paper is dedicated to the identification of geochemical characteristics of hydrothermal alterations, focusing on the determination of the mass gains and losses of REEs, to gain significant insights regard...     

粤北油洞岩体SHRIMP锆石U-Pb年龄、地球化学特征及其成因研究

油洞岩体位于诸广南部岩体中部,是一个重要的产铀岩体,岩性为中粒小斑状二云母花岗岩。SHRIMP锆石U-Pb年龄为232±4Ma(MSWD=3.2),属于印支早期岩浆活动产物。该岩体在主量元素方面,具有富硅(SiO2平均为72.65%)、富铝(A/CNK值平均为1.12)和高的K2O/Na2O比值(平均为1.79);在微量元素方面,大离子元素富集,Ba、Sr、P、Ti、Nb、Ta亏损明显,具有高的Rb/Sr(平均为8.08)和Rb/Nb比值(平均为20.96);在稀土元素方面,轻稀土明显富集,配分模式呈右倾型,Eu亏损明显;在同位素方面,εNd(t)值低(平均为-11.9),(87Sr/86Sr)i高(平均为0.72330),Nd模式年龄古老(平均为1954 Ma)。这些特征一致表明,油洞岩体属于典型的壳源型花岗岩范畴,是在华南地块和印支地块碰撞结束后不久形成的伸展构造环境中,位于中—下地壳部位的古—中元古代地壳组分由于在地壳缩短之后的伸展、减薄环境下产生的减压、导水和地幔上涌等因素的综合影响下,由泥质岩和砂质岩混合组成的源区发生部分熔融而形成。     

基于三维裂隙网络的裂隙岩体表征单元体研究

岩体表征单元体（REV）的研究一直是岩土工程中的重要问题。裂隙岩体REV是其内部复杂结构在尺寸效应上的体现,归因于岩体内裂隙系统的随机性。在三维网络模型的基础上通过计算体积节理数建立了确定岩体REV的方法。这种方法考虑了岩体内裂隙的三维空间密度的特征,采用体积节理数这一参数的收敛状况确定岩体REV的大小。建立了实际岩体的模拟裂隙系统,用一定大小的立方体来分割裂隙系统,统计不同立方体内的体积节理数,并将均值作为最终裂隙岩体的参数。将此参数作为裂隙岩体的属性,通过判别其收敛状况获得了岩体REV的尺寸。考虑了岩体不连续面的特征,建立了REV的分析模型,得出工程岩体的REV的尺寸是裂隙平均迹长的4倍的结论。     

面向对象的土石坝参数随机反演程序设计

采用面向对象的软件PowerBuilder作为开发平台 ,运用动态链接库的原理 ,建立了一套面向对象的反演分析法 ,定义了对象类之间的关系 ,确立了有限元的消息流 ,使程序的界面更加友好 ,程序操作更为直观。并针对岩土工程反分析中使用的传统贝叶斯法中存在的缺陷 ,提出了扩展贝叶斯法 ,从概率论和数理统计的原理出发 ,建立了基于决策信息论中AIC准则和最大熵准则的岩土工程随机反演的准则函数。所建立的计算模型和分析方法 ,通过了数值算例的验证 ,从而也表明本文提出的模型和分析方法基础可靠、仿真性强 ,具有广阔的应用前景     

Pore orientation of granular materials during biaxial compression

The local pore spaces in granular materials tend to be aligned parallel to the major principal stress direction upon particle mobilization. Manifestation of this response has been numerically validated in our previous studies with the aid of discrete element method modeling and image processing techniques during creep and shearing. We now extend the modeling of pore geometry, constructed with spherical particles, to assemblies of particle clumps. Two-dimensional simulations are performed for both loose and dense assemblies of spherical particles and particle clumps. Each particle packing is bound by rigid or flexible walls and subjected to biaxial compression and the particle mobilization effect on the evolution of pore orientation is explored. Randomly shaped pores surrounded by adjacent particles are geometrically quantified by Delaunay tessellation and fitted with ellipses. Results show that localization is apparent in dense assemblies, in particular for clumped particle packing, while loose assemblies exhibit diffusive failure. Small pores within well-defined shear bands tend to align either parallel to the direction of the shear band or perpendicular to the major principal stress. On the other hand, small pores within the blocks and large pores have a tendency to become elongate towards the major principal stress direction. This study reveals for the first time that pore orientation is dependent upon particle shape, pore size, and assembly conditions on the pore and global scales.