土-结构动力相互作用问题分析中地震动输入的一种新方法

为实现地震作用下土-结构动力相互作用问题的有限元模拟,需要在人工边界上完成地震动的有效输入,目前工程和科研中常用的地震动输入方法有两种:波动输入方法和振动输入方法。波动输入方法的模拟精度高,但实施上相对复杂且耗时,而振动输入方法处理简单,但模拟精度较低。针对应力型人工边界提出一种在人工边界上实现地震动输入的新方法,该方法通过对土-结构有限元模型中由人工边界节点及相邻节点组成的局部子结构施加自由波场位移时程并进行动力分析,从而直接获得可实现地震波动有效输入的等效地震荷载,然后在土-结构有限元模型的人工边界节点上施加等效输入地震荷载并完成动力计算,由此完成土-结构动力相互作用问题的地震动输入和地震反应计算。与原有波动输入方法相比,新方法避免了原方法需分别计算人工边界上自由场应力和由引入人工边界条件引起的附加力,以及需要根据不同人工边界面的外法线方向确定荷载作用方向等较为复杂的处理过程,具有等效地震荷载计算简便、地震动输入过程更易于实施的特点。采用均匀弹性半空间和成层弹性半空间一维地震反应算例初步验证新方法的正确性和可靠性。     

ABAQUS模拟土-结构相互作用时人工边界的选取

利用有限元数值模拟方法研究半无限空间中的动力响应问题时,人工边界的选取对于计算结果的可靠性与准确性有直接的影响。本文利用ABAQUS有限元软件计算土-结构地震响应,分析了选用不同人工边界对计算结果的影响。通过与远置边界自由场和土-结构模型的数值结果对比,讨论了底边界分别采用固定边界和黏性边界;侧边界分别采用固定边界、绑定边界、竖向约束边界和黏性边界时的计算误差。研究表明,利用ABAQUS软件计算土-结构动力相互作用影响时,黏性边界适合作为非刚性地基问题的底边界,固定边界适合作为刚性地基问题底边界。对于侧向边界当采用相对计算宽度较小模型时,固定边界和黏性边界无法再现自由场和土-结构响应;绑定边界和竖向约束边界可以模拟精确的自由场响应,并在动力土-结构作用分析中具有良好的性能。     

用于地下结构地震反应分析的改进反应加速度法

为考虑结构的存在对反应加速度法中地震输入荷载的影响,基于《城市轨道交通结构抗震设计规范》（GB 50909—2014）中的反应加速度法,提出将土-结构模型转换为自由场模型的等效方式进行地震反应分析。采用等效自由场进行地震输入荷载计算,近似地考虑了结构的存在对地震输入荷载的影响,以期提高反应加速度法的计算精度。使用有限元软件ABAQUS,采用改进前后反应加速度法对日本大开车站不同工况下的地震反应进行计算。研究结果表明,采用等效自由场计算的地震输入荷载有效提高了反应加速度法的计算精度。在验证等效方式有效性的基础上,分析等效模型宽度的选取对反应加速度法计算结果的影响,建议等效模型边界至结构侧边的距离取为1～3倍结构宽度。     

考虑河谷场地效应的拱坝-地基地震响应分析方法研究

正河谷地形对场地地震响应有很大影响,但一般进行拱坝-地基体系动力分析时为方便计算,常常将三维拱坝-地基模型的河谷四周侧边界简化为水平成层半空间来求解自由场(本文称为自由场近似解),这与拱坝-地基体系实际受到的地震作用差别是显而易见的。在拱坝-地基动力相互作用地震反应分析中,包括自由场在内的波动输入处理的合理与否,将直接影响     

城市地下变电站非线性地震反应分析

为研究城市地下变电站结构及其内部电气设备的动力特性及地震反应规律,以西安市快速轨道交通二号线张家堡地下变电站为工程背景,依据黏弹性人工边界理论及有限单元法,建立考虑土体-结构-设备动力相互作用的三维有限元计算模型,首先通过半空间自由场模型验证了地震动输入方法的有效性,继而开展地下变电站系统的模态分析和动力时程分析,得到了该类结构的动力特性及地震响应规律。研究结果表明,结构上覆土体及电气设备在地震作用下产生的惯性力是结构构件应力增长的主要原因,8度罕遇地震作用下地下变电站内电气设备的加速度放大系数明显超出规范要求,需要采取相应保护措施。     

