首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 211 毫秒
1.
The contribution of the (linear) unbounded soil to the basic equation of motion of a non-linear analysis of soil-structure interaction consists of convolution integrals of the displacement-force relationship in the time domain and the history of the interaction forces. The former is calculated using the indirect boundary-element method, which is based on a weighted-residual technique and involves Green's functions. As an example of a non-linear soil-structure-interaction analysis, the partial uplift of the basemat of a structure is examined. As the convolution integrals have to be recalculated for each time step, the computational effort in this rigorous procedure is substantial. A reduction can be achieved by simplifying the Green's function by ‘concentrating’ the region of influence. Alternatively, assuming a specified wave pattern, a coupled system of springs and dashpots with frequency-independent coefficients can be used as an approximation.  相似文献   

2.
Starting from a weighted-residual formulation, the various boundary-element methods, i.e. the weighted-residual technique, the indirect boundary-element method and the direct boundary-element method, are systematically developed for the calculation of the dynamic-stiffness matrix of an embedded foundation. In all three methods, loads whose analytical response in the unbounded domain can be determined are introduced acting on the continuous soil towards the region to be excavated. In the weighted-residual technique and in the indirect boundary-element method, a weighting function is used; in the latter case, it is selected as the Green's function for the surface traction. In the direct boundary-element method, the surface traction along the structure-soil interface is interpolated. The same type of boundary matrices which have a clear physical interpretation are identified in the three formulations, each of which is illustrated with a simple static example. The indirect boundary-element method leads to the most accurate results. The guaranteed symmetry and the fact that the displacement arising from the applied loads can easily be calculated and compared to the prescribed displacement makes the indirect boundary-element method especially attractive for calculating the dynamic-stiffness matrix of the soil. Instead of calculating the dynamic-stiffness matrix of the embedded foundation with the boundary-element method, it can be determined as the difference of those of the regular free field and of the excavated part. The calculation of the former does not require the Green's function for the surface traction. The dynamic stiffness of the excavated part can be calculated by the finite-element method.  相似文献   

3.
The weighted-residual technique, the indirect boundary element method, the truncated indirect boundary element method and the direct boundary element method can be used to analyse nonlinear soil-structure interaction in the time domain. They are illustrated and compared by using the one-dimensional dynamic problem of the spherical cavity in an infinite space. For realistic time steps, all formulations lead to accurate results, but the weighted-residual technique and the truncated indirect boundary element method are much more efficient than the direct boundary element method in the time domain. Hysteretic damping leads to noncausal behaviour, which can, however, be neglected from a practical point of view.  相似文献   

4.
Green's influence functions are derived for a linearly distributed load acting on part of a layered elastic halfplane on a line which is inclined to the horizontal. Using these Green's functions as fundamental solutions in the boundary-element method, the dynamic-stiffness matrices of the unbounded soil with excavation, of the excavated part and of the free field are calculated. The indirect boundary-element method using distributed loads and no offset leads to more accurate results than the weighted-residual technique and the direct boundary-element method. At the natural frequencies of the undamped excavated part built-in along the structure-soil interface, the spring coefficients associated with the dynamic-stiffness matrices of the excavated part and of the free field will become infinite. If the dynamic-stiffness matrix of the soil with excavation is calculated as the difference of that of the free field and that of the excavated part, the difference of two large numbers will arise in the vicinity of these frequencies. A consistent discretization must then be used. In particular, the dynamic-stiffness matrix of the embedded part cannot be determined by the finite-element method in this case. A parametric study is performed for the dynamic-stiffness matrix of the free field for a rectangular foundation embedded in a halfplane and in a layer built-in at its base; the aspect ratio and the damping of the soil are varied.  相似文献   

5.
Ground vibrations induced by railway traffic at grade and in tunnels are often studied by means of two-and-half dimensional (2.5D) models that are based on a Fourier transform of the coordinate in the longitudinal direction of the track. In this paper, the need for 2.5D coupled finite element-boundary element models is demonstrated in two cases where the prediction of railway induced vibrations is considered. A recently proposed novel 2.5D methodology is used where the finite element method is combined with a boundary element method, based on a regularized boundary integral equation. In the formulation of the boundary integral equation, Green's functions of a layered elastic halfspace are used, so that no discretization of the free surface or the layer interfaces is required. In the first case, two alternative models for a ballasted track on an embankment are compared. In the first model, the ballast and the embankment are modelled as a continuum using 2.5D solid elements, whereas a simplified beam representation is adopted in the second model. The free field vibrations predicted by both models are compared to those measured during a passage of the TGVA at a site in Reugny (France). A very large difference is found for the free field response of both models that is due to the fact that the deformation of the cross section of the embankment is disregarded in the simplified representation. In the second case, the track and free field response due to a harmonic load in a tunnel embedded in a layered halfspace are considered. A simplified methodology based on the use of the full space Green's function in the tunnel–soil interaction problem is investigated. It is shown that the rigorous finite element-boundary element method is required when the distance between the tunnel and the free surface and the layer interfaces of the halfspace is small compared to the wavelength in the soil.  相似文献   

