Nonlinear finite elements/Axial bar weak form

Axially loaded bar: Weak FormEdit

Instead of deriving the differential equation using a balance of forces on a differential element, we may arrive at the same problem description via a different route - the principle of virtual work (also called the variational approach).

If we imagine there were forces (virtual forces) inside and outside of the bar, then the virtual work generated by these 'virtual forces' should conserve energy. For the bar, this principle can be stated as

 

where   is the virtual work of the internal forces,   is the virtual work of the external forces, and   is the virtual work of the body forces.

The virtual internal work is given by

 

The virtual work done by the external forces is given by

 

The virtual work done by the body forces is given by

 

The principle of virtual work for the bar can then be expressed as

 

Now, the stress and the virtual strain expressed in terms of the displacements are

 

Therefore, we have

 

Note that the virtual displacement   is zero at points on the boundary where displacements are prescribed.

The above equation is called the variational form or weak form of the problem.

Why is it called a variational form?Edit

Let us start with the weak form. Using the formula for the first variation (from variational calculus)

 

we have

 

Therefore, the weak form can be written as

 

This is equivalent to the following variational statement of the problem

 

where   means variation in and   is an arbitrary variation on   subject to the condition that  .


Principle of minimum potential energyEdit

Note that the first term of   is the strain energy stored and the next two terms are the potential energy of the external forces (  ). The above variational statement can also be interpreted as statement of minimum potential energy: the total potential of the system (elastic body and external forces) attains a stationary value when equilibrium is satisfied. Further more, this stationary value is minimum. This can be reinterpreted as of all the possible displacements, the actual displacement   minimizes the potential energy of the system.

Equivalence of the strong and weak formsEdit

The strong and weak statements of the problem are equivalent. Since finite elements uses the weak form of the problem statement, we can either derive the weak form from the strong form or using the principle of virtual work. The preferred approach is chosen on the basis of convenience.

We can derive the weak form from the strong form as follows. We multiply the strong form with an arbitrary virtual displacement (or weighting function) and integrate over the length of the bar. Thus

 

Next, we integrate the first term of the equation by parts (to get rid of the higher order derivatives of  ) and get

 

Recall that   and  . Therefore, we have

 

This is the same as the weak form derived from the principle of virtual work.