单点子域积分与差分

通过稳定性分析、显式与隐式积分，表明了单点子域积分相对于差分法的优越性．     

精细积分的非线性动力学积分方程及其解法

给出了非线性动力学积分方程的表达式,针对该方程提出了一个显式预测－校正的单步四阶精度的积分算法,适用于多自由度、强非线性,非保守系统。算例表明该方法精度高、计算量较少。     

结构非线性动力分析混合时间积分并行算法

综合隐式和显式时间积分技术,对结构非线性动力反应分析提出一种并行混合时间积分算法.该算法采用区域分解技术.将并发性引入到算法中,即利用显式时间积分技术进行界面节点积分而利用隐式算法求解局部子区域.为实现并行混合时间积分算法,设计了灵活的并行数据信息流.编写了该算法的程序,在工作站机群实现了数值算例,验证了算法的精度和性能.计算结果表明该算法具有良好的并行性能,优于隐式算法.     

一类奇异积分计算方法及其在断裂力学中的应用

一类奇异积分计算方法及其在断裂力学中的应用任传波，云大真（大连理工大学，大连１１６０２３）在力学及其它的工程计算中，常常遇到奇异积分，不同奇异程度的积分可以采用不同的方法来解决。本文提出的方法可以解决如下一类的奇异积分问题其中１求解方法对于式（１）的...     

边界元法中计算几乎奇异积分的一种无奇异算法

边界元法中存在几乎奇异积分的计算困难。引起边界单元上几乎奇异积分的因素是源点到其邻近单元的最小距离δ。本文拓展文[1]的思想,进一步采用分部积分将δ移出奇异积分式中积分核之外,转换后积分核是δ的正则函数。所以几乎强奇异和超奇异积分被化为无奇异的规则积分与解析积分的和,可由通常的Gauss数值积分解。文中应用此正则化技术求解了弹性力学平面问题的近边界点位移和应力。     

三维间断位移法及强奇异和超奇异积分的处理方法

从积分方程Somigliana等式出发,导出三维状态下单位位错集度的基本解.在此基础上,建立了边界积分方程,并给出了其离散形式.对强奇异和超奇异积分,采用了Hadamard定义的有限部分积分来处理.最后,给出了计算裂纹应力强度因子的算例,并与解析解进行了比较,证实了该方法的有效性.     

积分方程特征值的一个下界估计式及其应用

导出了具有正定对称核的积分方程最小特征值的一个下界估计式，为此为基础，获得了相应特征值的一种算法。     

反平面扩展裂纹的J积分与COD

给出一个以任意速率扩展的反平面裂纹与路径无关的Ｊ积分，证明Ｊ积分扩展裂纹尖端的张开位移（动态ＣＯＤ）之间有的简单的关系，Ｊ积分与能量释放率，动应力强度因子之间也有简单关系，利用这些关系，给出了动态ＣＯＤ与动应力强度因子之间的关系式。     

瑞典圆弧法积分模型的边坡稳定性解析计算

以积分模型代替条分模型可以提高边坡稳定性系数的计算精度,参考坡脚圆和坡底圆的两类滑动面形式,通过水平和竖直积分得到瑞典圆弧法的边坡稳定性解析式,并与不同方法比较。结果表明:忽略条间力的水平积分模型的解要小于竖直积分解,得到更为保守的稳定性系数,是稳定性系数的下限解;通过数学归纳法得到了边坡水平积分模型的广义形式,匀质土层、异形边坡,成层土、含有稳定水位的边坡都可适用于水平积分方法;竖直积分模型逼近稳定性系数的上限解,应用时要注意稳定性系数的放大,避免不安全的评价。竖直积分和水平积分两种方法可分别界定稳定性系数的上下限,为合理选用边坡稳定性系数提供思路。     

非线性动力学积分方程分块积分解法

对于非线性动力学方程组分块地应用精细积分算法,使其化成积分方程表达式,求解的表达式中具有相对低阶的转换矩阵,从而使精细积分更适用于多自由度、强非线性、变系数、非保守系统,针对积分方程提出了一个显示预测-校正的单步四阶精度自起步的精细积分算法。算例表明本方法是有效的。     

RETHINKING TO FINITE DIFFERENCE TIME-STEP INTEGRATIONS

ＲＥＴＨＩＮＫＩＮＧＴＯＦＩＮＩＴＥＤＩＦＦＥＲＥＮＣＥＴＩＭＥ－ＳＴＥＰＩＮＴＥＧＲＡＴＩＯＮＳＺｈｏｎｇＷａｎｘｉｅ（钟万勰）（ＲｅｓｅａｒｔｃｈＩｎｓｔｉｔｕｔｅｏｆＥｎｇｉｎｅｅｒｉｎｇＭｅｃｈａｎｉｃｓ，ＤａｌｉａｎＵｎｉｖｅｒｓｉｔｙｏｆ...     

结构动力分析中时间积分方法进展

叙述了结构动力分析中时间积分方法的最新发展情况,对这一领域的基本原理和思想进行了总结,重点介绍一些新型计算方法的基本性质,为时间积分方法的进一步研究奠定基础。     

Exponential time integration using Krylov subspaces

The application of exponential integrators based on Krylov techniques to large‐scale simulations of complex fluid flows with multiple time‐scales demonstrates the efficiency of these schemes in reducing the associated time‐step restrictions due to numerical stiffness. Savings of approximately 50% can be achieved for simulations of the three‐dimensional compressible Navier–Stokes equations while still maintaining a truncation error typical of explicit time‐stepping schemes. Exponential time integration techniques of this type are particularly advantageous for fluid flows with a wide range of temporal scales such as low‐Mach number, reactive or acoustically dominated flows. Copyright © 2008 John Wiley & Sons, Ltd.     

