求解非线性伪抛物方程的重心Lagrange插值配点法

提出了重心Lagrange插值配点法求解一类非线性伪抛物方程。首先,介绍了重心Lagrange插值并给出了微分矩阵表达式。其次,构造了求解非线性伪抛物方程的直接线性化迭代格式、部分线性化迭代格式、Newton线性化迭代格式。再次,未知函数和初边值条件利用重心Lagrange插值函数来近似,利用配点法得到离散方程,获得了方程的矩阵表达式。最后,数值算例表明,重心Lagrange插值配点法具有高精度和高效率的优点。     

重心插值配点法求解Allen-Cahn方程

利用重心插值配点法(重心Lagrange插值配点法和重心有理插值配点法)构造包含时间、空间变量的近似函数,给定Chebyshev点族;将重心插值配点法代入Allen-Cahn方程及定解条件,得到离散方程组,并采用Newton迭代格式求解方程组.数值算例表明:文中的配点法具有较高精度;利用2种配点法的能量函数满足能量递减规律.     

时间分数阶Black-Scholes方程的重心Lagrange插值配点法

针对欧式期权定价的时间分数阶Black-Scholes模型,设计一种重心Lagrange插值配点法格式.首先,采用Laplace变换近似Caputo型分数阶导数,将分数阶方程转化为整数阶方程;然后,在时-空方向上均采用重心Lagrange插值配点法进行离散,构造重心Lagrange插值配点法格式.结果表明：时间分数阶Black-Scholes方程的重心Lagrange插值配点法具有高精度和有效性.     

重心插值配点法求解Volterra积分方程

文章提出了求解Volterra积分方程的一种高精度数值方法：重心插值配点法（包括重心Lagrange插值配点法和重心有理插值配点法）。该方法分为两步：首先对Volterra积分方程采用两种重心插值配点法进行离散，构造出Volterra积分方程的数值求解格式；然后，依次选取第二类Chebyshev节点和等距节点进行数值计算。文章主要研究积分项中含有未知函数的一阶导函数的Volterra积分方程的离散格式构造及数值实现。数值实验结果表明：在使用第二类Chebyshev节点时，用重心Lagrange插值配点法较好；在使用等距节点时，使用重心有理插值配点法较好。     

极坐标系下弹性问题的重心插值配点法

针对岩土工程中的孔洞及曲梁问题,提出一种在极坐标系下求解二维弹性问题的重心插值配点法.该方法分别在r和θ方向分别布置m和n个节点,生成求解区域上的节点.以一维重心Lagrange插值的张量积插值形式近似二维弹性问题的位移函数,代入位移表达的平衡方程和边界条件,平衡方程和边界条件分别在所有的计算节点和边界节点上精确成立,得到极坐标下弹性力学平衡方程和边界条件的离散代数表达式.利用一维重心Lagrange插值微分矩阵,将离散的平衡方程和边界条件表达为矩阵形式.利用置换法施加边界条件,求得在计算节点处的位移,进而通过微分矩阵直接求得计算节点处的应力.数值算例表明:极坐标下重心插值配点法具有计算格式简单、程序实施容易和计算精度高的特点.     

重心插值配点法求解分数阶Fredholm积分方程

重心插值配点法是插值法和配点法的结合和推广,它具有稳定性好、高精度和计算效率高等优点.主要运用高精度无网格重心插值配点法求解分数阶Fredholm积分方程.首先推导了基于分数阶Fredholm积分方程重心插值配点法的离散公式,然后通过理论分析得出其解的存在唯一性与误差分析,最后利用数值算例通过对等距节点与第二类Chebyshev节点的对比,验证了所用方法的高精度和可靠性,并得出影响精度的条件.     

时间分数阶Allen-Cahn方程的重心插值配点法

采用Laplace变换近似Caputo型分数阶导数,将分数阶方程转换成整数阶方程;然后,在时-空方向均采用重心插值配点法离散,非线性项采用Newton迭代格式求解,并给出配点格式的相容性误差分析.数值实验表明:该配点法格式具有较高精度,能满足能量递减规律.     

波动问题的高精度重心有理插值配点法

提出一种数值求解波动问题的高精度重心有理插值配点法。对于给定的时间和空间上的计算节点,采用重心有理插值近似未知函数,建立未知函数关于时间和空间变量导数的微分矩阵。将未知函数的重心有理插值近似函数代入波动问题的控制方程,得到波动问题方程和定解条件的离散代数方程组。利用微分矩阵的记号,将离散后的代数方程组写成简洁的矩阵形式。通过置换法施加边界条件和初始条件,求解代数方程组,得到波动问题在计算节点处的位移值。数值算例表明,重心有理插值配点法具有计算公式简单、计算节点适应性好、程序实施方便和计算精度高的优点。     

对流扩散方程最优控制问题的重心插值配点格式

为了讨论对流扩散方程最优控制问题的重心插值配点格式，首先，借助Lagrange乘子法，推导出由状态方程、伴随方程、最优性方程构成的最优性条件.其次，在空间x,y方向均运用重心插值配点格式离散方程组，并给出该配点格式的相容性分析.最后，数值实验验证格式的有效性，与经典有限差分格式比较，重心插值配点格式用较少的节点数就能具有很高的精度.     

重心插值配点法分析矩形薄板弯曲问题

采用重心Lagrange插值近似未知函数建立未知函数各阶导数的微分矩阵.采用微分矩阵近似未知函数的导数,利用配点法将矩形薄板的控制方程和边界条件离散为代数方程组,通过求解代数方程组,求得矩形薄板的各个离散点的挠度,进而利用微分矩阵求得矩形薄板的内力.给出详细的控制方程和边界条件的离散公式.数值算例表明,重心插值配点法具有原理简单,易于程序实现和数值计算精度高的优点.     

