动力反应分析的显式积分方法及其稳定性

Newm ark-更新精细积分法是动力方程求解的隐式的时域逐步积分法,其稳定性条件非常容易满足。与隐式方法相比较,显式积分方法不需要求解耦联的方程组,可以有效地减少内存占用和机时耗费。因此,根据显式积分方法的特点和优点,基于Newm ark-更新精细积分法的基本思想,提出其显式积分格式。对显式积分方法的精度与稳定性进行了初步的分析,指出该显式积分方法具有极好的稳定性,其精度比隐式积分方法的精度稍低。随着时间步长的增加,其精度优于传统的方法。     

基于极大初始时间步法的非线性静动力联合分析

静动力联合分析是考虑非线性条件下结构响应分析的关键技术问题之一。本文通过分析非线性有限元静、动力平衡方程的内在联系,并结合Newmark积分算法的隐式求解特点,提出了一种基于极大初始时间步法的非线性静动力耦合算法,可保证结构在静动荷载的耦合激励下获得合理的非线性地震响应。通过对Koyna重力坝的地震响应数值分析算例,验证了该算法的合理性与有效性。另外,该算法在各类隐式求解的通用有限元程序中均易于实现,可操作性强,便于工程应用。     

显式分析方法在高层建筑弹塑性地震反应分析中的适用性研究

目前隐式方法是动力弹塑性时程分析最常用的分析方法,然而隐式方法在强非线性分析中常常存在迭代不收敛的问题,并且刚度矩阵求解消耗的存储空间随结构自由度增加呈几何级数增长。因此,在求解高层建筑这种大规模问题中,极易遭遇计算瓶颈。显式分析方法直接求解解耦的方程组,不需要迭代。本文对隐式方法和显式方法进行了对比分析,研究了显式分析方法在高层建筑弹塑性地震反应分析中适用性。实例分析表明,从计算精度来讲,隐式方法和显式方法在稳定条件下都能得到较好的精度。从计算效率来讲,对于自由度较少的结构,隐式方法计算效率较高;对于自由度庞大的结构,显式方法计算效率较高。建议在进行自由度多的高层建筑弹塑性地震反应分析时,采用显式分析方法。     

求解振动方程的一种显式积分格式及其精度与稳定性

介绍了一种求解有阻尼体系振动方程的显式积分格式及其逐步求解的过程，并对其计算精度和稳定性进行了分析。该方法不但能同时求得体系的位移、速度和加速度反应，而且所得到的加速度反应的精度能满足工程需要。     

一种全局优化的隐式交错网格有限差分算法及其在弹性波数值模拟中的应用

在数值模拟中,隐式有限差分具有较高的精度和稳定性.然而,传统隐式有限差分算法大多由于需要求解大型矩阵方程而存在计算效率偏低的局限性.本文针对一阶速度-应力弹性波方程,构建了一种优化隐式交错网格有限差分格式,然后将改进格式由时间-空间域转换为时间-波数域,利用二范数原理建立目标函数,再利用模拟退火法求取优化系数.通过对均匀模型以及复杂介质模型进行一阶速度-应力弹性波方程数值模拟所得单炮记录、波场快照分析表明:这种优化隐式交错网格差分算法与传统的几种显式和隐式交错网格有限差分算法相比不但降低了计算量,而且能有效的压制网格频散,使弹性波数值模拟的精度得到有效的提高.     

结构动力学方程的显式积分格式

本文从空间解耦有限元常微分方程组出发，探讨了结构动力学方程的高精度显式积分格式。通过被积函数的拉格朗日多项式内插和分部积分导出了波动数值模拟的一组显式时步积分公式。这组公式是时间和空间解耦的，即波场内任一离散节点在任一时刻的波动数据可以用这组公式依据该节点及其邻近节点在该时刻之前的n＋1个时刻的波动数据显式地算出（n为非负整数），阐明了这组公式的如下特点：第一，其截断误差的量级不超过0（Δt^n＋3），Δt为时间步距。第二，它不仅可用于线性波动的数值模拟，而且可用于本构方程具有强非线性情形。第三，这组公式也可推广应用于一系列数学物理暂态问题的数值求解。针对一个简单的时不变系统初步分析了此组积分格式的稳定性。但是，对其稳定性尚需作进一步研究。     

粘弹性场地地形对地震动影响分析的显式有限元-有限差分方法

推导出了分析二维粘及弹性场地地形对地震动影响的显式有限元-有限差分方法．这一方法中，首先利用人工边界及有限元离散方法，给出问题分析的有限元离散网格计算力学模型，并利用一种类似于差分方法的有限元方法，建立局部网格节点的动力方程，而后利用笔者提出的有阻尼体系动力方程求解的显式差分格式，及推广的多次透射边界公式，给出网格节点运动量计算的时域显式逐步积分公式．利用计算机程序实现这一方法的计算具有所需计算机内存量小及计算时间量小的优势，而且，这一方法适用于任意地形情况，具有较高的计算精度及较好的计算稳定性．      

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

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

