共查询到20条相似文献，搜索用时 15 毫秒
1.
A. BurmanParsuram Nayak P. AgrawalDamodar Maity 《Soil Dynamics and Earthquake Engineering》2012,34(1):62-68
A time domain transient analysis of a concrete gravity dam and its foundation has been carried out in a coupled manner using finite element technique and the effect of Soil-Structure Interaction (SSI) has been incorporated using a simplified direct method. A two dimensional plane strain dam-foundation model has been used for the time history analysis to compute the stresses and displacements against earthquake loading considering the effect of soil-structure interaction. An effective boundary condition has been implemented by attaching dashpots to the vertical boundaries. The material damping effects have also been considered and the dam and foundation have both been modeled as linear, elastic materials. To achieve a greater degree of accuracy, the displacements and stresses calculated in the free-field analysis have also been added to those developed in the complete dam-foundation analysis. The proposed algorithm has been simulated for the case of two published problems and in both the cases the results have been found to be in close agreement. The proposed technique is quite simple and easy to implement in the computer code. The outcomes of the results show the efficacy of the developed method. 相似文献
2.
3.
为减少直接分析三维大尺度复杂土-结构动力相互作用问题的计算量,提高计算效率,本文直接从波动方程出发,提出了较常规子结构法更简单的两步简化计算过程,即第一步简化上部复杂结构体系为集中质量杆系模型,并求基础处等效输入,第二步通过等效输入求上部结构各位置的动力反应.其中第一步计算主要采用集中质量显式有限单元法结合局部透射人工... 相似文献
4.
Edward H. Stehmeyer III Dimitris C. Rizos 《Soil Dynamics and Earthquake Engineering》2006,26(5):421-434
This work presents a simplified, yet accurate model of rigid foundation-soil systems for the dynamic analysis of structures including SSI effects. The simplified model is based on closed form solutions that reproduce the characteristic B-spline impulse response functions (BIRF) of 3D continuous soil-foundation systems, as obtained from rigorous boundary element method (BEM) analysis. The proposed simplified model is used within the framework of B-spline impulse response techniques and can be coupled directly to other solution techniques, such as the finite element method (FEM). Validation and application studies demonstrate the accuracy and versatility of the simplified model for the direct time domain solution of dynamic SSI problems involving rigid square surface foundations of any size. The proposed model, although simplified, demonstrates similar high accuracy to that of more rigorous solutions based on domain discretizations. 相似文献
5.
A local time-domain transmitting boundary for simulating cylindrical elastic wave propagation in infinite media 总被引：3，自引：0，他引：3
Finite element simulation of the time-dependent wave propagation in infinite media requires enforcing the transmitting boundary to replace the truncated far-field infinite domain so as to model the effect of the wave radiation towards infinity. This paper proposed a novel local time-domain transmitting boundary for simulating the cylindrical elastic wave radiation problem. This boundary is a mechanical model consisting of the spring, dashpot and mass elements, with the auxiliary degrees of freedom introduced, which is dynamically stable and easily implemented into the commercial finite element codes. Numerical analysis of the cylindrical elastic wave radiation problem indicates that the proposed transmitting boundaries with the order N=3 for cylindrical P and SV waves and with the order N=4 for cylindrical SH wave have very high accuracy, even when the artificial boundary at wave source. The proposed transmitting boundary with order N=0 can be applied approximately to the general two-dimensional infinite elastic wave problems that contain the more complex outgoing wave fields at artificial boundary than the cylindrical waves. The plane-strain Lamb problem is analyzed with the acceptable engineering accuracy achieved. On the other hand, the proposed transmitting boundary with higher order can be a tool to localize the temporal convolution that appears in an exact time-domain transmitting boundary for the general infinite wave problems. This potential applicability is mentioned. 相似文献
6.
We present a new finite element (FE) method for magnetotelluric modelling of three-dimensional conductivity structures. Maxwell's equations are treated as a system of first-order partial differential equations for the secondary fields. Absorbing boundary conditions are introduced, minimizing undesired boundary effects and allowing the use of small computational domains. The numerical algorithm presented here is an iterative, domain decomposition procedure employing a nonconforming FE space. It does not use global matrices, therefore allowing the modellization of large and complicated structures. The algorithm is naturally parallellizable, and we show results obtained in the IBM SP2 parallel supercomputer at Purdue University. The accuracy of the numerical method is verified by checking the computed solutions with the results of COMMEMI, the international project on the comparison of modelling methods for electromagnetic induction. 相似文献
7.
本用震源力学理论和方法研究了徐淮地区从１９７０年以来构造应力场的方向和强度的时空变化过程。结果表明：以唐山地震为分界线，本区的应力场Ｐ轴取向由震前平均６１．８°变为震后平均７７．７°。如果将本区以宿北断裂为界分为南区和北区两个部分，则北区的Ｐ轴取向从６８．１°变７１．２°，而南区的Ｐ轴取向由６２．５°变到８３．６°，南区的变化明显于北区。 相似文献
8.
《地震工程与工程振动(英文版)》2021,20(3)
The direct finite element method is a type commonly used for nonlinear seismic soil-structure interaction(SSI) analysis. This method introduces a truncated boundary referred to as an artificial boundary meant to divide the soilstructure system into finite and infinite domains. An artificial boundary condition is used on a truncated boundary to achieve seismic input and simulate the wave radiation effect of infinite domain. When the soil layer is particularly thick, especially for a three-dimensional problem, the computational efficiency of seismic SSI analysis is very low due to the large size of the finite element model, which contains an whole thick soil layer. In this paper, an accurate and efficient scheme is developed to solve the nonlinear seismic SSI problem regarding thick soil layers. The process consists of nonlinear site response and SSI analysis. The nonlinear site response analysis is still performed for the whole thick soil layer. The artificial boundary at the bottom of the SSI analysis model is subsequently relocated upward from the bottom of the soil layer(bedrock surface) to the location nearest to the structure as possible. Finally, three types of typical sites and underground structures are adopted with seismic SSI analysis to evaluate the accuracy and efficiency of the proposed efficient analysis scheme. 相似文献
9.
A procedure which involves a non‐linear eigenvalue problem and is based on the substructure method is proposed for the free‐vibration analysis of a soil–structure system. In this procedure, the structure is modelled by the standard finite element method, while the unbounded soil is modelled by the scaled boundary finite element method. The fundamental frequency, and the corresponding radiation damping ratio as well as the modal shape are obtained by using inverse iteration. The free vibration of a dam–foundation system, a hemispherical cavity and a hemispherical deposit are analysed in detail. The numerical results are compared with available results and are also verified by the Fourier transform of the impulsive response calculated in the time domain by the three‐dimensional soil–structure–wave interaction analysis procedure proposed in our previous paper. The fundamental frequency obtained by the present procedure is very close to that obtained by Touhei and Ohmachi, but the damping ratio and the imaginary part of modal shape are significantly different due to the different definition of damping ratio. This study shows that although the classical mode‐superposition method is not applicable to a soil–structure system due to the frequency dependence of the radiation damping, it is still of interest in earthquake engineering to evaluate the fundamental frequency and the corresponding radiation damping ratio of the soil–structure system. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
10.
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. 相似文献
11.
This paper presents a time-dependent semi-analytical artificial boundary for numerically simulating elastic wave propagation problems in a two-dimensional homogeneous half space. A polygonal boundary is considered in the half space to truncate the semi-infinite domain, with an appropriate boundary condition imposed. Using the concept of the scaled boundary finite element method, the wave equation of the truncated semi-infinite domain is represented by the partial differential equation of non-constant coefficients. The resulting partial differential equation has only one spatial coordinate variable and time variable. Through introducing a few auxiliary functions at the truncated boundary, the resulting partial differential equations are further transformed into linear time-dependent equations. This allows an artificial boundary to be derived from the time-dependent equations. The proposed artificial boundary is local in time, global at the truncated boundary and semi-analytical in the finite element sense. Compared with the scaled boundary finite element method, the main advantage in using the proposed artificial boundary is that the requirement for solving a matrix form of Lyapunov equation to obtain the unit-impulse response matrix is avoided, so that computer efforts are significantly reduced. The related numerical results from some typical examples have demonstrated that the proposed artificial boundary is of high accuracy in dealing with time-dependent elastic wave propagation in two-dimensional homogeneous semi-infinite domains. 相似文献
12.
A direct time domain boundary element method is presented based on the Stokes fundamental solutions, discretized in both time and space, and an efficient time step-by-step solution that minimizes the accumulation of errors. A non-singular numerical integration procedure, in the Cauchy sense, is proposed for the generation of the associated influence matrices. This methodology is shown to be efficient for the solution of a number of computationally intensive problems in the area of soil–structure interaction. In addition, an algorithm for the direct calculation of the response of massive foundations to externally applied forces and/or obliquely incident seismic waves is introduced. The accuracy and computational efficiency of the proposed methodologies is established through a number of comparison studies. 相似文献
13.
Using reciprocal theorems for dynamic and static boundary value problems, boundary integral equations are presented for wave propagation in elastic, isotropic media and compressible, inviscid fluids in the time domain as well as in the frequency domain. For the analysis of fluid–soil and fluid–structure systems, suitable coupling conditions are prescribed along the interfaces. The numerical treatment of the boundary integral equations consists of a point collocation and of a discretization of the boundary, in which constant and linear approximation functions are assumed. Step-by-step integration is applied to the time-dependent equations, where again the states are taken to be linear and constant over each time interval. These boundary element procedures are used to analyse the response of dams due to horizontal and vertical ground motions considering dam–water interaction and absorption of hydrodynamic pressure waves at the reservoir bottom or at the far end into the soil medium. Both the frequency response and the impulse generated transient response are investigated. 相似文献
14.
《Soil Dynamics and Earthquake Engineering》2003,23(4):249-253
A parallel soil–structure interaction (SSI) model is presented for applications on distributed computer systems. Substructring method is applied to the SSI system and a coupled finite–infinite element based parallel computer program is developed. In the SSI system, infinite elements are used to represent the soil which extends to infinity. In this case, a large finite element mesh is required to define the near field for reliable predictions. The resulting large-scale problems are solved on distributed computer systems in this study. The domain is represented by separated substructures and an interface. The number of substructures are determined by the available processors in the parallel platform. To avoid the formation of large interface equations, smaller interface equations are distributed to processors while substructure contributions are performed. This saves a lot of memory storage and computational effort. Direct solution techniques are used for the solution of interface and substructure equation systems. The program is investigated through some example problems. The example problems exposed the need for solving large-scale problems in order to reach better results. The results of the example problems demonstrated the benefits of the parallel SSI algorithm. 相似文献
15.
从Maxwell方程组出发,推导了探地雷达(GPR)有限元波动方程.阐述了透射边界条件和Sarma边界条件的原理,推导了这两种边界条件的理论公式;通过在衰减层内加入过渡带优化了Sarma边界条件的加载方法,压制了介质区和衰减层交界面处的人为反射.考虑到透射边界条件与Sarma边界条件不同的理论机制,提出了一种结合透射边界条件和Sarma边界条件的混合边界条件,它利用Sarma边界条件对到达边界区域的GPR波能量衰减功能和透射边界对GPR波能量的透射功能,使GPR波经过Sarma边界条件的衰减吸收后,再通过透射边界条件将剩余能量透射出去,集成了二者的优势.并以二维均匀模型中的中心脉冲激励源方式为例,通过Matlab程序实现,以GPR的全波场快照的直观方式,对比了有、无边界条件及不同边界条件对人工截断边界的处理效果,说明了该混合边界条件对到达截断边界处的GPR波的处理优于单一边界条件.最后,以基于混合边界条件的有限单元法对两个典型的GPR地电模型进行了正演模拟,指导了GPR数据处理与工程实践. 相似文献
16.
The study of wave propagation in finite/infinite media has many applications in geotechnical and structural earthquake engineering and has been a focus of research for the past few decades. This paper presents an analysis of 2D anti- plane problems （Love waves） and 2D in-plane problems （Rayleigh waves） in the frequency domain in media consisting of a near-field irregular and a far-field regular part. The near field part may contain structures and its boundaries with the far-field can be of any shape. In this study, the irregular boundaries of the near-field are treated as consistent boundaries, extending the concept of Lysmer＇s vertical consistent boundaries. The presented technique is called the Condensed Hyperelements Method （CHM）. In this method, the irregular boundary is limited to a vertical boundary at each end that is a consistent boundary at the far-field side. Between the two ends, the medium is discretized with hyperelements. Using static condensation, the stiffness matrix of the far-field is derived for the nodes on the irregular boundary. Examples of the application of the CHM illustrate its excellent accuracy and efficiency. 相似文献
17.
S. H. Ju 《地震工程与结构动力学》2003,32(9):1431-1442
This paper discusses how to use the three‐dimensional (3D) time‐domain finite‐element method incorporating the least‐squares method to calculate the equivalent foundation mass, damping and stiffness matrices. Numerical simulations indicate that the accuracy of these equivalent matrices is acceptable when the applied harmonic force of 1+sine is used. Moreover, the accuracy of the least‐squares method using the 1+sine force is not sensitive to the first time step for inclusion of data. Since the finite‐element method can model problems flexibly, the equivalent mass, damping and stiffness matrices of very complicated soil profiles and foundations can be established without difficulty using this least‐squares method. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
18.
3D viscous-spring artificial boundary in time domain 总被引：7，自引：1，他引：7
After a brief review of studies on artificial boundaries in dynamic soil-structure interaction, a three-dimensional viscous-spring artificial boundary （VSAB） in the time domain is developed in this paper. First, the 3D VSAB equations in the normal and tangential directions are derived based on the elastic wave motion theory. Secondly, a numerical simulation technique of wave motion equations along with the VSAB condition in the time domain is studied. Finally, numerical examples of some classical elastic wave motion problems are presented and the results are compared with the associated theoretical solutions, demonstrating that high precision and adequate stability can be achieved by using the proposed 3D VSAB. The proposed 3D VSAB can be conveniently incorporated in the general finite element program, which is commonly used to study dynamic soil-structure interaction problems. 相似文献
19.
The boundary element method has been successfully applied in the past to the analysis of hydrodynamic forces in two- and three-dimensional finite water reservoirs subjected to seismic ground motions. In extending the method to an infinite reservoir, the loss of energy due to pressure waves moving away towards infinity must be taken into account. In addition, for both finite and infinite reservoirs, energy is lost owing to partial absorption of the waves incident on a flexible bottom consisting of alluvial deposits. This paper presents the results of more recent research on the application of the boundary element method to the analysis of 2D reservoir vibration. Two different formulations are used: a constant boundary element formulation and a linear boundary element formulation. Special boundary conditions to treat infinite radiation and foundation damping have been incorporated in both formulations. Numerical results have been obtained for each of the two alternative formulations and compared against each other as well as with classical solutions and results obtained by other researchers. 相似文献
20.
Simple formulas are derived for the dynamic stiffness of pile group foundations subjected to horizontal and rocking dynamic loads. The formulations are based on the construction of a general model of impedance matrices as the condensation of matrices of mass, damping, and stiffness, and on the identification of the values of these matrices on an extensive database of numerical experiments computed using coupled finite element–boundary element models. The formulations obtained can be readily used for the design of both floating piles on homogeneous half‐space and end‐bearing piles and are applicable for a wide range of mechanical and geometrical parameters of the soil and piles, in particular for large pile groups. For the seismic design of a building, the use of the simple formulas rather than a full computational model is shown to induce little error on the evaluation of the response spectra and time histories. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献