一维非线性对流扩散方程特征有限元的两重网格算法

针对一维非线性对流扩散方程，构造了特征有限元两重网格算法。该算法只需要在粗网格上进行非线性迭代运算，而在所需要求解的细网格上进行一次线性运算即可。对于非线性对流占优扩散方程，不仅可以消除因对流占优项引起的数值振荡现象，更重要的是可以加快收敛速度、提高计算效率。误差估计表明，只要选取粗网格步长与细网格步长的平方根同阶，就可以使两重网格解与有限元解保持同样的计算精度。     

非线性反应扩散问题两网格混合有限元法的数值分析

对于一类非线性反应扩散问题,给出了一种两网格混合有限元解法.首先在粗网格上求解非线性方程组,然后在细网格上采用了牛顿迭代求解.从数值分析的角度对两网格混合有限元算法进行了研究.数值算例结果表明,与混合有限元方法相比,两网格混合有限元方法在不降低解的精度阶数的条件下,提高了计算速度.     

非线性双曲型方程的混合有限元两层网格算法

针对一类非线性双曲型方程, 利用混合有限元法,构造了1种混合有限元两层网格算法, 给出了两网格方法的误差分析. 结果表明, 当两层网格算法所选取的粗网格和细网格步长满足H=б(h^1/2)时,能获得渐近最优的离散逼近解. 并用数值例子验证了该混合有限元两层网格算法的有效性.     

结构非线性有限元分析现状综述

简要介绍非线性有限元概念及基本算法,浅谈结构非线性有限元分析之现状,例举数值算法和网格划分技术对结构非线性有限元分析的精度和效率的影响。     

非线性抛物方程的自适应有限元算法

本文讨论了一维空间中解非线性抛物方程的自适应有限元算法。对此非线性抛物方程,我们选择3层有限元格式,并以等分布原则为基础,对每一时间层选择非一致网格以保证误差具有较高的精度。此算法对大梯度的问题有较强的实用性。应用此算法的计算实例在文末给出。     

一维非线性弦平衡方程的有限元两重网格算法

针对一维非线性弦的平衡方程,构造了有限元两重网格算法,该算法只需要在粗网格上进行非线性迭代,而在所需要求解的细网格上进行一次线性运算即可。与非线性迭代直接求解结果进行对比可知,有限元两重网格算法在保持了计算精度的前提下,所用的时间更短,从而证明了该算法是一种求解非线性问题的高效方法。     

半线性椭圆方程的一个新的双重网格差分算法

用所提出的双重网格算法研究了半线性椭圆方程,其对粗网格（可以很粗）的非线性解在细网格上进行了几次线性修正．无需求解细网格上的非线性解,且重复算法最后两步,可使解的误差估计达到任意阶精度,并提出了相应的数值算例．     

非对称不定椭圆方程的两网格内罚间断有限元方法

针对一类非对称或不定椭圆方程的内罚间断有限元方法,设计和分析了相应的两网格求解算法.首先给出了内罚间断有限元解的适定性,及其在L2和间断H1范数下的先验估计;其次设计了相应的两网格求解算法,并给出算法的误差分析;最后,数值实验结果验证了算法的高效性.     

四次B样条Galerkin有限元法求解KdV方程

采用有限元方法进行空间离散,提出了解一维非线性KdV方程的四次B样条Galerkin方法.通过两个数值算例来体现这种算法的精确度, 对该方法得到的数值解与精确解以及二次B样条Galerkin有限元解进行比较,结果表明所求得的数值解与精确解符合得很好.     

Burgers方程的高精度多步显式格式

为提高Burgers方程的数值计算精度和效率，提出了一种新的高精度多步显式格式．在空间坐标上按差分法离散，在时间方向上将差分改为积分，应用显式指数时程差分法构造出了不同精度的计算格式．对不同初边值Burgers方程进行了数值模拟，并与显式交替分组法、交替Crank-Nicolson并行算法和小波法等算法进行了比较．结果表明，当新方法的网格比是参考算法网格比的2．5～20倍时，新方法数值解的绝对误差仍然小于参考算法数值解的绝对误差．该方法为数值求解非线性偏微分方程提供了一族不同精度计算格式，扩大了指数时程差分法的应用领域．     

结构动力分析的高精度有限单元法

基于哈密顿变分原理,建立了增加内部自由度的结构动力分析高精度非协调有限元法计算列式,揭示了它与动态有限元法的内在联系;同时,对增加单元结点的高精度动力有限元也进行了评述与讨论.通过实例分析,对高精度动力有限元与常规动力有限元进行了比较.算例表明,与常规有限元相比,常规动态有限元、本文的动态有限元均能给出更好的结果;高阶动力有限元能给出甚至优于动态有限元的计算结果.     

