非饱和土地基中复合多层波阻板对S波的隔离效应

基于弹性波在非饱和多孔介质与单相弹性介质中的传播理论,考虑在非饱和土地基中设置一定厚度的复合多层波阻板（复合多层波阻板以3层为例）,利用Helmholtz矢量分解定理,推导了非饱和土地基中S波通过复合多层波阻板的透射、反射振幅比的解析解。通过数值算例,分析了层间波阻板剪切模量和密度等物理力学参数对非饱和土地基中S波通过复合多层波阻板时传播特性的影响规律。结果表明：复合多层波阻板中层间波阻板材料的剪切模量和密度对透反射系数影响显著。复合多层波阻板是一种有效的隔振屏障,严格控制层间波阻板的剪切模量和密度可以获得最佳隔振效果,这为复合多层波阻板在地基振动控制领域中的应用提供理论指导。     

基于薄层法和移动坐标系法的地铁竖向振动在成层土层中传播

针对地铁竖向振动在成层地基中的传播,提出了移动荷载作用下层状地基的分析模型。基于该模型,利用狄拉克函数及三重傅里叶变换将时空域内单个移动简谐荷载转换为频率-波数域内的荷载,结合薄层法和移动坐标系法推导了单个移动简谐荷载和移动简谐线荷载作用下三维层状地基动力响应的解析解,并给出了移动简谐线荷载动力响应解析解中参数n的经验取值范围;分析了荷载的移动速度对层状地基动力响应的影响以及弹性模量、泊松比、阻尼比和荷载频率对土体临界速度的影响。结果表明：荷载的移动速度对不同频带的动力响应的影响范围不同;移动速度对低频响应的影响程度大于对高频响应的影响程度;相比于泊松比和阻尼比,弹性模量对土体临界速度的影响最大;频率越接近荷载振动频率的振动响应,其幅值越大,临界速度越小。     

Vertical deformation of rigid foundations of arbitrary shape on layered soil media

A numerical procedure is described for the analysis of vertical deformation of smooth, rigid foundations of arbitrary shape on homogeneous and layered soil media. The contact area at the interface of the foundation and soil medium is approximated by square subdivisions. The response of the system is then obtained from the superposition of the influence of the individual subdivisions. The flexibility influence coefficients are based on equivalent smooth, rigid circular areas with the same contact area as the square subdivisions. For foundations on a homogeneous, isotropic elastic half-space, the flexibility coefficients are given analytically by the integrated forms of the Boussinesq's solution. For a layered soil medium, the flexibility coefficients are determined from an axisymmetric finite element analysis which is essentially two dimensional. Thus, there is no necessity for a full three-dimensional finite element analysis. Comparison with solutions obtained using the integral transform technique for smooth, rigid rectangular foundations on a homogeneous, isotropic elastic half-space shows good agreement. Parametric solutions are presented for the response of rectangular foundations on some ‘typical’ soil profiles. The use of a simplified method to estimate the settlement of rectangular foundations on a layered soil medium by superposing solutions for homogeneous, elastic strata is discussed.     

Monitoring in situ stress/strain behaviour during plastic yielding in polymineralic rocks using neutron diffraction

Attempts to use rock deformation experiments to examine the elastic and plastic behaviour of polymineralic rocks are hampered by the fact that usually only whole sample properties can be monitored as opposed to the separate contribution of each phase. To circumvent this difficulty, room-temperature, uniaxial compression experiments were performed in a neutron beam-line on a suite of calcite + halite samples with different phase volume proportions. By collecting diffraction data during loading, the elastic strain and hence stress in each phase was determined as a function of load to bulk strains of 1–2%. In all cases, the calcite behaved elastically while the halite underwent plastic yielding. During the fully elastic part of the deformation, the composite elastic properties and the within-phase stresses are well-described both by recent shear lag models and by analyses based on Eshelby's solution for the elastic field around an ellipsoidal inclusion in a homogeneous medium. After the onset of yielding, the halite in situ stress/total strain curve may be reconstructed using the rule of mixtures. At calcite contents of greater than 30%, the in situ halite response may be significantly weaker or stronger than that obtained at lesser calcite contents. The results highlight the potential that such techniques offer for developing an explicitly experimental approach for determining the influence of microstructural variables on the mechanical properties of polymineralic rocks.     

Effect of rock discontinuities on certain rock strength and fracture energy parameters under uniaxial compression