重心插值配点法求解Cahn-Hilliard方程

对Cahn-Hilliard方程中的时、空方向均采用重心插值配点格式(重心Lagrange插值配点格式和重心有理插值配点格式)进行离散,非线性项采用一般迭代法,导出离散的线性代数方程组,并给出重心Lagrange插值的逼近误差估计.数值算例表明:两种重心插值配点格式均具有高精度,且满足能量递减规律.     

多点Newton-Raphson迭代的新几何解释

正如"线性化"揭示了Newton迭代的构造思想一样,本文给出的一种几何解释揭示了多点Newton-Raphson迭代的构造思想,由此我们能够给出它的4阶收敛速度的一个简单证明,以及相关的一些重要结果.此外,我们还将多点Newton-Raphson迭代与Olver迭代、Newton迭代进行了综合比较,结论是:多点Newton-Raphson迭代更实用.     

基于重心型插值的数值计算方法

对以重心型插值作为近似函数,数值求解微分方程初值问题和边值问题的数值计算方法作了介绍。给出了重心Lagrange插值和重心有理插值的计算公式和插值节点类型。归纳了求解微分方程初边值问题的重心插值配点法、重心插值Galerkin法和重心插值单元法的计算公式、边界条件/初始条件的离散和施加方法。     

二维Allen-Cahn方程的有限差分法/配点法求解

对二维Allen-Cahn方程中的时间方向采用有限差分法,空间方向采用重心插值配点法,非线性项采用牛顿迭代法,导出离散的线性代数方程组.最后,通过数值算例验证配点法格式的精度及能量递减规律.     

SAV方法求解Allen-Cahn方程的数值比较

该文研究基于标量辅助变量(SAV)格式下求解Allen-Cahn方程的数值比较.首先给出1维Allen-Cahn方程的SAV格式; 然后,对方程的时间方向采用2阶向后差分(BDF2)格式和Crank-Nicolson(CN)格式离散,对方程的空间方向采用重心Lagrange插值配点法和2阶中心差分法离散,用离散正弦变换(DST)、快速傅里叶变换(FFT)求解差分导出的线性代数方程组; 最后,通过数值算例验证重心Lagrange插值配点法是指数收敛,与差分格式比较,配点格式用较少的点就能达到较高的精度且耗时少,并进一步验证几种SAV离散格式都满足能量递减规律.     

处理优化计算力学结果中出现拐点的方法

处理优化计算力学结果中出现拐点的方法顾冲时，吴中如，大坝观测与土工测试，１９９３年第３期针对在利用计算力学结果求解大坝变形方程（如大坝的梁挠曲线方程）时所出现的拐点问题，提出了用插值补点回归法、动态迭代回归法和动态插值迭代回归法，以求得该问题合理解决...     

极坐标下薄板弯曲问题的重心有理插值法

利用重心有理插值配点法(BRICM)研究了极坐标下薄板的弯曲问题,该方法是以重心有理插值近似未知函数强迫微分方程在离散节点处成立,得到微分方程的离散代数方程组,进而采用重心有理插值的微分矩阵将离散代数方程组表达为矩阵的形式。利用置换法施加边界条件,求解微分方程组。数值算例结果表明,该方法在解决极坐标下薄板弯曲问题上公式简单,程序实施方便且计算精度高。     

基于等几何配点法的柔索非线性动力学分析

针对柔索结构非线性动力学问题,提出了基于非均匀有理B样条(NURBS)的等几何配点法.利用哈密顿原理建立了耦合三个平动自由度和一个扭转自由度的柔索非线性运动方程,基于等参元的思想,利用NURBS基函数的高阶连续性和配点离散强形式的控制偏微分方程组,采用Newmark-β时间积分与修正的Newton-Raphson非线性迭代法求解离散动力学方程.通过与传统有限元法计算结果对比,验证了该数值方法的准确性和高效性.进一步分析了等几何配点法典型的收敛特性,表明该方法能有效处理非线性动力学问题.     

利用重心有理插值配点法求解一、二维对流扩散方程

对流扩散是自然界中一种最为常见的物理现象,在气液固中均可发生,该方程已被广泛应用于飞行器设计、热磁辐射、天气预报、化工反应、生物斑点生长等重要领域.为了进一步提高该微分方程的逼近精度,可通过改进基底函数或者调整离散点分布来实现.借助于重心有理插值数值逼近稳定性好,离散矩阵具有稀疏性等优势,求解了一、二维对流扩散方程.将该数值方法与传统的FDM以及Meshfree等方法进行比较,得到的结论是:重心有理插值配点法在求解一、二维对流扩散问题上具有精度高、条件数小、收敛快等优点.从插值节点的分布效果上看, Chebyshev点比等距网格点更稳定,逼近精度略高,且能有效地抑制"龙格"现象的发生.最后,给出了相应的误差估计与收敛性分析,并使用软件画出了热流密度的分布云图,该图有利于分析对流扩散方程的数值解变化趋势问题.     

再生核质点法在非线性冲击动力学中的应用

利用再生核方法来构造计算域函数的插值函数,通过区域内的质点务数来实现数值计算,满足一致性要求;在大变形、大转动、大应变的情况下,引入Jaumann应力和应变速率张量来表征实际应力和应变速率,进而可以表征出材料冲击条件下的本构关系;将本构关系引入冲击过程虚功原理,采取Lagrangian描述的方法,并通过再生核质点方法将能量积分方程离散,推导了非线性冲击过程的再生核质点法计算控制方程,给出了方程求解的Newmark递推公式和Newton-Raphson迭代方法．通过具体实例来说明了再生核质点法在非线性动力学中的应用过程,并验证了计算方法的正确性．