自然单元法在Winkler地基薄板计算中的应用

自然单元法是一种基于Voronoi图及Delaunay三角形剖分图,以自然邻接点插值为试函数的一种无网格数值方法.本文以目前该方法中自然邻接点的Laplace插值形函数为基础,求出了其一阶及二阶导函数,建立了Winkler地基上薄板弯曲挠度的自然单元法求解控制方程,并编制了相应的计算程序,通过算例分析表明了本文方法的可行性和有效性.     

基于局部搜索算法的自然邻接点方法

自然邻接点方法(NNM)采用自然邻接点形函数进行插值,其插值形函数具有严格定义,且与有限元形函数一样形式简洁、性能优良,因而避免了EFG法里难以准确施加位移边界条件和材料不连续条件等诸多主要困难. 但是从形式上看自然邻接点方法仍然属于有网格的方法,其研究和应用受到了较大的限制. 为了克服这个缺点,对于任意给定的数值积分点,提出了一种基于局部搜索自然邻接点的寻找算法对NNM进行改进. 改进后的NNM与无单元伽辽金法(EFG)的插值和求解过程类似,兼具有EFG的真正无网格特性及NNM的便于处理边界和材料不连续条件等优点. 所得计算结果表明,改进后的NNM的计算精度和计算时间与NNM相当,是一种比较理想的数值求解方法.     

自然单元法研究进展

自然单元法是一种基于Voronoi图和Delaunay三角化几何结构,以自然邻点插值为试函数的一种新型数值方法.其既具有无网格方法和经典有限元方法的优点,又克服了两者的一些缺陷,是一种发展前景广阔的求解微分方程的数值方法.自然单元法的形函数满足插值性质,可以像有限元法一样直接施加本质边界条件,不存在基于移动最小二乘拟合的无网格方法不能直接施加本质边界条件的难题.由于自然单元法是无网格方法,可以方便处理有限元方法较难处理的一些问题,例如移动边界和大变形等问题.自然单元法与其他数值方法的最根本区别于其插值格式的不同.将自然邻点插值用于Galerkin过程,就得到基于Voronoi结构的自然单元Galerkin法.自然邻点插值有自然邻点Sibson插值和Laplace插值(非Sibson插值)两种.Laplace插值比Sibson插值在计算上要简单的多,并且不论对凸的或非凸的区域都能精确施加本质边界条件.以Laplace插值为试函数的自然单元法在数值实施上比以Sibson插值为试函数的自然单元法简单.本文对基于Voronoi结构的自然邻点插值和自然单元法的基本思想作了介绍,综述了国内外关于自然单元法的研究成果,总结了自然单元法的优点和尚需解决的问题.      

动力弹塑性分析的无网格自然单元法

基于无网格自然单元法,提出了结构动力弹塑性响应分析的一条新途径.自然单元法是一种新兴的无网格数值计算方法,其实质是基于自然邻近插值的伽辽金法.自然单元法在本质边界条件的施加上较采用移动最小二乘法的无网格法具有明显的优势.在空间域上采用自然单元法离散,并运用加权余量法推导了动力弹塑性分析的离散控制方程.然后,采用预校正形式的Newmark法在时间域上进行求解.最后给出了数值算例,并验证了所提方法的有效性和正确性.     

基于Lasserre算法的自然单元法形函数计算

基于Lasserre体积算法推导了两种插值方案下自然单元法形函数及其导数的具体计算方法,特别是对计算点处于某些特殊位置时可能造成计算失败的原因和处理方法进行了较为深入的研究. 算例结果验证了Sibson与non-Sibson插值形函数在三角形外接圆的圆周上具有不同的连续性,自然单元法形函数在凸区域的边界结点间是线性变化的,因而可以方便地施加本质边界条件.     

自然单元法在线弹性断裂力学中的应用研究

采用Laplace插值构造近似场函数,代入弹性力学边值问题的变分弱形式,由变分原理推导出自然单元法的离散控制方程。运用线弹性断裂力学理论,分析了有限板单边Ⅰ型裂纹的应力强度因子和有限板拉剪复合裂纹的扩展,给出了分析线弹性断裂力学问题的自然单元法。数值计算结果表明了方法的正确性和有效性。     

自然单元法的自适应研究

自然单元法(NEM)是一种新兴的无网格数值计算方法,具有前处理简单和易于准确施加本质边界条件等优点.本文基于Z-Z后验误差估计方法,给出了一种自然单元法的误差估计因子和自适应分析细化方案.采用结点处的光滑应变计算相应的恢复应力,并用于构造全域上的恢复应力场.通过结点Voronoi单胞内的能量范数相对误差指示需要进行结点...     

基于自然单元法的极限上限分析