Summary Five series of test blocks of Pendeli marble with artificially created discontinuities of different crack densities (simulating three mutually orthogonal joint sets) were tested in uniaxial compression in order to study the effect of discontinuities on: (a) the compressive strength and the modulus of elasticity, and (b) certain fracture energy parameters expressed by the ratio W _A/W _V, where W _A is the surface energy and W _V the volume elastic strain energy. Mathematical relationships are derived similar to those suggested by other authors relating strength parameters to crack densities. Such relationships clearly show a reduction in strength with increased crack density. The experimental results obtained permit the extension of Persson's relation (which refers to ideal intact rock) to the more realistic case of discontinuous rock mass by introducing the appropriate term that takes into consideration the effect of rock mass discontinuities on the energy ratio W _A/W _V. A comparison between laboratory results and field observations was subsequently carried out assuming the rock mass to behave as a linearly elastic material, obeying the Hoek and Brown failure criterion. This comparison showed that laboratory results can be extended to larger scale. Furthermore, in order to predict the in situ strength and stability of a rock mass in uniaxial compression (which is of major importance in underground excavations) certain concepts are proposed based on laboratory tests, in situ investigations and first principles of linear elastic fracture mechanics.     

A semi-discrete procedure for dynamic response analysis of saturated soils

Biot's equations of wave propagation through fluid-saturated porous elastic media are discretized spatially using the finite element method in conjunction with Galerkin's procedure. Laplace transformation of the discretized equations is used to suppress the time variable. Introducing Laplace transforms of constituent velocities at nodal points as additional variables, the quadratic set of equations in the Laplace transform parameter is reduced to a linear form. The solution in the Laplace transform space is inverted, term by term, to get the complete time history of the solid and fluid displacements and velocities. Since the solution is exact in the time domain, the error in the calculated response is entirely due to the spatial approximation. The procedure is applied to one-dimensional wave propagation in a linear elastic material and in a fluid-saturated elastic soil layer with ‘weak’, ‘strong’ as well as ‘moderate’ coupling. With refinement of the spatial mesh, convergence to the exact solution is established. The procedure can provide a useful benchmark for validation of approximate temporal discretization schemes and estimation of errors due to spatial discretization.     

Transient response of a circular cavity in a poroelastic medium

This paper considers the transient response of a pressurized long cylindrical cavity in an infinite poroelastic medium. To obtain transient solutions, Biot's equations for poroelastodynamics are specialized for this problem. A set of exact general solutions for radial displacement, stresses, pore pressure and discharge are derived in the Laplace transform space by using analytical techniques. Solutions are presented for three different types of prescribed transient radial pressures acting on the surface of a permeable as well as an impermeable cavity surface. Time domain solutions are obtained by inverting Laplace domain solutions using a reliable numerical scheme. A detailed parametric study is presented to illustrate the influence of poroelastic material parameters and hydraulic boundary conditions on the response of the medium. Comparisons are also presented with the corresponding ideal elastic solutions to portray the poroelastic effects. It is noted that the maximum radial displacement and hoop stress at the cavity surface are substantially higher than the classical static solutions and differ considerably from the transient elastic solutions. Time histories and radial variations of displacement, hoop stress, pore pressure and fluid discharge corresponding to a cavity in two representative poroelastic materials are also presented.     

Load transfer from an elastic pile to a saturated porous elastic soil

The quasi-static development of the force, pore pressure and displacement is obtained in the system of a circular elastic pile partially embedded in a saturated porous elastic soil, and loaded axially on the top. The porous elastic soil is governed by Biot's theory. The problem is decomposed into two systems, namely, an extended porous elastic half-space in the absence of the pile characterized by the material constants of the medium, and a fictitious pile represented by a Young's modulus equal to the difference between the Young's moduli of the real pile and the medium. The problem is found to be governed by a Fredholm integral equation of the second kind. Laplace transforms are applied to time functions involved, and Hankel transforms to the radial coordinate of which the origin is at the centre of the pile. Numerical solutions are obtained for final and initial solutions for various practical values of the parameters involved.     

Analysis of stone-column reinforced foundations

A numerical model is proposed to analyse elastic as well as elastoplastic behaviour of stone-column reinforced foundations. The stone-columns are assumed to be dispersed within the in situ soil and a homogenization technique is invoked to establish equivalent material properties for in situ soil and stone-column composite. The difficulties encountered in carrying out elastoplastic analyses of composite materials are overcome by adopting a separate yield function for each of the constituent materials and a sub-iteration procedure within an implicit backward Euler stress integration scheme. In the proposed procedure, equilibrium as well as kinematic conditions implied in the homogenization procedure are satisfied for both elastic as well as elastoplastic stress states. The proposed model is implemented in an axi-symmetric finite element code and numerical prediction is made for the behaviour of model circular footings resting on stone-column reinforced foundations. This prediction indicates good agreement with experimental observation. Finally, a new scheme in which the length of stone-column is variable is proposed and its behaviour is examined through a numerical example. © 1998 John Wiley & Sons, Ltd.     