用积分方程法模拟各向异性地层中三维电性异常体的电磁响应

本文研究并建立了一种模拟各向异性地层中三维电性异常体电磁响应的积分方程算法.首先讨论了并矢Green函数及其相关积分的计算,将水平层状各向异性地层中的电场并矢Green函数分解成含有奇异项的直达波与非奇异的来自各个层界面的反射和透射波两个部分,再应用等效体积单元和表面积分技术对积分方程的奇异核进行离散化处理以便提高离散方程的精度.然后为了节省计算机内存以及计算时间,引入基于Krylov子空间的迭代算法求解积分方程的离散化矩阵方程.最后通过与现有文献中的结果作对比从而检验了所述算法的有效性,并结合具体算例考察分析了地层的各向异性对三维电性异常体电磁响应的影响特征和规律.     

地震波动方程的局部间断有限元方法数值模拟

近年来,随着地震波数值模拟对计算精度和效率的要求越来越高,间断有限元方法开始受到越来越多的关注.本文中,针对具有吸收边界条件的二维地震声波波动方程,作者提出了一种基于局部间断有限元方法的数值模拟算法.该算法在空间上使用局部间断有限元方法进行离散,在时间上采用了显式蛙跳格式.在这种时空离散的组合方式下,每个时间步上,此算法在空间剖分的每个单元上的求解计算是相互独立的,因而具有极高的并行性.通过数值算例,我们将该算法与连续有限元方法进行了比较.结果表明,本算法不仅具有对起伏构造的良好适应性,而且在计算效率和计算精度等方面,都具有优越性.     

Stability analysis of SDOF real‐time hybrid testing systems with explicit integration algorithms and actuator delay

Real‐time hybrid testing is a method that combines experimental substructure(s) representing component(s) of a structure with a numerical model of the remaining part of the structure. These substructures are combined with the integration algorithm for the test and the servo‐hydraulic actuator to form the real‐time hybrid testing system. The inherent dynamics of the servo‐hydraulic actuator used in real‐time hybrid testing will give rise to a time delay, which may result in a degradation of accuracy of the test, and possibly render the system to become unstable. To acquire a better understanding of the stability of a real‐time hybrid test with actuator delay, a stability analysis procedure for single‐degree‐of‐freedom structures is presented that includes both the actuator delay and an explicit integration algorithm. The actuator delay is modeled by a discrete transfer function and combined with a discrete transfer function representing the integration algorithm to form a closed‐loop transfer function for the real‐time hybrid testing system. The stability of the system is investigated by examining the poles of the closed‐loop transfer function. The effect of actuator delay on the stability of a real‐time hybrid test is shown to be dependent on the structural parameters as well as the form of the integration algorithm. The stability analysis results can have a significant difference compared with the solution from the delay differential equation, thereby illustrating the need to include the integration algorithm in the stability analysis of a real‐time hybrid testing system. Copyright © 2007 John Wiley & Sons, Ltd.     

Improving explicit time integration by modal truncation techniques

This paper describes a modal weighting technique that improves the stability characteristics of explicit time-integration schemes used in structural dynamics. The central difference method was chosen as the trial algorithm because of its simplicity, both in terms of formulation and ease of numerical stability and convergence analysis. It is shown how explicit algorithms may be reformulated in order to make them stable for any integration time by attenuating high-frequency oscillation modes that are generated by mesh geometry rather than generic dynamical features. We discuss results from trial calculations obtained from mathematical models that represent hysteretic restoring force elements and an application on a physical, four-degree-of-freedom, base-isolated structure using the pseudodynamic technique. © 1998 John Wiley & Sons, Ltd.     

求解运动方程的一种不等时间步长的显式数值积分方法

在波动有限元模拟中, 若采用传统的显式数值积分方法求解运动方程, 计算时间步长需采用计算区内满足稳定条件要求的最小时间步长. 然而, 对于大部分计算区域, 这一时间步长过小, 是不必要的. 本文提出了一种不等时间步长的显式数值积分方法, 其基本思想是不同的计算区域采用满足各自稳定条件的计算时间步长. 最后, 本文通过数值试验检验了这一方法的可行性及其对数值计算精度的影响.      

Real‐time hybrid testing using the unconditionally stable explicit CR integration algorithm

