Stability conditions of explicit integration algorithms when using 3D viscoelastic artificial boundaries

Viscoelastic artificial boundaries are widely adopted in numerical simulations of wave propagation problems. When explicit time-domain integration algorithms are used, the stability condition of the boundary domain is stricter than that of the internal region due to the influence of the damping and stiffness of an viscoelastic artificial boundary. The lack of a clear and practical stability criterion for this problem, however, affects the reasonable selection of an integral time step when using viscoelastic artificial boundaries. In this study, we investigate the stability conditions of explicit integration algorithms when using three-dimensional (3D) viscoelastic artificial boundaries through an analysis method based on a local subsystem. Several boundary subsystems that can represent localized characteristics of a complete numerical model are established, and their analytical stability conditions are derived from and further compared to one another. The stability of the complete model is controlled by the corner regions, and thus, the global stability criterion for the numerical model with viscoelastic artificial boundaries is obtained. Next, by analyzing the impact of different factors on stability conditions, we recommend a stability coefficient for practically estimating the maximum stable integral time step in the dynamic analysis when using 3D viscoelastic artificial boundaries.

     

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

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

粘弹性人工边界单元及地震动输入方法比较研究

本文引用一致粘弹性边界的概念,推导了与一致粘弹性边界等效的粘弹性边界单元,得到了粘弹性边界单元的刚度矩阵和阻尼矩阵。通过对比粘弹性边界单元阻尼系数不同取值的计算精度,给出了阻尼系数推荐取值。在输入地震动时,采用一种基于局部人工边界的子结构输入方法,实现了和波动法相同的输入效果,并且避免了复杂的应力计算和荷载作用方向判断等问题,应用更加方便。采用自由场的算例,比较了人工边界和不同输入方式的组合在SV波垂直入射时的模拟精度。为了说明边界单元的变形对计算结果的影响,提出了使用底边界输入位移幅值u_max和边界单元宽度b的比值ε=u_max/b作为变形指标,调节边界单元边长尺寸。通过算例给出了ε的推荐取值。最后,选取合适的ε后建立一个圆形孔洞散射场地模型,对散射场地模型的模拟结果验证了粘弹性边界单元和子结构输入方式应用于散射场地时仍具有较高精度。     

基于黏弹性边界的地震动输入方法和边界条件选择研究

针对地下工程抗震分析中边界条件和地震动输入方式选择问题,对一地下半无限体模型在一剪切波垂直入射条件下进行不同边界条件和不同地震动输入方法的算例试算。结果发现:两侧采用竖直向位移约束、底部采用黏弹性边界、并在模型底部用等效节点力输入地震动的方法,与三边均采用黏弹性边界条件和三边均采用等效节点荷载输入地震动的方法可以获得同样的结果。该方法易于在ANSYS中通过APDL编程实现,可为地下工程抗震分析中黏弹性边界的使用提供便利。     

弱形式时域完美匹配层——滞弹性近场波动数值模拟

本文旨在构建适用于滞弹性近场时域波动有限元模拟的高精度人工边界条件：完美匹配层（Perfectly Matched Layer：PML）,其中阻尼介质时域本构基于广义标准线性体建立.与以往研究不同,本文采用复坐标延拓技术变换弱形式波动方程构建了可直接用有限元离散的弱形式时域PML,规避以往独立对无限域内波动方程及界面条件进行延拓可导致的PML场方程和界面条件匹配不合理引发数值失稳、计算精度低下等问题.其次,针对PML中多极点有理分式与频域函数乘积的傅里叶反变换难以计算的问题,利用PML精度对复坐标延拓函数中延拓参数微调不敏感这一特点,明确给出了参数微调准则以规避多重极点,进而利用有理分式分解给出了一种普适、简便的计算方法,极大地简化了PML计算.基于该方法可实现任意高阶PML.最后,将本文构建滞弹性PML与高阶勒让德谱元（高精度集中质量有限元）结合得到滞弹性近场波动谱元离散方案.基于算例验证了滞弹性PML的计算效率、精度及新离散方案的长持时稳定特性.新离散方案可应用于计入实际介质阻尼的地震波动正、反问题数值模拟,提高波形模拟的精度以及地下波速结构反演的精度和可靠性.     

Design formulations for supplemental viscous dampers to highway bridges

