The reflection of plane waves in a poroelastic half-space saturated with inviscid fluid

This paper discusses surface displacements, surface strain, rocking, and energy partitioning during reflection-of-plane waves in a fluid-saturated poroelastic half-space. The medium is modeled by Biot's theory, and is assumed to be saturated with inviscid fluid. A linear porosity-modulus relation based on experimental data on sandstones is used to determine the material parameters for Biot's model. Numerical results in terms of angle of incident waves and Poisson's ratio are illustrated for various porosities and degrees of solid frame stiffness. The results show that the amount of solid frame stiffness controls the response of a fluid-saturated porous system. A poroelastic medium with essentially dry-frame stiffness behaves like an elastic medium, and the influence of pore fluid increases as dry-frame stiffness is reduced. The effects of a second P-wave become noticeable in poroelastic media with low dry-frame stiffness.     

Analytical exploration of γ-function explicit method for pseudodynamic testing of nonlinear systems

It has been well studied that the γ-function explicit method can be effective in providing favorable numerical dissipation for linear elastic systems. However, its performance for nonlinear systems is unclear due to a lack of analytical evaluation techniques. Thus, a novel technique is proposed herein to evaluate its efficiency for application to nonlinear systems by introducing two parameters to describe the stiffness change. As a result, the numerical properties and error propagation characteristics of the γ-function explicit method for the pseudodynamic testing of a nonlinear system are analytically assessed. It is found that the upper stability limit decreases as the step degree of nonlinearity increases; and it increases as the current degree of nonlinearity increases. It is also shown that this integration method provides favorable numerical dissipation not only for linear elastic systems but also for nonlinear systems. Furthermore, error propagation analysis reveals that the numerical dissipation can effectively suppress the severe error propagation of high frequency modes while the low frequency responses are almost unaffected for both linear elastic and nonlinear systems.     

Simplified Method for the Assessment of the Stiffness Anisotropy of Rocks at Small Strains

Summary A new, practically applicable method for characterizing the stiffness anisotropy of rocks is presented. The anisotropy of geo-materials is often ignored in engineering applications, with potentially serious ramifications, because of the number of parameters required for characterization. The elastic anisotropy has often been considered to be a function of mathematical symmetry, and the restrictions due to layering, microcracking and granularity of the materials have not been considered in the assessment of the anisotropy. The practicality of the method proposed here is achieved by rationally reducing the number of independent anisotropy parameters, typically 9 for orthotropic anisotropy, to a system of 4 independent parameters through a systematic theoretical and experimental analysis of these structural restrictions. These 4 parameters are shown to be sufficient for describing the anisotropy of some rocks and sands at small strains, and parameter determination by back-analysis is demonstrated to be stable using appropriate measurement systems involving 9 elastic wave velocities even when the directions of anisotropic axes are unknown and the velocity data contains appreciable error.     

On size and boundary effects in scaled model blasts spatial problems

This paper addresses size and boundary effects on wave propagation, fracture pattern development and fragmentation in small scale laboratory-size specimens for model blasting. Small block type specimens are centre-line loaded by linear explosive charges and supersonically detonated. Using elastic wave propagation theory and fracture mechanics it is shown that the type of boundary conditions which prevail at the outer boundary of the cylinder control the extension of bore-hole cracking and fragmentation within the body of the cylinder. In the case of a composite block where a cylindrical core of different material is embedded, the level of fracturing and fragmentation is controlled by the separation of the interface which in turn depends on the relative dimensions of the core and the block. The most important parameter is the ratio between the length of the pulse (space-wise or time-wise) and the characteristic dimensions of the models, i.e. in this case the dimensions of the core and the mantel. Stress wave superposition effects occur in the corner sections of the mantel. Theoretical results are in good agreement with recent experimental findings.     

Prediction of creep characteristic of rock under varying environment

