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

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

基于无网格法的弹性地基加肋斜板频率数值解

为了给实际地基工程中的经济效益提供技术指标参考,应用无网格法计算加肋斜板在地基上的自由振动行为,应用移动最小二乘近似MLS和一阶剪切变形原理描述加肋斜板的位移场,分别建立斜板与肋条的势能泛函,使用Winkler弹性地基模型处理加肋板与地基之间的接触势能。通过斜板与肋条的位移协调关系寻找斜板和肋条节点参数转换公式,确立加肋斜板的自由振动控制方程。本文运用了无网格法的优势,即使肋条位置改变也不需重置网格。与有限元解的对比验证了本文方法的有效性和精确性。     

双重点最小二乘配点无网格法

配点类无网格法需要计算近似函数的二阶导数,因而在移动最小二乘(MLS)近似中至少要采用二次基函数。本文利用Voronoi图对双重点移动最小二乘近似法进行了改进,建立了基于Voronoi图的双重点移动最小二乘近似(VDG),并利用加权最小二乘法离散微分方程,导出了双重点最小二乘配点无网格法(MD GLS)。该方法将求解域用节点离散,并以节点为生成点建立Voronoi图,取Voronoi多边形的顶点为辅助点。近似函数及其二阶导数的计算过程可分解为两个步骤:首先用场函数节点值拟合辅助点处近似函数的一阶导数,再以辅助点处近似函数的一阶导数值拟合节点处近似函数的二阶导数。由于在每一步中只需计算MLS形函数及其一阶导数,这种近似方法需要较少的影响点和较小的影响域。同时借助于Voronoi结构的优良几何性质,可以快速地搜索影响点。研究表明,与基于MLS的加权最小二乘无网格法(MWLS)相比,这种方法可以显著提高计算效率,并且在精度和收敛性方面也有所改善。     

加权最小二乘无网格法

在最小二乘法和移动最小二乘近似的基础上提出了加权最小二乘无网格法．该方法除节点外又引入了一些辅助点，控制方程在所有节点和辅助点处的残差用最小二乘法予以消除，边界条件用罚函数法引入．另外对移动最小二乘近似进行了改进，并给出了最小二乘法中泛函的简化格式，因而提高了计算效率．与配点法相比，新方法精度高，稳定性好，并且系数矩阵是对称正定矩阵．与Galerkin法相比，该方法不需要进行高斯积分，因而计算量小．算例表明该方法具有效率高、精度高和稳定性好等优点，并且易于实现．     

弹性地基上矩形加肋板自由振动分析的无网格法

应用移动最小二乘无网格法研究弹性地基上矩形加肋板的自由振动问题。假设弹性地基与加肋板紧密接触,以弹簧模拟弹性地基,将弹性地基上的加肋板视为板与肋条组合的结构。基于一阶剪切理论,用无网格伽辽金法推出了板和肋条各自的动能与势能;再通过位移协调条件将两者的能量叠加,得到了弹性地基上整个加肋板的动能与势能。由Hamilton原理导出了弹性地基上加肋板自由振动的控制方程。采用完全转换法引入边界条件,求解自由振动方程,并编制了计算程序,给出了算例。将算例与ABAQUS有限元解及已有文献结果进行了比较分析,其相对误差均在5%以内,验证了该方法计算弹性地基上矩形加肋板结构自振频率的有效性。     

瞬态热传导问题的加权最小二乘无网格法

加权最小二乘无网格法是一种基于节点信息的纯无网格法,该方法使用最小二乘法建立系统的变分原理,通过移动最小二乘法构造近似函数,控制方程在节点处的残量使用最小二乘法予以消除,边界条件通过罚函数法引入。本文推导了瞬态热传导问题的加权最小二乘无网格计算格式,编制了相应的计算程序,算例结果表明,该方法具有精度高、前后处理简单的优点,是一种高效的的新型无网格法。     

基于径向基函数的加权最小二乘无网格法

径向基函数插值是一种新型的无网格插值方法,具有形式简单、空间维数无关等优点.这种插值方法具有δ函数的性质,易于满足本质边界条件,且插值函数的导数求解过程比通常的移动最小二乘插值(MLS)简单,精度也较高.另一方面,通过加权最小二乘法离散控制方程不需要积分,具有效率高,精度高等优点.本文试图将两者的优点结合起来,发展一种...     

无网格法在弹塑性问题中的应用

在弹塑性问题中引入无网格法,得到了弹塑性无网格法的控制方法;计算了中心裂纹板的塑性区分布以及极限荷载;数值结果表明了该方法的可行性和有效性。     

弹性地基加肋功能梯度板自由振动分析的无网格法

基于一阶剪切变形理论和移动最小二乘近似研究Winkler弹性地基上加肋功能梯度板的固有频率。假设功能梯度板的材料性质沿厚度方向按幂函数连续变化,基于物理中面和移动最小二乘近似分别推导功能梯度板和肋条的动能和势能,再通过引入位移协调条件,建立板和肋条节点参数转换关系,叠加两者的总能量,然后利用Hamilton原理推导加肋功能梯度板自由振动控制方程。采用完全转换法施加边界条件。通过将本文的计算结果与有限元以及文献的结果对比,验证方法的收敛性以及准确性。     

基于非线性分析的加肋板肋条位置无网格优化

