A novel method for simulating nuclear explosion with chemical explosion to form an approximate plane wave: Field test and numerical simulation

A nuclear explosion in the rock mass medium can produce strong shock waves,seismic shocks,and other destructive effects,which can cause extreme damage to the underground protection infrastructures.With the increase in nuclear explosion power,underground protection engineering enabled by explosion-proof impact theory and technology ushered in a new challenge.This paper proposes to simulate nuclear explosion tests with on-site chemical explosion tests in the form of multi-hole explosions.First,the mechanism of using multi-hole simultaneous blasting to simulate a nuclear explosion to generate approximate plane waves was analyzed.The plane pressure curve at the vault of the underground protective tunnel under the action of the multi-hole simultaneous blasting was then obtained using the impact test in the rock mass at the site.According to the peak pressure at the vault plane,it was divided into three regions:the stress superposition region,the superposition region after surface reflection,and the approximate plane stress wave zone.A numerical simulation approach was developed using PFC and FLAC to study the peak particle velocity in the surrounding rock of the underground protective cave under the action of multi-hole blasting.The time-history curves of pressure and peak pressure partition obtained by the on-site multi-hole simultaneous blasting test and numerical simulation were compared and analyzed,to verify the correctness and rationality of the formation of an approximate plane wave in the simulated nuclear explosion.This comparison and analysis also provided a theoretical foundation and some research ideas for the ensuing study on the impact of a nuclear explosion.

Scattering of plane P waves by circular-arc layered alluvial valleys: An analytical solution

An analytical solution for scattering of plane P waves by circular-arc layered alluvial valleys was derived by Fourier-Bessel series expansion technique, and the solution was utilized to analyze the effects of alluvial sequence and their relative stiffness on the scattering of incident waves.

Diffraction of plane SV waves by a cavity in poroelastic half-space

This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green＇s functions of compressive and shear wave sources are derived based on Biot＇s theory. The scattered waves are constructed using fictitious wave sources close to the boundary of the cavity, and their magnitudes are determined by the boundary conditions. Verification of the accuracy is performed by： （1） checking the satisfaction extent of the boundary conditions, （2） comparing the degenerated solutions of a single-phased case with well- known solutions, and （3） examining the numerical stability of the solutions. The nature of diffraction of plane SV waves around a cavity in a poroelastic half-space is investigated by numerical examples.

Diffraction of plane P waves by a canyon of arbitrary shape in poroelastic half-space (Ⅰ): Formulation

This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green＇s functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot＇s theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.

Scattering of plane SH waves by a circular-arc hill with a circular tunnel

No abstract available

Anti-plane (SH) waves diffraction by an underground semi-circular cavity:analytical solution

Diffraction of a two-dimensional （2D） semi-circular cavity in a half-space under incident SH-waves is studied using the classic wave function expansion method with a new de-coupling technique. This so-called ＂improved cosine half- range expansion＂ algorithm exhibits an excellent performance in reducing displacement residual errors at two rim points of concern. The governing equations are developed in a manner that minimizes the residues of the boundary conditions. Detailed derivation and analysis procedures as well as truncation of infinite linear governing equations are presented. The semi-circular cavity model presented in this paper, due to its simple profile, is expected to be used in seismic wave propagation studies as a benchmark for examining the accuracies of various analytical or numerical methods for mixed-boundary wave propagation problems.

2.5D scattering of incident plane SV waves by a canyon in layered half-space

This paper presents a 2.5D scattering of incident plane SV waves by a canyon in a layered half-space by using the indirect boundary element method （IBEM）. A free field response analysis is performed to provide the displacements and stresses on the boundary of the canyon where fictitious uniform moving loads are applied to calculate the Green＇s fi.mctions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The free field displacements are added to the fictitious uniform moving loads induced displacements and the total response is obtained. Numerical calculations are performed for a canyon with homogenous and in one layer over bedrock. The effects of the thickness and stiffness of the layer on the amplification are studied and discussed.

Surface motion of a semi-elliptical hill for incident plane SH waves

A closed-form analytical solution of surface motion of a semi-elliptical cylindrical hill for incident plane SH waves is presented. Although some previous analytical work had already dealt with hill topography of semi-circular and shallow circular, our work aims at calculating surface motion of very prolate hill for high incident frequency, and explaining the special vibrating is checked by boundary conditions, numerical results for and some conclusions are obtained. properties of very prolate hill. Accuracy of the solution surface motion of oblate and prolate hills are calculated,

Diffraction of plane P waves by a canyon of arbitrary shape in poroelastic half-space (Ⅱ): Numerical results and discussion

