首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 117 毫秒
1.
基于单位分解法的无网格数值流形方法   总被引:20,自引:1,他引:19  
李树忱  程玉民 《力学学报》2004,36(4):496-500
在数值流形方法和单位分解法的基础上,提出了无网格数值流形方法. 无网格数值流形 方法在分析时采用了双重覆盖系统,即数学覆盖和物理覆盖. 数学覆盖提供的节点形成求解 域的有限覆盖和单位分解函数;而物理覆盖描述问题的几何区域及其域内不连续性. 与原有 的数值流形方法相比,无网格数值流形方法的数学覆盖形状更加灵活,可以用一系列节点的 影响域来建立数学覆盖和单位分解函数,具有无网格方法的特性,从而摆脱了传统的数值流 形方法中网格所带来的困难. 与无网格方法相比,由于采用了有限覆盖技术,试函数的构造 不受域内不连续的影响,克服了原有的无网格方法在处理不连续问题时所遇到的困难. 详细推导了无网格数值流形方法的试函数和求解方程,最后给出了算例,验证了该方法的正 确性.  相似文献   

2.
借鉴流形方法思想,引入广义节点的概念,对传统的无网格法进行了改进,建立了可具有任意高阶多项式插值函数的广义节点无网格方法。同时采用径向插值函数构造具有插值特性的逼近函数;采用配点法建立系统的离散方程。在阐述了这种方法基本原理的同时,针对线弹性力学问题给出了这种方法的数值计算列式。与传统无网格方法相比,这种方法更具有一般性;同时由于采用了配点法而不需要背景积分网格,所以可以认为这种方法是某种真正意义上的无网格法。当选取0阶广义节点位移插值函数时便可得到传统的无网格法;在不增加支持域内节点数目的条件下,通过选取高阶广义节点位移插值函数可以提高计算精度。最后通过算例分析,对0阶、1阶及2阶广义节点无网格法与现有的有关解答进行了对比,论证了其合理性。  相似文献   

3.
数值流形方法在流固耦合谐振分析中的应用   总被引:2,自引:1,他引:1  
数值流形方法(流形法)是石根华博士利用现代数学中流形分析的有限覆盖技术建立起来的新的数值分析方法,统一解决了连续和非连续变形的力学问题,具有广阔的应用前景。本文将流形法应用于交界面耦合的流固振动分析,采用平面矩形数学网格,针对无粘、无旋、不可压缩流体和无阻尼的固体结构,提出分析流固耦合系统简谐振动的高阶流形法公式,其中,采用拉格朗日乘子法引入流场的已知边界条件。本文还初步研究了在特殊的无限远流场中采用解析解覆盖函数的实现技术。文中算例体现了流形法网格划分的方便性和计算的高精度,显示出流形法在数值解和解析解联合运用上的优势。  相似文献   

4.
数值流形方法研究及应用进展   总被引:2,自引:0,他引:2  
基于有限覆盖技术的数值流形方法是一种新的广义的数值方法.该方法的场函数近似原理和有限元、无网格、单位分解等方法相似,但在网格划分、覆盖形式、近似函数等方面有其自身的特点和优势.对该方法近年来在理论研究和应用方面取得的重要进展进行了综述.在理论研究方面,日前已对不同形式物理覆盖流形单元的性能进行了研究,结果表明流形单元的精度较有限单元高,且提高覆盖函数的阶次能提高单元的精度;同时理论研究已由二维低阶流形方法推广到三维高阶流形方法,由线性流形方法推广到非线性流形方法,由基于能量原理的流形方法推广到基于加权余量的流形方法,非协调流形方法、无网格流形方法等也已开展了研究;此外,覆盖系统的自动生成、覆盖函数的形式以及边界条件的处理方法等流形方法相关理论的研究也取得了进展.在应用方面,开展了有关岩石破坏和裂纹扩展等非连续变形分析更深入的研究,并已逐步推广到金属塑性变形分析、多孔介质变形分析以及温度场的数值分析等多个领域.针对日前流形方法的研究和应用现状,该文展望了流形方法理论及实现方法的研究方向、及其在计算流体力学、金属成形等大变形问题、多物理场分析等领域的应用前景.  相似文献   

