To address the problem of subdividing inflexible rectangular grid models and their poor definition of velocity interfaces,we propose a complex structure triangular net for a minimum traveltime ray tracing global algor...To address the problem of subdividing inflexible rectangular grid models and their poor definition of velocity interfaces,we propose a complex structure triangular net for a minimum traveltime ray tracing global algorithm.Our procedure is：（1） Subdivide a triangle grid based on the Delaunay triangular subdivision criterion and the relationships of the points,lines,and the surfaces in the subdividing area.（2） Define the topology relationships and related concepts of triangular unit ray tracing.（3） The source point and wave arrival points at any time compose the propagating plane wave and the minimum traveltime and secondary source positions are calculated during the plane wave propagation.We adopt the hyperbolic approximation global algorithm for secondary source retrieving.（4） By minimum traveltime ray tracing,collect the path from receiver to source points with the neighborhood point＇s traveltime and the direction of the secondary source.Numerical simulation examples are given to test the algorithm.The results show that the triangular net ray tracing method demonstrates model subdivision flexibility,precise velocity discontinuity interfaces,and accurate computations.展开更多
The model of equations of micropolar generalized magneto-thermoelasticity is introduced within the context of the theory of two temperatures generalized thermoelasticity and we consider a problem of an isotropic homog...The model of equations of micropolar generalized magneto-thermoelasticity is introduced within the context of the theory of two temperatures generalized thermoelasticity and we consider a problem of an isotropic homogeneous micropolar medium taking into account the heat effects and allowing the magnetic field effects. A plane wave analysis is employed to obtain the exact formulas of the two temperatures （conductive and mechanical）, displacement components, micro-rotation components, stresses, couple stresses, induced electric current, electric field and magnetic field. Arbitrary application is chosen to enable us to get the complete solution. The considered variables are presented graphically and discussions are made for the results.展开更多
For the last few decades,the parameter estimation of electromagnetic plane waves i.e.,far field sources,impinging on antenna array geometries has attracted a lot of researchers due to their use in radar,sonar and unde...For the last few decades,the parameter estimation of electromagnetic plane waves i.e.,far field sources,impinging on antenna array geometries has attracted a lot of researchers due to their use in radar,sonar and under water acoustic environments.In this work,nature inspired heuristics based on the flower pollination algorithm(FPA)is designed for the estimation problem of amplitude and direction of arrival of far field sources impinging on uniform linear array(ULA).Using the approximation in mean squared error sense,a fitness function of the problem is developed and the strength of the FPA is utilized for optimization of the cost function representing scenarios for various number of sources non-coherent located in the far field.The worth of the proposed FPA based nature inspired computing heuristic is established through assessment studies on fitness,histograms,cumulative distribution function and box plots analysis.The other worthy perks of the proposed scheme include simplicity of concept,ease in the implementation,extendibility and wide range of applicability to solve complex optimization problems.These salient features make the proposed approach as an attractive alternative to be exploited for solving different parameter estimation problems arising in nonlinear systems,power signal modelling,image processing and fault diagnosis.展开更多
We construct the two-flux colliding plane wave solutions in higher-dimensional gravity theory with dilaton,and two complementary fluxes. Two kinds of solutions have been obtained: Bell-Szekeres (BS) type and homogeneo...We construct the two-flux colliding plane wave solutions in higher-dimensional gravity theory with dilaton,and two complementary fluxes. Two kinds of solutions have been obtained: Bell-Szekeres (BS) type and homogeneous type. After imposing the junction condition, we find that only the BS type solution is physically well-defined. Furthermore, we show that the future curvature singularity is always developed for our solutions.展开更多
Extremely large-scale array(XL-array)communications can significantly improve the transmission rate,spectral efficiency,and spatial resolution,and has great potential in next-generation mobile communication networks.A...Extremely large-scale array(XL-array)communications can significantly improve the transmission rate,spectral efficiency,and spatial resolution,and has great potential in next-generation mobile communication networks.A crucial problem in XLarray communications is to determine the boundary of applicable regions of the plane wave model(PWM)and spherical wave model(SWM).In this paper,we propose new PWM/SWM demarcations for XL-arrays from the viewpoint of channel gain and rank.Four sets of results are derived for four different array setups.First,an equi-power line is derived for a point-touniform linear array(ULA)scenario,where an inflection point is found at±π6 central incident angles.Second,an equi-power surface is derived for a point-touniform planar array(UPA)scenario,and it is proved that cos2(ϕ)cos2(φ)=12 is a dividing curve,where ϕ andφdenote the elevation and azimuth angles,respectively.Third,an accurate and explicit expression of the equi-rank surface is obtained for a ULA-to-ULA scenario.Finally,an approximated expression of the equirank surface is obtained for a ULA-to-UPA scenario.With the obtained closed-form expressions,the equirank surface for any antenna structure and any angle can be well estimated.Furthermore,the effect of scatterers is also investigated,from which some insights are drawn.展开更多
In this paper,the governing equations of linear,isotropic,homogeneous and generalized micropolar thermoelasticity are specialized in a plane.The governing equations are solved for plane harmonic wave solutions.Two sep...In this paper,the governing equations of linear,isotropic,homogeneous and generalized micropolar thermoelasticity are specialized in a plane.The governing equations are solved for plane harmonic wave solutions.Two separate velocity equations are obtained which indicate the existence of four plane waves with distinct speeds.A problem on reflection of plane waves from a thermally insulated/isothermal surface is considered with impedance boundary conditions.Appropriate potentials for incident and reflected waves are formulated which satisfy the boundary conditions at a plane surface.Relations between reflection coefficients as well as the expressions of energy ratios for various reflected waves are obtained.For illustration,the reflection coefficients and energy ratios of reflected waves are computed for relevant material parameters of an aluminium-epoxy composite.Effect of impedance parameters on all reflected waves is shown graphically at each angle of incidence.展开更多
Solutions to the equation of waves motion are derived for homogeneous and transversely isotropic media such as fiber-reinforced composites, and three dimensional slowness surfaces are shown as well. A brief discussion...Solutions to the equation of waves motion are derived for homogeneous and transversely isotropic media such as fiber-reinforced composites, and three dimensional slowness surfaces are shown as well. A brief discussion on the propagation of plane waves is given.Elastic plane waves are characterized by slowness vectors, wave vectors, polarization vectors and group velocity vectors, etc. The results obtained are presented in a coordinate-free form due to the introduction of the crystal axis' orieniation vector which specifies the anisotropy of the media. Therefore, the results are the most general and convenient for further application展开更多
Numerical stability when integrating plane waves of cubic SchrSdinger equation is thor- oughly analysed for some explicit exponential methods. We center on the following second- order methods： Strang splitting and La...Numerical stability when integrating plane waves of cubic SchrSdinger equation is thor- oughly analysed for some explicit exponential methods. We center on the following second- order methods： Strang splitting and Lawson method based on a one-parameter family of 2-stage 2nd-order explicit Runge-Kutta methods. Regions of stability are plotted and numerical results are shown which corroborate the theoretical results. Besides, a tech- nique is suggested to avoid the possible numerical instabilities which do not correspond to continuous ones.展开更多
We study the numerical identification of an unknown portion of the boundary on which either the Dirichlet or the Neumann condition is provided from the knowledge of Cauchy data on the remaining,accessible and known pa...We study the numerical identification of an unknown portion of the boundary on which either the Dirichlet or the Neumann condition is provided from the knowledge of Cauchy data on the remaining,accessible and known part of the boundary of a two-dimensional domain,for problems governed by Helmholtz-type equations.This inverse geometric problem is solved using the plane wavesmethod(PWM)in conjunction with the Tikhonov regularizationmethod.The value for the regularization parameter is chosen according toHansen’s L-curve criterion.The stability,convergence,accuracy and efficiency of the proposedmethod are investigated by considering several examples.展开更多
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 ＂impro...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.展开更多
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.展开更多
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 a...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.展开更多
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...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,展开更多
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, subj...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.展开更多
This paper presents 2.5D scattering of incident plane SH waves by a canyon in layered half-space by the indirect boundary element method (IBEM). The free field response is carried out to give the displacements and str...This paper presents 2.5D scattering of incident plane SH waves by a canyon in layered half-space by the indirect boundary element method (IBEM). The free field response is carried out to give the displacements and stresses on the line which forms boundary of the canyon. The fictitious uniform moving loads are applied to the same line to calculate the Green's functions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The displacements due to the free field and from the fictitious uniform moving loads have to be added to obtain the whole motion. The numerical results are carried out for the cases of a canyon in homogenous and in one layer over bedrock. The results show that the 2.5D wave scattering problem is essentially different from the 2D case, and there exist distinct differences between the wave amplification by a canyon in layered half-space and that in homogeneous half-space. The reasons for the distinct difference are explored, and the effects of the thickness and stiffness of the layer on the amplification are discussed.展开更多
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 ...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.展开更多
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 diffract...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.展开更多
In this paper we extend the standard Ultra Weak Variational Formulation （UWVF） of Maxwell＇s equations in an isotropic medium to the case of an anisotropic medium. We verify that the underlying theoretical framework...In this paper we extend the standard Ultra Weak Variational Formulation （UWVF） of Maxwell＇s equations in an isotropic medium to the case of an anisotropic medium. We verify that the underlying theoretical framework carries over to anisotropic media （however error estimates are not yet available） and completely describe the new scheme. We then consider TM mode scattering, show how this results in a Helmholtz equation in two dimensions with an anisotropic coefficient and demonstrate how to formulate the UWVF for it. In one special case, convergence can be proved. We then show some numerical results that suggest that the UWVF can successfully simulate wave propagation in anisotropic media.展开更多
基金the National Natural Science Foundation of China(Grant Nos.50974081,50774051, 51034003)the Foundation for the Author of National Excellent Doctoral Dissertation of PR China(Grant No.200958)+1 种基金the Program for Changjiang Scholars and Innovative Research Team in University(Grant No.IRT0843)the Open Research Fund Program of the State Key Laboratory of Coal Resources and Safe Mining(China University of Mining and Technology, Beijing)(No.2007-07)
文摘To address the problem of subdividing inflexible rectangular grid models and their poor definition of velocity interfaces,we propose a complex structure triangular net for a minimum traveltime ray tracing global algorithm.Our procedure is：（1） Subdivide a triangle grid based on the Delaunay triangular subdivision criterion and the relationships of the points,lines,and the surfaces in the subdividing area.（2） Define the topology relationships and related concepts of triangular unit ray tracing.（3） The source point and wave arrival points at any time compose the propagating plane wave and the minimum traveltime and secondary source positions are calculated during the plane wave propagation.We adopt the hyperbolic approximation global algorithm for secondary source retrieving.（4） By minimum traveltime ray tracing,collect the path from receiver to source points with the neighborhood point＇s traveltime and the direction of the secondary source.Numerical simulation examples are given to test the algorithm.The results show that the triangular net ray tracing method demonstrates model subdivision flexibility,precise velocity discontinuity interfaces,and accurate computations.
文摘The model of equations of micropolar generalized magneto-thermoelasticity is introduced within the context of the theory of two temperatures generalized thermoelasticity and we consider a problem of an isotropic homogeneous micropolar medium taking into account the heat effects and allowing the magnetic field effects. A plane wave analysis is employed to obtain the exact formulas of the two temperatures （conductive and mechanical）, displacement components, micro-rotation components, stresses, couple stresses, induced electric current, electric field and magnetic field. Arbitrary application is chosen to enable us to get the complete solution. The considered variables are presented graphically and discussions are made for the results.
基金the Deanship of Scientific Research at Majmaah University for supporting this work under Project Number No.R-2021-27.
文摘For the last few decades,the parameter estimation of electromagnetic plane waves i.e.,far field sources,impinging on antenna array geometries has attracted a lot of researchers due to their use in radar,sonar and under water acoustic environments.In this work,nature inspired heuristics based on the flower pollination algorithm(FPA)is designed for the estimation problem of amplitude and direction of arrival of far field sources impinging on uniform linear array(ULA).Using the approximation in mean squared error sense,a fitness function of the problem is developed and the strength of the FPA is utilized for optimization of the cost function representing scenarios for various number of sources non-coherent located in the far field.The worth of the proposed FPA based nature inspired computing heuristic is established through assessment studies on fitness,histograms,cumulative distribution function and box plots analysis.The other worthy perks of the proposed scheme include simplicity of concept,ease in the implementation,extendibility and wide range of applicability to solve complex optimization problems.These salient features make the proposed approach as an attractive alternative to be exploited for solving different parameter estimation problems arising in nonlinear systems,power signal modelling,image processing and fault diagnosis.
文摘We construct the two-flux colliding plane wave solutions in higher-dimensional gravity theory with dilaton,and two complementary fluxes. Two kinds of solutions have been obtained: Bell-Szekeres (BS) type and homogeneous type. After imposing the junction condition, we find that only the BS type solution is physically well-defined. Furthermore, we show that the future curvature singularity is always developed for our solutions.
基金supported in part by the National Natural Science Foundation of China (NSFC) under Grants 62271310 and 62125108in part by the Fundamental Research Funds for the Central Universities of Chinain part by the NSFC under Grant 62431014
文摘Extremely large-scale array(XL-array)communications can significantly improve the transmission rate,spectral efficiency,and spatial resolution,and has great potential in next-generation mobile communication networks.A crucial problem in XLarray communications is to determine the boundary of applicable regions of the plane wave model(PWM)and spherical wave model(SWM).In this paper,we propose new PWM/SWM demarcations for XL-arrays from the viewpoint of channel gain and rank.Four sets of results are derived for four different array setups.First,an equi-power line is derived for a point-touniform linear array(ULA)scenario,where an inflection point is found at±π6 central incident angles.Second,an equi-power surface is derived for a point-touniform planar array(UPA)scenario,and it is proved that cos2(ϕ)cos2(φ)=12 is a dividing curve,where ϕ andφdenote the elevation and azimuth angles,respectively.Third,an accurate and explicit expression of the equi-rank surface is obtained for a ULA-to-ULA scenario.Finally,an approximated expression of the equirank surface is obtained for a ULA-to-UPA scenario.With the obtained closed-form expressions,the equirank surface for any antenna structure and any angle can be well estimated.Furthermore,the effect of scatterers is also investigated,from which some insights are drawn.
文摘In this paper,the governing equations of linear,isotropic,homogeneous and generalized micropolar thermoelasticity are specialized in a plane.The governing equations are solved for plane harmonic wave solutions.Two separate velocity equations are obtained which indicate the existence of four plane waves with distinct speeds.A problem on reflection of plane waves from a thermally insulated/isothermal surface is considered with impedance boundary conditions.Appropriate potentials for incident and reflected waves are formulated which satisfy the boundary conditions at a plane surface.Relations between reflection coefficients as well as the expressions of energy ratios for various reflected waves are obtained.For illustration,the reflection coefficients and energy ratios of reflected waves are computed for relevant material parameters of an aluminium-epoxy composite.Effect of impedance parameters on all reflected waves is shown graphically at each angle of incidence.
文摘Solutions to the equation of waves motion are derived for homogeneous and transversely isotropic media such as fiber-reinforced composites, and three dimensional slowness surfaces are shown as well. A brief discussion on the propagation of plane waves is given.Elastic plane waves are characterized by slowness vectors, wave vectors, polarization vectors and group velocity vectors, etc. The results obtained are presented in a coordinate-free form due to the introduction of the crystal axis' orieniation vector which specifies the anisotropy of the media. Therefore, the results are the most general and convenient for further application
文摘Numerical stability when integrating plane waves of cubic SchrSdinger equation is thor- oughly analysed for some explicit exponential methods. We center on the following second- order methods： Strang splitting and Lawson method based on a one-parameter family of 2-stage 2nd-order explicit Runge-Kutta methods. Regions of stability are plotted and numerical results are shown which corroborate the theoretical results. Besides, a tech- nique is suggested to avoid the possible numerical instabilities which do not correspond to continuous ones.
文摘We study the numerical identification of an unknown portion of the boundary on which either the Dirichlet or the Neumann condition is provided from the knowledge of Cauchy data on the remaining,accessible and known part of the boundary of a two-dimensional domain,for problems governed by Helmholtz-type equations.This inverse geometric problem is solved using the plane wavesmethod(PWM)in conjunction with the Tikhonov regularizationmethod.The value for the regularization parameter is chosen according toHansen’s L-curve criterion.The stability,convergence,accuracy and efficiency of the proposedmethod are investigated by considering several examples.
文摘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.
基金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.
基金support from the Program for New Century Excellent Talents in University (NCET-05-0248)the Key Program for Applied Basic Research of Tianjin Municipality (07JCZDJC10100)
文摘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.
基金supported by National Natural Science Foundation of China under grant No.50978183
文摘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,
基金National Natural Science Foundation of China under Grant 51378384Key Project of Natural Science Foundation of Tianjin Municipality under Grant 12JCZDJC29000
文摘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.
基金supported by National Natural Science Foundation of China (Nos. 50908156 and 50978183)
文摘This paper presents 2.5D scattering of incident plane SH waves by a canyon in layered half-space by the indirect boundary element method (IBEM). The free field response is carried out to give the displacements and stresses on the line which forms boundary of the canyon. The fictitious uniform moving loads are applied to the same line to calculate the Green's functions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The displacements due to the free field and from the fictitious uniform moving loads have to be added to obtain the whole motion. The numerical results are carried out for the cases of a canyon in homogenous and in one layer over bedrock. The results show that the 2.5D wave scattering problem is essentially different from the 2D case, and there exist distinct differences between the wave amplification by a canyon in layered half-space and that in homogeneous half-space. The reasons for the distinct difference are explored, and the effects of the thickness and stiffness of the layer on the amplification are discussed.
基金National Natural Science Foundation of China Under Grant No.50378063Excellent Young Teacher Program of Ministry of Education
文摘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.
基金support from the Program for New Century Excellent Talents in University (NCET-05-0248)the Key Program for Applied Basic Research of Tianjin Municipality (07JCZDJC10100)
文摘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.
文摘In this paper we extend the standard Ultra Weak Variational Formulation （UWVF） of Maxwell＇s equations in an isotropic medium to the case of an anisotropic medium. We verify that the underlying theoretical framework carries over to anisotropic media （however error estimates are not yet available） and completely describe the new scheme. We then consider TM mode scattering, show how this results in a Helmholtz equation in two dimensions with an anisotropic coefficient and demonstrate how to formulate the UWVF for it. In one special case, convergence can be proved. We then show some numerical results that suggest that the UWVF can successfully simulate wave propagation in anisotropic media.
