Finite Element Method: Weak form to Algebraic Equations?

Click For Summary
SUMMARY

The discussion focuses on the application of the Finite Element Method (FEM) in engineering, specifically addressing the mathematical foundation using a one-dimensional heat transfer problem. The governing equation derived is k(d²T/dx²) + Q = 0, where k represents thermal conductivity, T is temperature, and Q is the heat generated. The discretization process involves the weighted integral form ∫₀ᴸ w_e(k(d²T/dx²) + Q) = 0, where w is an arbitrary linear weighting function. The necessity of the weighting term is emphasized, as it facilitates the transition to a variational expression, allowing for a matrix representation that simplifies programming.

PREREQUISITES
  • Understanding of the Finite Element Method (FEM)
  • Familiarity with partial differential equations (PDEs)
  • Knowledge of variational principles in mathematical modeling
  • Basic programming skills for implementing FEM algorithms
NEXT STEPS
  • Study the derivation of variational principles in FEM
  • Learn about the implementation of weighted residual methods
  • Explore numerical methods for solving PDEs
  • Investigate software tools for FEM analysis, such as ANSYS or COMSOL Multiphysics
USEFUL FOR

Engineers, researchers, and students interested in computational mechanics, particularly those focusing on heat transfer analysis and the application of the Finite Element Method in engineering problems.

ramzerimar
Messages
178
Reaction score
23
Okay, I'm following a series of video lectures on applications of finite element method to engineering, and the tutor started by demonstrating the mathematical background of FEM using a simple heat transfer problem. He derived the governing equation (in just one dimension):

(1) k\frac{d^2 T}{dx^2} + Q = 0

Where K is a constant, T is temperature and Q is the heat generated. The next step was discretizing the domain (in this case, a bar with length L). For this, we used the weighed integral form, which is:

(2) \int_{0}^{L}w_e(k\frac{d^2 T}{dx^2} + Q) = 0

Where w is a arbitrary linear weighting function. I understand that we can't solve (1) by using this discretization, because the temperatures would be discontinuous at the nodes and the second derivative wouldn't be defined, and that's why we integrate it to get only first derivative terms. But I didn't understand the weighting term. Just integrating it wouldn't be enough? What's the purpose of it?
 
Last edited:
Physics news on Phys.org
I think the author is backing his way into a variational expression for the equation of motion. The equation written is valid for all reasonable choices of the test function. This allows for a simple matrix form of the equation being modeled that is easy to program.
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 6 ·
Replies
6
Views
4K
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 10 ·
Replies
10
Views
5K
  • · Replies 3 ·
Replies
3
Views
1K
  • · Replies 1 ·
Replies
1
Views
3K
Replies
1
Views
1K
  • · Replies 2 ·
Replies
2
Views
3K