Mode-based equivalent multi-degree-of-freedom system for one-dimensional viscoelastic response analysis of layered soil deposit

Discrete models such as the lumped parameter model and the finite element model are widely used in the solution of soil amplification of earthquakes. However, neither of the models will accurately estimate the natural frequencies of soil deposit, nor simulate a damping of frequency independence. This research develops a new discrete model for one-dimensional viscoelastic response analysis of layered soil deposit based on the mode equivalence method. The new discrete model is a one-dimensional equivalent multi-degree-of-freedom (MDOF) system characterized by a series of concentrated masses, springs and dashpots with a special configuration. The dynamic response of the equivalent MDOF system is analytically derived and the physical parameters are formulated in terms of modal properties. The equivalent MDOF system is verified through a comparison of amplification functions with the available theoretical solutions. The appropriate number of degrees of freedom (DOFs) in the equivalent MDOF system is estimated. A comparative study of the equivalent MDOF system with the existing discrete models is performed. It is shown that the proposed equivalent MDOF system can exactly present the natural frequencies and the hysteretic damping of soil deposits and provide more accurate results with fewer DOFs.     

高温气冷堆核电厂建筑结构采用不同有限元模型的模态分析对比研究

建立有限元模型是核电厂建筑结构模态分析的重要前提。本文以某高温气冷堆核电厂建筑结构为原型，在分析方法相同的前提下，建立2种不同模型（Solid模型和Shell模型），并对这2种模型进行模态分析。重点分析、对比2种模型的自振频率和振型图，计算分析表明:Solid模型与Shell模型相比，计算得到的结构自振频率值较高，但两者的差异很小，前30阶自振频率相对误差小于3.4%；2种模型的计算结构振型基本一致。研究结果可为核电厂抗震性能分析和设计提供参考。     

Asymptotic formulae for correcting finite element predicted natural frequencies of gravity and embankment dams

Extraction of natural frequencies of a gravity dam or an embankment dam plays an important role in the seismic design of the dam because the seismic response of a dam is dependent largely on the dynamic characteristics of the dam. Owing to the lack of exact solutions and the geometry of a dam, numerical methods such as finite element methods have been often used to extract the natural frequencies of the dam. Since the finite element method is an approximate one, the resulting finite element solution to the natural frequency of a dam cannot be safely used unless its accuracy is evaluated within the acceptable range for the seismic design of the dam. To solve this problem, some asymptotic formulae for correcting the finite element predicted natural frequencies of a gravity dam and an embankment dam have been developed in this paper. Since the present asymptotic formulae are derived from the fact that the finite element solution tends to the exact one if the finite element size used approaches zero, they provide a corrected solution of higher accuracy for the natural frequency of a dam so that the accuracy of a finite element solution can be evaluated against this corrected solution. After the correctness and usefulness of the present formulae are assessed, two practical examples have been given to show how the asymptotic formulae can be used to correct and evaluate the discretization error for the finite element predicted natural frequencies of gravity dams and embankment dams.     

Mass lumping models of the linear diffusion equation

The nodal domain integration method is used to develop a numerical model of the linear diffusion equation. The nodal domain integration approach is shown to represent an infinity of finite element mass matrix lumping schemes including the Galerkin and subdomain integration versions of the weighted residual method and an integrated finite difference method. Neumann, Dirichlet and mixed boundary conditions are accommodated analogous to the Galerkin finite element method. In order to reduce the overall integrated approximation relative error, a mass matrix lumping formulation is developed which is based on the Crank-Nicolson time advancement approximation. The optimum mass lumping factors are found to be strongly related to the model timestep size.     

Determination of scoured bridge natural frequencies with soil–structure interaction

