用无网格局部Petrov-Galerkin法分析非线性地基梁

利用无网格局部Petrov-Galerkin法求解了非线性地基梁。在Petrov-Galerkin方法中，采用移动最小二乘（MLS）近似函数作为场主量挠度的试函数并取移动最小二乘近似函数中的体验函数作为近似场函数的加权函数，采用罚因子法施加本质边界条件。文末给出了两个计算实例，算例的结果表明，Petrov-galerkin法不仅能成功地分析线性地基梁，而且也适用于解非线性地基梁，在分析非线性地基梁时具有收敛快，稳定性好的优点。     

壳结构的无网格局部Petrov-Galerkin方法

无网格近似函数具有高度光滑性,能够很好的逼近曲壳表面及其位移场。无网格局部Petrov-Galerkin方法不论插值还是离散都不需要单元,是一种真正的无网格方法。本文基于无网格局部Petrov-Galerkin方法的基本原理,采用移动最小二乘插值,利用控制微分方程弱形式,建立了Mindlin壳结构的无网格局部Petrov-Galerkin分析方法,用屋顶壳、受夹圆柱壳、几何非线性圆柱壳作为计算实例分析了求解精度、收敛性和稳定性,并与精确解和有限元计算结果进行了对比,表明该方法计算精度高及收敛性好。     

无单元法分析薄板自由振动问题

提出了分析四边简支和四边固定薄板自由振动问题的无单元法，推导了无单元法的插值函数，从变分原理出发导出薄板自由振动的质量矩阵和刚度矩阵，编制了相应的计算程序，通过计算实例与其它方法的结果进行比较．数值结果表明无单元法具有一系列优点。     

无网格局部Petrov-Galerkin方法的并行计算研究

研究了无网格局部Petrov-Galerkin方法MLPG(Meshless Local Petrov-Galerkin Method)的并行算法与并行实现过程。将MLPG方法推广到弹性动力学问题,研究了MLPG方法中节点搜索、积分点搜索、数值积分及方程组求解等过程的并行算法,并给出了MLPG方法并行计算的具体实现过程。两个数值算例验证了MLPG并行算法的有效性;计算结果表明,MLPG方法的并行计算具有很好的并行性能和可扩展性。     

无网格局部Petrov-Galerkin法求解板壳弹塑性大变形

无网格局部Petrov-Galerkin法构造的高阶光滑的形函数非常适合建立板壳结构场函数的逼近函数,是一种比较理想的研究板壳问题的方法。基于Mindlin板壳理论,采用更新拉格朗日原理和大变形条件下场量的无网格表达形式,实现了率型无网格局部Petrov-Galerkin方法对板壳弹塑性大变形的求解,算例分析表明了方法的有效性和较高的分析精度。     

边界元分析薄板自由振动问题的一个近似方法

本文提出边界元法分析域内具有支承及集中质量的薄板自由振动问题的近似方法，该方法在利用基本解的基础上，将域内积分化为边界积分来处理，节省了工作量，文中计算实例结果表明，该方法的精度满足实际工程的要求。．     

改进广义移动最小二乘近似的无网格法

无网格法是求解微分方程定解问题的一种新数值方法.移动最小二乘近似只要求近似函数在各节点处的误差的平方和最小,对近似函数导数的误差没有任何约束.而广义移动最小二乘近似要求近似函数及其导数在所有节点处的误差的平方和最小.为了降低计算工作量,本文构造了要求近似函数在全部节点处和任意阶导数在部分节点处误差的平方和最小的改进广义移动最小二乘近似.数值计算显示本文提供的方法关于函数值和各阶导数值都具有很高的精度.     

用杂交边界点法分析简支薄板的热弯曲问题

该文利用杂交边界点法对简支薄板的热弹性弯曲进行了分析计算.采用薄板的热弹性理论,通过薄板的修正变分原理建立了各向同性薄板的边界局部积分方程,域内变量使用基本解插值,而边界上的变量则用移动最小二乘法近似.计算时仅需边界上离散点的信息,无论变量近似还是数值积分都不需要网格,因此该方法是一种纯边界类型无网格方法.数值算例表明,杂交边界点法在分析薄板的热弯曲问题时具有效率高、精度高和收敛性好等优点.     

无网格方法中本质边界条件的处理

基于新的离散模定义，采用移动最小二乘近似方法，给出了场变量的近似表达形式；从约束变分原理出发，给出了无网格方法中新的本质边界条件的处理方案——罚函数法；对权函数、罚函数的选择对计算精度的影响进行了研究。数值计算结果表明该方法不仅简单合理，而且具有较高的精度和稳定性。     

弹性力学问题的局部边界积分方程方法

