Smoothed finite element method

Last updated

Smoothed finite element methods (S-FEM) [1] are a particular class of numerical simulation algorithms for the simulation of physical phenomena. It was developed by combining meshfree methods [2] with the finite element method. S-FEM are applicable to solid mechanics as well as fluid dynamics problems, although so far they have mainly been applied to the former.

Contents

Description

The essential idea in the S-FEM is to use a finite element mesh (in particular triangular mesh) to construct numerical models of good performance. This is achieved by modifying the compatible strain field, or construct a strain field using only the displacements, hoping a Galerkin model using the modified/constructed strain field can deliver some good properties. Such a modification/construction can be performed within elements but more often beyond the elements (meshfree concepts): bring in the information from the neighboring elements. Naturally, the strain field has to satisfy certain conditions, and the standard Galerkin weak form needs to be modified accordingly to ensure the stability and convergence. A comprehensive review of S-FEM covering both methodology and applications can be found in [3] ("Smoothed Finite Element Methods (S-FEM): An Overview and Recent Developments").

History

The development of S-FEM started from the works on meshfree methods, where the so-called weakened weak (W2) formulation based on the G space theory [4] were developed. The W2 formulation offers possibilities for formulate various (uniformly) "soft" models that works well with triangular meshes. Because triangular mesh can be generated automatically, it becomes much easier in re-meshing and hence automation in modeling and simulation. In addition, W2 models can be made soft enough (in uniform fashion) to produce upper bound solutions (for force-driving problems). Together with stiff models (such as the fully compatible FEM models), one can conveniently bound the solution from both sides. This allows easy error estimation for generally complicated problems, as long as a triangular mesh can be generated. Typical W2 models are the Smoothed Point Interpolation Methods (or S-PIM). [5] The S-PIM can be node-based (known as NS-PIM or LC-PIM), [6] edge-based (ES-PIM), [7] and cell-based (CS-PIM). [8] The NS-PIM was developed using the so-called SCNI technique. [9] It was then discovered that NS-PIM is capable of producing upper bound solution and volumetric locking free. [10] The ES-PIM is found superior in accuracy, and CS-PIM behaves in between the NS-PIM and ES-PIM. Moreover, W2 formulations allow the use of polynomial and radial basis functions in the creation of shape functions (it accommodates the discontinuous displacement functions, as long as it is in G1 space), which opens further rooms for future developments.[ citation needed ]

The S-FEM is largely the linear version of S-PIM, but with most of the properties of the S-PIM and much simpler. It has also variations of NS-FEM, ES-FEM and CS-FEM. The major property of S-PIM can be found also in S-FEM. [11]

List of S-FEM models

Applications

S-FEM has been applied to solve the following physical problems:

  1. Mechanics for solid structures and piezoelectrics; [24] [25]
  2. Fracture mechanics and crack propagation; [26] [27] [28] [29]
  3. Nonlinear and contact problems; [30] [31]
  4. Stochastic analysis; [32]
  5. Heat transfer; [33] [34]
  6. Structural acoustics; [35] [36] [37]
  7. Adaptive analysis; [38] [18]
  8. Limited analysis; [39]
  9. Crystal plasticity modeling. [40]

Basic Formulation of S-FEM

The fundamental problem addressed by SFEM is typically the solution of Poisson's equation with Dirichlet boundary conditions, given as follows:

Δu+f=0 in Ω, u=g on ΓD

where Ω is the domain and Γ is its boundary, consisting of ΓD=Γ. Here, u: Ω→R is the trial solution, f: Ω→R is a given function, and g represents Dirichlet boundary conditions.

S-FEM involves discretizing the domain Ω using finite element meshes, which can be global or local. The global mesh represents the entire domain, while the local mesh is used to discretize regions requiring high resolution within the global domain. The local domain is assumed to be included in the global domain (ΩL​⊆ΩG).

Weak Formulation

The weak form of the problem is derived by multiplying the equation by suitable test functions and integrating over the domain. In SFEM, the weak form is expressed as follows: Given f and g, find u∈U such that for all w∈V,

aΩ​(w,u)=LΩ​(w)

where aΩ is a bilinear form, and LΩ is a linear functional.

S-FEM Formulation

In S-FEM, the trial solution u and test functions w are defined separately for the global (ΩG) and local (ΩL) domains. The trial solution spaces UG, UL and test function spaces VG, VL are defined accordingly. The weak form in the S-FEM formulation becomes:

aΩ′​(w,u)=LΩ′​(w)

where aΩ′​(⋅,⋅) and LΩ′​(⋅) are modified bilinear forms and linear functionals, respectively, to accommodate the S-FEM approach.