三维不规则地形河谷场地地震响应分析方法研究

本文提出了不规则地形地震响应分析的自由场计算方法,并以小湾拱坝场址作为三维场地模型,比较了三维不规则场地侧边界自由场的精确解和近似解;同时,采用大型有限元软件Ansys/Ls-Dyna计算了自由场对场地内部场点地震响应的影响.计算结果表明,三维场地模型的二维自由场近似解与精确解有很大差异,其结果对三维场地内部场点地震响应影响显著,因此,要准确求解三维不规则地形的地震响应必须采用真实的二维侧边界模型.     

基于UPFs的粘弹性人工边界单元及地震动输入方法研究

基于粘弹性人工边界相关理论和地震动输入方法的有关理论,推导并建立了既具有粘弹性人工边界特性,同时可以将地震动转化为施加在人工边界上等效荷载的三维粘弹性人工边界,并利用二次开发工具UPFs将其成功地嵌入至通用有限元分析软件ANSYS中。通过Lamb表面源问题以及内源问题对上述二次开发单元进行了验证。算例表明,该单元可以很好地对无限地基辐射阻尼进行模拟。随后通过对自由场波动输入问题的验证,表明该单元可以实现在有限元单元级别上准确有效地施加地震动等效荷载。与传统的粘弹性边界施加方式相比,该单元有着更简单的使用性以及更强的操作性,同时又不增加计算模型的规模等优势,因而计算效率更高。同时,结合通用有限元分析软件ANSYS强大的建模和计算分析能力,使粘弹性边界有着更广泛的应用,并可适用于各种结构-地基动力相互作用问题的分析。     

相对位置的地上结构对地铁地下结构地震反应的影响

应用有限元方法和粘弹性人工边界,建立了土-结构动力相互作用体系。对比分析了相对位置上不同的地上结构对地铁地下结构地震反应的影响。分析过程采用平面应变条件下的二维动力时程分析,选择的输入地震时程为Kobe时程。结果表明,地上结构的存在对地铁地下结构的层间相对位移、节点内力等地震反应指标产生了不可忽视的影响;且随着两者相对位置的改变,呈规律性变化。     

SV波在饱和土体自由表面的反射

基于饱和两相介质弹性波动方程分析SV波在饱和土体自由表面的反射问题,引入波动方程的势函数解答,求解出二维问题中SV波入射情况下饱和土体自由场的位移、速度、加速度和应力响应。在饱和土体自由场响应解析解基础上,建立SV波入射下饱和土体自由场静、动力有限元模型。建模中考虑了如下几方面因素:(1)在不同分析步,对土体单元赋予不同材料本构。通过*model change命令进行单元生死设定,从而实现在初始应力场平衡的静力状态下采用DuncanChang本构模型,而地震波动输入时采用Davidenkov动力本构模型;(2)采用多孔介质黏弹性人工边界条件,在人工边界上分别施加固相和液相介质的弹簧和阻尼来模拟饱和土体中能量的传播;(3)将地震波转化为作用在人工边界上的等效地震荷载,施加到人工边界节点上;(4)土体单元采用4结点平面应变孔压单元(CPE4P)。有限元计算与解析解比较结果表明:SV波在垂直入射和掠入射时,竖向位移响应为零;在45°左右入射时,水平位移响应最大;60°左右入射时,竖向位移响应最大。这些结论与解析解吻合较好,本文模型为建立土-结构动力相互作用模型打下良好的基础。     

考虑土-结构相互作用的高层建筑抗震分析

本文采用通用有限元程序ANSYS,针对上海地区一例土-箱基-高层建筑结构进行了三维有限元分析,计算中土体的本构模型采用等效线性模型,利用粘一弹性人工边界作为土体的侧向边界,并研究了土体边界位置、土性、基础埋深、基础形式以及上部结构刚度等参数对动力相互作用体系动力特性及地震反应的影响。     

