Galerkin's Method: APL705 Finite Element Method
Galerkin's Method: APL705 Finite Element Method
Galerkin's Method: APL705 Finite Element Method
Lecture 5: Weighted Residual Methods:
Galerkin’s Method
APL705 Finite Element Method
Weighted Residual Methods
• Here we start with a set of governing differential
equations.
• Consider a set of governing equations on region V as
Lu = P
• Here L is an operator operating on the displacements u. In the
previous example of one‐dimensional problem this will be
d ⎛ du ⎞
⎜ EA ⎟ = 0
dx ⎝ dx ⎠
• Here the operator L is replaced by d d
EA ( ) = 0
operating on u dx dx
1
1/24/2017
Galerkin’s Method
• This is a weighted residual method
To find an exact solution forLu = P , it should be satisfied at all x;
To find an exact solution for , it should be satisfied at all x;
however, for an approximate solution , it will introduce an error ε(x),
called the residual which can be expressed as follows:
The approximation methods attempt to make the residual zero relative
to a weighting function Wi as
Depending on the choice of a weighting function Wi gives rise to
various methods. In Galerkin’s method, weighting function Wi is
chosen from the basis function used to construct .
Galerkin’s Method
• The method is outlined below
Choose a basis function G
Ch b i f ti Gi and determine coefficients Q
dd t i ffi i t Qi in
i
n
~
u = QG ∑
i =1
i i
Such that
∫ φ (Lu~ − P)dV = 0
V
n
For every φ
For of the type φ = ∑φiGi where coefficients φ
every φ of the type where coefficients φi are
are
i =1
arbitrary except that φ has to satisfy the boundary homogeneous
boundary conditions. The solution of the resulting equations Qi
then gives the approximate solution .
2
1/24/2017
Galerkin’s Method: Illustrative example
Assuming a polynomial function u = a + a x + a x 2
1 2 3
This function should be satisfied at the two boundaries: u=0 at
This function should be satisfied at the two boundaries: u=0 at
x=0 and x=2. Therefore 0 = a1 and 0 = +2a2 + 4a3
∴ a2 = −2a3
u = a3 (−2x + x 2 ) ⇒ u1 = −a3
Next we will find he potential energy
du
= a3 (−1+
( 1+ x))
dx
2
π = ∫ 4a32 (−1+ x)2 dx − 2(−a3 )
0
2
= 2a32 ∫ (1− 2x + x 2 )dx + 2a3 = 2a32 ( 23 ) + 2a3
0
Rayleigh‐Ritz Method: Illustrative example
∂π
Now to find minimum potential energy we se = 2a3 ( 23 ) + 2 = 0
∂a3
a3 = −0.75
0 u1 = −a3 = 0.75
0
Next let us find the stress in the bar σ = E du = 1.5(1− x)
dx