共查询到20条相似文献，搜索用时 569 毫秒
1.
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. 相似文献
2.
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 相似文献
3.
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 相似文献
4.
A systematic procedure to construct the (symmetric) static-stiffness, damping and mass matrices representing the unbounded medium is presented addressing the unit-impulse response matrix corresponding to the degrees of freedom on the structure–medium interface. The unit-impulse response matrix is first diagonalized which then permits each term to be modelled independently from the others using expansions in a series of Legendre polynomials in the time domain. This leads to a rational approximation in the frequency domain of the dynamic-stiffness coefficient. Using a lumped-parameter model which provides physical insight the property matrices are constructed. 相似文献
5.
Starting from the unit-impulse response matrix of the unbounded medium, a discrete-time formulation permitting the recursive evaluation of the interaction forces and a continuous-time formulation yielding property matrices corresponding to a model with a finite number of degrees of freedom are discussed. This is achieved using the balancing approximation method which is easily automated, guarantees stability and leads to highly accurate results. © 1998 John Wiley & Sons Ltd. 相似文献
6.
Formulation of a matrix‐valued force–displacement relationship which can take radiation damping into account is of major importance when modelling unbounded domains. This can be done by means of fundamental solutions in space and time in connection with convolution integrals or by means of a frequency dependent boundary element representation, but for discrete frequencies Ω only. In this paper a method for interpolating discrete values of dynamic stiffness matrices by a continuous matrix valued rational function is proposed. The coupling between interface degrees of freedom is fully preserved. Another crucial point in soil–structure interaction analysis is how to implement an approximation in the spectral domain into a time‐domain analysis. Well‐known approaches for the scalar case are based on the partial‐fraction expansion of a scalar rational function. Here, a more general procedure, applicable to MDOF‐systems, for the transformation of spectral rational approximations into the time‐domain is introduced. Evaluation of the partial‐fraction expansion is avoided by using the so‐called mixed variables. Thus, unknowns in the time‐domain are displacements as well as forces. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
7.
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. 相似文献
8.
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. 相似文献
9.
Consideration of structure-foundation-soil dynamic interaction is a basic requirement in the evaluation of the seismic safety of nuclear power facilities. An efficient and accurate dynamic interaction numerical model in the time domain has become an important topic of current research. In this study, the scaled boundary finite element method （SBFEM） is improved for use as an effective numerical approach with good application prospects. This method has several advantages, including dimensionality reduction, accuracy of the radial analytical solution, and unlike other boundary element methods, it does not require a fundamental solution. This study focuses on establishing a high performance scaled boundary finite element interaction analysis model in the time domain based on the acceleration unit-impulse response matrix, in which several new solution techniques, such as a dimensionless method to solve the interaction force, are applied to improve the numerical stability of the actual soil parameters and reduce the amount of calculation. Finally, the feasibility of the time domain methods are illustrated by the response of the nuclear power structure and the accuracy of the algorithms are dynamically verified by comparison with the refinement of a large-scale viscoelastic soil model. 相似文献
10.
Exact representation of unbounded soil contains the single output–single input relationship between force and displacement in the physical or transformed space. This relationship is a global convolution integral in the time domain. Rational approximation to its frequency response function (frequency‐domain convolution kernel) in the frequency domain, which is then realized into the time domain as a lumped‐parameter model or recursive formula, is an effective method to obtain the temporally local representation of unbounded soil. Stability and identification for the rational approximation are studied in this paper. A necessary and sufficient stability condition is presented based on the stability theory of linear system. A parameter identification method is further developed by directly solving a nonlinear least‐squares fitting problem using the hybrid genetic‐simplex optimization algorithm, in which the proposed stability condition as constraint is enforced by the penalty function method. The stability is thus guaranteed a priori. The infrequent and undesirable resonance phenomenon in stable system is also discussed. The proposed stability condition and identification method are verified by several dynamic soil–structure‐interaction examples. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
11.
John P. Wolf 《地震工程与结构动力学》2002,31(1):15-32
The scaled boundary finite‐element method is a powerful semi‐analytical computational procedure to calculate the dynamic stiffness of the unbounded soil at the structure–soil interface. This permits the analysis of dynamic soil–structure interaction using the substructure method. The response in the neighbouring soil can also be determined analytically. The method is extended to calculate numerically the response throughout the unbounded soil including the far field. The three‐dimensional vector‐wave equation of elasto‐dynamics is addressed. The radiation condition at infinity is satisfied exactly. By solving an eigenvalue problem, the high‐frequency limit of the dynamic stiffness is constructed to be positive definite. However, a direct determination using impedances is also possible. Solving two first‐order ordinary differential equations numerically permits the radiation condition and the boundary condition of the structure–soil interface to be satisfied sequentially, leading to the displacements in the unbounded soil. A generalization to viscoelastic material using the correspondence principle is straightforward. Alternatively, the displacements can also be calculated analytically in the far field. Good agreement of displacements along the free surface and below a prism foundation embedded in a half‐space with the results of the boundary‐element method is observed. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
12.
13.
The seismic response of a dam is strongly influenced by its interaction with the water reservoir and the foundation. The hydrodynamic forces in the reservoir are in turn affected by radiation of waves towards infinity, wave absorption at the reservoir bottom, and cross-coupling between the foundation below the dam and the reservoir bottom. The fluid–foundation interaction effect, i.e. the wave absorption along the reservoir bottom, can be accounted for by using either an approximate one-dimensional (1D) wave propagation model or a rigorous analysis of interaction between the flexible soil along the base and the water. The rigorous approach requires enormous computational effort because of (a) cross-coupling between the foundation of the dam and the soil below the reservoir and (b) frequency dependence of the boundary condition along the fluid-foundation interface. The analysis can be simplified by ignoring the cross-coupling and by using the approximate 1D wave propagation model. The effects of each of these two simplifications on the accuracy and computational efficiency of the procedure used for the seismic response analysis of a dam are examined. Analytical results are presented for the complex frequency-response functions as well as the time histories of the response of Pine Flat dam to Taft and E1 Centro ground motions. 相似文献
14.
A new predictor–corrector (P–C) method for multi‐site sub‐structure pseudo‐dynamic (PSD) test is proposed. This method is a mixed time integration method in which computational components separable from experimental components are solved by implicit time integration method (Newmark β method). The experiments are performed quasi‐statically based on explicit prediction of displacement. The proposed P–C method has an important advantage as it does not require the determination of the initial stiffness values of experimental components and is thus suitable for representing elastic and inelastic systems. A parameter relating to quality of displacement prediction at boundaries nodes is introduced. This parameter is determined such that P–C method can be applicable to many practical problems. Error‐propagation characteristics of P–C method are also presented. A series of examples including linear and non‐linear soil–foundation–structure interaction problem demonstrate the performance of the proposed method. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
15.
基础动力刚度的精确数值解及集中参数模型 总被引：2，自引：1，他引：2
土－结构相互作用分析的关键是建立以土－结构界面定义的无限半空间的动力刚度矩阵。本文介绍了一种求解半无限地基动力刚度的新方法，通过两个算例验证了该方法的精度，并给出了一种利用频域刚性基础动力刚度计算基础时域荷载响应的实用方法，该研究为刚性基础设计提供了一种新的，可靠的理论方法。 相似文献
16.
《地震工程与结构动力学》2018,47(5):1291-1308
Evaluating the inelastic seismic response of structures accurately is of great importance in earthquake engineering and generally requires refined simulation, which is a time‐consuming process. Because the material nonlinearity generally occurs in a small part of the whole structure, many researches focus on taking advantage of this characteristic to improve the computational efficiency and the inelasticity‐separated finite element method (IS‐FEM) proposed recently provide a generic finite element formulation for solving this kind of problems efficiently. Although the fiber beam‐column element is widely used for the simulation of reinforced concrete (RC) framed structures, the inelastic deformation is often detected in a large part of the numerical model under earthquake excitation so that it is hard to achieve high efficient computation when applying the IS‐FEM to the inelastic response analysis of RC fiber models directly. In this paper, a new numerical scheme for seismic response analysis of RC framed structures model by fiber beam‐column element is proposed based on the IS‐FEM. To implement the RC fiber model for use in IS‐FEM and improve the computational performance of proposed scheme, a method of identifying the local domains with severe section inelasticity level is proposed and a modified Kent‐Park concrete material model is developed. Because the Woodbury formula is adopted as the solver, the global stiffness matrix can keep unchanged throughout the analysis and the main computational effort is only invested on a small matrix representing local inelastic behavior. The numerical examples demonstrate the validity and efficiency of the proposed scheme. 相似文献
17.
P. Dangla 《地震工程与结构动力学》1988,16(8):1115-1128
A plane strain model for dynamic soil-structure interaction problems under harmonic state is presented. The boundary element method is used to study the response of a homogeneous isotropic linear elastic soil. The far field displacement at the free surface is approximated by an outgoing Rayleigh wave. The finite element method is used to describe the response of the building, of the foundation and possibly of a finite part of the inhomogeneous non-linear soil. Two coupling procedures are described. The model is applied to a problem previously studied in the antiplane case. Incident P, SV and Rayleigh waves are considered. The results show an amplification and an attenuation of the structure motion with frequency when incident Rayleigh waves and P, SV body waves are respectively considered. 相似文献
18.
阻尼矩阵的选择对土层时域内地震反应计算的准确性有重要意义。本文构造了质量比例阻尼、刚度比例阻尼、Rayleigh阻尼、Caughey阻尼和Clough阻尼矩阵,研究各阻尼对土层地震反应分析精度的影响。以苏通大桥5号桥塔基础处的工程场地为例,分别以人工合成基岩波,Northridge波、Parkfield波和汶川波为输入,得到不同阻尼模型下土层的地震反应,结果表明Clough阻尼在计算精度和误差稳定性上要优于其他阻尼模型,针对Clough阻尼的不足之处提出了改进方法,算例表明改进后的Clough阻尼能在保证计算精度的前提下减少计算机的存储量和计算量。 相似文献
19.
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. 相似文献
20.
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. 相似文献