A Generalized Backward Euler algorithm for the numerical integration of an isotropic hardening elastoplastic model for mechanical and chemical degradation of bonded geomaterials

An extended version of the classical Generalized Backward Euler (GBE) algorithm is proposed for the numerical integration of a three‐invariant isotropic‐hardening elastoplastic model for cemented soils or weak rocks undergoing mechanical and non‐mechanical degradation processes. The restriction to isotropy allows to formulate the return mapping algorithm in the space of principal elastic strains. In this way, an efficient and robust integration scheme is developed which can be applied to relatively complex yield surface and plastic potential functions. Moreover, the proposed algorithm can be linearized in closed form, thus allowing for quadratic convergence in the global Newton iteration. A series of numerical experiments are performed to illustrate the accuracy and convergence properties of the algorithm. Selected results from a finite element analysis of a circular footing on a soft rock layer undergoing chemical weathering are then presented to illustrate the algorithm performance at the boundary value problem level. Copyright © 2002 John Wiley & Sons, Ltd.     

A relook at NEH‐4 curve number data and antecedent moisture condition criteria

This paper investigates the variation of the popular curve number (CN) values given in the National Engineering Hand Book–Section 4 (NEH‐4) of the Soil Conservation Service (SCS) with antecedent moisture condition (AMC) and soil type. Using the volumetric concept, involving soil, water, and air, a significant condensation of the NEH‐4 tables is achieved. This leads to a procedure for determination of CN for gauged as well as ungauged watersheds. The rainfall‐runoff events derived from daily data of four Indian watersheds exhibited a power relation between the potential maximum retention or CN and the 5‐day antecedent rainfall amount. Including this power relation, the SCS‐CN method was modified. This modification also eliminates the problem of sudden jumps from one AMC level to the other. The runoff values predicted using the modified method and the existing method utilizing the NEH‐4 AMC criteria yielded similar results. Copyright © 2006 John Wiley & Sons, Ltd.     

A bimodal four‐parameter lognormal linear model of soil water repellency persistence

Soil water repellency may be characterized in terms of the delayed infiltration time of a water droplet resting on the soil surface, which is, water drop penetration time (WDPT), or repellency persistence. Such repellency persistence varies nonlinearly with soil water content (θ_g), although no models have been proposed to reproduce the variation of WDPT with θ_g in soils. Dynamic factor analysis (DFA) is used to identify two common patterns of unexplained variability in a scattered dataset of WDPT versus θ_g measurements. A four‐parameter lognormal distribution was fitted to both common patterns obtained by DFA, and these were combined additively in a weighted multiple linear bimodal model. We show how such an empirical model is capable of reproducing a large variety of WDPT versus θ_g curve shapes (N = 80) both within a wide range of measured WDPTs (0–17 000 s) and for samples with organic matter content ranging from 21·7 to 80·6 g (100 g)^?1. Copyright © 2009 John Wiley & Sons, Ltd.     

Experimental study on the dynamic modulus of compacted loess under bidirectional dynamic load

The dynamic characteristics of compacted loess are of great significance to the seismic construction of the Loess Plateau area in Northwest China, where earthquakes frequently occur. To study the change in the dynamic modulus of the foundation soil under the combined action of vertical and horizontal earthquakes, a hollow cylindrical torsion shear instrument capable of vibrating in four directions was used to perform two-way coupling of compression and torsion of Xi'an compacted loess under different dry density and deviator stress ratios. The results show that increasing the dry density can improve the initial dynamic compression modulus and initial dynamic shear modulus of compacted loess. With an increase in the deviator stress ratio, the initial dynamic compression modulus increases, to a certain extent, but the initial dynamic shear modulus decreases slightly. The dynamic modulus gradually decreases with the development of dynamic strain and tends to be stable, and the dynamic modulus that reaches the same strain increases with an increasing dry density. At the initial stage of dynamic loading, the attenuation of the dynamic shear modulus with the strain development is faster than that of the dynamic compression modulus. Compared with previous research results, it is determined that the dynamic modulus of loess under bidirectional dynamic loading is lower and the attenuation rate is faster than that under single-direction dynamic loading. The deviator stress ratio has a more obvious effect on the dynamic compression modulus. The increase in the deviator stress ratio can increase the dynamic compression modulus, to a certain extent. However, the deviator stress ratio has almost no effect on the dynamic shear modulus, and can therefore be ignored.     

Inverse Analysis of Deep Excavation Using Differential Evolution Algorithm

This paper presents the applications of the differential evolution (DE) algorithm in back analysis of soil parameters for deep excavation problems. A computer code, named Python‐based DE, is developed and incorporated into the commercial finite element software ABAQUS, with a parallel computing technique to run an FE analysis for all trail vectors of one generation in DE in multiple cores of a cluster, which dramatically reduces the computational time. A synthetic case and a well‐instrumented real case, that is, the Taipei National Enterprise Center (TNEC) project, are used to demonstrate the capability of the proposed back‐analysis procedure. Results show that multiple soil parameters are well identified by back analysis using a DE optimization algorithm for highly nonlinear problems. For the synthetic excavation case, the back‐analyzed parameters are basically identical to the input parameters that are used to generate synthetic response of wall deflection. For the TNEC case with a total of nine parameters to be back analyzed, the relative errors of wall deflection for the last three stages are 2.2, 1.1, and 1.0%, respectively. Robustness of the back‐estimated parameters is further illustrated by a forward prediction. The wall deflection in the subsequent stages can be satisfactorily predicted using the back‐analyzed soil parameters at early stages. Copyright © 2014 John Wiley & Sons, Ltd.     

