With the three dimensional(3D)oblique incident waves exactly determined for the free field,the soil seismic responses in both frequency and time domains are studied by the 2.5 dimension(2.5D)finite/infinite element me...With the three dimensional(3D)oblique incident waves exactly determined for the free field,the soil seismic responses in both frequency and time domains are studied by the 2.5 dimension(2.5D)finite/infinite element method.First,the free-field responses in frequency domain are solved exactly for 3D arbitrary incident P and SV waves,which requires no coordinate conversion or extra effort for SV waves with super-critical incident angles.Next,the earthquake spectra are incorporated by the concept of equivalent seismic forces on the near-field boundary,based only on the displacements input derived for unit ground accelerations of each frequency using the 2.5D approach.For the asymmetric 2.5D finite/infinite element model adopted,the procedure for soil seismic analysis is presented.The solutions computed by the proposed method are verified against those of Wolf’s and de Barros and Luco’s and for inversely calculated ground motions.Of interest is that abrupt variation in soil response occurs around the critical angle on the wave propagation plane for SV waves.In addition,the horizontal displacements attenuate with increasing horizontal incident angle,while the longitudinal ones increase inversely for 3D incident P and SV waves.展开更多
This paper presents,for the first time,the consideration of three-dimensional(3D)oblique incident P and SV waves in calculating the 3D seismic response of a lined tunnel embedded in a half-space by the 2.5D finite/inf...This paper presents,for the first time,the consideration of three-dimensional(3D)oblique incident P and SV waves in calculating the 3D seismic response of a lined tunnel embedded in a half-space by the 2.5D finite/infinite element method(FIEM).Firstly,the applicability of the 2.5D FIEM for 3D seismic analysis is summarized.With the exact solutions obtained for the free field in the Appendix,the equivalent seismic forces are rationally computed for the near-field boundary,considering the horizontal and vertical excitations of the Chi-Chi Earthquake.By performing seismic analysis of the half space embedded with a tunnel using the 2.5D FIEM,the time-domain responses of the tunnel are obtained.The accuracy of the present solutions is verified against those of de Barros and Luco.Conclusions drawn from the parametric study include:(1)Stress concentration for the principal stress under oblique incident seismic waves occurs at the polar angles of 0(vault),90,180(inverted arch),and 270of the lining wall.(2)The vault and inverted arch are the weakest parts of the tunnel during earthquakes.(3)The accelerations of the tunnel during earthquakes can be regarded as of the rigid body type.(4)The responses of the tunnel lining caused by SV waves of an earthquake are much more critical than those by P waves.(5)For arbitrary seismic waves,the maximum longitudinal acceleration azmax is of the same order of magnitude as the maximum horizontal acceleration axmax.展开更多
The interaction of oblique incident water waves with a small bottom deformation on a porous ocean-bed is examined analytically here within the framework of linear water wave theory. The upper surface of the ocean is a...The interaction of oblique incident water waves with a small bottom deformation on a porous ocean-bed is examined analytically here within the framework of linear water wave theory. The upper surface of the ocean is assumed to be covered by an infinitely extended thin uniform elastic plate, while the lower surface is bounded by a porous bottom surface having a small deformation. By employing a simplified perturbation analysis, involving a small parameter c^（〈〈l ）, which measures the smallness of the deformation, the governing Boundary Value Problem （BVP） is reduced to a simpler BVP for the first-order correction of the potential function. This BVP is solved using a method based on Green＇s integral theorem with the introduction of suitable Green＇s function to obtain the first-order potential, and this potential function is then utilized to calculate the first-order reflection and transmission coefficients in terms of integrals involving the shape function c（x） representing the bottom deformation. Consideration of a patch of sinusoidal ripples shows that when the quotient of twice the component of the incident field wave number propagating just below the elastic plate and the ripple wave number approaches one, the theory predicts a resonant interaction between the bed and the surface below the elastic plate. Again, for small angles of incidence, the reflected wave energy is more as compared to the other angles of incidence. It is also observed that the reflected wave energy is somewhat sensitive to the changes in the flexural rigidity of the elastic plate, the porosity of the bed and the ripple wave numbers. The main advantage of the present study is that the results for the values of reflection and transmission coefficients obtained are found to satisfy the energy-balance relation almost accurately.展开更多
The seismic analysis of a viscoelastic half-space under two-dimensional(2D)oblique incident waves is carried out by the finite/infinite element method(FIEM).First,the frequency-domain exact solutions for the displacem...The seismic analysis of a viscoelastic half-space under two-dimensional(2D)oblique incident waves is carried out by the finite/infinite element method(FIEM).First,the frequency-domain exact solutions for the displacements and stresses of the free field are derived in general form for arbitrary incident P and SV waves.With the present formulation,no distinction needs to be made for SV waves with over-critical incident angles that make the reflected P waves disappear,while no critical angle exists for P waves.Next,the equivalent seismic forces of the earthquake(Taft Earthquake 1952)imposed on the near-field boundary are generated by combining the solutions for unit ground accelerations with the earthquake spectrum.Based on the asymmetric finite/infinite element model,the frequency-domain motion equations for seismic analysis are presented with the key parameters selected.The results obtained in frequency and time domain are verified against those of Wolf’s,Luco and de Barros’and for inversely computed ground motions.The parametric study indicated that distinct phase difference exists between the horizontal and vertical responses for SV waves with over-critical incident angles,but not for under-critical incident angles.Other observations were also made for the numerical results inside the text.展开更多
The Goda's method of separating the frequency spectrum of the unidirectional incident and reflected waves is improved. The proposed method can be applied to the separation of oblique incident and reflected waves a...The Goda's method of separating the frequency spectrum of the unidirectional incident and reflected waves is improved. The proposed method can be applied to the separation of oblique incident and reflected waves and the two wave gauges can be arranged in an arbitrary angle in front of a structure. When the projected distance of the two probes on the incident wave direction is the multiple ofthe half length of the incident waves, the singular problem will emerge by using the method. It is advised that when the projected distance of the two measured points on the incident wave direction is 0.05～0.45 times the wave length of peak frequency wave, good results can be obtained. The simulated resultant waves are separated by the method of numerical simulation and the separated wave spectra are basically corresponding to the target spectra input. The wave trains calculated by the separated incident and reflected wave frequency spectrum are approximated to the input wave trains and the reflected coefficient can be derived correctly. Therefore, the method proposed in this paper is reliable.展开更多
The analysis technology of Amplitude Variation with Offset(AVO)is one of the important methods for oil and gas reservoir prediction.Zoeppritz equation and its approximations are the theoretical basis of AVO analysis,w...The analysis technology of Amplitude Variation with Offset(AVO)is one of the important methods for oil and gas reservoir prediction.Zoeppritz equation and its approximations are the theoretical basis of AVO analysis,which assumes that the upper and lower media of a horizontal interface are single-phase media.Limited by this assumption,AVO analysis has limited prediction and identification accuracy for complex porous reservoirs.In view of this,the first-order approximate analytical expressions of oblique elastic wave at an interface of porous media are derived.Firstly,the incident and scattering characteristics of various waves at the interface of porous media are analyzed,and the displacement vectors generated by these elastic waves are described by exponential function.Secondly,the kinematic and dynamic boundary conditions at the interface of porous media are discussed.Thirdly,by substituting the displacement vectors of incident and scattered waves into boundary conditions,the exact analytical equation is derived.Then,considering the symmetry of scattering matrix in the equation,the exact analytical expressions of each scattered wave are obtained.Furthermore,under the assumptions of small incident angle,weak elasticity at an interface of porous media,and ignoring the second-and higherorder terms,the first-order approximate analytical expressions are derived.Establishing a model of sandstone porous media with different porosity in upper and lower media,the correctness of the approximate analytical expressions is verified,and the elastic wave response characteristics of lithology and pore fluids are analyzed.展开更多
In order to investigate the seismic performance of geosynthetic reinforced and pile supported(GRPS) embankment under seismic loads, an input method for three-dimensional oblique incidence of P wave was proposed. This ...In order to investigate the seismic performance of geosynthetic reinforced and pile supported(GRPS) embankment under seismic loads, an input method for three-dimensional oblique incidence of P wave was proposed. This method is based on the explicit finite element method while considering the viscous-spring artificial boundary(VSAB) condition. Using the proposed method, a numerical study was conducted, and the influence of oblique incidence on the seismic response of GRPS embankment under the oblique incident P waves was analyzed. The results indicate that in comparison with vertical incidence, the oblique incidence can significantly increase the displacement, velocity and acceleration of key locations in the GRPS embankment. The existence of geosynthetics can alleviate the impact of seismic load on the response of the embankment to a certain degree. Moreover, the number of reinforcement layers and modulus of geogrid also greatly influence the seismic performance of GRPS embankment.展开更多
Precision flat surfaces finished by hand-scraping are extensively used in the sliding parts of precision positioning systems or the surface plates of precision measuring instruments.These hand-scraped surfaces are emp...Precision flat surfaces finished by hand-scraping are extensively used in the sliding parts of precision positioning systems or the surface plates of precision measuring instruments.These hand-scraped surfaces are employed to improve the flatness of the machined surfaces and the motion accuracy of the instruments.The numerous depressions of micrometer-scale depth on a scraped surface are thought to prevent wringing and improve the lubrication properties by forming oil pockets.Form measurement of the micrometer-scale features of a hand-scraped surface is one of the reasonable methods for investigating its effect on the motion accuracy of the positioning systems or the lubrication properties of the sliding parts.Because a hand-scraped surface is a rough surface composed of micrometric complex features,it is difficult to conduct precision measurement of the three-dimensional surface form using conventional measurement instruments,such as vertical-incident interferometers or stylus-type profilometers.In this study,an oblique-incident interferometer based on Abramson interferometry is introduced to measure the form of a rough surface produced by hand-scraping.Because this technique can detect high-intensity reflected light on a rough surface,it allows the measurement of the form of the hand-scraped surface in noncontact condition.In order to shorten the measurement time,the form calculation based on a five-step phase shifting method was introduced.展开更多
The reflection of oblique incident waves from breakwaters with a partially-perforated front wall is investigated. The fluid domain is divided into two sub-domains and the eigenfunction expansion method is applied to e...The reflection of oblique incident waves from breakwaters with a partially-perforated front wall is investigated. The fluid domain is divided into two sub-domains and the eigenfunction expansion method is applied to expand velocity potentials in each domain. In the eigen-expansion of the velocity potential, evanescent waves are included. Numerical results of the present model are compared with experimental data. The effect of porosity, the relative chamber width, the relative water depth in the wave absorbing chamber and the water depth in front of the structure are discussed.展开更多
基金National Natural Science Foundation of China(Grant Nos.52078082,52008057)Chongqing Science and Technology Commission(Nos.cstc2021yszx-jscxX0001,2022YSZX-JSX0004CSTB).
文摘With the three dimensional(3D)oblique incident waves exactly determined for the free field,the soil seismic responses in both frequency and time domains are studied by the 2.5 dimension(2.5D)finite/infinite element method.First,the free-field responses in frequency domain are solved exactly for 3D arbitrary incident P and SV waves,which requires no coordinate conversion or extra effort for SV waves with super-critical incident angles.Next,the earthquake spectra are incorporated by the concept of equivalent seismic forces on the near-field boundary,based only on the displacements input derived for unit ground accelerations of each frequency using the 2.5D approach.For the asymmetric 2.5D finite/infinite element model adopted,the procedure for soil seismic analysis is presented.The solutions computed by the proposed method are verified against those of Wolf’s and de Barros and Luco’s and for inversely calculated ground motions.Of interest is that abrupt variation in soil response occurs around the critical angle on the wave propagation plane for SV waves.In addition,the horizontal displacements attenuate with increasing horizontal incident angle,while the longitudinal ones increase inversely for 3D incident P and SV waves.
文摘This paper presents,for the first time,the consideration of three-dimensional(3D)oblique incident P and SV waves in calculating the 3D seismic response of a lined tunnel embedded in a half-space by the 2.5D finite/infinite element method(FIEM).Firstly,the applicability of the 2.5D FIEM for 3D seismic analysis is summarized.With the exact solutions obtained for the free field in the Appendix,the equivalent seismic forces are rationally computed for the near-field boundary,considering the horizontal and vertical excitations of the Chi-Chi Earthquake.By performing seismic analysis of the half space embedded with a tunnel using the 2.5D FIEM,the time-domain responses of the tunnel are obtained.The accuracy of the present solutions is verified against those of de Barros and Luco.Conclusions drawn from the parametric study include:(1)Stress concentration for the principal stress under oblique incident seismic waves occurs at the polar angles of 0(vault),90,180(inverted arch),and 270of the lining wall.(2)The vault and inverted arch are the weakest parts of the tunnel during earthquakes.(3)The accelerations of the tunnel during earthquakes can be regarded as of the rigid body type.(4)The responses of the tunnel lining caused by SV waves of an earthquake are much more critical than those by P waves.(5)For arbitrary seismic waves,the maximum longitudinal acceleration azmax is of the same order of magnitude as the maximum horizontal acceleration axmax.
基金Partially Supported by a Research from Department of Science and Technology(DST),India under Grant No.SB/FTP/MS-003/2013
文摘The interaction of oblique incident water waves with a small bottom deformation on a porous ocean-bed is examined analytically here within the framework of linear water wave theory. The upper surface of the ocean is assumed to be covered by an infinitely extended thin uniform elastic plate, while the lower surface is bounded by a porous bottom surface having a small deformation. By employing a simplified perturbation analysis, involving a small parameter c^（〈〈l ）, which measures the smallness of the deformation, the governing Boundary Value Problem （BVP） is reduced to a simpler BVP for the first-order correction of the potential function. This BVP is solved using a method based on Green＇s integral theorem with the introduction of suitable Green＇s function to obtain the first-order potential, and this potential function is then utilized to calculate the first-order reflection and transmission coefficients in terms of integrals involving the shape function c（x） representing the bottom deformation. Consideration of a patch of sinusoidal ripples shows that when the quotient of twice the component of the incident field wave number propagating just below the elastic plate and the ripple wave number approaches one, the theory predicts a resonant interaction between the bed and the surface below the elastic plate. Again, for small angles of incidence, the reflected wave energy is more as compared to the other angles of incidence. It is also observed that the reflected wave energy is somewhat sensitive to the changes in the flexural rigidity of the elastic plate, the porosity of the bed and the ripple wave numbers. The main advantage of the present study is that the results for the values of reflection and transmission coefficients obtained are found to satisfy the energy-balance relation almost accurately.
基金sponsored by the following agencies:National Natural Science Foundation of China(Grant No.52078082)Chongqing Science and Technology Commission(No.cstc2019yszx-jcyjX0001,cstc2020yszx-jscxX0002,and cstc2021yszxjscxX0001).
文摘The seismic analysis of a viscoelastic half-space under two-dimensional(2D)oblique incident waves is carried out by the finite/infinite element method(FIEM).First,the frequency-domain exact solutions for the displacements and stresses of the free field are derived in general form for arbitrary incident P and SV waves.With the present formulation,no distinction needs to be made for SV waves with over-critical incident angles that make the reflected P waves disappear,while no critical angle exists for P waves.Next,the equivalent seismic forces of the earthquake(Taft Earthquake 1952)imposed on the near-field boundary are generated by combining the solutions for unit ground accelerations with the earthquake spectrum.Based on the asymmetric finite/infinite element model,the frequency-domain motion equations for seismic analysis are presented with the key parameters selected.The results obtained in frequency and time domain are verified against those of Wolf’s,Luco and de Barros’and for inversely computed ground motions.The parametric study indicated that distinct phase difference exists between the horizontal and vertical responses for SV waves with over-critical incident angles,but not for under-critical incident angles.Other observations were also made for the numerical results inside the text.
文摘The Goda's method of separating the frequency spectrum of the unidirectional incident and reflected waves is improved. The proposed method can be applied to the separation of oblique incident and reflected waves and the two wave gauges can be arranged in an arbitrary angle in front of a structure. When the projected distance of the two probes on the incident wave direction is the multiple ofthe half length of the incident waves, the singular problem will emerge by using the method. It is advised that when the projected distance of the two measured points on the incident wave direction is 0.05～0.45 times the wave length of peak frequency wave, good results can be obtained. The simulated resultant waves are separated by the method of numerical simulation and the separated wave spectra are basically corresponding to the target spectra input. The wave trains calculated by the separated incident and reflected wave frequency spectrum are approximated to the input wave trains and the reflected coefficient can be derived correctly. Therefore, the method proposed in this paper is reliable.
基金financially supported by the National Natural Science Foundation of China(Grant No.42104131)the Natural Science Foundation of Sichuan Province of China(Grant No.2022NSFSC1140)Open Fund(PLC20211101)of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation
文摘The analysis technology of Amplitude Variation with Offset(AVO)is one of the important methods for oil and gas reservoir prediction.Zoeppritz equation and its approximations are the theoretical basis of AVO analysis,which assumes that the upper and lower media of a horizontal interface are single-phase media.Limited by this assumption,AVO analysis has limited prediction and identification accuracy for complex porous reservoirs.In view of this,the first-order approximate analytical expressions of oblique elastic wave at an interface of porous media are derived.Firstly,the incident and scattering characteristics of various waves at the interface of porous media are analyzed,and the displacement vectors generated by these elastic waves are described by exponential function.Secondly,the kinematic and dynamic boundary conditions at the interface of porous media are discussed.Thirdly,by substituting the displacement vectors of incident and scattered waves into boundary conditions,the exact analytical equation is derived.Then,considering the symmetry of scattering matrix in the equation,the exact analytical expressions of each scattered wave are obtained.Furthermore,under the assumptions of small incident angle,weak elasticity at an interface of porous media,and ignoring the second-and higherorder terms,the first-order approximate analytical expressions are derived.Establishing a model of sandstone porous media with different porosity in upper and lower media,the correctness of the approximate analytical expressions is verified,and the elastic wave response characteristics of lithology and pore fluids are analyzed.
基金Projects(41202220,41472278,51478438)supported by the National Natural Science Foundation of ChinaProjects(20120022120003,20134101120009)supported by the Research Fund for the Doctoral Program of Higher Education,China+1 种基金Project(2652012065)supported by the Fundamental Research Funds for the Central Universities,ChinaProject(14A560014)supported by the Key Project of Education Department in Henan Province,China
文摘In order to investigate the seismic performance of geosynthetic reinforced and pile supported(GRPS) embankment under seismic loads, an input method for three-dimensional oblique incidence of P wave was proposed. This method is based on the explicit finite element method while considering the viscous-spring artificial boundary(VSAB) condition. Using the proposed method, a numerical study was conducted, and the influence of oblique incidence on the seismic response of GRPS embankment under the oblique incident P waves was analyzed. The results indicate that in comparison with vertical incidence, the oblique incidence can significantly increase the displacement, velocity and acceleration of key locations in the GRPS embankment. The existence of geosynthetics can alleviate the impact of seismic load on the response of the embankment to a certain degree. Moreover, the number of reinforcement layers and modulus of geogrid also greatly influence the seismic performance of GRPS embankment.
基金This research is supported by the JSPS,JapanJapan Society for the Promotion of Science,Japan KAKENHI Grant Number JP 20K04195 and research grant of the Mitutoyo Association for Science and Technology Grant Number R1901.
文摘Precision flat surfaces finished by hand-scraping are extensively used in the sliding parts of precision positioning systems or the surface plates of precision measuring instruments.These hand-scraped surfaces are employed to improve the flatness of the machined surfaces and the motion accuracy of the instruments.The numerous depressions of micrometer-scale depth on a scraped surface are thought to prevent wringing and improve the lubrication properties by forming oil pockets.Form measurement of the micrometer-scale features of a hand-scraped surface is one of the reasonable methods for investigating its effect on the motion accuracy of the positioning systems or the lubrication properties of the sliding parts.Because a hand-scraped surface is a rough surface composed of micrometric complex features,it is difficult to conduct precision measurement of the three-dimensional surface form using conventional measurement instruments,such as vertical-incident interferometers or stylus-type profilometers.In this study,an oblique-incident interferometer based on Abramson interferometry is introduced to measure the form of a rough surface produced by hand-scraping.Because this technique can detect high-intensity reflected light on a rough surface,it allows the measurement of the form of the hand-scraped surface in noncontact condition.In order to shorten the measurement time,the form calculation based on a five-step phase shifting method was introduced.
基金by Joint Fund of the National Natural Science Foundation of China the Hong Kong Science Research Bureau (49910161985)+1 种基金the National Natural Science Foundation of China (50025924,50179004)the Research Fund for the Development of harbor engineeri
文摘The reflection of oblique incident waves from breakwaters with a partially-perforated front wall is investigated. The fluid domain is divided into two sub-domains and the eigenfunction expansion method is applied to expand velocity potentials in each domain. In the eigen-expansion of the velocity potential, evanescent waves are included. Numerical results of the present model are compared with experimental data. The effect of porosity, the relative chamber width, the relative water depth in the wave absorbing chamber and the water depth in front of the structure are discussed.