Parallel, similar and constrained folds

Theoretical analysis of folding of viscous multilayers with free slip or bonding at layer contacts indicates that folds in such multilayers can be described in terms of three end-members:parallel, in which orthogonal thicknesses of layers are largely constant;similar, in which vertical thicknesses of layers and shapes of successive interfaces are essentially constant; andconstrained, in which amplitudes of anticlines and synclines decrease to zero at upper and lower boundaries. Constrained,internal folds form if the multilayer is confined by rigid media; parallel,concentric-like folds form if the multilayer is confined by soft media, provided soft interbeds are sufficiently thin for the stiff layers to fold as an ensemble. Similar,sinusoidal orchevron folds form throughout much of the thickness of a multilayer, for any stiffness of confining media, provided wavelengths of folds are short relative to the thickness of the multilayer or soft interbeds are sufficiently soft and thick for the stiff layers to act independently. The analysis shows that multilayer folds may have the same form regardless of whether the layer contacts are freely slipping or bonded.

The forms of folds in multilayers confined by media with different viscosities above and below depend on the viscosity contrast of the media. For no medium above and a rigid medium below, the forms are concentric-like in the upper part and internal in the lower part of the multilayer. For no medium above and a soft medium below, the folds are concentric-like throughout the multilayer.

The theory indicates that a useful way to analyze forms of folds in rocks or in experiments is in terms of component waveforms, as defined, for example, by Fourier series. The distributions of amplitudes of component waveforms throughout the multilayer appears to be diagnostic, reflecting contrasts in properties of the multilayer and its confining media. Analysis of a large fold in the central Appalachians, Pennsylvania, and of a smaller fold in the Huasna syncline, California, indicates that at least three component waveforms are required to produce the gross forms of those folds.

The theory closely predicts wavelengths and shapes of folds produced in analogous elastic multilayers, indicating that nonlinearities in material behavior, which are inherent in the elastic material but are absent in the viscous material, are less significant than nonlinearities in the boundary conditions, which are the same in elastic and viscous materials.     

Analysis of pavement response to subsurface deformations

The paper deals with modeling the response of pavement systems to subsurface deformations. An analytic formulation is developed in the context of layered elastic theory wherein the bottom boundary is distorted in a prescribed manner – forcing the medium to deform. In the basic case considered, an axisymmetric displacement field is imposed, consisting of vertical and radial components. A version of the model, dealing with a single elastic layer deformed only in the vertical direction (with frictionless bottom), is parametrically interrogated to study the influence of layer thickness and modulus on the resulting surface profile and vertical contact stresses. Furthermore, it is shown that by spatially superposing several axisymmetric solutions, general displacement shapes can be reproduced and analyzed. Such a scheme is demonstrated in an inverse type of analysis to a case study dealing with a sagging runway; using measured surface elevations, an attempt is made to backcalculate subsurface deformations and stress changes at the bottom interface of the structure.     

Derivation of material variability from settlement measurements

In pavement evaluation using non-destructive testing (NDT), a large amount of deflection bowls are analysed in terms of the elastic moduli of the layers. The results are used to evaluate the material variability, which could serve in an overlay design procedure based on the concept of reliability. The model currently used for interpreting deflection bowls is based on the random variable theory which neglects the spatial distribution of the elastic modulus of the material. Since the subgrade and pavement materials have a spatial distribution, the analysis of NDT could lead to an underestimate of the material variability. The random field theory, which is more adequate than the random variable theory, is presented and used to correct the NDT analysis. The pavement is modelled as a multilayer random elastic medium with, for each layer, a constant Poisson's ratio and a random shear modulus characterized by three statistical moments: average value, standard deviation and autocorrelation function. The stochastic integral formulation presented in the previous publications is generalized here for a multilayer system. The multilayer system is analysed with the random field theory and the covariance matrix of the deflection bowl is obtained and used to generate deflection bowls corresponding to the properties of the random field. These bowls are then interpreted with the current procedure and elastic modulus variabilities are computed. It is found that the current procedure for interpreting deflection bowls underestimates the variability of the subgrade, by a factor of 0.4–1.0. It is interesting to note that the average moduli of the Boussinesq layer and of the two layers are not affected by the type of theory used. The variability of the upper layer in the two-layer system is also unaffected for a small variation range.     