Challenges

One of the primary challenges of S-FEM is the difficulty in exact integration of the submatrices representing the relationship between global and local meshes (KGL and KLG). Additionally, the matrix K can become singular, posing numerical challenges in solving the resulting linear algebraic equations.

These challenges and potential solutions are discussed in detail in the literature, aiming to improve the efficiency and accuracy of S-FEM for various applications. [41]

B-spline S-FEM (BFSEM)

S-FEM can reasonably model an analytical domain by superimposing meshes with different spatial resolutions, it has intrinsic advantages of local high accuracy, low computation time, and simple meshing procedure. However, it has disadvantages such as accuracy of numerical integration and matrix singularity. Although several additional techniques have been proposed to mitigate these limitations, they are computationally expensive or ad-hoc, and detract from the method’s strengths. These issues can be address by incorporating cubic B-spline functions with C squared continuity across element boundaries as the global basis function. To avoid matrix singularity, applying different basis functions to different meshes. In a recent study Lagrange basis functions were used as local basis functions. With this method the numerical integration can be calculated with sufficient accuracy without any additional techniques used in conventional S-FEM. Furthermore, the proposed method avoids matrix singularity and is superior to conventional methods in terms of convergence for solving linear equations. Therefore, the proposed method has the potential to reduce computation time while maintaining a comparable accuracy to conventional S-FEM. [41]

See also

Related Research Articles

<span class="mw-page-title-main">Computational fluid dynamics</span> Analysis and solving of problems that involve fluid flows

Computational fluid dynamics (CFD) is a branch of fluid mechanics that uses numerical analysis and data structures to analyze and solve problems that involve fluid flows. Computers are used to perform the calculations required to simulate the free-stream flow of the fluid, and the interaction of the fluid with surfaces defined by boundary conditions. With high-speed supercomputers, better solutions can be achieved, and are often required to solve the largest and most complex problems. Ongoing research yields software that improves the accuracy and speed of complex simulation scenarios such as transonic or turbulent flows. Initial validation of such software is typically performed using experimental apparatus such as wind tunnels. In addition, previously performed analytical or empirical analysis of a particular problem can be used for comparison. A final validation is often performed using full-scale testing, such as flight tests.

Topology optimization is a mathematical method that optimizes material layout within a given design space, for a given set of loads, boundary conditions and constraints with the goal of maximizing the performance of the system. Topology optimization is different from shape optimization and sizing optimization in the sense that the design can attain any shape within the design space, instead of dealing with predefined configurations.

Numerical methods for partial differential equations is the branch of numerical analysis that studies the numerical solution of partial differential equations (PDEs).

The boundary element method (BEM) is a numerical computational method of solving linear partial differential equations which have been formulated as integral equations, including fluid mechanics, acoustics, electromagnetics, fracture mechanics, and contact mechanics.

<span class="mw-page-title-main">Mesh generation</span> Subdivision of space into cells

Mesh generation is the practice of creating a mesh, a subdivision of a continuous geometric space into discrete geometric and topological cells. Often these cells form a simplicial complex. Usually the cells partition the geometric input domain. Mesh cells are used as discrete local approximations of the larger domain. Meshes are created by computer algorithms, often with human guidance through a GUI, depending on the complexity of the domain and the type of mesh desired. A typical goal is to create a mesh that accurately captures the input domain geometry, with high-quality (well-shaped) cells, and without so many cells as to make subsequent calculations intractable. The mesh should also be fine in areas that are important for the subsequent calculations.

<span class="mw-page-title-main">Fluid–structure interaction</span>

Fluid–structure interaction (FSI) is the interaction of some movable or deformable structure with an internal or surrounding fluid flow. Fluid–structure interactions can be stable or oscillatory. In oscillatory interactions, the strain induced in the solid structure causes it to move such that the source of strain is reduced, and the structure returns to its former state only for the process to repeat.

<span class="mw-page-title-main">Meshfree methods</span> Methods in numerical analysis not requiring knowledge of neighboring points

In the field of numerical analysis, meshfree methods are those that do not require connection between nodes of the simulation domain, i.e. a mesh, but are rather based on interaction of each node with all its neighbors. As a consequence, original extensive properties such as mass or kinetic energy are no longer assigned to mesh elements but rather to the single nodes. Meshfree methods enable the simulation of some otherwise difficult types of problems, at the cost of extra computing time and programming effort. The absence of a mesh allows Lagrangian simulations, in which the nodes can move according to the velocity field.

The finite-difference frequency-domain (FDFD) method is a numerical solution method for problems usually in electromagnetism and sometimes in acoustics, based on finite-difference approximations of the derivative operators in the differential equation being solved.

