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.     

关于岩石圈有效弹性厚度的地质理解

简要回顾了岩石圈均衡理论的发展及岩石圈区域均衡和挠曲理论在岩石圈动力学研究中起的作用,阐述岩石圈有效弹性厚度(Te)的概念和特征。强调Te的研究是地质学和地球物理学的紧密结合,即通过岩石圈挠曲理论和区域均衡原理,对地形和重力资料进行谱分析计算,来获取岩石圈的物理性质信息。计算的Te(和相应的挠曲刚度)是岩石圈等效的强度,与爆破地震、地震层析成像和大地电磁测深等方法观测到的岩石圈和地壳厚度不同,它们之间只有通过岩石圈的屈服刚度包络面(YSE)才能比较。大洋和大陆岩石圈YSE的理论计算,表明Te值显著小于地震学的地壳和岩石圈厚度。尤其对于大陆岩石圈,地壳厚度、热年龄和应变率均可显著影响岩石圈的强度。本文还以滇西为实例介绍了对相干值曲线计算的新认识和当前岩石圈Te研究的最新趋势。     

Analytical prediction of aggregates' effects on the ITZ volume fraction and Young's modulus of concrete

In order to evaluate analytically the ITZ volume fraction (f_ITZ) in concrete, a three phase model is proposed for the random concrete microstructure using the Voronoï tessellation. Within this model, the ITZ local thickness is a statistical variable depending on the local paste thickness available between each couple of neighbouring aggregates. The f_ITZ is found to not exceed 7% for typical concretes. Then, the concrete Young's modulus is predicted analytically using a four‐phase generalized self consistent model but in which the proposed f_ITZ is considered. It is found that the concrete Young's modulus increases when increasing aggregates volume fraction, aggregates maximum size and the proportion of coarse aggregates and when decreasing the ITZ thickness and Young's modulus. Finally, the validity of the proposed model is discussed based on a comparison between its predictions and three sets of experimental results related to normal and high strength concretes taken from literature. Copyright © 2016 John Wiley & Sons, Ltd.     

Elastic lateral dynamic impedance functions for a rigid cylindrical shell type foundation

Elastic lateral dynamic impedance functions are defined as the ratio of the lateral dynamic force/moment to the corresponding lateral displacement/rotation at the top ending of a foundation at very small strains. Elastic lateral dynamic impedance functions have a defining influence on the natural frequencies of offshore wind turbines supported on cylindrical shell type foundations, such as suction caissons, bucket foundations, and monopiles. This paper considers the coupled horizontal and rocking vibration of a cylindrical shell type foundation embedded in a fully saturated poroelastic seabed in contact with a seawater half‐space. The formulation of the coupled seawater–shell–seabed vibration problem is simplified by treating the shell as a rigid one. The rigid shell vibration problem is approached by the integral equation method using ring‐load Green's functions for a layered seawater‐seabed half‐space. By considering the boundary conditions at the shell–soil interface, the shell vibration problem is reduced to Fredholm integral equations. Through an analysis of the corresponding Cauchy singular equations, the intrinsic singular characteristics of the problem are rendered explicit. With the singularities incorporated into the solution representation, an effective numerical method involving Gauss–Chebyshev method is developed for the governing Fredholm equations. Selected numerical results for the dynamic contact load distributions, displacements of the shell, and lateral dynamic impedance functions are examined for different shell length–radius ratio, poroelastic materials, and frequencies of excitation. Copyright © 2016 John Wiley & Sons, Ltd.     

Inverse modelling using PS‐InSAR data for improved land subsidence simulation in Las Vegas Valley,Nevada

Las Vegas Valley has had a long history of groundwater development and subsequent surface deformation. InSAR interferograms have revealed detailed and complex spatial patterns of subsidence in the Las Vegas Valley area that do not coincide with major pumping regions. This research represents the first effort to use high spatial and temporal resolution subsidence observations from InSAR and hydraulic head data to inversely calibrate transmissivities (T), elastic and inelastic skeletal storage coefficients (S_ke and S_kv) of the developed‐zone aquifer and conductance (CR) of the basin‐fill faults for the entire Las Vegas basin. The results indicate that the subsidence observations from InSAR are extremely beneficial for accurately quantifying hydraulic parameters, and the model calibration results are far more accurate than when using only groundwater levels as observations, and just a limited number of subsidence observations. The discrepancy between distributions of pumping and greatest levels of subsidence is found to be attributed to spatial variations in clay thickness. The Eglington fault separates thicker interbeds to the northwest from thinner interbeds to the southeast and the fault may act as a groundwater‐flow barrier and/or subsidence boundary, although the influence of the groundwater barrier to this area is found to be insignificant. Copyright © 2016 John Wiley & Sons, Ltd.     

Displacement analytical prediction of shallow tunnel based on unified displacement function under slope boundary

For slope condition of ground surface, the asymmetrical deformation about the vertical center line and the horizontal center line of the tunnel cross section can be formed. A unified displacement function expressed by the Fourier series is presented to express the asymmetrical deformation of the tunnel cross section. Five basic deformation modes corresponding to the expansion order 2 are a complete deformation mode to reflect deformation behaviors of the tunnel cross section under slope boundary. Such this complete displacement mode is implemented into the complex variable solution for analytically predicting tunneling-induced ground deformation under slope boundary. All of these analytical solutions are verified by good agreements of the comparison between the analytical solutions and finite element method results. A parameter study is carried out to investigate the influence of deformation modes of the tunnel cross section, geometrical conditions of the tunnel and the slope angle, and “Buoyancy effect” on the displacement field. Finally, the proposed method is consistent with measured data of the Hejie tunnel in China qualitatively. The presented solution can provide a simplified indication for evaluating the ground deformation under slope condition of ground surface.     

A multidirectional semi‐analytical method for analysis of laterally loaded pile groups in multi‐layered elastic strata

A semi‐analytical method for calculating the response of single piles and pile groups subjected to lateral loading is developed in this paper. Displacements anywhere in the soil domain are tied to the displacements of the piles through decay functions. The principle of virtual work and the calculus of variations are used to derive the governing differential equations that describe the response of the piles and soil. The eigenvalue method and the finite difference technique are used to solve the system of coupled differential equations for the piles and soil, respectively. The proposed method takes into account the soil surface displacement along and perpendicular to the loading direction and produces displacement fields that are very close to those produced by the finite element method but at lower computational effort. Compared with the previous method that considered only the soil displacement along the loading direction, accounting for the multi‐directional soil displacement field produces responses for the piles and soil that are closer to those approximated by the finite element method. Copyright © 2016 John Wiley & Sons, Ltd.