非线性抛物型方程小波Galerkin逼近

对非线性抛物型方程的初边值问题,用具有紧支集的Ｄａｕｂｃｈｉｅｓ 小波基,给出小波Ｇａｌｅｒｋｉｎ 逼近方法．由于小波基函数的特性,使所得数值方法计算量小、精度高． 数值实验表明,即使在奇异点,也能高精度求解．     

基于双正交小波基的热传导方程数值解法

对给定的热传导方程构造一个双正交小波基,从而用Ｇａｌｅｒｋｉｎ方法给出方程数值近似解的小波表示。由于小波基的双正交性、使得用Ｇａｌｅｒｋｉｎ方法得到的矩阵大大简化,因而得到一个数值求解热传导方程的简单、有效方法。     

Hermiter多项式在双参数弹性地基上圆板非线性稳定分析中的应用

双参数弹性地基上圆板的非线性稳定问题在本文中得到了研究，控制方程为Ｋａｒｍａｎ型大挠度方程，试函数由Ｈｅｒｍｉｔｅｒ多项式构成，方法采用Ｇａｌｅｒｋｉｎ加权残数法，边界条件为周边转动弹性约束，并给出了一些数值结果。     

B样条小波及在钢构件稳定分析中的应用

利用B样条小波的性质,以B样条小波尺度函数作为插值基函数,构造了基于钢构件稳定分析问题的偏微分方程数值求解方法.根据样条函数的数值特征,得到了简支梁的截面转角数值解.在不同的阶数和分解尺度下,该方法能得到不同精度、不同分辨率的结果.数值算例分析表明:与解析解相比,文中构造的B样条小波求解方法具有计算精度高、收敛比较快的特点.     

关于无网格法的几点研究

无网格伽辽金法（EFGM）是近几年发展起来的与有限元相似的一种数值算法，它采用移动的最小二乘法构造形函数，从能量泛函的弱变分形式中得到控制方程．本文讨论了无网格的两种处理本征边界条件的方法：拉格朗日乘子法和引入罚参数的方法．讨论了用不同的基函数对插值函数及对无单元法的计算精度的影响，并用算例说明了处理本征边界条件和基函数不同时的影响．     

用Galerkin—有限条法分析矩形板在非保守力作用下的稳定性

提出了一种求解板的颤振失稳问题的Ｇａｌｅｒｋｉｎ－有限条法。该方法利用加权残值法中的Ｇａｌａｒｋｉｎ方法直接从定解方程出发，建立有限条模式对问题进行求解，通过对ＧＳＦＳ板的计算结果同Ｓ．Ａｄａｌｌ的结果比较，证明了Ｇａｌｅｒｋｉｎ－有限条法是解决此类问题的一种行之有效的方法，利用Ｇａｌｅｒｋｉｎ－有限条法对边界分别为ＣＳＦＳ，ＳＳＦＳ，ＣＣＦＳ，ＣＣＦＣ的４种边界条件的矩形板在面内力和切向跟随力联     

应变梯度弹性〖WTBX〗理论C 1自然单元法

将Sibson插值作为三次单纯形Bernstein-Bézier多项式的自然邻近坐标,得到C ¹连续插值函数.将C ¹连续插值函数应用于应变梯度弹性理论.由于该函数对结点函数值和梯度值具有插值特性,应变梯度理论C ¹自然单元法可以直接施加本质边界条件.数值算例分析了双材料系统界面边界层问题和中心圆孔无限大板双轴拉伸问题,数值解与理论解吻合得较好,表明C ¹自然单元法能够用来分析应变梯度弹性理论问题.     

无网格法计算中厚板弯曲问题

无网格数值方法用于中厚平板结构的数值计算,与传统的数值方法不同,无网格数值方法只需准备节点数据,不必划分单元,减少了数据准备工作．采用最小二乘法构造形函数,从能量泛函的弱变分形式中得到控制方程,并用拉氏乘子法满足本质边界条件,从而得到偏微分方程的数值解．通过计算可见,无网格数值方法和有限元相比具有精度好、后处理方便等优点．     

层合板大挠度分析的一种新方法

提出了一种用于求解薄板结构大挠度问题的缩减基方法。分析了复合材料层合板结构在简单铰支和可移夹紧两种不同边界条件下的静态响应。数值结果表明本文方法克服了摄动方法和Ｇａｌｅｒｋｉｎ方法的不足之处，拓展了正则摄动法的使用范围，提高了Ｇａｌｅｒｋｉｎ法的有效性。     

C1自然邻近迦辽金法在偶应力理论中的应用

以non-Sibsonian插值函数作为三次单纯形Bernstein-Bézier的基坐标,构建了C1自然邻近插值函数,该函数具有二次完备性以及对结点函数值和梯度值的插值特性等性质.将C1插值函数应用于Toupin-Mindlin偶应力弹性理论,由于C1形函数的插值特性,偶应力理论迦辽金法可以直接施加本质边界条件,克服了其它无网格法施加本质边界条件的困难.具体算例包括单剪问题和中心圆孔无限大板单轴拉伸问题,数值解与理论解吻合得较好,表明C1自然邻近迦辽金法能够用来分析偶应力理论问题.     

