Method of Manufactured Solutions
Nowadays, we often use the tools of computational analysis in engineering design. For the results of such analyses to be believable, the tools that are used have to pass rigorous tests. There are two categories of tests involved:
- Verification tests where one checks whether the mathematical model has been implemented according to design.
- Validation tests where one checks how closely the correctly implemented mathematical model mimics reality.
A method of verification of a numerical code that has received some attention is is the method of manufactured solutions. A systematic application of the method for complex nonlinear constitutive models in solid mechanics is yet to be found in the literature. The fluid mechanics community, however, uses the method on a regular basis.
Last year, at the World Conference on Computational Mechanics in Los Angeles, Prof. R. C. Batra told me that he had been using the method to verify his codes for a while. However, he had never felt the need to publish the approach because of its seemingly trivial nature. In this post I will give you an example of the method of fictitious body forces used by Prof. Batra [1,2].
The method of fictitious body forces involves assuming a displacement field over the body; finding the body force, initial conditions, and boundary conditions that fit that solution; and then using that body force and the computed boundary conditions and initial conditions in the numerical algorithm to arrive at a solution. If the computed solution matches the assumed displacement, then we're good to go.
The example below was meant to provide a starting point for a student. Please let me know if you find it useful or whether the equations are confusing rather that clarifying.
A 1-D example:
Consider a body in which the material point moves to during a motion. The motion is governed by the equations
If we now consider a one-dimensional bar of the material, then
The left Cauchy-Green deformation tensor is
The Cauchy stress is given by (using a simplified constitutive model):
where is an elastic constant. Also, the Cauchy stress and the 1st Piola-Kirchhoff stress are related by
The one-dimensional momentum equation is
The method of manufactured solutions:
Let us try to apply the method of fictitious body forces to the bar problem.
The first step is to assume that the displacement field is given by (you can choose any other appropriate functional form)
and the deformation gradient is
The left Cauchy-Green deformation is
The first Piola-Kirchhoff stress is
Plugging and into the momentum equation, we get
Therefore, the body force is
If the bar is of length , the boundary conditions are
The initial conditions are
When you apply the body force and the initial and boundary conditions, you should get a solution that matches the chosen function.
Note that all these have been done in a Lagrangian configuration. You could alternatively do the same assuming a function and transforming the equations accordingly.
- R. C. Batra and X. Q. Liang.
Finite dynamic deformations of smart structures.
Computational Mechanics, 20:427-438, 1997.
- R. C. Batra and B. M. Love.
Multiscale analysis of adiabatic shear bands in tungsten heavy alloy particulate composites.
Int. J. Multiscale Computational Engineering, 4(1):95-114, 2006.