基于S-R和分解定理的几何非线性问题的数值计算分析

为了探究几何非线性问题的数值求解方法,采用理论推导、MATLAB编程计算、有限元模拟相结合的方法,基于S-R和分解定理及更新拖带坐标描述法,运用插值型无单元Galerkin方法对几何非线性问题的增量变分方程进行了推导,并通过四点Gauss积分法和不动点迭代法对其进行求解.最后以平面悬臂梁的大变形问题为例进行求解计算,发现与ANSYS的计算结果拟合相似度很高,说明了所采用的几何非线性力学理论及数值计算方法的正确性和合理性,为求解几何非线性问题提供了一种新的依据.     

考虑地基耦合效应含裂纹中厚矩形板的非线性振动分析

基于Reissner板理论和Hamilton变分原理,建立了双参数地基上具有表面横向贯穿裂纹的中厚矩形板的非线性运动控制方程.在周边自由的条件下,提出了一组满足问题全部边界条件和裂纹处连续条件的试函数.且利用Galerkin法和谐波平衡法对方程进行求解,分析了考虑地基耦合效应的中厚矩形裂纹板的非线性振动特性.数值计算中,讨论了不同裂纹位置、裂纹深度、板的结构参数和地基物理参数对弹性地基上具裂纹的四边自由中厚矩形板的非线性幅频响应的影响.     

无单元Galerkin法在热传导中的应用

对热传导问题的微分方程采用无单元Galerkin法进行数值求解.首先,将微分方程用Galerkin加权残量法转化为等效的积分形式.然后,先将时间变量看作参数,对空间变量进行离散化,得到方程的半离散形式,接着,对时间采用向后Euler—Galerkin格式进行离散,得到方程的全离散形式最后,编制MATLAB程序,上机计算.列举了两个热传导算例,通过计算说明EFG法适用于热传导问题,且其计算速度快,精确度高、前后处理也十分方便,是一种具有潜力的温度场数值计算的新方法.     

二维瞬态热传导问题的无单元Galerkin法分析

采用无单元Galerkin(element-free Galerkin,EFG)法求解具有混合边界条件的二维瞬态热传导问题.首先采用二阶向后微分公式离散热传导方程的时间变量,将该问题转化为与时间无关的混合边值问题;然后采用罚函数法处理Dirichlet边界条件,建立了二维瞬态热传导问题的无单元Galerkin法;最后基于移动最小二乘近似的误差结果,详细推导了无单元Galerkin法求解二维瞬态热传导问题的误差估计公式.给出的数值算例表明计算结果与解析解或已有数值解吻合较好,该方法具有较高的计算精度和较好的收敛性.     

求解Bratu型方程的径向基函数逼近法

基于径向基函数可以逼近几乎所有函数的强大逼近功能,借鉴弹塑性静力学的处理方法,提出位移、速度、加速度联合插值的径向基函数表达式,结合MATLAB数值软件进行计算机编程,成功求解了Bratu型强非线性方程,并给出相应的相对误差.通过分析几种典型的算例,并将计算结果与一些现有的数值分析法得到的数值解进行对比,表明了该方法的可行性和精确性,为求解强非线性Bratu型方程提供了一种新思路.     

Hamilton-Jacobi方程的小波Galerkin方法

本文选择Daubechies小波尺度函数空间作为Galerkin方法的测试函数空间,并将其应用于Hamilton-Jacobi方程,得到了求解Hamilton-Jacobi方程的小波Galerkin方法的数值格式．由于小波在时间和频率上的局部性,本方法适用于处理具有奇异解的问题,可以有效地防止数值振荡．数值试验显示,本方法是有效的．     

一种抑制激波计算中数值振荡现象的双重小波收缩方法

在激波数值计算中,容易出现数值振荡的问题,振荡激烈时会掩盖真实解,为此提出了许多高精度复杂计算格式或采用人工粘性抑制数值振荡.从信号处理的角度,提出双重小波收缩方法,它能自适应提取激波数值振荡解中的真实物理解.先用局部微分求积法求解浅水波方程和理想流体Euler运动方程中的激波问题,发现其数值振荡现象严重,然后采用双重小波收缩方法对其处理,获得了无数值振荡解,它能准确捕捉激波的位置并且保持激波结构.相比于复杂的Riemann(黎曼)求解格式,借助小波收缩方法,可以采用相对简单的计算格式如微分求积法求解激波问题.     

求解非线性Black-Scholes模型的自适应小波精细积分法