边界点方法解Reissner型板

采用Reissner型板基本解来构建一系列特解,再通过边界点法确定边界元方程系数矩阵的全部元素,解算中不涉及具体插值,不用数值积分,避免了奇性处理,而任意点物理量的计算不依赖于待解的边界未知量,算效高,精度好。本法还可用来分析其它各类板壳问题,无论是各向同性还是各向异性的。不同的只是应按各自的基本解来构造全特解场矩阵。     

功能梯度材料的点插值无网格法

建立了一种新的求解功能梯度材料问题的点插值无网格法,这种无网格方法将径向基函数和多项式基函数耦合构造具有插值特性的近似函数,并将其应用于弹性力学问题Galerkin形式的无网格方法。在计算过程中,取高斯点的材料参数模拟功能梯度材料特性的变化,由于形函数及其导数的构造相对简单,并且满足Delta函数性质,所以该方法具有计算量小、精度高、可以像有限元法一样直接施加边界条件的优点。最后通过数值算例证明了该方法的有效性。     

自然邻近无网格Petrov-Galerkin法 求解稳态热传导问题

自然邻近无网格Petrov-Galerkin法采用自然邻近插值构造试函数,并且在由Delaunay三角形构成的多边形局部子域上采用局部Petrov-Galerkin方法建立整体求解的平衡控制方程,是一种真正的无网格法。该方法能够方便准确地施加本质边界条件,而且得到的系统矩阵是带状稀疏矩阵。对该方法在稳态热传导问题中的应用进行了研究,算例结果表明该方法具有良好的数值精度和稳定性。     

基于紧支径向基函数的局部径向点插值方法

文章对紧支径向基函数进行完备性修正,利用完备性修正的紧支径向基函数,并结合局部残差的思想,建立了局部径向点插值方法.由于该方法中的插值函数满足Delta函数性质,因此本质边界条件可以像传统的有限元方法一样容易施加,在计算过程中不需要积分网格,是一种"纯无网格方法".将该方法用于二维弹性静力问题的求解,导出其相应的离散方程.数值算例初步验证了该方法的有效性与合理性.     

周期子波在二维声辐射和声散射中的应用

提出了一种新的求解二维Helmholtz积分方程的方法。它通过将边界量用周期子波展开,将Helmhotz积分方程化为一组代数方程求解。即可求解Dirichlet、Neumann问题,也可求解合边值问题,方程的系数形成可用快速子波变换。用该方法形成的Helmholtz积分方程的系数矩阵是一稀疏矩阵,这样大大提高了计算效率,本算例表明：该方法收敛快,精度高,相同的精度下,本方法求解的未知量大大少于边界元所用未知量。     

对流扩散方程的点插值无网格算法

介绍了点插值原理,针对一维对流扩散方程构造带有多项式基的径向点插值无网格方法,证明了解的存在唯一性,并对具体对流扩散方程沿特征线进行无网格计算,计算结果表明,新算法结构简单,计算精度高。     

重心插值配点法分析梁屈曲问题

重心Lagrange插值具有数值稳定性好、计算精度高的优点。采用重心Lagrange插值多项式建立未知函数的微分矩阵,采用配点法将梁的控制方程表示为代数方程组。通过求解代数方程组系数矩阵的特征值和特征向量,求得梁的临界应力和屈曲模态。数值算例表明,文中所提出的方法具有原理简单,易于程序实现和数值计算精度高的优点。     

连续深梁的自重应力

本文讨论连续深梁在自重作用下的应力分析问题.利用体积应变函数将耦联的控制偏微分方程组解耦成为非耦联的,然后利用有限积分变换法求解所得非耦联的控制微分方程,从而得到给定问题的解答。本文所述方法可用以求解具有其他边界条件和载荷的深梁的应力分析问题.文末对连续深梁的自重应力进行了数值计算,并与有限元法结果进行了比较,二者吻合良好。     

两端固定深梁的自重应力

本文讨论两端固定深梁在自重作用下的应力分析问题,利用体积应变函数将耦联的控制偏微分方程组解耦成为非耦联的,然后利用有限积分变换法求解所得非耦联的控制偏微分方程,从而得到给定问题的解答。文中所述方法可用以求解具有其它边界条件和载荷的深梁的应力分析问题。文末对两端固定深梁的自重应力进行了数值计算,并与有限元素法的结果进行了比较,二者吻合良好。     

基于增量本构关系弹塑性分析的无网格伽辽金法

在无网格伽辽金法的基础上，采用增量形式的无网格伽辽金法的插值方法，并利用应力应变增量形式表征材料的弹塑性本构关系，在小变形假设的前提下，提出了基于增量本构关系的弹塑性分析的无网格伽辽金法；采用罚参数修正了能量变分方程式，方便地实现了无网格伽辽金法的本质边界条件．对二维平面问题及轴对称问题的计算结果表明：运用增量形式的无网格伽辽金法进行弹塑性问题的数值分析，具有较高的计算精度，且适应性强，从而为材料的弹塑性分析计算提供了一种实用方法。