This study developed a finite element method with the effect of soil–fluid–structure interaction to calculate bridge natural frequencies. The finite element model includes bridge girders, piers, foundations, soil, and water. The effective mass above the soil surface was then used to find the first natural frequency in each direction. A field experiment was performed to validate that the natural frequencies calculated using the proposed finite element method had acceptable accuracy. The calculated natural frequencies with the fluid–structure interaction effect are always smaller than those without this effect. However, the frequency change due to the fluid effect is not obvious, so using the soil–structure interaction model is accurate enough in the bridge natural frequency analysis. The trend of the frequency decreases with the increase of the scour depth, but the curve is not smooth because of non-uniform foundation sections and layered soils. However, when the scour depth is such that pile cap is exposed, the changes in natural frequency with the scour depth are more obvious, and this is useful for measurement of the depth using bridge natural frequencies.     

Analysis of the observed seismic response of a highway bridge

Strong-motion accelerograms obtained on the San Juan Bautista 156/101 Separation Bridge during the 6 August 1979 Coyote Lake, California, earthquake are used to examine the response of this multiple-span bridge to moderate levels of earthquake loading. Although the bridge was not damaged, the records are of significant engineering interest as they are the first to be recorded on a highway bridge structure in North America. A technique of system identification is used to determine optimal modal parameters for linear models which can closely replicate the observed time-domain seismic response of the bridge. Time variations in frequency and damping in the horizontal response are identified using a moving-window analysis. A three-dimensional finite element model is developed to study the bridge response in detail. The first two horizontal modal frequencies computed from this model are in excellent agreement with information obtained during the system identification analysis provided the finite element model's expansion joints are locked, preventing relative translational motions from occurring across the joints. Locking is confirmed by the observed seismic deformations of the structure in the fundamental mode. Fundamental vertical frequencies of the individual spans, predicted by the finite element model, are in very good agreement with ambient vibration test data.     

Vibration tests and analysis of a model arch dam

Vibration tests were conducted on a 1/24-scale model of the North Fork Dam, a double-curvature arch dam, to determine natural frequencies, mode shapes and hydrodynamic pressures. The mode shapes and natural frequencies were determined from tests using two vibrators mounted on the crest of the dam. Hydrodynamic pressures at the dam/reservoir interface were determined from tests in which the vibrator was attached to the downstream foundation of the dam. The hydrodynamic pressures calculated using Westergaard's theory and a theory for arch dams developed by Perumalswami and Kar accurately predicted the measured pressure at frequencies below the first mode frequency of the dam. The differences in the two theories were insignificant. The Structural Analysis Program (SAP), a linear three-dimensional (3-D) finite element code, was used to compute mode shapes and frequencies for the dam with its base fixed and with a foundation. Numerical solution schemes used in the finite element analysis consisted of a Ritz analysis and a subspace iteration method. Calculations were conducted for both full and empty reservoir conditions. The accuracy of the Ritz analysis improved considerably as more nodes in flexible regions of the dam were loaded. However, the lowest eigenvalues were computed using the subspace iteration method. For the full reservoir, the natural frequencies decreased by 20-30 per cent when the foundation was included in the finite element model. The difference was less when the reservoir was empty. The calculations using the subspace iteration scheme and including the foundation agreed closely with experimental mode shapes and corresponding natural frequencies.     

高墩抗震梁桥单墩模型地震时程响应分析

直接基于经典的欧拉梁理论,推导针对规则高墩梁桥单墩模型弹性阶段地震时程响应的半解析计算步骤,考虑梁墩铰接的形式,在墩底引入等效基础弹簧变形协调条件,解析的获得各阶实模态和频率方程,用牛顿法搜索各阶频率．应用振型叠加法求解体系的地震时程响应。最后应用该方法对一高墩梁桥单墩模型的地震响应进行分析,与离散模型的有限元时程积分的结果比较表明了方法的有效性。     

The use of reduced finite element models in system identification

