微分求积时间单元方法

基于微分求积法则,提出了一种求解动力学常微分方程的高效高精度微分求积时间单元方法(DQTEM).给出了DQTEM施加初始位移和初始速度的方法,其结果相当于构造了C1时间单元.与递推格式的直接积分方法不同,对于考虑的时间域通常只需用一个微分求积时间单元.与RK法和Newmark法相比,用少量时间结点的DQTEM结果就与精确解吻合.稳定性分析表明,DQTEM通常是条件稳定的.     

A unified set of single step algorithms part 3: The beta-m method,a generalization of the Newmark scheme

Introduced herein is a generalization of Newmark's time marching integration scheme, called the β-m method. Like the SSpj method (introduced in Parts 1 and 2 of this series), the β-m method provides a gcneral single-step scheme applicable to any set of initial value problems. The method is specialized by specifying the method order m along with rn integration parameters, β0, β1, …, β_m?1. For a particular choice of m, the integration parameters provide a subfamily of methods which control accuracy and stability, as well as offering options for explicit or implicit algorithms. For the most part, attention is focused on the application to structural dynamic equations. Most well-known methods (e.g. Newmark, Wilson, Houbolt, etc.) are shown to be special cases within the β-m family. Hence, one computationally efficient and surprisingly simple algorithm unifies old and new methods. Stability and accuracy analyses are presented for method orders m = 2, 3 and 4 to determine optimal parameters for implicit and explicit schemes, along with numerical verification.     

On the equivalence of the time domain differential quadrature method and the dissipative Runge–Kutta collocation method

Numerical solutions for initial value problems can be evaluated accurately and efficiently by the differential quadrature method. Unconditionally stable higher order accurate time step integration algorithms can be constructed systematically from this framework. It has been observed that highly accurate numerical results can also be obtained for non‐linear problems. In this paper, it is shown that the algorithms are in fact related to the well‐established implicit Runge–Kutta methods. Through this relation, new implicit Runge–Kutta methods with controllable numerical dissipation are derived. Among them, the non‐dissipative and asymptotically annihilating algorithms correspond to the Gauss methods and the Radau IIA methods, respectively. Other dissipative algorithms between these two extreme cases are shown to be B‐stable (or algebraically stable) as well and the order of accuracy is the same as the corresponding Radau IIA method. Through the equivalence, it can be inferred that the differential quadrature method also enjoys the same stability and accuracy properties. Copyright © 2001 John Wiley & Sons, Ltd.     

Effects of time integration schemes on discontinuous deformation simulations using the numerical manifold method

Numerical stability by using certain time integration scheme is a critical issue for accurate simulation of discontinuous deformations of solids. To investigate the effects of the time integration schemes on the numerical stability of the numerical manifold method, the implicit time integration schemes, ie, the Newmark, the HHT‐α, and the WBZ‐α methods, and the explicit time integration algorithms, ie, the central difference, the Zhai's, and Chung‐Lee methods, are implemented. Their performance is examined by conducting transient response analysis of an elastic strip subjected to constant loading, impact analysis of an elastic rod with an initial velocity, and excavation analysis of jointed rock masses, respectively. Parametric studies using different time steps are conducted for different time integration algorithms, and the convergence efficiency of the open‐close iterations for the contact problems is also investigated. It is proved that the Hilber‐Hughes‐Taylor‐α (HHT‐α), Wood‐Bossak‐Zienkiewicz‐α (WBZ‐α), Zhai's, and Chung‐Lee methods are more attractive in solving discontinuous deformation problems involving nonlinear contacts. It is also found that the examined explicit algorithms showed higher computational efficiency compared to those implicit algorithms within acceptable computational accuracy.     

动力学方程的解析逐步积分法

