Electromagnetic scattering from targets situated in half space is solved by applying fast inhomogeneous plane wave algorithm combined with a tabulation and interpolation method. The integral equation is set up based o...Electromagnetic scattering from targets situated in half space is solved by applying fast inhomogeneous plane wave algorithm combined with a tabulation and interpolation method. The integral equation is set up based on derivation of dyadic Green's functions in this environment. The coupling is divided into nearby region and well-separated region by grouping. The Green's function can be divided into two parts: primary term and reflected term. In the well-separated region, the two terms are both expressed as Sommerfeld integral, which can be accelerated by deforming integral path and taking interpolation and extrapolation. For the nearby region, the direct Sommerfeld integral makes the filling of impedance matrix time-expensive. A tabulation and interpolation method is applied to speed up this process. This infinite integral is pre-computed in sampling region, and a two-dimensional table is then set up. The impedance elements can then be obtained by interpolation. Numerical results demonstrate the accuracy and efficiency of this algorithm.展开更多
In this paper,we introduce an incident angle based fusion method for radar and infrared sensors to improve the recognition rate of complex targets under half space scenarios,e.g.,vehicles on the ground in this paper.F...In this paper,we introduce an incident angle based fusion method for radar and infrared sensors to improve the recognition rate of complex targets under half space scenarios,e.g.,vehicles on the ground in this paper.For radar sensors,convolutional operation is introduced into the autoencoder,a“winner-take-all(WTA)”convolutional autoencoder(CAE)is used to improve the recognition rate of the radar high resolution range profile(HRRP).Moreover,different from the free space,the HRRP in half space is more complex.In order to get closer to the real situation,the half space HRRP is simulated as the dataset.The recognition rate has a growth more than 7%com-pared with the traditional CAE or denoised sparse autoencoder(DSAE).For infrared sensor,a convolutional neural network(CNN)is used for infrared image recognition.Finally,we com-bine the two results with the Dempster-Shafer(D-S)evidence theory,and the discounting operation is introduced in the fusion to improve the recognition rate.The recognition rate after fusion has a growth more than 7%compared with a single sensor.After the discounting operation,the accuracy rate has been improved by 1.5%,which validates the effectiveness of the proposed method.展开更多
The present paper contributes in studying the phase velocities of P- and S-waves in a half space subjected to a compressive initial stress and gravity field. The density and acceleration due to gravity vary quadratica...The present paper contributes in studying the phase velocities of P- and S-waves in a half space subjected to a compressive initial stress and gravity field. The density and acceleration due to gravity vary quadratically along the depth. The dispersion equation is derived in a closed form. It is shown that the phase velocities depend not only on the initial stress, gravity, and direction of propagation but also on the inhomogeneity parameter associated with the density and acceleration due to gravity. Various particular cases are obtained, and the results match with the classical results. Numerical investigations on the phase velocities of P- and S-waves against the wave number are made for various sets of values of the material parameters, and the results are illustrated graphically. The graphical user interface model is developed to generalize the effect.展开更多
A novel cloaking scheme to hide an object in a half space from electromagnetic (EM) detection without reflection is firstly presented. The proposed cloaking scheme contains a couple of matching strips, which consist...A novel cloaking scheme to hide an object in a half space from electromagnetic (EM) detection without reflection is firstly presented. The proposed cloaking scheme contains a couple of matching strips, which consist of an isotropic material layer and an anisotropic UPML layer, located right under the bottom surface of a semi-cylindrical cloaking shell. Simple expressions for the material parameters of the cloaking scheme are derived. Numerical simulations are also performed, and a good cloaking effect is achieved. The cloaking scheme is effective to hide the local object with strong scattering characters placed on mobile carders, such as the radar antenna system on an aircraft.展开更多
The vibration analysis of a plate on an elastic foundation is an important problem in engineering. It is the interaction of a plate with the three-dimensional half space and the plate is usually loaded from both the u...The vibration analysis of a plate on an elastic foundation is an important problem in engineering. It is the interaction of a plate with the three-dimensional half space and the plate is usually loaded from both the upper and lower surfaces. The contact pressure from the soil can not be predefined. According to Lamb's solution for a single oscillating force acting on a point on the surface of an elastic half space, and the relevant approximation formulae, a relation between the local pressure and the deflection of the plate has been proposed. Based on this analysis, the reaction of the soil can be represented as the deformation of the plate. Therefore, the plate can be separated from the soil and only needs to be divided by a number of elements in the analysis. The following procedure is the same as the standard finite element method. This is a semi-analytical and semi-numerical method. It has been applied to the dynamic analysis of circular or rectangular plates on the elastic half space, at low or high frequency vibration, and on rigid, soft or flexible foundations. The results show that this method is versatile and highly accurate.展开更多
The dual integral equations of vertical forced vibration of elastic plate on an elastic half space subject to harmonic uniform distribution loading are established according to the mixed boundary-value condition. By a...The dual integral equations of vertical forced vibration of elastic plate on an elastic half space subject to harmonic uniform distribution loading are established according to the mixed boundary-value condition. By applying Abel transformation the dual integral equations are reduced to Fredholm integral equation of the second kind which is solved numerically.展开更多
The uniqueness for the solutions mentioned in the subject is proved by using the uniqueness of the solution for the internal boundary problem of Laplace and bi-Laplace equations of the first kind as well as of the sec...The uniqueness for the solutions mentioned in the subject is proved by using the uniqueness of the solution for the internal boundary problem of Laplace and bi-Laplace equations of the first kind as well as of the second.展开更多
In this paper, some thermoelastic problems in the half space are studied by using the general solutions of the elastic equations. The method presented here is extremely effective for the axisymmetric problems of the h...In this paper, some thermoelastic problems in the half space are studied by using the general solutions of the elastic equations. The method presented here is extremely effective for the axisymmetric problems of the half space as well as the half plane problems.展开更多
The precise L^(p) norm of a class of Forelli-Rudin type operators on the Siegel upper half space is given in this paper.The main result not only implies the upper L^(p) norm estimate of the Bergman projection,but also...The precise L^(p) norm of a class of Forelli-Rudin type operators on the Siegel upper half space is given in this paper.The main result not only implies the upper L^(p) norm estimate of the Bergman projection,but also implies the precise L^(p) norm of the Berezin transform.展开更多
In this paper,an analytic method is developed to address steady SH-wave scattering and perform dynamic analysis of multiple circular cavities in half space.The scattered wave function used for scattering of SH-waves b...In this paper,an analytic method is developed to address steady SH-wave scattering and perform dynamic analysis of multiple circular cavities in half space.The scattered wave function used for scattering of SH-waves by multiple circular cavities,which automatically satisfies the stress-free condition at the horizontal surface,is constructed by applying the symmetry of the SH-wave scattering and the method of multi-polar coordinates system.Applying this scattered wave function and method of moving coordinates,the original problem can be transformed to the problem of SH-wave scattering by multiple circular cavities in the full space.Finally,the solution of the problem can be reduced to a series of algebraic equations and solved numerically by truncating the infinite algebraic equations to the finite ones.Numerical examples are provided for case with two cavities to show the effect of wave number,and the distances between the centers of the cavities and from the centers to the ground surface on the dynamic stress concentration around the cavity impacted by incident steady SH-wave.展开更多
General solution of stresses solved from the two dimensiona l system of equilibrium equations in Cartesian coordinates is characterized by the presence ...General solution of stresses solved from the two dimensiona l system of equilibrium equations in Cartesian coordinates is characterized by the presence of two families of characteristic lines along which initial stresses and discontinuities in them are transmitted intact far down to infinity.This is against our intuition and not verifiable by experimental findings. For the fundamental case of infinite uniform pressure on the upper surface,a comparison between solutions from equilibrium equations in Cartesian coordinates and from those in polar coordinates is carried out in details.The semi infinite characteristic lines in the former are bent up to exponential spirals with both ends on the upper surface in the latter.Thus,the transmission pattern from solution in polar coordinates comes closer to actual situation.However,in polar reference frame,the solution for distribution of stresses in particulate half space under surface strip pressure or so can then only be obtained from boundary value problem of second order partial differential equation.展开更多
We present in this paper a numerical method for hypersingular boundary integral equations. This method was developed for planar crack problems: additional edge singularities are known to develop in that case. This pa...We present in this paper a numerical method for hypersingular boundary integral equations. This method was developed for planar crack problems: additional edge singularities are known to develop in that case. This paper includes a rigorous error analysis proving the convergence of our numerical scheme. Three types of examples are covered: the Laplace equation in free space, the linear elasticity equation in free space, and in half space.展开更多
We consider the asymptotic behavior of solutions to a model of hyperbolicelliptic coupled system on the half-line R+ = (0, ∞),with the Dirichlet boundary condition u(0, t) = 0. S. Kawashima and Y. Tanaka [Kyushu...We consider the asymptotic behavior of solutions to a model of hyperbolicelliptic coupled system on the half-line R+ = (0, ∞),with the Dirichlet boundary condition u(0, t) = 0. S. Kawashima and Y. Tanaka [Kyushu J. Math., 58(2004), 211-250] have shown that the solution to the corresponding Cauchy problem behaviors like rarefaction waves and obtained its convergence rate when u_ u+. Our main concern in this paper is the boundary effect. In the case of null-Dirichlet boundary condition on u, asymptotic behavior of the solution (u, q) is proved to be rarefaction wave as t tends to infinity. Its convergence rate is also obtained by the standard L2-energy method and Ll-estimate. It decays much lower than that of the corresponding Cauchy problem.展开更多
The expression of the equivalent stiffness of the saturated poro-elastic half space interacting with an infinite beam to harmonic moving loads is obtained via the Fourier transformation method in the frequency wave nu...The expression of the equivalent stiffness of the saturated poro-elastic half space interacting with an infinite beam to harmonic moving loads is obtained via the Fourier transformation method in the frequency wave number domain. Based on the obtained equivalent stiffness, the frequency wave number domain solution of the beam-half-space system is obtained by the compatibility condition between the beam and the half space. Critical velocity of harmonic moving loads along an infinite Euler-Bernoulli elastic beam is determined. The time domain solutions for the beam and the saturated poro-elastic half space are derived by means of the inverse Fourier transformation method. Also, the influences of the load speed, frequency and material parameters of the poro-elastic half space on the responses of the beam are investigated. Numerical results show that the frequency corresponding to the maximum deflection and bending moment increases with increasing load speed. Moreover, it can be seen that at higher frequencies, the dynamic response is independent of the load speed. The present results also show that for a beam overlying a saturated poro-elastic half space, there still exist critical velocities even when the load velocity is larger than the shear wave speed of the medium.展开更多
We study positive solutions of the following polyharmonic equation with Hardy weights associated to Navier boundary conditions on a half space:where rn is any positive integer satisfying 0 〈 2m 〈 n. We first prove ...We study positive solutions of the following polyharmonic equation with Hardy weights associated to Navier boundary conditions on a half space:where rn is any positive integer satisfying 0 〈 2m 〈 n. We first prove that the positive solutions of (0.1) are super polyharmonic, i.e.,where x* = (x1,... ,Xn-1, --Xn) is the reflection of the point x about the plane Rn-1. Then, we use the method of moving planes in integral forms to derive rotational symmetry and monotonicity for the positive solution of (0.3), in which α can be any real number between 0 and n. By some Pohozaev type identities in integral forms, we prove a Liouville type theorem--the non-existence of positive solutions for (0.1).展开更多
In this paper, we investigate the Cahn-Hilliard equation defined on the half space and subject to the Neumann boundary and initial condition. For given initial data in some sobolev space, we prove the existence and an...In this paper, we investigate the Cahn-Hilliard equation defined on the half space and subject to the Neumann boundary and initial condition. For given initial data in some sobolev space, we prove the existence and analytic smoothing effect of the solution.展开更多
Reducing the linear system of two first order equilibrium equations involving normal stress σ(ρ,θ) and shearing stress v(ρ,θ), by elimination, to two decoupled second order equations in σ and v, ...Reducing the linear system of two first order equilibrium equations involving normal stress σ(ρ,θ) and shearing stress v(ρ,θ), by elimination, to two decoupled second order equations in σ and v, we find that, for pressure only case, v(ρ,θ) vanishes in the half space. Consequently, the second order equation in σ can be simplified. In the language of linear system analysis, the medium(system) function, characterizing the mechanical behavior of a particulate medium in pressure only case, is obtained from the simplified second order equation ( 2 ρ+ 2 θ)σ(ρ,θ)=0 and can be inverted to give impulse reponse explicitly. Thus, response σ α(ρ,θ) may be computed directly from input, i.e., the surface pressure φ α(ρ) , by integration. Some explicit formulas for transmission problems, including response to input of strip linearly increasing pressure, are given in the paper.展开更多
The wave propagation behavior in an elastic wedge-shaped medium with an arbitrary shaped cylindrical canyon at its vertex has been studied.Numerical computation of the wave displacement field is carried out on and nea...The wave propagation behavior in an elastic wedge-shaped medium with an arbitrary shaped cylindrical canyon at its vertex has been studied.Numerical computation of the wave displacement field is carried out on and near the canyon surfaces using weighted-residuals(moment method).The wave displacement fields are computed by the residual method for the cases of elliptic,circular,rounded-rectangular and flat-elliptic canyons,The analysis demonstrates that the resulting surface displacement depends,as in similar previous analyses,on several factors including,but not limited,to the angle of the wedge,the geometry of the vertex,the frequencies of the incident waves,the angles of incidence,and the material properties of the media.The analysis provides intriguing results that help to explain geophysical observations regarding the amplification of seismic energy as a function of site conditions.展开更多
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 deriv...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.展开更多
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 d...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.展开更多
文摘Electromagnetic scattering from targets situated in half space is solved by applying fast inhomogeneous plane wave algorithm combined with a tabulation and interpolation method. The integral equation is set up based on derivation of dyadic Green's functions in this environment. The coupling is divided into nearby region and well-separated region by grouping. The Green's function can be divided into two parts: primary term and reflected term. In the well-separated region, the two terms are both expressed as Sommerfeld integral, which can be accelerated by deforming integral path and taking interpolation and extrapolation. For the nearby region, the direct Sommerfeld integral makes the filling of impedance matrix time-expensive. A tabulation and interpolation method is applied to speed up this process. This infinite integral is pre-computed in sampling region, and a two-dimensional table is then set up. The impedance elements can then be obtained by interpolation. Numerical results demonstrate the accuracy and efficiency of this algorithm.
基金supported by the National Natural Science Foundation of China(61571022,61971022).
文摘In this paper,we introduce an incident angle based fusion method for radar and infrared sensors to improve the recognition rate of complex targets under half space scenarios,e.g.,vehicles on the ground in this paper.For radar sensors,convolutional operation is introduced into the autoencoder,a“winner-take-all(WTA)”convolutional autoencoder(CAE)is used to improve the recognition rate of the radar high resolution range profile(HRRP).Moreover,different from the free space,the HRRP in half space is more complex.In order to get closer to the real situation,the half space HRRP is simulated as the dataset.The recognition rate has a growth more than 7%com-pared with the traditional CAE or denoised sparse autoencoder(DSAE).For infrared sensor,a convolutional neural network(CNN)is used for infrared image recognition.Finally,we com-bine the two results with the Dempster-Shafer(D-S)evidence theory,and the discounting operation is introduced in the fusion to improve the recognition rate.The recognition rate after fusion has a growth more than 7%compared with a single sensor.After the discounting operation,the accuracy rate has been improved by 1.5%,which validates the effectiveness of the proposed method.
基金supported by the Research Fellow of Indian School of Mines in Dhanbad (No. 2010DR0016)
文摘The present paper contributes in studying the phase velocities of P- and S-waves in a half space subjected to a compressive initial stress and gravity field. The density and acceleration due to gravity vary quadratically along the depth. The dispersion equation is derived in a closed form. It is shown that the phase velocities depend not only on the initial stress, gravity, and direction of propagation but also on the inhomogeneity parameter associated with the density and acceleration due to gravity. Various particular cases are obtained, and the results match with the classical results. Numerical investigations on the phase velocities of P- and S-waves against the wave number are made for various sets of values of the material parameters, and the results are illustrated graphically. The graphical user interface model is developed to generalize the effect.
基金Project supported by the Program for New Century Excellent Talents in University of China (Grant No. NCET-10-0894)
文摘A novel cloaking scheme to hide an object in a half space from electromagnetic (EM) detection without reflection is firstly presented. The proposed cloaking scheme contains a couple of matching strips, which consist of an isotropic material layer and an anisotropic UPML layer, located right under the bottom surface of a semi-cylindrical cloaking shell. Simple expressions for the material parameters of the cloaking scheme are derived. Numerical simulations are also performed, and a good cloaking effect is achieved. The cloaking scheme is effective to hide the local object with strong scattering characters placed on mobile carders, such as the radar antenna system on an aircraft.
文摘The vibration analysis of a plate on an elastic foundation is an important problem in engineering. It is the interaction of a plate with the three-dimensional half space and the plate is usually loaded from both the upper and lower surfaces. The contact pressure from the soil can not be predefined. According to Lamb's solution for a single oscillating force acting on a point on the surface of an elastic half space, and the relevant approximation formulae, a relation between the local pressure and the deflection of the plate has been proposed. Based on this analysis, the reaction of the soil can be represented as the deformation of the plate. Therefore, the plate can be separated from the soil and only needs to be divided by a number of elements in the analysis. The following procedure is the same as the standard finite element method. This is a semi-analytical and semi-numerical method. It has been applied to the dynamic analysis of circular or rectangular plates on the elastic half space, at low or high frequency vibration, and on rigid, soft or flexible foundations. The results show that this method is versatile and highly accurate.
文摘The dual integral equations of vertical forced vibration of elastic plate on an elastic half space subject to harmonic uniform distribution loading are established according to the mixed boundary-value condition. By applying Abel transformation the dual integral equations are reduced to Fredholm integral equation of the second kind which is solved numerically.
基金theResearchFoundationofEducationalCommitteeofYunnanProvince China
文摘The uniqueness for the solutions mentioned in the subject is proved by using the uniqueness of the solution for the internal boundary problem of Laplace and bi-Laplace equations of the first kind as well as of the second.
文摘In this paper, some thermoelastic problems in the half space are studied by using the general solutions of the elastic equations. The method presented here is extremely effective for the axisymmetric problems of the half space as well as the half plane problems.
基金supported by the National Natural Science Foundation of China(11801172,11771139,12071130)supported by the Natural Science Foundation of Zhejiang Province(LQ21A010002)supported by the Natural Science Foundation of Zhejiang Province(LY20A010007).
文摘The precise L^(p) norm of a class of Forelli-Rudin type operators on the Siegel upper half space is given in this paper.The main result not only implies the upper L^(p) norm estimate of the Bergman projection,but also implies the precise L^(p) norm of the Berezin transform.
文摘In this paper,an analytic method is developed to address steady SH-wave scattering and perform dynamic analysis of multiple circular cavities in half space.The scattered wave function used for scattering of SH-waves by multiple circular cavities,which automatically satisfies the stress-free condition at the horizontal surface,is constructed by applying the symmetry of the SH-wave scattering and the method of multi-polar coordinates system.Applying this scattered wave function and method of moving coordinates,the original problem can be transformed to the problem of SH-wave scattering by multiple circular cavities in the full space.Finally,the solution of the problem can be reduced to a series of algebraic equations and solved numerically by truncating the infinite algebraic equations to the finite ones.Numerical examples are provided for case with two cavities to show the effect of wave number,and the distances between the centers of the cavities and from the centers to the ground surface on the dynamic stress concentration around the cavity impacted by incident steady SH-wave.
文摘General solution of stresses solved from the two dimensiona l system of equilibrium equations in Cartesian coordinates is characterized by the presence of two families of characteristic lines along which initial stresses and discontinuities in them are transmitted intact far down to infinity.This is against our intuition and not verifiable by experimental findings. For the fundamental case of infinite uniform pressure on the upper surface,a comparison between solutions from equilibrium equations in Cartesian coordinates and from those in polar coordinates is carried out in details.The semi infinite characteristic lines in the former are bent up to exponential spirals with both ends on the upper surface in the latter.Thus,the transmission pattern from solution in polar coordinates comes closer to actual situation.However,in polar reference frame,the solution for distribution of stresses in particulate half space under surface strip pressure or so can then only be obtained from boundary value problem of second order partial differential equation.
文摘We present in this paper a numerical method for hypersingular boundary integral equations. This method was developed for planar crack problems: additional edge singularities are known to develop in that case. This paper includes a rigorous error analysis proving the convergence of our numerical scheme. Three types of examples are covered: the Laplace equation in free space, the linear elasticity equation in free space, and in half space.
基金The research was supported by the National Natural Science Foundation of China #10625105 and #10431060, the Program for New Century Excellent Talents in University #NCET-04-0745. Acknowledgement Authors would like to thank the anonymous referee for his/her helpful suggestions and comments.
文摘We consider the asymptotic behavior of solutions to a model of hyperbolicelliptic coupled system on the half-line R+ = (0, ∞),with the Dirichlet boundary condition u(0, t) = 0. S. Kawashima and Y. Tanaka [Kyushu J. Math., 58(2004), 211-250] have shown that the solution to the corresponding Cauchy problem behaviors like rarefaction waves and obtained its convergence rate when u_ u+. Our main concern in this paper is the boundary effect. In the case of null-Dirichlet boundary condition on u, asymptotic behavior of the solution (u, q) is proved to be rarefaction wave as t tends to infinity. Its convergence rate is also obtained by the standard L2-energy method and Ll-estimate. It decays much lower than that of the corresponding Cauchy problem.
基金the National Natural Science Foundatio of China (No. 50679041)the Foundation of Jiangx Educational Committee (No. GJJ09367)
文摘The expression of the equivalent stiffness of the saturated poro-elastic half space interacting with an infinite beam to harmonic moving loads is obtained via the Fourier transformation method in the frequency wave number domain. Based on the obtained equivalent stiffness, the frequency wave number domain solution of the beam-half-space system is obtained by the compatibility condition between the beam and the half space. Critical velocity of harmonic moving loads along an infinite Euler-Bernoulli elastic beam is determined. The time domain solutions for the beam and the saturated poro-elastic half space are derived by means of the inverse Fourier transformation method. Also, the influences of the load speed, frequency and material parameters of the poro-elastic half space on the responses of the beam are investigated. Numerical results show that the frequency corresponding to the maximum deflection and bending moment increases with increasing load speed. Moreover, it can be seen that at higher frequencies, the dynamic response is independent of the load speed. The present results also show that for a beam overlying a saturated poro-elastic half space, there still exist critical velocities even when the load velocity is larger than the shear wave speed of the medium.
文摘We study positive solutions of the following polyharmonic equation with Hardy weights associated to Navier boundary conditions on a half space:where rn is any positive integer satisfying 0 〈 2m 〈 n. We first prove that the positive solutions of (0.1) are super polyharmonic, i.e.,where x* = (x1,... ,Xn-1, --Xn) is the reflection of the point x about the plane Rn-1. Then, we use the method of moving planes in integral forms to derive rotational symmetry and monotonicity for the positive solution of (0.3), in which α can be any real number between 0 and n. By some Pohozaev type identities in integral forms, we prove a Liouville type theorem--the non-existence of positive solutions for (0.1).
文摘In this paper, we investigate the Cahn-Hilliard equation defined on the half space and subject to the Neumann boundary and initial condition. For given initial data in some sobolev space, we prove the existence and analytic smoothing effect of the solution.
文摘Reducing the linear system of two first order equilibrium equations involving normal stress σ(ρ,θ) and shearing stress v(ρ,θ), by elimination, to two decoupled second order equations in σ and v, we find that, for pressure only case, v(ρ,θ) vanishes in the half space. Consequently, the second order equation in σ can be simplified. In the language of linear system analysis, the medium(system) function, characterizing the mechanical behavior of a particulate medium in pressure only case, is obtained from the simplified second order equation ( 2 ρ+ 2 θ)σ(ρ,θ)=0 and can be inverted to give impulse reponse explicitly. Thus, response σ α(ρ,θ) may be computed directly from input, i.e., the surface pressure φ α(ρ) , by integration. Some explicit formulas for transmission problems, including response to input of strip linearly increasing pressure, are given in the paper.
文摘The wave propagation behavior in an elastic wedge-shaped medium with an arbitrary shaped cylindrical canyon at its vertex has been studied.Numerical computation of the wave displacement field is carried out on and near the canyon surfaces using weighted-residuals(moment method).The wave displacement fields are computed by the residual method for the cases of elliptic,circular,rounded-rectangular and flat-elliptic canyons,The analysis demonstrates that the resulting surface displacement depends,as in similar previous analyses,on several factors including,but not limited,to the angle of the wedge,the geometry of the vertex,the frequencies of the incident waves,the angles of incidence,and the material properties of the media.The analysis provides intriguing results that help to explain geophysical observations regarding the amplification of seismic energy as a function of site conditions.
基金Program for New Century Excellent Talents in University Under Grant No. NCET-05-0248the Key Program for Applied Basic Research of Tianjin Municipality Under Grant No. 07JCZDJC10100
文摘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.
基金National Natural Science Foundation of China Under Grant No.50908156 and 50978183Tianjin Natural Science Foundation Under Grant No. 07JCZDJC10100
文摘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.