在加肋板无网格模型中,肋条的位置对各种工况下加肋板受力性能的影响至关重要.文章基于一阶剪切变形和移动最小二乘法理论提出一种考虑非线性影响的加肋板无网格模型,并利用遗传算法优化肋条位置.首先,采用离散节点分别对平板和肋条进行离散,得到加肋板的无网格离散模型;其次,通过冯·卡门大挠度理论得到非矩形板几何非线性问题的弯曲控制方程;再次,通过哈密顿原理得到加肋非矩形板自由振动问题的控制方程;最后引入遗传算法,以肋条的位置为设计变量、非矩形加肋板中心点挠度最小或自振频率最大为目标函数,对肋条位置进行优化.在考虑了几何非线性影响的肋条位置优化过程中,肋条位置改变时只需重新计算位移转换矩阵,避免了网格重构.本文以全局荷载下单肋条菱形板为例与理论解进行对比,进行有效性验证.再以板的中点挠度最小和自振频率最大为优化目标,对局部荷载作用下不同形状、不同肋条布置方式的加肋板进行优化,分析方法的收敛性及稳定性.     

基于小波微分求积法的薄板弯曲分析

利用小波微分求积法(WDQM)对任意荷载作用下的薄板弯曲问题进行了求解分析。数值算例表明,小波微分求积法与一般的DQ法相比具有很好的适用性,特别是薄板受集中荷载或不连续分布荷载作用时,由于小波基函数的紧支撑特性与其对突变信号良好的描述能力,WDQ法的精度明显优于一般的DQ法,具有良好的应用前景。     

薄板弯曲分析的节点积分高阶无网格法

采用高阶无网格法求解薄板弯曲问题,在已发展的线性曲率光顺方案的基础上,通过引入泰勒展开技术,建立了能够精确再现纯弯曲和线性弯曲模式的节点积分方法。与之相比,目前无网格薄板分析主要采用的节点积分方法仅能精确再现纯弯曲模式。数值结果表明,本文方法可精确通过纯弯曲和线性弯曲试验,且能得到光滑、无振荡的弯矩场。与标准的高斯积分方法和目前已存在的节点积分方法相比,本文方法在计算精度、效率以及弯矩分布等方面均展现出显著优势。     

任意边界条件下双模量矩形薄板的弯曲

将双模量板等效为两个各向同性小矩形板组成的层合板,假定该层合板的中性面即为两个小矩形板的交界面。根据中性面上应力为零且薄板全厚度上应力的代数和为零,推导了双模量矩形薄板的中性面位置。本文采用严宗达提出的带补充项的双重正弦傅里叶级数通解,该通解可以适用于任意边界条件的矩形薄板且不需要叠加或者重新构造。联立边界条件和控制方程,求得通解中的待定系数并代入到通解中,即可得到任意边界条件下双模量矩形薄板的弯曲解析解。与有限元结果比较,本文结果符合工程精度要求。     

Large deflection analysis of rectangular plates by combined perturbation and finite strip method

The perturbation method and finite strip method are combined to solve the largedeflection bending problems of rectangular plates.Perturbation method is used to reducethe nonlinear differential equations into a series of linear differential equations.The finitestrip method is then employed to tackle these linear equations.Some calculation examplesare compared with those got by other methods.     

Bending of cantilever rectangular plates with the effect of transeverse shear deformation

On the basis of Reissner's theory, the exact solutions of the bending of cantilever rectangular plates are obtained by means of the concept of generalized simply-supported boundary. From the results obtained, it can be found that the method is valid.     

混合边界各向异性矩形板的静力分析

采用一般解析解和配点法相结合的方法,求解混合边界各向异性矩形板的弯曲问题.先由弯曲挠度的微分方程求出各种类型的齐次解和特解,然后组成一般解析解,再将板的每个边等分 为很多微小的段,仅对每一微段的中点建立应满足的边界条件,由全部边界条件方程式即可求得全部积分常数.以每边一半边界为平夹、另一半边界为简支或自由的方板为例进行了计算,并与四边均为简支的方板进行了对比,表明理论简单,结果实用.     

几何非线性分析的高效高阶无网格法

准确高效地处理几何非线性对于材料破坏等大变形过程的数值分析至关重要。考虑到无网格法具有易于形成高阶近似函数等诸多优点,本文发展了几何非线性分析的高阶无网格法。采用上一载荷步收敛的构形作为计算的参考构形,位移本质边界条件由罚函数法施加。为提高计算效率,将针对线性问题发展的二阶一致三点积分格式QC3（Quadratically Consistent 3-point integration scheme）拓展到考虑构形变化的几何非线性分析,大幅度减少了所需的积分点数目。数值结果表明,本文发展的高阶无网格法能够准确有效地处理几何非线性问题,而且在计算效率、精度以及应力场光滑性等方面均表现出显著优势。     

Bending of cantilever rectangular plate with concentrated load

In this paper the solutions for the bending of cantilever rectangular plates with concentrated load acting at any point of the middle line perpendicular to the clamped edge are given by means of a conception named modified simply supported edges and the method of superposition. Some numerical examples are presented. The total bending moment checks very well with the value determined statically.     

THE METHOD OF THE RECIPROCAL THEOREM OF FORCED VIBRATION FOR THE ELASTIC THIN RECTANGULAR PLATES (III)—CANTILEVER RECTANGULAR PLATES

In this paper, applying the method of the reciprocal theorem, we give the stationary solutions of the forced vibration of cantilever rectangular plates under uniformly distributed harmonic load and concentrated harmonic load acting at any point of the plates, the figures and tables of number value of bending moment and the deflection amplitudes as well.     

NONLINEAR BENDING OF SIMPLY SUPPORTED RECTANGULAR SANDWICH PLATES

In this paper,fundamental equations and boundary conditions of the nonlinearbending theory for a rectangular sandwicl plate with a soft core are derived by meansof the method of calculus of variations.Then the nonlinear bending for a simplysupported rectangular sandwich plate under the uniform lateral load is investigated byuse of the perturbation method and a quite accurate analytic solution is obtained.