5.
用径向基函数构造无网格点插值法的形函数,插值函数具有Kronecker delta函数性质,因此可以很方便地施加本质边界条件.利用无网格局部径向点插值方法分别对一个对边固支另对边简支中厚板和一个悬臂中厚板的弯曲进行了分析计算.该方法不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行,是一种真正的无网格方法.算例表明:将无网格局部径向点插值法应用于计算中厚板的弯曲问题,所求得的位移场和应力场都是光滑的;在径向基函数的基础上,附加多项式大大提高了插值精度;所得结果与弹性力学理论解以及有限元解都十分吻合.  相似文献   

6.
数值流形方法 (numerucal manifold method,NMM)通过引入数学覆盖和物理覆盖两套系统来统一处理连续和非连续问题.通过用移动最小二乘插值(moving least squares interpolation,MLS)中的节点影响域构造数学覆盖,得到了基于数值流形方法的无网格伽辽金法(element free Galerkin,EFG).该方法在保证前处理简单的同时,又能方便处理如裂纹等不连续问题.建立了适用于小变形和大变形的裂纹扩展计算格式,并通过对曲折裂纹(kinked crack)的处理,在不加密的情况下实现了任意小步长的裂纹扩展,大大提高了在固定网格中模拟裂纹扩展的实用性.大小变形的结果对比表明,按照不考虑构型变化的小变形计算,结果可能偏于危险.  相似文献   

7.
弹性力学中的一种非协调数值流形方法   总被引:1,自引:0,他引:1  
魏高峰  冯伟 《力学学报》2006,38(1):79-88
通过引入数学和物理双重网格,将插值域与 积分域分别定义在不同的覆盖上,即在数学网格上进行插值函数的构造,物理网格上完成 系统能量泛函积分运算,最后通过覆盖权函数将二者联结在一起. 它的优点是单元网格划 分随意,不受复杂边界形状和二相材料界面的限制,单元可以是任意形状,是较之于有限 元方法更一般的数值模拟方法. 在4节点四边形数值流形方法中,由于单元总体位移函数 包含的完全多项式不完全,使得计算精度不够精确,为此,在单元总体位移函数上附 加非协调位移基本项,使之趋于完全,提出了弹性力学问题的一种改进的数值流形 方法------非协调数值流形方法. 通过内部自由度静力凝聚处理,导出了消除内参后的单元应变矩阵 和单元刚度矩阵,使得在不增加广义节点自由度的前提下,大大提高了数值流形方法的计 算精度和计算效率. 同时对非协调项进行了显式处理,可以对工程实践起到更切实的帮助. 数值试验表明,它们能够保证收敛,有较高的精度,对畸变不敏感,从而证明了该方法的 可行性.  相似文献   

8.
数值流形方法是一种非常灵活的数值计算方法,连续体的有限单元方法和块体系统的非连续变形分析方法只是这一数值方法的特例.数值流形方法中高阶位移函数的构造可通过提高权函数的阶次来实现,这种方法往往需要沿单元边界配置适当的边内节点,这些结点的出现增加了前处理的复杂性,特别是对于大型复杂的空间问题.另一方面,在数值流形方法中可通过缩小单元尺寸(h加密)来提高求解精度.当模拟裂纹扩展时,这种细化策略可用来克服裂纹尖端的奇异性.一个传统的解决方案是细化整个网格,但这会导致计算效率的显著降低.将适合分析的T样条(analysis-suitable T-spline,AST)引入数值流形方法中来建立高阶数值流形方法的分析格式,有效的避免了该问题的出现.AST样条基函数具有线性无关,单位分解,局部加密等许多重要性质,使得其非常适合用于工程设计及分析.在引入AST样条后,可通过改变数学覆盖的构造形式建立不同阶次的数值流形方法分析格式;AST样条自身的局部加密性质也使得数值流形方法中的数学网格局部加密更容易实现.算例结果表明:随着AST样条基函数阶次的提高,数值流形方法的计算结果有了明显的改善;基于AST样条基函数的数值流形方法在保持计算精度的前提下降低了自由度的数量.  相似文献   

9.
无网格法因为不需要划分网格, 可以避免网格畸变问题,使得其广泛应用于大变形和一些复杂问题. 径向基函数配点法是一种典型的强形式无网格法,这种方法具有完全不需要任何网格、求解过程简单、精度高、收敛性好以及易于扩展到高维空间等优点,但是由于其采用全域的形函数, 在求解高梯度问题时 存在精度较低和无法很好地反应局部特性的缺点. 针对这个问题,本文引入分区径向基函数配点法来求解局部存在高梯度的大变形问题. 基于完全拉格朗日格式,采用牛顿迭代法建立了分区径向基函数配点法在大变形分析中的增量求解模式.这种方法将求解域根据其几何特点划分成若干个子域, 在子域内构建径向基函数插值, 在界面上施加所有的界面连续条件,构建分块稀疏矩阵统一求解. 该方法仍然保持超收敛性, 且将原来的满阵转化成了稀疏矩阵, 降低了存储空间,提高了计算效率. 相比较于传统的径向基函数配点法和有限元法, 这种方法能够更好地反应局部特性和求解高梯度问题.数值分析表明该方法能够有效求解局部存在高梯度的大变形问题.   相似文献   