A comparative study of ANN and Neuro-fuzzy for the prediction of dynamic constant of rockmass

Physico-mechanical properties of rocks have great significance in all operational parts in mining activities, from exploration to final dispatch of material. Compressional wave velocity (p-wave velocity) and anisotropic behaviour of rocks are two such properties which help to understand the rock response under varying stress conditions. They also influence the breakage mechanism of rock. There are different methods to determine thep-wave velocity and anisotropyin situ and in the laboratory. These methods are cumbersome and time consuming. Fuzzy set theory, Fuzzy logic and Neural Networks techniques seem very well suited for typical geotechnical problems. In conjunction with statistics and conventional mathematical methods, hybrid methods can be developed that may prove to be a step forward in modeling geotechnical problems. Here, we have developed and compared two different models, Neuro-fuzzy systems (combination of fuzzy and artificial neural network systems) and Artificial neural network systems, for the prediction of compressional wave velocity.     

Dynamic response of concrete gravity dams including dam–water–foundation interaction

In this paper, the gravity dam–water–foundation system including the physical and mechanical properties of the sediment at the reservoir bottom is modelled using a finite element and infinite element coupling model. The sediment at the reservoir bottom has been assumed to be a viscoelastic solid medium. The effects of thickness, elastic modulus, Poisson's ratio and material damping of the sediment on the response of the dam have been studied. The related numerical results from this study illustrated that the existence of the sediment at the bottom of the reservoir has significant effects on the response of concrete gravity dams since the soft layer of the sediment plays two main roles in the dam–water–foundation system, the energy dissipation of the system and the amplification of the incident wave on the water–sediment interface. It is the amplified acceleration on the water–sediment interface that results in different mechanisms of the effect of the sediment on the response of the dam. Therefore, apart from the incident wave, the thickness, the softness and the damping ratio of the sediment can also affect the response of the dam.     

Transverse seismic kinematics of single piles by a modified Vlasov model

Within the framework of soil–pile interaction, a novel displacement scheme for the transverse kinematic response of single piles to vertically propagating S waves is proposed on the basis of the modified Vlasov foundation model. The displacement model contains a displacement function along the pile axis and an attenuation function along the radial direction. The governing equations and boundary conditions of the two undetermined functions are obtained in a coupled form by using Hamilton's principle. An iterative algorithm is adopted to decouple and solve the two unknown functions. In light of the governing equation of the pile kinematics, a mechanical model is proposed to evaluate the present method on a physical basis considering material damping. The coefficient of the equivalent Winkler spring is derived explicitly as function of the displacement decay parameter γ and soil Poisson's ratio. A parametric study is performed to investigate the effects of the soil–pile system properties on the kinematic response of single piles. The results show that the dimensionless pile length controls the transverse kinematics of piles. In terms of the theory of beams on elastic foundation, the classification limits of the dimensionless pile length may be π ∕ 4 and π, respectively. Copyright © 2014 John Wiley & Sons, Ltd.     

Vertical vibration of a rigid strip footing on viscoelastic half-space

This paper presents an analytical method for modeling the dynamic response of a rigid strip footing subjected to vertical-only loads. The footing is assumed to rest on the surface of a viscoelastic half-space; therefore, effects of hysteretic soil damping on the impedance of the foundation and the generated ground vibrations are considered in the solution. To solve the mixed boundary value problem, we use the Fourier transform to cast a pair of dual integral equations providing contact stresses, which are solved by means of Jacobi orthogonal polynomials. The resulting soil and footing displacements and stresses are obtained by means of the Fourier inverse transform. The solution provides more realistic estimates of footing impedance, compared to existing solutions for elastic soil, as well as of the attenuation of ground vibrations with distance of the footing. The latter is important for the estimation of machine vibration effects on nearby structures and installations.     

Poroelastic model for pile–soil interaction in a half‐space porous medium due to seismic waves