Design formulas for supplemental viscous dampers to building structures are readily available in FEMA provisions and MCEER research reports. However, for the design of supplemental viscous dampers corresponding to a desired system damping ratio of highway bridges, there exist, if any, few design guidelines. This is particularly true if the bridge components such as elastomeric bearings, piers and abutment possess different damping ratios, stiffnesses, and lumped masses. In this paper, the design formulas for supplemental viscous dampers to highway bridges have been derived based on the concept of ‘composite damping ratio’. The design formulas can be used to determine the damping coefficients of the dampers corresponding to a desired system damping ratio of the bridge in which different component damping ratios may be assumed for the elastomeric bearings, piers and abutments. The proposed design formulas are numerically validated by comparing the seismic responses of a three‐span bridge equipped with viscous dampers with those of the same bridge without viscous dampers but with an assigned inherent system damping ratio equal to the target system damping ratio. Copyright © 2005 John Wiley & Sons, Ltd.     

基于Marc二次开发实现粘弹性人工边界和振型分解反应谱法

新版《水电工程水工建筑物抗震设计规范》（NB35047-2015）中规定:"抗震设防类别为甲类的混凝土坝应考虑远域地基的辐射阻尼效应"。针对通用商业有限元软件Marc在振型分解反应谱法中的不足,且缺少时程分析法中模拟地基辐射阻尼效应的粘弹性人工边界,采用Fortran语言编制相应的独立程序及二次开发程序,以便在Mrac软件中精确实现新规范要求下的混凝土坝抗震安全评价。数值算例和工程实例分析结果验证了基于Marc二次开发实施思路和自编程序的正确性;重力坝抗震薄弱部位主要为坝体断面突变处,考虑无限地基的辐射阻尼效应后,坝体地震动力响应明显下降。     

Seismic response of bridge piers on pile groups for different soil damping models and lumped parameter representations of the foundation

This paper presents a wide parametric study aimed at elucidating the influence, on the computed seismic response of bridge piers, of two related aspects of the model: (1) the adoption of the classical hysteretic or the causal Biot's damping models for the soil and (2) the use of two different lumped parameter models of different complexity and accuracy to approximate the impedances of the pile foundation. A total of 2072 cases, including different superstructures, pile foundations, soil deposits, and seismic input signals, are studied. The results are presented so that the influence of the different parameters involved in the analysis can be assessed. From an engineering point of view, both lumped parameter models provide, in general, sufficiently low errors. The choice of the most adequate model for each case will depend not only on the configuration of the structure and the soil-foundation system but also on the assumed soil damping model, whose influence on the computed seismic responses is relevant in many cases. The nonphysical behaviour provided by the classical hysteretic damping model for the soil at zero frequency generates issues in the process of fitting the impedance functions. It is also found that larger deck displacements are predicted by Biot's model due to the higher damping at low frequencies provided by the classical hysteretic damping model.     

Damping identification of a full‐scale passively controlled five‐story steel building structure

A series of large‐scale dynamic tests was conducted on a passively controlled five‐story steel building on the E‐Defense shaking table facility in Japan to accumulate knowledge of realistic seismic behavior of passively controlled structures. The specimen was tested by repeatedly inserting and replacing each of four damper types, that is, the buckling restrained braces, viscous dampers, oil dampers, and viscoelastic dampers. Finally, the bare steel moment frame was tested after removing all dampers. A variety of excitations was applied to the specimen, including white noise, various levels of seismic motion, and shaker excitation. System identification was implemented to extract dynamic properties of the specimen from the recorded floor acceleration data. Damping characteristics of the specimen were identified. In addition, simplified estimations of the supplemental damping ratios provided by added dampers were presented to provide insight into understanding the damping characteristics of the specimen. It is shown that damping ratios for the specimen equipped with velocity‐dependent dampers decreased obviously with the increasing order of modes, exhibiting frequency dependency. Damping ratios for the specimen equipped with oil and viscoelastic dampers remained constant regardless of vibration amplitudes, whereas those for the specimen equipped with viscous dampers increased obviously with an increase in vibration amplitudes because of the viscosity nonlinearity of the dampers. In very small‐amplitude vibrations, viscous and oil dampers provided much lower supplemental damping than the standard, whereas viscoelastic dampers could be very efficient. Copyright © 2012 John Wiley & Sons, Ltd.     

Dynamic stiffness of deep foundations with inclined piles

The influence of inclined piles on the dynamic response of deep foundations and superstructures is still not well understood and needs further research. For this reason, impedance functions of deep foundations with inclined piles, obtained numerically from a boundary element–finite element coupling model, are provided in this paper. More precisely, vertical, horizontal, rocking and horizontal–rocking crossed dynamic stiffness and damping functions of single inclined piles and 2 × 2 and 3 × 3 pile groups with battered elements are presented in a set of plots. The soil is assumed to be a homogeneous viscoelastic isotropic half‐space and the piles are modeled as elastic compressible Euler–Bernoulli beams. The results for different pile group configurations, pile–soil stiffness ratios and rake angles are presented. Copyright © 2010 John Wiley & Sons, Ltd.     