10.
原有数值流形方法通过积累每一时步的小变形而得到结构最终的大变形,然而,当结构发生大变形、大转动时往往产生较大计算误差.针对该问题,从动量守恒方程以及应力边界条件的积分弱形式出发,引入流形方法的插值函数,建立了基于有限变形理论的数值流形方法.通过对比改进前后流形方法的计算迭代格式,指出了原有流形方法计算大变形问题时的误差来源.最后,通过大变形悬臂梁和旋转块体算例对有限变形流形方法进行了验证.数值结果表明,改进后的流形方法能够很好地处理大变形大转动问题,消除了转动所带来的计算误差,其计算结果与解析解及ABAQUS软件求得的数值解相吻合.  相似文献   

11.
将比例边界法与无单元伽辽金法相结合,建立了反平面断裂分析的无单元伽辽金比例边界法。这是一种边界型无网格法,在环向方向上采用无单元伽辽金法进行离散,因此计算时仅需要边界上的节点信息,不需要边界元所要求的基本解。为了便于施加本质边界条件,通过建立节点值和虚拟节点值之间的关系给出了修正的移动最小二乘形函数。在径向方向上,该方法利用解析的方法求解,因此是一种半解析的数值方法。最后,给出了数值算例,并验证了所提方法后处理简单和计算精度高的特点,适合于求解反平面断裂问题。  相似文献   

12.
王涛  陈莘莘 《力学季刊》2021,42(3):507-516
作为一种最近发展起来的半解析数值方法,插值型无单元伽辽金比例边界法不仅无需基本解,且在处理应力奇异性问题和无限域问题时十分有效.为了更有效地求解粘弹性问题,对插值型无单元伽辽金比例边界法应用于此类问题进行了研究,并发展了相应的算法. 通过时域分段展开,将时空耦合的初边值问题转化为一系列递推形式的边值问题,然后采用插值型无单元伽辽金比例边界法进行自适应计算.在径向保持解析特性的基础上,环向采用无单元伽辽金法离散可简化前处理和后处理工作量.此外,改进的插值型移动最小二乘法形函数具有插值性,有效地解决了本质边界条件不能直接施加的困难.最后给出了数值算例,并验证了所提方法的有效性和正确性.  相似文献   

13.
宋彦琦  周涛 《力学学报》2018,50(4):853-862
S-R(strain-rotation)和分解定理克服了经典有限变形理论的一些缺点, 使其可以为几何非线性数值分析提供可靠的理论基础. 对于大变形问题, 由于无网格法(element-free method)避免了对单元网格的依赖, 从而从根本上避免了有限单元法(finite element method, FEM)的单元畸变问题, 保证了求解精度. 因此, 将无网格法和S-R和分解定理结合起来势必能建立一套更加合理可靠的几何非线性数值计算方法. 目前基于S-R 定理的无网格数值方法研究较少并且只能用于二维平面问题的求解, 但实际上绝大多数问题都必须以三维模型来进行处理, 因此建立适用于三维情况的S-R无网格法是非常有必要的. 本文给出了适用于三维情况的S-R 无网格法: 采用由更新拖带坐标法和势能率原理推导出来的增量变分方程, 利用基于全局弱式的无网格Galerkin 法(EFG)得到了用于求解三维空间问题的离散格式. 利用MATLAB编制三维S-R 无网格法程序, 对受均布载荷的三维悬臂梁和四边简支矩形板结构的非线性弯曲问题进行了计算. 最后将所得的数值结果与已有文献进行了比较, 验证了本文的三维S-R无网格数值算法的合理性、有效性和准确性. 本文的三维S-R无网格数值算法可以作为一种可靠的三维几何非线性数值分析方法.   相似文献   