Taylor-Galerkin-based spectral element methods for convection-diffusion problems

Several explicit Taylor-Galerkin-based time integration schemes are proposed for the solution of both linear and non-linear convection problems with divergence-free velocity. These schemes are based on second-order Taylor series of the time derivative. The spatial discretization is performed by a high-order Galerkin spectral element method. For convection-diffusion problems an operator-splitting technique is given that decouples the treatment of the convective and diffusive terms. Both problems are then solved using a suitable time scheme. The Taylor-Galerkin methods and the operator-splitting scheme are tested numerically for both convection and convection-diffusion problems.     

Provably stable and time‐accurate extensions of Runge–Kutta schemes for CFD computations on moving grids

Extending fixed‐grid time integration schemes for unsteady CFD applications to moving grids, while formally preserving their numerical stability and time accuracy properties, is a nontrivial task. A general computational framework for constructing stability‐preserving ALE extensions of Eulerian multistep time integration schemes can be found in the literature. A complementary framework for designing accuracy‐preserving ALE extensions of such schemes is also available. However, the application of neither of these two computational frameworks to a multistage method such as a Runge–Kutta (RK) scheme is straightforward. Yet, the RK methods are an important family of explicit and implicit schemes for the approximation of solutions of ordinary differential equations in general and a popular one in CFD applications. This paper presents a methodology for filling this gap. It also applies it to the design of ALE extensions of fixed‐grid explicit and implicit second‐order time‐accurate RK (RK2) methods. To this end, it presents the discrete geometric conservation law associated with ALE RK2 schemes and a method for enforcing it. It also proves, in the context of the nonlinear scalar conservation law, that satisfying this discrete geometric conservation law is a necessary and sufficient condition for a proposed ALE extension of an RK2 scheme to preserve on moving grids the nonlinear stability properties of its fixed‐grid counterpart. All theoretical findings reported in this paper are illustrated with the ALE solution of inviscid and viscous unsteady, nonlinear flow problems associated with vibrations of the AGARD Wing 445.6. Copyright © 2011 John Wiley & Sons, Ltd.     

Exact integration of constitutive equations of kinematic hardening materials and its extended applications

In this paper, an exact formula for the integration of the constitutive equations of kinematic hardening material is presented. Its algorithms are simple and clear. For isotropic hardening or mixed hardening material, the formula is still an exact solution for the case of radial loading, and it is an approximate solution with reasonable accuracy for the case of non-radial loading. The computation results show that the procedure proposed in this paper improves both accuracy and efficiency of numerical integration schemes adopted widely in elastic-plastic finite element analysis.     

Nonlinear Dynamics of Shells: Theory,Finite Element Formulation,and Integration Schemes

The paper is concerned with a dynamical formulation of a recently established shell theory capable to catch finite deformations and falls within the class of geometrically exact shell theories. A basic aspect is the design of time integration schemes which preserve specific features of the continuous system such as conservation of momentum, angular momentum, and energy when the applied forces allow to. The integration method differs from the one recently proposed by Simo and Tarnow in being applicable without modifications to shell formulations with linear as well as nonlinear configuration spaces and in being independent of the nonlinearities involved in the strain-displacement relations. A finite element formulation is presented and various examples of nonlinear shell dynamics including large overall and chaotic motions are considered.     

On Fractional Adaptive Control

Introducing fractional operators in the adaptive control loop, and especially in Model Reference Adaptive Control (MRAC), has proven to be a good mean for improving the plant dynamics with respect to response time and disturbance rejection. The idea of introducing fractional operators in adaptation algorithms is very recent and needs to be more established, that is why many research teams are working on the subject. Previously, some authors have introduced a fractional model reference in the adaptation scheme, and then fractional integration has been used to deal directly with the control rule. Our original contribution in this paper is the use of a fractional derivative feedback of the plant output, showing that this scheme is equivalent to the fractional integration, one with a certain benefit action on the system dynamical behaviour and a good robustness effect. Numerical simulations are presented to show the effectiveness of the proposed fractional adaptive schemes.     

On unconditionally positive implicit time integration for the DG scheme applied to shallow water flows

We present a new unconditionally positivity‐preserving (PP) implicit time integration method for the DG scheme applied to shallow water flows. This novel time discretization enhances the currently used PP DG schemes, because in the majority of previous work, explicit time stepping is implemented to deal with wetting and drying. However, for explicit time integration, linear stability requires very small time steps. Especially for locally refined grids, the stiff system resulting from space discretization makes implicit or partially implicit time stepping absolutely necessary. As implicit schemes require a lot of computational time solving large systems of nonlinear equations, a much larger time step is necessary to beat explicit time stepping in terms of CPU time. Unfortunately, the current PP implicit schemes are subject to time step restrictions due to a so‐called strong stability preserving constraint. In this work, we hence give a novel approach to positivity preservation including its theoretical background. The new technique is based on the so‐called Patankar trick and guarantees non‐negativity of the water height for any time step size while still preserving conservativity. In the DG context, we prove consistency of the discretization as well as a truncation error of the third order away from the wet–dry transition. Because of the proposed modification, the implicit scheme can take full advantage of larger time steps and is able to beat explicit time stepping in terms of CPU time. The performance and accuracy of this new method are studied for several classical test cases. Copyright © 2014 John Wiley & Sons, Ltd.