Micromechanical viscoelasto‐plastic models and finite element implementation for rate‐independent and rate‐dependent permanent deformation of stone‐based materials

This paper presents parallel and serial viscoelasto‐plastic models to simulate the rate‐independent and the rate‐dependent permanent deformation of stone‐based materials, respectively. The generalized Maxwell viscoelastic and Chaboche's plastic models were employed to formulate the proposed parallel and serial viscoelasto‐plastic constitutive laws. The finite element (FE) implementation of the parallel model used a displacement‐based incremental formulation for the viscoelastic part and an elastic predictor—plastic corrector scheme for the elastoplastic component. The FE framework of the serial viscoelasto‐plastic model employed a viscoelastic predictor—plastic corrector algorithm. The stone‐based materials are consisted of irregular aggregates, matrix and air voids. This study used asphalt mixtures as an example. A digital sample was generated with imaging analysis from an optically scanned surface image of an asphalt mixture specimen. The modeling scheme employed continuum elements to mesh the effective matrix, and rigid bodies for aggregates. The ABAQUS user material subroutines defined with the proposed viscoelasto‐plastic matrix models were employed. The micromechanical FE simulations were conducted on the digital mixture sample with the viscoelasto‐plastic matrix models. The simulation results showed that the serial viscoelasto‐plastic matrix model generated more permanent deformation than the parallel one by using the identical material parameters and displacement loadings. The effect of loading rates on the material viscoelastic and viscoelasto‐plastic mixture behaviors was investigated. Permanent deformations under cyclic loadings were determined with FE simulations. The comparison studies showed that the simulation results correctly predicted the rate‐independent and rate‐dependent viscoelasto‐plastic constitutive properties of the proposed matrix models. Overall, these studies indicated that the developed micromechanical FE models have the abilities to predict the global viscoelasto‐plastic behaviors of the stone‐based materials. Copyright © 2009 John Wiley & Sons, Ltd.     

Description of deformation process in inherently anisotropic granular materials

The main focus in this work is on modeling of mechanical response of granular materials that display inherent anisotropy. Both the experimental and numerical investigations are described. First, the results of direct shear as well as drained/undrained triaxial tests that involve crushed limestone with elongated angular‐shaped particles are reviewed. Afterward, a mathematical framework is presented for modeling of elastic/ inelastic deformation that incorporates the multi‐laminate approach. The deformation is monitored on a set of randomly oriented planes, and the formulation incorporates the thickness of the shear band that is associated with sliding/separation process. A systematic procedure for identification of material functions/ parameters is outlined that is based on the results of direct shear tests, and the framework is later applied to simulate the behavior under triaxial conditions. The results of numerical simulations are compared with the experimental data. Copyright © 2011 John Wiley & Sons, Ltd.     

Garnet crystal plasticity in the continental crust,new example from south Madagascar

Garnet (10 vol.%; pyrope contents 34–44 mol.%) hosted in quartzofeldspathic rocks within a large vertical shear zone of south Madagascar shows a strong grain‐size reduction (from a few cm to ～300 μm). Electron back‐scattered diffraction, transmission electron microscopy and scanning electron microscope imaging coupled with quantitative analysis of digitized images (PolyLX software) have been used in order to understand the deformation mechanisms associated with this grain‐size evolution. The garnet grain‐size reduction trend has been summarized in a typological evolution (from Type I to Type IV). Type I, the original porphyroblasts, form cm‐sized elongated grains that crystallized upon multiple nucleation and coalescence following biotite breakdown: biotite + sillimanite + quartz = garnet + alkali feldspar + rutile + melt. These large garnet grains contain quartz ribbons and sillimanite inclusions. Type I garnet is sheared along preferential planes (sillimanite layers, quartz ribbons and/or suitably oriented garnet crystallographic planes) producing highly elongated Type II garnet grains marked by a single crystallographic orientation. Further deformation leads to the development of a crystallographic misorientation, subgrains and new grains resulting in Type III garnet. Associated grain‐size reduction occurs via subgrain rotation recrystallization accompanied by fast diffusion‐assisted dislocation glide. This plastic deformation of garnet is associated with efficient recovery as shown by the very low dislocation densities (10¹⁰ m^?3 or lower). The rounded Type III garnet experiences rigid body rotation in fine‐grained matrix. In the highly deformed samples, the deformation mechanisms in garnet are grain‐size‐ and shape‐dependent: dislocation creep is dominant for the few large grains left (>1 mm; Type II garnet), rigid body rotation is typical for the smaller rounded grains (300 μm or less; Type III garnet) whereas diffusion creep may affect more elliptic garnet (Type IV garnet). The P–T conditions of garnet plasticity in the continental crust (≥950 °C; 11 kbar) have been identified using two‐feldspar thermometry and GASP conventional barometry. The garnet microstructural and deformation mechanisms evolution, coupled with grain‐size decrease in a fine‐grained steady‐state microstructure of quartz, alkali feldspar and plagioclase, suggests a separate mechanical evolution of garnet with respect to felsic minerals within the shear zone.     