The strain developed due to creep is mainly proportional to the logarithm of the time under load, and is mostly proportional to the stress and temperature. At higher temperature the creep rate falls slowly with respect to time, and the creep strain is proportional to a fractional power of time, with the exponent increasing as the temperature increases and reaching a value approximately one-third at temperatures of about 0.5°C. At these temperatures, the creep increases with stress according to a power greater than unity and possibly exponentially. It increases with temperature as (−U/kT), where U is an activation energy and k is Boltzman’s constant. There are different methods to determine the creep strain and the energy of Jog (B) including experimental methods, multivariate regression analysis, and by numerical simulation. These methods are less cumbersome and time consuming. In the present investigation, artificial neural network technique has been used for prediction of the creep strain and energy of Jog (B). Two different networks have been tested and validated. Both the networks have four input neurons and one hidden layer with five neurons, and one output neuron. The data for different rocks at temperatures up to 750°C under conditions of compressive or tortional stress are taken from the literatures. The training and testing data sets used were 163 and 14, respectively. To deal with the problem of overfitting of data, Bayesian regulation has been used and network is trained with suitable training epochs. The coefficients of correlation among the predicted and observed values are found high and they improve the confidence of the users. The mean absolute percentage error obtained are also very low.     

Snow loads on wire mesh and cable net rockfall slope protection systems

Snow load on mesh systems is complicated by many factors. This paper presents field instrumentation data on snow load variation with temperature, snowfall and snow depth on a mesh system. It was found that snow load pattern on mesh systems changed with temperature even without variation in snow depth. It reached its maximum value when the temperature rose just above freezing to melt the interface. The field data was used to formulate appropriate snow load models for the various conditions of temperature in the field. The snow load models were used to study the performance of a number of mesh systems in North America and estimate the interface friction that was prevalent for the different surface conditions.     

中国陆地净初级生产力的季节变化研究

了解不同季节陆地净初级生产力(NPP)的变化及与气候的相互关系以及在不同地类的差异对深刻理解我国陆地生态系统对全球气候变化的响应和陆地碳循环研究具有重要意义。本文使用1981～2000年间GLO-PEM模型模拟的我国陆地NPP数据和同期气温、降水以及土地利用数据, 研究不同季节我国陆地植被NPP的变化。结果表明,在1981～2000年期间,四个季节的NPP都呈显著增加趋势,春季是NPP增加速率最快的季节,夏季是NPP增加量最大的季节。耕地在春、夏和秋季NPP增长和增长率最高,林地冬季NPP增长最多而水域冬季NPP增长率最高。夏季NPP增长最高的区域分布于我国东部的多数地区、内蒙古东部、四川盆地、贵州东部、藏南和新疆西部;夏季NPP降低最多的区域分布于在呼伦贝尔高原、鄂尔多斯高原、黄土高原、青藏高原东部和新疆西北部。     

Validating the ability of the discontinuous deformation analysis method to model normal P‐wave propagation across rock fractures

The effects of fractures on wave propagation problems are increasingly abstracting the attention of scholars and engineers in rock engineering field. This study aims to fully validate the ability of discontinuous deformation analysis (DDA) to model normal P‐wave propagation across rock fractures. The effects of a single fracture and multiple parallel fractures are all tested. The results indicate that DDA can accurately reflect the fracture effects, including the fractures stiffness, the fracture spacing and the fracture number, and the effects of incident wave frequency on one‐dimensional P‐wave propagation problems. Thus, DDA is able to deal well with normal incident P‐wave propagation problems. Copyright © 2017 John Wiley & Sons, Ltd.     

Discrete element model for crack propagation in brittle materials

We propose a discrete element model for brittle rupture. The material consists of a bidimensional set of closed‐packed particles in contact. We explore the isotropic elastic behavior of this regular structure to derive a rupture criterion compatible to continuum mechanics. We introduce a classical criterion of mixed mode crack propagation based on the value of the stress intensity factors, obtained by the analysis of two adjacent contacts near a crack tip. Hence, the toughness becomes a direct parameter of the model, without any calibration procedure. We verify the consistency of the formulation as well as its convergence by comparison with theoretical solutions of tensile cracks, a pre‐cracked beam, and an inclined crack under biaxial stress. Copyright © 2015 John Wiley & Sons, Ltd.