基于微分求积法的边界值方法

探究了边界值方法与微分求积法两者之间的关系。利用经典的微分求积公式,系统地构造了三类不同的边界值方法;当采用均匀网格点时,本文所导出的边界值方法与已有的边界值方法是一致的。研究结果揭示了微分求积法与边界值方法两者之间的内在关系,也建立了线性多步法与单步多级方法之间的联系。     

微分求积法的特性及其改进

微分求积法已在科学和工程计算中得到了广泛应用。然而,有关时域微分求积法的数值稳定性、计算精度即阶数等基本特性,仍缺乏系统性的分析结论。依据微分求积法的基本原理,推导证明了微分求积法的权系数矩阵满足V-变换这一重要特性;利用微分求积法和隐式Runge-Kutta法的等值性,证明了时域微分求积法是A-稳定、s级s阶的数值方法。在此基础上,为进一步提高传统微分求积法的计算精度,利用待定系数法和Padé逼近,推导出了一类新的s级2s阶的微分求积法。数值计算对比结果验证了所提出的新微分求积法比传统的微分求积法具有更高的计算精度。     

基于微分求积法的逐步积分法与常用时间积分法的比较

微分求积法具有数学概念简单、精度高和计算时间少等优点,是一种不断受到重视的数值方法.目前微分求积法在方法本身的研究已经相当充分和成熟,而应用方面的研究则大多集中在边值问题的求解,本文的研究集中在采用微分求积法求解动力学初值问题方面.先介绍了一种新近提出的逐步积分方法,该方法基于一种特殊的节点分布的微分求积法.然后通过理论分析与几种常用的时间积分方法进行了稳定性、精确性和计算量的比较.最后计算了一个双质点系在强迫力下的瞬态响应.比较结果表明新近提出的逐步积分方法具有无条件稳定、高精度和计算量少的优点.     

求解具有弹性接头桩基动力学大变形的微分—代数方法

本文采用弧坐标首先建立了求解具有弹性接头的桩基大变形分析的非线性动力学微分方程,其中, 广义Winkler模型用来模拟土对桩基的抗力.其次,在空间域内应用微分求积单元法来离散非线性微分方程组,并给出了处理弹性接头处连接条件的微分求积单元公式,得到了时间域内的一组微分-代数方程,采用二阶向后差分来代替二阶时间导数离散微分-代数方程组,得到一组离散化的非线性代数方程,应用Newton-Raphson方法求解了该非线性代数方程组.最后给出了数值算例,得到了桩基在顶部处受到组合动载荷作用时的响应,考察了弹性接头的刚度、位置对桩基动力学行为的影响.     

轴向均布载荷下压杆稳定问题的DQ解

叙述了微分求积法(differential quadrature method)的一般方法，研究用微分求积法求解在均布轴向载荷下细长杆的稳定问题．通过Newton-Raphson法求解非线性方程组，以及对问题进行线性假设后求解广义特征值方程，得到了精度很高的后屈曲挠度数值和临界载荷数值．与解析解和其他近似解相比，微分求积法具有较高的精度和简便性．     

用微分求积方法计算二维薄板在超音速流中的非线性颤振

采用了一种微分求积方法将二维薄板在超音速气流作用下的非线性动力学方程离散为常微分方程,并用Runge-Kutta数值方法进行了计算.为验证微分求积方法的结果,与伽辽金方法计算结果进行了比较,取得了一致的结果.微分求积法的计算结果用分叉图、相平面、时域曲线以及功率谱进行了描述,结果表明在特定的参数区间存在混沌运动,而通向混沌的道路是经过一系列周期倍化分叉产生的.     

微分求积单元法在结构工程中的应用

微分求积法（Differential Quadrature Method）是求鳃偏微分方程和积分-微分方程的一种数值方法,该法具有计算简便、精度较高和易于实现等优点。微分求积单元法（Differential Quadrature Element Method）是在微分求积法的基础上结合区域分割和集成规则而形成的一种新的数值计算方法,能通过自适应地选取微分求积网点数目正确模拟构件的刚度和荷载性质,其精度可通过细分单元或增加离散点数目加以提高。微分求积单元法是一种可供选择的、性能优越的数值计算方法。本文将详细论述这一数值方法的基本原理,并通过数值算例说明该方法的应用过程及其优越性,为这一方法在结构工程中的推广应用提供参考。     

集中载荷作用下梯度复合材料梁的小波-微分求积法分析

将梯度复合材料梁作为平面应力问题处理,采用小波和微分求积混合法,对集中荷载作用下结构的响应进行了分析.考虑材料特性参数沿高度方向呈梯度分布,在该方向上采用广义微分求积法进行离散;鉴于广义微分求积法求解集中荷载问题精度不高的缺点,在梁的长度方向上引入对突变信号敏感的小波插值函数.数值计算表明,小波-微分求积混合法不仅保留了广义微分求积法高效的优点,而且能够很好地模拟结构局部化特征.     

非齐次动力学方程的一种精细积分单步方法

