立式储液罐地震破坏快速评估方法

本文采用大型有限元程序ANSYS对立式储液罐的地震反应进行了非线性数值分析,考虑了罐体和内部液体的相互作用,得到了3个不同容积的储罐,在不同峰值的地震动时程输入下的反应结果。结合地震灾害中立式储罐的地震破坏现象和特点,以及现有的立式储罐地震破坏预测方法,提出了立式储罐的地震破坏快速评估方法。     

基础隔震结构的耗能分析

采用Bouc-wen模型,利用状态空间迭代法,对基础隔震结构进行了多质点的弹塑性时程分析,并根据此结果,利用能量方程,求得隔震结构的各项能量,绘制了各能量项时程曲线。以一实际工程为例,求得隔震结构的各项耗能情况,说明了基础隔震结构以减少地震输入和隔震层滞回耗能来减小对上部结构的损坏。     

自适应时间步长法及三维堤坝地震液化数值模拟

三维大模型数值计算因巨大的单元和结点数目而非常耗时,在地震响应分析中受计算时间步长的限值则更加耗时。在饱和砂土动力液化计算平台上开发时域离散误差评估方法和时间步长自适应调整的计算程序,并成功应用于三维堤坝地震液化响应分析。时域离散误差包括土骨架的位移误差和单元孔压误差,通过定义孔压误差影响系数计算出混合误差,根据混合误差和设定的误差允许值进行计算步长的自适应调整。在三维堤坝地震液化数值模拟中,采用自适应时间步长法有效避免小步长精确但耗时、大步长省时而不精确的缺点。在大模型和超大模型计算中,最优调整每一步的计算时间步长,完美实现既节省时间又不失精度的时域离散策略。     

Optimal design of structures for earthquake loads by a hybrid RBF-BPSO method

The optimal seismic design of structures requires that time history analyses （THA） be carried out repeatedly. This makes the optimal design process inefficient, in particular, if an evolutionary algorithm is used. To reduce the overall time required for structural optimization, two artificial intelligence strategies are employed. In the first strategy, radial basis function （RBF） neural networks are used to predict the time history responses of structures in the optimization flow. In the second strategy, a binary particle swarm optimization （BPSO） is used to find the optimum design. Combining the RBF and BPSO, a hybrid RBF-BPSO optimization method is proposed in this paper, which achieves fast optimization with high computational performance. Two examples are presented and compared to determine the optimal weight of structures under earthquake loadings using both exact and approximate analyses. The numerical results demonstrate the computational advantages and effectiveness of the proposed hybrid RBF-BPSO optimization method for the seismic design of structures.     

非线性粘滞阻尼器消能结构减振效果分析

应用并完善了非线性粘滞阻尼器消能结构地震反应预测的反应谱方法,通过与时程分析计算结果对比证明了方法的可行性。利用本方法研究了支撑刚度及阻尼器参数对非线性粘滞阻尼器减振效果的影响。通过数值分析,给出了位移降低率达到最佳时支撑刚度取值的建议式。提出了为保证剪力降低率不大于1时非线性粘滞阻尼器参数的控制方法。     

基于重启多项式Krylov子空间模型降阶的全波形瞬变电磁三维正演

瞬变电磁三维正演在内存占用和求解时间上都是大的挑战.瞬变电磁响应受到发射波形的影响, 正演计算中还需考虑发射波形.基于重启多项式Krylov子空间模型降阶方法, 本文实现了一种新的全波形瞬变电磁三维正演算法.采用规则六面体网格的拟态有限体积方法对全波形控制方程进行空间离散, 引入一阶φ函数推导得到全时段响应的统一表达式, 利用重启多项式Krylov子空间算法求解φ函数表达式, 对于给定的重启子空间维度, 基于残差公式得到满足给定精度的任意时刻的正演响应.该算法内存占用小, 不需要求解大型线性方程组, 内存占用主要由空间离散网格数量和重启子空间维度决定.数值模型正演结果验证了本文算法的有效性.

     