This paper investigates in detail the nature of diffraction of plane P waves around a canyon in poroelastic half-space, and studies the effects of incident frequency, drainage condition, porosity, etc, on the diffraction of waves. It is shown that the surface displacement amplitudes of the drained case are close to those of the undrained case, however, the surface displacement amplitudes of the dry case are very different from those of the saturated （either drained or undrained） cases. There are large phase shift between the dry case and the saturated cases, as well as slightly longer resultant wavelengths for the undrained case than those for the drained case and longer resultant wavelengths for the drained case than those for the dry case. For small porosity the surface displacement amplitudes for the saturated cases are almost identical to those for the dry case; while for large porosity, the effect of drainage condition becomes significant, and the surface displacement amplitudes for the undrained case are larger than those for the drained case. As the incident frequency increases, the effect of porosity becomes significant, and more significant for the undrained case than that for the drained case. As the porosity increases, the pore pressures increase significantly but their oscillations become smoother. As the incident frequency increases, the pore pressures become more complicated.

Scattering of plane P waves by a semi-cylindrical hill: analytical solution

An analytical solution for scattering of plane P waves by a semi-cylindrical hill was derived by using the wave function expansion method, and convergence of the solution and accuracy of truncation were verified. The effect of incident frequency and incident angle on the surface motion of the hill was discussed, and it was shown that a hill greatly amplifies incident plane P waves, and maximum horizontal displacement amplitudes appear mostly at the inclined incidence of waves, which are located at the half-space; and maximum vertical displacement amplitudes emerge mostly at the vertical incidence of waves, which are situated at the hill.

Dynamic soil-tunnel interaction in layered half-space for incident plane SH waves

The dynamic soil-tunnel interaction is studied by indirect boundary element method （IBEM）, using the model of a rigid tunnel in layered half-space, which is simplified to a single soil layer on elastic bedrock, subjected to incident plane SH waves. The accuracy of the results is verified through comparison with the analytical solution. It is shown that soil-tunnel interaction in layered half-space is larger than that in homogeneous half-space and this interaction mechanism is essentially different from that of soil-foundation-superstructure interaction.

Investigation of a silicon-based one-dimensional phononic crystal plate via the super-cell plane wave expansion method

The super-cell plane wave expansion method is employed to calculate band structures for the design of a siliconbased one-dimensional phononic crystal plate with large absolute forbidden bands. In this method, a low impedance medium is introduced to replace the free stress boundary, which largely reduces the computational complexity. The dependence of band gaps on structural parameters is investigated in detail. To prove the validity of the super-cell plane wave expansion, the transmitted power spectra of the Lamb wave are calculated by using a finite element method. With the detailed computation, the band-gap of a one-dimensional plate can be designed as required with appropriate structural parameters, which provides a guide to the fabrication of a Lamb wave phononic crystal.

Tunnel effect of fractal fault and transient S-wave velocity rupture (TSVR) of in-plane shear fault

Transient S wave velocity rupture (TSVR) means the velocity of fault rupture propagation is between S wave velocity α and P wave velocity β . Its existing in the rupture of in plane ( i.e . strike slip) fault has been proved, but in 2 dimensional classical model, there are two difficulties in transient S wave velocity rupture, i.e ., initialization difficulty and divergence difficulty in interpreting the realization of TSVR. The initialization difficulty means, when v ↑ v R (Rayleigh wave velocity), the dynamic stress strength factor K 2(t) →+0, and changes from positive into negative in the interval ( v R, β ). How v transit the forbidden of ( v R, β )? The divergence difficulty means K 2(t) →+∞ when v ↓ β . Here we introduce the concept of fractal and tunnel effect that exist everywhere in fault. The structure of all the faults is fractal with multiple cracks. The velocity of fault rupture is differentiate of the length of the fault respect to time, so the rupture velocity is also fractal. The tunnel effect means the dynamic rupture crosses over the interval of the cracks, and the coalescence of the intervals is slower than the propagation of disturbance. Suppose the area of earthquake nucleation is critical or sub critical propagation everywhere, the arriving of disturbance triggers or accelerates the propagation of cracks tip at once, and the observation system cannot distinguish the front of disturbance and the tip of fracture. Then the speed of disturbance may be identified as fracture velocity, and the phenomenon of TSVR appears, which is an apparent velocity. The real reason of apparent velocity is that the mathematics model of shear rupture is simplified of complex process originally. The dual character of rupture velocity means that the apparent velocity of fault and the real velocity of micro crack extending, which are different in physics, but are unified in rupture criterion. Introducing the above mentioned concept to the calculation of K 2 (t) , the difficulty of initialization can be overcome, and the integral equation of triggering the initialization of TSVR is given quantitatively. By solving this integral equation, the lower limit of TSVR is 1.105 3 β , not β , and the divergence difficulty is overcome. TSVR is unstable solution, and may degenerate to sub Rayleigh wave velocity rupture immediately where the non critical condition can be measured. The results of this paper show that the initialization and continuum depends on the condition of earthquake nucleation in seismogenic area.

Reflection of plane seismic waves at the surface of double-porosity dual-permeability materials