针对非线性Black-Scholes方程,基于quasi-Shannon小波函数给出了一种求解非线性偏微分方程的自适应多尺度小波精细积分法.该方法首先利用插值小波理论构造了用于逼近连续函数的多尺度小波插值算子,利用该算子可以将非线性Black-Scholes方程自适应离散为非线性常微分方程组;然后将用于求解常微分方程组的精细积分法和小波变换的动态过程相结合,并利用非线性处理技术(如同伦分析技术)可有效求解非线性Black-Scholes方程.数值结果表明了该方法在数值精度和计算效率方面的优越性.     

非线性动力系统中两鞍-结分岔点间非稳定曲线的确定

采用将伪弧长延拓法与Poincaré映射法相结合的方法,确定非自治动力系统中两鞍-结分岔点间非稳定曲线,并对采用一般延拓法时出现的奇异性进行了证明。该方法引入了一正则化方程,避免了在求解过程中出现的奇异问题,并给出了相应的迭代格式。在曲线的延拓过程中,由于存在两个延拓方向,为保证将曲线延拓出来,给出了一种确定切线方向的方法,该方法在分析非线性振动系统中的双稳态现象等问题是很有效的。     

一种直接间断Galerkin方法的误差分析

提出一种求解非线性扩散方程的直接间断Galerkin方法.在空间离散上,采用Stirling插值公式构造了一种含有高阶导数项的数值流;在时间离散上,采用通常的Runge-Kutta方法.在理论上,证明了由这种直接间断Galerkin方法得到的计算格式是非线性稳定的,相应数值解不仅具有最优的能量误差估计,而且具有最优的L~2误差估计.数值实验验证了方法的有效性.     

Secondary Buckling Analysis of Thin Rectangular Plates Based on the Wavelet Galerkin Method北大核心CSCD

Application of the wavelet Galerkin method (WGM) to numerical solution of nonlinear buckling problems was studied with classical elastic thin rectangular plates. First, the discretized scheme of the von Kármán equation were introduced, then a simple calculation approach to the Jacobian and Hessian matrices based on the WGM was proposed, and the wavelet discretized scheme-based eigenvalue equation method, the extended equation method and the pseudo arc-length method for nonlinear buckling analysis were discussed. Second, the secondary post-buckling equilibrium paths of elastic thin rectangular plates and the effects of aspect ratios, boundary conditions and bi-directional compression on the mode jumping behaviors, were discussed in detail. Numerical results show that, the WGM possesses good convergence for solving buckling loads on rectangular plates, and the obtained equilibrium paths are in good agreement with those from the stability experiments, the 2-step perturbation method and the nonlinear finite element method. Given the feasibility of combination with different bifurcation computation methods, the WGM makes an efficient spatial discretization method for complex nonlinear stability problems of typical plates and shells. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Finite dimensional approximation of two-parameter nonlinear problems

Summary. In this paper we study a general theory for the numerical approximation of functional nonlinear two-parameter problems in a neighbourhood of an isola center. The results are also valid for a certain class of perturbed bifurcation points. The abstract theory is applied to the Galerkin approximation of nonlinear variational posed problems. In this case, as a consequence of the error being orthogonal to the approximating space, we prove the superconvergence of the perturbation parameter, whereas for the bifurcation parameter and the solution we obtain the same order as in the linear problem. Numerical results are given for the one-dimensional Brussellator model. Received June 10, 1992 / Revised version received May 16, 1994     

弹性约束浅拱的非线性动力响应与分岔分析

浅拱采用竖向、转动方向弹性约束时,自振频率和模态与理想的铰支/固结边界存在差异,不同约束刚度将改变外激励下的非线性响应及各种分岔产生的参数域．由浅拱基本假定建立无量纲动力学方程, 采用在频率和模态中考虑约束刚度大小的方法,通过Galerkin全离散和多尺度摄动分析导出极坐标、直角坐标形式的平均方程, 其中方程系数与约束刚度一一对应．用数值方法分析了周期激励下竖向弹性约束系统最低两阶模态之间1∶2内共振时的动力行为, 所得结果与有限元的对比以及平均方程系数的收敛性证明了所采用方法是可行的．随着激励幅值、频率的变化存在若干分岔点,分岔发生时的参数分布与约束刚度值有关,在由分岔点连接的不稳定区或共振区附近,存在一系列稳态解、周期解、准周期解和混沌解窗口,且随参数的变化可观测到倍周期分岔．     