自然单元法是一种基于离散点集的Voronoi图和Delaunay三角化几何信息,以自然邻近插值为试函数的新型数值方法.相对于一般无网格法中常采用的移动最小二乘近似而言,自然邻近插值不涉及到复杂的矩阵求逆运算,更不需要任何人为的参数,可以提高计算效率.采用该方法构造的形函数满足Delta函数的性质,可以像有限元一样准确地施加边界条件,可以方便处理场函数及其导数的不连续性的问题.论文将自然单元法应用到极限上限分析中,编制了相应的计算程序,通过极限分析的几个经典算例进行了验证,同时采用类似于分片应力磨平的方式,编制相应的磨平程序,由计算点上的塑性耗散功外推得到了节点上的塑性耗散功的值,从而画出了极限状态下结构的塑性耗散功的分布云图.计算结果表明采用自然单元法求解极限上限分析具有稳定性好,精度高,收敛快等优点.     

基于自然单元法的极限下限分析

自然单元法是一种基于自然邻近插值的无网格方法,其形函数构造简单,且易于施加本质边界条件．基于极限分析理论和降维思想,建立了用自然单元法进行二维理想弹塑性结构极限下限分析的整套求解算法．极限下限分析中的弹性应力场由自然单元法直接求出,利用自然单元法弹塑性迭代的结果得到自平衡应力基矢量并构造出所需的自平衡应力场．极限下限问题被转化为一系列数学规划子问题,并采用复合形法直接求解．数值算例结果表明,将自然单元法应用于极限下限分析是可行有效的．     

Helmholtz方程的无网格自然邻接点Petrov-Galerkin法

采用无网格自然邻接点Petrov-Galerkin法求解Helmholtz方程。通过自然邻接点插值构造试函数,并采用有限元法的三角形线性单元的形函数作为加权残值法的权函数,基于局部Petrov-Galerkin法建立了Helmholtz问题的离散方程。由于所构造的形函数满足KroneckerDelta性质,因此本质边界条件的施加十分方便。数值算例表明,基于无网格自然邻接点Petrov-Galerkin法的计算结果非常接近精确解,且随着节点的增加,其精确度越来越高,验证了本文方法具有良好的收敛性。     

加锚体自然单元法分析与程序设计

建立了加锚体的自然单元法分析模型.结合自然相邻插值和一维轴力杆单元的劲度矩阵,推导了锚杆支护数值模拟的单元列式,给出了自然单元法分析加锚体的实施步骤和主要程序设计.结构的离散和锚杆的离散完全独立,因此可以自由地加设或拆除锚杆.算例分析表明自然单元法分析加锚体的精度与有限元的相当,但比有限元的实施要方便得多.     

高度不规则网格多边形单元的有理函数插值格式

借鉴自然邻点插值法,提出了基于高度不规则网格多边形单元的有理函数插值格式一多边形有理函数插值.给出了多边形有理函数插值形函数的计算表达式.该插值格式以多边形的顶点作为插值点,插值形函数为有理函数形式,克服了传统有限方法中构造边数大于4单元多项式形式位移插值的困难.     

基于局部Petrov-Galerkin离散方案的无网格法

基于局部Petrov-Galerkin离散方案,选用自然邻近插值构造试函数,用Shepard函数作为权函数,提出了一种无网格方法(MNNPG),这种方法充分发挥了局部Petrov-Galerkin法的优势,并且结合了自然邻近插值的特点,方便引入边界条件,由于以Shepard函数的圆形支集作为积分子域,用分片中点插值来完成区域积分,无需额外背景网格,是一种真正的无网格法。本文将该无网格方法用于求解二维弹性力学边值问题,算例结果很好地吻合了精确解,表明该方法具有良好的数值精度和稳定性。     

多边形有限元研究进展

有限元法是数值求解偏微分方程边值问题的重要方法,采用不规则多边形单元网格, 可以方便有效地模拟材料的力学性能, 又使得区域网格剖分变得灵活方便. 特别是对于复杂的几何形状, 多边形单元网格具有更大的优势. 本文对国内外有关多边形有限元法的最新进展作了初步的总结和评述, 主要以基于位移法的多边形有限元为主.论述了多边形有限元的发展历史, 给出了多边形单元上的Wachspress插值、Laplace插值和重心坐标的一些最新研究成果. 与经典有限元法形函数为多项式形式不同, 多边形单元的形函数为有理函数或者无理函数形式. 多边形单元插值形函数满足线性完备性, 可以再现线性位移场, 像经典有限元法一样直接施加本质边界条件; 插值函数在多边形的边界上是线性的,确保不同单元间的自动协调. 不同单元的插值形函数表达公式形式统一, 方便混合单元网格计算的程序编写. 提出了多边形有限元法今后需要研究的问题.