Preconditioned IDR(s) iterative solver for non‐symmetric linear system associated with FEM analysis of shallow foundation

Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix K_ep (ILU(0)‐K_ep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐K_ep. To further reduce computational cost, the global stiffness matrix K_ep is divided into two parts. The first part is the linear elastic stiffness matrix K_e, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐K_e preconditioner is more time effective than IDR(s = 6) with ILU(0)‐K_ep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd.     

Explicit integration of bounding surface model for the analysis of earthquake soil liquefaction

This paper presents a new plasticity model developed for the simulation of monotonic and cyclic loading of non‐cohesive soils and its implementation to the commercial finite‐difference code FLAC, using its User‐Defined‐Model (UDM) capability. The new model incorporates the framework of Critical State Soil Mechanics, while it relies upon bounding surface plasticity with a vanished elastic region to simulate the non‐linear soil response. Stress integration of constitutive relations is performed using a recently proposed explicit scheme with automatic error control and substepping, which so far has been employed in the literature only for constitutive models aiming at monotonic loading. The overall accuracy of this scheme is evaluated at element level by simulating cyclic loading along complex stress paths and by using iso‐error maps for paths involving change of the Lode angle. The performance of the new constitutive model and its stress integration scheme in complex boundary value problems involving earthquake‐induced liquefaction is evaluated, in terms of accuracy and computational cost, via a number of parametric analyses inspired by the successful simulation of the VELACS centrifuge Model Test No. 2 studying the lateral spreading response of a liquefied sand layer. Copyright © 2009 John Wiley & Sons, Ltd.     

岩体质量分级的风险分析方法

提出了基于RMR岩体质量分级系统的岩体质量研究风险分析方法。该方法分析步骤如下 :(1)通过岩芯样品的现场观测和实验室试验获得分类所需的变量 ;(2 )统计分析拟合得出各变量的分布函数及参数 ;(3)运用Monte Carlo模拟方法获得 2万个RMR值 ,并将结果绘成岩体质量描述图 ;(4)利用以上结果作出岩体质量风险分析评价。该方法用于润扬长江公路大桥岩体质量评价研究获得良好效果     

利用水汽图象分析中尺度云团

利用ＮＯＡＡ卫星提供的水汽图（６．７μｍ），配合ＧＭＳ云图、多普勒雷达照片，分析了中尺度暴雨云团发生的位置、演变及其结构，最后给出了有关方法的预报框图。     

震级特征图特点与a,b值关系

本文主要研究了a、b值的变化对震级特征图的影响。当a、b值的变化幅度不太大时,仍然可从震级特征图上识别出最小均匀震级等震级阈值;当a、b值的变化达到相当的幅度后,震级阈值点变模糊了。     

用储层岩石抽提物的饱和烃色谱指纹识别油气层

利用储层岩石抽提物的饱和烃色 质总离子图的指纹特征判识储层的流体（油、气）类型。油层或油饱和的储层,饱和烃色 质总离子流图的特征是正构烷烃碳数分布宽,碳数分布在C₁₅～C₃₈之间,与原油样品的正构烷烃碳数分布相似;凝析气层的正构烷烃碳数分布略窄一些,碳数分布在C₁₅～C₃₅之间,低碳数（小于C₂₁）的正构烷烃相对于油层富集,高碳数（大于C₂₁）的正构烷烃丰度明显低于油层;干气层的正构烷烃碳数分布最窄,碳数分布在C₁₅～ C₂₈之间,只有低碳数的正构烷烃,高碳数的正构烷烃丰度极低。据此可用来识别油层、凝析气层和干气层,也可用于一些测井资料不全或测井质量差的老井的油气层复查、测井难以识别的火成岩油气层的识别,还可用于地层评价,为油气酸化层位的优选提供科学依据。     

“地气图”方法用于省级旱涝预测的实践

简要介绍用“地气图”方法进行降水预测的基本思路,总结此方法的预报实践,指出省级尺度旱涝预测需要解决的几个关键问题,强调除加密观测外,上游地热涡增强的能量频散作用是下游地热涡发展的重要原因