Finite Element Methods 1 (NMNV405) — Practicals (Winter Semester 2023/2024)
Lectures
Monday 09:00 – 10:30, K7 Sokolovská 83 Karlín
- Lecture 1 Notes
- Lecture 2 Notes
- Lecture 3 Notes
- Lecture 4 Notes
- Lecture 5 Notes
- Lecture 6 Notes
- Lecture 7 Notes
- Lecture 8 Notes
- Lecture 9 Notes
- Lecture 10 Notes
- Lecture 11 Notes
- Lecture 12 Notes
- Lecture 13 Notes
Practicals
Monday 10:40 – 12:10, K7 Sokolovská 83 Karlín
- Practical 2 Notes
- Practical 3 Notes
- Practical 4 Notes
- Practical 5 Notes
- Practical 6 Notes
- Practical 7 Notes
- Practical 8 Notes
- Practical 9 Notes
- Practical 10 Notes
- Practical 11 Notes
- Practical 12 Notes
- Practical 13 Notes
Homework
There will be four homeworks during the course of the year. Obtaining credit for this course will involve obtaining at least 50% of the marks available from these homeworks.
- 23.10.2023: Homework (Deadline: 06.11.2023); Solutions
- 06.11.2023: Homework (Deadline: 20.11.2023); Solutions
- 27.11.2023: Homework (Deadline: 11.12.2023); Solutions
- 18.12.2023: Homework (Deadline: 08.01.2024); Solutions
Exam
The course examination will consist of a 30 minute oral examination on the syllabus covered in the lecture.
Finite Elements and Basis Functions
Some examples of finite elements can be found in these notes
Examples of various basis functions on triangles and rectanges can be found here.
The displayed basis function can be changed via the selection boxes in the top left corner. The bottom left corner allows switching the view between 3D (orthogonal and perspective) and 2D views. In 3D view the basis function can be rotated by dragging. In the bottom left corner the principal lattice (white lines) can be toggled on/off.
The following website also lists details about all the basis functions we will discuss (plus a lot more): defelement.com
Notes on Gaussian Quadrature
The following journal article has information on the choice of points and weights for Gaussian quadrature on a triangle:
Dunavant, D.A. (1985), High degree efficient symmetrical Gaussian quadrature rules for the triangle. Int. J. Numer. Meth. Engng., 21:1129-1148. https://doi.org/10.1002/nme.1620210612
Example FEM Code
An example MATLAB finite element code which solves the problem
-Δu = f,
on the unit square [0, 1]2 with forcing function
f = 2π2 sin πx sin πy,
and homogeneous Dirichlet boundary conditions is available in the following zip file
or below:- Main Function: fem.m
- Convergence Analysis: fem_convergence.m
- Quadrature: quadrature_tri.m, quadrature_rect.m
- Lagrange Basis Functions (1D Interval): basis_lagrange_interval.m, deriv_basis_lagrange_interval.m
- Lagrange Basis Functions (Triangle): basis_lagrange_tri.m, grad_basis_lagrange_tri.m
- Lagrange Basis Functions (Rectangle): basis_lagrange_quad.m, grad_basis_lagrange_quad.m
- Reduced Lagrange Basis Functions (Triangle): basis_reduced_tri.m, grad_basis_reduced_tri.m
- Reduced Lagrange Basis Functions (Rectangle): basis_reduced_quad.m, grad_basis_reduced_quad.m
- Crouzeix-Raviart Basis Functions: basis_crouzeix_raviart.m, grad_basis_crouzeix_raviart.m
- Rotated Bilinear Basis Functions: basis_rotated_bilinear.m, grad_basis_rotated_bilinear.m
- Plotting Utility Functions: plot_mesh.m, plot_soln.m