首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 553 毫秒
1.
Based on the up formulation of Biot equation with an assumption of zero permeability coefficient, a high-order transmitting boundary is derived for cylindrical elastic wave propagation in infinite saturated porous media. By this transmitting boundary the total stresses on the truncated boundaries of a numerical model, such as a finite element model, are replaced by a set of spring, dashpot and mass elements, with some additionally introduced auxiliary degrees of freedom. The transmitting boundaries are incorporated into the DIANA SWANDYNE II program and an unconditionally stable implicit time integration algorithm is adopted. Despite the assumption made in the derivation of the transmitting boundary, numerical examples show that it can provide highly accurate results for cylindrical elastic wave propagation problems in infinite saturated porous medium in case the up formulation is applicable. Although the direct applications of the proposed transmitting boundary to general two dimensional wave problems in infinite saturated porous media are not highly accurate, acceptable accuracy can still be achieved by placing the transmitting boundary at relatively large distance from the wave source.  相似文献   

2.
Based on the uU formulation of Biot equation and the assumption of zero permeability coefficient, a viscous-spring transmitting boundary which is frequency independent is derived to simulate the cylindrical elastic wave propagation in unbounded saturated porous media. By this viscous-spring boundary the effective stress and pore fluid pressure on the truncated boundary of the numerical model are replaced by a set of spring, dashpot and mass elements, and its simplified form is also given. A uU formulation FEA program is compiled and the proposed transmitting boundaries are incorporated therein. Numerical examples show that the proposed viscous-spring boundary and its simplified form can provide accurate results for cylindrical elastic wave propagation problems with low or intermediate values of permeability or frequency content. For general two dimensional wave propagation problems, spuriously reflected waves can be greatly suppressed and acceptable accuracy can still be achieved by placing the simplified boundary at relatively large distance from the wave source.  相似文献   

3.
In this paper, a study on the transient response of an elastic structure embedded in a homogeneous, isotropic and linearly elastic half-plane is presented. Transient dynamic and seismic forces are considered in the analysis. The numerical method employed is the coupled Finite-Element–Boundary-Element technique (FE–BE). The finite element method (FEM) is used for discretization of the near field and the boundary element method (BEM) is employed to model the semi-infinite far field. These two methods are coupled through equilibrium and compatibility conditions at the soil–structure interface. Effects of non-zero initial conditions due to the pre-dynamic loads and/or self-weight of the structure are included in the transient boundary element formulation. Hence, it is possible to analyse practical cases (such as dam–foundation systems) involving initial conditions due to the pre-seismic loads such as water pressure and self-weight of the dam. As an application of the proposed formulation, a gravity dam has been analysed and the results for different foundation stiffness are presented. The results of the analysis indicate the importance of including the foundation stiffness and thus the dam–foundation interaction.  相似文献   

4.
A seismic free field input formulation of the coupling procedure of the finite element (FE) and the scaled boundary finite-element (SBFE) is proposed to perform the unbounded soil-structure interaction analysis in time domain. Based on the substructure technique, seismic excitation of the soil-structure system is represented by the free-field motion of an elastic half-space. To reduce the computational effort, the acceleration unit-impulse response function of the unbounded soil is decomposed into two functions; linear and residual. The latter converges to zero and can be truncated as required. With the prescribed tolerance parameter, the balance between accuracy and efficiency of the procedure can be controlled. The validity of the model is verified by the scattering analysis of a hemi-spherical canyon subjected to plane harmonic P, SV and SH wave incidence. Numerical results show that the new procedure is very efficient for seismic problems within a normal range of frequency. The coupling procedure presented herein can be applied to linear and nonlinear earthquake response analysis of practical structures which are built on unbounded soil. Supproted by: the National Key Basic Research and Development Program under Grant No. 2002CB412709  相似文献   

5.
An axisymmetric infinite element and a three-dimensional infinite element are developed to solve three-dimensional elastic wave propagation problems in unbounded media. The elements are capable of transmitting Rayleigh, shear and compressional waves in the frequency domain. A scheme to integrate numerically the characteristic matrices of the elements is formulated based upon Gauss—Laguerre quadrature. Finally, the axisymmetric infinite element is used to find the compliance functions of a rigid circular plate subjected to harmonic loading on a semi-infinite medium. By using infinite elements, the size of the near field may be kept small. Consequently, the system is characterized by relatively few degrees of freedom, thus providing the analyst with an inexpensive solution.  相似文献   