提出了求解动力学方程的一个新型的逐步积分法。基于动力学方程的解析齐次解,构造出动力学方程解的一般积分表达式,借助于显式、自起动、预测－校正的单步四阶精度的积分算法,离散方程右端的等价荷载项,给出了一个新的解析逐步积分方法格式。如果用分块求解,其刚度阵、质量阵等将有较小的规模,将使计算效率更高。算例表明本文方法比中心差分法、Ｎｅｗｍａｒｋ、Ｗｉｌｓｏｎ－θ、Ｈｏｕｂｏｌｔ法等有较高的精度,本文结果更接近解析解。本文方法也适用于非线性,因为本计算格式是显示,因此不需要迭代求解。     

基于微分求积法及V-变换的大规模动力系统快速数值计算方法

针对大规模动力系统动态响应的数值计算,传统的微分求积法通常在时间域上逐步离散、整体求解,存在“维数灾”问题。在多级高阶时域微分求积法的基础上,提出了基于V-变换的大规模动力系统动态响应的快速数值计算方法。利用微分求积法的加权系数矩阵满足V-变换这一重要特性,将离散后的雅可比矩阵方程进行解耦分块,推导形成了多级分块递推计算方法。数值算例表明,即使采用相当于Newmark方法2s倍的步长,微分求积法的计算精度仍比Newmark方法要高出2~3个数量级。进一步对3个不同规模的算例系统进行了测试,结果表明：相对于传统的数值计算方法,多级分块递推计算方法可以获得较大的加速比,能够显著提高大规模动力系统动态响应的计算效率。     

A coupled numerical scheme of dual reciprocity BEM with DQM for the transient elastodynamic problems

The two‐dimensional transient elastodynamic problems are solved numerically by using the coupling of the dual reciprocity boundary element method (DRBEM) in spatial domain with the differential quadrature method (DQM) in time domain. The DRBEM with the fundamental solution of the Laplace equation transforms the domain integrals into the boundary integrals that contain the first‐ and the second‐order time derivative terms. Thus, the application of DRBEM to elastodynamic problems results in a system of second‐order ordinary differential equations in time. This system is then discretized by the polynomial‐based DQM with respect to time, which gives a system of linear algebraic equations after the imposition of both the boundary and the initial conditions. Therefore, the solution is obtained at any required time level at one stroke without the use of an iterative scheme and without the need of very small step size in time direction. The numerical results are visualized in terms of graphics. Copyright © 2008 John Wiley & Sons, Ltd.     

Time‐domain integration methods of exponentially damped linear systems

In this paper, three new kinds of time‐domain numerical methods of exponentially damped systems are presented, namely, the simplified Newmark integration method, the precise integration method, and the simplified complex mode superposition method. Based on the traditional Newmark integration method and transforming the equation of motion with exponentially damping kernel functions into an equivalent second‐order equation of motion by using the internal variables technique, the simplified Newmark integration method is developed by using a decoupling technique to reduce the computer run time and storage. By transforming the equation of motion with exponentially damping kernel functions into a first‐order state‐space equation, the precise integration technique is used to numerically solve the state‐space equation. Based on a symmetric state‐space equation and the complex mode superposition method, a delicate and simplified general solution of exponentially damped linear systems, completely in real‐value form, is developed. The accuracy and efficiency of the developed numerical methods are compared and discussed by two benchmark examples.     

基于微分求积的结构随机振动时域分析方法

将微分求积法应用于结构动力学方程的逐步时程积分时存在计算效率低的问题.为此,从数值积分角度出发,采用复化微分求积公式计算Duhamel积分项,并将其和精细积分法结合,可形成一种计算任意随机激励下结构随机振动时域分析的显式求解方法.该方法无需对系数矩阵求逆,能够减小在一个积分步长内载荷量线性化所造成的误差,同时也提高数值...     

A simple explicit single step time integration algorithm for structural dynamics

In this article, a new single-step explicit time integration method is developed based on the Newmark approximations for the analysis of various dynamic problems. The newly proposed method is second-order accurate and able to control numerical dissipation through the parameters of the Newmark approximations. Explicitness and order of accuracy of the proposed method are not affected in velocity-dependent problems. Illustrative linear and nonlinear examples are used to verify performances of the proposed method.     