针对非齐次动力学方程■,结合精细积分法和微分求积法,利用同阶的显式龙格-库塔法对计算过程中待求的v_(k+i/s)(i=1,2,…,s)进行预估,提出了一种避免状态矩阵求逆的高效精细积分单步方法。该方法采用精细积分法计算e~(Ht),而Duhamel积分项采用s级s阶的时域微分求积法,计算格式统一且易于编程,可灵活实现变阶变步长。仿真结果表明,与其他单步法及预估校正-辛时间子域法进行数值比较,该方法具有高精度、高效率及良好的稳定性,在求解大规模动力系统时间响应问题中具有较大的优势。     

微分求积方法对于桩基大变形分析的应用

本文基于大变形的理论,采用弧坐标首先建立了具有初始位移的桩基的非线性数学模型,一组强非线性的微分-积分方程,其中,地基的抗力采用了Winkeler模型;其次,引入变数变换将微分-积分方程转化为一组非线性微分方程,并用微分求积方法离散了方程组,得到一组离散化的非线性代数方程;最后用Newton-Raphson迭代方法对离散化方程进行了求解,得到了桩基变形前后的构形、弯矩和剪力.计算中选取了两种不同类型的初始位移,并考察了它们对桩基大变形力学行为的影响.     

Solution to transient Navier–Stokes equations by the coupling of differential quadrature time integration scheme with dual reciprocity boundary element method

The two‐dimensional time‐dependent Navier–Stokes equations in terms of the vorticity and the stream function are solved numerically by using the coupling of the dual reciprocity boundary element method (DRBEM) in space with the differential quadrature method (DQM) in time. In DRBEM application, the convective and the time derivative terms in the vorticity transport equation are considered as the nonhomogeneity in the equation and are approximated by radial basis functions. The solution to the Poisson equation, which links stream function and vorticity with an initial vorticity guess, produces velocity components in turn for the solution to vorticity transport equation. The DRBEM formulation of the vorticity transport equation results in an initial value problem represented by a system of first‐order ordinary differential equations in time. When the DQM discretizes this system in time direction, we obtain a system of linear algebraic equations, which gives the solution vector for vorticity at any required time level. The procedure outlined here is also applied to solve the problem of two‐dimensional natural convection in a cavity by utilizing an iteration among the stream function, the vorticity transport and the energy equations as well. The test problems include two‐dimensional flow in a cavity when a force is present, the lid‐driven cavity and the natural convection in a square cavity. The numerical results are visualized in terms of stream function, vorticity and temperature contours for several values of Reynolds (Re) and Rayleigh (Ra) numbers. Copyright © 2008 John Wiley & Sons, Ltd.     

Boundary element solution of unsteady magnetohydrodynamic duct flow with differential quadrature time integration scheme

A numerical scheme which is a combination of the dual reciprocity boundary element method (DRBEM) and the differential quadrature method (DQM), is proposed for the solution of unsteady magnetohydrodynamic (MHD) flow problem in a rectangular duct with insulating walls. The coupled MHD equations in velocity and induced magnetic field are transformed first into the decoupled time‐dependent convection–diffusion‐type equations. These equations are solved by using DRBEM which treats the time and the space derivatives as nonhomogeneity and then by using DQM for the resulting system of initial value problems. The resulting linear system of equations is overdetermined due to the imposition of both boundary and initial conditions. Employing the least square method to this system the solution is obtained directly at any time level without the need of step‐by‐step computation with respect to time. Computations have been carried out for moderate values of Hartmann number (M?50) at transient and the steady‐state levels. As M increases boundary layers are formed for both the velocity and the induced magnetic field and the velocity becomes uniform at the centre of the duct. Also, the higher the value of M is the smaller the value of time for reaching steady‐state solution. Copyright © 2005 John Wiley & Sons, Ltd.     

求解非线性方程组的混合遗传算法

非线性方程组的求解是数值计算领域中最困难的问题。大多数的数值求解算法例如牛顿法的收敛性和性能特征在很大程度上依赖于初始点。但是对于很多非线性方程组,选择好的初始点是一件非常困难的事情。本文结合遗传算法和经典算法的优点,提出了一种用于求解非线性方程组的混合遗传算法。该混合算法充分发挥了遗传算法的群体搜索和全局收敛性,有效地克服了经典算法的初始点敏感问题;同时在遗传算法中引入经典算法(Powell法、拟牛顿迭代法)作局部搜索,克服了遗传算法收敛速度慢和精度差的缺点。选择了几个典型非线性方程组,从收敛可靠性、计算成本和适用性等指标分析对比了不同算法。计算结果表明所设计的混合遗传算法有着可靠的收敛性和较高的收敛速度和精度,是求解非线性方程组的一种成功算法。     

用DQM求解成层地基一维非线性固结问题