On the effective seismic input for non-linear soil-structure interaction systems

Two equivalent semi-discrete formulations are presented for the problem of the transient response of soil-structure interaction systems to seismic excitation, considering linear behaviour of the soil material and arbitrary non-linear structural properties. One formulation results in a direct method of analysis in which the motion in the structure and the entire soil medium, rendered finite by an artificial absorbing boundary, is determined simultaneously. The other represents a substructuring technique in which the structure and the soil are analysed separately. The forces induced in the discretized system by the incident seismic motion are obtained as part of the general formulation by using the free-field motion of the unaltered soil as the earthquake input. It is shown that these forces act within the soil region in the direct method, but only on the soil-structure interface in the substructure formulation. Both sets of forces, however, involve only the displacements and tractions acting on the fictitious surface in the unaltered (linear) soil which coincides with the soil-structure interface of the complete system. It is shown, further, that the free-field displacements alone define a minimal set of data for evaluating the seismic response of the structure, since the tractions and displacements on that surface are interrelated. In practice, the minimal set must be obtained by extrapolating the available information, as the free-field ground motion at a site is usually specified at a single reference point.     

A seismic free field input model for FE-SBFE coupling in time domain

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     

基于混合波场地震动输入技术的近海场地地震反应分析方法

基于人工边界子结构模型,提出一种利用混合波场实现近海场地中地震P波和SV波垂直输入的方法.该方法中用于波动输入的混合波场由计算模型两侧截断边界的自由波场和底面边界的入射波场构成,避免了不规则近海场地的自由波场求解.同时采用基于声流体单元的流固耦合算法模拟场地-海水动力相互作用,利用流体介质人工边界和黏弹性人工边界单元模...     

Efficient seismic analysis for nonlinear soil-structure interaction with a thick soil layer

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.     

极限安全地震动下考虑土-结构相互作用的核电站安全壳响应分析

根据黏弹性人工边界的基本原理,结合有限元分析软件ABAQUS和MATLAB辅助程序,在地基有限区域上添加黏弹性人工边界并实现极限安全地震动的输入。基于ABAQUS软件平台,对CPR1000安全壳构建了土-结构相互作用体系的数值模拟模型,分析其在极限地震动下的动力响应,并将计算结果与考虑刚性基础的安全壳结构响应数据进行对比。结果表明:核电站CPR1000安全壳结构在极限安全地震动下仍能保持良好的密闭性。考虑土-结构相互作用后分析所得安全壳结构受到的应力、加速度峰值和相对位移均有所增大,使用刚性地基模型要偏于危险。     

A simple boundary element formulation and its application to wavefield excited soil-structure interaction

A simple boundary element formulation which is based directly on the point load solutions for an elastic full-space is presented. It is integrated in a finite element program to calculate dynamic soil-structure interaction problems. The combined boundary and finite element method is applied to structures which are excited by horizontally propagating waves in the soil. For three different types of flexible structure-elastic beams, low and high (square) shear walls-and the corresponding rigid structures the vibration modes and the soil-structure transfer functions have been investigated. The flexible foundations display the same wave pattern as the exciting free-field of the soil, but the amplitudes are reduced with increasing frequency, depending on the stiffness or wave resistance of the structure. Rigid structures show, in part, quite different behaviour, giving free-field reductions caused by kinematic and inertial soil-structure interaction.     

Correlation of predicted seismic response using hybrid modelling with EPRI/TPC lotung experimental data

The hybrid modelling method is presented herein along with the equivalent linearization method to take account of the strain-dependent non-linearity of soils in a soil-structure interaction (SSI) seismic analysis. A refined substructuring of the soil-structure system is utilized and two separate analyses are made to determine the soil free-field and SSI motions induced by earthquake excitation. This method is used to predict the seismic response of a 1/4-scale containment model built in the seismically active area of Lotung, Taiwan. The results obtained show excellent correlation with the field test results.