An implicit multi-time step integration method for structural dynamics problems

A multi-time step integration algorithm is developed based on the trapezoidal rule time integration method for finite element equations of motion. This algorithm uses nodal groups to partition the mesh into subdomains that are updated with different time steps. A␣stability analysis of the method shows that the scheme retains the unconditionally stable behavior of the trapezoidal rule and conserves the same pseudo energy as the parent algorithm. Several numerical examples are used to verify the stability of the method and to investigate the accuracy of the scheme.     

On time integration of viscoplastic constitutive models suitable for creep

Creep of critical components such as electrical solder connections may occur over long periods of time. Efficient numerical simulations of such problems generally require the use of quasi‐static formulations with conjugate‐gradient techniques for solving the large number of algebraic equations. Implicit in the approach is the need to solve the constitutive equation several times for large time steps and for loading directions that may have no resemblance to the actual solution. Therefore, an unconditionally stable and efficient algorithm for solving the constitutive equation is essential for the overall efficiency of the solution procedure. Unfortunately, constitutive equations suitable for simulating the materials of interest are notoriously difficult to solve numerically and most existing algorithms have a stability limit on the time step which may be several orders of magnitude smaller than the desired time step. Here an algorithm is proposed which is a combination of the use of a trapezoidal rule and an iterative Newton–Raphson method for solving implicitly the non‐linear equations. The key to the success of the proposed approach is to always use an initial guess based on the steady‐state solution to the constitutive equation. A representative viscoplastic constitutive equation is used as a model for illustrating the approach. The algorithm is developed and typical numerical results are provided to substantiate the claim that stability has been achieved. Copyright © 2001 John Wiley & Sons, Ltd.     

A time integration algorithm for elasto-viscoplastic cubic crystals applied to modelling polycrystalline deformation

An algorithm for integrating the constitutive equations for an elasto-viscoplastic cubic crystal is presented which is shown to be easily employed in a polycrystalliné analysis. Anisotropic elastic behaviour is incorporated into the standard constitutive equations for ductile single crystals. The algorithm is shown to be efficient, robust and general. The primary advantage of this algorithm is that is provides an implicit integration of the plastic deformation gradient while including the elastic response. This permits taking large time steps while maintaining accuracy and stability. Several polycrystalline examples are presented to demonstrate the effect of the time step on the solution. Examples also are presented which compare the algorithm described herein to an algorithm which neglects the elastic part of the deformation. In addition, the effect of the anisotropic component of the elasticity is investigated by comparing the results with those obtained assuming isotropic elasticity.     

Analysis and implementation of a new constant acceleration subcycling algorithm

An algorithm for explicit integration of structural dynamics problems with multiple time steps is proposed that averages accelerations to obtain subcycle states at a nodal interface between regions integrated with different time steps. With integer time step ratios, the resulting subcycle updates at the interface sum to give the same effect as a central difference update over a major cycle. The algorithm is shown to have good accuracy, and stability properties in linear elastic analysis similar to those of constant velocity subcycling algorithms. The implementation of a generalised form of the algorithm with non-integer time step ratios is presented. © 1997 John Wiley & Sons, Ltd.     

Time‐dependent fractional advection–diffusion equations by an implicit MLS meshless method

Recently, many new applications in engineering and science are governed by a series of fractional partial differential equations (FPDEs). Unlike the normal partial differential equations (PDEs), the differential order in an FPDE is with a fractional order, which will lead to new challenges for numerical simulation, because most existing numerical simulation techniques are developed for the PDE with an integer differential order. The current dominant numerical method for FPDEs is finite difference method (FDM), which is usually difficult to handle a complex problem domain, and also difficult to use irregular nodal distribution. This paper aims to develop an implicit meshless approach based on the moving least squares (MLS) approximation for numerical simulation of fractional advection–diffusion equations (FADE), which is a typical FPDE The discrete system of equations is obtained by using the MLS meshless shape functions and the meshless strong‐forms. The stability and convergence related to the time discretization of this approach are then discussed and theoretically proven. Several numerical examples with different problem domains and different nodal distributions are used to validate and investigate the accuracy and efficiency of the newly developed meshless formulation. It is concluded that the present meshless formulation is very effective for the modeling and simulation of the FADE. Copyright © 2011 John Wiley & Sons, Ltd.     

An implicit RBF meshless approach for time fractional diffusion equations