Real‐time hybrid testing combines experimental testing and numerical simulation, and provides a viable alternative for the dynamic testing of structural systems. An integration algorithm is used in real‐time hybrid testing to compute the structural response based on feedback restoring forces from experimental and analytical substructures. Explicit integration algorithms are usually preferred over implicit algorithms as they do not require iteration and are therefore computationally efficient. The time step size for explicit integration algorithms, which are typically conditionally stable, can be extremely small in order to avoid numerical stability when the number of degree‐of‐freedom of the structure becomes large. This paper presents the implementation and application of a newly developed unconditionally stable explicit integration algorithm for real‐time hybrid testing. The development of the integration algorithm is briefly reviewed. An extrapolation procedure is introduced in the implementation of the algorithm for real‐time testing to ensure the continuous movement of the servo‐hydraulic actuator. The stability of the implemented integration algorithm is investigated using control theory. Real‐time hybrid test results of single‐degree‐of‐freedom and multi‐degree‐of‐freedom structures with a passive elastomeric damper subjected to earthquake ground motion are presented. The explicit integration algorithm is shown to enable the exceptional real‐time hybrid test results to be achieved. Copyright © 2008 John Wiley & Sons, Ltd.     

Implementation and application of the unconditionally stable explicit parametrically dissipative KR‐α method for real‐time hybrid simulation

In real‐time hybrid simulations (RTHS) that utilize explicit integration algorithms, the inherent damping in the analytical substructure is generally defined using mass and initial stiffness proportional damping. This type of damping model is known to produce inaccurate results when the structure undergoes significant inelastic deformations. To alleviate the problem, a form of a nonproportional damping model often used in numerical simulations involving implicit integration algorithms can be considered. This type of damping model, however, when used with explicit integration algorithms can require a small time step to achieve the desired accuracy in an RTHS involving a structure with a large number of degrees of freedom. Restrictions on the minimum time step exist in an RTHS that are associated with the computational demand. Integrating the equations of motion for an RTHS with too large of a time step can result in spurious high‐frequency oscillations in the member forces for elements of the structural model that undergo inelastic deformations. The problem is circumvented by introducing the parametrically controllable numerical energy dissipation available in the recently developed unconditionally stable explicit KR‐α method. This paper reviews the formulation of the KR‐α method and presents an efficient implementation for RTHS. Using the method, RTHS of a three‐story 0.6‐scale prototype steel building with nonlinear elastomeric dampers are conducted with a ground motion scaled to the design basis and maximum considered earthquake hazard levels. The results show that controllable numerical energy dissipation can significantly eliminate spurious participation of higher modes and produce exceptional RTHS results. Copyright © 2014 John Wiley & Sons, Ltd.     

A discontinuous finite element baroclinic marine model on unstructured prismatic meshes

We describe the time discretization of a three-dimensional baroclinic finite element model for the hydrostatic Boussinesq equations based upon a discontinuous Galerkin finite element method. On one hand, the time marching algorithm is based on an efficient mode splitting. To ensure compatibility between the barotropic and baroclinic modes in the splitting algorithm, we introduce Lagrange multipliers in the discrete formulation. On the other hand, the use of implicit–explicit Runge–Kutta methods enables us to treat stiff linear operators implicitly, while the rest of the nonlinear dynamics is treated explicitly. By way of illustration, the time evolution of the flow over a tall isolated seamount on the sphere is simulated. The seamount height is 90% of the mean sea depth. Vortex shedding and Taylor caps are observed. The simulation compares well with results published by other authors.     

Energy-consistent integration method and its application to hybrid testing

The time-integration algorithm is an indispensable element to determine response of the boundary of the numerical as well as physical parts in a hybrid test. Instability of the time-integration algorithm may directly lead to failure of the test, so stability of an integration algorithm is particularly important for hybrid testing. The explicit algorithms are very popular in hybrid testing, because iteration is not needed. Many unconditionally stable explicit-algorithms have been proposed for hybrid testing. However, the stability analysis approaches used in all these methods are valid only for linear systems. In this paper, a uniform formulation for energy-consistent time integrations, which are unconditionally stable, is proposed for nonlinear systems. The solvability and accuracy are analyzed for typical energy-consistent algorithms. Some numerical examples and the results of a hybrid test are provided to validate the effectiveness of energy-consistent algorithms.     

Extensions of nonlinear error propagation analysis for explicit pseudodynamic testing

Two important extensions of a technique to perform a nonlinear error propagation analysis for an explicit pseudodynamic algorithm (Chang, 2003) are presented. One extends the stability study from a given time step to a complete step-by-step integration procedure. It is analytically proven that ensuring stability conditions in each time step leads to a stable computation of the entire step-by-step integration procedure. The other extension shows that the nonlinear error propagation results, which are derived for a nonlinear single degree of freedom (SDOF) system, can be applied to a nonlinear multiple degree of freedom (MDOF) system. This application is dependent upon the determination of the natural frequencies of the system in each time step, since all the numerical properties and error propagation properties in the time step are closely related to these frequencies. The results are derived from the step degree of nonlinearity. An instantaneous degree of nonlinearity is introduced to replace the step degree of nonlinearity and is shown to be easier to use in practice. The extensions can be also applied to the results derived from a SDOF system based on the instantaneous degree of nonlinearity, and hence a time step might be appropriately chosen to perform a pseudodynamic test prior to testing.