Weakened weak form

Last updated

Weakened weak form (or W2 form) [1] 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.

Contents

Description

For simplicity we choose elasticity problems (2nd order PDE) for our discussion. [2] Our discussion is also most convenient in reference to the well-known weak and strong form. In a strong formulation for an approximate solution, we need to assume displacement functions that are 2nd order differentiable. In a weak formulation, we create linear and bilinear forms and then search for a particular function (an approximate solution) that satisfy the weak statement. The bilinear form uses gradient of the functions that has only 1st order differentiation. Therefore, the requirement on the continuity of assumed displacement functions is weaker than in the strong formulation. In a discrete form (such as the Finite element method, or FEM), a sufficient requirement for an assumed displacement function is piecewise continuous over the entire problems domain. This allows us to construct the function using elements (but making sure it is continuous a long all element interfaces), leading to the powerful FEM.

Now, in a weakened weak (W2) formulation, we further reduce the requirement. We form a bilinear form using only the assumed function (not even the gradient). This is done by using the so-called generalized gradient smoothing technique, [3] with which one can approximate the gradient of displacement functions for certain class of discontinuous functions, as long as they are in a proper G space. [4] Since we do not have to actually perform even the 1st differentiation to the assumed displacement functions, the requirement on the consistence of the functions are further reduced, and hence the weakened weak or W2 formulation.

History

The development of systematic theory of the weakened weak form started from the works on meshfree methods. [2] It is relatively new, but had very rapid development in the past few years.[ when? ]

Features of W2 formulations

  1. 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. This is very important for our long-term goal of development of fully automated computational methods.
  2. 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. This is important for producing so-called certified solutions.
  3. W2 models can be built free from volumetric locking, and possibly free from other types of locking phenomena.
  4. W2 models provide the freedom to assume separately the displacement gradient of the displacement functions, offering opportunities for ultra-accurate and super-convergent models. It may be possible to construct linear models with energy convergence rate of 2.
  5. W2 models are often found less sensitive to mesh distortion.
  6. W2 models are found effective for low order methods

Existing W2 models

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. 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] The S-FEM models are:

Applications

Some of the applications of W2 models are:

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

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.

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">Computational electromagnetics</span> Branch of physics

Computational electromagnetics (CEM), computational electrodynamics or electromagnetic modeling is the process of modeling the interaction of electromagnetic fields with physical objects and the environment using computers.

<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.

<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.

<span class="mw-page-title-main">Discontinuity layout optimization</span>

Discontinuity layout optimization (DLO) is an engineering analysis procedure which can be used to directly establish the amount of load that can be carried by a solid or structure prior to collapse. Using DLO the layout of failure planes, or 'discontinuities', in a collapsing solid or structure are identified using mathematical optimization methods. It is assumed that failure occurs in a ductile or 'plastic' manner.

Smoothed finite element methods (S-FEM) are a particular class of numerical simulation algorithms for the simulation of physical phenomena. It was developed by combining meshfree methods 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.

In numerical mathematics, the regularized meshless method (RMM), also known as the singular meshless method or desingularized meshless method, is a meshless boundary collocation method designed to solve certain partial differential equations whose fundamental solution is explicitly known. The RMM is a strong-form collocation method with merits being meshless, integration-free, easy-to-implement, and high stability. Until now this method has been successfully applied to some typical problems, such as potential, acoustics, water wave, and inverse problems of bounded and unbounded domains.

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.

Fluid motion is governed by the Navier–Stokes equations, a set of coupled and nonlinear partial differential equations derived from the basic laws of conservation of mass, momentum and energy. The unknowns are usually the flow velocity, the pressure and density and temperature. The analytical solution of this equation is impossible hence scientists resort to laboratory experiments in such situations. The answers delivered are, however, usually qualitatively different since dynamical and geometric similitude are difficult to enforce simultaneously between the lab experiment and the prototype. Furthermore, the design and construction of these experiments can be difficult, particularly for stratified rotating flows. Computational fluid dynamics (CFD) is an additional tool in the arsenal of scientists. In its early days CFD was often controversial, as it involved additional approximation to the governing equations and raised additional (legitimate) issues. Nowadays CFD is an established discipline alongside theoretical and experimental methods. This position is in large part due to the exponential growth of computer power which has allowed us to tackle ever larger and more complex problems.

MoFEM is an open source finite element analysis code developed and maintained at the University of Glasgow. MoFEM is tailored for the solution of multi-physics problems with arbitrary levels of approximation, different levels of mesh refinement and optimised for high-performance computing. MoFEM is the blend of the Boost MultiIndex containers, MOAB and PETSc. MoFEM is developed in C++ and it is open-source software under the GNU Lesser General Public License (GPL).

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.

Simulation governance is a managerial function concerned with assurance of reliability of information generated by numerical simulation. The term was introduced in 2011 and specific technical requirements were addressed from the perspective of mechanical design in 2012. Its strategic importance was addressed in 2015. At the 2017 NAFEMS World Congress in Stockholm simulation governance was identified as the first of eight “big issues” in numerical simulation.

<span class="mw-page-title-main">Magd Abdel Wahab</span> Belgian academic

Magd Abdel Wahab is a Belgian academic, researcher, author and Imam of Islam. He is full professor and chair of applied mechanics at Ghent University, Belgium, where he is also the Head of Finite Element Modelling Research Group of Laboratory Soete.

References

  1. ↑ 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
  2. 1 2 Liu, G.R. 2nd edn: 2009 Mesh Free Methods, CRC Press. 978-1-4200-8209-9
  3. ↑ Liu GR, "A Generalized Gradient Smoothing Technique and the Smoothed Bilinear Form for Galerkin Formulation of a Wide Class of Computational Methods", International Journal of Computational Methods , Vol. 5, Issue: 2, 199–236, 2008
  4. ↑ Liu GR, "On G Space Theory", International Journal of Computational Methods, Vol. 6 Issue: 2, 257–289, 2009
  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". International Journal for Numerical Methods in Engineering. 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. ↑ 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
  20. ↑ 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
  21. ↑ 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
  22. ↑ 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
  23. ↑ 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
  24. ↑ 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
  25. ↑ 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
  26. ↑ Liu GR, Nourbakhshnia N, Chen L, et al. "A Novel General Formulation for Singular Stress Field Using the Es-Fem Method for the Analysis of Mixed-Mode Cracks", International Journal of Computational Methods, Vol. 7, Issue: 1, 191–214, 2010
  27. ↑ 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
  28. ↑ 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
  29. ↑ 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
  30. ↑ 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
  31. ↑ 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
  32. ↑ 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
  33. ↑ 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
  34. ↑ 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
  35. ↑ 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
  36. ↑ 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
  37. ↑ Li E, Liu GR, Tan V, et al. "An efficient algorithm for phase change problem in tumor treatment using alpha FEM", International Journal of Thermal Sciences, Vol. 49, Issue: 10, 1954–1967, 2010
  38. ↑ 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
  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
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.