6.
We present and evaluate the formulation of a reduced‐dimension (one‐dimensional) finite element for the nonlinear analysis of a vibrating disk in a two‐dimensional unbounded domain. As this problem is relevant in studies of the dynamic response of laterally loaded piles, numerous spring‐and‐dashpot representations of the disk undergoing displacement in an unbounded material domain have been developed to date: static and dynamic, linear and nonlinear. With the focus on material nonlinearity, the present simplified formulation circumvents the complications associated with nonlinear springs and dashpots. Indeed, the continuum‐based treatment described herein accounts for the interaction between the two modes of energy dissipation, due to wave propagation in the unbounded domain and loss associated with inelastic behavior. The formulation is a good compromise between the competing desires for realistic representation and efficient computation. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

7.
The performance of a coupled Finite Element—Analytic method in solving elastic wave scattering in infinite spaces is studied for plane SH-waves. The influence of several parameters on the accuracy of the results is investigated. These parameters are the number of terms in the series solution, the number of compatibility points across an imaginary circular interface (i.e. the points for which the continuity of both stresses and displacements are enforced), and the ratio of element dimension to wavelength in the direction of propagation. It was found that the method converges rapidly with the number of terms. Furthermore, both element dimension to wavelength ratio and the number of compatibility points significantly influence the accuracy of the approximate solution.  相似文献   

8.
Transient wave propagation in three-dimensional unbounded domains is studied. An efficient numerical approach is proposed, which is based on using the displacement unit-impulse response matrix representing the interaction force–displacement relationship on the near field/far field interface. Spatially, an approximation is used to reduce the computational effort associated with the large size of three-dimensional problems. It is based on subdividing the fully coupled unbounded domain into multiple subdomains. The displacement unit-impulse response matrices of all subdomains are calculated separately. The error associated with this spatial decoupling can be reduced by placing the near field/far field interface further away from the domain of interest. Detailed parameter studies have been conducted using numerical examples, in order to provide guidelines for the proposed spatially local schemes, and to demonstrate the accuracy and high efficiency of the proposed method for three-dimensional soil–structure interaction problems.  相似文献   

9.
To predict the earthquake response of saturated porous media it is essential to correctly simulate the generation, redistribution, and dissipation of excess pore water pressure during and after earthquake shaking. To this end, a reliable numerical tool requires a dynamic, fully coupled formulation for solid–fluid interaction and a versatile constitutive model. Presented in this paper is a 3D finite element framework that has been developed and utilized for this purpose. The framework employs fully coupled dynamic field equations with a upU formulation for simulation of pore fluid and solid skeleton interaction and a SANISAND constitutive model for response of solid skeleton. After a detailed verification and validation of the formulation and implementation of the developed numerical tool, it is employed in the seismic response of saturated porous media. The study includes examination of the mechanism of propagation of the earthquake-induced shear waves and liquefaction phenomenon in uniform and layered profiles of saturated sand deposits.  相似文献   

10.
A study on the seismic response of massive flexible strip-foundations embedded in layered soils and subjected to seismic excitation is presented. Emphasis is placed on the investigation of the system response with the aid of a boundary element–finite element formulation proper for the treatment of such soil–structure interaction problems. In the formulation, the boundary element method (BEM) is employed to overcome the difficulties that arise from modeling the infinite soil domain, and the finite element method (FEM) is applied to model the embedded massive flexible strip-foundation. The numerical solution for the soil–foundation system is obtained by coupling the FEM with the BEM through compatibility and equilibrium conditions at the soil–foundation and soil layer interfaces. A parametric study is conducted to investigate the effects of foundation stiffness and embedment on the seismic response.  相似文献   

11.
The scaled boundary finite-element method is a powerful tool used to analyse far-field boundary soil–structure interaction problems. In this paper, the method is extended to include Biot's coupled consolidation in order to deal with fully saturated soil as a two-phase medium. The advantages of this method are explained in this paper. The detailed formulation considers the general two-dimensional (2D) analysis case, accounting for body forces and surface tractions in both bounded and unbounded media.  相似文献   

