An approach to model the mechanical behavior of transversely isotropic materials

The present work proposes an approach to adapt existing isotropic models to transversely isotropic materials. The main idea is to introduce equivalence relations between the real material and a fictitious isotropic one on which one can take all the advantages of the well‐established isotropic theory. Two applications of this approach are presented here: a failure criterion and a damage model that takes into account the load‐induced anisotropy. In both cases, theoretical predictions are in agreement with the experimental data. In the present paper, the developed approach is applied to sedimentary rock materials; nevertheless, it can be generalized to any material that exhibits transverse isotropy. Copyright © 2015 John Wiley & Sons, Ltd.     

Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

One-dimensional consolidation of soil under multistage load based on continuous drainage boundary

In engineering practice, a rapid loading rate can result in ground failure when the strength of soft soils is relatively low, and a multistage loading scheme is always utilized to deal with this situation. Firstly, under a multistage load and the continuous drainage boundary, an analytical solution of excess pore-water pressure and consolidation degree is obtained by virtue of the superposition formula of excess pore-water pressure, and a more general continuous drainage boundary under arbitrary time-dependent load is developed. Then, a comparison with existing analytical solutions is conducted to verify the present solution. A preliminary attempt on applying the continuous drainage boundary into the finite element model is made, and the feasibility of the numerical model for the one-dimensional consolidation under the continuous drainage boundary is verified by comparing the results calculated by FEM with that from present analytical solution. Finally, the consolidation behavior of soil is investigated in detail for different int erface parameters or loading scheme. The results show that, in land reclamation projects, a horizontal drain should be placed close to the boundary with a smaller interface parameter to improve the consolidation efficiency. The degree of consolidation is also related to the applied time-dependent load and interface parameters.     

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.     

On the relationship between historical land‐use change and water availability: the case of the lower Tarim River region in northwestern China

Natural ecosystems in the region of the lower Tarim River in northwestern China strongly deteriorated since the 1950s due to an expanding desertification. As a result, the downstream Tarim River reaches became permanently dry land. This historical evolution in land‐use change is typically the result of the anthropogenic impact on natural ecosystems. On the basis of a spatially distributed hydrological catchment model bidirectionally linked with a fully hydrodynamic MIKE11 river model, land‐use changes characterized by historical changes in leaf area index (LAI) of vegetation, as well as the evolution of irrigated surface areas, can be causally related to changes in water resources (groundwater storage and surface water resources). An increased surface area of irrigated (agricultural) land, together with a majority of inefficient irrigation methods, did lead to a strong increase of water resources consumption of the farmlands located in the upper Tarim River area. Evidently, this evolution influenced available water resources downstream in the Tarim basin. As a result, farmland has been gradually relocated to the upstream regions. This has led to reduced flows from the upper Tarim stream, which subsequently accelerated the dropping of the groundwater level downstream in the basin. This study moreover demonstrates that land surface biomass changes (cumulative LAI) along the lower Tarim River are strongly related to the changes in groundwater storage. Copyright © 2012 John Wiley & Sons, Ltd.     

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

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

2000～2015年会理县地表覆盖变化监测

以LandSat多时相影像为数据源,采用基于规则的面向对象分类方法,监测会理县2000～ 2015年地表覆盖变化.结合1∶50 000 DEM数据,分析了会理县地表覆盖构成与转移、变化动态度及潜在退耕地情况.监测分析结论能为会理县加强空间管制、促进土地资源集约利用、指导生态文明建设等提供决策信息支撑.