Development of fragility curves using high‐dimensional model representation

Fragility curves represent the conditional probability that a structure's response may exceed the performance limit for a given ground motion intensity. Conventional methods for computing building fragilities are either based on statistical extrapolation of detailed analyses on one or two specific buildings or make use of Monte Carlo simulation with these models. However, the Monte Carlo technique usually requires a relatively large number of simulations to obtain a sufficiently reliable estimate of the fragilities, and it is computationally expensive and time consuming to simulate the required thousands of time history analyses. In this paper, high‐dimensional model representation based response surface method together with the Monte Carlo simulation is used to develop the fragility curve, which is then compared with that obtained by using Latin hypercube sampling. It is used to replace the algorithmic performance‐function with an explicit functional relationship, fitting a functional approximation, thereby reducing the number of expensive numerical analyses. After the functional approximation has been made, Monte Carlo simulation is used to obtain the fragility curve of the system. Copyright © 2012 John Wiley & Sons, Ltd.     

Conversion between solid and beam element solutions of finite element method based on meta-modeling theory: development and application to a ramp tunnel structure

In this study, a new method for conversion of solid finite element solution to beam finite element solution is developed based on the meta-modeling theory which constructs a model consistent with continuum mechanics. The proposed method is rigorous and efficient compared to a typical conversion method which merely computes surface integration of solid element nodal stresses to obtain cross-sectional forces. The meta-modeling theory ensures the rigorousness of proposed method by defining a proper distance between beam element and solid element solutions in a function space of continuum mechanics. Results of numerical verification test that is conducted with a simple cantilever beam are used to find the proper distance function for this conversion. Time history analysis of the main tunnel structure of a real ramp tunnel is considered as a numerical example for the proposed conversion method. It is shown that cross-sectional forces are readily computed for solid element solution of the main tunnel structure when it is converted to a beam element solution using the proposed method. Further, envelopes of resultant forces which are of primary importance for the purpose of design, are developed for a given ground motion at the end.     

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

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

基于数值模拟和地震活动性统计的混合地震预测: 在中国地震科学实验场的应用

地震数值预测是近年来最为前沿的科学研究问题之一,但是大部分地震数值预测的探索都受困于数值模拟步长的问题.本文引入混合预测方法,探索了较短预测时长的地震数值预测方法,并在中国地震科学实验场区中进行应用.该混合预测方法主要建立在活动块体理论和构造力学分析的基础上,利用数值模拟推测未来长时间的应力变化量,并引入上一次7级地震之后数十年间,中等强度地震同震应力的影响,最后通过重分类和加权,在应力变化量上叠加年尺度小地震活动异常的统计,从而实现区域地震危险程度的混合评估.基于该方法,我们将地震数值预测的预测时间尺度缩短到年尺度,以年为单位滚动预测了中国地震科学实验场自2021—2030年或更长时间段内可能出现中强地震的区域.2022年芦山6.1级地震、2021年漾濞6.4级地震、2021年泸州6.0级地震和2022年泸定6.8级地震均发生在危险性较高的区域,而2022年马尔康6.0级地震则不在预测区域内.本文展示的混合预测尝试为地震数值预测研究提供了一个可供参考的例证,也为解决数值预测步长控制难的问题提供了一个可能的探索途径.

     

仿古建筑的隔震加固研究

为对具有传统文化的仿古建筑在不破坏其原有建筑风貌的基础上进行抗震加固,将隔震技术应用于仿古建筑。采用数值模拟的方法研究仿古建筑隔震,建立隔震与非隔震两种有限元计算分析模型,并对加固前、后模型进行动力时程分析。通过输入不同频谱的El Centro Site波、Taft波和RH1TG055波,对比分析两种模型在地震作用下的位移、加速度时程曲线。结果表明：采用隔震技术加固方案后,结构的地震响应得到明显的降低,从而显著地提高了结构的抗震能力。利用该方法加固,节约造价,施工方便,可供同类工程参考。     

