What Are Shape Funct Derived? ns In FEA - And How Are They Many of us have utilized the power of FEA solutions, such as Abaqus, to analyze all manner of products. But how much do we really know about exactly how it works? In this article, we will learn how to discretize a simple 1D domain into linear and quadratic elements and derive shape functions by approximating the unknown variables over the elements. Mathematical Modeling ‘As we know, any physical problem can be expressed in terms of a mathematical governing equation. These governing equations can be simple algebraic equations (heat conduction between two points) or complex differential and functional descriptions (spring-mass system, Navier-Stokes). The exact numerical solution to these governing equations cannot be always obtained using standard solution methods like substitutions and separable equations. Finite element analysis (FEA) is a procedure to obtain approximate numerical solution to governing equations. This is achieved by converting the mathematical model to set of solvable algebraic equations and attaining the numerical solution at set of discrete locations in the ‘computational domain. Deriving Weak Form Once the governing diferential equation (strong form) is obtained by considering the physics, kinematics and dynamics of a physical problem, the weak form can be obtained using different approaches like virtual work principle and Galerkin weighted residual method. For example, the weak form of 1D elastic problem is given below, j ® a (#)ax = | Sub(x) dx ‘The figure below summarized different techniques for obtaining weak form equations. Potential Energy Euler equation Virtual Work ‘Functional statement Strong Form. ‘Vantinalcaeulue Derivation of Shape FunctionsThe above weak form equation is solved in a discretized computational domain. The figure below shows the discretization of 1D domain into 3 elements. Elements Nodes Let us assume the unknown displacement u at the nodes are u1, u2, u3 and u4. This unknown variable u is approximated using a polynomial function in x as given below, UN) a5. 14 09.4402. as. 0 +s x where 1, x, x2, x3 are called monomial basis functions and a1, a2, a3, a4 are coefficients. 1 el 2 x For linear elements, a linear function is assumed for u. ux) Since there are two unknown coefficients, we need two boundary conditions to solve them. (xs) (xz) ‘Substituting the boundary conditions in the linear approximation gives, 1 1402.x A Us= Gr +02. X Up = ay +42. X2 This can be written in matric form as, t et3-(3The above equations are solved for unknown coefficients a1, and a2 which gives, uf) = Na (x) us + Nox) uz where, N,Q) = 22 NG) ==* Here N1 (x), N2(x) are called linear shape functions. Z~ 1 et 2 x Similarly, if we have three node elements as shown in the below figure, a quadratic polynomial is assumed for u. Here, we will need three nodal displacements to solve for the unknown coefficients. el The coefficients are obtained by solving the equations below, u(x) [email protected]+a2.x+03.x° Expressing u in terms of unknown coefficients gives, Tox d)]imy cy 1 x, xf ]\42) Ms. &—me—m) NO = Gane =x) @ Oe where, Nz) = _ = Hex) MO) = Tega)Here N4(x), N2(x) and N3{x) are called quadratic shape functions. 7 a ~~ The shape functions are to be calculated for all the elements in the discretized domain. They are substituted back into the weak form governing equations to perform finite element approximation and to derive element matrices and vectors. The above procedure can be extended to derive shape functions for 2D, 3D linear and quadratic elements. Ni(e2)Nale2)p¥i(e3)_ Male Properties of Shape Functions 1, Ata node, the value of shape function is either 0 or 1 m@={9 i2) 2, Summation of all the shape functions at any point is equal to 1 Yueo=t 3. Continuity Within the element boundary, the polynomial approximation for u is C~ continuous, The continuity across the element boundaries depends on the polynomial approximation for the variable. For example, the variable has only CO continuity across the linear and quadratic element boundaries (.e., its first derivative is piecewise continuous).Us ua a | x continuity for linear elements 4. Completeness This is an important requirement for shape functions to ensure the convergence of finite element approximate solution to the exact solution of governing equation Final Thoughts. Shape functions are the backbone of finite element method as they play a crucial part in converting the weak form equation to set of computer solvable algebraic equations. Hopefully this article has provided some basic understanding of discretization of a 1D computational domain and derivation of shape functions.