Chemo‐mechanics of cemented granular solids subjected to precipitation and dissolution of mineral species

This paper studies the chemo‐mechanics of cemented granular solids in the context of continuum thermodynamics for fluid‐saturated porous media. For this purpose, an existing constitutive model formulated in the frame of the Breakage Mechanics theory is augmented to cope with reactive processes. Chemical state variables accounting for the reactions between the solid constituents and the solutes in the pore fluid are introduced to enrich the interactions among the microstructural units simulated by the model (i.e., grains and cement bonds). Two different reactive processes are studied (i.e., grain dissolution and cement precipitation), using the chemical variables to describe the progression of the reactions and track changes in the size of grains and bonds. Finally, a homogenization strategy is used to derive the energy potentials of the solid mixture, adopting probability density functions that depend on both mechanical and chemical indices. It is shown that the connection between the statistics of the micro‐scale attributes and the continuum properties of the solid enables the mathematical capture of numerous mechanical effects of lithification and chemical deterioration, such as changes in stiffness, expansion/contraction of the elastic domain, and development of inelastic strains during reaction. In particular, the model offers an interpretation of the plastic strains generated by aggressive environments, which are here interpreted as an outcome of chemically driven debonding and comminution. As a result, the model explains widely observed macroscopic signatures of geomaterial degradation by reconciling the energetics of the deformation/reaction processes with the evolving geometry of the microstructural attributes. Copyright © 2015 John Wiley & Sons, Ltd.     

Anisotropic damage of rock joints during cyclic loading: constitutive framework and numerical integration

This work describes a constitutive framework for modeling the behavior of rough joints under cyclic loading. Particular attention is paid to the intrinsic links between dilatancy, surface degradation, and mobilized shear strength. The framework also accounts for the important effect of shear‐induced anisotropy. The resulting approach is fully three‐dimensional and is not restricted to plane‐displacement kinematics. Both the governing formulation and an algorithm for implicit numerical integration are presented. While the proposed methods are general, we also postulate a specific model that is compared with experimental data. It employs relatively few free parameters but shows good agreement with laboratory tests. Copyright © 2013 John Wiley & Sons, Ltd.     

Enhancement of backtracking search algorithm for identifying soil parameters

In this paper, an enhanced backtracking search algorithm (so-called MBSA-LS) for parameter identification is proposed with two modifications: (a) modifying the mutation of original backtracking search algorithm (BSA) considering the contribution of current best individual for accelerating convergence speed and (b) novelly incorporating an efficient differential evolution (DE) as local search for improving the quality of population. The proposed MBSA-LS is first validated with better performance than the original BSA and some other typical state-of-the-art optimization algorithms on a benchmark of soil parameter identification in terms of effectiveness, efficiency, and robustness. Then, the efficiency of the MBSA-LS is further illustrated by two representative cases: identifying soil parameters from both laboratory tests and field measurements. All comparisons demonstrate that the proposed MBSA-LS algorithm can give accurate results in a short time. Finally, to conveniently solve the problems of parameter identification, a practical tool ErosOpt for parameter identification is developed by integrating the proposed MBSA-LS and some other efficient algorithms for readers to conduct the parameter identification using optimisation algorithms.     

An original semi‐discrete approach to assess gas conductivity of concrete structures

For civil engineering structures with a tightness role, structural permeability is a key issue. In this context, this paper presents a new proposition of a numerical modelling of leakage rate through a cracked concrete structure undergoing mode I cracking. The mechanical state of the material, considered in the framework of continuum mechanics based on finite element modelling, is described by means of the stress‐based nonlocal damage model which takes into account the stress state and provides realistic local mechanical fields. A semi‐discrete method based on the strong discontinuity approach to estimate crack opening is then considered in the post‐treatment phase. Using a Poiseuille's like relation, the coupling between the mechanical state of the material and its dry gas conductivity is performed. For validation purposes, an original experimental campaign is conducted on a dry concrete disc loaded in a splitting setup. During the loading, gas conductivity and digital image correlation analysis are performed. The comparison with the 3D experimental mechanical global response highlights the performance of the mechanical model. The comparison between crack openings measured by digital image correlation and estimated by the strong discontinuity method shows a good agreement. Finally, the results of the semi‐discrete approach coupled with the gas conductivity compared with experimental data show a good estimation of the structural conductivity. Consequently, if the mechanical problem is well modelled at the global scale, then the proposed approach provides good estimation of gas conductivity. Copyright © 2016 John Wiley & Sons, Ltd.     

Review and enhancement of 3D concrete models for large‐scale numerical simulations of concrete structures

The present paper focuses on selected plasticity and damage‐plasticity models for describing the 3D material behavior of concrete. In particular, a plasticity model and a damage‐plasticity model are reviewed and evaluated. Based on the results of the evaluation, enhancements are proposed, aiming at improving the correspondence between predicted and observed material behavior and aiming at implementing a robust and efficient stress update algorithm in a finite element program for performing large‐scale 3D numerical simulations of concrete structures. The capabilities of the concrete models are demonstrated by 3D numerical simulations of benchmark tests with combined bending and torsional loading and combined compression and shear loading and by a large‐scale 3D finite element analysis of a model test of a concrete arch dam. Copyright © 2011 John Wiley & Sons, Ltd.     

雅安芦山4.20震后县级单元生态系统风险评价