微分求积法在结构动力分析中的应用

微分求积法（DQM）是1种求解微分方程初（边）值问题的数值方法,通常以较小的计算工作量即可获得较高的数值精度。这种方法应用于工程领域时多用来解决梁、板等结构的静力分析或结构特征值分析等问题,即对边值问题的微分方程的求解。结构动力分析属于初值问题,荷载和结构反应都具有特殊性,直接套用DQM求解边值问题并不能获得问题的解。本文尝试利用微分求积原理建立求解结构动力反应的具体方法。借鉴单元法的思想,将荷载持时划分为若干个时步,在每个时步内对动态荷载和结构反应进行离散,然后用DQM对时步逐个进行求解,得到体系在整个时域内的反应过程。通过对3种不同自振周期的线弹性单自由度体系在不同频率简谐激励下反应的计算,阐释了本文方法的可行性以及高精度、高效率的特点,通过数值试验确定了时步内相对较优的节点数,并为时步长度的选取提供了建议。     

A new analytical approach to reconstruct the acceleration time history at the bedrock base from the free surface signal records

Acceleration time histories of earthquake events are typically measured in seismic stations that are placed close to the soil top surface. These acceleration records are often used as input data for seismic analysis. It may be used for base excitation in seismic analysis of above ground structures with shallow foundations.. However it may not be used for seismic analysis of underground structures, or even for above ground buildings with deep foundations and several underground stories. The required base excitation data of the latter should have been measured below the top surface, at a level that may be determined according to the specific analyzed building geometry or at the bedrock below. If the acceleration time history at the bedrock would have been known, the seismic wave propagation through the soil medium, from the bedrock towards the top surface, could have been carried out and the base excitation of the buried structure could be determined. Since there is no data on the acceleration time history at the bedrock, and the only given data is the acceleration records at the top surface, the goal of this paper is to provide an exact reverse analysis procedure to determine the unknown acceleration time history at the bedrock that would exactly produce the measured acceleration time history at the top surface. Once this goal is achieved, seismic analysis of buried structures may be carried out with the determined acceleration record at the bedrock as input. This paper presents an analytical exact solution of the inverse problem for determination of the acceleration, velocity and displacement time histories at the bedrock base of a layered geological medium that are compatible with the given acceleration record at the soil top surface. This new proposed method is based on analytical solutions of the initial-boundary value problems of the linear wave equation in the case of a layered medium. The relationship between waves in one layer and waves in another adjacent layer is derived considering the continuity of stresses and displacements at the common interface between the layers. The efficiency and accuracy of the proposed method is demonstrated through several examples involving the nonstationary response of the free surface. The case of the San Fernando Earthquake is studied. Excellent agreement is achieved between the recorded free surface time history and the reconstructed signal. This excellent agreement is obtained due to the exact analytical method used in deriving the inverse problem solution. This exact analytical method allows one to obtain an acceleration (velocity/displacement) distribution along all the layers at any time.     

Homotopy perturbation method for solving the space–time fractional advection–dispersion equation

In this paper we present a reliable algorithm, the homotopy perturbation method, to construct numerical solutions of the space–time fractional advection–dispersion equation in the form of a rapidly convergent series with easily computable components. Fractional advection–dispersion equations are used in groundwater hydrology to model the transport of passive tracers carried by fluid flow in a porous medium. The fractional derivatives are described in the Caputo sense. Some examples are given. Numerical results show that the homotopy perturbation method is easy to implement and accurate when applied to space–time fractional advection–dispersion equations.     

Displacement amplification factors for steel eccentrically braced frames