一种多自由度阻尼体系动力问题的高阶分析方法

传统动力时程直接积分法多采用低阶数值格式,需要选择非常小的时间步距才能获得满足精度要求的动力分析结果.该文将结构动力时程分析的积分求微法推广至多自由度情形,发展了一种具有较高计算效率的多自由度阻尼体系的动力时程高阶分析方法.将相邻的ρ个时步组成一个待求解时段,基于多自由度体系动力响应积分解,以精细积分法结合秦九韶算法计...     

任意阶显式精细积分多步法在刚性方程中的应用研究

将作者提出的高精度任意阶显式精细积分多步法应用于刚性方程中。本文的方法可方便地进行不同阶次的运算。将本文的方法与精确值和其它数值计算方法进行比较,数值计算结果表明本文方法是一种高精度、高效率的方法。     

On the accuracy of direct integration methods for structural dynamics

Abstract

This paper presents a method of eigen analysis to evaluate the accuracy of the time integration method in the structural dynamic analysis of both transient and steady state responses. Accuracy measure for transient response is evaluated by period elongation and amplitude decay after a complete cycle of response of an undamped system. The proposed method is a unique measure, independent of the initial conditions. Accuracy measure for steady state response is evaluated by the amplitude decay and phase angle of a frequency response. The accuracy measure proposed by this method can be applied to all ranges of ?t/T. The significance of the large range ?t/T is discussed. The proposed method is used to evaluate the accuracy of three commonly used time integration methods: the Wilson method, the Newmark method and the Houbolt method.     

An improved non-classical method for the solution of fractional differential equations

A procedure to construct temporally local schemes for the computation of fractional derivatives is proposed. The frequency-domain counterpart (iω) ^α of the fractional differential operator of order α is expressed as an improper integral of a rational function in iω. After applying a quadrature rule, the improper integral is approximated by a series of partial fractions. Each term of the partial fractions corresponds to an exponential kernel in the time domain. The convolution integral in a fractional derivative can be evaluated recursively leading to a local scheme. As the arguments of the exponential functions are always real and negative, the scheme is stable. The present procedure provides a convenient way to evaluate the quality of a given algorithm by examining its accuracy in fitting the function (iω) ^α . It is revealed that the non-classical solution methods for fractional differential equations proposed by Yuan and Agrawal (ASME J Vib Acoust 124:321–324, 2002) and by Diethelm (Numer Algorithms 47:361–390, 2008) can also be interpreted as applying specific quadrature rules to evaluate the improper integral numerically. Over a wider range of frequencies, Diethelm’s algorithm provides a more accurate fitting than the YA algorithm. Therefore, it leads to better performance. Further exploiting this advantage of the proposed derivation, a novel quadrature rule leading to an even better performance than Diethelm’s algorithm is proposed. Significant gains in accuracy are achieved at the extreme ends of the frequency range. This results in significant improvements in accuracy for late time responses. Several numerical examples, including fractional differential equations of degree α = 0.3 and α = 1.5, demonstrate the accuracy and efficiency of the proposed method.     

On the elastodynamic solution of frictional contact problems using variational inequalities

This article is concerned with the development, implementation and application of variational inequalities to treat the general elastodynamic contact problem. The solution strategy is based upon the iterative use of two subproblems. Quadratic programming and Lagrange multipliers are used to solve the respective first and second subproblems and to identify the candidate contact surface and contact stresses. This approach guarantees the imposition of the active kinematic contact constraints, avoids the use of special contact elements and the interference of the user in dictating the accuracy of the solution. A modified Newmark formulation is developed to integrate the resulting time‐dependent variational inequality. This newly devised implicit time integration scheme is unconditionally stable, second‐order accurate, avoids numerical oscillations present in the traditional Newmark method, and does not cause numerical dissipation. To demonstrate the versatility and accuracy of the newly proposed algorithm, several examples are examined and compared with existing solutions where the penalty method has been employed. Copyright © 2001 John Wiley & Sons, Ltd.     

