finite element method

 

  • “[15] Generalized finite element method[edit] The generalized finite element method (GFEM) uses local spaces consisting of functions, not necessarily polynomials, that reflect
    the available information on the unknown solution and thus ensure good local approximation.

  • A discretization strategy is understood to mean a clearly defined set of procedures that cover (a) the creation of finite element meshes, (b) the definition of basis function
    on reference elements (also called shape functions), and (c) the mapping of reference elements onto the elements of the mesh.

  • The global system of equations has known solution techniques and can be calculated from the initial values of the original problem to obtain a numerical answer.

  • Although the approaches used by these pioneers are different, they share one essential characteristic: mesh discretization of a continuous domain into a set of discrete sub-domains,
    usually called elements.

  • Extended finite element methods enrich the approximation space to naturally reproduce the challenging feature associated with the problem of interest: the discontinuity, singularity,
    boundary layer, etc.

  • So we now have to solve a linear system in the unknown where most of the entries of the matrix , which we need to invert, are zero.

  • The finite element method formulation of a boundary value problem finally results in a system of algebraic equations.

  • Typical work out of the method involves: 1. dividing the domain of the problem into a collection of subdomains, with each subdomain represented by a set of element equations
    to the original problem 2. systematically recombining all sets of element equations into a global system of equations for the final calculation.

  • In applying FEA, the complex problem is usually a physical system with the underlying physics such as the Euler–Bernoulli beam equation, the heat equation, or the Navier-Stokes
    equations expressed in either PDE or integral equations, while the divided small elements of the complex problem represent different areas in the physical system.

  • As we refine the triangulation, the space of piecewise linear functions must also change with .

  • A piecewise linear function in two dimensions For problem P2[edit] We need to be a set of functions of .

  • Basic concepts The subdivision of a whole domain into simpler parts has several advantages:[2] • Accurate representation of complex geometry • Inclusion of dissimilar material
    properties • Easy representation of the total solution • Capture of local effects.

  • For instance, for a fourth-order problem such as , one may use piecewise quadratic basis functions that are .

  • (b) The sparse matrix L of the discretized linear system(c) The computed solution, The primary advantage of this choice of basis is that the inner products and will be zero
    for almost all .

  • In contrast, ordinary differential equation sets that occur in the transient problems are solved by numerical integration using standard techniques such as Euler’s method
    or the Runge-Kutta method.

  • In the one-dimensional case, for each control point we will choose the piecewise linear function in whose value is at and zero at every , i.e., for ; this basis is a shifted
    and scaled tent function.

  • However, this method of solving the boundary value problem (BVP) works only when there is one spatial dimension.

  • FEA may be used for analyzing problems over complicated domains (like cars and oil pipelines) when the domain changes (as during a solid-state reaction with a moving boundary),
    when the desired precision varies over the entire domain, or when the solution lacks smoothness.

  • For vector partial differential equations, the basis functions may take values in .

  • When the errors of approximation are larger than what is considered acceptable, then the discretization has to be changed either by an automated adaptive process or by the
    action of the analyst.

  • The weak form of P2[edit] If we integrate by parts using a form of Green’s identities, we see that if solves P2, then we may define for any by where denotes the gradient and
    denotes the dot product in the two-dimensional plane.

  • [3] History While it is difficult to quote the date of the invention of the finite element method, the method originated from the need to solve complex elasticity and structural
    analysis problems in civil and aeronautical engineering.

  • A clear, detailed, and practical presentation of this approach can be found in The Finite Element Method for Engineers.

  • CutFEM[edit] The Cut Finite Element Approach was developed in 2014.

  • [14] The approach is “to make the discretization as independent as possible of the geometric description and minimize the complexity of mesh generation, while retaining the
    accuracy and robustness of a standard finite element method.

  • More advanced implementations (adaptive finite element methods) utilize a method to assess the quality of the results (based on error estimation theory) and modify the mesh
    during the solution aiming to achieve an approximate solution within some bounds from the exact solution of the continuum problem.

  • The FEM is a general numerical method for solving partial differential equations in two or three space variables (i.e., some boundary value problems).

  • One hopes that as the underlying triangular mesh becomes finer and finer, the solution of the discrete problem (3) will, in some sense, converge to the solution of the original
    boundary value problem P2.

  • P1 is a one-dimensional problem where is given, is an unknown function of , and is the second derivative of with respect to .

  • The process, in mathematical language, is to construct an integral of the inner product of the residual and the weight functions and set the integral to zero.

  • This is achieved by a particular space discretization in the space dimensions, which is implemented by the construction of a mesh of the object: the numerical domain for the
    solution, which has a finite number of points.

  • XFEM[edit] Main article: Extended finite element method The extended finite element method (XFEM) is a numerical technique based on the generalized finite element method (GFEM)
    and the partition of unity method (PUM).

  • However, the derivative exists at every other value of , and one can use this derivative for integration by parts.

  • The function is the unique function of whose value is at and zero at every .

  • After this second step, we have concrete formulae for a large but finite-dimensional linear problem whose solution will approximately solve the original BVP.

  • Illustrative problems P1 and P2[edit] The following two problems demonstrate the finite element method.

  • We define a new operator or map by using integration by parts on the right-hand-side of (1): where we have used the assumption that .

  • To solve a problem, the FEM subdivides a large system into smaller, simpler parts called finite elements.

  • [1] The simple equations that model these finite elements are then assembled into a larger system of equations that models the entire problem.

  • To explain the approximation in this process, the finite element method is commonly introduced as a special case of Galerkin method.

  • Hence the convergence properties of the GDM, which are established for a series of problems (linear and nonlinear elliptic problems, linear, nonlinear, and degenerate parabolic
    problems), hold as well for these particular FEMs.

  • General form of the finite element method[edit] In general, the finite element method is characterized by the following process.

  • If this condition is not satisfied, we obtain a nonconforming element method, an example of which is the space of piecewise linear functions over the mesh, which are continuous
    at each edge midpoint.

  • Under specific hypotheses (for instance, if the domain is convex), a piecewise polynomial of order method will have an error of order .

  • Main article: Applied element method A-FEM[edit] Yang and Lui introduced the Augmented-Finite Element Method, whose goal was to model the weak and strong discontinuities without
    needing extra DoFs, as PuM stated.

  • [24] Finite Element Model of a human knee joint[25] This powerful design tool has significantly improved both the standard of engineering designs and the design process methodology
    in many industrial applications.

  • Comparison to the finite difference method The finite difference method (FDM) is an alternative way of approximating solutions of PDEs.

  • The name virtual derives from the fact that knowledge of the local shape function basis is not required and is, in fact, never explicitly calculated.

  • • There are reasons to consider the mathematical foundation of the finite element approximation more sound, for instance, because the quality of the approximation between
    grid points is poor in FDM.

  • Various numerical solution algorithms can be classified into two broad categories; direct and iterative solvers.

  • [11] The method has since been generalized for the numerical modeling of physical systems in a wide variety of engineering disciplines, e.g., electromagnetism, heat transfer,
    and fluid dynamics.

  • [12][13] Technical discussion The structure of finite element methods[edit] A finite element method is characterized by a variational formulation, a discretization strategy,
    one or more solution algorithms, and post-processing procedures.

  • It includes the use of mesh generation techniques for dividing a complex problem into small elements, as well as the use of software coded with a FEM algorithm.

  • For higher-order partial differential equations, one must use smoother basis functions.

  • • The quality of a FEM approximation is often higher than in the corresponding FDM approach, but this is highly problem-dependent, and several examples to the contrary can
    be provided.

  • Hrennikoff’s work discretizes the domain by using a lattice analogy, while Courant’s approach divides the domain into finite triangular subregions to solve second order elliptic
    partial differential equations that arise from the problem of torsion of a cylinder.

  • Choosing a basis[edit] Interpolation of a Bessel function 16 scaled and shifted triangular basis functions (colors) used to reconstruct a zeroeth order Bessel function J0
    (black) The linear combination of basis functions (yellow) reproduces J0 (black) to any desired accuracy.

  • Generally, the higher the number of elements in a mesh, the more accurate the solution of the discretized problem.

  • Post-processing procedures are designed to extract the data of interest from a finite element solution.

  • The practical application of FEM is known as finite element analysis (FEA).

  • It is increasingly being adopted by other commercial finite element software, with a few plugins and actual core implementations available.

  • In step (2) above, a global system of equations is generated from the element equations by transforming coordinates from the subdomains’ local nodes to the domain’s global
    nodes.

  • In contrast, computational fluid dynamics (CFD) tend to use FDM or other methods like finite volume method (FVM).

  • If instead of making h smaller, one increases the degree of the polynomials used in the basis function, one has a p-method.

  • High-order methods with large uniform p are called spectral finite element methods (SFEM).

  • Mesh adaptivity may utilize various techniques; the most popular are: • moving nodes (r-adaptivity) • refining (and unrefined) elements (h-adaptivity) • changing order of
    base functions (p-adaptivity) • combinations of the above (hp-adaptivity).

 