14.
张赞  程玉民 《力学季刊》2007,28(2):333-339
无网格方法与有限元法或边界元法耦合是无网格方法处理边界条件的方法之一,在无网格方法中研究无网格方法与有限元法或边界元法耦合的研究显得非常重要.本文在无单元Galerkin法和边界元法的基础上,基于无单元Galerkin法子域和边界元法子域的界面上位移连续和面力平衡条件,提出了一种新的无单元Galerkin法和边界元法的直接耦合方法,对弹性力学问题详细推导了在整个求解域上的耦合公式.与以往的耦合法相比,这种方法简单直观,不需要增加新的耦合区域,也不需要建立新的逼近函数来保证界面位移的连续性.算例结果表明,该方法具有较好的计算精度.  相似文献   

15.
孟智娟  迟晓菲 《力学季刊》2022,43(2):355-365
本文以求解三维波动方程为例,介绍了改进的插值型维数分裂无单元Galerkin方法,推导了方程的弱形式,构造了具有插值特性的逼近函数,建立了可直接施加本质边界条件的离散方程组,研究不同本质边界条件施加方法对计算结果的影响.本文列举了三种常用的处理本质边界条件的方法:直接配点法、对角元素置大数法和对角元素化一法.选取了三个数值算例,分别采用不同的本质边界条件施加方法,分析计算结果,证明了三种施加方法的有效性,讨论了每种施加方法的优缺点,并针对问题需求选出合适的施加本质边界条件的方法.与改进的无单元Galerkin方法相比,改进的插值型维数分裂无单元Galerkin方法具有更高的计算精度和更快的计算速度.  相似文献   

16.
One of major difficulties in the implementation of meshfree methods using the moving least square (MLS) approximation, such as element-free Galerkin method (EFG), is the imposition of essential boundary conditions as the approximations do not pass through the nodal parameter values. Another class of meshfree methods based on the radial basis point interpolation can satisfy the essential boundary conditions exactly since its approximation function passes through each node in an influence domain and thus its shape functions possess the properties of delta function. In this paper, a coupled element-free Galerkin(EFG)-radial point interpolation method (RPIM) is proposed to enhance their advantages and avoid their disadvantages. Discretized equations of equilibrium are obtained in the RPIM region and the EFG region, respectively. Then a collocation approach is introduced to couple the RPIM and the EFG method. This method satisfies the linear consistency exactly and can maintain the stiffness matrix symmetric. Numerical tests show that this method gives reasonably accurate results consistent with the theory.  相似文献   

17.
The finite element method is used to find the elastic strain (and thus the stress) for given velocity fields of the Leonov model fluid. With a simple linearization technique and the Galerkin formulation, the quasi-linear coupled first-order hyperbolic differential equations together with a non-linear equality constraint are solved over the entire domain based on a weighted residual scheme. The proposed numerical scheme has yielded efficient and accurate convective integrations for both the planar channel and the diverging radial flows for the Leonov model fluid. Only the strain in the inflow plane is required to be prescribed as the boundary conditions. In application, it can be conveniently incorporated in an existing finite element algorithm to simulate the Leonov viscoelastic fluid flow with more complex geometry in which the velocity field is not known a priori and an iterative procedure is needed.  相似文献   

18.
基于径向基函数的无单元法求解力学问题误差分析   总被引:1,自引:0,他引:1  
径向基函数形状参数的选择在无单元法数值计算中一直是一个热门的问题,现在已总结出许多确定形状参数的经验公式. 但还没有相关研究表明这些形状参数是如何随着影响域尺寸而变化的. 本文研究了MQ(multi-quadrics) 径向基函数中形状参数对无单元法计算误差的影响. 首先,从理论上分析了形函数导数随着形状参数值的变化趋势,和以计算点为中心节点对称布置与不对称布置的形函数导数的变化规律;然后分析了影响域尺寸对误差的影响,得到了在不同影响域尺寸下,误差随形状参数值变化的规律;在此基础上,给出了影响域范围值.  相似文献   

19.
混合变换法在无网格伽辽金方法中的应用   总被引:1,自引:0,他引:1  
移动最小二乘近似的非插值特征给无网格伽辽金方法的应用带来了一定的的困难,本文将再生核质点法中的混合变换法与无网格伽辽金方法相结合,通过对移动最小二乘近似进行局部修正,实现了无网格伽辽金方法中本质边界条件的精确施加。对权函数、影响半径、积分阶等对计算精度的影响进行了有益的探讨。数值计算结果表明了方法的可行性。  相似文献   

20.
IntroductionMeshless methods, as a special numerical method, originated from1970s. Since thediffuse element method was proposed by Nayroleset al.[1]in1992, the meshless methodshave received wide attentions in the mechanics area, and have shown obvious adv…  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号