基于第二代小波的自适应有限元构造研究

为了有效地实现有限元自适应分析中函数平滑逼近和细节信息的分离，根据提出的第二代小波有限元多分辨分析方法构造了相应的自适应算法．采用第二代小波作为基函数构造了逐级嵌套的有限元逼近空间，并引入消失矩提出有限元求解方程刚度矩阵解耦的条件．函数在低分辨空间的逼近可以通过增加细节空间提升到高分辨空间，通过使刚度矩阵解耦，可以消除低分辨逼近空间和细节空间之间以及不同尺度的细节空间之间的耦合项．通过构造轴力杆的自适应分析，验证了方法的有效性，同时为实现复杂第二代小波自适应有限元分析提供了较好的研究思路、     

有限元整体刚度矩阵的一种快速解算方法

根据有限元整体刚度矩阵的特点，提出了快速解算方程的一种方法，通过对程序的运行，可以实现运用常规内存快速计算４０００个节点的有限元计算，并可实现快速解算大型稀疏对称线性方程。     

一类Poisson-Nernst-Planck方程的边平均有限元计算

针对一类三维Poisson-Nernst-Planck方程, 给出一种边平均有限元离散形式. 在适当的网格条件下, 该离散形式得到的总刚度矩阵为M-矩阵, 从而保证了数值解的非负性. 数值实验结果表明, 边平均有限元方法相比于标准有限元的CPU时间更短, 且误差较小.     

一种平面问题的加权Nitsche间断伽辽金有限元法

以经典间断伽辽金有限元法求解弹性力学界面问题,存在着由于稳定系数取值不当引起的数值不稳定问题,而加权Nitsche间断伽辽金有限元法可以缓解这种问题,但仅应用于常量单元离散的情况。为解决上述问题,基于加权Nitsche间断伽辽金有限元法,针对平面弹性力学问题,推导了四节点四边形单元离散情况下的加权系数和稳定参数的计算公式,建立了权重与稳定参数间的定性依赖关系。通过建立和求解广义特征值问题,实现了加权系数和稳定参数的自动计算,使得高阶单元的使用成为可能。通过数值试验检验了方法的收敛性和稳定性。结果表明:在求解均匀或材料分区不均匀介质问题时,加权Nitsche间断伽辽金有限元法均表现出良好的稳定性,且计算结果具有较高的精度。所提出的方法在一定程度上无须人工干预,具有高效率、高精度和良好的稳定性,可以应用于复杂界面问题。     

带源项二维浅水方程正问题和反问题的求解

考虑到求解二维浅水方程的正问题及反问题的复杂性,建立高分辨率的有限元格式对正问题进行模拟,运用最佳摄动量法对反问题进行了研究,通过对两类溃坝问题进行数值模拟及文献对比,验证了有限元方法具有较高的离散精度同时又避免了数值解的伪震荡,并运用最佳摄动量法对正问题算例进行糙率率定,模拟结果较为理想.     

有限单元法的非结构化网格及其自动生成

本文介绍了有限单元法非结构化网格的基本原理及其自动形成方法。由于有限元法是一种离散的数值求解方法，其近似求解方法的精确度，很大程度上取决于所形成网格的质量。另外，对于工程中一些形状复杂的问题，一般的网格生成方法很难对其进行离散。非结构化网格及其自动生成，使复杂形状的工程问题容易地进行离散，改善所形成网格的质量，提高近似计算的精确度，并且在有限元网格修正自适应分析中具有重要作用。     

一类二维非线性发展方程的交替方向有限体积元方法

采用交替方向有限元方法的离散思想,利用差分方法中的Douglas格式,给出了求解一类二维非线性发展方程的基于双线性插值的交替方向有限体积元法.文中对该方法进行了误差估计,并用数值例子验证了方法的有效性.     

聚合物熔体轴对称共挤出流动的有限元计算

建立数学模型，利用有限元法作数值解，确定界面，速度场及压力场，利用界面处法向速度为零，通过迭代求解界面位置，在界面处采用双结点法以保证速度和应力的连续性并满足压力的不连续性，使用ｅｎｎｅｒ迭代法求解圆形通道中两个共轴层流动的非线性方程。求出共挤出轴对称棒材的流量压降关系，计算结果与解析解比较，有很好的一致性。该数值方法可以运用于无法得到解析解的复杂流道。