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.     

Effects of temporal scales and space mismatches on the TRMM 3B42 v7 precipitation product in a remote mountainous area

The accurate measurement of precipitation is essential to understanding regional hydrological processes and hydrological cycling. Quantification of precipitation over remote regions such as the Tibetan Plateau is highly unreliable because of the scarcity of rain gauges. The objective of this study is to evaluate the performance of the satellite precipitation product of tropical rainfall measuring mission (TRMM) 3B42 v7 at daily, weekly, monthly, and seasonal scales. Comparison between TRMM grid precipitation and point‐based rain gauge precipitation was conducted using nearest neighbour and bilinear weighted interpolation methods. The results showed that the TRMM product could not capture daily precipitation well due to some rainfall events being missed at short time scales but provided reasonably good precipitation data at weekly, monthly, and seasonal scales. TRMM tended to underestimate the precipitation of small rainfall events (less than 1 mm/day), while it overestimated the precipitation of large rainfall events (greater than 20 mm/day). Consequently, TRMM showed better performance in the summer monsoon season than in the winter season. Through comparison, it was also found that the bilinear weighted interpolation method performs better than the nearest neighbour method in TRMM precipitation extraction.     

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.     

Efficient estimation of flood forecast prediction intervals via single‐ and multi‐objective versions of the LUBE method

Prediction intervals (PIs) are commonly used to quantify the accuracy and precision of a forecast. However, traditional ways to construct PIs typically require strong assumptions about data distribution and involve a large computational burden. Here, we improve upon the recent proposed Lower Upper Bound Estimation method and extend it to a multi‐objective framework. The proposed methods are demonstrated using a real‐world flood forecasting case study for the upper Yangtze River Watershed. Results indicate that the proposed methods are able to efficiently construct appropriate PIs, while outperforming other methods including the widely used Generalized Likelihood Uncertainty Estimation approach. Copyright © 2016 John Wiley & Sons, Ltd.     

Granular element method for three‐dimensional discrete element calculations

This paper endows the recently‐proposed granular element method (GEM) with the ability to perform 3D discrete element calculations. By using non‐uniform rational B‐Splines to accurately represent complex grain geometries, we proposed an alternative approach to clustering‐based and polyhedra‐based discrete element methods whereby the need for complicated and ad hoc approaches to construct 3D grain geometries is entirely bypassed. We demonstrate the ability of GEM in capturing arbitrary‐shaped 3D grains with great ease, flexibility, and without excessive geometric information. Furthermore, the applicability of GEM is enhanced by its tight integration with existing non‐uniform rational B‐Splines modeling tools and ability to provide a seamless transition from binary images of real grain shapes (e.g., from 3D X‐ray CT) to modeling and discrete mechanics computations.Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

单片卫星影像在航测测图中的利用

在边境测图过程中,经常遇到较大面积云层问题,使航测测图无法顺利进行。能充分利用的影像资源是第一次全国地理国情普查的影像,有国产高分2号、SPOT-6、资源3号等单片卫星影像,无立体像对影像。通过对国情普查的影像分析选取,选取一试验区域进行试验,用国产高分2号卫星影像与SPOT-6卫星影像组成立体像对测制地物地貌,用Worldview卫星立体像对影像测制同一区域的地物地貌,将两种影像源组成立体模型测制的地物地貌进行比较及精度统计,说明单片卫星影像在航测测图中达到的精度。     

GPS精密卫星钟差估计研究

随着精密单点定位技术的发展,对于精确的卫星坐标以及卫星钟差改正精度的要求越来越高,精密卫星星历以及精密卫星钟差的求解成为制约精密单点定位技术发展的瓶颈。本文基于修复周跳的载波相位观测值与相位平滑伪距观测值,采用无电离层延迟星间单差精密卫星钟差估计模型,在先估计出整周模糊度后,进行了精密卫星钟差的估计,并采用与IGS事后精密钟差作二次差的方法进行精度分析,这对于提高精密单点定位精度具有一定的意义。     

高分二号卫星影像融合方法探析

图像融合技术是卫星影像处理过程中的一个重要步骤,它可以将全色影像的高空间分辨率特点与多光谱影像的高光谱分辨率特点有机地结合起来,使得融合后的影像更利于后续解译。由于融合处理方法的原理以及获取影像的传感器特性不同,不同的卫星影像适合于不同的融合方法。针对高分二号卫星影像的特点,本文采用Gram-Schmidt融合法、FIHS变换方法、HPF方法、Pansharp方法进行试验。实验结果表明:HPF方法所得融合结果在光谱特性保持方法性能最优,而Pansharp方法融合结果在信息量的保持以及清晰度方面效果更优,因此,在实际应用中,应根据具体需求进行方法选择。