This paper aims to develop an implicit meshless approach based on the radial basis function (RBF) for numerical simulation of time fractional diffusion equations. The meshless RBF interpolation is firstly briefed. The discrete equations for two-dimensional time fractional diffusion equation (FDE) are obtained by using the meshless RBF shape functions and the strong-forms of the time FDE. The stability and convergence of this meshless approach are discussed and theoretically proven. Numerical examples with different problem domains and different nodal distributions are studied to validate and investigate accuracy and efficiency of the newly developed meshless approach. It has proven that the present meshless formulation is very effective for modeling and simulation of fractional differential equations.     

Multiple-parameter reduced basis technique for bifurcation and post-buckling analyses of composite plates

A multiple-parameter reduced basis technique and a problem-adaptive computational algorithm are presented for the bifurcation and post-buckling analyses of composite plates subjected to combined loadings. The computational algorithm can be conveniently divided into three distinct stages. The first stage is that of determining the stability boundary. The plate is discretized by using displacement finite element models and the analysis region is reduced by exploiting the special symmetries exhibited by the response of the plate. The vector of unknown nodal displacements is expressed as a linear combination of a small number of path derivatives (derivatives of the nodal displacements with respect to path parameters), and a Rayleigh-Ritz technique is used to approximate the finite element equations by a small system of algebraic equations. The reduced equations are used to determine the stability boundary of the plate. In the second stage, a nonllnear solution in the vicinity of the stability boundary is obtained by using a bifurcation buckling mode as a predictor, and a set of reduced equations is generated. In the third stage, the reduced equations are used to trace post-buckling paths corresponding to various combinations of the load parameters. The potential of the proposed approach is discussed and its effectiveness is demonstrated by means of a numerical example of laminated composite plate subjected to combined compressive and shear loadings.     

Reactive power transfer allocation method with the application of artificial neural network

A novel method to identify the reactive power transfer between generators and load using modified nodal equations is proposed. On the basis of the solved load flow results, the method partitions the Y-bus matrix to decompose the current of the load buses as a function of the generators' current and voltage. Then it uses the load voltages from the load flow results and decomposed load currents to determine reactive power contribution from each generator to loads. The validation of the proposed methodology is demonstrated by using a simple 3-bus system and the 25-bus equivalent system of south Malaysia. Next part here focuses on creating an appropriate artificial neural network (ANN) to solve the same problem in a simpler and faster manner. The basic idea is to use supervised learning paradigm to train the ANN. Most commonly used feedforward architecture has been chosen for the proposed ANN reactive power transfer allocation technique. Almost all system variables obtained from load flow solutions are utilised as an input to the neural network. Moreover, tan-sigmoid activation functions are incorporated in the hidden layer to realise the nonlinear nature of the reactive power transfer allocation. The targets of the ANN corresponding to the previously developed reactive power transfer allocation method. The 25-bus equivalent system of south Malaysia is utilised as a test system to illustrate the effectiveness of the ANN output compared with that of the modified nodal equations method. The ANN output provides promising results in terms of accuracy and computation time.     

A FETI‐based domain decomposition technique for time‐dependent first‐order systems based on a DAE approach

We present a novel partitioned coupling algorithm to solve first‐order time‐dependent non‐linear problems (e.g. transient heat conduction). The spatial domain is partitioned into a set of totally disconnected subdomains. The continuity conditions at the interface are modeled using a dual Schur formulation where the Lagrange multipliers represent the interface fluxes (or the reaction forces) that are required to maintain the continuity conditions. The interface equations along with the subdomain equations lead to a system of differential algebraic equations (DAEs). For the resulting equations a numerical algorithm is developed, which includes choosing appropriate constraint stabilization techniques. The algorithm first solves for the interface Lagrange multipliers, which are subsequently used to advance the solution in the subdomains. The proposed coupling algorithm enables arbitrary numeric schemes to be coupled with different time steps (i.e. it allows subcycling) in each subdomain. This implies that existing software and numerical techniques can be used to solve each subdomain separately. The coupling algorithm can also be applied to multiple subdomains and is suitable for parallel computers. We present examples showing the feasibility of the proposed coupling algorithm. Copyright © 2008 John Wiley & Sons, Ltd.     

Multi-time-step integration using nodal partitioning

An algorithm is presented which integrates different groups of nodes of a finite element mesh with different time steps and different integrators. Since the nodal groups are updated independently no unsymmetric systems need be solved. Stability is demonstrated by showing that an energy norm of the solution decreases after every update if the time step is less than a given critical value. The element eigenvalue inequality theorem is used to give the critical time step in terms of element eigenvalues.     