Zircon U-Pb Geochronology and Geochemical Characteristics of the Volcanic Host Rocks from the Tongyu VHMS Copper Deposit in the Western North Qinling Orogen and Their Geological Significance

Precise in situ zircon U-Pb dating and Lu–Hf isotopic measurement using an LA-ICP-MS system, whole-rock major and trace element geochemistry and Sr–Nd isotope geochemistry were conducted on the volcanic host rocks of the Tongyu copper deposit on the basis of further understanding of its geological characteristics. Three zircon samples from the volcanic host rocks yielded 206Pb/238 U weighted average ages ranging from 436±4 Ma to 440±5 Ma, which are statistically indistinguishable and coeval with the ca. 440 Ma northward subduction event of the Paleo-Qinling oceanic slab. The volcanic host rocks were products of magmatic differentiation that evolved from basalt to andesite to dacite to rhyolite, forming an integrated tholeiitic island arc volcanic rock suite. The primitive mantle-normalized trace element patterns for most samples show characteristics of island arc volcanic rocks, such as relative enrichment of LILE(e.g. Th, U, Pb and La) and depletion of HFSE(e.g. Nb, Ta, Ti, Zr and Hf). Discrimination diagrams of Ta/Yb vs Th/Yb, Ta vs Th, Yb vs Th/Ta, Ta/Hf vs Th/Hf, Hf/3 vs Th vs Nb/16, La vs La/Nb and Nb vs Nb/Th all suggest that both the volcanic host rocks from the Tongyu copper deposit and the volcanic rocks from the regional Xieyuguan Group were formed in an island arc environment related to subduction of an oceanic slab. Values of ISr(0.703457 to 0.708218) and εNd(t)(-2 to 5.8) indicate that the source materials of volcanic rocks from the Tongyu copper deposit and the Xieyuguan Group originated from the metasomatised mantle wedge with possible crustal material assimilation. Most of the volcanic rock samples show good agreement with the values of typical island arc volcanic rocks in the ISr-εNd(t) diagram. The involvement of crustal-derived material in the magma of the volcanic rocks from the Tongyu copper deposit was also reflected in the zircon εHf(t) values, which range from-3.08 to 10.7, and the existence of inherited ancient xenocrystic zircon cores(2616±39 Ma and 1297±22 Ma). The mineralization of the Tongyu copper deposit shows syn-volcanic characteristics such as layered orebodies interbedded with the volcanic rock strata, thus, the zircon U-Pb age of the volcanic host rocks can approximately represent the mineralization age of the Tongyu copper deposit. Both the Meigou pluton and the volcanic host rocks were formed during the ca. 440 Ma northward subduction of the Paleo-Qinling Ocean when high oxygen fugacity aqueous hydrothermal fluid released by dehydration of the slab and the overlying sediments fluxed into the mantle wedge, triggered partial melting of the mantle wedge, and activated and extracted Cu and other ore-forming elements. The magma and ore-bearing fluid upwelled and erupted, and consequently formed the island arc volcanic rock suite and the Tongyu VHMS-type copper deposit.     

高度计测距精度对沿轨迹重力异常反演的影响

应用求解沿轨迹重力异常的垂线偏差法以及求解空间分辨率的交叉谱分析法,建立了高度计测距精度与沿轨迹重力异常反演精度以及空间分辨率的关联性模型。首先依据卫星测高原理,给出了沿轨迹重力异常的误差传播公式,然后以此为基础通过推导交叉谱分析中一致性系数与信噪比的数学表达式,建立了高度计测距精度与空间分辨率的解析关系。数值仿真结果表明:雷达高度计测距精度与沿轨迹重力异常反演精度成正比关系,与空间分辨率成幂函数关系,即高度计测距精度提高m倍,沿轨迹重力异常反演精度提高m倍,全球海域平均空间分辨率提高m0.464 4倍。将数值仿真结果与相关文献中对实际测高数据的处理结果进行比较,验证了理论分析及模型的正确性。     

加速迭代硬阈值ISAR稀疏成像

加速迭代硬阈值(AIHT)作为一种压缩感知(CS)算法,具有计算复杂度低、内存占用小的特点,在较小数目的观测下,经过有限次迭代便能获得较好的重建效果。因此,本文将该算法应用到逆合成孔径雷达(ISAR)稀疏成像,对该方法的具体特征和适用性进行了详细分析。仿真结果表明,在相同信噪比和回波脉冲数条件下,基于AIHT的ISAR成像方法,相比于稀疏贝叶斯算法(SBL),在重建时间上缩短了至少2个量级;相比于迭代加权最小二乘法(IRLS)和平滑l0范数算法,在重建质量上提高了约3倍。经实测数据处理结果验证,在1/4原始脉冲数条件下,该方法仍具有良好的聚焦性能和压低旁瓣效果。