An a priori, analytical model in system identification of vibrating structures is improved by input and output measurements with least square fitting. The structure is modelled by the finite element method. Finite element models usually need a large number of degrees of freedom to simulate a small number of lower spectrum eigenfrequencies with accuracy. The large finite element model is reduced to a subspace of significant eigenfrequencies. The size of the subspace is chosen with regard to the frequency content of the measured data and the accuracy of the large analytical model. An identification method is formulated for the large analytical model. This procedure improves system parameters in the matrices of the large model by measured input and output data with a least square functional. The objective function is consistently reduced, so that the whole identification procedure can be performed in the small subspace. The proposed reduction method permits very large and accurate analytical models to be used, and it decreases the computational cost of the identification procedure significantly. The computational efficiency is demonstrated on an in situ experiment of a radar tower.     

Energy-based dynamic parameter identification for Pasternak foundation model

Parameter identification of Pasternak foundation models(PFM) is never satisfactory, which discourages the application and popularization of PFM. In the present study, an energy-based model to predict the dynamic foundation coefficients was proposed using the vibration kinetic energy and potential energy of a Pasternak foundation-rigid plate system. On the basis of the Pasternak foundation, the relationship among the natural frequency, dynamic foundation coefficients, rigid plate configuration, and vibrating soil equivalent mass per unit area was considered. To obtain the natural frequencies of the Pasternak foundation-rigid plate system, dynamic tests were performed. Using two or more dynamic test results of various rigid plates on a foundation, a set of equations of dynamic foundation coefficients was set up to directly identify the foundation coefficients and equivalent mass per unit area of vibrating soil. The feasibility of the proposed method was verified by comparing it with the outdoor and indoor test results and finite element analysis results. When the proposed method is used to obtain the dynamic parameters, PFM can be generalized and applied more widely in engineering practice.     

Rapid evaluation and damage assessment of instrumented highway bridges

This study focuses on the use of strong motion data recorded during earthquakes and aftershocks to provide a preliminary assessment of the structural integrity and possible damage in bridges. A system identification technique is used to determine dynamical characteristics and high‐fidelity first‐order linear models of a bridge from low level earthquake excitations. A finite element model is developed and updated using a genetic algorithm optimization scheme to match the frequencies identified and to simulate data from a damaging earthquake for the bridge. Here, two criteria are used to determine the state of the structure. The first criteria uses the error between the data recorded or simulated by the calibrated nonlinear finite element model and the data predicted by the linear model. The second criteria compares relative displacements of the structure with displacement thresholds identified using a pushover analysis. The use of this technique can provide an almost immediate, yet reliable, assessment of the structural health of an instrumented bridge after a seismic event. Copyright © 2011 John Wiley & Sons, Ltd.     

Development of basic technique to improve seismic response accuracy of tributary area-based lumped-mass stick models

Although a detailed finite element(FE) model provides more precise results, a lumped-mass stick(LMS) model is preferred because of its simplicity and rapid computational time. However, the reliability of LMS models has been questioned especially for structures dominated by higher modes and seismic inputs. Normally, the natural frequencies and dynamic responses of a LMS model based on tributary area mass consideration are different from the results of the FE model. This study proposes a basic updating technique to overcome these discrepancies; the technique employs the identical modal response, D(t), to the detailed FE model. The parameter D(t) is a time variable function in the dynamic response composition and it depends on frequency and damping ratio for each mode, independent of the structure's mode shapes. The identical response D(t) for each mode is obtained from the frequency adaptive LMS model; the adaptive LMS model which can provide identical modal frequencies as the detailed FE model. Theoretical backgrounds and formulations of the updating technique are proposed. To validate the updating technique, two types of structures(a symmetric straight column and an unsymmetric T-shaped structure) are considered. From the seismic response results including base shear and base moment, the updating technique considerably improves the seismic response accuracy of the tributary area-based LMS model.     

土石坝地震响应的非线性剪切条模型与对比分析