Works Cited

[‘Daryl L. Logan (2011). A first course in the finite element method. Cengage Learning. ISBN 9780495668275.
2. ^ Jump up to:a b Reddy, J. N. (2006). An Introduction to the Finite Element Method (Third ed.). McGraw-Hill. ISBN 9780071267618.
3. ^ Huebner,
Kenneth H. (2001). The Finite Element Method for Engineers. Wiley. ISBN 978-0-471-37078-9.
4. ^ Liu, Wing Kam; Li, Shaofan; Park, Harold S. (2022). “Eighty Years of the Finite Element Method: Birth, Evolution, and Future”. Archives of Computational
Methods in Engineering. 29 (6): 4431–4453. arXiv:2107.04960. doi:10.1007/s11831-022-09740-9. ISSN 1134-3060. S2CID 235794921.
5. ^ Hrennikoff, Alexander (1941). “Solution of problems of elasticity by the framework method”. Journal of Applied Mechanics.
8 (4): 169–175. Bibcode:1941JAM…..8A.169H. doi:10.1115/1.4009129.
6. ^ Courant, R. (1943). “Variational methods for the solution of problems of equilibrium and vibrations”. Bulletin of the American Mathematical Society. 49: 1–23. doi:10.1090/s0002-9904-1943-07818-4.
7. ^
“СПб ЭМИ РАН”. emi.nw.ru. Archived from the original on 30 September 2015. Retrieved 17 March 2018.
8. ^ Hinton, Ernest; Irons, Bruce (July 1968). “Least squares smoothing of experimental data using finite elements”. Strain. 4 (3): 24–27. doi:10.1111/j.1475-1305.1968.tb01368.x.
9. ^
“SAP-IV Software and Manuals”. NISEE e-Library, The Earthquake Engineering Online Archive. Archived from the original on 2013-03-09. Retrieved 2013-01-24.
10. ^ Gard Paulsen; Håkon With Andersen; John Petter Collett; Iver Tangen Stensrud (2014).
Building Trust, The history of DNV 1864-2014. Lysaker, Norway: Dinamo Forlag A/S. pp. 121, 436. ISBN 978-82-8071-256-1.
11. ^ Strang, Gilbert; Fix, George (1973). An Analysis of The Finite Element Method. Prentice Hall. ISBN 978-0-13-032946-2.
12. ^
Olek C Zienkiewicz; Robert L Taylor; J.Z. Zhu (31 August 2013). The Finite Element Method: Its Basis and Fundamentals. Butterworth-Heinemann. ISBN 978-0-08-095135-5.
13. ^ Bathe, K.J. (2006). Finite Element Procedures. Cambridge, MA: Klaus-Jürgen
Bathe. ISBN 978-0979004902.
14. ^ celledoni (2023-02-27). “CutFEM: Discretizing Partial Differential Equations and Geometry”. ECMI. Retrieved 2023-10-13.
15. ^ Burman, Erik; Claus, Susanne; Hansbo, Peter; Larson, Mats G.; Massing, André (2015-11-16).
“CutFEM: Discretizing geometry and partial differential equations”. International Journal for Numerical Methods in Engineering. 104 (7): 472–501. Bibcode:2015IJNME.104..472B. doi:10.1002/nme.4823. ISSN 0029-5981.
16. ^ Babuška, Ivo; Banerjee, Uday;
Osborn, John E. (June 2004). “Generalized Finite Element Methods: Main Ideas, Results, and Perspective”. International Journal of Computational Methods. 1 (1): 67–103. doi:10.1142/S0219876204000083.
17. ^ P. Solin, K. Segeth, I. Dolezel: Higher-Order
Finite Element Methods, Chapman & Hall/CRC Press, 2003
18. ^ Song, Chongmin; Wolf, John P. (5 August 1997). “The scaled boundary finite-element method – alias consistent infinitesimal finite-element cell method – for elastodynamics”. Computer Methods
in Applied Mechanics and Engineering. 147 (3–4): 329–355. Bibcode:1997CMAME.147..329S. doi:10.1016/S0045-7825(97)00021-2.
19. ^ “Spectral Element Methods”. State Key Laboratory of Scientific and Engineering Computing. Archived from the original
on 2017-08-10. Retrieved 2017-07-28.
20. ^ Beirão da Veiga, L.; Brezzi, F.; Cangiani, A.; Manzini, G.; Marini, L. D.; Russo, A. (2013). “Basic principles of Virtual Element Methods”. Mathematical Models and Methods in Applied Sciences. 23 (1): 199–214.
doi:10.1142/S0218202512500492.
21. ^ Jump up to:a b c Topper, Jürgen (January 2005). “Option pricing with finite elements”. Wilmott. 2005 (1): 84–90. doi:10.1002/wilm.42820050119 (inactive 2024-04-07). ISSN 1540-6962.
22. ^ “What’s The Difference
Between FEM, FDM, and FVM?”. Machine Design. 2016-04-18. Archived from the original on 2017-07-28. Retrieved 2017-07-28.
23. ^ Kiritsis, D.; Eemmanouilidis, Ch.; Koronios, A.; Mathew, J. (2009). “Engineering Asset Management”. Proceedings of the
4th World Congress on Engineering Asset Management (WCEAM): 591–592.
24. ^ “Finite Element Analysis: How to create a great model”. Coventive Composites. 2019-03-18. Retrieved 2019-04-05.[permanent dead link]
25. ^ Naghibi Beidokhti, Hamid; Janssen,
Dennis; Khoshgoftar, Mehdi; Sprengers, Andre; Perdahcioglu, Emin Semih; Boogaard, Ton Van den; Verdonschot, Nico (2016). “A comparison between dynamic implicit and explicit finite element simulations of the native knee joint” (PDF). Medical Engineering
& Physics. 38 (10): 1123–1130. doi:10.1016/j.medengphy.2016.06.001. PMID 27349493. Archived (PDF) from the original on 2018-07-19. Retrieved 2019-09-19.
26. ^ Jump up to:a b c Hastings, J. K., Juds, M. A., Brauer, J. R., Accuracy and Economy of
Finite Element Magnetic Analysis, 33rd Annual National Relay Conference, April 1985.
27. ^ McLaren-Mercedes (2006). “McLaren Mercedes: Feature – Stress to impress”. Archived from the original on 2006-10-30. Retrieved 2006-10-03.
28. ^ Peng Long;
Wang Jinliang; Zhu Qiding (19 May 1995). “Methods with high accuracy for finite element probability computing”. Journal of Computational and Applied Mathematics. 59 (2): 181–189. doi:10.1016/0377-0427(94)00027-X.
29. ^ Haldar, Achintya; Mahadevan,
Sankaran (2000). Reliability Assessment Using Stochastic Finite Element Analysis. John Wiley & Sons. ISBN 978-0471369615.
Photo credit: https://www.flickr.com/photos/aidanmorgan/3556123230/’]