All courses › Finite Element Method › Master element and Jacobian

Master element and Jacobian

All elements are computed on the same master element ξ∈[−1,1]\xi\in[-1, 1] and mapped to the physical element. The Jacobian JJ is how much the element is stretched. Integrals get the factor JJ, derivatives are divided by JJ.

φ1=1−ξ2,φ2=1+ξ2\varphi_1 = \frac{1 - \xi}{2},\quad \varphi_2 = \frac{1 + \xi}{2}linear shape functions
x(ξ)=x1φ1+x2φ2,J=h2x(\xi) = x_1\varphi_1 + x_2\varphi_2,\quad J = \frac{h}{2}mapping and Jacobian
Bije=∫−11EA dφidξdφjdξ1J dξB^e_{ij} = \int_{-1}^{1} EA\,\frac{d\varphi_i}{d\xi}\frac{d\varphi_j}{d\xi}\frac{1}{J}\,d\xielement matrix
Fie=∫−11Fx φi J dξF^e_i = \int_{-1}^{1} F_x\,\varphi_i\,J\,d\xielement vector

Symbols

ξ\xicoordinate on the master element
JJJacobian dx/dξdx/d\xim
hhelement lengthm

Example

h=2h = 2 m, EA=600EA = 600 kN: J=1J = 1, dφ1dξ=−12\frac{d\varphi_1}{d\xi} = -\tfrac12, so B11e=∫−11600⋅14 dξ=300B^e_{11} = \int_{-1}^{1} 600\cdot\tfrac14\,d\xi = 300 kN/m =EA/h= EA/h.

Chain rule: dφdx=1Jdφdξ\frac{d\varphi}{dx} = \frac1J\frac{d\varphi}{d\xi}. Forget JJ and the stiffness is off by a factor h/2h/2.
Practise weak form and shape functions for free →

← The weak form for the bar · Shape functions →

Part of Finite Element Method: Weak form and shape functions.