|
=====================================
最近无网格方法发展很快,大有取代有限元方法的势头。可是我们的论坛上关于这个话题的讨论还很少。这种方法对传统的有限元方法进行了改造,特别是在插值函数构造上。
无网格方法一改有限元方法采用单元进行计算区域离散,而是采用离散点进行区域离散——不用划分网格!!这是无网格方法的最大卖点。
我这里先抛砖引玉,如果大家有兴趣,我们可以继续讨论。
=====================================
80年代初,Liszka和Orkisz等提出广义有限差分法(Generalized Finite Difference Method, 简写为GFDM,Lucy以及R.A.Gingold和J.J.Monaghan等各自独立地提出Smoothed particle Hydrodynamics(SPH),成为最早的无网格方法。此后,相继出现了The diffuse element method(Nayroles,1992; Krongauz, 1997), The element free Galerkin method (Belytschko, 1994;Lu,1994;Belytschko ,1995;Beissel,1996), Wavelet Galerkin Method(Mattos, 2003), Reproducing kernel particle method(Liu,1996; Chen, 1996; Li, 1996; Hulbert, 1996; Aluru, 1999), Hp-clouds method(Duarte, 1995), Partition of unity method(Melenk,1996; 刘欣,2001), finite point method(Onate,1996),natural element method(Braun,1995;Sukumar, 1998), free mesh method(Yagawa,1996)等,人们对无网格方法的认识逐渐深入。
=====================================
无 网 格 法
近几十年来,有限元法已成为计算力学中解决工程问题的主要数值手段,然而随着其应用范围的扩展,其固有的一些缺陷也日益突出。在金属冲压成形、高速碰撞、流固耦合等涉及特大变形的领域中,基于拉格朗日法的有限元网格可能产生严重的扭曲,甚至使得单元的雅可比行列式为负值,不仅在计算中需要网格重构,而且严重地影响解的精度;对高速冲击等动态问题,显式时间积分的步长取决于有限元网格的最小尺寸,因而网格的扭曲将使得时间积分步长过小,大幅度地增加了计算工作量;对裂纹的动态扩展问题,由于裂纹的扩展方向不能事先确定,因而在计算过程中需要不断地重新划分网格以模拟裂纹的动态扩展过程。由于有限元近似基于网格,因此必然难于处理与原始网格线不一致的不连续性和大变形。网格重构不仅计算费用昂贵,而且会损害计算精度。鉴于这种缺陷,近几年来国际上许多著名的计算力学学者,如 T. Belytschko, O.C. Zienkiewicz, S.N. Atluri, J.T. Oden, W.K. Liu 等都对无网格方法表现出了极大的兴趣,并进行了大量的研究工作。无网格方法采用基于点的近似,可以彻底或部分地消除网格,不需要网格的初始划分和重构,不仅可以保证计算的精度,而且可以大大减小计算的难度。然而,由于目前的无网格近似一般没有解析表达式,且大都基于伽辽金原理,因此计算量很大,要超出传统的有限元法;另外,无网格近似大都是拟合,因此对于位移边界的处理比较困难,多采用拉格朗日乘子法处理。
目前已提出了十余种无网格法,其主要区别在于离散微分方程的方法(如伽辽金法、配点法、最小二乘法、彼得洛夫-伽辽金法等)和建立近似函数的方法(移动最小二乘近似、核近似、重构核质点近似、单位分解法、hp云团法、径向基函数法、点插值法等)。
=====================================
有关书籍
《无网格法》(张雄,刘岩著,清华大学出版社/Springer出版社,2004年8月)以紧支试函数加权残量法为主线,系统地论述了目前现有的各种无网格方法的基本原理以及它们之间的区别与联系,在此基础上建立了一些新型有效的无网格方法,如最小二乘配点无网格法、加权最小二乘无网格法、伽辽金最小二乘无网格法、伽辽金配点无网格法等,其中许多内容都是作者课题组在国家自然科学基金资助下取得的成果。书中给出了作者用C++研制的面向对象的无网格法程序OMLL,另外各章都给出了相应的MATLAB 程序,以帮助读者理解各种无网格方法的程序实现过程。
=====================================
clili,你好。很高兴能认识你。
我觉得在有限元上我们国家已经落后太多了,商业数值计算基本上被美国的有限元公司如MSC、ANSYS等垄断。要想在这方面与他们竞争是很难的。
但是,现在无网格的商业程序还没有出现,大家都处在同一个起跑线上,这是我们在数值计算领域占有一席之地的好机会。因此,希望各位同行同心协力,希望能在中国率先推出性能卓越的三维无网格商业软件。
=====================================
我正在做无网格方法在断裂力学中的应用,其中要用到不连续性的判断,这个东东虽然原理很简单,可是涉及到自动布点和计算权函数时dI的求解。不知道有哪位大虾能提供点思路不~~
=====================================
无网格方法在断裂力学方面的应用这个题目国内已经有人做过了,你可以查一下相关的博士论文。
我不知道你说得dI是什么意思,我想模拟不连续性就是几种方法:可视性准则、透射准则、衍射准则,这在文献中也是可以找到的。
=====================================
感谢yongjq兄的指点。dI指高斯型权函数中高斯点跟节点之间的距离,即高斯型权函数的自变量。我现在就是用衍射法则作的,可是计算效果不好,在裂纹尖端部分计算的位移振荡很明显。我怀疑是否在用衍射法则时有一些其他的处理技术在里边?我手中的资料都是只讲理论,没有讲其具体实现步骤或者在实现过程出现的一些情况。
你是如何看这个问题的?
=====================================
Star_R_A,是的,我看过的资料里都没有提到采用衍射法则后位移场会出现什么情况,也没有很细致的描述。这些问题我看只能慢慢摸索了。
看来你是要用滑动最小二乘法?希望能看到你的位移到底是怎么震荡的。我想大家一起想办法,应该可以解决这些问题。
=====================================
我就是用伽辽金方法计算的,其核心是最小二乘法。用这个方法计算没有不连续线(我计算的是裂纹)时,在裂纹尖端很靠近的几个点上位移出现了振荡而其它绝大多数点上的位移还是非常好的。我不知道别人在计算时是所有点都很好还是也像我一样在个别点上不好。
很郁闷,如果可以解决这个问题我就可以完成论文可以早点回家了~~
=====================================
无网格法》(张雄,刘岩著,清华大学出版社/Springer出版社,2004年8月)这本书谁有?可以传上来........感谢
=====================================
台湾中原大学硕士论文
元素释放法于弹性动力之研究
注:元素释放法为无网格法的一种 |
|