Buoyancy potential conserving technique for dynamic analysis of offshore pipes and risers

The paper describes the midpoint algorithm for the solution of dynamic equilibrium equations in the marine environment for the non‐linear dynamic analysis of flexible pipes and riser systems in 3D configuration. It has been shown that proposed method has an advantage over traditional time marching method by means of long‐term stability and larger time steps. The mid‐point algorithm used for the time integration is considered ‘approximate buoyancy‐energy conserving algorithm’. Emphasis is placed on elements that use the Reissner–Simo beam theory or a related co‐rotational approach. Copyright © 2005 John Wiley & Sons, Ltd.     

饱和两相介质近场波动问题的一种时域全显式数值计算方法

基于 u -p形式的饱和两相介质弹性波动方程,开展了饱和两相介质近场波动问题时域显式数值计算方法的研究。通过对波动方程中的质量矩阵和孔隙流体压缩矩阵进行对角化处理,消除了方程中的动力耦联,实现了波动方程的解耦。分别应用中心差分法和Newmark常平均加速度法求解固相位移和速度,基于向后差分法求解孔隙流体压力,推导得到了饱和两相介质动力响应的时域显式逐步积分的计算列式,建立了饱和两相介质近场波动问题的一种新的时域全显式数值计算方法。进行了该文方法中矩阵对角化合理性的验证。将该方法的数值解与相应的解析解进行对比,二者符合良好,验证了该方法的正确性。将该文建立的时域数值计算方法与透射人工边界方法相结合,应用于饱和两相介质的近场波动问题,进行了饱和土场地地震响应的计算研究,计算结果符合弹性波动理论的基本规律,表明该方法对于饱和两相介质近场波动问题时域计算求解的适用性。基于该方法中时域递推计算格式的传递矩阵,进行了该方法稳定性特性的研究。该文建立的数值计算方法具有时域全显式算法的基本特征。方法中对动力响应的全部分量均采用递推和迭代的模式进行求解,避免了求解耦联的动力方程组。该方法具有较高的计算效率,是进行饱和两相介质近场波动问题时域计算求解的一种有效的算法。     

A simple algorithm for localized construction of non‐matching structural interfaces

A simple and effective algorithm for the modular construction of non‐matched interfaces is presented for the partitioned solution of large‐scale structural problems. The formulation is based on a recently developed four‐field variational principle, which introduces a connection frame between the interfaced partitions. A key result of the present study is a frame nodal placement criterion that uniquely determines the frame discretization into piecewise linear elements so that the interface patch test condition is satisfied a priori. The method is demonstrated with several 2D and 3D example problems. Copyright © 2002 John Wiley & Sons, Ltd.     

Singularity confinement for a class of m-th order difference equations of combinatorics

In a recent publication, it was shown that a large class of integrals over the unitary group U(n) satisfy nonlinear, non-autonomous difference equations over n, involving a finite number of steps; special cases are generating functions appearing in questions of the longest increasing subsequences in random permutations and words. The main result of the paper states that these difference equations have the discrete Painlevé property; roughly speaking, this means that after a finite number of steps the solution to these difference equations may develop a pole (Laurent solution), depending on the maximal number of free parameters, and immediately after be finite again ("singularity confinement"). The technique used in the proof is based on an intimate relationship between the difference equations (discrete time) and the Toeplitz lattice (continuous time differential equations); the point is that the Painlevé property for the discrete relations is inherited from the Painlevé property of the (continuous) Toeplitz lattice.     

Timoshenko梁能量守恒逐步积分算法

该文提出了Timoshenko梁非线性动力分析的能量守恒逐步积分算法。采用共旋技术考虑结构的几何非线性,空间离散采用相关插值形式,避免了剪切锁定现象。在时间离散时利用多参数修正方法对等效的节点动力平衡方程进行修正,实现了离散系统在保守荷载作用下的能量守恒。算法具备二阶局部精度,与已有的平均加速度方法和隐式中点方法相比,具有更好的数值稳定性。在二维情形下与Simo方法对比,指出了Simo方法在受保守外弯矩作用时系统能量不守恒。最后,通过三个数值模拟算例验证了算法的性能和能量守恒特性。     

A SPACE-TIME COUPLED p-VERSION LEAST SQUARES FINITE ELEMENT FORMULATION FOR UNSTEADY TWO-DIMENSIONAL NAVIER–STOKES EQUATIONS

