动力方程求解的显式积分格式及其稳定性与适用性

文献（１）给出了一种求解有阻尼体系动力方程的显式积分格式,文中以数值表格的形式给出了格式的稳定性条件,本文对该格式的稳定性问题作了进一步的分析,并给出了其计算稳定性条件的表示式。本文还着重讨论了基于这一显式积分格式的推导过程而派生出的另一形式的积分格式的稳定性,并指出了该派生格式的适用性问题。     

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

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

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

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

结构动力分析的非线性拟动力方程法

本文在时间步长△t内用二次曲线来拟合动力反应的速度和加速度时程曲线，并且得到了求解结构动力响应的拟动力方程。由于可直接由动力方程求解动力反应，因而物理概念明确，可克服计算中误差积累问题，并且与线性法相比，在保证计算精度的条件下，可减少时间步长的划分，提高运算速度。     

求解动力方程的三种显式积分格式的比较

比较了工程上常用于求解动力方程的三种显式积分格式的能耗特性,结果表明,李小军等人的显式差分格式在物理阻尼比为0到1的范围内均具有较强的能耗特性,而且随着物理阻尼比的增大,能耗特性也增强,这种特性优于另外两种显式积分格式,能有效地克服透射边界的高频失稳,因此在处理复杂场地地震反应方面是一种比较理想的显式积分格式。     

阻尼振动方程的一种显式直接积分方法

介绍了一种求解阻尼振动方程的条件稳定的直接积分方法及其逐步求解过程,该方法是显式的,适于分析阻尼系统的动力反应。     

有阻尼振动方程常用显式积分格式稳定性分析

本文主要介绍了用于求解有阻尼振动方程的四种常用显式积分方法，并针对其稳定性和精度进行了分析对比，讨论了其适用范围及阻尼对稳定性的影响。     

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

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

波动有限元方程显式逐步积分格式稳定性分析

对常用的波动有限元方程的两种显式格式的稳定性进行了分析,利用单元的最大频率乃是系统的特征频率的上界的概念结合模态分析方法,给出了便于实际应用的稳定性必要条件。同时,利用Von Neumann方法给出了这两种格式稳定性的充分条件,并通过算例对这两种稳定性条件进行验证。     

动力边界积分方程法及其应用

文中简要介绍了作者在国家地震局工程力学研究所学习和工作期间，该所动力边界积分方程法及其应用课题研究组这几年来的研究成果。     

Precise integration method without inverse matrix calculation for structural dynamic equations

The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.     

Explicit concomitance of implicit method to solve vibration equation

It has been proven that the implicit method used to solve the vibration equation can be transformed into an explicit method,which is called the concomitant explicit method.The constant acceleration met...     

Comparison of capability of time integration methods in capturing dynamic loading

Numerical properties of the time integration method proposed by the first author of this paper in 2007 are the same as those of the constant average acceleration method (AAM) for linear elastic systems, except that the capability to capture dynamic loading was not explored. It was found that there were different quadrature equations to predict the next step displacement increment. A modified quadrature equation of this method was derived so that the equation to determine the next step displacement was numerically equivalent to the equation used in the constant AAM. It was verified that the original form of this method, in general, had a better capability to capture dynamic loadings than the constant AAM. This excellent property, in addition to computational efficiency, will help to make this method competitive with general secondorder accurate integration methods.     

Stability of average acceleration method for structures with nonlinear damping

The energy approach is used to theoretically verify that the average acceleration method （AAM）, which is unconditionally stable for linear dynamic systems, is also unconditionally stable for structures with typical nonlinear damping, including the special case of velocity power type damping with a bilinear restoring force model. Based on the energy approach, the stability of the AAM is proven for SDOF structures using the mathematical features of the velocity power function and for MDOF structures by applying the virtual displacement theorem. Finally, numerical examples are given to demonstrate the accuracy of the theoretical analysis.     

散射波积分方程的Adomian分解解法

在背景模型基础上,求解模型扰动后的地震波散射场,这是目前地震反演中的一个关键步骤.本文将计算数学中求解非线性积分方程的Adomian分解方法,应用到求解标量波散射场的Lippmann-Schwinger积分方程和Ricatti积分方程中,分别得到了散射场的Born序列解和Rytov序列解.通过一维和二维数值算例说明：在满足一定的条件下,散射场的这两种序列解稳定收敛,与传统的Born和Rytov近似解相比,引入散射序列中的高阶项可以更精确地描述地震波散射场.

     

Analytical power series solutions to the two-dimensional advection–dispersion equation with distance-dependent dispersivities

As is frequently cited, dispersivity increases with solute travel distance in the subsurface. This behaviour has been attributed to the inherent spatial variation of the pore water velocity in geological porous media. Analytically solving the advection–dispersion equation with distance-dependent dispersivity is extremely difficult because the governing equation coefficients are dependent upon the distance variable. This study presents an analytical technique to solve a two-dimensional (2D) advection–dispersion equation with linear distance-dependent longitudinal and transverse dispersivities for describing solute transport in a uniform flow field. The analytical approach is developed by applying the extended power series method coupled with the Laplace and finite Fourier cosine transforms. The developed solution is then compared to the corresponding numerical solution to assess its accuracy and robustness. The results demonstrate that the breakthrough curves at different spatial locations obtained from the power series solution show good agreement with those obtained from the numerical solution. However, owing to the limited numerical operation for large values of the power series functions, the developed analytical solution can only be numerically evaluated when the values of longitudinal dispersivity/distance ratio e_L exceed 0·075. Moreover, breakthrough curves obtained from the distance-dependent solution are compared with those from the constant dispersivity solution to investigate the relationship between the transport parameters. Our numerical experiments demonstrate that a previously derived relationship is invalid for large e_L values. The analytical power series solution derived in this study is efficient and can be a useful tool for future studies in the field of 2D and distance-dependent dispersive transport. Copyright © 2008 John Wiley & Sons, Ltd.     

结构动力分析阻尼模型研究

提出了结构动力时程分析修正的Ｒａｙｌｅｉｇｈ阻尼模型和一个基于单元特性的阻尼参数方法,可用于合理确定复杂阻尼结构的阻尼参数,工程应用实例表明,阻尼参数的计算对时程响应计算结果的影响是十分明显的,必须在分析结构体系振动特性的基础上慎重确定阻尼参数,以避免为为放大或滤掉主要振型,造成虚假的结论。     

用于结构健康诊断的动静结合法及其试验验证

介绍了一种用于梁式结构的健康诊断方法,该方法以具有较高测试精度的低阶模态参数和静力测试参数为损伤识别量,适用于本构相同的多个结构或构件的健康诊断.具体实施中只选取少量结构作为样本,对其进行静、动力的联合测试,以获得静力、动力参数间的数学回归关系,依此关系对其余的结构仅通过动力测试便可获得等效于静力加载法的健康诊断效果.以12个不同损伤程度的钢桁架结构为试验对象,对其静、动参数的相关性进行了研究.结果表明,桁架的静力刚度和一阶固有频率与桁架的损伤程度密切相关,二者可作为标识桁架损伤的静、动参数.静、动参数对桁架的损伤程度指示明确,规律是损伤越大,静、动参数值越小.采用最小二乘法建立了桁架静、动参数间数学关系式,说明梁式结构的静、动参数密切相关,动静结合法可行.