共查询到20条相似文献，搜索用时 171 毫秒
1.
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. 相似文献
2.
The scaled boundary finite‐element method is extended to simulate time‐harmonic responses of non‐homogeneous unbounded domains with the elasticity modulus and mass density varying as power functions of spatial coordinates. The unbounded domains and the elasticity matrices are transformed to the scaled boundary coordinates. The scaled boundary finite‐element equation in displacement amplitudes are derived directly from the governing equations of elastodynamics. To enforce the radiation condition at infinity, an asymptotic expansion of the dynamic‐stiffness matrix for high frequency is developed. The dynamic‐stiffness matrix at lower frequency is obtained by numerical integration of ordinary differential equations. Only the boundary is discretized yielding a reduction of the spatial dimension by one. No fundamental solution is required. Material anisotropy is modelled without additional efforts. Examples of two‐ and three‐dimensional non‐homogeneous isotropic and transversely isotropic unbounded domains are presented. The results demonstrate the accuracy and simplicity of the scaled boundary finite‐element method. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
3.
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. 相似文献
4.
比例边界有限元方法是一种半解析的数值计算方法,具有降维、网格灵活、严格模拟无限域和无需基本解等特点。比例边界有限元方法的基本理论是在整体坐标与局部坐标的比例边界转换基础之上建立的,相似中心的选取是否合理对分析计算具有重要的影响,导致在模拟拱坝这种不规则的空间壳体结构时,具有一定的局限性。采用子结构方法,将坝体分为若干满足相似性要求的区域可解决上述问题,以某拱坝为例给出了合理的坝体子结构分区形式,验证了子结构方法的精确性,为建立基于比例边界有限元方法的坝体-库水-地基系统的计算模型奠定了基础。 相似文献
5.
采用显式有限元方法求解开放系统动力响应的主要环节和问题 总被引：3，自引：1，他引：2
分析了采用显式有限元方法求解开放系统动力响应的主要环节中影响精度和稳定性的因素及其相互关系,指出数值解的精度和稳定性是这些因素共同作用的结果。归纳了合理输入数字荷载应满足的条件以及系统离散化准则。基于有限元方法推导了吸收边界的一般表达式,并据此提出吸收边界和透射边界之间存在确定的函数关系,并不是两类独立人工边界的新见解。认为内域的阻尼稳定性和人工边界的近场稳定性是大型复杂开放系统动力响应数值方法研究和应用中的关键科学问题,并从应用的角度提出了解决问题的思路和看法。 相似文献
6.
时域高阶双渐近透射边界能够同时模拟层状介质中行波和快衰波的传播,具有很高的计算精度和计算效率.本文将高阶双渐近透射边界推广应用到多层层状地基系统弹性波传播问题的模拟,采用广义特征值分解分析该透射边界的数值稳定性,通过移谱法消除导致数值不稳定的虚假模态.将高阶双渐近透射边界以超单元的形式直接嵌入到近场有限元方程,建立了有限元-高阶双渐近透射边界时域耦合分析模型,并将其应用到重力坝-层状地基动力相互作用分析.数值算例分析结果表明,该时域耦合分析模型具有很高的精度和计算效率,适用于实际重力坝工程的地震响应分析. 相似文献
7.
Local transmitting boundaries for transient elastic analysis 总被引：1，自引：0，他引：1
The aim of this paper is to investigate and develop alternative methods of analyzing problems in dynamic soil–structure-interaction (SSI). The interaction means that the amplitude of structural response is effected by additional energy dissipation through radiation and material damping in the soil. The surrounding soft soil behaves as a natural damper for a massive and stiff structure supported or embedded in it. The main focus is the major difficulty posed by such an analysis — the phenomena of waves that radiate outward from the excited structures towards infinity. In numerical calculations only a finite region of the foundation medium is analyzed and something is done to prevent the outgoing radiation waves from reflecting at the boundary region.Development of a simple and efficient finite element (FE) procedure for the solution directly in the time domain of transient SSI problems is the main concern. The central feature of the procedure is local absorbing boundaries used to render the computational domain finite. These boundaries are local in both time and space and are completely defined by a pair of symmetric stiffness and damping matrices. As the effort for implementing them is the same as for the impedance boundary condition (BC) considering the angle of incidence, standard assembly procedure can be used. Due to the local nature they also preserve the overall structure of the global equations of motion. Even though the focus is in the time domain the same equations of motions can be used to determine the solution under time-harmonic excitation directly in the frequency domain. Explicit formulae for the element matrices are included in the paper and numerical examples for transient radiation model problems to illustrate the validity and accuracy of the new procedures, are given. 相似文献
8.
A new numerical procedure is proposed for the analysis of three-dimensional dynamic soil–structure interaction in the time domain. In this study, the soil is modelled as a linear elastic solid, however, the methods developed can be adapted to include the effects of soil non-linearities and hysteretic damping in the soil. A substructure method, in which the unbounded soil is modelled by the scaled boundary finite-element method, is used and the structure is modelled by 8–21 variable-number-node three-dimensional isoparametric or subparametric hexahedral curvilinear elements. Approximations in both time and space, which lead to efficient schemes for calculation of the acceleration unit-impulse response matrix, are proposed for the scaled boundary finite-element method resulting in significant reduction in computational effort with little loss of accuracy. The approximations also lead to a very efficient scheme for evaluation of convolution integrals in the calculation of soil–structure interaction forces. The approximations proposed in this paper are also applicable to the boundary element method. These approximations result in an improvement over current methods. A three-dimensional Dynamic Soil–Structure Interaction Analysis program (DSSIA-3D) is developed, and seismic excitations (S-waves, P-waves, and surface waves) and externally applied transient loadings can be considered in analysis. The computer program developed can be used in the analysis of three-dimensional dynamic soil–structure interaction as well as in the analysis of wave scattering and diffraction by three-dimensional surface irregularities. The scattering and diffraction of seismic waves (P-, S-, and Rayleigh waves) by various three-dimensional surface irregularities are studied in detail, and the numerical results obtained are in good agreement with those given by other authors. Numerical studies show that the new procedure is suitable and very efficient for problems which involve low frequencies of interest for earthquake engineering. Copyright © 1999 John Wiley & Sons Ltd 相似文献
9.
高阶双渐近时域透射边界能够同时模拟行波和快衰波的传播,并且能够在全频范围内迅速逼近准确解,具有优良的收敛性能和计算效率.本文将动水压力波高阶双渐近透射边界直接嵌入到近场有限元方程中,建立了大坝-库水动力相互作用的直接耦合分析模型.该模型的整体控制方程保留了近场有限元方程系数矩阵对称稀疏的优势,可以方便地利用现有的通用有限元求解器求解.基于有限元开源软件框架体系OpenSees（Open System for Earthquake Engineering Simulation）,编程实现了直接耦合分析模型,并将其应用于二维重力坝、三维拱坝与库水动力相互作用分析.数值算例表明,该直接耦合分析模型具有很高的精度和计算效率. 相似文献
10.
Kyung Hwan Cho Moon Kyum Kim Yun Mook Lim Seong Yong Cho 《Soil Dynamics and Earthquake Engineering》2004,24(11):839-852
The seismic response analysis of a base-isolated liquid storage tank on a half-space was examined using a coupling method that combines the finite elements and boundary elements. The coupled dynamic system that considers the base isolation system and soil–structure interaction effect is formulated in time domain to evaluate accurately the seismic response of a liquid storage tank. Finite elements for a structure and boundary elements for liquid are coupled using equilibrium and compatibility conditions. The base isolation system is modeled using the biaxial hysteretic element. The homogeneous half-space is idealized using the simple spring-dashpot model with frequency-independent coefficients. Some numerical examples are presented to demonstrate accuracy and applicability of the developed method.Consequently, a general numerical algorithm that can analyze the dynamic response of base-isolated liquid storage tanks on homogeneous half-space is developed in three-dimensional coordinates and dynamic response analysis is performed in time domain. 相似文献
11.
考虑土-结构相互作用和岩土参数不确定性的核电厂结构地震响应分析 总被引：2，自引：0，他引：2
针对核电厂结构,在考虑土-结构相互作用(SSI)的情况下进行随机地震反应分析,探讨地基岩土参数的不确定性对反应堆厂房楼层反应谱(FRS)的影响。运用ANSYS软件模块建立核电厂(NPP)结构有限元模型,通过设置边界弹簧单元和阻尼装置来考虑SSI效应;并且通过设置具有概率意义的弹簧刚度和阻尼系数,来模拟土特性参数的不确定性。随机响应分析与确定性分析的结果对比,揭示了岩性地基条件下SSI效应对核电厂FRS的影响以及地基岩土参数不确定性对FRS的影响程度。研究表明,在岩性地基条件下,亦不应忽略SSI效应;考虑SSI效应的随机分析模型同确定性模型相比,二者的分析结果较为接近,两方法都可用于NPP的FRS敏感性分析评估之中,并可进行相互比照。 相似文献
12.
考虑土-结构相互作用的高层建筑抗震分析 总被引：17，自引：0，他引：17
本文采用通用有限元程序ANSYS,针对上海地区一例土-箱基-高层建筑结构进行了三维有限元分析,计算中土体的本构模型采用等效线性模型,利用粘一弹性人工边界作为土体的侧向边界,并研究了土体边界位置、土性、基础埋深、基础形式以及上部结构刚度等参数对动力相互作用体系动力特性及地震反应的影响。 相似文献
13.
This paper explores dynamic soil–bridge interaction in high speed railway lines. The analysis was conducted using a general and fully three-dimensional multi-body finite element–boundary element model formulated in the time domain to predict vibrations caused by trains passing over the bridge. The vehicle was modelled as a multi-body system, the track and the bridge were modelled using finite elements and the soil was considered as a half-space by the boundary element method. The dynamic response of bridges to vehicle passage is usually studied using moving force and moving mass models. However, the multi-body system allows to consider the quasi-static and dynamic excitation mechanisms. Soil–structure interaction was taken into account by coupling finite elements and boundary elements. The paper presents the results obtained for a simply supported short span bridge in a resonant regime under different soil stiffness conditions. 相似文献
14.
本用震源力学理论和方法研究了徐淮地区从１９７０年以来构造应力场的方向和强度的时空变化过程。结果表明：以唐山地震为分界线，本区的应力场Ｐ轴取向由震前平均６１．８°变为震后平均７７．７°。如果将本区以宿北断裂为界分为南区和北区两个部分，则北区的Ｐ轴取向从６８．１°变７１．２°，而南区的Ｐ轴取向由６２．５°变到８３．６°，南区的变化明显于北区。 相似文献
15.
John L. Tassoulas 《地震工程与结构动力学》2011,40(5):531-550
A half‐space finite element and a transmitting boundary are developed for a water‐saturated layered half‐space using a paraxial boundary condition. The exact dynamic stiffness of a half‐space in plane strain is derived and a second‐order paraxial approximation of the stiffness is obtained. A half‐space finite element and a transmitting boundary are then formulated. The development is verified by comparison of the dynamic stiffness of impermeable and permeable rigid strip foundations with other published results. The advantage of using the paraxial boundary condition in comparison with the rigid boundary condition is examined. It is shown that the paraxial boundary condition offers significant gain and the resulting half‐space finite element and transmitting boundary can represent the effects of a water‐saturated layered half‐space with good accuracy and efficiency. In addition, the numerical method described herein maintains the strengths and advantages of the finite element method and can be easily applied to demanding problems of soil–structure interaction in a water‐saturated layered half‐space. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
16.
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. 相似文献
17.
An up to date literature survey on the dynamics of underground structures is presented briefly. The dynamic response of large three-dimensional underground structures to external or internal dynamic forces or to seismic waves is numerically determined by the frequency domain boundary element method. This method is used to model both the structure and the soil medium, which are assumed to behave as linear elastic or viscoelastic bodies. The full-space dynamic fundamental solution is employed in the formulation and this requires a free soil surface discretization, confined to a finite portion around the area of interest, in addition to soil—structure interface and free structural surface discretizations. The dynamic disturbances can have a harmonic or a transient time variation. The transient case is treated with the aid of numerical Laplace transforms with respect to time. Various numerical examples involving lined cavities and long lined tunnels buried in the full- or the half-space subjected to harmonic or transient external forces or seismic waves are presented to illustrate the method and demonstrate its advantages. 相似文献
18.
19.
20.
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 相似文献