The present work deals with the reflection of plane seismic waves at the stress-free plane surface of double-porosity dualpermeability material. The incidence of two main waves（i.e., P1 and SV） is considered. As a result of the incident waves,four reflected（three longitudinal and one shear） waves are found in the medium. The expressions of reflection coefficients for a given incident wave are obtained as a non-singular system of linear equations. The energy shares of reflected waves are obtained in the form of an energy matrix. A numerical example is considered to calculate the partition of incident energy for fully closed as well as perfectly open pores. Effect of incident direction on the partition of the incident energy is analyzed with the change in wave frequency, wave-induced fluid-flow, pore-fluid viscosity and double-porosity structure. It has been confirmed from the numerical interpretation that during the reflection process, conservation of incident energy is obtained at each angle of incidence.

Diffraction of plane P waves around an alluvial valley in poroelastic half-space

Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green＇s fimctions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of ficti- tious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, and the comparison between the degenerated solutions and available results in single-phase case. Finally, the nature of diffraction of plane P waves around an alluvial valley in poroelastic half-space is investigated in detail through nu- merical examples.

Reflection and refraction of plane waves at the boundary of an elastic solid and double-porosity dual-permeability materials

Phenomena of reflection and refraction of plane harmonic waves at a plane interface between an elastic solid and doubleporosity dual-permeability material are investigated. The elastic solid behaves non-dissipatively, while double-porosity dual-permeability materials behave dissipatively to wave propagation due to the presence of viscosity in pore fluids. All the waves(i.e., incident and reflected) in an elastic medium are considered as homogeneous(i.e., having the same directions of propagation and attenuation), while all the refracted waves in double-porosity dual-permeability materials are inhomogeneous(i.e., having different directions of propagation and attenuation). The coefficients of reflection and refraction for a given incident wave are obtained as a non-singular system of linear equations. The energy shares of reflected and refracted waves are obtained in the form of an energy matrix. A numerical example is considered to calculate the partition of incident energy among various reflected and refracted waves. The effect of incident direction on the partition of the incident energy is analyzed with a change in wave frequency, wave-induced fluid-flow, pore-fluid viscosity and double-porosity structure.It has been confirmed from numerical interpretation that during the reflection/refraction process, conservation of incident energy is obtained at each angle of incidence.

Surface motion of a semi-cylindrical hill for incident plane SV waves: Analytical solution

An analytical solution for surface motion of a semi-cylindrical hill for incident plane SV waves was derived by using the wave function expansion method and the auxiliary function technique, and convergence of the solution and accuracy of truncation were verified. The effect of incident frequency and angle as well as hill width on the surface motion of the hill was discussed by numerical examples. It was shown that, a hill greatly amplifies incident plane SV waves, and the maximal amplification may reach 4 times of that for free-field response; and for incident waves of low frequency, the maximal displacement amplitudes emerge mostly at the half-space, however, for incident waves of high frequency, the maximal displacement amplitudes appear mostly at the hill.

Scattering of plane P waves by canyons containing multiple circular-arc-shaped layers

By expressing the wave functions in the form of Fourier-Bessel series, the analytical solution for the two-dimension scattering problem of plane P waves by the cylindrical canyon topography that contains arbitrary number of circular-arc-shaped layers is presented firstly. And then, the convergence of the proposed series solution with the truncation number of terms is discussed, which demonstrates that the analytical solution can converge even for very high frequencies of the incident P wave. Finally, using this solution, the influences that are imposed on the stationary ground motion by the number and the sequence of alluvial layers, as well as the stiffness of soft interlayer contained in the canyon, are studied,

Calculation of Electronic Structure of Anatase TiO_2 Doped with Transition Metal V,Cr,Fe and Cu Atoms by the Linearized Augmented Plane Wave Method

The electronic density of states and band structures of doped and un-doped anatase TiO2 were studied by the Linearized Augmented Plane Wave method based on the density functional theory. The calculation shows that the band structures of TiO2 crystals doped with transition metal atoms become narrower. Interesting, an excursion towards high energy level with increasing atomic number in the same element period could be observed after doping with transition metal atoms.

Reflection of Plane Waves from a Free Surface of an Initially Stressed Transversely Isotropic Dissipative Medium

The governing equations of a transversely isotropic dissipative medium are solved analytically to obtain the speeds of plane waves. The appropriate solutions satisfy the required boundary conditions at the stress-free surface to obtain the expressions of the reflection coefficients of reflected quasi-P (qP) and quasi-SV (qSV) waves in closed form for the incidence of qP and qSV waves. A particular model is chosen for numerical computation of these reflection coefficients for a certain range of the angle of incidence. The numerical values of these reflection coefficients are shown graphically against the angle of incidence for different values of initial stress parameter. The impact of initial stress parameter on the reflection coefficients is observed significantly.