滑动最小二乘插值函数在板弯曲问题上的应用

 利用滑动最小二乘插值函数作为加权残值法的试函数，分析了 该试函数的拟合特性，对试函数中的基函数以及权函数的选取提出了 建议；采用最小二乘配点法求出试函数中的系数，进而可得到定解问 题的近似解；利用该试函数对薄板的挠曲、中厚板的弯曲两个例子进 行了数值计算，并与理论结果或其它数值结果进行对比，结果表明， 该试函数适用于多种边值问题，且精度高. 该法简化了选择试函数的 过程，尤其适用于工程中的各种数值计算.     

板分析的滑动最小二乘插值函数残值法

利用滑动最小二乘插值数作为加权残值法中试函数，对试函数中的基函数的以及权函数的选取提出了建议；分析了形函数的特性；对试函数拟合原函数的效果进行了分析，进而提出了权函数及相应的影响半径的取值；采用最小二北配点法求解定解问题的近似解；利用该试函数对矩形薄板和L形板的弯曲进行了数值计算，并与理论结果和有限元数值结果进行对比，结果表明，该试函数适用于多种边值问题，且精度高。该法简化了选择试验的过程，尤其适用于工程中的各种数值计算。     

紧支试函数加权残值法

将紧支函数引入加权残值法中，提出了紧支试函数加权残值法，其数值格式具有和有限元相似的窄带系数矩阵，提高了加权残值法的计算效率．在紧支试函数加权残值的基础上，导出了紧支试函数直接配点法、紧支试函数Hermite配点法和紧支试函数最小二乘配点法的具体格式，并且对几个典型算例进行了分析．与配点法相比，这些方法精度高，稳定性好，而与Galerkin法相比，这些方法效率高．     

无网格介点法求解Helmholtz方程

作为一种配点型无网格法，无网格介点MIP法具有数值实施简单、计算精度高、运算高效和适用范围广等优点。Helmholtz方程是科学与工程问题中广泛应用的一类特殊方程，因此对MIP法求解此类方程的适用性进行了验证。利用MIP法的d适应性，给出了MIP法求解该方程的两种计算格式。在数值算例中，分别对平面规则域和不规则域上的一般Helmholtz方程，以及轴对称Helmholtz方程进行了数值分析。结果表明，MIP法完全适用于求解Helmholtz方程。而且，MIP法的计算精度和收敛性都优于普通配点法。此外，MIP法的两种计算格式中，L2C0型通常具有更好的计算效果，故建议将该计算格式作为MIP法求解该类方程的标准形式。     

求解Helmholtz方程基于核重构思想的最小二乘配点法

基于核重构思想构造近似函数，将配点法和最小二乘原理相结合对微分方程进行离散， 建立了Helmholtz方程的最小二乘配点格式，并分别研究了Helmholtz方程的波传播问题和 边界层问题. 通过数值算例可以发现，给出的数值计算结果非常接近于精确解，计算精度明显高于SPH 法的数值结果，且随着节点数目的增加，其精确度越来越高，具有良好的收敛性.     

弹性薄板的广义多维微分求积分析方法

多维微分求积DC(differential cubature)法是近年迅速发展起来的一种处理多维微分方程的高精度数值方法。本文围绕多元函数插值适定性这一基本问题，探讨了实现DC法时插值结点组和插值空间应满足的条件，研究了离散结点组选择及其构造适定试函数的方法。提出了一种沿各坐标方向可以布置不同数量结点的广义交错网格DC方法，并给出了确定相应的插值空间及其选择试函数的具体方案。通过弹性薄板变形分析研究了广义交错网格DC法解决实际结构力学问题的可靠性及其潜力。研究结果表明，广义交错网格DC法较传统的交错网格DC格式具有更强的适用性和灵活性，采用数量少得多的离散结点即可达到传统交错网格DC法相同的数值精度，表现出精度高和计算工作量小的优点。     

直接增强自然单元法计算应力强度因子

