Verification and validation of a numerical simulation
Introduction
The main goal of this course is to show you how to deploy a finite element simulation. The FE simulations should be updated, deployed, and optimized depending of the purpose of the simulations.
​
We will present the concepts through a series of examples.
Difference between verification and validation
The objective of a simulation is to predict physical phenomena. But the question is how can you make sure that your numerical model works well? And the other question is that what does it mean “works well”?
The first question can be answered through verification and validation steps. With the verification step, you make sure that your model is mathematically correct. And with the validation step, you ensure that your model is physically correct.
For instance in verification, you compare your simulation predictions with the predictions of an analytical model. Imagine that you want to model a beam by Finite Element (FE) simulation. Your FE model predicts the deformation of a cantilever beam under a uniform and static loading. Then, you also need to compute the beam deformation by analytical model of the Euler-Bernoulli beam or Timoshenko beam. You compare the predicted deformation by your FE simulation with the predicted deformation by the Euler-bernoli beam. If the error is less than your expected error, then you can say that Yes my model works well! (congratulations!).
Besides that, you also need to make sure that you have a fine mesh. So you should do a convergence study. In convergence study, you reduce the element size to have a converged system response quantity. For instance, in our case you reduce the element size (whether it’s solid element or beam element) and then you plot the deformation at the tip of the beam with respect to the total degrees of freedom.
​In validation, you compare the model’s predictions with the experimental measurement. This time, it is much more difficult. First, an experimental setup should be designed and it should be defined what and how you want to measure? Let’s say like the verification step, you want to measure the vertical displacement at the tip of the beam (it would be called system response quantity). Then, you should define how you want to measure and what is the accuracy of the measurement. It would matter whether you want to measure it with a three thousand dollar camera or with a three dollar plastic ruler. You will decide. And of course, your decision will affect your validation results.
The verification results and validation results indicate how your model performs. It can work well or not. Now the second question that was asked at the beginning should be answered. What does it mean that the model works well? Well, it depends again on you. There are no good or bad results. There are just numbers that describe the model performance. The results can be acceptable or not. If you want to design a blade in a jet engine the error of simulation should be less and 1 millimeter. If you want to do an assignment, maybe 10 mm is still ok (it depends on your TA). So the application of the modeling in fact determines the expected accuracy of the model.
Example #1
Example 1 is modelling of a cantilever beam (the following figure). It is a very simple finite element simulation. You should write a python code to simulate a cantilever beam in Abaqus.
​
​
Write a code to simulate a cantiliver beam. The left side of the beam is fixed (all DOF of U1, U2 and UR1) are equal to zero. A concentrated force of F is applied on the right side of the beam in the tip. The lenth of the beam is 1000 mm. The Young's modulus of E=200 Gpa and Poisson's ration of v=0.3.
​
The force F=100 N is applied on the tip of the beam (the force is negative and it is in direction of y).
​
The objective of simulation is to compute the maximum deflection of the beam under this loading.
Verification
Verify you FE simulation by doing a convergence study. The convergence criterion is 1%. Reduce the element size from 100 mm to 1 mm. Plot the convergence study.
​
​
Validation
Validate the prediction of your simulation ​with the experimental measurement.
Compute the relative difference of the prediction agains the experimental measurement.
​
​
Sensitivity analysis
Do a sensitivity analysis on the length and the cross section of the beam. Plot the maximum stress with respect to the different length of the beam and also with respect to different aspect ratio of a/b. ​
​
Solution
We should reduce the element size from 100 to 1 mm and we stop the reduction until we satisfy the convergence criteria. For each mesh, the maximum displacement is computed and then we compute the relative different between the coarse mesh and the fine mesh.