Inelastic deformation capacity of links is a factor that significantly influences design of steel eccentrically braced frames (EBFs). The link rotation angle is used to describe inelastic link deformation. The link rotation angle is generally calculated by making use of design story drifts that in turn are calculated by modifying the elastic displacements by a displacement amplification factor. This paper presents a numerical study undertaken to evaluate the displacement amplification factor given in ASCE7‐10 for EBFs and the rigid‐plastic mechanism used for calculating link rotation angles. A total of 72 EBFs were designed by considering the number of stories, the bay width, the link length to bay width ratio, and the seismic hazard level as the prime variables. All structures were analyzed using elastic and inelastic time history analyses. The results indicated that the displacement amplification factor given in ASCE7‐10 provides unconservative estimates of the story drifts. On the other hand, the rigid‐plastic mechanism provides conservative estimates of link rotations. Based on the results of the numerical study, a new set of displacement amplification factors that vary along the height of the structure and modifications to the rigid‐plastic mechanism were developed. In light of the proposed modifications, the EBFs were redesigned and analyzed using inelastic time history analysis. The results indicated that the proposed modifications provide improvements for the displacement amplification factor and link rotation angle calculation procedures. Copyright © 2014 John Wiley & Sons, Ltd.     

Assessment of errors in the newmark method in structural dynamics

In the Newmark and other approximate step-by-step methods, having introduced assumptions in order to transform the differential equations, which are characteristic of response problems, into simultanéous equations, successive solutions lead to a response-time history. In this paper numerical results and formulae are given for the errors which are generated by this procedure. These errors are oscillatory in nature and, in general, the oscillations increase in magnitude as the number of time steps increases. Recommendations for upper limits on the time step, which will provide acceptable accuracy for a wide range of system and excitation parameters, are presented.     

一种显隐式积分方法及其应用

在实际工程结构动力反应分析中,往往由于结构型式十分复杂,常用的两种直接积分方法,即显式积分方法和隐式积分方法,在使用中都存在着一定的局限性,如何将这两种积分方法合理有效地结合起来,是一个十分有意义的研究课题。针对实际工程问题中整体结构计算时间步长的选择往往受局部区域的材料特性、尺寸大小等因素影响的这一现象,提出了一种对结构局部区域进行隐式积分、对其余区域进行显式积分的显隐式积分方法,这种积分格式相对于显式积分格式而言,能显著提高整体结构的计算速度。最后采用两个数值计算实例对这一方法进行验证。     

Extracting predictable components and forecasting techniques in extended-range numerical weather prediction

This paper refers to the CNOP-related algorithms and formulates the practical method and forecast techniques of extracting predictable components in a numerical model for predictable components on extended-range scales.Model variables are divided into predictable components and unpredictable chaotic components from the angle of model prediction error growth.The predictable components are defined as those with a slow error growth at a given range.A targeted numerical model for predictable components is established based on the operational dynamical extended-range forecast(DERF)model of the National Climate Center.At the same time,useful information in historical data are combined to find the fields for predictable components in the numerical model that are similar to those for the predictable components in historical data,reducing the variable dimensions in a similar judgment process and further correcting prediction errors of predictable components.Historical data is used to obtain the expected value and variance of the chaotic components through the ensemble forecast method.The numerical experiment results show that this method can effectively improve the forecast skill of the atmospheric circulation field in the 10–30 days extended-range numerical model and has good prospects for operational applications.     

Modified precise time step integration method of structural dynamic analysis

The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.     

R.C.框架结构的振动台试验和面向设计的时程分析方法

钢筋混凝土框架结构的弹塑性时程反应分析是抗震设计规范所要求的一项重要内容,但目前还滑有比较简单实用肯计算可靠的方法。本文进行了两个六层钢筋混凝土框架结构模型的地震模拟震动台试验,通过分析试验模型的自振特性、变形形式、加速度反应、位移反应和破坏开矿以及钢筋混凝土框架结构的实际震害,建立了面向设计的钢筋混凝土框架结构弹塑性里程反应分析的计算方法。在这一方法中,对钢筋混凝土框架整体结构采用层间剪切模型,