An extended Galerkin weak form and a point interpolation method with continuous strain field and superconvergence using triangular mesh

被引:28
作者
Liu, G. R. [2 ,3 ]
Xu, X. [1 ,3 ]
Zhang, G. Y. [3 ]
Gu, Y. T. [4 ]
机构
[1] Jilin Univ, Coll Math, Changchun 130012, Peoples R China
[2] Natl Univ Singapore, Ctr Adv Computat Engn Sci, Dept Mech Engn, Singapore 117576, Singapore
[3] SMA, Singapore 117576, Singapore
[4] Queensland Univ Technol, Sch Engn Syst, Brisbane, Qld 4001, Australia
基金
中国国家自然科学基金;
关键词
Finite element method; Meshfree methods; Point interpolation method; Solution bound; Superconvergence; x-Galerkin weak form; FINITE-ELEMENT-METHOD; SOLID MECHANICS PROBLEMS; METHOD LC-PIM; PARTICLE METHODS; METHOD SFEM; COMPUTATIONAL MECHANICS; ELASTICITY PROBLEMS; STABILITY ANALYSIS; CRACK INITIATION; PROPAGATION;
D O I
10.1007/s00466-008-0336-5
中图分类号
O1 [数学];
学科分类号
0701 ; 070101 ;
摘要
A point interpolation method (PIM) with continuous strain field (PIM-CS) is developed for mechanics problems using triangular background mesh, in which PIM shape functions are used to construct both displacement and strain fields. The strain field constructed is continuous in the entire problem domain, which is achieved by simple linear interpolations using locally smoothed strains around the nodes and points required for the interpolation. A general parameterized functional with a real adjustable parameter a are then proposed for establishing PIM-CS models of special property. We prove theoretically that the PIM-CS has an excellent bound property: strain energy obtained using PIM-CS lies in between those of the compatible FEM and NS-PIM models of the same mesh. Techniques and procedures are then presented to compute the upper and lower bound solutions using the PIM-CS. It is discovered that an extended Galerkin (x-Galerkin) model, as special case resulted from the extended parameterized functional with alpha = 1, is outstanding in terms of both performance and efficiency. Intensive numerical studies show that upper and lower bound solutions can always be obtained, there exist a values at which the solutions of PIM-CS are of superconvergence, and the x-Galerkin model is capable of producing superconvergent solutions of ultra accuracy that is about 10 times that of the FEM using the same mesh.
引用
收藏
页码:651 / 673
页数:23
相关论文
共 65 条
[1]  
Anderson - T.L., 1991, FRACTURE MECH FUNDAM
[2]  
[Anonymous], 2002, Cambridge Texts in Applied Mathematics, DOI [10.1017/CBO9780511791253, DOI 10.1017/CBO9780511791253]
[3]  
[Anonymous], 2003, Smoothed particle hydrodynamics: a meshfree particle method, DOI DOI 10.1007/S00466-004-0573-1
[4]   A new meshless local Petrov-Galerkin (MLPG) approach in computational mechanics [J].
Atluri, SN ;
Zhu, T .
COMPUTATIONAL MECHANICS, 1998, 22 (02) :117-127
[5]   ELEMENT-FREE GALERKIN METHODS [J].
BELYTSCHKO, T ;
LU, YY ;
GU, L .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 1994, 37 (02) :229-256
[6]   Three-dimensional crack initiation, propagation, branching and junction in non-linear materials by an extended meshfree method without asymptotic enrichment [J].
Bordas, Stephane ;
Rabczuk, Timon ;
Zi, Goangseup .
ENGINEERING FRACTURE MECHANICS, 2008, 75 (05) :943-960
[7]  
Braess Dietrich., 2007, Finite Elemente
[8]  
Brenner S.C., 2007, The mathematical theory of finite element methods, V3rd
[9]  
Chen JS, 2001, INT J NUMER METH ENG, V50, P435, DOI 10.1002/1097-0207(20010120)50:2<435::AID-NME32>3.0.CO
[10]  
2-A