共查询到20条相似文献，搜索用时 359 毫秒
1.
The various boundary-element methods, well established in the frequency domain, are developed in the time domain for a foundation embedded in a layered halfspace. They are the weighted-residual technique and the indirect boundary-element method, based on a weighted-residual equation, and the direct boundary-element method based on a reciprocity equation, both equations involving time and space. In the indirect approach, formulating the weighted-residual equation over the last time step only results in the truncated indirect boundary-element formulation which requires a reduced computational effort. In all cases, convolution integrals occur. The truncated indirect boundary-element method leads to a highly reliable algorithm, as is verified when a linear analysis in the time domain is compared to the corresponding one in the frequency domain. This boundary-element formulation, which is non-local in space and time, represents a rigorous generally applicable method taking into account a layered halfspace in a non-linear soil-structure interaction analysis. As an example, the non-linear soil-structure interaction analysis of a structure embedded in a halfspace with partial uplift of the basemat and separation of the side wall is investigated. 相似文献
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.
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. 相似文献
4.
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. 相似文献
5.
The scattering of plane SH-waves by topographic features in a layered transversely isotropic (TI) half-space is investigated by using an indirect boundary element method (IBEM). Firstly, the anti-plane dynamic stiffness matrix of the layered TI half-space is established and the free fields are solved by using the direct stiffness method. Then, Green’s functions are derived for uniformly distributed loads acting on an inclined line in a layered TI half-space and the scattered fields are constructed with the deduced Green’s functions. Finally, the free fields are added to the scattered ones to obtain the global dynamic responses. The method is verified by comparing results with the published isotropic ones. Both the steady-state and transient dynamic responses are evaluated and discussed. Numerical results in the frequency domain show that surface motions for the TI media can be significantly different from those for the isotropic case, which are strongly dependent on the anisotropy property, incident angle and incident frequency. Results in the time domain show that the material anisotropy has important effects on the maximum duration and maximum amplitudes of the time histories. 相似文献
6.
Elastic wave propagation in an irregularly layered medium 总被引：1，自引：0，他引：1
Rossana Vai Jos Manuel Castillo-Covarrubias Francisco J. Snchez-Sesma Dimitri Komatitsch Jean-Pierre Vilotte 《Soil Dynamics and Earthquake Engineering》1999,18(1):385
The indirect boundary element method (IBEM) is used to simulate wave propagation in two-dimensional irregularly layered elastic media for internal line sources. The method is based on the integral representation for scattered elastic waves using single layer boundary sources. Fulfillment of the boundary conditions leads to a system of integral equations. Results are obtained in the frequency domain and seismograins are computed through Fourier synthesis. In order to test and validate the method we present various comparisons between our results and the time series obtained analytically for a buried line source in a half-space and by using the recently developed spectral element method (SEM). 相似文献
7.
《地震工程与工程振动(英文版)》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. 相似文献
8.
模拟地震波传播的三维逐元并行谱元法 总被引：2，自引：0，他引：2
高效地震波场正演模拟对于复杂模型中地震波传播与成像研究至关重要.本文在谱元法原理框架内,对已有逐元谱元法改进,提出一种新的逐元并行谱元法求解三维地震波运动方程,并得到地震波场.逐元并行谱元法的核心思想在于在单元上进行质量矩阵与解向量的乘积运算,并将此运算平均分配至每一个CPU计算核心,此处理有利提升谱元法的并行计算效率.同时,根据Gauss-Lobatto-Legendre(GLL)数值积分点与插值点重合的特点,将稠密单元刚度矩阵的存储转化成单元雅克比矩阵行列式的值及其逆的存储,大幅减少谱元法计算内存开销.此外,在模型边界上利用逐元并行谱元法求解二阶位移形式完美匹配层(PML)吸收边界条件,消除边界截断而引入的虚假反射.通过逐元并行谱元法得到的数值解与解析解对比,以及实际地震波场模拟,数值结果证实了逐元并行谱元法用于地震波场模拟的高效性. 相似文献
9.
本用震源力学理论和方法研究了徐淮地区从１９７０年以来构造应力场的方向和强度的时空变化过程。结果表明：以唐山地震为分界线，本区的应力场Ｐ轴取向由震前平均６１．８°变为震后平均７７．７°。如果将本区以宿北断裂为界分为南区和北区两个部分，则北区的Ｐ轴取向从６８．１°变７１．２°，而南区的Ｐ轴取向由６２．５°变到８３．６°，南区的变化明显于北区。 相似文献
10.
<正>The objective of this paper is to develop a dynamic slip model for a shear crack under constant stress drop.This crack problem is formulated by a traction boundary integral equation(BIE) in the frequency domain and then solved by the hyper-singular boundary element method as well as the regularization technique proposed in this paper.Based on the spectral integral form of the kernel function,the unbounded term can be isolated and extracted from the hyper-singular kernel function by using the method of subtracted and added back in wave number domain.Finally,based on the inverse transformation from the frequency domain to the time domain,the time histories of crack opening displacement under constant stress drop can be determined.Three rupture models(simultaneous rupture model,symmetric bilaterally-propagating model and unilaterally propagating model) with specified time histories of stress drop are considered in this paper.Even though these three models will cause the same final slip shapes because of the same constant stress drop,the associated slip time functions differ significantly from each other during the rupture process. 相似文献
11.
A technique for modeling transient wave propagation in unbounded media is extended and applied to seismic soil–structure interaction analysis in the time domain. The technique, based on the discontinuous Galerkin method, requires lower computational cost and less storage than the boundary element method, and the time‐stepping scheme resulting from Newmark's method in conjunction with the technique is unconditionally stable, allowing for efficient and robust time‐domain computations. To extend the technique to cases characterized by seismic excitation, the free‐field motion is used to compute effective forces, which are introduced on the boundary of the computational domain containing the structure and the soil in the vicinity of the structure. A numerical example on a dam–foundation system subjected to seismic excitation demonstrates the performance of the method. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
12.
We propose a domain decomposition method based on the spectral element method(DDM-SEM)for elastic wave computation in frequency domain.It combines the high accuracy of the spectral element method and the high degree of parallelism of a domain decomposition technique,which makes this method suitable for accurate and efficient simulations of large scale problems in elastodynamics.In the DDM-SEM,the original large-scale problem is divided into a number of well designed subdomains.We use the spectral element method independently for each subdomain,and the neighboring subdomains are connected by a frequency-domain version of Riemann transmission condition(RTC)for elastic waves.For the proposed method,we can employ the non-conforming meshes and different interpolation orders in different subdomains to maximize the efficiency.By separating the internal and boundary unknowns of each subdomain,an efficient and naturally parallelizable block LDU direct solver is developed to solve the final system matrix.Numerical experiments verify its accuracy and efficiency,and show that the proposed DDM-SEM can be a promising numerical tool for accurately and effectively solving large and multi-scale problems of elastic waves.It is potentially valuable for the frequency domain seismic inversion where multiple source illuminations are required. 相似文献
13.
2.5D direct‐current resistivity forward modelling and inversion by finite‐element–infinite‐element coupled method 
下载免费PDF全文
下载免费PDF全文 Yuan Yuan Jianke Qiang Jingtian Tang Zhengyong Ren Xiao Xiao 《Geophysical Prospecting》2016,64(3):767-779
To reduce the numerical errors arising from the improper enforcement of the artificial boundary conditions on the distant surface that encloses the underground part of the subsurface, we present a finite‐element–infinite‐element coupled method to significantly reduce the computation time and memory cost in the 2.5D direct‐current resistivity inversion. We first present the boundary value problem of the secondary potential. Then, a new type of infinite element is analysed and applied to replace the conventionally used mixed boundary condition on the distant boundary. In the internal domain, a standard finite‐element method is used to derive the final system of linear equations. With a novel shape function for infinite elements at the subsurface boundary, the final system matrix is sparse, symmetric, and independent of source electrodes. Through lower upper decomposition, the multi‐pole potentials can be swiftly obtained by simple back‐substitutions. We embed the newly developed forward solution to the inversion procedure. To compute the sensitivity matrix, we adopt the efficient adjoint equation approach to further reduce the computation cost. Finally, several synthetic examples are tested to show the efficiency of inversion. 相似文献
14.
Amplification of in-plane seismic ground motion by group cavities in layered half-space (Ⅰ) 下载免费PDF全文
下载免费PDF全文 Amplification of in-plane seismic ground motion by underground group cavities in layered half-space is studied both in frequency domain and time domain by using indirect boundary element method (IBEM), and the effect of cavity interval and spectrum of incident waves on the amplification are studied by numerical examples. It is shown that there may be large interaction between cavities, and group cavities with certain intervals may have significant amplification to seismic ground motion. The amplification of PGA (peak ground acceleration) and its PRS (peak response spectrum) can be increased up to 45.2% and 84.4%, for an example site in Tianjin, under the excitation of Taft wave and El Centro wave; and group cavities may also affect the spectra of the seismic ground motion. It is suggested that the effect of underground group cavities on design seismic ground motion should be considered. 相似文献
15.
The impact of the incident angle of earthquake motion on the seismic response of the long lined tunnels is studied. Based on the time‐domain finite element method with the viscous‐spring artificial boundary condition, the earthquake motion of oblique incidence is transformed into the equivalent nodal forces acting on the truncated boundary of finite element model. In the present work, the formulas of equivalent nodal forces for the plane P wave with arbitrary incident angle are deduced and implemented into the commercial software abaqus 1 . The effectiveness of the formulas and its implementation are demonstrated by two numerical examples with the reference solutions. The proposed method is applied to investigate the seismic responses of the long lined tunnels under the obliquely incident P waves. The numerical results indicate that the seismic responses of the long lined tunnels are highly affected by the incident angles of P waves. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
16.
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. 相似文献
17.
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. 相似文献
18.
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 相似文献
19.
目前,对于可控源电磁法各向异性介质2.5维问题,主要采用一次场、二次场分离的方法消除场源奇异性并降低截断边界对计算区域的影响.该方法数值计算精度高,但是很难适用于复杂地形条件下的数值模拟.针对复杂地形问题,基于总场的有限元方法表现出一定的优越性,然而,这种方法存在场源奇异性问题和截断边界问题.本文采用基于总场计算的方法对带地形的可控源电磁法2.5维各向异性介质进行模拟研究,推导了考虑电导率和介电常数各向异性的2.5维控制方程;引入网格加密-收缩算法降低场源奇异性的影响范围,提升数值计算效率;引入行波分解边界条件降低截断边界的影响;提出任意采样反傅里叶变换方法,快速、高精度地计算出空间域电磁场分量.理论模型数值算例中:首先,验证了本文算法的有效性;其次,对任意各向异性倾角产生的可控源电磁响应规律进行研究;最后,采用山丘模型对各向异性介质电磁场的响应规律进行了模拟和分析. 相似文献
20.
平面SV波在层状半空间中沉积谷地周围的散射 总被引：2，自引：0，他引：2
采用间接边界元方法求解了入射平面SV波在层状半空间中沉积谷地周围的散射问题.问题的解答包含自由场和散射场两部分.自由场可由直接刚度法求得,散射场由层状半空间中斜线荷载动力格林函数来模拟.文中以入射平面SV波在基岩上单一土层中沉积谷地周围散射为例研究了土层和沉积谷地周围的影响.结果表明,由于考虑了土层的动力特性,平面SV... 相似文献