In the numerical solution of partial differential equations, a topic in mathematics, the spectral element method (SEM) is a formulation of the finite element method (FEM) that uses high-degree piecewise polynomials as basis functions. The spectral element method was introduced in a 1984 paper by A. T. Patera. Although Patera is credited with development of the method, his work was a rediscovery of an existing method

<span class="mw-page-title-main">Finite element method</span> Numerical method for solving physical or engineering problems

The finite element method (FEM) is a popular method for numerically solving differential equations arising in engineering and mathematical modeling. Typical problem areas of interest include the traditional fields of structural analysis, heat transfer, fluid flow, mass transport, and electromagnetic potential.

<span class="mw-page-title-main">Hermes Project</span>

Hermes2D is a C++/Python library of algorithms for rapid development of adaptive hp-FEM solvers. hp-FEM is a modern version of the finite element method (FEM) that is capable of extremely fast, exponential convergence.

hp-FEM is a generalization of the finite element method (FEM) for solving partial differential equations numerically based on piecewise-polynomial approximations. hp-FEM originates from the discovery by Barna A. Szabó and Ivo Babuška that the finite element method converges exponentially fast when the mesh is refined using a suitable combination of h-refinements and p-refinements .The exponential convergence of hp-FEM has been observed by numerous independent researchers.

OOFEM is a free and open-source multi-physics finite element code with object oriented architecture. The aim of this project is to provide efficient and robust tool for FEM computations as well as to offer highly modular and extensible environment for development.

Weakened weak form is used in the formulation of general numerical methods based on meshfree methods and/or finite element method settings. These numerical methods are applicable to solid mechanics as well as fluid dynamics problems.

In numerical mathematics, the boundary knot method (BKM) is proposed as an alternative boundary-type meshfree distance function collocation scheme.

The Kansa method is a computer method used to solve partial differential equations. Its main advantage is it is very easy to understand and program on a computer. It is much less complicated than the finite element method. Another advantage is it works well on multi variable problems. The finite element method is complicated when working with more than 3 space variables and time.

p-FEM or the p-version of the finite element method is a numerical method for solving partial differential equations. It is a discretization strategy in which the finite element mesh is fixed and the polynomial degrees of elements are increased such that the lowest polynomial degree, denoted by , approaches infinity. This is in contrast with the "h-version" or "h-FEM", a widely used discretization strategy, in which the polynomial degrees of elements are fixed and the mesh is refined such that the diameter of the largest element, denoted by approaches zero.

The finite point method (FPM) is a meshfree method for solving partial differential equations (PDEs) on scattered distributions of points. The FPM was proposed in the mid-nineties in, and with the purpose to facilitate the solution of problems involving complex geometries, free surfaces, moving boundaries and adaptive refinement. Since then, the FPM has evolved considerably, showing satisfactory accuracy and capabilities to deal with different fluid and solid mechanics problems.

<span class="mw-page-title-main">FEATool Multiphysics</span>

FEATool Multiphysics is a physics, finite element analysis (FEA), and partial differential equation (PDE) simulation toolbox. FEATool Multiphysics features the ability to model fully coupled heat transfer, fluid dynamics, chemical engineering, structural mechanics, fluid-structure interaction (FSI), electromagnetics, as well as user-defined and custom PDE problems in 1D, 2D (axisymmetry), or 3D, all within a graphical user interface (GUI) or optionally as script files. FEATool has been employed and used in academic research, teaching, and industrial engineering simulation contexts.

<span class="mw-page-title-main">Michel Bercovier</span> Professor of Scientific Computing

Michel Bercovier is a French-Israeli Professor (Emeritus) of Scientific Computing and Computer Aided Design (CAD) in The Rachel and Selim Benin School of Computer Science and Engineering at the Hebrew University of Jerusalem. Bercovier is also the head of the School of Computer Science at the Hadassah Academic College, Jerusalem.