提出了弹性力学平面问题的局部边界积分方程方法。这种方法是一种无网格方法,它采用移动最小二乘近似试函数,且只包含中心在所考虑节点的局部边界上的边界积分。它易于施加本质边界条件。所得系统矩阵是一个带状稀疏矩阵。它组合了伽辽金有限元法、整体边界元法和无单元伽辽金法的优点。该方法可以容易推广到求解非线性问题以及非均匀介质的力学问题。计算了两个弹性力学平面问题的例子,给出了位移和能量的索波列夫模,所得计算结果证明：该方法是一种具有收敛快、精度高、简便有效的通用方法。     

Local Petrov-Galerkin method for a thin plate

The meshless local Petrov-Galerkin (MLPG) method for solving the bendingproblem of the thin plate were presented and discussed. The method used the moving least-squares approximation to interpolate the solution variables, and employed a local symmetricweak form. The present method was a truly meshless one as it did not need a finite elementor boundary element mesh, either for purpose of interpolation of the solution, or for theintegration of the energy. All integrals could be easily evaluated over regularly shapeddomains ( in general, spheres in three-dimensional problems ) and their boundaries. Theessential boundary conditions were enforced by the penalty method. Several numericalexamples were presented to illustrate the implementation and performance of the presentmethod. The numerical examples presented show that high accuracy can be achieved forarbitrary grid geometries for clamped and simply-supported edge conditions. No postprocessing procedure is required to computer the strain and stress, since the originalsolution from the present method, using the moving least squares approximation, is already smooth enough.     

基于无网格局部Petrov-Galerkin(MLPG)法求解任意磁场方向下的定常MHD流动

磁流体流动在现代工业和科研中有着广泛的应用,但磁流体的流动受到磁场的影响,与一般流体区别较大,需要对其进行深入的研究。磁流体的流动受到流体力学流动方程和麦克斯韦方程的共同影响,其精确解在有限条件下才能得到,因此对磁流体的流动进行数值模拟具有重要的意义。本文采用移动最小二乘法计算形函数,利用无网格局部Petrov-Galerkin(MLPG)法得到控制方程的离散形式,在管壁为任意电导率及任意方向外加磁场的条件下,对方形直管道中定常流动的磁流体进行了数值计算。MLPG法的计算是基于节点的,不需要任何网格或单元,是一种真正的无网格方法。计算结果与Scheriff精确解进行了比较,表明该方法适用于中等以下哈特曼数的磁流体流动计算。     

粘弹性结构自由振动分析

提出一种分析粘弹性结构自振特性的数值方法,通过弹性等效空间特征值问题及一个代数方程的求解,得出粘弹性结构的频率及振型。此外建立了薄板的粘弹性动力方程并给出相应的求解方法。文中做了数值验证。     

带源参数热传导问题的基于滑动Kriging插值的MLPG法

利用基于滑动Kriging插值的无网格局部Petrov-Galerkin(MLPG)法来求解带源参数的二维热传导问题,推导了相应的离散方程。由于滑动Kriging插值法构造的形函数满足Kronecker Delta特性,因此可以直接施加本质边界条件。在离散过程中采用Heaviside分段函数作为局部弱形式的权函数,时间域则通过向后差分法进行离散,这一处理过程中刚度矩阵只涉及到边界积分,而没有涉及到区域积分。最后通过算例验证了本方法的有效性。     

薄板振动特征值问题的一种改进积分方程解法

本文根据一种改进的边界元/有限元混合法求解薄板振动固有频率问题,既避开了标准的边界元法所导致的求解非代数特征值方程的困难,亦能够基本上消除通常的边界元/有限元混合法结果精度受区域内部单元划分影响较大的弊端。文中讨论了迭代算法的收敛问题,并用于薄板固有频率分析。数值结果表明,即便是在域内单元很粗疏划分的情况下,本文的方法仍能给出相当满意的结果。     

无单元法求解正交各向异性板自由振动问题

无单元法是一种新出现的数值方法。本文对无单元法的数学基础—滑动最小二乘法进行了详细的研究,推导了无单元法的形函数,并对一些关键问题,如权函数的选取,正交基函数,边界条件,数值实现方法等得出了研究结论。用无单元法研究了正交各向异性板的自由振动问题,由动力学变分原理和滑动最小二乘法导出了正交各向异性板的无单元法质量矩阵和刚度矩阵,编制了相应的计算程序,通过计算实例验证了该方法的有效性。     

The stokes flow of an arbitrary prolate axisymmetric body towards an infinite plane wall

By distributing continuously the image Sampsonlets with respect to the plane and by applying the constant density, the linear and the parabolic approximation, the analytic expressions in closed form for flow field are obtained. The drag factor of the prolate spheroid and the Cassini oval are calculated for different slender ratios and different distances between the body and the plane. It is demonstrated that the proposed method is satisfactory both in convergence and accuracy. Comparison with existing results in the case of prolate spheroid shows that the coincidence is quite good.