Impedance functions for foundations embedded in a layered medium: An integral equation approach

An integral equation technique to calculate the dynamic response of foundations embedded in a layered viscoelastic half-space when subjected to external forces and moments is presented. The technique is based on representing the radiated field as resulting from a set of sources distributed over a surface internal to the actual boundary of the foundation and by imposing the boundary conditions in an integral sense. The resulting non-singular integral equation with symmetric kernel is solved by discretization and reduction to a system of linear algebraic equations. The technique is validated by comparison with previous results for cylindrical foundations with different embedment ratios.     

Vibration of three-dimensional rigid raft foundation on viscoelastic medium

Using the coupled model of finite and infinite elements, the compliance of a rigid raft (plate) on a viscoelastic medium is calculated in this paper. Further, the effects of different proportions of radiation damping and material damping in the total damping and the distribution of wave motion in the near field are also investigated. From the numerical results, it can be concluded that: (1) the radiation damping plays an important role in the total damping, compared to the material damping; (2) the material damping has some effect on the compliance of the raft; (3) under the action of harmonic concentrated loads, the displacement amplitudes in the near field gradually decrease with the increase in the distance from the raft and the material damping; (4) the strata characteristics of the soil medium have considerable influence on the compliance of the raft and the distribution of wave motion in the near field.     

Driving a lumped hydrological model with precipitation output from weather generators of different complexity

ABSTRACT

This paper deals with the question of whether a lumped hydrological model driven with lumped daily precipitation time series from a univariate single-site weather generator can produce equally good results compared to using a multivariate multi-site weather generator, where synthetic precipitation is first generated at multiple sites and subsequently lumped. Three different weather generators were tested: a univariate “Richardson type” model, an adapted univariate Richardson type model with an improved reproduction of the autocorrelation of precipitation amounts and a semi-parametric multi-site weather generator. The three modelling systems were evaluated in two Alpine study areas by comparing the hydrological output with respect to monthly and daily statistics as well as extreme design flows. The application of a univariate Richardson type weather generator to lumped precipitation time series requires additional attention. Established parametric distribution functions for single-site precipitation turned out to be unsuitable for lumped precipitation time series and led to a large bias in the hydrological simulations. Combining a multi-site weather generator with a hydrological model produced the least bias.     

大型河谷场地地震动特性研究

采用有限差分方法,通过算例研究了大型河谷场地地震动特性分析中的人工边界的选取方法,对比分析了不同人工边界的选取对数值模拟结果的影响,确定了散射场地震反应分析输入边界的地震动输入方法,认为在进行有限差分动力计算时,模型两侧施加自由场边界的模拟效果要优于两侧施加粘性边界的模拟效果;同时,对FLAC计算软件进行二次开发,对2个地形差异较大的河谷场地,采用线性和非线性摩尔-库仑模型进行了地震反应对比分析,研究了河谷场地地震动幅值和频谱特性随地形变化的规律。模拟数据表明,河谷场地谷底处地震动基本无放大作用,地势凸起处放大作用则较为明显;当考虑土体非线性时,随着地震动强度的增加,放大作用逐渐减小;谱分析结果表明,地势凸起处受高频地震动的影响显著,而地势平坦的谷底则受低频地震动的影响显著。     

Dynamic soil–structure interaction of coupled shear walls by boundary element method

For a class of civil engineering structures, that can be accurately represented by ‘coupled shear walls’ (CSWs), a discrete model for the analysis of the dynamic interaction with the underlying soil is proposed. The CSWs, with one or more rows of openings, rest on a rigid foundation embedded in the elastic or viscoelastic half-space. A hierarchical finite element model based on an equivalent continuum approach is adopted for the structure. A frequency-domain boundary element method is used to represent the half-space. Finally, the set of equations governing the response of the coupled soil-structure system to harmonic lateral loads acting on the structure is also given. The frequency deviation effect with respect to the fixed-base structure and the effects of radiation and material damping in the soil are presented for different characteristics of the structure and different soil properties.     

Dynamic response of partially sealed circular tunnel in viscoelastic saturated soil