结构瞬态响应数值积分法的稳定性

本文借模态分析方法,导出了多自由度系统数值积分中的判稳方法,并以常用的尤拉后差法、纽马克法、呼伯特法为例,阐述了这一通用的数值稳定性分析方法,并得到了一些对选择恰当的差分格式有用的结论.     

随机振动响应计算的精细积分时域平均法

发展了计算结构随机振动响应的精细积分时域平均法,详细讨论了各种激励作用下系统动力响应的精细积分时域平均,并给出算例。与传统的其它离散积分方法如随机中心差分方法、随机纽马克差分方法等比较,本方法具有非常高的精度和计算效率,特别是能给出速度响应方差计算的良好结果。精细积分时域平均法在各种随机振动系统中有广泛的应用前景     

DRBEM solution of exterior nonlinear wave problem using FDM and LSM time integrations

The nonlinear wave equation is solved numerically in an exterior region. For the discretization of the space derivatives dual reciprocity boundary element method (DRBEM) is applied using the fundamental solution of Laplace equation. The time derivative and the nonlinearity are treated as the nonhomogenity. The boundary integrals coming from the far boundary are eliminated using rational and exponential interpolation functions which have decay properties far away from the region of interest. The resulting system of ordinary differential equations in time are solved using finite difference method (FDM) with a relaxation parameter and least squares method (LSM). The proposed methods are examined with numerical test problems in which the behaviours of solutions are known. Although it gives almost the same accuracy with the DRBEM+FDM procedure, DRBEM+LSM solution procedure is preferred, since it is a direct method without the need of a parameter.     

A Computational Analysis to Burgers Huxley Equation

The efficiency of solving computationally partial differential equations can be profoundly highlighted by the creation of precise, higher-order compact numerical scheme that results in truly outstanding accuracy at a given cost. The objective of this article is to develop a highly accurate novel algorithm for two dimensional non-linear Burgers Huxley (BH) equations. The proposed compact numerical scheme is found to be free of superiors approximate oscillations across discontinuities, and in a smooth flow region, it efficiently obtained a high-order accuracy. In particular, two classes of higher-order compact finite difference schemes are taken into account and compared based on their computational economy. The stability and accuracy show that the schemes are unconditionally stable and accurate up to a two-order in time and to six-order in space. Moreover, algorithms and data tables illustrate the scheme efficiency and decisiveness for solving such non-linear coupled system. Efficiency is scaled in terms of L₂ and L_∞ norms, which validate the approximated results with the corresponding analytical solution. The investigation of the stability requirements of the implicit method applied in the algorithm was carried out. Reasonable agreement was constructed under indistinguishable computational conditions. The proposed methods can be implemented for real-world problems, originating in engineering and science.     

The time-marching method of fundamental solutions for wave equations

In this paper, a meshless numerical algorithm is developed for the solution of multi-dimensional wave equations with complicated domains. The proposed numerical method, which is truly meshless and quadrature-free, is based on the Houbolt finite difference (FD) scheme, the method of the particular solutions (MPS) and the method of fundamental solutions (MFS). The wave equation is transformed into a Poisson-type equation with a time-dependent loading after the time domain is discretized by the Houbolt FD scheme. The Houbolt method is used to avoid the difficult problem of dealing with time evolution and the initial conditions to form the linear algebraic system. The MPS and MFS are then coupled to analyze the governing Poisson equation at each time step. In this paper we consider six numerical examples, namely, the problem of two-dimensional membrane vibrations, the wave propagation in a two-dimensional irregular domain, the wave propagation in an L-shaped geometry and wave vibration problems in the three-dimensional irregular domain, etc. Numerical validations of the robustness and the accuracy of the proposed method have proven that the meshless numerical model is a highly accurate and efficient tool for solving multi-dimensional wave equations with irregular geometries and even with non-smooth boundaries.