12.
A three dimensional numerical model is presented capable of modelling the propagation and transmission of ground vibration in the vicinity of high speed railways. It is used to investigate the effect of embankment constituent material on ground borne vibration levels at various distances from the track.The model is a time domain explicit, dynamic finite element model capable of simulating non-linear excitation mechanisms. The entire model, including the wheel/rail interface is fully coupled. To account for the unbounded nature of the soil structure an absorbing boundary condition (infinite element) is placed at the truncated interfaces. To increase boundary absorption performance, the soil structure is modelled using an elongated spherical geometry.The complex geometries associated with the track components are modelled in detail thus allowing a highly realistic simulation of force transmission from vehicle to embankment. Lastly, quasi-static and dynamic excitation mechanisms of the vehicle locomotives are described using a multi-body approach which is fully coupled to the track using non-linear Hertzian contact theory.The resulting model is verified using experimental ground borne vibration data from high speed trains, gathered through field trials. It is then used to investigate the role of embankments in the transmission of vibration. It is found that soft embankments exhibit large deflections and act as a waveguide for railway vibrations which are trapped within the structure. This results in increased vibration levels both inside the embankment and in the surrounding soil. In contrast it is found that embankments formed from stiffer material reduce vibrations in the near and far fields.  相似文献   

13.
FiniteelementsimulationofsteadystateSHwavemotionZhen-PengLIAO(廖振鹏)andGuangYANG(杨光)(InstituteofEngineeringMechanics,StateSeism...  相似文献   

14.
The scaled boundary finite‐element method has been developed for the dynamic analysis of unbounded domains. In this method only the boundary is discretized resulting in a reduction of the spatial dimension by one. Like the finite‐element method no fundamental solution is required. This paper extends the scaled boundary finite‐element method to simulate the transient response of non‐homogeneous unbounded domains with the elasticity modulus and mass density varying as power functions of spatial coordinates. To reduce the number of degrees of freedom and the computational cost, the technique of reduced set of base functions is applied. The scaled boundary finite‐element equation for an unbounded domain is reformulated in generalized coordinates. The resulting acceleration unit‐impulse response matrix is obtained and assembled with the equation of motion of standard finite elements. Numerical examples of non‐homogeneous isotropic and transversely isotropic unbounded domains demonstrate the accuracy of the scaled boundary finite‐element method. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

15.
The physical background of singular memory models and in particular the Cole–Cole model is discussed. Three models of anisotropic linear viscoelasticity with frequency-dependent stiffness coefficients are considered. The models are constructed in such a way that anisotropic properties are separated from anelastic effects. Two of these models represent finite-speed wave propagation with singularities at the wavefronts (the exponential relaxation model) and without singularities at the wavefronts (the Cole–Cole model), while a third model called the fractional model is related to the constant Q with unbounded propagation speed. The Cole–Cole and fractional models belong to the class of singular memory models studied earlier because of their applications in polymer rheology, poroelasticity, poroacoustics, seismic wave propagation and other applications. Well-posedness of initial boundary value problems with mixed Dirichlet–Neumann boundary conditions is established for the three models. Regularity properties of the three models are examined.  相似文献   

16.
Two prediction models for calculating vibration from underground railways are developed: the pipe-in-pipe model and the coupled periodic finite element–boundary element (FE–BE) model.The pipe-in-pipe model is a semi-analytical three-dimensional model that accounts for the dynamic interaction between the track, the tunnel and the soil. The continuum theory of elasticity in cylindrical coordinates is used to model two concentric pipes: an inner pipe to represent the tunnel wall and an outer pipe to represent the surrounding soil. The tunnel and soil are coupled accounting for equilibrium of stresses and compatibility of displacements at the tunnel–soil interface. This method assumes that the tunnel is invariant in the longitudinal direction and the problem is formulated in the frequency–wavenumber domain using a Fourier transformation. A track, formulated as an Euler–Bernoulli beam, is then coupled to this model. Results are transformed to the space domain using the inverse Fourier transform.The coupled periodic FE–BE model is based on a subdomain formulation, where a boundary element method is used for the soil and a finite element method for the tunnel. The Craig–Bampton substructuring technique is used to efficiently incorporate the track in the tunnel. The periodicity of the tunnel is exploited using the Floquet transformation to formulate the track–tunnel–soil interaction problem in the frequency–wavenumber domain and to compute the wave field radiated into the soil.An invariant concrete tunnel, embedded in a homogeneous full space is analyzed using both approaches. The pipe-in-pipe model offers an exact solution to this problem, which is used to validate the coupled periodic FE–BE model. The free field response due to a harmonic load in the tunnel is predicted and results obtained with both models are compared. The advantages and limitations of both models are highlighted. The coupled periodic FE–BE model has a greater potential as it can account for the complex periodic geometry of the tunnel and the layering in a soil medium. The effect of coupling a floating slab to the tunnel–soil system is also studied with both models by calculating the insertion gain.  相似文献   

17.
Starting from two-dimensional wave equations and making use of Galerkin weighted residual approximations, discretized formulations for wave problems of a visco-elastic foundation have been derived. With considerations of geometrical and mechanical characteristics of a semi-infinite domain, a frequency-dependent compatible infinite element is also presented. Finally, by coupling the infinite elements with ordinary finite elements the system is used for simulation of propagating waves in a semi-infinite foundation. This model is not only suitable for simulations of complicated variations of geometrical conditions, but also for describing the unbounded behaviour of arbitrary multiple layers. Examples given indicate the model has excellent computational accuracy and feasibility for analysing the effects on foundation response due to the existence of faults or any other soft layers.  相似文献   

18.
An efficient method for modelling the propagation of elastic waves in unbounded domains is developed. It is applicable to soil–structure interaction problems involving scalar and vector waves, unbounded domains of arbitrary geometry and anisotropic soil. The scaled boundary finite element method is employed to derive a novel equation for the displacement unit-impulse response matrix on the soil–structure interface. The proposed method is based on a piecewise linear approximation of the first derivative of the displacement unit-impulse response matrix and on the introduction of an extrapolation parameter in order to improve the numerical stability. In combination, these two ideas allow for the choice of significantly larger time steps compared to conventional methods, and thus lead to increased efficiency. As the displacement unit-impulse response approaches zero, the convolution integral representing the force–displacement relationship can be truncated. After the truncation the computational effort only increases linearly with time. Thus, a considerable reduction of computational effort is achieved in a time domain analysis. Numerical examples demonstrate the accuracy and high efficiency of the new method for two-dimensional soil–structure interaction problems.  相似文献   

19.
高毅超  徐艳杰  金峰  王翔 《地球物理学报》2013,56(12):4189-4196
高阶双渐近时域透射边界能够同时模拟行波和快衰波的传播,并且能够在全频范围内迅速逼近准确解,具有优良的收敛性能和计算效率.本文将动水压力波高阶双渐近透射边界直接嵌入到近场有限元方程中,建立了大坝-库水动力相互作用的直接耦合分析模型.该模型的整体控制方程保留了近场有限元方程系数矩阵对称稀疏的优势,可以方便地利用现有的通用有限元求解器求解.基于有限元开源软件框架体系OpenSees(Open System for Earthquake Engineering Simulation),编程实现了直接耦合分析模型,并将其应用于二维重力坝、三维拱坝与库水动力相互作用分析.数值算例表明,该直接耦合分析模型具有很高的精度和计算效率.  相似文献   

20.
流体饱和多孔隙介质弹性波方程边界元解法研究   总被引:4,自引:2,他引:2       下载免费PDF全文
基于流体饱和多孔隙各向同性介质模型,本文首先推导了流体饱和多孔隙介质中弹性波传播的频率域系统动力方程及边界积分方程,然后给出了流体饱和多孔隙介质弹性波方程的基本解,最后,利用本文给出的边界元方法对流体饱和多孔隙各向同性介质中的弹性波传播进行了数值模拟.结果表明:不论是从固相位移,还是液相位移的地震合成记录都能看到明显的慢速P波,本文提出的流体饱和多孔隙介质弹性波边界元法是有效可行的.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号