对一维剪切条计算模型进行改进,提出了土石坝非线性地震反应的简化计算方法。首先将坝体沿坝高离散为一系列的具有不同剪切刚度与阻尼比等参数特性的层状体系,建立了各层的振动控制方程及其边值条件,进而采用数学物理方程方法进行了求解,确定了体系的振动特性,并根据振型叠加原理和Duhamel积分确定了坝体地震反应的线弹性解。采用等价线性化方法考虑坝料的动力非线性性质,通过对线弹性地震响应的反复迭代计算,使得各层土的模量和阻尼比与其相应的剪应变水平相协调,确定出与非线性坝体系统相等效的线性解答,并将所得到的地震响应作为非线性地震响应的近似解。最后,以均质坝和心墙坝作为算例进行了具体的数值计算,将所得结果与有限元数值解进行对比分析,论证了所提方法的适用性和合理性。     

有限元离散模型中的出平面运动形式

波动在有限元离散模型中传播时,会出现许多不同于连续介质中波动的现象.本文引用经典物理学中晶格动力学的分析方法,由有限元离散形式的出平面波动方程得到频散关系,在此基础上,讨论了集中质量有限元模型中各种可能的出平面运动形式,并分析了考虑时域离散化影响下各种形式可能存在的频率范围.     

运行状态下海上单桩风机系统自振频率分析

系统自振频率限制是海上风机结构设计中的一个关键因素。运行状态下风机动力荷载会引起基础的水平侧移,较大的水平侧移会导致基础刚度的降低,进一步影响风机系统的自振频率。该文基于有限元软件ABAQUS平台,建立单桩式海上风机结构系统的自振频率数值模型,并讨论运行状态下基础水平侧移对大直径海上风机系统自振频率的影响。模型中考虑了塔筒的变截面特性;桩-土相互作用通过p-y曲线方法模拟;桩和塔采用梁单元模拟;通过Pushover分析汇总出水平侧移引起的桩顶水平刚度。研究结果表明:桩基侧向位移会降低风机结构体系的自振频率;桩基侧向位移对基频的影响较小,对高阶频率的影响显著;大直径海上风机的频率计算中可忽略风机运行状态对体系自振频率的影响。     

CAP1400核岛结构缩尺模型对地震反应影响分析

针对核电厂CAP1400核岛结构地震反应问题,构建了核岛屏蔽厂房和辅助厂房整体结构的3个分析模型:原型和1/16、1/40缩尺模型,并在AP000谱和RG1.60谱地震动输入下进行了有限元模拟对比分析,探讨了振动台试验模型缩尺处理的合理性和精确性。研究表明,基于缩尺模型得到的结构自振频率相对于原型结构模型有所降低,降低幅度在8.5%以内;结构模型的缩尺对结构反应峰值加速度和高频（大于3Hz）加速度反应谱的影响较为显著,但对较低频（小于3Hz）的加速度反应谱影响较小;模型缩尺对结构不同方向反应的影响中,刚度越大的方向其影响越大。进一步将结构模型数值模拟结果与1/16缩尺模型的振动台试验结果进行了比较分析,试验给出的结构自振频率远低于模型数值模拟结果,但原型和1/16缩尺模型数值模拟得到的结构反应均与试验结果较为接近。基于模型数值模拟和振动台试验研究,认为对于缩尺比1/16或更大的模型可以忽略模型的缩尺效应。     

基于卷积完全匹配层的非规则网格时域有限元探地雷达数值模拟

复频移完全匹配层(Complex Frequency-Shifted PML,CFS-PML)在长时间时域计算中对凋落波、倏失波具有好的吸收效果,并被广泛应用于时域有限差分模拟中.而本文采用卷积方法将CFS-PML应用于时域有限元求解GPR波动方程的数值模拟中.论文以TM波为例,推导了基于CPML(Convolutional PML)边界的时域有限元GPR波动方程求解公式,采用Newmark-β方法对时间导数进行离散,有效改善了时域有限元GPR数值计算程序的稳定性.并以狭长模型为例,开展了CPML边界中关键参数m、R和κ的选取实验,通过对比反射误差大小确定了综合最优参数组合.相同时刻UPML与CPML波场快照、3个检测点的反射误差比较,说明CPML较UPML具有更好的吸收效果.最后,采用非规则四边形网格对1个复杂GPR模型进行剖分,应用加载CPML边界条件的FETD程序对该模型进行了正演,得到了二维剖面法、宽角法正演GPR剖面图,说明非规则四边形对复杂模型的良好适应性,基于CPML边界条件的FETD可有效减少边界反射误差,能实现对任意复杂不规则模型的正演模拟.     

