Numerical Integration
Numerical Integration
1
Finite Element Method by G. R. Liu and S. S. Quek
•
1D Jacobian
1/17/2017 Dr.Yilma T. 2
Finite Element Method by G. R. Liu and S. S. Quek
1/17/2017 Dr.Yilma T. 3
Finite Element Method by G. R. Liu and S. S. Quek
2D Jacobian
• For a continuous 1-to-1 transformation from
(x,y) to (u,v), then
2D Jacobian
Where Region (in the xy plane) maps onto region in the uv plane
maps areas dxdy to areas dudv
1/17/2017 Dr.Yilma T. 4
Finite Element Method by G. R. Liu and S. S. Quek
1/17/2017 Dr.Yilma T. 5
Finite Element Method by G. R. Liu and S. S. Quek
Relation between Jacobians
1/17/2017 Dr.Yilma T. 6
Finite Element Method by G. R. Liu and S. S. Quek
3D Jacobian
1/17/2017 Dr.Yilma T. 7
Finite Element Method by G. R. Liu and S. S. Quek
Jacobian Matrix
It is used to transformer the derivatives of shape function
from global to local coordinate system
1/17/2017 Dr.Yilma T. 8
Inverse of Jacobian Matrix
1/17/2017 Dr.Yilma T. 9
Elements of Jacobian Matrix
1/17/2017 Dr.Yilma T. 10
,
1/17/2017 Dr.Yilma T. 11
FEM Numerical Integration
Two types of elements
• Shape functions
• Spatial derivatives of
defined here
shape functions
needed here (in
matrices K ) • Numerical integration
performed on local
element (using Gauss-
• Integration
Legendre quadrature)
boundaries given
1/17/2017
here. Dr.Yilma T. 12
Transformation of spatial
derivatives
Global element Local element
• Spatial derivatives of • Spatial derivatives of
shape functions
shape functions
Transformation
1/17/2017 Dr.Yilma T. 13
Transformation of integration boundaries
Local element
Global element
• Integration boundaries
• Integration
boundaries
Tranasformation
1/17/2017 Dr.Yilma T. 14
Numerical integration
1/17/2017 Dr.Yilma T. 15
Gauss-Legendre quadrature
• The integeration is carried between 1 and -1
• In 1D (n=3):
• In 2D (nip=9):
1/17/2017 Dr.Yilma T. 16
Gauss-Legendre quadrature
1/17/2017 Dr.Yilma T. 17
Integration points and weights
1/17/2017 Dr.Yilma T. 18
More on Gauss-Legendre-Quadrature
1/17/2017 Dr.Yilma T. 19
Distorted vs. idealized element
1/17/2017 Dr.Yilma T. 20
Two transformations are necessary
1/17/2017 Dr.Yilma T. 21
First transformation in 1D
1/17/2017 Dr.Yilma T. 22
First transformation in 1D
→ Derivatives of shape fcts. from local to global
1/17/2017 Dr.Yilma T. 23
Second transformation in 1D
1/17/2017 Dr.Yilma T. 24
First transformation in 2D
1/17/2017 Dr.Yilma T. 25
First transformation in 2D
1/17/2017 Dr.Yilma T. 26
First transformation in 2D
→ Derivatives of shape fcts. from local to global
1/17/2017 Dr.Yilma T. 27
Second transformation in 2D
1/17/2017 Dr.Yilma T. 28