In this paper, frequency domain dynamic response of a pile embedded in a half‐space porous medium and subjected to P, SV seismic waves is investigated. According to the fictitious pile methodology, the problem is decomposed into an extended poroelastic half‐space and a fictitious pile. The extended porous half‐space is described by Biot's theory, while the fictitious pile is treated as a bar and a beam and described by the conventional 1‐D structure vibration theory. Using the Hankel transformation method, the fundamental solutions for a half‐space porous medium subjected to a vertical or a horizontal circular patch load are established. Based on the obtained fundamental solutions and free wave fields, the second kind of Fredholm integral equations describing the vertical and the horizontal interaction between the pile and the poroelastic half‐space are established. Solution of the integral equations yields the dynamic response of the pile to plane P, SV waves. Numerical results show the parameters of the porous medium, the pile and incident waves have direct influences on the dynamic response of the pile–half‐space system. Significant differences between conventional single‐phase elastic model and the poroelastic model for the surrounding medium of the pile are found. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical simulation of fully saturated porous materials

Fully coupled, porous solid–fluid formulation, implementation and related modeling and simulation issues are presented in this work. To this end, coupled dynamic field equations with u?p?U formulation are used to simulate pore fluid and soil skeleton (elastic–plastic porous solid) responses. Present formulation allows, among other features, for water accelerations to be taken into account. This proves to be useful in modeling dynamic interaction of media of different stiffnesses (as in soil–foundation–structure interaction). Fluid compressibility is also explicitly taken into account, thus allowing excursions into modeling of limited cases of non‐saturated porous media. In addition to these features, present formulation and implementation models in a realistic way the physical damping, which dissipates energy. In particular, the velocity proportional damping is appropriately modeled and simulated by taking into account the interaction of pore fluid and solid skeleton. Similarly, the displacement proportional damping is physically modeled through elastic–plastic processes in soil skeleton. An advanced material model for sand is used in present work and is discussed at some length. Also explored in this paper are the verification and validation issues related to fully coupled modeling and simulations of porous media. Illustrative examples describing the dynamical behavior of porous media (saturated soils) are presented. The verified and validated methods and material models are used to predict the behavior of level and sloping grounds subjected to seismic shaking. Copyright © 2008 John Wiley & Sons, Ltd.     

Dynamic responses of a poroelastic half‐space from moving trains caused by vertical track irregularities

The vibrations of railway tracks on a poroelastic half‐space generated by moving trains are investigated through a vehicle–track–ground coupling model. The theoretical model incorporates a vehicle, a track, and a fully saturated poroelastic half‐space soil medium. The source of vibration excitation is divided into two components: the quasi‐static loads and the dynamic loads. The quasi‐static loads are related to the static component of the axle loads, whereas the dynamic loads are due to the dynamic wheel–rail interaction. A linear Hertizian contact spring is introduced between each wheelset and the rail to consider the dynamic loads. Biot's dynamic theory is used to characterize the poroelastic half‐space soil medium. Using the Fourier transform, the governing equations for the track–ground system are solved and the numerical results are presented for a single axle vehicle model. The different dynamic characteristics of the elastic soil medium and the saturated poroelastic medium are investigated. In addition, the different roles of the moving axle loads and the roughness‐induced dynamic loads are identified. It is concluded that the vibration level of the free field off the track predicted by the poroelastic soil medium is smaller than that predicted by the elastic soil medium for vehicle speed below the Rayleigh wave speed of the poroelastic half‐space, whereas it is larger for vehicle speed above the Rayleigh wave speed. The dynamic loads play an important role in the dynamic responses of the track–ground system. Copyright © 2010 John Wiley & Sons, Ltd.     

Finite element analysis of consolidation of layered clay soils using an elastic visco-plastic model

This paper presents a general one-dimensional (1-D) finite element (FE) procedure for a highly non-linear 1-D elastic visco-plastic (1-D EVP) model proposed by Yin and Graham for consolidation analysis of layered clay soils. In formulating the 1-D FE procedure, a trapezoidal formula is used to avoid the unsymmetry of the stiffness matrix for a Newton (modified Newton) iteration scheme. Unlike many other 1-D FE approaches in which the initial in situ stresses (or stress/strain states) are considered indirectly or even not considered, the initial in situ stress/strain states are taken into account directly in this paper. The proposed FE procedure is used for analysis of 1-D consolidation of a clay with published test results in the literature. The FE modelling results are in good agreement with the measured results. The FE model and procedure is then used to analyse the consolidation of a multi-layered clay soils with a parametric study on the effects of the variations of creep parameters in Yin and Graham's 1-D EVP model. It is found that the creep parameters ψ/V and t₀ have significant influence on the compression and porewater pressure dissipation. For some boundary conditions, changes of parameters in one layer will have some effects on the consolidation behaviour of another layer due to the different consolidation rates. Finally, the importance of initial stress/strain states is illustrated and discussed. Copyright © 1999 John Wiley & Sons, Ltd.