6.
The frequency-independent foundation impedances, commonly used in soil-structure dynamic interaction problems, are developed for a circular footing resting on a homogeneous halfspace. As they ignore the structure attached to the foundation, the error introduced in the structural response may be 50 per cent or more in the neighbourhood of the fundamental frequency of the soil-structure system. The present study proposes a new method developed for most dynamic soil-structure interaction problems. The key idea is to retain for the frequency-independent impedances values computed for the fundamental frequency of the soil-structure system; thus these values include the dynamic characteristics of the whole soil-structure system and lead to a satisfactory approximation of the exact solution over a wide frequency range. The method is developed here for the horizontal and rocking modes of a structure with a circular base resting on a homogeneous halfspace. Numerical applications are given for a simple linear oscillator in order to make possible a thorough parametric study. The response of some idealized building-foundation systems to harmonic excitation or to a seismic input is next examined in order to illustrate the efficiency of the proposed model.  相似文献   

7.
For a reservoir with an arbitrary shape of the upstream dam face and of the bottom including an adjacent regular part of constant depth extending to infinity, the hydrodynamic-stiffness matrix in the frequency domain for a displacement formulation is derived using the boundary-element method. The fundamental solution takes the boundary condition at the free surface into account. The analytical solution of the semi-infinite reservoir is used to improve the accuracy. To be able to transform the hydrodynamic-stiffness matrix from the frequency to the time domain, the singular part consisting of its asymptotic value of ω ∞ is split off. It consists of an imaginary linear term in ω which can be interpreted as a damper with a coefficient per unit area equal to the product of the mass density and the wave velocity. This also applies for a reservoir bottom of arbitrary shape. The remaining regular part of the stiffness matrix is transformed numerically. The corresponding interaction force-displacement relationship involves convolution integrals. This boundary-element solution agrees well with analytical results and with those of other numerical procedures based on a time-stepping method. The method is also applied to an actual earthquake acting on a reservoir with an irregular part with an inclined bottom and a regular part extending to infinity. The results of the analysis in the time domain coincide with those determined in the frequency domain.  相似文献   

8.
Few studies of wave propagation in layered saturated soils have been reported in the literature. In this paper, a general solution of the equation of wave motion in saturated soils, based on one kind of practical Biot‘s equation,was deduced by introducing wave potentials. Then exact dynamic-stiffness matrices for a poroelastic soil layer and halfspace were derived, which extended Wolf‘s theory for an elastic layered site to the case of poroelasticity, thus resolving a fundamental problem in the field of wave propagation and soil-structure interaction in a poroelastic layered soil site. By using the integral transform method, Green‘s functions of horizontal and vertical uniformly distributed loads in a poroelastic layered soil site were given. Finally, the theory was verified by numerical examples and dynamic responses by comparing three different soil sites. This study has the following advantages: all parameters in the dynamic-stiffness matrices have explicitly physical meanings and the thickness of the sub-layers does not affect the precision of the calculation which is very convenient for engineering applications. The present theory can degenerate into Wolf‘s theory and yields numerical results approaching those for an ideal elastic layered site when porosity tends to zero.  相似文献   

9.
本文研究了土-结构动力相互作用对采取不同控制措施的结构控制效果的影响。文中首先建立了主动调谐质量阻尼器(ATMD)、半主动磁流变阻尼器(MR)和被动多重调谐质量阻尼器(MTMD)等三种结构控制措施在时域中的控制算法和控制律,然后基于子结构法,采用间接边界元方法,通过傅里叶变换,推导了分别安装三种结构控制措施的受控结构在频域中的运动方程,数值仿真分析了某36层高层建筑的地震反应及其控制效果。结果表明,当采用ATMD或MTMD控制时,考虑土-结构动力相互作用后结构地震反应有所减小;当采用MR控制时,考虑土-结构动力相互作用后结构地震反应有很大程度的减小。由此看来,在设计软土地基上高层结构的结构控制措施时,不考虑土-结构动力相互作用对结构控制效果的影响是偏于安全的。  相似文献   