Based on Biot's wave equation, dynamic response of a circular tunnel with partially sealed liner in viscoelastic saturated soil is investigated. By introducing two scalar potential functions, the analytical solutions of stresses, displacements and pore pressure induced by axisymmetric gradually applied step load are derived in Laplace transform domain. Numerical results are obtained by inverting Laplace transform presented by Durbin and used to analyze the influences of partial permeable property of boundary and viscoelastic damping coefficient of soil on dynamic response of the tunnel. It is shown that the attenuation of radial displacement appeared with the increase of viscoelastic damping coefficient of soil, and relative rigidity of liner and soil, and the influence of partial sealing property of boundary on stresses, displacements and pore pressure is remarkable. The available solutions of permeable and impermeable boundary conditions are only two extreme cases of this paper.     

起伏地表复杂介质波动方程有限元数值模拟方法

波动方程数值模拟是深入研究地震波传播规律的有效方法.有限差分法因其方法简单、精度高而得到了广泛的应用.但其缺点是不能准确模拟具有复杂几何形态的物性界面.因而当遇到起伏地表或复杂构造时,求解精度低.为了准确模拟起伏地形、复杂构造和复杂介质条件下的地震波场,本文采用有限元法模拟二维声波方程.用三角形单元模拟地形和速度界面;把单元内的场和波速均看作单元上的线性函数,以适应复杂介质压制边角散射;采用吸收边界条件去除来自截断边界上的反射;采用集中质量矩阵和集中阻尼矩阵使得显式时间递推无需对矩阵求逆,提高了计算效率.对模型的计算表明该方法正确有效.     

水平成层场地地震反应的集中质量切比雪夫谱元分析方法

提出了一种用于水平成层场地地震反应分析的时域高阶显式算法. 首先,将覆盖土层和基岩划分为若干个切比雪夫谱单元,在模型底部设置多次透射人工边界;其次,以切比雪夫正交多项式构建高阶单元位移模式,通过高斯?洛巴托积分严格导出对角形式的切比雪夫谱单元集中质量矩阵,结合中心差分时域逐步积分格式,建立了高效的集中质量切比雪夫谱元波动模拟方法;最后,利用日本Kik-net强震台网提供的不同类型场地上获得的实际地震观测记录检验了本文方法的有效性. 该方法避免了传统切比雪夫谱元法由于具有一致质量矩阵形式而造成的计算效率不高的问题。数值结果表明,本文方法能够较好地预测Ⅰ₁,Ⅱ和Ⅳ类场地在较弱地震和中等强度地震作用下的地面运动特征,每个波长内仅需布置少量单元即可取得较高精度的计算结果。      

Structure-soil-structure interaction of adjacent buildings subjected to seismic loading

The problem of the through-soil coupling of structures has puzzled the researchers in the field for a long while, especially regarding the varied performance of identical, adjacent buildings in earthquakes. The phenomenon of structure-soil-structure interaction (SSSI) that has often been overlooked is recently being recognized: The possible effects in urban regions are yet to be thoroughly quantified. In this respect, the goal of this work was to rigorously investigate the interacting effects of adjacent buildings in a two-dimensional setting. Detailed finite element models of 5-, 15-, and 30-story structures, realistically designed, were used in forming building clusters on the viscoelastic half-space. Perfectly matched layers were used to properly define the half-space boundaries. The interaction of the structure and the soil medium because of the presence of spatially varying ground motion on the boundary of excavated region was considered. The effects of the foundation material and the distance between adjacent buildings on the structural behavior of the neighboring buildings were investigated using drift ratios and base shear quantities as the engineering demand parameters of interest. The effects of SSSI, first investigated in the frequency domain, was then quantified in the time domain using suites of appropriate ground motions in accordance with the soil conditions, and the results were compared with the counterpart SSI solution of a single building. The results showed that, for identical low-rise structures, the effects of SSSI were negligible. Yet, neglecting SSSI for neighboring closely spaced high-rise structures or building clusters with a large stiffness contrast was shown to lead to a considerable underestimation of the true seismic demands even compared with solutions obtained using the rigid base assumption.     

正常工况下大型岸桥地震瑞利阻尼系数计算方法及试验研究

针对大型岸桥地震响应分析中不同的瑞利阻尼系数会导致计算结果存在较大偏差,基于相似模型试验数据与仿真计算结果的对比分析,并综合考虑岸桥边界约束等因素,借鉴工程中常用的瑞利阻尼系数计算方法改进了岸桥结构的阻尼系数计算方法及频率选取方式。利用仿真方法和振动台模型地震试验对比分析不同瑞利阻尼系数下岸桥结构的动力响应,结果表明：瑞利阻尼系数计算方法和频率选取对岸桥结构地震响应的影响随着地震加强而提高,当加速度峰值为0.4g、0.62g,现有的阻尼系数计算方法得到的仿真结果与试验值之间的平均误差超过20%,而采用改进后的阻尼系数计算方法,岸桥地震仿真结果与试验值误差都在10%以内。