In-plane and anti-plane strong shaking of soil systems and structures

The concept of in-plane and anti-plane shaking is introduced with a rigid block on a plane surface with Coulomb friction. Using a hypoplastic constitutive relation to model the mechanical behaviour of the soil, numerical solutions for a rigid block on a thin dry or saturated soil layer are obtained. The coupled nature of dynamic problems involving granular materials is shown, i.e. the motion of the block changes the soil state—skeleton stresses and density—which in turn affects the block motion. Motions of the block as well as soil response can be more realistically calculated by the new model. The same constitutive equation is applied to the numerical simulation of the propagation of plane waves in homogeneous and layered level soil deposits induced by a wave coming from below. Experiments with a novel laminar shake box as well as real seismic records from well-documented sites during strong earthquakes are used to verify the adequacy of the hypoplasticity-based numerical model for the prediction of soil response during strong earthquakes. The response of a homogeneous earth dam subjected to in-plane and anti-plane shaking is investigated numerically. In-plane and anti-plane shaking is shown to cause nearly the same spreading of a sand dam under drained conditions, whereas under undrained conditions anti-plane shaking causes stronger spreading of the dam. The dynamic behaviour of a breakwater founded on rockfill and soft clay during the 1995 Kobe earthquake is back-calculated to show the good performance of the proposed numerical model also with a structure. Section 9 deals with buildings on mattresses of densified cohesionless soils or fine-grained soils with granular columns, slopes with ‘hidden’ dams and structures on piles traversing clayey slopes to show the suitability of hypoplasticity-based models for the earthquake-resistant design and safety assessment of geotechnical systems.     

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.     

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.     

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.     

南海东北部内潮的模态结构和传播

The evolution of energy, energy flux and modal structure of the internal tides(ITs) in the northeastern South China Sea is examined using the measurements at two moorings along a cross-slope section from the deep continental slope to the shallow continental shelf. The energy of both diurnal and semidiurnal ITs clearly shows a～14-day spring-neap cycle, but their phases lag that of barotropic tides, indicating that ITs are not generated on the continental slope. Observations of internal tidal energy flux suggest that they may be generated at the Luzon Strait and propagate west-northwest to the continental slope in the northwestern SCS. Because the continental slope is critical-supercritical with respect to diurnal ITs, about 4.6 kJ/m~2 of the incident energy and 8.7 kW/m of energy flux of diurnal ITs are reduced from the continental slope to the continental shelf. In contrast, the semidiurnal internal tides enter the shelf because of the sub-critical topography with respect to semidiurnal ITs.From the continental slope to the shelf, the vertical structure of diurnal ITs shows significant variation, with dominant Mode 1 on the deep slope and dominant higher modes on the shelf. On the contrary, the vertical structure of the semidiurnal ITs is stable, with dominant Mode 1.