自然单元法是一种新兴的无网格数值计算方法,但应用于裂纹问题计算时,其近似函数并不能准确反映裂纹尖端渐进应力场的奇异性,为获得足够的计算精度,需要在缝尖附近增大结点的布置密度。针对裂纹问题提出一种增强的自然单元法,将缝尖渐近位移场函数嵌入到自然单元法近似函数中,给出了增强试函数的构造方法,推导了总体刚度矩阵和荷载列阵的相关列式。应力强度因子可以作为附加未知量直接算得,也可用J积分或相互作用能量积分方法进行计算,对增强区域的选择和影响进行了分析。算例结果表明,基于增强自然单元法采用围线积分方法计算应力强度因子具有很高的精度,但直接以附加结点自由度形式计算则精度有所降低。     

基于参数曲面描述的层合壳体结构分析的分片Ritz法

为了对层合壳体结构进行准确和有效的分析，给出了基于分片参数曲面描述的适用于层合壳体结构分析的分片Ritz法。该方法首先将壳体结构的参考曲面进行分片处理，利用参数曲面来描述分片的壳体参考曲面并将分片按照C^0连续性进行拼接。构造壳体参考曲面分片上的具有参数曲面特征的Ritz法的试解函数，并得到相应问题的泛函和控制方程。最后给出了两个典型的层合壳体结构的数值算例，结果表明该方法准确有效。     

物理信息依赖核函数配点法的研究进展

在过去几十年里, 尽管有限元等传统计算方法已被成功用于众多科学与工程领域, 但是其在数值模拟无限域波传播、大尺寸比结构、工程反演和移动边界问题时仍面临计算量大、计算效率低、网格生成困难等计算难题. 本文介绍一类基于物理信息依赖核函数的无网格配点法及其在上述难点问题中的应用. 物理信息依赖核函数配点法的关键在于构建能反映问题微分控制方程物理信息的基函数. 基于这些物理信息依赖核函数, 该方法无需/仅需少量配点对所求微分控制方程进行离散, 即可有效提高计算效率. 本文首先介绍满足常见齐次微分方程的基本解、调和函数、径向Trefftz函数以及T完备函数等典型物理信息依赖核函数. 接着依次介绍非齐次、非均质、非稳态以及隐式微分方程构造物理信息依赖核函数的方法. 随后, 根据所求问题特点, 选用全域配点或局部配点技术, 建立相应的物理信息依赖核函数配点法. 最后, 通过几个典型算例验证所提物理信息依赖核函数配点法的有效性.      

冲击波的小波数值计算方法

基于自适应小波配点法和人工黏性技术，构造出一种简单稳定的冲击波数值计算方法。采用小波阈值滤波，生成适应流场分布的多尺度自适应网格，并利用密度场最细尺度的小波系数构造幂函数形式的冲击波定位函数，用以判断冲击波位置。联合人工黏性与冲击波定位函数，自动根据流场梯度严格控制人工黏性的大小和分布。对强/弱冲击波管问题进行计算，结果表明，该方法能够准确捕捉冲击波和有效抑制数值振荡，并且使用简单、分辨率高、计算量小。     

薄板弯曲问题的虚边界元-最小二乘法

本文提出求解任意形状的薄板弯曲问题的虚边界元-最小二乘法。本法首先利用薄板弯曲平衡方程的格林函数和离开实际边界上分布的未知的横向荷载和法向弯矩函数建立满足实际边界条件的积分方程;然后采用最小二乘法和沿虚边界分段离散化的待定的分布横向荷载和法向弯矩函数得到求上述积分方程离散化数值解的线性代数方程组。导出了一系列的数值积分的公式,并求解了许多例题,数值结果说明本法完全避免了奇异积分及其复杂的处理方法和耗时的运算,而且在边界及其附近区域解的精度比普通边界元(以后简称边界元)法大大地提高了。     

Boundary element calculation of shallow water waves

