ME 428: Finite Element Method: Dr. S. M. Kamal

ME 428: Finite Element Method
Lecture 4
Dr. S. M. Kamal
1. Introduction
In the last lecture, you have learnt how to obtain the variational form of a given boundary value
problem. We also discussed the virtual work principle in reference to obtain the variational form
particularly for structural problems. In the process of obtaining the variational form of the axial
bar problem you have seen that at the end of the procedure we arrived at the condition of
extremizing a functional along with the satisfaction of the boundary conditions. That functional
is the variational form or the integral form of the boundary value problem. The integral form
needs to be extremized to obtain the solution of the boundary value problem. This formulation is
called variational formulation. In this lecture, we will discuss the physical interpretation of the
variational form. We also discuss some problems.
2. Principle of minimum potential energy

We obtained the variational form of the axial bar problem as
L
 1  du 
2

I = ∫  EA   − qu  dx − P u x = L . (1)
0 2  dx  
Equation (1) has a physical interpretation. The quantity I in Eq. (1) represents the total potential
energy of the bar.

The total potential energy of the bar consists of (a) the strain energy of the elastic
deformation and (b) the potential possessed by the applied loads. The strain energy of a linear
elastic body is given by
1
U = ∫ σε dV , (2)
V
2
where σ is the stress and ε is the strain. For the problem of axial bar
du
ε= , (3)
dx
du
σ =E . (4)
dx
1
Thus, in the present case,
2
1  du 
U = ∫ E   dV ,
V
2  dx 
L 2
1  du 
U = ∫ EA   dx, (5)
0
2  dx 
where A is the cross sectional area and L is the length of the rod.
The work potential due to the axial distributed load q(x) is given by
L
Vq = − ∫ qudx. (6)
0
The work potential due to the concentrated load at the extreme end of the bar is
VP = − P u x = L . (7)
Thus, the total potential energy of the rod is given by

L 2 L
1  du 
∏ = U + Vq + VP = ∫ EA   dx − ∫ qudx − P u x = L . (8)
0
2  dx  0
Equation (8) is the same expression that we obtained for the variational form I in Eq. (1).
Therefore the variational formulation states the principle of minimum potential energy, which
can be stated as
“For conservative structural systems, of all the kinematically admissible deformations, those
corresponding to the equilibrium state, extremize the total potential energy. If the extremum is a
minimum, the equilibrium state is stable’.
3. Strong and weak form of a problem representing a physical system
We can represent a physical system either as a boundary value problem consisting of a
differential equation with boundary conditions or by an integral form known as the variational
form if obtained by using the variational calculus. The boundary value problem, i.e., the
differential form of the physical system is called the strong form. It is called as the strong form,
because it contains higher order derivatives. However, the integral form containing the lower
order derivatives is called the weak form.
2
4. Exercise problems
(i) The functional governing static buckling of the column in Fig. 4.1 is given as
2 2
1
L
 d2w  L
P  dw  1 2
∏ = ∫ EI  2  dx − ∫   dx + kwL (9)
2 0  dx  2 0  dx  2
dw
where w x = 0 = 0 and = 0. (10)
dx x = 0
Derive the governing differential equation of the problem.
Fig. 4.1 Static buckling of a column

(ii) Consider the following boundary value problem
d 2u
+ u = 1, 0 ≤ x ≤1 (11)
dx 2
du
with u (0) = 0 and = 0 at x = 1. Convert this problem into variational form.
dx
(iii) The bending of a beam is governed by the following differential equation:
d2  d 2 v  dmz
 zz 2  +
EI − qy = 0 (12)
dx 2  dx  dx
with the following essential boundary conditions:

dv
at x = 0, v = v∗ and = θ z∗ (13)
dx
and the following natural boundary conditions:
3
d2v ∗ d  d2v  ∗
at x = L, EI = M and  zz 2  + mz = −Vy
EI (14)
dx 2 z
dx  dx 
where mz is the distributed moment per unit length about z-axis and qy is the distributed force per
unit length in y-direction. Obtain the variational form of the problem.

ME 428: Finite Element Method: Dr. S. M. Kamal

Uploaded by

Copyright:

Available Formats

ME 428: Finite Element Method: Dr. S. M. Kamal

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

ME 428: Finite Element Method: Dr. S. M. Kamal

Uploaded by

Copyright:

Available Formats

ME 428: Finite Element Method

2. Principle of minimum potential energy

energy of the bar.

Thus, the total potential energy of the rod is given by

Fig. 4.1 Static buckling of a column

with the following essential boundary conditions:

You might also like