区域生态风险评价是对各种生态风险及环境问题进行评价和管理的重要手段。针对雅安地震灾区特殊的自然地理及生态环境特征,选取芦山县为研究对象,采用遥感、GIS及SPSS统计分析的方法,通过风险源、风险受体、暴露和易损性分析,建立生态风险综合评价模型,划分生态风险区类型,进而提出生态风险管理对策。结果表明:1)微度和低度生态风险区集中分布在高海拔的森林及草地生态系统,该区生物多样性丰富,抗干扰能力较强,地质灾害及人类活动影响较小;2)中度和高度生态风险区具有沿农田及建设用地生态系统集聚分布的特征,该区地质灾害频繁,地壳活动性较强,生态系统抵抗灾害的能力较差。研究结果可为地震灾区防御、规避风险及安全选址提供科学依据。     

X射线衍射技术应用于宝石鉴定-合成及晶体结构研究进展

因宝石检测常附带无损要求、宝石评价具科学性和不确定性等特点,X射线衍射(XRD)在宝石学这门矿物学的分支学科中的应用具有特殊性,其应用领域主要集中在三个方面。①宝石鉴定领域:对单晶宝石和多晶宝石进行物相鉴定、宝石产地特征研究以及宝石矿物多型的种类划分等。例如根据XRD图谱,可将品种繁多的微晶石英隐晶质变种玉髓中不同的SiO₂质矿物相精确区分;通过黏土矿物成分及其含量可将鸡血石中"地"的种类进行划分。②宝石矿物的晶体结构研究:对成分复杂的单晶宝石进行晶体化学分析、获取矿物类质同象替代情况、利用结晶度划分宝石品质。例如晶格参数中的c/a比值可将绿柱石的Al八面体类质同象替代和Be四面体类质同象替代进行区分;软玉中通过XRD数据计算出的结晶度与软玉的品质关系密切。③宝石学研究技术的开发:对宝石的优化处理过程进行监测,在宝石合成过程中检验矿物晶形大小、结晶质量及内应力等,以有机宝石的晶体结构为基础的人工养殖技术研究。例如尖晶石在热处理过程中,Mg离子由四面体中的T位迁移至八面体中的M位,导致M-O和T-O键长随温度变化,并反映在尖晶石的色调中;在托帕石合成过程中,XRD数据显示770~800℃时形成含氟托帕石(与天然托帕石结构相同),而在950℃和1000℃出现含氟托帕石离解成刚玉和莫来石。目前,我国在宝石学的XRD应用方面起步较晚,国外在宝石鉴定、优化与合成方面的应用相对成熟。本文认为,该技术与扫描电镜、质子探针等分析技术的联用仍有很大发展空间。     

Modeling discontinuous rock mass based on smoothed finite element method

This paper aims at developing a method for modeling rock mass with preexisting multiple discontinuities within the framework of the smoothed finite element method (SFEM). The discontinuity is simulated by an interface element with zero thickness, the stiffness matrix of which are derived explicitly based on the SFEM. An elastic damage constitutive relation with residual strength is introduced in order to describe the nonlinear mechanical behavior of the discontinuities. The computation codes of the present method were developed. The present method has been verified to be a sound approach for modeling discontinuous rock mass, inheriting the advantages of the SFEM.     

风暴分类识别技术在人工防雹中的应用

利用新一代多普勒天气雷达资料,在风暴跟踪识别算法的基础上,发展了风暴分类技术,以提高人工防雹作业指挥的效率。首先以SCIT算法为基础,结合风暴的结构特征,综合利用雷达、探空资料,自动提取风暴结构特征指数;其次采用基于决策树模型的风暴自动分类技术,将风暴按强度分为雷雨云、单体风暴、多单体风暴和强风暴;最后根据风暴强度、高度和位置等属性,对有可能产生冰雹的单体,结合GIS,自动对下游方向或附近作业点进行预警或输出作业参数。通过对2006—2014年期间重庆、辽宁大连和河南三门峡三地发生的较为典型的31次冰雹天气过程、182站次冰雹样本的检验来看:该方法通过对风暴按强度、垂直结构等综合属性进行分类,能有效提高冰雹识别的命中率、降低空报率,其中强风暴的命中率能达到100%,空报率仅为11.4%。能有效提高人工防雹作业的自动化程度,对防雹作业的科学决策有着重要参考作用。     

利用双线偏振雷达分析人工防雹作业效果

以北京市平谷区2011年6月11日防雹作业为例,利用双线偏振雷达资料,选取与作业区条件十分接近的云块为对比区,根据回波移动方向和速度,跟踪分析在不同高度上防雹作业前后云体宏观结构特征和粒子相态等微物理变化过程。作业后云体所呈现特征为:1作业区的云顶高度、强回波中心高度迅速下降,对比区变化不明显;2作业区水平反射率Zh减小,差分反射率Zdr、零相关系数ρhv增大,单位差分传播相移Kdp小范围内波动,对比区Zh、Zdr、Kdp变化不明显;ρhv增大;3作业区对流减弱,高层较大冰雹粒子、大雨滴下沉明显,最终以霰粒子为主;而对比区域则对流仍然旺盛,冰雹粒子有增多趋势。以上特征表明防雹作业可有效抑制冰雹胚胎成长为冰雹的过程,通过偏振雷达观测参量可对防雹作业效果进行较好的验证。