A wavelet method for solving a class of nonlinear boundary value problems

We develop an approximation scheme for a function defined on a bounded interval by combining techniques of boundary extension and Coiflet-type wavelet expansion. Such a modified wavelet approximation allows each expansion coefficient being explicitly expressed by a single-point sampling of the function, and allows boundary values and derivatives of the bounded function to be embedded in the modified wavelet basis. By incorporating this approximation scheme into the conventional Galerkin method, the interpolating property makes the solution of boundary value problems with strong nonlinearity to be very effective and accurate. As an example, we have applied the proposed method to the solution of the Bratu-type equations. Results demonstrate a much better accuracy than most methods developed so far. Interestingly, unlike most existing methods, numerical errors of the present solutions are not sensitive to the nonlinear intensity of the equations.     

Some modifications of the quasilinearization method with higher-order convergence for solving nonlinear BVPs

In this paper, modifications of the quasilinearization method with higher-order convergence for solving nonlinear differential equations are constructed. A general technique for systematically obtaining iteration schemes of order m (?>?2) for finding solutions of highly nonlinear differential equations is developed. The proposed iterative schemes have convergence rates of cubic, quartic and quintic orders. These schemes were further applied to bifurcation problems and to obtain critical parameter values for the existence and uniqueness of solutions. The accuracy and validity of the new schemes is tested by finding accurate solutions of the one-dimensional Bratu and Frank-Kamenetzkii equations.     

拟弧长延拓法在静电激励MEMS吸合特性研究中的应用

在静电激励微机电系统MEMS(micro-electro-mechanical systems)吸合特性研究中,基于应变梯度理论的微梁结构的控制方程是非线性高阶微分方程,给方程的求解带来了困难.由于该问题的数学模型本质上是分叉问题,方程的解支上出现奇异点,而运用局部延拓法无法通过奇异点.因此,通过运用广义微分求积法将控制方程降阶离散,结合拟弧长延拓法使迭代顺利通过奇异点,求出了整个解曲线.结果表明,拟弧长延拓法能有效并准确地求解具有分叉现象的高阶微分方程问题,为精确预测静电激励MEMS的吸合电压提供有力帮助.     

Thin plate spline Galerkin scheme for numerically solving nonlinear weakly singular Fredholm integral equations

The present work proposes a numerical method to obtain an approximate solution of non-linear weakly singular Fredholm integral equations. The discrete Galerkin method in addition to thin-plate splines established on scattered points is utilized to estimate the solution of these integral equations. The thin-plate splines can be regarded as a type of free shape parameter radial basis functions which create an efficient and stable technique to approximate a function. The discrete Galerkin method for the approximate solution of integral equations results from the numerical integration of all integrals in the method. We utilize a special accurate quadrature formula via the non-uniform composite Gauss-Legendre integration rule and employ it to compute the singular integrals appeared in the scheme. Since the approach does not need any background meshes, it can be identified as a meshless method. Error analysis is also given for the method. Illustrative examples are shown clearly the reliability and efficiency of the new scheme and confirm the theoretical error estimates.     

An Alternating-Direction Sinc–Galerkin method for elliptic problems

A new Alternating-Direction Sinc–Galerkin (ADSG) method is developed and contrasted with classical Sinc–Galerkin methods. It is derived from an iterative scheme for solving the Lyapunov equation that arises when a symmetric Sinc–Galerkin method is used to approximate the solution of elliptic partial differential equations. We include parameter choices (derived from numerical experiments) that simplify existing alternating-direction algorithms. We compare the new scheme to a standard method employing Gaussian elimination on a system produced using the Kronecker product and Kronecker sum, as well as to a more efficient algorithm employing matrix diagonalization. We note that the ADSG method easily outperforms Gaussian elimination on the Kronecker sum and, while competitive with matrix diagonalization, does not require the computation of eigenvalues and eigenvectors.     

Bifurcation in the neighbourhood of a non-isolated singular point

The Lyapunov-Schmidt method for bifurcation problems has, until recently, been applied only to operator equations whose singular points are isolated in the solution set of the equation. For bifurcation at a multiple eigenvalue involving several parameters, however, singular points are often non-isolated. In this paper, the case of intersecting curves of singular points is considered. Under natural hypotheses on these curves, and assuming suitable transversality conditions on the first order nonlinearity of the operator, it is shown that the solution set of the equation may be completely determined locally in terms of the solutions of associated finite dimensional polynomial equations.