Êtes-vous un étudiant de l'EPFL à la recherche d'un projet de semestre?
Travaillez avec nous sur des projets en science des données et en visualisation, et déployez votre projet sous forme d'application sur Graph Search.
Concept
Stiffness matrix
Résumé
In the finite element method for the numerical solution of elliptic partial differential equations, the stiffness matrix is a matrix that represents the system of linear equations that must be solved in order to ascertain an approximate solution to the differential equation.
For simplicity, we will first consider the Poisson problem
on some domain Ω, subject to the boundary condition u = 0 on the boundary of Ω. To discretize this equation by the finite element method, one chooses a set of basis functions {φ_1, ..., φ_n} defined on Ω which also vanish on the boundary. One then approximates
The coefficients u_1, u_2, ..., u_n are determined so that the error in the approximation is orthogonal to each basis function φ_i:
The stiffness matrix is the n-element square matrix A defined by
By defining the vector F with components the coefficients u_i are determined by the linear system Au = F. The stiffness matrix is symmetric, i.e. A_ij = A_ji, so all its eigenvalues are real. Moreover, it is a strictly positive-definite matrix, so that the system Au = F always has a unique solution. (For other problems, these nice properties will be lost.)
Note that the stiffness matrix will be different depending on the computational grid used for the domain and what type of finite element is used. For example, the stiffness matrix when piecewise quadratic finite elements are used will have more degrees of freedom than piecewise linear elements.
Determining the stiffness matrix for other PDEs follows essentially the same procedure, but it can be complicated by the choice of boundary conditions. As a more complex example, consider the elliptic equation
where is a positive-definite matrix defined for each point x in the domain. We impose the Robin boundary condition
where ν_k is the component of the unit outward normal vector ν in the k-th direction. The system to be solved is
as can be shown using an analogue of Green's identity. The coefficients u_i are still found by solving a system of linear equations, but the matrix representing the system is markedly different from that for the ordinary Poisson problem.
Source officielle
À propos de ce résultat
Cette page est générée automatiquement et peut contenir des informations qui ne sont pas correctes, complètes, à jour ou pertinentes par rapport à votre recherche. Il en va de même pour toutes les autres pages de ce site. Veillez à vérifier les informations auprès des sources officielles de l'EPFL.