Smoothed finite element method

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 ("Smoothed Finite Element Methods (S-FEM): An Overview and Recent Developments", <http://link.springer.com/article/10.1007/s11831-016-9202-3/fulltext.html>).

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 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). The S-PIM can be node-based (known as NS-PIM or LC-PIM), edge-based (ES-PIM), and cell-based (CS-PIM). The NS-PIM was developed using the so-called SCNI technique. It was then discovered that NS-PIM is capable of producing upper bound solution and volumetric locking free. 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.

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.

List of S-FEM models

Applications

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

1. Mechanics for solid structures and piezoelectrics;
2. Fracture mechanics and crack propagation;
3. Nonlinear and contact problems;
4. Stochastic analysis;
5. Heat transfer;
6. Structural acoustics;
7. Adaptive analysis;
8. Limited analysis;
9. Crystal plasticity modeling.