Mesh Refinement Near Discontinuities PDF
Mesh Refinement Near Discontinuities PDF
Mesh Refinement Near Discontinuities PDF
In this study stresses are measured at 5 fixed locations in a simple shape as the mesh size is
changed. The stress errors predicted by the error function are compared with ultimate stress
predicted from mesh refinement. CosmosDesigner 2008 SP5.0 FEA software is used for this
report using 2nd order shell elements.
Sharp Corner
Sharp Radius
The test shape a simple flat plate modeled at thickness. The model has a sharp radius and a
sharp corner. Test data will be taken near the sharp radius and at the sharp corner.
Disclaimer: This document is provided for educational purposes only. Pressure Vessel Engineering
Ltd. is not liable for its use.
The surface is split to provide test areas 1 through 5 at fixed locations regardless of the mesh size
used.
The applied load is a 500 psi force at the right edge. The left edge is fixed. The shell model is
meshed at thick.
The error plot for the 1 Mesh size results for areas 2, 3, 4 and 5 should not be used. Elements
with errors greater than 5% are within 1 element of the test node locations. The result for test
node 1 can be used.
The overall stress pattern for the 1/16 mesh size shown as Tresca stress intensity (P1-P3).
How accurate are these stress values?
Reported Error
Err Point Err Point Err Point Err Point Err Point
Mesh 1 2 3 4 5
1" 1.141 2.53 2.021 2.304 31.742
0.75" 1.113 0.538 1.954 1.14 28.884
0.5" 0.739 0.445 1.291 1.337 13.377
0.375" 0.568 0.532 0.982 0.594 46.642
0.25" 0.379 0.589 0.65 0.641 45.822
0.1875" 0.287 0.604 0.513 0.498 51.904
0.125" 0.185 0.487 0.335 0.326 62.154
0.0625" 0.089 0.324 0.156 0.179 79.544
Stress and Error results for the 5 test points. All data for points 1-4 meet the 5% acceptance
criteria, but points 2-5 were previously disqualified at the 1 mesh size due to being within 1
element of areas of >5% error. Point 5 fails the criteria for all mesh sizes.
4.5
120
4
100 3.5
y = 1.0958x + 113.52
3
Error %
80
Stress
Stress
2.5
60 Error
2
Linear (Stress)
40 1.5
1
20
0.5
0 0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Mesh Size
Point 1 all stress results show an error of less than 5%. The ultimate stress is extrapolated to
113.52 psi.
4.5
30
y = -3.525x + 30.482 4
25 3.5
3
Error %
20
Stress
Stress 2.5
15
Error 2
Linear (Stress) 1.5
10
1
5
0.5
0 0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Mesh Size
Point 2 all stress results show an error of less than 5%. Stress at 1 mesh size has been
disqualified as being within 1 element of a node with greater than 5% error.
4.5
30
y = -7.0118x + 31.47 4
25 3.5
Error %
20
Stress
Stress 2.5
15
Error 2
Linear (Stress) 1.5
10
1
5
0.5
0 0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Mesh Size
Point 2 again - stress results extrapolated from mesh sizes 0.5 or smaller. The ultimate stress is
extrapolated to 31.47 psi
4.5
350
y = 10.917x + 351.87 4
300
Stress 3.5
250 Error 3
Error %
Stress
Linear (Stress)
200 2.5
2
150
1.5
100
1
50 0.5
0 0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Mesh Size
Point 3 all stress results show an error of less than 5%. Stress at 1 mesh size has been
disqualified as being within 1 element of a node with greater than 5% error. The ultimate stress
is extrapolated to 351.87 psi.
4.5
350
y = 11.946x + 350.38 4
300
3.5
250 Stress 3
Error %
Stress
Error
200 2.5
Linear (Stress)
150 2
1.5
100
1
50
0.5
0 0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Mesh Size
Point 4 all stress results show an error of less than 5%. Stress at 1 mesh size has been
disqualified as being within 1 element of a node with greater than 5% error. The ultimate stress
is extrapolated to 350.38psi.
1000 50
Error %
Stress
800 40
600 30
y = 455.74x -0.4074
400 20
200 10
0 0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Mesh Size
Point 5 all stress results show an error of more than 5%. Stresses at all sizes are disqualified as
being over 5%. The ultimate stress is extrapolated to infinity at 0 mesh size.
The ultimate stress allows the cosmos error estimate to be compared to the true error. The true
error is (actual stress - ultimate stress) / ultimate stress.
Point 3
2
Point 4
1.5
0.5
0
0.00% 5.00% 10.00% 15.00% 20.00% 25.00%
True Error
The graph of the predicted error from the Cosmos error functions vs. the true error. The
predicted error is not very good for point 2. Point 2 is a very low stress area of the model. The
results for point 2 change a lot as the mesh size changes and are not accurate.
Point 4
Point 2
Point 2 is a low stress area adjacent to a high stress area (the sharp radius area) on the model. As
the mesh is refined, small changes in the mesh in the radius area have large effects on point 2.
The error function can not predict the error for this situation. Areas of low stress are not usually
of interest in pressure vessel finite element studies.
Point 4
1.5
0.5
0
0.00% 0.50% 1.00% 1.50% 2.00% 2.50% 3.00% 3.50%
True Error
1.5
Error Comparison
0.5
Linear (Error
Comparison)
0
0.00% 0.50% 1.00% 1.50% 2.00% 2.50% 3.00% 3.50%
True Error %
The graph shows good correlation between the extrapolated true stress and the reported error
from the Cosmos error function. The true stress is approximately 1.49x the reported stress for
points 1, 3 and 4 for this shape.
Conclusion:
The Cosmos Error function will not work for all locations on a model. For this model,
the error results for Point 2 a low stress area adjacent to a high stress area - were found to be
not useable (the reported stress was too low vs. the real stress). Areas of low stress like this
would not normally be of interest in a pressure vessel study.
As a guideline, at Pressure Vessel Engineering we consider using caution when viewing
results at a node when there are elements within 1 node that have errors over 5%.
Laurence Brundrett
Pressure Vessel Engineering Ltd.
Dec 16 2008
Rev 1.