10.
IntroductionThe analysis of dynamic soil-structure interaction for important engineering project is still based on linear model (including equivalent linear model) with complex damping, and traditional frequency domain method (Lysmer, et al, 1975, 1981; DING, et al, 1999). Namely, first calculating frequency domain solution by Fourier transform, and then calculating time domain solution by Fourier inverse transform. The motion equation of a system in frequency domain is usually written as (…  相似文献   

11.
Seismic response of bridge pier on rigid caisson foundation in soil stratum   总被引:2,自引:0,他引:2  
An analytical method to study the seismic response of a bridge pier supported on a rigid caisson foundation embedded in a deep soil stratum underlain by a homogeneous half space is developed. The method reproduces the kinematic and inertial responses, using translational and rotational distributed Winkler springs and dashpots to simulate the soil-caisson interaction. Closed-form solutions are given in the frequency domain for vertical harmonic S-wave excitation. Comparison with results from finite element (FE) analysis and other available solutions demonstrates the reliability of the model. Results from parametric studies are given for the kinematic and inertial responses. The modification of the fundamental period and damping ratio of the bridge due to soil-structure interaction is graphically illustrated.  相似文献   

12.
An Erratum has been published for this article in Earthquake Engineering & Structural Dynamics 33(6) 2004, 793. The dynamic stiffness of a foundation embedded in a multiple‐layered halfspace is calculated postulating one‐dimensional wave propagation in cone segments. In this strength‐of‐materials approach the sectional property of the cone segment increases in the direction of wave propagation. Reflections and refractions with waves propagating in corresponding cone segments occur at layer interfaces. Compared to rigorous procedures the novel method based on cone segments is easy to apply, provides conceptual clarity and physical insight in the wave propagation mechanisms. This method postulating one‐dimensional wave propagation in cone segments with reflections and refractions at layer interfaces is evaluated, calculating the dynamic stiffness of a foundation embedded in a multiple‐layered halfspace. For sites resting on a flexible halfspace and fixed at the base, engineering accuracy (deviation of ±20%) is achieved for all degrees of freedom with a vast parameter variation. The behaviour below the cut‐off frequency in an undamped site fixed at its base is also reliably predicted. The accuracy is, in general, better than for the method based on cone frustums, which can lead to negative damping. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

13.
A time-domain formulation is proposed for the transient response analysis of general, three-dimensional structures resting on a homogeneous, elastic halfspace subjected to either external loads or seismic motions. The formulation consists of two parts: (a) the time domain formulation of the soil behaviour and (b) the coupling of the corresponding soil algorithms to the Finite Element Code ANSYS. As far as the structure is concerned, this coupling opens the way for the analysis of non-linear soil–structure interaction. The approach is based on halfspace Green's functions for displacements elicited by Heaviside time-dependent surface point loads. Hence, the spatial discretisation can be confined to the contact area between the foundation and the soil, i.e. no auxiliary grid beyond the foundation as for conventional boundary element formulations is required. The method is applied to analyse the dynamic response of a railway track due to a moving wheel set by demonstrating the influence of ‘through-the-soil coupling’.  相似文献   

14.
Vertical vibration of an embedded rigid foundation in a poroelastic soil   总被引:4,自引:0,他引:4  
This paper considers time-harmonic vertical vibration of an axisymmetric rigid foundation embedded in a homogeneous poroelastic soil. The soil domain is represented by a homogeneous poroelastic half space that is governed by Biot's theory of poroelastodynamics. The foundation is subjected to a time-harmonic vertical load and is perfectly bonded to the surrounding half space. The contact surface can be either fully permeable or impermeable. The dynamic interaction problem is solved by employing an indirect boundary integral equation method. The kernel functions of the integral equation are the influence functions corresponding to vertical and radial ring loads, and a ring fluid source applied in the interior of a homogeneous poroelastic half space. Analytical techniques are used to derive the solution for influence functions. The indirect boundary integral equation is solved by using numerical quadrature. Selected numerical results for vertical impedance of rigid foundations are presented to demonstrate the influence of poroelastic effect, foundation geometry, hydraulic boundary condition along the contact surface and frequency of excitation.  相似文献   

15.
A study on the dynamic characteristics of rigid foundations with special geometries such as square or circular with concentric internal holes, is presented. The foundations are resting on a homogeneous, linear elastic halfspace and are subjected to external forces or seismic wave excitation. Both ‘relaxed’ and ‘non-relaxed’ boundary conditions at the interface between the foundation and the halfspace are considered, and several parametric studies are conducted to assess the influence of either type of boundary conditions upon each of the possible modes of vibration. Results for massive and massless foundations are presented in time and frequency domains for impulsive and harmonic excitations, respectively. A time domain boundary element method (BEM) developed by the authors for the solution of a class of 3-D soil-structure interaction (SSI) problems is used for all the analyses reported in this work. The accuracy and efficiency of the method and the BEM models developed in this work are assessed on the basis of comparison studies with published results.  相似文献   

16.
A new model named double-shear model based on Pasternak foundation and Timoshenko beam theory is developed to evaluate the effect of a forced harmonic vibration pile to its adjacent pile in multilayered soil medium. The double-shear model takes into account the shear deformation and the rotational inertia of piles as well as the shear deformation of soil. The piles are simulated as Timoshenko beams, which are embedded in a layered Pasternak foundation. The differential equation of transverse vibration for a pile is solved by the initial parameter method. The dynamic interaction factors for the layered soil medium are obtained by the transfer matrix method. The formulation and the implementation have been verified by means of several examples. The individual shear effects of soil and piles on the interaction factors are evaluated through a parametric study. Compared to Winkler model with Euler beam, the present model gives much better results for the dynamic interaction of piles embedded in stiff soil with small slenderness ratios. Finally, the effect of a forced long pile to a short pile embedded in multilayered soil medium is studied in detail.  相似文献   

17.
To calculate the dynamic-stiffness matrix in the time domain (unit-impulse response functions) of the unbounded medium, the infinitesimal finite element cell method based solely on the finite element formulation and working exclusively in the time domain is developed. As in the cloning algorithm, the approach is based on similarity of the unbounded media corresponding to the interior and exterior boundaries of the infinitesimal finite element cell. The derivation can be performed exclusively in the time domain, or alternatively in the frequency domain. At each time station a linear system of equations is solved. The consistent-boundary method to analyse a layered medium in the frequency domain and the viscous-dashpot boundary method are special cases of the infinitesimal finite element cell method. The error is governed by the finite element discretization in the circumferential direction, as the width of the finite-element cell in the radial direction is infinitesimal. The infinitesimal finite element cell method is thus ‘exact in the finite-element sense’. This method leads to highly accurate results for a vast class of problems, ranging from a one-dimensional spherical cavity to a rectangular foundation embedded in a half-plane.  相似文献   

18.
The basic equation of motion to analyse the interaction of a non-linear structure and an irregular soil with the linear unbounded soil is formulated in the time domain. The contribution of the unbounded soil involves convolution integrals of the dynamic-stiffness coefficients in the time domain and the corresponding motions. Alternatively, a flexibility formulation for the contribution of the unbounded soil using the dynamic-flexibility coefficients in the time domain, together with the direct-stiffness method for the structure and the irregular soil can be applied. The dynamic-stiffness or flexibility coefficient in the time domain is calculated as the inverse Fourier transform of the corresponding value in the frequency domain. The dynamic-stiffness coefficient's asymptotic behaviour for high frequencies determines the singular part whose transformation exists only in the sense of a distribution. As the dynamic-flexibility coefficient converges to zero for the frequency approaching infinity, the corresponding coefficient in the time domain is simpler to calculate, as no singular part exists. The salient features of the dynamic-stiffness and flexibility coefficients in the time domain are illustrated using a semi-infinite rod with exponentially increasing area. The dynamic-flexibility coefficients in the time domain are calculated for a rigid circular disc resting on the surface of an elastic halfspace and of a layer built-in at its base. Material damping is also introduced using the three-parameter Kelvin and the Voigt models.  相似文献   

19.
An impedance matrix is derived for the relationship between displacements and external excitations of a rigid or flexible foundation embedded in a layered soil medium. The unknown contact distributed force between the foundation and soil is expanded in the frequency domain as a twofold series of azimuthal and radial components; each term represents a basic or fundamental distribution. As a result, the total response of the soil, either of displacements or stresses, has the same type of series expression except for the fundamental distributions replaced by influence functions. The coefficients of the series expansion, appearing in both equilibrium conditions of the foundation and compatibility conditions on the contact surface, relate the foundation displacements and excitations, and, therefore, result in the impedance matrix. Avoidance of integral equations in the soil-structure interaction analysis is the merit of the present approach.  相似文献   

20.
通过方位角的Fourier变换及关于径向的Hankel变换,将柱坐标系下一般形式的弹性力学基本方程和本构方程转化为波数域内非齐次的二阶常微分方程组,利用初参数法,建立了介质层传递矩阵.利用层间完全接触条件,给出了任意埋藏源层状弹性半空间频域奇异解,时域奇异解可通过关于频率的Fourier积分得到,该方法可应用到固体、流体层的情况.  相似文献   

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

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