共查询到20条相似文献，搜索用时 15 毫秒
1.
Dionyssia-Pinelopi N. Kontoni Dimitri E. Beskos 《Soil Dynamics and Earthquake Engineering》1987,6(4):227-238
In this paper the applicability of an approximate Boundary Element Method to uniform half-plane elastodynamic problems is investigated. This method employs the concept of images to construct approximate fundamental solutions for the half-plane and does not require any half-plane surface discretization. The method is formulated in the frequency domain for the case of harmonic disturbances or the Laplace transform domain for the case of transient disturbances. In the latter case a numerical inversion of the transformed solution is necessary to obtain the time domain response. The proposed method can be used as an alternative to boundary element methods that either utilize the infinite plane fundamental solution and thus require a half-plane surface discretization, or employ the exact half-plane fundamental solution, which even though leads to no surface discretization, is of a very lengthy and complicated form. Two characteristics numerical examples are used to illustrate the proposed method and study its advantages and disadvantages. 相似文献
2.
A boundary element formulation having discontinuous curved quadratic elements is presented for 2D elastodynamics. The first fundamental solution for static case is subtracted from and added to the first fundamental solution for dynamic case. As both kernels have the same order of singularity, the integral involving the regular expression arising from the subtraction can be calculated. matrix is calculated by employing the well-known rigid-body motion technique. The formulation is performed in Fourier transform space. Based on the formulation presented in this study, a general purpose computer program is developed for elastic or visco-elastic 2D elastodynamic problems. The program performs the analysis in Fourier transform space and can also be used for static analysis by assigning a very small value close to zero for the frequency. The results of some elastodynamic and dynamic soil–structure interaction problems obtained using the present study are compared with those in the literature. 相似文献
3.
By using integral transform methods, the Green's functions ofhorizontal harmonic force applied at the interior of the saturated half-space soil are obtained in the paper. The general solutions of the Biot dynamic equations in frequency domain are established through the use of Hankel integral transforms technique. Utilizing the above- mentioned general solutions, and the boundary conditions of the surface of the half-space and the continuous con-ditions at the plane of the horizontal force, the solutions of the boundary value problem can be determined. By the numerical inverse Hankel transforms method, the Green's functions of the harmonic horizontal force are obtainable. The degenerate case of the results deduced from this paper agrees well with the known results. Two numerical examples are given in the paper. 相似文献
4.
By using integral transform methods, the Green’s functions of horizontal harmonic force applied at the interior of the saturated
half-space soil are obtained in the paper. The general solutions of the Biot dynamic equations in frequency domain are established
through the use of Hankel integral transforms technique. Utilizing the above-mentioned general solutions, and the boundary
conditions of the surface of the half-space and the continuous conditions at the plane of the horizontal force, the solutions
of the boundary value problem can be determined. By the numerical inverse Hankel transforms method, the Green’s functions
of the harmonic horizontal force are obtainable. The degenerate case of the results deduced from this paper agrees well with
the known results. Two numerical examples are given in the paper.
Foundation item: State Natural Science Foundation (59879012) and Doctoral Foundation from State Education Commission (98024832). 相似文献
5.
IntroductionThe wave propagation problems in saturated soil are very important for the civil engineering, geophysics and seismology. Biot (1956,1962) established the theory of wave propagation in saturated soil firstly, and hereafter many researchers have used Biot theory to study wave propagation problems in saturated soil. By using integral transform and potential function method, Philippacopoulos (1988) studied the Lamb(s problem of a vertical point force applied to the surface of saturate… 相似文献
6.
The transient dynamic response of saturated soil under suddenly applied normal and horizontal concentrated loading is studied in this paper. The behavior of saturated soil is governed by Biot's consolidation theory. The general solutions for Biot equations of equilibrium are derived in terms of displacements and variations of fluid volume, using Laplace–Hankel integral transforms. The solutions in the time domain can be evaluated by numerical inverse Laplace–Hankel transforms. Selected numerical results for displacements, stresses, and pore pressures are presented. Comparisons with existing closed-form solutions for the elastic half-space are made to confirm the accuracy of the present solutions. The solutions can be used to study a variety of transient wave propagation problems and dynamical interactions between saturated soil and structures. 相似文献
7.
根据积分变换方法得出了半空间内部作用简谐水平力时的Gree函数．首先，利用Hankel积分变换方法，直接对频域内的Biot波动方程进行求解，得出Biot波动方程的通解；利用通解和半空间内部作用水平力时边界上的边界条件，以及力作用面上的连续性条件，可以得出上述边值问题的解；对于边值问题在变换域内的解进行相应的逆变换，就可以得出频域内的Gree函数．本文得到的线弹性退化解与文献中的结果吻合．最后，文中给出了两个算例． 相似文献
8.
In this paper the coupled equations governing the dynamic behavior of unsaturated soils are derived based on the poromechanics theory within the framework of the suction-based mathematical model presented by Gatmiri (1997) [Gatmiri B. Analysis of fully coupled behavior of unsaturated porous medium under stress, suction and temperature gradient. Final report of CERMES-EDF, 1997] and Gatmiri et al. (1998) [Gatmiri B, Delage P, Cerrolaza M, UDAM: a powerful finite element software for the analysis of unsaturated porous media. Adv Eng Software 1998; 29(1): 29–43]. In this formulation, the solid skeleton displacements, water pressure and air pressure are presumed to be independent variables. The Boundary Integral formulations as well as fundamental solutions for such a dynamic u–pw–pa theory are presented in this paper for the first time. The boundary integral equations are derived via the use of the weighted residuals method in a way that permits an easy discretization and implementation in a Boundary Element code. Also, the associated two dimensional (2D) fundamental solutions for such deformable porous medium with linear elastic behavior are derived in Laplace transform domain using the method of Hörmander. Finally, some numerical results are presented to show the accuracy of the proposed solutions. The derived results are verified analytically by comparison with the previously introduced corresponding fundamental solutions in elastodynamic limiting case. 相似文献
9.
Starting with the governing equations of motion and the constitutive equations of transversely isotropic elastic body, and based on the corresponding algebraic operations and the Hankel transform, the analytical layer-elements of a finite layer and a half-space are obtained in the transformed domain. According to the continuity conditions between adjacent layers, the global stiffness matrix equation is obtained by assembling the analytical layer-element of each single layer. The solutions in the transformed domain are acquired by introducing the boundary conditions into the global stiffness matrix equation, and thus, the corresponding solutions in frequency domain are achieved by taking the inversion of Hankel transform. Finally, some numerical examples are given to illustrate the accuracy of the proposed method, and to study the influence of properties and the frequency of excitation on the dynamic response of the medium. 相似文献
10.
The Laplace domain solutions have been obtained for three-dimensional groundwater flow to a well in confined and unconfined wedge-shaped aquifers. The solutions take into account partial penetration effects, instantaneous drainage or delayed yield, vertical anisotropy and the water table boundary condition. As a basis, the Laplace domain solutions for drawdown created by a point source in uniform, anisotropic confined and unconfined wedge-shaped aquifers are first derived. Then, by the principle of superposition the point source solutions are extended to the cases of partially and fully penetrating wells. Unlike the previous solution for the confined aquifer that contains improper integrals arising from the Hankel transform [Yeh HD, Chang YC. New analytical solutions for groundwater flow in wedge-shaped aquifers with various topographic boundary conditions. Adv Water Resour 2006;26:471–80], numerical evaluation of our solution is relatively easy using well known numerical Laplace inversion methods. The effects of wedge angle, pumping well location and observation point location on drawdown and the effects of partial penetration, screen location and delay index on the wedge boundary hydraulic gradient in unconfined aquifers have also been investigated. The results are presented in the form of dimensionless drawdown-time and boundary gradient-time type curves. The curves are useful for parameter identification, calculation of stream depletion rates and the assessment of water budgets in river basins. 相似文献
11.
D. PATELLA 《Geophysical Prospecting》1980,28(5):744-749
The Hankel transform theorem can be applied to the inversion of gravity data for the buried sphere, the horizontal cylinder, and the vertical rod. This new approach leads to exact solutions of the transforms for the assumed bodies. A comparison with the classical procedure by Fourier transform reveals that for the sphere and the vertical rod, the Hankel transform is preferable. 相似文献
12.
A transversely isotropic material in the sense of Green is considered. Using a series of potential functions proposed in [Eskandari-Ghadi M. A complete solution of the wave equations for transversely isotropic media. J Elasticity 2005; 81:1–19], the solutions of the transient wave equations within a half-space under surface load are obtained in the Laplace–Hankel domain for axisymmetric problems. The solutions are investigated in detail in the special case of a surface point force pulse varying with time as Heaviside function. Using Cagniard–De Hoop method, the inverse Laplace transform and inverse Hankel transform of the solutions are then obtained in the form of integrals with finite limits. For validity of the analytical results, the final formulations for surface waves are degenerated for an isotropic material and compared with the existing formulation obtained by Pekeris [The seismic surface pulse. Proc Natl Acad Sci USA 1955;41:469–80], to show that they are exactly the same. The numerical evaluations of the integrals for some transversely isotropic materials as well as an isotropic one are obtained. The present approach is then numerically verified by comparing a particular case of displacements for the surface of an isotropic half-space subjected to a point load of Heaviside function with the solutions obtained by Pekeris [The seismic surface pulse. Proc Natl Acad Sci USA 1955;41:469–80]. In addition, the wave equations for the mentioned medium are obtained on the vertical line directly under the applied surface load. The final formulations are degenerated for an isotropic material and compared with the existing formulation given in Graff [Wave motion in elastic solids. New York: Dover Publications Inc; 1975 [New Ed edition, November 1991]], to show that they are also exactly the same. Then equations are presented in graphical forms using an appropriate numerical evaluation. 相似文献
13.
The numerical evaluation of certain integral transforms is required for the interpretation of some geophysical exploration data. Digital linear filter operators are widely used for carrying out such numerical integration. It is known that the method of Wiener–Hopf minimization of the error can be used to design very efficient, short digital linear filter operators for this purpose. We have found that, with appropriate modifications, this method can also be used to design longer filters. Two filters for the Hankel J0 transform (61-point and 120-point operators), and two for the Hankel J1 transform (47-point and 140-point operators) have been designed. For these transforms, the new filters give much lower errors compared to all other known filters of comparable, or somewhat longer, size. The new filter operators and some results of comparative performance tests with known integral transforms are presented. These filters would find widespread application in many numerical evaluation problems in geophysics. 相似文献
14.
A. V. Vostroukhov S. N. Verichev A. W. M. Kok C. Esveld 《Soil Dynamics and Earthquake Engineering》2004,24(6):449-459
This paper outlines a method of analysis of the steady-state dynamic response of a stratified soil to a horizontal time-harmonic loading applied at a circular area. It is assumed the load to be uniformly distributed over the contact area and embedded at an arbitrary depth. It is shown that by means of application of integral Hankel transform, the problem can be reduced to a form of a single integral, which can be taken numerically. Results of the numerical simulations for an embedded horizontal impulse load are presented to demonstrate the efficiency of the developed procedure. 相似文献
15.
An analytical approach is developed to study the dynamic response of a flexible plate on single-layered saturated soil. The analysis is based on Biot's two-phased theory of poroelasticity and also on the classical thin-plate theory. First, the governing differential equations for saturated soil are solved by the use of Hankel transform. The general solutions of the skeleton displacements, stresses, and pore pressures, derived in the transformed domain, are subsequently incorporated into the imposed boundary conditions, which leads to a set of dual integral equations describing the corresponding mixed boundary value problem. These governing integral equations are finally reduced to the Fredholm integral equations of the second kind and solved by standard numerical procedures. The accuracy of the present solution is validated via comparisons with existing solutions for an ideal elastic half-space. Furthermore, some numerical results are presented to show the influences of the layer depth, the plate flexibility, and the soil porosity on the dynamic compliances. 相似文献
16.
Franois Chapel 《地震工程与结构动力学》1987,15(7):815-829
An extension of the boundary element method to heterogeneous domains composed of horizontal layers is here proposed. It includes a numerical computation of the corresponding Green's functions, thanks to an inverse Hankel transform of the closed form solutions obtained in the spectral domain with suitable variables derived from displacements and stress vectors to obtain the decoupling between P–SV and SH waves. Transmission and reflection operators are introduced to avoid the problem of overflowing exponentials met with in Thomson–Haskell matrices. Applications are given in the soil–structure interaction field to compute the impedances of surface and embedded circular foundations resting on a viscoelastic halfspace. 相似文献
17.
An exact stiffness matrix method is presented to evaluate the dynamic response of a multi-layered poroelastic medium due to time-harmonic loads and fluid sources applied in the interior of the layered medium. The system under consideration consists of N layers of different properties and thickness overlying a homogeneous half-plane or a rigid base. Fourier integral transform is used with respect to the x-co-ordinate and the formulation is presented in the frequency domain. Fourier transforms of average displacements of the solid matrix and pore pressure at layer interfaces are considered as the basic unknowns. Exact stiffness (impedance) matrices describing the relationship between generalized displacement and force vectors of a layer of finite thickness and a half-plane are derived explicitly in the Fourier-frequency space by using rigorous analytical solutions for Biot's elastodynamic theory for porous media. The global stiffness matrix and the force vector of a layered system is assembled by considering the continuity of tractions and fluid flow at layer interfaces. The numerical solution of the global equation system for discrete values of Fourier transform parameter together with the application of numerical quadrature to evaluate inverse Fourier transform integrals yield the solutions for poroelastic fields. Numerical results for displacements and stresses of a few layered systems and vertical impedance of a rigid strip bonded to layered poroelastic media are presented. The advantages of the present method when compared to existing approximate stiffness methods and other methods based on the determination of layer arbitrary coefficients are discussed. 相似文献
18.
In this article, we have derived a new and more general formulation of Love waves in arbitrarily irregular multi-layered media by using the global generalized reflection/transmission (abbreviated to R/T thereafter) matrices method developed earlier by Chen [17~20]. From the basic principle that the modal solutions are the non-trivial solutions of the free elastodynamic equation under appropriate boundary conditions, we naturally derived the characteristic frequencies and the corresponding distorted modes of Love wave in irregular multi-layered media. Moreover, we have derived the corresponding excitation formulation of Love waves in such laterally heterogeneous media by using the general solution of elastodynamic equation [17~20]. Similar to the result for laterally homogeneous layered structure, the Love waves radiated from a point source in irregular multi-layered media can be expressed as a superposition of distorted modes. Since the structure model used here is quite arbitrary, it can be used for so 相似文献
19.
By virtue of the precise integration method (PIM) and the technique of mixed variable formulations, solutions for the dynamic response of the multi-layered transversely isotropic medium subjected to the axisymmetric time-harmonic forces are presented. The planes of cross anisotropy are assumed to be parallel to the horizontal surface of the stratified media. Four kinds of vertically acting axisymmetric loads are prescribed either at the external surface or in the interior of the soil system. Thicknesses and number of the medium strata are not limited. Employing the Hankel integral transform in cylindrical coordinate, the axisymmetric governing equations in terms of displacements of the multi-layered media are uncoupled. Applying mixed variable formulations, more concise first-order ordinary differential matrix equations from the uncoupled motion equations can be obtained. Solutions of the ordinary differential matrix equations in the transformed domain are acquired by utilizing the approach of PIM. Since PIM is highly accurate to solve the sets of first-order ordinary differential equations, any desired accuracy of the solutions can be achieved. All calculations are based on the corresponding algebraic operations and computational efforts can be reduced to a great extent. Comparisons with the existing numerical solutions are made to confirm the accuracy of the present solutions proposed by this procedure. Several examples are illustrated to explore the influences of the type and degree of material anisotropy, the frequency of excitation and loading positions on the dynamic response of the stratified medium. 相似文献
20.
采用刚度矩阵方法结合Hankel积分变换,求解了层状黏弹性半空间中球面SH、P和SV波的自由波场.首先,在柱坐标系下建立层状黏弹性半空间的反轴对称(柱面SH波)和轴对称(柱面P-SV波)情况精确动力刚度矩阵.进而由Hankel变换将空间域内的球面波展开为波数域内柱面波的叠加,然后将球面波源所在层的上下端面固定,求得固定层内的动力响应和固定端面反力,将固端反力反向施加到层状黏弹性半空间上,采用直接刚度法求得固端反力的动力响应,叠加固定层内和固端反力动力响应,求得波数域内球面波源动力响应.最后由Hankel积分逆变换求得频率-空间域内球面波源自由场,时域结果由傅里叶逆变换求得.文中验证了方法的正确性,并以均匀半空间和基岩上单一土层中球面SH、P和SV波为例分别在频域和时域内进行了数值计算分析.研究表明基岩上单一土层中球面波自由场与均匀半空间情况有着本质差异;基岩上单一土层中球面波位移频谱峰值频率与场地固有频率相对应,基岩面的存在使得基岩上单一土层地表点的位移时程非常复杂,振动持续时间明显增长;阻尼的增大显著降低了动力响应的峰值,同时也显著减少了波在土层的往复次数. 相似文献