References

  1. ↑ Liu, G.R., 2010 Smoothed Finite Element Methods, CRC Press, ISBN   978-1-4398-2027-8.
  2. ↑ Liu, G.R. 2nd edn: 2009 Mesh Free Methods, CRC Press. 978-1-4200-8209-9
  3. ↑ W. Zeng, G.R. Liu. Smoothed finite element methods (S-FEM): An overview and recent developments. Archives of Computational Methods in Engineering, 2016, doi: 10.1007/s11831-016-9202-3
  4. ↑ G.R. Liu. A G space theory and a weakened weak (W2) form for a unified formulation of compatible and incompatible methods: Part I theory and Part II applications to solid mechanics problems. International Journal for Numerical Methods in Engineering, 81: 1093-1126, 2010
  5. ↑ Liu, G.R. 2nd edn: 2009 Mesh Free Methods, CRC Press. 978-1-4200-8209-9
  6. ↑ Liu GR, Zhang GY, Dai KY, Wang YY, Zhong ZH, Li GY and Han X, A linearly conforming point interpolation method (LC-PIM) for 2D solid mechanics problems, International Journal of Computational Methods , 2(4): 645-665, 2005.
  7. ↑ G.R. Liu, G.R. Zhang. Edge-based Smoothed Point Interpolation Methods. International Journal of Computational Methods, 5(4): 621-646, 2008
  8. ↑ G.R. Liu, G.R. Zhang. A normed G space and weakened weak (W2) formulation of a cell-based Smoothed Point Interpolation Method. International Journal of Computational Methods, 6(1): 147-179, 2009
  9. ↑ Chen, J. S., Wu, C. T., Yoon, S. and You, Y. (2001). A stabilized conforming nodal integration for Galerkin mesh-free methods. Int. J. Numer. Meth. Eng. 50: 435–466.
  10. ↑ G. R. Liu and G. Y. Zhang. Upper bound solution to elasticity problems: A unique property of the linearly conforming point interpolation method (LC-PIM). International Journal for Numerical Methods in Engineering, 74: 1128-1161, 2008.
  11. ↑ Zhang ZQ, Liu GR, Upper and lower bounds for natural frequencies: A property of the smoothed finite element methods, International Journal for Numerical Methods in Engineering Vol. 84 Issue: 2, 149-178, 2010
  12. ↑ Liu GR, Nguyen-Thoi T, Nguyen-Xuan H, Lam KY (2009) A node-based smoothed finite element method (NS-FEM) for upper bound solutions to solid mechanics problems. Computers and Structures; 87: 14-26.
  13. ↑ Liu GR, Nguyen-Thoi T, Lam KY (2009) An edge-based smoothed finite element method (ES-FEM) for static, free and forced vibration analyses in solids. Journal of Sound and Vibration; 320: 1100-1130.
  14. ↑ Nguyen-Thoi T, Liu GR, Lam KY, GY Zhang (2009) A Face-based Smoothed Finite Element Method (FS-FEM) for 3D linear and nonlinear solid mechanics problems using 4-node tetrahedral elements. International Journal for Numerical Methods in Engineering; 78: 324-353
  15. ↑ Liu GR, Dai KY, Nguyen-Thoi T (2007) A smoothed finite element method for mechanics problems. Computational Mechanics; 39: 859-877
  16. ↑ Dai KY, Liu GR (2007) Free and forced vibration analysis using the smoothed finite element method (SFEM). Journal of Sound and Vibration; 301: 803-820.
  17. ↑ Dai KY, Liu GR, Nguyen-Thoi T (2007) An n-sided polygonal smoothed finite element method (nSFEM) for solid mechanics. Finite Elements in Analysis and Design; 43: 847－860.
  18. 1 2 Li Y, Liu GR, Zhang GY, An adaptive NS/ES-FEM approach for 2D contact problems using triangular elements, Finite Elements in Analysis and Design Vol.47 Issue: 3, 256-275, 2011
  19. ↑ Jiang C, Zhang ZQ, Liu GR, Han X, Zeng W, An edge-based/node-based selective smoothed finite element method using tetrahedrons for cardiovascular tissues, Engineering Analysis with Boundary Elements Vol.59, 62-77, 2015
  20. ↑ Liu GR, Nguyen-Thoi T, Lam KY (2009) A novel FEM by scaling the gradient of strains with factor α (αFEM). Computational Mechanics; 43: 369-391
  21. ↑ Liu GR, Nguyen-Xuan H, Nguyen-Thoi T, Xu X (2009) A novel weak form and a superconvergent alpha finite element method (SαFEM) for mechanics problems using triangular meshes. Journal of Computational Physics; 228: 4055-4087
  22. ↑ Zeng W, Liu GR, Li D, Dong XW (2016) A smoothing technique based beta finite element method (βFEM) for crystal plasticity modeling. Computers and Structures; 162: 48-67
  23. ↑ Zeng W, Liu GR, Jiang C, Nguyen-Thoi T, Jiang Y (2016) A generalized beta finite element method with coupled smoothing techniques for solid mechanics. Engineering Analysis with Boundary Elements; 73: 103-119
  24. ↑ Cui XY, Liu GR, Li GY, et al. A thin plate formulation without rotation DOFs based on the radial point interpolation method and triangular cells, International Journal for Numerical Methods in Engineering Vol.85 Issue: 8, 958-986, 2011
  25. ↑ Liu GR, Nguyen-Xuan H, Nguyen-Thoi T, A theoretical study on the smoothed FEM (S-FEM) models: Properties, accuracy and convergence rates, International Journal for Numerical Methods in Engineering Vol. 84 Issue: 10, 1222-1256, 2010
  26. ↑ Liu GR, Nourbakhshnia N, Zhang YW, A novel singular ES-FEM method for simulating singular stress fields near the crack tips for linear fracture problems, Engineering Fracture Mechanics Vol.78 Issue: 6 Pages: 863-876, 2011
  27. ↑ Liu GR, Chen L, Nguyen-Thoi T, et al. A novel singular node-based smoothed finite element method (NS-FEM) for upper bound solutions of fracture problems, International Journal for Numerical Methods in Engineering Vol.83 Issue: 11, 1466-1497, 2010
  28. ↑ Zeng W, Liu GR, Kitamura Y, Nguyen-Xuan H. "A three-dimensional ES-FEM for fracture mechanics problems in elastic solids", Engineering Fracture Mechanics Vol. 114, 127-150, 2013
  29. ↑ Zeng W, Liu GR, Jiang C, Dong XW, Chen HD, Bao Y, Jiang Y. "An effective fracture analysis method based on the virtual crack closure-integral technique implemented in CS-FEM", Applied Mathematical Modelling Vol. 40, Issue: 5-6, 3783-3800, 2016
  30. ↑ Zhang ZQ, Liu GR, An edge-based smoothed finite element method (ES-FEM) using 3-node triangular elements for 3D non-linear analysis of spatial membrane structures, International Journal for Numerical Methods in Engineering, Vol. 86 Issue: 2 135-154, 2011
  31. ↑ Jiang C, Liu GR, Han X, Zhang ZQ, Zeng W, A smoothed finite element method for analysis of anisotropic large deformation of passive rabbit ventricles in diastole, International Journal for Numerical Methods in Biomedical Engineering , Vol. 31 Issue: 1,1-25, 2015
  32. ↑ Liu GR, Zeng W, Nguyen-Xuan H. Generalized stochastic cell-based smoothed finite element method (GS_CS-FEM) for solid mechanics, Finite Elements in Analysis and Design Vol.63, 51-61, 2013
  33. ↑ Zhang ZB, Wu SC, Liu GR, et al. Nonlinear Transient Heat Transfer Problems using the Meshfree ES-PIM, International Journal of Nonlinear Sciences and Numerical Simulation Vol.11 Issue: 12, 1077-1091, 2010
  34. ↑ Wu SC, Liu GR, Cui XY, et al. An edge-based smoothed point interpolation method (ES-PIM) for heat transfer analysis of rapid manufacturing system, International Journal of Heat and Mass Transfer Vol.53 Issue: 9-10, 1938-1950, 2010
  35. ↑ He ZC, Cheng AG, Zhang GY, et al. Dispersion error reduction for acoustic problems using the edge-based smoothed finite element method (ES-FEM), International Journal for Numerical Methods in Engineering Vol. 86 Issue: 11 Pages: 1322-1338, 2011
  36. ↑ He ZC, Liu GR, Zhong ZH, et al. A coupled ES-FEM/BEM method for fluid-structure interaction problems, Engineering Analysis with Boundary Elements Vol. 35 Issue: 1, 140-147, 2011
  37. ↑ Zhang ZQ, Liu GR, Upper and lower bounds for natural frequencies: A property of the smoothed finite element methods, International Journal for Numerical Methods in Engineering Vol.84 Issue: 2,149-178, 2010
  38. ↑ Nguyen-Thoi T, Liu GR, Nguyen-Xuan H, et al. Adaptive analysis using the node-based smoothed finite element method (NS-FEM), International Journal for Numerical Methods in Biomedical Engineering Vol. 27 Issue: 2, 198-218, 2011
  39. ↑ Tran TN, Liu GR, Nguyen-Xuan H, et al. An edge-based smoothed finite element method for primal-dual shakedown analysis of structures, International Journal for Numerical Methods in Engineering Vol.82 Issue: 7, 917-938, 2010
  40. ↑ Zeng W, Larsen JM, Liu GR. Smoothing technique based crystal plasticity finite element modeling of crystalline materials, International Journal of Plasticity Vol.65, 250-268, 2015
  41. 1 2 Magome, Nozomi; Morita, Naoki; Kaneko, Shigeki; Mitsume, Naoto (January 2024). "Higher-continuity s-version of finite element method with B-spline functions". Journal of Computational Physics. 497: 112593. doi:10.1016/j.jcp.2023.112593. ISSN   0021-9991.
This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.