Finite element model updating for structures with parametric constraints

This paper presents a finite element (FE) model updating procedure applied to complex structures using an eigenvalue sensitivity‐based updating approach. The objective of the model updating is to reduce the difference between the calculated and the measured frequencies. The method is based on the first‐order Taylor‐series expansion of the eigenvalues with respect to some structural parameters selected to be adjusted. These parameters are assumed to be bounded by some prescribed regions which are determined according to the degrees of uncertainty that exist in the parameters. The changes of these parameters are found iteratively by solving a constrained optimization problem. The improvement of the current study is in the use of an objective function that is the sum of a weighted frequency error norm and a weighted perturbation norm of the parameters. Two weighting matrices are introduced to provide flexibility for individual tuning of frequency errors and parameters' perturbations. The proposed method is applied to a 1/150 scaled suspension bridge model. Using 11 measured frequencies as reference, the FE model is updated by adjusting ten selected structural parameters. The final updated FE model for the suspension bridge model is able to produce natural frequencies in close agreement with the measured ones. Copyright © 2000 John Wiley & Sons, Ltd.     

三维陆地可控源电磁法有限元模型降阶快速正演

三维陆地可控源电磁法有限元快速正演的主要瓶颈在于多频率大型稀疏方程组求解问题.本文引入一种基于模型降阶的Krylov子空间投影算法,推导了有限元刚度矩阵的模型降阶形式,构建了频率域传递函数;采用标准正交向量序列,构建一个远远小于有限元刚度矩阵维度的矩阵,该矩阵与频率无关,通过一次模型降阶即可实现多频点有限元方程快速求解.采用基于电场的变分方程,加入散度校正条件,以消除伪解;引入伪δ函数,消除了源点的奇异性,可适用于复杂背景模型三维有限元数值模拟,并为多源的求解奠定了基础;以层状介质模型解析解为标准,通过和基于Pardiso直接求解器的有限元算法（3DFEM）进行比较,模型降阶法计算时间小于前者的1/10,平均相对误差在1.72%,在满足精度要求下,实现了高效率三维有限元数值求解;分别设计了横向高低阻模型和纵向高低阻模型,分析了从近区到远区电场和卡尼亚视电阻率的变化规律,假极值的表现特征,阴影效应的影响等,从而也验证了该算法的正确性.最后,建立了一个地层陷落柱模型,通过模型降阶有限元正演模拟,发现视电阻率断面图在陷落柱上方出现"凹陷",与模型设计吻合,表明该算法对复杂地层模拟具有同样的适用性.     

面向目标自适应海洋可控源电磁三维矢量有限元正演

本文实现了一种面向目标自适应海洋可控源电磁三维矢量有限元方法.为满足三维复杂电性结构模拟的需求,网格剖分采用非结构化六面体.在组装刚度矩阵之后,形成的大型复数线性方程组分解为等价的实数形式,利用带预条件的广义最小残差法进行求解.在获得微分方程的解之后,为提高解的准确性,通过面向目标的自适应误差估计来指示网格细化,重点加密能使观测点数值模拟精度提高的网格.对于大规模三维数据,为了使模型空间的并行计算达到均衡负载的效果,我们使用METIS函数库来进行网格计算任务量的划分.最后,通过对比一维解析解与三维自适应矢量有限元计算结果,验证了程序的正确性;通过自适应过程中误差指示子的分布,验证了面向目标自适应的有效性;通过对三维复杂模型进行均衡负载下的并行计算,测试了程序的可扩展性.