由于多层地基的一维非线性固结问题求解的复杂性,其解析解很难求得。本文基于Davis和Raymond一维非线性固结理论,利用DQM(Differential Quadrature Method)导了初始有效应力沿深度变化、任意边界条件、任意荷载作用下成层地基一维非线性固结的统一表达式,求得了孔压、有效应力和平均固结度的解答。通过解的收敛性分析讨论了DQM解的有效性。由于DQM解对于固结间题各种复杂条件具有统一的矩阵表达式,更便于编程计算和工程应用。最后,用本文解答对三层地基一维非线性固结问题进行了讨论。     

High order multiplication perturbation method for singular perturbation problems

This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singularly perturbed two-point boundary value problems are first transformed into the singularly perturbed initial value problems. With the variable coefficient dimensional expanding, the non-homogeneous ordinary dif- ferential equations （ODEs） are transformed into the homogeneous ODEs, which are then solved by the high order multiplication perturbation method. Some linear and nonlinear numerical examples show that the proposed method has high precision.     

On singular perturbation boundary-value problem of coupling type system of convection-diffusion equations

In this paper we consider the singular perturbation boundary-value problem of thefollowing coupling type system of convection-diffusion equationsWe advance two methods:the first one is the initial value solving method,by which theoriginal boundary-value problem is changed into a series of unperturbed initial-valueproblems of the first order ordinary differential equation or system so that an asymptoticexpansion is obtained;the second one is the boundary-value solving method,by which theoriginal problem is changed into a few boundary-value problems having no phenomenon ofboundary-layer so that the exact solution can be obtained and any classical numericalmethods can be used to obtain the numerical solution of consismethods can be used to obtainthe numerical solution of consistant high accuracy with respect to the perturbationparameterε     

Analysis of velocity equation of steady flow of a viscous incompressible fluid in channel with porous walls

Steady flow of a viscous incompressible fluid in a channel, driven by suction or injection of the fluid through the channel walls, is investigated. The velocity equation of this problem is reduced to nonlinear ordinary differential equation with two boundary conditions by appropriate transformation and convert the two‐point boundary‐value problem for the similarity function into an initial‐value problem in which the position of the upper channel. Then obtained differential equation is solved analytically using differential transformation method and compare with He's variational iteration method and numerical solution. These methods can be easily extended to other linear and nonlinear equations and so can be found widely applicable in engineering and sciences. Copyright © 2009 John Wiley & Sons, Ltd.     

基于精细积分技术的非线性动力学方程的同伦摄动法

将精细积分技术（PIM）和同伦摄动方法（HPM）相结合,给出了一种求解非线性动力学方程的新的渐近数值方法。采用精细积分法求解非线性问题时,需要将非线性项对时间参数按Taylor级数展开,在展开项少时,计算精度对时间步长敏感;随着展开项的增加,计算格式会变得越来越复杂。采用同伦摄动法,则具有相对筒单的计算格式,但计算精度较差,应用范围也限于低维非线性微分方程。将这两种方法相结合得到的新的渐近数值方法则同时具备了两者的优点,既使同伦摄动方法的应用范围推广到高维非线性动力学方程的求解,又使精细积分方法在求解非线性问题时具有较简单的计算格式。数值算例表明,该方法具有较高的数值精度和计算效率。     

康托洛维奇-里茨杂交法及其应用

本文将Kantorovich法与Ritz法进行适当组合,吸收了二者的主要优点,提出了康托洛维奇法的一种改进方法.以二维问题为例,Kantorovich法在一个方向(例如y方向)的分布完全预先选定,这含有很大的主观任意性,因而限制了近似解的精度,改进法则在y方向仿Ritz法改进为一个含有若干个自由参数的分布函数,由于增加了近似解的自由度,故可改善解的精度.对Kantorovich法的另一改进是在计算高阶近似解时,通过逐项求解待定函数避免了求解更高阶微分方程或含更多方程的方程组,减少了计算量,降低了计算难度.用改进法求解了固体力学里的矩形截面柱体扭转问题和四边固支矩形板的弯曲问题,通过算例充分说明了此方法的特点和优越性.     

A simple and accurate mixed FE-DQ formulation for free vibration of rectangular and skew Mindlin plates with general boundary conditions

A simple and accurate mixed finite element-differential quadrature formulation is proposed to study the free vibration of rectangular and skew Mindlin plates with general boundary conditions. In this technique, the original plate problem is reduced to two simple bar (or beam) problems. One bar problem is discretized by the finite element method (FEM) while the other by the differential quadrature method (DQM). The mixed method, in general, combines the geometry flexibility of the FEM and high accuracy and efficiency of the DQM and its implementation is more easier and simpler than the case where the FEM or DQM is fully applied to the problem. Moreover, the proposed formulation is free of the shear locking phenomenon that may be encountered in the conventional shear deformable finite elements. A simple scheme is also presented to exactly implement the mixed natural boundary conditions of the plate problem. The versatility, accuracy and efficiency of the proposed method for free vibration analysis of rectangular and skew Mindlin plates are tested against other solution procedures. It is revealed that the proposed method can produce highly accurate solutions for the natural frequencies of rectangular and skew Mindlin plates with general boundary conditions.