共查询到19条相似文献,搜索用时 171 毫秒
1.
基于局部搜索算法的自然邻接点方法 总被引:9,自引:0,他引:9
自然邻接点方法(NNM)采用自然邻接点形函数进行插值,其插值形函数具有严格定义,且与
有限元形函数一样形式简洁、性能优良,因而避免了EFG法里难以准确施加位移边界条件和
材料不连续条件等诸多主要困难. 但是从形式上看自然邻接点方法仍然属于有网格的方法,
其研究和应用受到了较大的限制. 为了克服这个缺点,对于任意给定的数值积分点,提出了
一种基于局部搜索自然邻接点的寻找算法对NNM进行改进. 改进后的NNM与无单元伽辽金法
(EFG)的插值和求解过程类似,兼具有EFG的真正无网格特性及NNM的便于处理边界和材料
不连续条件等优点. 所得计算结果表明,改进后的NNM的计算精度和计算时间与NNM相当,
是一种比较理想的数值求解方法. 相似文献
2.
基于Voronoi结构的无网格局部Petrov-Galerkin方法 总被引:24,自引:2,他引:24
基于自然邻结点近似位移函数提出了一种用于求解弹性力学平面问题的无网格局部局部Petrov-Galerkin方法。这种方法在结构求解域Ω内任意布置离散的结点,并且利用需求结点的自然邻结点和Voronoi结构来构造整腐朽 求解的近似位移函数,对于构造好的近似位移函数,在局部Petrov-Galerkin方法建立整体求解的平控制方程,这样平衡方程的积分可在背景三角积分网格的形心上解析计算得到,而采用标准Galerkin方法的自然单元法需要三个数值积分点。该方法能够准确地施加边界条件,得到的系统矩阵是带状稀疏矩阵,对软件用户来说,这它学是一种安全的,真正的无网格方法,所得计算结果表明,该方法的计算精度与有限元四边界单元相当,但计算和形成系统平衡方程的时间比有限元法四边界单元提高了将近一倍,是一种理想的数值求解方法。 相似文献
3.
结合了扩展有限元法(extended finite elementmethods,XFEM)和比例边界有限元法(scaled boundary finite elementmethods,SBFEM)的主要优点,提出了一种改进型扩展比例边界有限元法(improvedextended scaled boundary finite elementmethods,$i$XSBFEM),为断裂问题模拟提供了一条新的途径.类似XFEM,采用两个正交的水平集函数表征材料内部裂纹面,并基于水平集函数判断单元切割类型;将被裂纹切割的单元作为SBFE的子域处理,采用SBFEM求解单元刚度矩阵,从而避免了XFEM中求解不连续单元刚度矩阵需要进一步进行单元子划分的缺陷;同时,借助XFEM的主要思想,将裂纹与单元边界交点的真实位移作为单元结点的附加自由度考虑,赋予了单元结点附加自由度明确的物理意义,可以直接根据位移求解结果得出裂纹与单元边界交点的位移;对于含有裂尖的单元,选取围绕裂尖单元一圈的若干层单元作为超级单元,并将此超级单元作为SBFE的一个子域求解刚度矩阵,超级单元内部的结点位移可通过SBFE的位移模式求解得到,应力强度因子可基于裂尖处的奇异位移(应力)直接获得,无需借助其他的数值方法.最后,通过若干数值算例验证了建议的$i$XSBFEM的有效性,相比于常规XFEM,$i$XSBFEM的基于位移范数的相对误差收敛性较好;采用$i$XSBFEM通过应力法和位移法直接计算得到的裂尖应力强度因子均与解析解吻合\较好. 相似文献
4.
求解不连续中厚板自由振动的微分容积单元法 总被引:2,自引:0,他引:2
基于区域叠加原理和微分容积法,发展了一种新型的数值方法——微分容积单元法,用以分析具有不连续几何特征的中厚板的自由振动。根据板的不连续情况将其划分为若干单元,在每个单元内用微分容积法将控制微分方程离散成为一组线性代数方程.在相邻的单元连接处应用位移连续条件和平衡条件,引入边界约束条件后得到一套关于各配点位移的齐次线性代数方程,由此可导出求解系统固有频率的特征方程。本文用子空间迭代法求解特征方程,并以开孔板、混合边界条件板和突变厚度板为例研究了方法的收敛性和计算精度。 相似文献
5.
本文采用了一种基于不连续场修正权函数的无网格方法来处理二维平面问题中的有限长裂纹。相较于目前常用的无网格裂纹不连续性处理方案,采用修正权函数处理裂纹附近不连续场时只需要对原权函数进行修正,算法简便易实现。本文采用基于不连续场修正权函数的无单元Galerkin方法(EFGM),对在边界上施加I-II混合型裂纹位移场的斜裂纹板进行了数值分析。并与可视性准则、衍射法和透射法等不连续准则对比了裂尖位移场、应力场和应力强度因子解的数值精度。另外,本文还对这四种不连续准则形函数的计算效率进行了分析和比较。 相似文献
6.
无网格方法与有限元法或边界元法耦合是无网格方法处理边界条件的方法之一,在无网格方法中研究无网格方法与有限元法或边界元法耦合的研究显得非常重要.本文在无单元Galerkin法和边界元法的基础上,基于无单元Galerkin法子域和边界元法子域的界面上位移连续和面力平衡条件,提出了一种新的无单元Galerkin法和边界元法的直接耦合方法,对弹性力学问题详细推导了在整个求解域上的耦合公式.与以往的耦合法相比,这种方法简单直观,不需要增加新的耦合区域,也不需要建立新的逼近函数来保证界面位移的连续性.算例结果表明,该方法具有较好的计算精度. 相似文献
7.
8.
9.
10.
11.
针对平面孔洞问题提出了一种新的数值模拟方法。本文通过水平集方法引入孔洞边界、力边界和位移边界水平集函数,利用边界水平集函数来构造边界试探项,将试探空间表示为二元幂级数与边界试探项的线性组合;同时提出一种基于水平集方法的位移边界条件施加方法,利用位移边界水平集函数来构造满足位移边界条件的近似位移场,并给出了相应的刚度矩阵和载荷矩阵表达式。与FEM、XFEM、无网格法等方法相比,该方法无需将求解域离散,具有较低的计算成本、特性良好的刚度矩阵和较为广泛的适用性。数值算例验证了该方法的有效性。 相似文献
12.
稳态蠕变是高温环境下材料的一个重要的考虑问题,它也是材料破坏的一种重要形式。由于存在划分网格,利用传统有限元法模拟稳态蠕变有一定的不足之处。作为一种新兴的数值模拟方法,无网格法不需要划分单元,只需将求解问题离散成独立的结点,计算过程可以局部细化。利用连续介质的虚功原理,将无网格伽辽金法应用于稳态蠕变的数值模拟,推导了稳态蠕变的无网格伽辽金法控制方程。利用罚参数来实现稳态蠕变的不可压缩条件和本质边界条件,能够保证求解过程中刚度矩阵的对称正定性。通过实例的计算结果表明,无网格伽辽金法在求解稳态蠕变时具有较高的计算精度,结果与理论解结果吻合,而且前后处理较为简单。 相似文献
13.
基于局部Petrov-Galerkin离散方案的无网格法 总被引:2,自引:0,他引:2
基于局部Petrov-Galerkin离散方案,选用自然邻近插值构造试函数,用Shepard函数作为权函数,提出了一种无网格方法(MNNPG),这种方法充分发挥了局部Petrov-Galerkin法的优势,并且结合了自然邻近插值的特点,方便引入边界条件,由于以Shepard函数的圆形支集作为积分子域,用分片中点插值来完成区域积分,无需额外背景网格,是一种真正的无网格法。本文将该无网格方法用于求解二维弹性力学边值问题,算例结果很好地吻合了精确解,表明该方法具有良好的数值精度和稳定性。 相似文献
14.
15.
对任意位移边界条件下的旋转周期对称结构,由拉格朗日乘子法建立有限元方程。在对称适应的坐标系下,由结构刚度矩阵的块循环性质,利用群变换给出一种新的求解方法。数值验证给出令人满意的结果。 相似文献
16.
17.
自然单元法研究进展 总被引:15,自引:2,他引:13
自然单元法是一种基于Voronoi图和Delaunay三角化几何结构,以自然邻点插值为试函数的一种新型数值方法.其既具有无网格方法和经典有限元方法的优点,又克服了两者的一些缺陷,是一种发展前景广阔的求解微分方程的数值方法.自然单元法的形函数满足插值性质,可以像有限元法一样直接施加本质边界条件,不存在基于移动最小二乘拟合的无网格方法不能直接施加本质边界条件的难题.由于自然单元法是无网格方法,可以方便处理有限元方法较难处理的一些问题,例如移动边界和大变形等问题.自然单元法与其他数值方法的最根本区别于其插值格式的不同.将自然邻点插值用于Galerkin过程,就得到基于Voronoi结构的自然单元Galerkin法.自然邻点插值有自然邻点Sibson插值和Laplace插值(非Sibson插值)两种.Laplace插值比Sibson插值在计算上要简单的多,并且不论对凸的或非凸的区域都能精确施加本质边界条件.以Laplace插值为试函数的自然单元法在数值实施上比以Sibson插值为试函数的自然单元法简单.本文对基于Voronoi结构的自然邻点插值和自然单元法的基本思想作了介绍,综述了国内外关于自然单元法的研究成果,总结了自然单元法的优点和尚需解决的问题. 相似文献
18.
不可压缩Stokes流动的PSPG无网格法 总被引:2,自引:0,他引:2
将应用于有限元法的Pressure-Stabilizing/Petrov-Galerkin(PSPG)稳定化机制引入到无网格法中,有效消除了由于速度和压力的插值模式违反LBB条件而导致的压力场的虚假振荡。采用与有限元法耦合的连续掺混法(Continuous Blending Method)施加本质边界条件,使得边界条件不仅在边界节点上而且在整条边界上都得到严格满足。给出了两个典型算例的数值模拟结果,表明了所建议无网格法模拟不可压缩Stokes流动的有效性。 相似文献
19.
一种高效的局部径向基点插值无网格方法 总被引:1,自引:0,他引:1
提出了一种弹性动力分析的高效局部径向基点插值无网格方法(MLRPI).该方法采用径向基点插值形函数近似解变量,运用局部Petrov-Galerkin法推导出了相应的离散方程,并根据波动模拟的精度要求,得到某一结点的动力方程.然后采用Newmark常平均加速度法和中心差分法相结合的显式积分格式进行时域积分,得到每个自由度的一种解耦递推格式.最后,对一平面应变问题进行了求解,比较了该文提出的解耦MI.RPI方法、常规MLRPI方法和ANSYS有限元方法的精度和计算时间,结果表明解耦MLRPI方法与常规MLRPI方法的精度相当,但计算效率大大提高. 相似文献