A boundary element method is proposed for studying periodic shallow water problems. The numerical model is based on the shallow water equation. The key feature of this method is that the boundary integral equations are derived using the weighted residual method and the fundamental solutions for shallow water wave problems are obtained by solving the simultaneous singular equations. The accuracy of this method is studied for the wave reflection problem in a rectangular tank. As a result of this test, it has been shown that the number of element divisions and the distribution of nodes are significant to the accuracy. For numerical examples of external problems, the wave diffraction problems due to single cylindrical, double cylindrical and plate obstructions are analysed and compared with the exact and other numerical solutions. Relatively accurate solutions are obtained.     

Static analysis of point-supported super-elliptical plates

This paper reports the static analysis of point-supported super-elliptical plates of uniform thickness subjected to a uniformly distributed lateral load. The plate perimeter was defined by a super-elliptic function with a power, corresponding to shapes ranging from an ellipse to a rectangle. The analysis was based on the Kirchhoff–Love plate theory and the computations were carried out by the Ritz method. Lagrange multipliers were used to satisfy the boundary conditions. Isotropic and homogeneous plates with 20 different shapes were examined for two distinct aspect ratios. Convergence studies were performed for the central deflection and the central bending moments. The results were checked against those of a corner-supported square plate and good agreement was obtained.     

升阶试函数族在矩形板大挠度问题中的应用

为解决加权残值法求近似解的计算精度问题，将摄动法与加权残值法相结合，首先以板中心挠度为摄动参数进行摄动，将矩形板大挠度非线性偏微分方程组分解为线性偏微分方程组，然后用最小二乘法求解．求解中构造并应用了可以由控制参数，调节的升阶试函数族，计算结果与实验结果基本一致，与以前的研究比较，计算精度明显提高．该方法对于寻求最佳试函数和最佳近似值是一种有效的方法．     

弹性静力问题的无网格弱-强形式结合法

基于改进的移动最小二乘(MLS)二阶导数近似,建立了一种求解弹性静力问题的无网格弱-强形式结合法(MLS-MWS)。该方法采用节点离散求解域,通过MLS构造形函数,将求解域划分为边界域和内部域,并分别使用控制方程的局部弱形式和强形式来建立离散系统方程。对强形式中涉及的近似函数二阶导数计算,提出了一种将其转化为求两次一阶导数的方法,与传统方法相比,该方法计算简单、精度高。MLS-MWS法结合了弱、强形式无网格法的优点,Neumann边界条件容易满足,并且只需在边界区域进行积分。文中应用该方法分析了两个弹性力学平面问题,分析结果表明本文方法具有良好的精度和收敛性。     

无单元法及其工程应用

无单元法可以求解复杂边界条件的边值问题，它只需结点信息而不需单元信息，故信息简单，特别适用于岩土工程数值分析．它的理论基础是滑动最小二乘法．本文对无单元法的基本理论作了研究，并用算例说明了研究成果．     

不连续Reissner矩形板的自由振动

用分区加权残值法研究Ressner矩形板在几何形状,边界条件等有突变时的自由振动问题,将研究对象按照结构几何形状和边界条件的具体情况划分为若干子域,在每个子域内用不同的试函数代入该域的内的控制方程到内部残值,并代入板的边界条件和各子域的协调条件得到边界残值和连续性残值,然后用最小二乘法消除残值,得到特征方程,文中讨论了该方法的收敛性和计算精度,求解了开孔矩形板的固有频率,并与已有结果进行了比较,结果表明：？该方法收敛性好,精度较高,适用范围广。     

非均质中厚板的无网格LRPIM动力学分析

用局部加权残值法建立了非均质中厚板的局部径向点插值离散系统方程,采用无网格局部径向点插值法分析了非均质中厚板的自由振动和强迫振动问题。用径向基函数耦合多项式基函数来近似试函数,用四次样条函数做为加权残值法中的权函数。所构造的形函数具有Kronecker delta性质,可以很方便地施加本质边界条件。该方法不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行。在计算过程中,取积分中的高斯点的材料参数来模拟问题域材料特性的变化。计算结果表明,利用该方法计算非均质中厚板的自由振动和强迫振动问题可以得到具有较高精度的解。