This paper presents a p-version least squares finite element formulation for two-dimensional unsteady fluid flow described by Navier–Stokes equations where the effects of space and time are coupled. The dimensionless form of the Navier–Stokes equations are first cast into a set of first-order differential equations by introducing auxiliary variables. This permits the use of C⁰ element approximation. The element properties are derived by utilizing the p-version approximation functions in both space and time and then minimizing the error functional given by the space–time integral of the sum of squares of the errors resulting from the set of first-order differential equations. This results in a true space–time coupled least squares minimization procedure. The application of least squares minimization to the set of coupled first-order partial differential equations results in finding a solution vector {δ} which makes gradient of error functional with respect to {δ} a null vector. This is accomplished by using Newton's method with a line search. A time marching procedure is developed in which the solution for the current time step provides the initial conditions for the next time step. Equilibrium iterations are carried out for each time step until the error functional and each component of the gradient of the error functional with respect to nodal degrees of freedom are below a certain prespecified tolerance. The space–time coupled p-version approximation functions provide the ability to control truncation error which, in turn, permits very large time steps. What literally requires hundreds of time steps in uncoupled conventional time marching procedures can be accomplished in a single time step using the present space–time coupled approach. The generality, success and superiority of the present formulation procedure is demonstrated by presenting specific numerical examples for transient couette flow and transient lid driven cavity. The results are compared with the analytical solutions and those reported in the literature. The formulation presented here is ideally suited for space–time adaptive procedures. The element error functional values provide a mechanism for adaptive h, p or hp refinements.     

A novel accurate and computationally efficient integration approach to viscoplastic constitutive model

A novel, accurate, and computationally efficient integration approach is developed to integrate small strain viscoplastic constitutive equations involving nonlinear coupled first-order ordinary differential equations. The developed integration scheme is achieved by a combination of the implicit backward Euler difference approximation and the implicit asymptotic integration. For the uniaxial loading case, the developed integration scheme produces accurate results irrespective of time steps. For the multiaxial loading case, the accuracy and computational efficiency of the developed integration scheme are better than those of either the implicit backward Euler difference approximation or the implicit asymptotic integration. The simplicity of the developed integration scheme is equivalent to that of the implicit backward Euler difference approximation since it also reduces the solution of integrated constitutive equations to the solution of a single nonlinear equation. The algorithm tangent constitutive matrix derived for the developed integration scheme is consistent with the integration algorithm and preserves the quadratic convergence of the Newton–Raphson method for global iterations.     

A finite element/nodal volume technique for flow simulation of injection pultrusion

In this paper, a finite element/nodal volume technique is developed to simulate the resin flow through the reinforcement during injection pultrusion processes. The governing equations for the pressure distribution and the conservation of resin mass during injection pultrusion are first derived. The solution algorithm and its numerical implementation are described. In particular, an algorithm is developed to advance the flow front by taking into account both the resin flow relative to the reinforcement and the movement of the pultruded part as a whole. The numerical technique developed is validated against a one-dimensional analytical solution derived. A number of numerical tests are then conducted to investigate the numerical performance and capability of the technique.     

结构动力微分方程的一种高精度摄动解

将结构的位移及速度响应作为状态变量,把结构动力方程转化为状态方程,采用摄动方法求解状态方程,推出一种级数形式的摄动解,同时给出了该文算法的迭代格式和计算步骤。该算法无需对转换矩阵求逆,也无须作指数矩阵运算,仅做矩阵向量相乘及向量求和运算,计算稳定而且效率高,收敛速度快,解的级项数及精度可由允许误差参数直接控制,很容易达到任意精度要求,该方法兼具线性加速法的高效率和精细积分法的高精度,可应用于结构大型稀疏线性动力方程组的求解。最后通过典型算例进一步验证了该文算法的精度和效率。     

结构动力方程一种新的级数形式的解析解

将结构的位移及速度响应作为状态变量,采用Lyapunov(李雅普诺夫)人工小参数法求解状态方程,导出状态方程的一个新的级数形式的解析解,该解析解还可以推广到非线性动力方程的计算。将秦九韶算法引入级数解的计算,提高了计算的效率和稳定性,同时给出了算法的计算格式和步骤。该算法无需对转换矩阵H求逆,仅使用矩阵向量相乘,计算稳定,精度仅由收敛项数控制,很容易达到任意精度要求,而且适合并行计算及压缩存储。最后通过算例进一步证实了该算法的精度和效率。